HR: 1330h
AN: H32A-0530 [PDF]
TI: A Statistical Model for Predicting Unsaturated Hydraulic Properties of Deep Sediments at the Idaho
National Engineering and Environmental Laboratory
AU: * Winfield, K A
EM: kawinfie@usgs.gov
AF: U.S. Geological Survey, 345 Middlefield Road, Menlo Park, CA 94025 United States
AU: Nimmo, J R
EM: jrnimmo@usgs.gov
AF: U.S. Geological Survey, 345 Middlefield Road, Menlo Park, CA 94025 United States
AB:
The development and application of property-transfer functions is an important approach for predicting unsaturated hydraulic
properties from more easily measured bulk properties. At the Idaho National Engineering and Environmental Laboratory (INEEL),
the unsaturated zone is comprised of thick basalt flow sequences interbedded with thinner sedimentary layers. Buried
hazardous waste in the surficial soil is a possible source of contamination to the underlying Snake River Plain aquifer,
which can be as deep as 200 m below land surface. Determining the unsaturated hydraulic properties of the sedimentary layers
is one step in understanding water flow and solute transport processes through this complex unsaturated system. This study
uses multiple linear regression analysis to construct simple property-transfer functions for estimating the water retention
curve for deep sediments at the INEEL. The regression models were developed using laboratory measurements on 109 sediment
core samples collected at depths of 9 m to 175 m at two facilities within the southwestern portion of the INEEL. These data
included water retention measurements, the curve fit parameters for which are the dependent variables of the
property-transfer functions, and bulk properties (such as bulk density and various representations of the particle-size
distribution), which are the potential independent variables.
The Rossi-Nimmo junction model was used to represent the water retention measurements. Three parameters define this retention
curve model: 1) saturated water content ($\theta$$_{sat}$), 2) a scaling parameter for matric pressure ($\psi$$_{o}$), and
3) a curve shape parameter ($\lambda$). The bulk property data and optimized hydraulic parameter values were used to develop
a separate regression model for each parameter. The predicted parameters were then used to calculate the water retention
curve from saturation to oven dryness. A selection process for the independent variables, referred to as "all possible
subsets regression," was used to determine the best predictive model for each hydraulic parameter.
Preliminary regression results show that textural class percentages were consistently better able to explain the hydraulic
parameters than were other potential representations of the particle-size distribution. The adjusted coefficient of
determination (adjusted R$^{2}$) for the best models, which consisted of some linear combination of textural class
percentages and bulk density, ranged between 0.2 and 0.5 when all observations were included in the regression analyses. The
residuals were close to normally distributed and were fairly homoscedastic when plotted versus the predicted dependent
variable values. The low adjusted R$^{2}$ values may indicate that the bulk property data used in calibrating the models are
not sufficient to completely predict the hydraulic parameters or may indicate significant measurement errors in the dependent
or independent variables. Other bulk property data not available for calibrating the property-transfer functions, such as
mineralogy, specific surface areas, or adsorption capacities, might correlate more strongly with the hydraulic parameters,
and thus may be useful in future regression analyses. The property-transfer functions from this study provide a basis for
development of a theoretical model that relies on physical relationships between the pore-size distribution and the bulk
properties of the media and that should be more universal in its application throughout the INEEL and other geographic
locations.
DE: 1829 Groundwater hydrology
DE: 1866 Soil moisture
DE: 1875 Unsaturated zone
SC: Hydrology [H]
MN: 2003 Fall Meeting