Analytical solutions for flow and advective solute transport in unconfined watershed aquifers with depth-dependent hydraulic conductivity

Вячеслав Гениевич Румынин, Антон Михайлович Никуленков, Леонид Наумович Синдаловский

Research output: Contribution to journalArticlepeer-review


We present analytical and semi-analytical flow and transport models for 2D (x–z vertical plane) flow and particle transport in stratified rock of sedimentary type or vertically non-uniform bedrock aquifers, which are essentially phreatic. The main focus is on the flow induced by the areal discharge in watershed aquifers in hard-rock formations (e.g. granite, gneiss), the weathering profiles of which comprise two or more stratiform zones differing in hydraulic conductivity, k, which thus can be represented by a decay function of the vertical coordinate, i.e., the depth, k(z). Several models of depth-dependent hydraulic conductivity (stepwise, power-law, exponential) have been tested to show the sensitivity of groundwater flow paths and ages (transit or residence time distributions) to various geological and hydraulic simplifications related to rock heterogeneity. The mathematical setup handles the Dupuit–Forchheimer approximation the validity of which, as well as the suggested solutions, were verified via comparison with numerical results obtained using MODFLOW–PMPATH packages.
Original languageEnglish
Article number127116
Number of pages17
JournalJournal of Hydrology
Issue number127116
StatePublished - Dec 2021

Scopus subject areas

  • Water Science and Technology


  • Depth-dependent hydraulic conductivity
  • Recharge
  • Streamline pattern
  • TTD
  • Unconfined aquifer

Cite this