Effects of non‐Gaussian copula‐based hydraulic conductivity fields on macrodispersion

[1] This is an application of a spatial copula model that is fitted to a real world data set. The copula model allows modeling of pure spatial dependence independently of the marginal distribution. Using non-Gaussian copula models it is demonstrated that the spatial dependence structure of the Borden aquifer is significantly non-Gaussian—despite the fact that the Borden aquifer is commonly thought of as a relatively homogeneous porous medium with a small variance of hydraulic conductivity. In addition to evaluating the spatial dependence structure of the Borden hydraulic conductivity data set using copulas, goal of this study is to explore if the structure of the hydraulic conductivity field influences a physical property, such as plume evolution as evaluated by second spatial central moments of concentration fields. For this comparison, two types of hydraulic conductivity fields were fitted to the Borden hydraulic conductivity data set: one with a Gaussian and the other with a non-Gaussian type of dependence. These two types of hydraulic conductivity fields were constructed such that their second-order spatial moments are identical, and hence they cannot be distinguished by semivariogram-based geostatistics. This paper illustrates that the spatial dependence structure of the Borden hydraulic conductivity data set is significantly non-Gaussian. Despite the fact that Borden is a relatively homogeneous porous medium, and despite the fact that both types of spatial fields are not distinguishable by their variograms, the solute transport characteristics based on these two types of isotropic fields differ significantly in two-dimensional settings. The difference is less pronounced in three-dimensions with anisotropy. It is postulated that non-Gaussian spatial dependence of hydraulic conductivity and a more skewed marginal distribution of hydraulic conductivity will have significant implications in the other more heterogeneous aquifers.

[1]  Gedeon Dagan,et al.  Analysis of flow through heterogeneous random aquifers: 2. Unsteady flow in confined formations , 1982 .

[2]  J. Gómez-Hernández,et al.  To be or not to be multi-Gaussian? A reflection on stochastic hydrogeology , 1998 .

[3]  Y. Rubin.  Stochastic modeling of macrodispersion in heterogeneous porous media , 1990 .

[4]  E. Sudicky.  A natural gradient experiment on solute transport in a sand aquifer: Spatial variability of hydraulic conductivity and its role in the dispersion process , 1986 .

[5]  L. Gelhar,et al.  Plume scale-dependent dispersion in heterogeneous aquifers: 1. Lagrangian analysis in a stratified aquifer , 1993 .

[6]  Gedeon Dagan,et al.  Transport in heterogeneous porous formations: Spatial moments, ergodicity, and effective dispersion , 1990 .

[7]  A. Bárdossy.  Copula‐based geostatistical models for groundwater quality parameters , 2006 .

[8]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[9]  R. L. Naff,et al.  High‐resolution Monte Carlo simulation of flow and conservative transport in heterogeneous porous media: 1. Methodology and flow results , 1998 .

[10]  Charles F. Harvey,et al.  When good statistical models of aquifer heterogeneity go bad: A comparison of flow, dispersion, and mass transfer in connected and multivariate Gaussian hydraulic conductivity fields , 2003 .

[11]  R. L. Naff,et al.  Nonreactive and reactive solute transport in three-dimensional heterogeneous porous media: Mean displacement, plume spreading, and uncertainty , 1994 .

[12]  R. Allen‐King,et al.  Why did Sudicky [1986] find an exponential‐like spatial correlation structure for hydraulic conductivity at the Borden research site? , 2006 .

[13]  S. P. Neuman,et al.  A quasi-linear theory of non-Fickian and Fickian subsurface dispersion , 1990 .

[14]  G. Fogg,et al.  Transition probability-based indicator geostatistics , 1996 .

[15]  E. A. Sudicky,et al.  Micro-scale modelling in the study of plume evolution in heterogeneous media , 1987 .

[16]  Allan D. Woodbury,et al.  The geostatistical characteristics of the borden aquifer , 1991 .

[17]  Franklin W. Schwartz,et al.  Mass transport: 1. A stochastic analysis of macroscopic dispersion , 1980 .

[18]  David L. Freyberg,et al.  A natural gradient experiment on solute transport in a sand aquifer: 2. Spatial moments and the advection and dispersion of nonreactive tracers , 1986 .

[19]  M. Sklar.  Fonctions de repartition a n dimensions et leurs marges , 1959 .

[20]  E. G. Vomvoris,et al.  Stochastic analysis of the concentration variability in a three‐dimensional heterogeneous aquifer , 1990 .

[21]  H. Joe.  Multivariate models and dependence concepts , 1998 .

[22]  A. Bárdossy,et al.  Geostatistical interpolation using copulas , 2008 .

[23]  C. Axness,et al.  Three‐dimensional stochastic analysis of macrodispersion in aquifers , 1983 .

[24]  Sebastien Strebelle,et al.  Conditional Simulation of Complex Geological Structures Using Multiple-Point Statistics , 2002 .

[25]  E. Sudicky,et al.  Three-dimensional analysis of variably-saturated flow and solute transport in discretely-fractured porous media , 1996 .

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

[27]  John A. Cherry,et al.  Migration of contaminants in groundwater at a landfill: A case study: 4. A natural-gradient dispersion test , 1983 .

[28]  A. Bárdossy,et al.  Application of Copulas in Geostatistics , 2010 .

[29]  G. Dagan.  Time‐dependent macrodispersion for solute transport in anisotropic heterogeneous aquifers , 1988 .