Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Unsaturated soil”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Snodgrass Mountain in East River Watershed, Colorado 2020-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors at Snodgrass Mountain. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format SG-X-Y, where SG refers to Snodgrass Mountain, X is the location identifier, and Y is the depth identifier at the same X (shallow Y=1). Specifically, SG-EHS is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and SG-ERTn (upslope n=1) are points along the Snodgrass electrical resistivity tomography transect not associated with the existing site names in the directory. The sample and location information can be found in metadata.csv. Sampling and Measurements Each sample falls into one of the three sampling methods – (1) intact cores, (2) repacked samples, or (3) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. Both intact cores and repacked samples were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores and in-situ soil sensors could be prioritized because these sampling methods are less destructive. While the repacked samples were packed to the target bulk density (estimated post-sampling, when sample volume was known), these samples had altered pore structures. Nevertheless, intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method also has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2026

Soil Water Retention and Hydraulic Conductivity Data and Model at Trail Creek in Taylor River Watershed, Colorado 2024-2025

This data package includes soil water retention and hydraulic conductivity data and model fitting results from measurements of ex-situ soil samples and in-situ soil sensors near Trail Creek. Soil water retention curves (SWRC) characterize soil water content as a function of soil water potential. SWRC depends on soil texture and pore structure and can be used to describe the constraints on biogeochemical processes in terms of soil water availability. In this data package, the sample identification follows the format TR-X-Y, where TR refers to Trail Creek, X is the treatment block identifier, and Y is the location identifier. Specifically, TR-ASCC1 is the control treatment block under the Adaptive Silviculture for Climate Change (ASCC) project, and TR-ASCC2 is the clear-cut treatment block. TR-ASCC-EHSn is associated with ecohydrology sites under the East-Taylor Watershed Community Observatory Sites directory, and TR-ASCC-ERTn (upslope n=1) are ecohydrology sites along the electrical resistivity tomography transects. The sample and location information can be found in metadata.csv, and the data from the soil sensors will be included in a future data version when the observation period becomes sufficiently long for data analysis. Sampling and Measurements Each sample falls into one of the two sampling methods – (1) intact cores or (2) soil sensors – and one of the two measurement methods – (a) laboratory or (b) in-situ. The intact cores were measured using the laboratory methods, which include measurements of soil water potential (HYPROP & WP4C, METER), saturated (KSAT, METER) and unsaturated hydraulic conductivity (HYPROP). The in-situ method uses a pair of co-located soil sensors to measure volumetric water content (TEROS12, METER) and soil water potential (TEROS21, METER), and the hydraulic conductivity was not measured. In comparison, the laboratory methods progress from full saturation to dry conditions, and the in-situ method includes both dry-to-wet and wet-to-dry cycles. The sampling and measurement methods for each sample can be found in metadata.csv, and more information about the measurements is detailed in the Methods section below. Models Retention and hydraulic conductivity data were fitted with four van-Genuchten-type models (specified by “model_name” column in the files): (1) traditional constrained van Genuchten model (“vG_constrained”), (2) traditional unconstrained van Genuchten model (“vG_unconstrained”), (3) PDI-variant of the constrained van Genuchten model (“vG_constrained_PDI”), and (4) PDI-variant of the unconstrained van Genuchten model (“vG_unconstrained_PDI”). The difference between the constrained (1: n) and the unconstrained (2: n, m) van Genuchten models is the number of pore-size distribution parameters in the model equations, giving the unconstrained model more degrees of freedom when fitting the data. Between the traditional and the PDI-variant models, model fitting differs the most at the dry end of the measurements. The traditional models allow infinite suction at the residual water content (water content does not drop below residual water content), and the PDI-variant models enforce a soil water potential value of pF=6.8 (~ -630 MPa) at oven-dryness (water content reaches 0). The inclusion of the van-Genuchten-type models is due to their common application. If other retention models are required, users can access the data in data.csv for further data fitting. More information about the models can be found in the Methods section below. Fitting Tasks The model fitting can be categorized into three levels of tasks (specified by “fitting_task” column in the files). Level 1 (“fit_retention”) only includes retention data fitting (the only level available for the in-situ method). Level 2 (“fit_retention_conductivity”) includes both retention and hydraulic conductivity data fitting, and the saturated hydraulic conductivity (Ks, a parameter of the hydraulic conductivity functions) is fixed by the measurements from KSAT. Level 3 (“fit_retention_conductivity_Ks”) also includes both retention and hydraulic conductivity data fitting, but Ks is a fitted parameter without the constraints from KSAT measurements. Among the same retention models (e.g. vG_constrained models of the same sample), level 1 should produce the best retention data fitting. Level 2 should have the highest misfit of the retention and hydraulic conductivity data, because the retention and hydraulic conductivity functions share common model parameters, and the unsaturated hydraulic conductivity (HYPROP) data fitting is subject to Ks measured independently by KSAT. Level 3 should have mid-level misfits of the retention and hydraulic conductivity data. While level 3 fits the hydraulic conductivity data better than level 2, the fitted Ks value might be unreasonable due to the lack of constraints at the wet end of the measurements. General recommendation when using this data package: (1) Choice of sampling methods: Intact cores might suffer from sample gaps that would lead to overestimation of Ks (sample gaps can be inferred from the “soil_sample_volume” column in metadata.csv when the value is < 249). In-situ method has higher uncertainty in characterizing the wet end of the SWRC because of sensor limitations and the difficulty in reaching full saturation under natural conditions. (2) Choice of fitting tasks: When only retention data is needed, level 1 (“fit_retention”) should be prioritized. When both retention and hydraulic conductivity data are needed, level 2 (“fit_retention_conductivity”) could be prioritized. (3) Choice of models: This could depend on what the downstream models call for. If no specific model is required, model misfit could be used as a ranking criterion. Model misfit values in terms of RMSE can be found in model_parameters.csv. The following files are included in this data package: (1) metadata.csv – This file includes the general information of each sample, including location (description, geocoordinates, elevation), sampling and measurements details (method, depth, time or period, volume, instruments), and soil physical properties (bulk density, saturated hydraulic conductivity, only applicable to physical soil samples). (2) data.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity data of each sample. Column “instrument” specifies the instrument (HYPROP, WP4C, or TEROS) used to perform the measurements. (3) model_fit.csv – This file includes soil water potential, volumetric water content, and unsaturated hydraulic conductivity fitted from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the variable does not apply to that fitting task. (4) model_parameters.csv – This file includes the fitted model parameters, model misfits, and conventional water content thresholds (field capacity and wilting point) from the four models and three fitting tasks. Column “model_name” specifies the retention model used, and “fitting_task” specifies the level of data fitting. Missing values indicate that the parameter does not apply to that model and/or that fitting task. (5) data_Ks.csv – This file includes the saturated hydraulic conductivity measurements from KSAT. (6) /figure/*.png – This folder includes three quick visualizations of the data, retention model fitting results and misfits, and hydraulic conductivity model fitting results, misfits, and parameters. The model fitting results are separated by samples and fitting tasks and colored by models. Zoom-in required. (7) /hyprop/*.bdhx – This folder includes proprietary hyprop files that require the free Labros SoilView-Analysis (METER) to open. Users can explore data fitting using other retention models (i.e. Brooks-Corey, Fredlund-Xing, Kosugi, bimodal models). Be aware that Ks value is pre-entered under “Fitting tab, Conductivity functions parameters” for level 2 fitting. If the value is lost, please refer to metadata.csv under “Ks” column. (8) Six file-level metadata that summarize file, header, column, and variable information of all files. This work was supported by the Watershed Function Science Focus Area at Lawrence Berkeley National Laboratory funded by the US Department of Energy, Office of Science, Biological and Environmental Research under Contract No. DE-AC02-05CH11231.

EARTH SCIENCE > LAND SURFACE > SOILS↗
DOE OSTI2024

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL National Security Test Range

This assessment uses screening level models to calculate potential environmental impacts from proposed tests at two locations at the Idaho National Laboratory (INL) National Security Test Range (NSTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K 2 O, LaBr 3 , KBr, Cu, Zr, F, Ga, Ga 2 O 3 , NaNO 2 , Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface soil impacts were calculated with the computer program Mixing-Cell Model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated based on the assumption that 12 tests are conducted annually using all 11 material types for 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer code input, output, executable files, batch files, scripts, and spreadsheet files, are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS).

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Assessment of Potential Dose and Environmental Impacts from Proposed Testing at the INL Radiological Response Training Range

This assessment uses screening level models to calculate potential environmental impacts from proposed tests at the Idaho National Laboratory (INL) Radiological Response Training Range (RRTR) site. Proposed tests could be conducted using 11 different radioactive material types that include K 2 O, LaBr 3 , KBr, Cu, Zr, F, Ga, Ga 2 O 3 , NaNO 2 , Ga-68, and Tc-99m. The tests could potentially release radioactive material to the atmosphere and radionuclides and other contaminants to the soil, which could leach into the unsaturated zone and migrate to the aquifer. Atmospheric transport of radionuclides to potential human receptors and time-integrated air concentrations were calculated with a Gaussian plume model and three years of hourly meteorological data. Potential surface-soil impacts were calculated with the computer program Mixing-Cell Model (MCM). Groundwater impacts were calculated with the computer programs MCM and GWSCREEN. Radiological doses from potential atmospheric releases were calculated for public receptors off the INL Site and for workers at nearby INL facilities. Results were compared to regulatory dose limits. Maximum potential groundwater concentrations were estimated in the aquifer below the NSTR site and compared to drinking water standards or risk-based screening levels for resident tap water. Soil concentrations were calculated and compared to risk-based screening levels for workers and potential future residents. All impacts were estimated based on the assumption that 12 tests are conducted annually using all 11 material types for a period of 15 years. This document provides the resources to enable a subject matter expert in the field of environmental assessments to replicate the modeling and calculations. The methodology and parameters are presented in the text. All electronic files, including computer code input, output, executable files, batch files, scripts, and spreadsheet files are contained in a zip file that can be accessed by selecting “Additional Information” (select Native File) in the INL Electronic Document Management System (EDMS).

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2021

A Rapid Decision Support Tool for Estimating Impacts of a Vadose Zone Volatile Organic Compound Source on Groundwater and Soil Gas

Diminishing rates of subsurface volatile contaminant removal by soil vapor extraction (SVE) oftentimes warrants an in‐depth performance assessment to guide remedy decision‐making processes. Such a performance assessment must include quantitative approaches to better understand the impact of remaining vadose zone contamination on soil gas and groundwater concentrations. The spreadsheet‐based Soil Vapor Extraction Endstate Tool (SVEET) software functionality has recently been expanded to facilitate quantitative performance assessments. The updated version, referred to as SVEET2, includes expansion of the input parameter ranges for describing a site (site geometry, source characteristics, etc.), an expanded list of contaminants, and incorporation of elements of the Vapor Intrusion Estimation Tool for Unsaturated‐zone Sources software to provide soil gas concentration estimates for use in vapor intrusion evaluation. As part of the update, SVEET2 was used to estimate the impact of a tetrachloroethene (PCE) vadose zone source on groundwater concentrations, comparing SVEET2 results to field‐observed values at an undisclosed site where SVE was recently terminated. PCE concentrations from three separate monitoring wells were estimated by SVEET2 to be within the range of 6.0–6.7 μg/L, as compared to actual field concentrations that ranged from 3 to 11 μg/L PCE. These data demonstrate that SVEET2 can rapidly provide representative quantitative estimates of impacts from a vadose zone contaminant source at field sites. Finally, in the context of the SVE performance assessment, such quantitative estimates provide a basis to support remedial and/or regulatory decisions regarding the continued need for vadose zone volatile organic compound remediation or technical justification for SVE termination, which can significantly reduce the cost to complete for a site.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1988

Estimating water flow through a hillslope using the massively parallel processor

A new two-dimensional model of water flow in a hillslope has been implemented on the Massively Parallel Processor at the Goddard Space Flight Center. Flow in the soil both in the saturated and unsaturated zones, evaporation and overland flow are all modelled, and the rainfall rates are allowed to vary spatially. Previous models of this type had always been very limited computationally. This model takes less than a minute to model all the components of the hillslope water flow for a day. The model can now be used in sensitivity studies to specify which measurements should be taken and how accurate they should be to describe such flows for environmental studies.

Devaney, Judy E.↗
NASA NTRS2019

Constraints on the Metabolic Activity of Microorganisms in Atacama Surface Soils Inferred from Refractory Biomarkers: Implications for Martian Habitability and Biomarker Detection

Dryness is one of the main environmental challenges to microbial survival. Understanding the threshold of microbial tolerance to extreme dryness is relevant to better constrain the environmental limits of life on Earth and critically evaluate long-term habitability models of Mars. Biomolecular proxies for microbial adaptation and growth were measured in Mars-like hyper-arid surface soils in the Atacama Desert that experience only a few millimeters of precipitation per decade, and in biologically active soils a few hundred kilometers away that experience two- to fivefold more precipitation. Diversity and abundance of lipids and other biomolecules decreased with increasing dryness. Cyclopropane fatty acids (CFAs), which are indicative of adaptive response to environmental stress and growth in bacteria, were only detected in the wetter surface soils. The ratio of trans to cis isomers of an unsaturated fatty acid, another bacterial stress indicator, decreased with increasingly dry conditions. Aspartic acid racemization ratios increased from 0.01 in the wetter soils to 0.1 in the driest soils, which is indicative of racemization rates comparable to de novo biosynthesis over long timescales (10,000 years). The content and integrity of stress proteins profiled by immunoassays were additional indicators that biomass in the driest soils is sin situ microbial growth in the driest surface soils of the Atacama, and any metabolic activity is likely to be basal for cellular repair and maintenance only. Our data add to a growing body of evidence that the driest Atacama surface soils represent a threshold for long-term habitability (i.e., growth and reproduction). These results place constraints on the potential for extant life on the surface of Mars, which is 100–1000 times drier than the driest regions in the Atacama.

Wilhelm, Mary Beth↗
DOE OSTI2022

Large-scale paleo water-table rise in a deep desert aquifer recorded by dissolved noble gases

Illuminating past hydrological and climatological conditions in arid regions may provide insights into future trends in groundwater availability. The main objective of this study is to explore paleorecharge processes in the deep regional Nubian Sandstone Aquifer (NSA), which stretches below the hyperarid deserts of the Sinai Peninsula (Egypt) and the Negev (Israel), using dissolved noble gas data. Extremely high amounts of dissolved excess air (Ne concentrations up to 4 × solubility equilibrium) were observed in the ancient ( 81 Kr-depleted) groundwater that was recharged during past pluvial epochs in the basin. Here, the observed unique excess air signal and the clear spatiotemporal structure of the noble gas compositions in the aquifer are hypothesized to reflect two major characteristics of the groundwater system: (a) large-scale, long-term rises in the water table that facilitate the entrapment and dissolution of air bubbles under increased hydrostatic pressure, and (b) the geological settings in the southern recharge area in Sinai, including an almost horizontal position of the layers and intercalation of low-permeability rock formations within the aquifer, which inefficiently reject air bubbles from groundwater, thereby allowing for substantial entrapment of bubbles that dissolve as the water table rises. Enhanced dissolution of entrapped air was also observed in other paleo-groundwater-containing regional sandstone aquifers across North Africa, which, together with the recent findings from the NSA, suggest that largescale water table fluctuations have likely extensively occurred during past pluvial periods over these (and possibly other) arid regions. A preliminary assessment of noble gas recharge temperatures (NGTs) indicates an apparent decoupling of surface and water table temperatures in the case of deep aquifers with a thick (hundreds of m) unsaturated zone. This observed decoupling calls for a re-evaluation of previously obtained NGT records and a need for future work to consider the modification of NGTs due to soil air fractionation and geothermal heating of the underlying deep unsaturated zones.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1993

Moisture Content and Migration Dynamics in Unsaturated Porous Media

Fundamental studies of fluid mechanics and transport in partially saturated soils are presented. Solution of transient diffusion problems in support of the development of probes for the in-situ measurement of moisture content is given. Numerical and analytical methods are used to study the fundamental problem of meniscus and saturation front propagation in geometric models of porous media.

Homsy, G. M.↗
DOE OSTI2024

Decoupling of redox processes from soil saturation in Arctic tundra

Permafrost thaw in warming Arctic landscapes alters hydrology and saturation-driven biogeochemical processes. Models assume that aerobic respiration occurs in drained soils while saturated soils support methanogenesis; however, saturated soils maintain redox gradients that host a range of anaerobic metabolisms. We evaluated how redox potential and redox-active solutes vary with soil moisture in the active layer of permafrost-affected acidic and non-acidic tundra hillslopes. Oxidizing conditions persisted in highly permeable organic horizons of both unsaturated tussock tundra and saturated wet sedge meadows. Redox potential decreased with depth in all soils as increasing soil bulk density restricted groundwater flow and oxygen diffusion. High concentrations of dissolved iron, phosphate, and organic carbon coincided with redox boundaries below the soil surface in acidic tundra, indicating active iron redox cycling and potential release of adsorbed phosphate during iron (oxyhydr)oxide dissolution. In non-acidic tundra, weatherable minerals affected nutrient dynamics more than redox-driven iron cycling, especially in low-lying, saturated areas where thaw reached mineral soils. The role of thaw depth and the ability of saturated soils to maintain oxidizing conditions in organic surface layers highlight the importance of soil physical properties and hydrology in predicting biogeochemical processes and greenhouse gas emissions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Electrical Resistivity Tomography Monitoring of In Situ Soil Flushing at the Hanford 100-K East Area: 100KE Soil Flushing Monitoring

Former operations in the 100-KE Area at the Hanford Site resulted in vadose zone hexavalent chromium that provides a source of groundwater contamination. Previous efforts to remediate vadose zone chromium involved excavation and offsite treatment of contaminated soils. Although much of the chromium was removed, contamination still exists in native soil between the water table and the bottom of the pit, which has since been backfilled. In-situ soil flushing was tested at the 100-KE Area to accelerate the removal of the remaining hexavalent chromium in the vadose zone. Soil flushing works by applying clean flush water at the surface, which mobilizes and transports chromium as it migrates downward to the water table. Once the chromium reaches the water table, it is hydraulically contained, extracted, and treated through pump and treat operations. The efficacy of soil flushing is directly related to the volume of clean water that infiltrates through contaminated soils. Therefore, performance monitoring requires and understanding of flush water migration through the vadose zone. It is challenging to comprehensively assess soil flushing performance via borehole access due to the limited volume of investigation afforded by a given wellbore. As an alternative, 3D time-lapse electrical resistivity tomography (ERT) was tested as a method of monitoring the distribution of flush water over time within the vadose zone. ERT is a method of imaging the bulk electrical conductivity of the subsurface, which is highly dependent on water saturation levels in unconsolidated and unsaturated sediments. Thus, the timing and location of changes in bulk conductivity caused by soil flushing can be used to infer flush water migration pathways. In this report we present results of 3D time-lapse ERT imaging during two separate soil flushing campaigns conducted during the spring and summer months of 2022 and 2023. Results show generally that: 1) Pit backfill materials appear to nominally have larger permeability than native soils. Consequently, the boundary between pit backfill and native soil had a significant impact on flush water migration, causing some flush water to migrate along the pit boundary to the bottom of the pit. 2) Non-uniform flows, likely caused by variations in hydrogeologic properties, developed in the pit backfill materials, resulting in uneven flush water distribution on the southern margin of the soil flushing zone. 3) Redistribution of water at the interface between backfill materials and Hanford formation materials likely facilitated elevated and uniform distribution of flush water within native soils beneath the deeper parts of the pit boundary beneath the center of the flush zone, which presumably overlie soils with elevated chromium contamination. These areas appear to have been infiltrated by higher volumes of flush water than the northern and southern margins of the flush zone. 4) In comparison to 2022, high flush water application rates significantly improved flush water distribution throughout the target flushing zone. 5) Imaging resolution was limited to a depth of approximate 20 meters below the land surface, due primarily to limitations on the lateral extent of the surface ERT array. The ~10 m region of the vadose zone between approximately 130 meters elevation and the water table at approximately 120 meters elevation was unresolved.

100 Area↗
NASA NTRS2016

Xeropreservation of Functionalized Lipid Biomarkers in Hyperarid Soils in the Atacama Desert

Our understanding of long-term organic matter preservation comes mostly from studies in aquatic systems. In contrast, taphonomic processes in extremely dry environments are relatively understudied and are poorly understood. We investigated the accumulation and preservation of lipid biomarkers in hyperarid soils in the Yungay region of the Atacama Desert. Lipids from seven soil horizons in a 2.5 m vertical profile were extracted and analyzed using GC-MS and LC-MS. Diagnostic functionalized lipids and geolipids were detected and increased in abundance and diversity with depth. Deeper clay units contain fossil organic matter (radiocarbon dead) that has been protected from rainwater since the onset of hyperaridity. We show that these clay units contain lipids in an excellent state of structural preservation with functional groups and unsaturated bonds in carbon chains. This indicates that minimal degradation of lipids has occurred in these soils since the time of their deposition between greater than 40,000 and 2 million years ago. The exceptional structural preservation of biomarkers is likely due to the long-term hyperaridity that has minimized microbial and enzymatic activity, a taphonomic process we term xeropreservation (i.e., preservation by drying). The degree of biomarker preservation allowed us to reconstruct major changes in ecology in the Yungay region that reflect a shift in hydrological regime from wet to dry since the early Quaternary. Our results suggest that hyperarid environments, which comprise 7.5 percent of the continental landmass, could represent a rich and relatively unexplored source of paleobiological informationon Earth.

Hyperarid↗
NASA NTRS2016

Xeropreservation of Functionalized Lipid Biomarkers in Hyperarid Soils in the Atacama Desert

Our understanding of long-term organic matter preservation comes mostly from studies in aquatic systems. In contrast, taphonomic processes in extremely dry environments are relatively understudied and are poorly understood. We investigated the accumulation and preservation of lipid biomarkers in hyperarid soils in the Yungay region of the Atacama Desert. Lipids from seven soil horizons in a 2.5 meter vertical profile were extracted and analyzed using GC-MS ( Gas Chromatography-Mass Spectrometry) and LC-MS (Liquid Chromatography-Mass Spectrometry). Diagnostic functionalized lipids and geolipids were detected and increased in abundance and diversity with depth. Deeper clay units contain fossil organic matter (radiocarbon dead) that has been protected from rainwater since the onset of hyperaridity. We show that these clay units contain lipids in an excellent state of structural preservation with functional groups and unsaturated bonds in carbon chains. This indicates that minimal degradation of lipids has occurred in these soils since the time of their deposition between more than 40,000 and up to 2 million years ago. The exceptional structural preservation of biomarkers is likely due to the long-term hyperaridity that has minimized microbial and enzymatic activity, a taphonomic process we term xeropreservation (i.e. preservation by drying). The degree of biomarker preservation allowed us to reconstruct major changes in ecology in the Yungay region that reflect a shift in hydrological regime from wet to dry since the early Quaternary. Our results suggest that hyperarid environments, which comprise 7.5 percent of the continental landmass, could represent a rich and relatively unexplored source of paleobiological information on Earth.

Xeropreservation↗
DOE OSTI2025

Biogeochemical controls on iron speciation and cycling across upland to shoreline gradients in freshwater and estuarine coastal soils (Lake Erie and Chesapeake Bay, United States)

Coastal environments are dynamic interfaces that mediate carbon and nutrient exchanges between terrestrial landscapes and open waters, and understanding the biogeochemical factors controlling these exchanges, particularly iron (Fe) redox transformations, is crucial for predicting coastal ecosystem functions. Here, we investigated the mechanisms controlling Fe speciation changes across upland-to-shoreline gradients in freshwater and estuarine soils using Fe K-edge X-ray absorption spectroscopy, solid and porewater composition analysis, and 16S rRNA sequencing analysis. We show that Fe transformations depend primarily on inundation patterns. In unsaturated uplands, Fe occurs as Fe(III) oxyhydroxides, mainly goethite (9–35 %), Fe(II,III)-phyllosilicates (39–89 %), and Fe(III)-organic species (0–61 %). Soils influenced by estuarine waters exhibit porewater sulfide concentrations reaching up to 221 μM, Fe- and S-cycling bacteria, and up to 81 % pyrite (FeS 2 ), indicating that sulfur-driven redox dynamics control Fe transformations. In lacustrine wetlands, Fe(III) reduction is indicated by porewater Fe(II) concentrations increasing to 1.0–2.1 mM, and ~10–15 % of Fe as Fe(II,III)-(hydr)oxides (green rust), vivianite (Fe 3 (PO 4 ) 2 ·8H 2 O), and/or adsorbed Fe(II) species. EXAFS data also indicate reduction of structural Fe(III) to Fe(II) in phyllosilicates. The presence of Fe- and S-cycling bacteria, as well as sulfide (0–10 μM), suggests that Fe-cycling is microbially driven and potentially coupled with cryptic S-cycling. Fe(II) oxidation was indicated above/near the water table by the presence of Fe(III) oxyhydroxides (ferrihydrite, lepidocrocite). Furthermore, negligible Fe(III) or sulfate reduction was observed at some water-saturated sites located at the upland-wetland transition, likely due to oxic (sub-)surface water inputs. Overall, our results highlight the importance of considering both Fe-speciation and hydro-biogeochemical dynamics when predicting Fe-cycling at coastal interfaces.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Benzene homologues contaminants in a former herbicide factory site: distribution, attenuation, risk, and remediation implication

Benzene homologues often used as organic raw materials or as detergents in chemical industry are prone to accidental release into the environment which can cause serious long-term soil pollutions. In a large former herbicide factory site, we investigated 43 locations for benzene homologues contaminations in soil, soil gas, and groundwater; and studied the hydrogeological conditions. An inverse distance weighted interpolation method was employed to determine the pollutants three-dimensional spatial distribution in the soils. Results showed that benzene homologues residues were mainly originated from the herbicide production workshop; and that the pollution had horizontally expanded at the deeper soil layer. Contaminants had already migrated 15m downward from ground surface. Contaminant phase distribution study showed that NAPL was the primary phase (>99%) for the pollutants accumulated in the unsaturated zone; while it had not migrated to groundwater. The primary mechanism for contaminant transport and attenuation included dissolution of “occluded” NAPL into pore water and pollutant volatilization into soil pore space. Risk assessment revealed that the pollutants brought unacceptable high carcinogenic and non-carcinogenic risks to public health. In order to convert this former chemical processing factory site into a residential area, a remediation to the polluted production workshop sites is urgently required.

Yang, Shuo↗
NASA NTRS1998

A Catchment-Based Land Surface Model for GCMs and the Framework for its Evaluation

A new GCM-scale land surface modeling strategy that explicitly accounts for subgrid soil moisture variability and its effects on evaporation and runoff is now being explored. In a break from traditional modeling strategies, the continental surface is disaggregated into a mosaic of hydrological catchments, with boundaries that are not dictated by a regular grid but by topography. Within each catchment, the variability of soil moisture is deduced from TOP-MODEL equations with a special treatment of the unsaturated zone. This paper gives an overview of this new approach and presents the general framework for its off-line evaluation over North-America.

Ducharen, A.↗
DOE OSTI2020

Organic amendments change soil organic C structure and microbial community but not total organic matter on sub-decadal scales

Organic C has many benefits for soil, but it is depleted by tillage and crop harvest, and especially so for biofuel crops. Accordingly, strategies such as partially retaining stover or planting a cover crop can help ameliorate the negative effect of C removal. We used a long-term field experiment to study the impacts of stover retention and planting cover crop on soil organic matter (SOM), its extractable components and the soil microbial community. SOM chemical composition characterization was determined by electrospray ionization (ESI) coupled with Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR-MS) in sequential water, methanol (MeOH), and chloroform (CHCl3) extracts. The characteristics of the soil bacterial community were measured by phospholipid fatty acid (PLFA), real-time quantitative PCR, and 16S rRNA gene sequence. The variations in total SOM content, total microbial biomass, and bacterial population were slight among treatments, but SOM chemical compounds, arbuscular mycorrhizal fungi (AMF) biomass, and bacterial structure changed significantly, and especially so in the coupled application of stover retention and cover crop. Specifically, stover retention enriched more lignin-like compounds in soil, whereas cover crop enriched more condensed hydrocarbons, and had more compounds with an aromaticity index (AI) >0.5. The bacterial community was not altered by the cover crop, but the corn stover retention increased the relative abundances of Myxococcales (Deltaproteobacteria) and decreased that of Actinobacteria. Redundancy analysis (RDA) further unveiled that the bacterial community in the stover treatments had significantly positive association with CHCl3-extracted chemical classes, i.e. unsaturated hydrocarbon and lipid with the coupled application (stover and cover crop), and lignin and protein with the single corn stover treatment. Taken together, our study shows how different C addition practices influence the molecular composition of SOM and the structure of soil microbial communities.

Xu, Jiangbin↗
DOE OSTI2020

Changes in water vapor adsorption and water film thickness in clayey materials as a function of relative humidity

Abstract Soil water, with adsorbed water being an important component, plays an important role in fluid flow and chemical movement in the unsaturated zone. The adsorbed water consists of water film and interlayer water, and its content is related to the relative humidity (RH). This study theoretically and experimentally focuses on the variation in adsorbed water content in clayey materials as the RH changes. Based on a slit‐pore model, three types of water (water film, interlayer water, and capillary water) could be present in the water vapor adsorption process. The variation of water film is obtained from the analysis of interfacial forces, and then the total water content can be quantitatively subdivided into these three types of water for different RH values. The slit‐pore adsorption model was used with experimental measurements (water vapor adsorption and N 2 physisorption) to quantify the amounts and interrelationships of these three different types of water to six clayey materials, namely three clay minerals (kaolinite, montmorillonite, and illite–smectite mixed layer) and three clay‐rich sediments from the Jianghan Plain in China. The main results are that (a) the volumes of water vapor adsorption are much greater than those of N 2 adsorption except for kaolinite; (b) interlayer water largely dominates the growth of total water content in montmorillonite and results in a concentration of 0.090 ml g −1 at 95% RH; and (c) the minimum thickness of water film is calculated to be 0.23 nm, and the maximum value is one‐third of the pore width by considering the interfacial forces.

54 ENVIRONMENTAL SCIENCES↗