Modeling contaminant transport in a two-aquifer system with an intervening aquitard

Chih Tien Liu, Hund-Der Yeh*, Li-Ming Yeh

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

4 Scopus citations


This study deals with the issue of one-dimensional solute transport in a two-aquifer system, where an aquitard lies between two aquifers. Different from previous studies on analysis of the contaminant transport affected by the presence of an aquitard, we developed a mathematical transport model in an aquifer-aquitard-aquifer system with considering transport of solutes in the aquitard governed by both advection and diffusion. The Laplace-domain solution of the model for concentration distributions is obtained by the Laplace transform technique and its corresponding time-domain results are computed numerically by using Laplace numerical inversion. An explicit finite difference model is also developed to simulate two-dimensional contaminant transport process in the system. The simulated depth-averaged concentrations in the lower and upper aquifers slightly differ from those predicted by the present solution. The results show that the movement of contaminant in the upper aquifer is slowed down considerably due to the advective transport in aquitard. When neglecting the aquitard advection (a zero Peclet number), the concentration level in the lower aquifer will be underestimated, especially at late times. In addition, the contaminant concentration in the lower aquifer increases significantly with aquitard's Peclet number.

Original languageEnglish
Pages (from-to)200-209
Number of pages10
JournalJournal of Hydrology
StatePublished - 30 Aug 2013


  • Analytical solution
  • Aquitard
  • De Hoog et al.'s method
  • Finite difference solution
  • Laplace transform
  • Tauberian theorem

Fingerprint Dive into the research topics of 'Modeling contaminant transport in a two-aquifer system with an intervening aquitard'. Together they form a unique fingerprint.

Cite this