Formulation of the Governing Partial Differential Equation for Horizontal Velocity as a Function of Vertical Depth with Depth-Dependent Hydraulic Conductivity in Confined Aquifer Systems

References

[1] W. M. Alley, R. W. Healy, J. W. LaBaugh, and T. E. Reilly, “Flow and storage in groundwater systems,” Science, vol. 296, no. 5575, pp. 1985–1990, 2002.

[2] A. C. Ekwe, N. N. Onu, and K. M. Onuoha, “Estimation of aquifer hydraulic characteristics from electrical sounding data: The case of middle Imo River basin aquifers, south-eastern Nigeria,” Journal of Spatial Hydrology, vol. 6, no. 2, pp. 121–132, 2006.

[3] R. A. Freeze and J. A. Cherry, Groundwater. Englewood Cliffs, NJ: Prentice-Hall, 1979.

[4] N. H. Ibragimov, Symmetries, Exact Solutions and Conservation Laws. Springer International Publishing, 2016.

[5] J. Jiang, Y. Feng, Z. Jiang, and S. Li, “Simultaneous identification of aquifer parameters and aquitard hydraulic conductivity using a multi-objective optimization approach,” Hydrogeology Journal, vol. 18, no. 5, pp. 1147–1157, 2010.

[6] Y. Jiang, C. Zheng, C. M. Dong, and K. F. Nordstrom, “Quantitative assessment of depth-dependent groundwater flow patterns in beach aquifers,” Journal of Hydrology, vol. 557, pp. 185–199, 2018.

[7] J. J. Kashaigili, “Ground water governance and management in Tanzania,” International Journal of Water Resources Development, vol. 33, no. 4, pp. 655–672, 2017.

[8] D. O. Olago, “Constraints and solutions for groundwater development, supply and governance in urban areas in Kenya,” Hydrogeology Journal, vol. 27, no. 4, pp. 1315–1340, 2019.

[9] J.-C. Comte, R. Cassidy, J. Obando, and others, “Depletion of groundwater resources under rapid urbanisation in Africa: recent and future trends in the Nairobi Aquifer System, Kenya,” Hydrogeology Journal, vol. 29, no. 2, pp. 765–782, 2021.

[10] S. Oiro, H. Comte, C. Soulsby, and D. Tetzlaff, “Using stable water isotopes to identify spatio-temporal controls on groundwater recharge in two contrasting East African aquifer systems,” Hydrological Sciences Journal, vol. 63, no. 6, pp. 862–877, 2018.

[11] T. Patel, M. N. Mehta, and V. H. Pradhan, “A mathematical model of groundwater recharge through inclined porous media,” International Journal of Applied and Computational Mathematics, vol. 3, no. 4, pp. 3579–3599, 2017.

[12] K. Rotzoll and O. A. El-Kadi, “Estimating hydraulic conductivity from specific capacity for Hawaii aquifers, USA,” Hydrogeology Journal, vol. 22, no. 4, pp. 987–1001, 2014.

[13] X. Wang, J. Liu, S. Li, Y. Song, and Y. Wang, “Analytical solutions for non-Darcian flow in 3D, triple-porosity, multiphase systems with nonlinear interaction coefficients,” Journal of Hydrology, vol. 588, p. 125043, 2020.

[14] H. Zhu, L. Song, and Q. Zhang, “Symmetry analysis of Darcy’s law and continuity equation under depth-dependent hydraulic conductivity,” Applied Mathematics and Computation, vol. 167, no. 2, pp. 1110–1120, 2005.

@article{1718825,
      author = {Oyombe Aluala, Vincent Marani},
      title = {Formulation of the Governing Partial Differential Equation for Horizontal Velocity as a Function of Vertical Depth with Depth-Dependent Hydraulic Conductivity in Confined Aquifer Systems},
      journal = {Iconic Research And Engineering Journals},
      year = {2026},
      volume = {9},
      number = {12},
      pages = {1226-1231},
      issn = {2456-8880},
      url = {https://www.irejournals.com/formatedpaper/1718825.pdf},
      abstract = {In confined aquifer systems, hydraulic conductivity varies with depth due to compaction and lithological changes, but the partial differential equations governing horizontal groundwater velocity under such conditions have not been expressed in a form amenable to systematic analytical treatment. This paper formulates the governing PDEs for two commonly used conductivity-depth relationships: the exponential model K(z)=K_0 e^(-αz) and the power-law model K(z)=K_0 z^(-β). By inserting each model into the mass conservation equation coupled with Darcy’s law and carrying out the required differentiation, both equations are reduced to the common form h_xx+h_zz-f(z) h_z=g(z) h_t. For exponential conductivity, f(z)=α is constant and g(z)=(S_s/K_0 ) e^αz grows without bound; for power-law conductivity, f(z)=β/z is singular at the origin and g(z)=(S_s/K_0 ) z^β grows algebraically. Setting the time derivative to zero produces two elliptic equations: one with constant coefficients (exponential) and one retaining the β/z singularity (power-law). A term-by-term comparison reveals three structural contrasts—in the convective coefficient, the diffusivity growth rate, and the steady-state coefficient type—that govern the applicability of different analytical solution methods. Recovery of the horizontal velocity from the head field via Darcy’s law is stated explicitly.},
      keywords = {Confined Aquifer, Depth-Dependent Hydraulic Conductivity, Governing PDE, Exponential Decay; Power-Law Decay, Darcy’s Law, Velocity-Depth Relationship},
      month = {June},
      doi = {https://doi.org/10.64388/IREV9I12-1718825}
  }