Spatial connectivity in a highly heterogeneous aquifer: From cores to preferential flow paths

[1] This study investigates connectivity in a small portion of the extremely heterogeneous aquifer at the Macrodispersion Experiment (MADE) site in Columbus, Mississippi. A total of 19 fully penetrating soil cores were collected from a rectangular grid of 4 m by 4 m. Detailed grain size analysis was performed on 5 cm segments of each core, yielding 1740 hydraulic conductivity (K) estimates. Three different geostatistical simulation methods were used to generate 3-D conditional realizations of the K field for the sampled block. Particle tracking calculations showed that the fastest particles, as represented by the first 5% to arrive, converge along preferential flow paths and exit the model domain within preferred areas. These 5% fastest flow paths accounted for about 40% of the flow. The distribution of preferential flow paths and particle exit locations is clearly influenced by the occurrence of clusters formed by interconnected cells with K equal to or greater than the 0.9 decile of the data distribution (10% of the volume). The fraction of particle paths within the high-K clusters ranges from 43% to 69%. In variogram-based K fields, some of the fastest paths are through media with lower K values, suggesting that transport connectivity may not require fully connected zones of relatively homogenous K. The high degree of flow and transport connectivity was confirmed by the values of two groups of connectivity indicators. In particular, the ratio between effective and geometric mean K (on average, about 2) and the ratio between the average arrival time and the arrival time of the fastest particles (on average, about 9) are consistent with flow and advective transport behavior characterized by channeling along preferential flow paths.

[1]  C. Zheng,et al.  Relative importance of dispersion and rate‐limited mass transfer in highly heterogeneous porous media: Analysis of a new tracer test at the Macrodispersion Experiment (MADE) site , 2010 .

[2]  Y. Rubin,et al.  Spatial correlation of permeability in cross‐stratified sediment with hierarchical architecture , 2004 .

[3]  Brian Berkowitz,et al.  Theory of anomalous chemical transport in random fracture networks , 1998 .

[4]  R. Grayson,et al.  Toward capturing hydrologically significant connectivity in spatial patterns , 2001 .

[5]  Mauro Giudici,et al.  Assessment of the role of facies heterogeneity at the fine scale by numerical transport experiments and connectivity indicators , 2010 .

[6]  Georg Teutsch,et al.  Heterogeneity patterns of Quaternary glaciofluvial gravel bodies (SW-Germany): application to hydrogeology , 2003 .

[7]  C. Zheng,et al.  Evaluation of the applicability of the dual‐domain mass transfer model in porous media containing connected high‐conductivity channels , 2007 .

[8]  C. Zheng,et al.  Numerical Simulation of a Natural Gradient Tracer Experiment for the Natural Attenuation Study: Flow and Physical Transport , 2001, Ground water.

[9]  Numerical simulation of tracer tests in heterogeneous aquifer , 1998 .

[10]  Graham E. Fogg,et al.  Groundwater Flow and Sand Body Interconnectedness in a Thick, Multiple-Aquifer System , 1986 .

[11]  J. Jaime Gómez-Hernández,et al.  Modeling tracer transport at the MADE site: The importance of heterogeneity , 2007 .

[12]  F. Alabert,et al.  Non-Gaussian data expansion in the Earth Sciences , 1989 .

[13]  E. Eric Adams,et al.  Field study of dispersion in a heterogeneous aquifer: 1. Overview and site description , 1992 .

[14]  M. Anderson Subsurface Flow and Transport: Characterization of geological heterogeneity , 1997 .

[15]  David A. Benson,et al.  Lagrangian simulation of multidimensional anomalous transport at the MADE site , 2008 .

[16]  Steven M. Gorelick,et al.  Fractional packing model for hydraulic conductivity derived from sediment mixtures , 1995 .

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

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

[19]  Heidi Christiansen Barlebo,et al.  Investigating the Macrodispersion Experiment (MADE) site in Columbus, Mississippi, using a three‐dimensional inverse flow and transport model , 2004 .

[20]  V. Tidwell,et al.  Permeability Upscaling Measured on a Block of Berea Sandstone: Results and Interpretation , 1999 .

[21]  Xavier Sanchez-Vila,et al.  Point-to-point connectivity, an abstract concept or a key issue for risk assessment studies? , 2008 .

[22]  Rina Schumer,et al.  Fractal mobile/immobile solute transport , 2003 .

[23]  Arlen W. Harbaugh,et al.  MODFLOW-2000, The U.S. Geological Survey Modular Ground-Water Model - User Guide to Modularization Concepts and the Ground-Water Flow Process , 2000 .

[24]  Clayton V. Deutsch,et al.  GSLIB: Geostatistical Software Library and User's Guide , 1993 .

[25]  C. Zheng,et al.  Quantifying mass transfer in permeable media containing conductive dendritic networks , 2005 .

[26]  S. Gorelick,et al.  Heterogeneity in Sedimentary Deposits: A Review of Structure‐Imitating, Process‐Imitating, and Descriptive Approaches , 1996 .

[27]  Robert W. Ritzi,et al.  Behavior of indicator variograms and transition probabilities in relation to the variance in lengths of hydrofacies , 2000 .

[28]  G. Fogg,et al.  Role of Molecular Diffusion in Contaminant Migration and Recovery in an Alluvial Aquifer System , 2001 .

[29]  M. Anderson Hydrogeologic facies models to delineate large-scale spatial trends in glacial and glaciofluvial sediments , 1989 .

[30]  Z. Dai,et al.  Modeling Multiscale Heterogeneity and Aquifer Interconnectivity , 2004, Ground water.

[31]  W. D. Carrier Goodbye, Hazen; Hello, Kozeny-Carman , 2003 .

[32]  Chunmiao Zheng,et al.  Analysis of Solute Transport in Flow Fields Influenced by Preferential Flowpaths at the Decimeter Scale , 2003, Ground water.

[33]  Stephen E. Silliman,et al.  Stochastic analysis of paths of high hydraulic conductivity in porous media , 1988 .

[34]  Jack Eggleston,et al.  Identification of large‐scale hydraulic conductivity trends and the influence of trends on contaminant transport , 1998 .

[35]  Jesús Carrera,et al.  Scale effects in transmissivity , 1996 .

[36]  A. Journel,et al.  Entropy and spatial disorder , 1993 .

[37]  Air‐Based Measurements of Permeability in Pebbly Sands , 2008, Ground water.

[38]  C. Deutsch Fortran programs for calculating connectivity of three-dimensional numerical models and for ranking multiple realizations , 1998 .

[39]  Zhenxue Dai,et al.  Improving permeability semivariograms with transition probability models of hierarchical sedimentary architecture derived from outcrop analog studies , 2005 .

[40]  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 .

[41]  Steven F. Carle,et al.  Geologically based model of heterogeneous hydraulic conductivity in an alluvial setting , 1998 .

[42]  C. Zheng,et al.  Investigation of Small‐Scale Preferential Flow with a Forced‐Gradient Tracer Test , 2011, Ground water.

[43]  Jesús Carrera,et al.  On the relationship between indicators of geostatistical, flow and transport connectivity , 2005 .

[44]  Mary P. Anderson,et al.  Simulation of Preferential Flow in Three-Dimensional, Heterogeneous Conductivity Fields with Realistic Internal Architecture , 1996 .

[45]  Timothy D. Scheibe,et al.  Scaling of flow and transport behavior in heterogeneous groundwater systems , 1998 .

[46]  S. Gorelick,et al.  Rate‐limited mass transfer or macrodispersion: Which dominates plume evolution at the macrodispersion experiment (MADE) site? , 2000 .

[47]  David W. Pollock,et al.  User's guide for MODPATH/MODPATH-PLOT, Version 3; a particle tracking post-processing package for MODFLOW, the U.S. Geological Survey finite-difference ground-water flow model , 1994 .

[48]  C. Zheng,et al.  Delineating Alluvial Aquifer Heterogeneity Using Resistivity and GPR Data , 2005, Ground water.

[49]  Carlos Llopis-Albert,et al.  Gradual conditioning of non-Gaussian transmissivity fields to flow and mass transport data: 3. Application to the Macrodispersion Experiment (MADE-2) site, on Columbus Air Force Base in Mississippi (USA) , 2009 .

[50]  L. Gelhar,et al.  Field study of dispersion in a heterogeneous aquifer: 2. Spatial moments analysis , 1992 .

[51]  E. Eric Adams,et al.  Field study of dispersion in a heterogeneous aquifer , 1992 .

[52]  C. Zheng,et al.  A dual‐domain mass transfer approach for modeling solute transport in heterogeneous aquifers: Application to the Macrodispersion Experiment (MADE) site , 2000 .

[53]  Steven M. Gorelick,et al.  Geological modeling of submeter scale heterogeneity and its influence on tracer transport in a fluvial aquifer , 2010 .

[54]  Steven F. Carle,et al.  Connected-network paradigm for the alluvial aquifer system , 2000 .

[55]  C. Tsang,et al.  Flow channeling in heterogeneous fractured rocks , 1998 .

[56]  Graham E. Fogg,et al.  Role of Molecular Diffusion in Contaminant Migration and Recovery in an Alluvial Aquifer System , 2001 .

[57]  Thomas B. Stauffer,et al.  Database for the Second Macrodispersion Experiment (MADE-2). , 1993 .

[58]  Carlos Llopis-Albert,et al.  Gradual conditioning of non-Gaussian transmissivity fields to flow and mass transport data: 1. Theory , 2009 .

[59]  Yoram Rubin,et al.  Simulation of non‐Gaussian space random functions for modeling transport in groundwater , 1991 .

[60]  Steven F. Carle,et al.  Geologic heterogeneity and a comparison of two geostatistical models: Sequential Gaussian and transition probability-based geostatistical simulation , 2007 .

[61]  G. Fogg,et al.  Modeling Spatial Variability with One and Multidimensional Continuous-Lag Markov Chains , 1997 .

[62]  Jesús Carrera,et al.  Transport upscaling in heterogeneous aquifers: What physical parameters control memory functions? , 2008 .

[63]  David A. Benson,et al.  Subordinated advection‐dispersion equation for contaminant transport , 2001 .

[64]  Warren W. Wood,et al.  Large-scale natural gradient tracer test in sand and gravel, Cape Cod, Massachusetts. 1. Experimental design and observed tracer movement , 1991 .

[65]  Alberto Guadagnini,et al.  Relative importance of geostatistical and transport models in describing heavily tailed breakthrough curves at the Lauswiesen site. , 2008, Journal of contaminant hydrology.

[66]  Marco Bianchi,et al.  Lessons Learned from 25 Years of Research at the MADE Site , 2011, Ground water.

[67]  D. Freyberg,et al.  A natural gradient experiment on solute transport in a sand aquifer: 1. Approach and overview of plume movement , 1986 .

[68]  T. Harter Finite-size scaling analysis of percolation in three-dimensional correlated binary Markov chain random fields. , 2005, Physical review. E, Statistical, nonlinear, and soft matter physics.

[69]  Gaisheng Liu,et al.  Limits of applicability of the advection‐dispersion model in aquifers containing connected high‐conductivity channels , 2004 .

[70]  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 .

[71]  R. Ritzi,et al.  Studying the effect of correlation and finite‐domain size on spatial continuity of permeable sediments , 2008 .

[72]  W. W. Wood,et al.  Large-Scale Natural Gradient Tracer Test in Sand and Gravel, , 1991 .