An analytical model for solute transport in an infiltration tracer test in soil with a shallow groundwater table

Summary It is recommended that an in-situ infiltration tracer test is considered for simultaneously determining the longitudinal and transverse dispersion coefficients in soil. Analytical solutions have been derived for two-dimensional advective–dispersive transport in a radial geometry in the literature which can be used for interpreting the result of such a tracer test. However, these solutions were developed for a transport domain with an unbounded-radial extent and an infinite thickness of vadose zone which might not be realistically manifested in the actual solute transport during a field infiltration tracer test. Especially, the assumption of infinite thickness of vadose zone should be invalid for infiltration tracer tests conducted in soil with a shallow groundwater table. This paper describes an analytical model for interpreting the results of an infiltration tracer test based on improving the transport domain with a bounded-radial extent and a finite thickness of vadose zone. The analytical model is obtained with the successive application of appropriate integral transforms and their corresponding inverse transforms. A comparison of the newly derived analytical solution against the previous analytical solutions in which two distinct sets of radial extent and thickness of vadose zone are considered is conducted to determine the influence of the radial and exit boundary conditions on the solute transport. The results shows that both the radial and exit boundary conditions substantially affect the trailing segment of the breakthrough curves for a soil medium with large dispersion coefficients. Previous solutions derived for a transport domain with an unbounded-radial and an infinite thickness of vadose zone boundary conditions give lower concentration predictions compared with the proposed solution at late times. Moreover, the differences between two solutions are amplified when the observation positions are near the groundwater table. In addition, we compare our solution against the approximate solutions that derived from the previous analytical solution and has been suggested to serve as fast tools for simultaneously estimating the longitudinal and transverse dispersion coefficients. The results indicate that the approximate solutions offer predictions that are markedly distinct from our solution for the entire range of dispersion coefficient values. Thus, it is not appropriate to use the approximate solution for interpreting the results of an infiltration tracer test.

[1]  Jeroen C. J. H. Aerts,et al.  Estimates of spatial variation in evaporation using satellite‐derived surface temperature and a water balance model , 2008 .

[2]  Jui-Sheng Chen,et al.  Generalized analytical solution for advection-dispersion equation in finite spatial domain with arbitrary time-dependent inlet boundary condition , 2011 .

[3]  Gour-Tsyh Yeh,et al.  On the computation of Darcian velocity and mass balance in the finite element modeling of groundwater flow , 1981 .

[4]  G. T. Yeh AT123D: analytical transient one-, two-, and three-dimensional simulation of waste transport in the aquifer system , 1981 .

[5]  Jui-Sheng Chen,et al.  An analytical model for simulating two-dimensional multispecies plume migration , 2016 .

[6]  I. N. Sneddon The use of integral transforms , 1972 .

[7]  Jui-Sheng Chen,et al.  Exact analytical solutions for two-dimensional advection–dispersion equation in cylindrical coordinates subject to third-type inlet boundary condition , 2011 .

[8]  Vedat Batu,et al.  A generalized three-dimensional analytical solute transport model for multiple rectangular first-type sources , 1996 .

[9]  Todd H. Skaggs,et al.  Analytical solutions of the one-dimensional advection–dispersion solute transport equation subject to time-dependent boundary conditions , 2013 .

[10]  Roberto Cianci,et al.  Some analytical solutions for two-dimensional convection-dispersion equation in cylindrical geometry , 2006, Environ. Model. Softw..

[11]  B. Fu,et al.  A new mobile‐immobile model for reactive solute transport with scale‐dependent dispersion , 2010 .

[12]  Bruce Hunt Contaminant Source Solutions with Scale-Dependent Dispersivities , 1998 .

[13]  David Andrew Barry,et al.  Application of the convection-dispersion model to solute transport in finite soil columns. , 1988 .

[14]  S. Zou,et al.  Two‐Dimensional Dispersivity Estimation Using Tracer Experiment Data , 1994 .

[15]  Chuen-Fa Ni,et al.  Analytical power series solution for contaminant transport with hyperbolic asymptotic distance-dependent dispersivity , 2008 .

[16]  Bojie Fu,et al.  Contaminant transport in soil with depth-dependent reaction coefficients and time-dependent boundary conditions. , 2013, Water research.

[17]  Vedat Batu,et al.  A generalized two-dimensional analytical solution for hydrodynamic dispersion in bounded media with the first-type boundary condition at the source , 1989 .

[18]  H. Vereecken,et al.  Review of Dispersivities for Transport Modeling in Soils , 2007 .

[19]  Yeh,et al.  FEMWATER: A finite element model of WATER flow through saturated-unsaturated porous media: First revision , 1980 .

[20]  Jui-Sheng Chen,et al.  Generalized analytical solutions to sequentially coupled multi-species advective-dispersive transport equations in a finite domain subject to an arbitrary time-dependent source boundary condition , 2012 .

[21]  Chuen-Fa Ni,et al.  A novel method for analytically solving multi-species advective–dispersive transport equations sequentially coupled with first-order decay reactions , 2012 .

[22]  T. P. Clement,et al.  Modeling Multispecies Reactive Transport in Ground Water , 1998 .

[23]  G. Moridis,et al.  The Laplace Transform Finite Difference Method for Simulation of Flow Through Porous Media , 1991 .

[24]  D. A. Barry,et al.  EXIT CONDITION FOR MISCIBLE DISPLACEMENT EXPERIMENTS , 1992 .

[25]  I. Klenk,et al.  Transverse vertical dispersion in groundwater and the capillary fringe. , 2002, Journal of contaminant hydrology.

[26]  Karsten Pruess,et al.  TOUGHREACT User's Guide: A Simulation Program for Non-isothermal Multiphase Reactive geochemical Transport in Variable Saturated Geologic Media , 2004 .

[27]  A two-well method to evaluate transverse dispersivity for tracer tests in a radially convergent flow field , 1999 .

[28]  J. Sauty An analysis of hydrodispersive transfer in aquifers , 1980 .

[29]  D. A. Barry,et al.  Migration of volatile organic contaminations (VOCs) through a deforming clay liner , 2013 .

[30]  O. Paladino,et al.  A New Method for Laboratory Estimation of the Transverse Dispersion Coefficient , 2007, Ground water.

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

[32]  S. Zou,et al.  Estimation of Dispersion Parameters for Two‐Dimensional Plumes , 1993 .

[33]  M. V. Genuchten,et al.  Analytical Solutions for Solute Transport in Three‐Dimensional Semi‐infinite Porous Media , 1991 .

[34]  Thomas Ptak,et al.  Tracer tests for the investigation of heterogeneous porous media and stochastic modelling of flow and transport—a review of some recent developments , 2004 .

[35]  G. Gao,et al.  An analytical solution of two‐dimensional reactive solute transport in an aquifer‐aquitard system , 2009 .

[36]  C. Zheng,et al.  Applied contaminant transport modeling , 2002 .

[37]  B. Hunt Scale-Dependent Dispersion from a Pit , 2002 .

[38]  Todd H. Skaggs,et al.  Analytical Solution for Multi-Species Contaminant Transport in Finite Media with Time-Varying Boundary Conditions , 2010 .

[39]  W. J. Alves,et al.  Analytical solutions of the one-dimensional convective-dispersive solute transport equation , 1982 .

[40]  Xiaoxian Zhang,et al.  An in situ method to measure the longitudinal and transverse dispersion coefficients of solute transport in soil , 2006 .

[41]  Vedat Batu,et al.  A generalized two‐dimensional analytical solute transport model in bounded media for flux‐type finite multiple sources , 1993 .

[42]  Wolfgang Nowak,et al.  Experiments on vertical transverse mixing in a large-scale heterogeneous model aquifer. , 2005, Journal of contaminant hydrology.

[43]  Eungyu Park,et al.  Analytical solutions of contaminant transport from finite one-, two-, and three-dimensional sources in a finite-thickness aquifer. , 2001, Journal of contaminant hydrology.

[44]  Paul P. Wang,et al.  MT3DMS: A Modular Three-Dimensional Multispecies Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems; Documentation and User's Guide , 1999 .

[45]  B. Fu,et al.  A mobile–immobile model with an asymptotic scale-dependent dispersion function , 2012 .

[46]  B. Hunt,et al.  Solutions and verification of a scale-dependent dispersion model. , 2001, Journal of contaminant hydrology.

[47]  Chuen-Fa Ni,et al.  Analytical power series solutions to the two‐dimensional advection–dispersion equation with distance‐dependent dispersivities , 2008 .

[48]  G. L. Moltyaner,et al.  Twin Lake Tracer Tests: Transverse dispersion , 1988 .

[49]  D. Barry,et al.  Solute transport in partially-saturated deformable porous media: Application to a landfill clay liner , 2012 .

[50]  M. N. Özişik Boundary value problems of heat conduction , 1989 .

[51]  Feike J. Leij,et al.  Analytical solutions for non-equilibrium solute transport in three-dimensional porous media , 1993 .