Conditional stochastic analysis of solute transport in heterogeneous, variably saturated soils

A method is developed for the conditional (Monte Carlo) simulation of steady stateflow and transient transport from point sources in variably saturated porous media. It combines the geostatistical method, a linearized approximation of the soil water tension perturbation solution, and afinite element numerical model. The method is used to investigate the usefulness of conditional simulation for predicting solute transport under a variety of sampling network designs applied to a number of hypothetical soils. Saturated hydraulic conductivity data yield the largest reduction of conditional uncertainty in relatively wet soils with mild heterogeneities. In highly heterogeneous soils or under dry conditions, soil water tension data by themselves, taken at a sampling density of one to two correlation scales along the expected mean travel path, can greatly reduce prediction uncertainty about solute concentration. Parameter uncertainty about statistical properties of the independent random variables becomes less important as the number of conditioning data increases. However, even with a very high number of sampling data, uncertainty of predicted concentration levels remains significant.

[1]  Fred C. Schweppe,et al.  Uncertain dynamic systems , 1973 .

[2]  G. Matheron The intrinsic random functions and their applications , 1973, Advances in Applied Probability.

[3]  A. Journel Geostatistics for Conditional Simulation of Ore Bodies , 1974 .

[4]  S. P. Neuman,et al.  A statistical approach to the inverse problem of aquifer hydrology: 1. Theory , 1979 .

[5]  J. P. Delhomme,et al.  Spatial variability and uncertainty in groundwater flow parameters: A geostatistical approach , 1979 .

[6]  Franklin W. Schwartz,et al.  Mass transport: 3. Role of hydraulic conductivity data in prediction , 1981 .

[7]  Franklin W. Schwartz,et al.  mass transport: 2. Analysis of uncertainty in prediction , 1981 .

[8]  D. Myers Matrix formulation of co-kriging , 1982 .

[9]  S. P. Neuman,et al.  Effects of kriging and inverse modeling on conditional simulation of the Avra Valley Aquifer in southern Arizona , 1982 .

[10]  G. Dagan Stochastic modeling of groundwater flow by unconditional and conditional probabilities: 2. The solute transport , 1982 .

[11]  E. G. Vomvoris,et al.  A geostatistical approach to the inverse problem in groundwater modeling (steady state) and one‐dimensional simulations , 1983 .

[12]  G. Dagan Solute transport in heterogeneous porous formations , 1984, Journal of Fluid Mechanics.

[13]  S. P. Neuman,et al.  Estimation of Aquifer Parameters Under Transient and Steady State Conditions: 1. Maximum Likelihood Method Incorporating Prior Information , 1986 .

[14]  J. Bahr,et al.  Consequences of spatial variability in aquifer properties and data limitations for groundwater modelling practice , 1988 .

[15]  Dennis McLaughlin,et al.  Stochastic analysis of nonstationary subsurface solute transport: 1. Unconditional moments , 1989 .

[16]  Maliha S. Nash,et al.  Soil physical properties at the Las Cruces trench site , 1989 .

[17]  D. McLaughlin,et al.  Stochastic analysis of nonstationary subsurface solute transport: 2. Conditional moments , 1989 .

[18]  Yoram Rubin,et al.  Prediction of tracer plume migration in disordered porous media by the method of conditional probabilities , 1991 .

[19]  Thomas Harter,et al.  An efficient method for simulating steady unsaturated flow in random porous media: Using an analytical perturbation solution as initial guess to a numerical model , 1993 .

[20]  G. Destouni Stochastic modelling of solute flux in the unsaturated zone at the field scale , 1993 .

[21]  T. Harter,et al.  A Numerical Model for Water Flow and Chemical Transport in Variably Saturated Porous Media , 1993 .

[22]  R. Srivastava,et al.  Numerical Simulation of the Wicking Effect in Liner Systems , 1994 .

[23]  Sean Hatch,et al.  Joint conditional simulations and the spectral approach for flow modeling , 1994 .

[24]  Thomas Harter Unconditional and conditional simulation of flow and transport in heterogeneous, variably saturated porous media , 1994 .

[25]  S. P. Neuman,et al.  Eulerian-Lagrangian Analysis of Transport Conditioned on Hydraulic Data: 4. Uncertain Initial Plume State and Non-Gaussian Velocities , 1995 .

[26]  S. P. Neuman,et al.  Eulerian‐Lagrangian Analysis of Transport Conditioned on Hydraulic Data: 1. Analytical‐Numerical Approach , 1995 .

[27]  S. P. Neuman,et al.  Eulerian-Lagrangian analysis of transport conditioned on hydraulic data 3. Spatial moments, travel time distribution, mass flow rate, and cumulative release across a compliance surface , 1995 .

[28]  Shlomo P. Neuman,et al.  Eulerian-Lagrangian Analysis of Transport Conditioned on Hydraulic Data: 2. Effects of Log Transmissivity and Hydraulic Head Measurements , 1995 .

[29]  Thomas Harter,et al.  Stochastic analysis of solute transport in heterogeneous, variably saturated soils , 1996 .