In nature, aquifers are usually composed of distinct kinds of media, i.e., heterogeneous domains rather than homogeneous domains. Groundwater level and flow changes in such domains are more complicated than those in homogeneous domains; thus, building a mathematical model for addressing groundwater flow in heterogeneous aquifers is the present research goal. In conventional research on similar topics, many one-dimensional (1D) analytical models have been presented, but it is challenging to simulate real-world scenarios. This study develops a two-dimensional (2D) analytical model for modeling groundwater flow in a conceptual sloping heterogeneous domain imposed by variable recharge. This model can consider distinct slope angles, medium heterogeneity, and any type of lateral recharge for a semi-infinite domain. The results indicate that groundwater level and flow discharge are greatly affected by the abovementioned factors. The recharge intensity significantly affects the peak of the groundwater level. For example, when the recharge rate increases by 30%, the peak water level increases by 50% as the groundwater flows from the sandy loam zone to the loam zone. The loops delineating the relationship between discharge and groundwater level for different bottom slopes cannot become close for heterogeneous aquifers. The presented 2D analytical model can simulate and better predict results of groundwater changes than previous 1D analytical models. Further, this model can simultaneously consider the effect of varying recharge over time and space on groundwater level change.
Citation: Ping-Cheng Hsieh, Po-Wen Yu, Ming-Chang Wu. Analytical modeling of 2D groundwater flow in a semi-infinite heterogeneous domain with variable lateral sources[J]. AIMS Mathematics, 2024, 9(4): 10121-10140. doi: 10.3934/math.2024495
In nature, aquifers are usually composed of distinct kinds of media, i.e., heterogeneous domains rather than homogeneous domains. Groundwater level and flow changes in such domains are more complicated than those in homogeneous domains; thus, building a mathematical model for addressing groundwater flow in heterogeneous aquifers is the present research goal. In conventional research on similar topics, many one-dimensional (1D) analytical models have been presented, but it is challenging to simulate real-world scenarios. This study develops a two-dimensional (2D) analytical model for modeling groundwater flow in a conceptual sloping heterogeneous domain imposed by variable recharge. This model can consider distinct slope angles, medium heterogeneity, and any type of lateral recharge for a semi-infinite domain. The results indicate that groundwater level and flow discharge are greatly affected by the abovementioned factors. The recharge intensity significantly affects the peak of the groundwater level. For example, when the recharge rate increases by 30%, the peak water level increases by 50% as the groundwater flows from the sandy loam zone to the loam zone. The loops delineating the relationship between discharge and groundwater level for different bottom slopes cannot become close for heterogeneous aquifers. The presented 2D analytical model can simulate and better predict results of groundwater changes than previous 1D analytical models. Further, this model can simultaneously consider the effect of varying recharge over time and space on groundwater level change.
[1] | E. I. Anderson, An analytical solution representing groundwater-surface water interaction, Water Resour. Res., 39 (2003), 1071. https://doi.org/10.1029/2002WR001536 doi: 10.1029/2002WR001536 |
[2] | C. H. Lee, W. P. Chen, R. H. Lee, Estimation of groundwater recharge using water balance coupled with base-flow-record estimation and stable-base-flow analysis, Environ. Geol., 51 (2006), 73–82. https://doi.org/10.1007/s00254-006-0305-2 doi: 10.1007/s00254-006-0305-2 |
[3] | K. Y. Ke, Application of an integrated surface water-groundwater model to multi-aquifers modeling in Choushui River alluvial fan, Taiwan, Hydrol. Process., 28 (2014), 1409–1421. https://doi.org/10.1002/hyp.9678 doi: 10.1002/hyp.9678 |
[4] | J. Kong, P. Xin, G. F. Hua, Z. Y. Luo, C. J. Shen, D. Chen, et al., Effects of vadose zone on groundwater table fluctuations in unconfined aquifers, J. Hydrol., 528 (2015), 397–407. https://doi.org/10.1016/j.jhydrol.2015.06.045 doi: 10.1016/j.jhydrol.2015.06.045 |
[5] | Z. Zomlot, B. Verbeiren, M. Huysmans, O. Batelaan, Spatial distribution of groundwater recharge and base flow: assessment of controlling factors, J. Hydrol., 4 (2015), 349–368. https://doi.org/10.1016/j.ejrh.2015.07.005 doi: 10.1016/j.ejrh.2015.07.005 |
[6] | J. F. Águila, J. Samper, B. Pisani, Parametric and numerical analysis of the estimation of groundwater recharge from water-table fluctuations in heterogeneous unconfined aquifers, Hydrogeol. J., 27 (2019), 1309–1328. https://doi.org/10.1007/s10040-018-1908-x doi: 10.1007/s10040-018-1908-x |
[7] | A. Mahdavi, Response of triangular-shaped leaky aquifers to rainfall-induced groundwater recharge: an analytical study, Water Resour. Manage., 33 (2019), 2153–2173. https://doi.org/10.1007/s11269-019-02234-7 doi: 10.1007/s11269-019-02234-7 |
[8] | S. Kar, J. P. Maity, J. S. Jean, C. C. Liu, B. Nath, H. J. Yang, et al., Arsenic-enriched aquifers: occurrences and mobilization of arsenic in groundwater of Ganges Delta Plain, Barasat, West Bengal, India, Appl. Geochem., 25 (2010), 1805–1814. https://doi.org/10.1016/j.apgeochem.2010.09.007 doi: 10.1016/j.apgeochem.2010.09.007 |
[9] | M. M. Sedghi, H. Zhan, Groundwater dynamics due to general stream fluctuations in an unconfined single or dual-porosity aquifer subjected to general areal recharge, J. Hydrol., 574 (2019), 436–449. https://doi.org/10.1016/j.jhydrol.2019.04.052 doi: 10.1016/j.jhydrol.2019.04.052 |
[10] | N. Pastore, C. Cherubini, A. Doglioni, C. I. Giasi, V. Simeone, A novel approach to model the hydrodynamic response of the surficial level of the Ionian multilayered aquifer during episodic rainfall events, Water, 12 (2020), 2916. https://doi.org/10.3390/w12102916 doi: 10.3390/w12102916 |
[11] | Y. Xin, Z. Zhou, M. Li, C. Zhuang, Analytical solutions for unsteady groundwater flow in an unconfined aquifer under complex boundary conditions, Water, 12 (2020), 75. https://doi.org/10.3390/w12010075 doi: 10.3390/w12010075 |
[12] | M. C. Wu, P. C. Hsieh, Variation of groundwater flow caused by any spatiotemporally varied recharge, Water, 12 (2020), 287. https://doi.org/10.3390/w12010287 doi: 10.3390/w12010287 |
[13] | Y. Zheng, M. Yang, H. Liu, Horizontal two-dimensional groundwater-level fluctuations in response to the combined actions of tide and rainfall in an unconfined coastal aquifer, Hydrogeol. J., 30 (2022), 2509–2518. https://doi.org/10.1007/s10040-022-02564-8 doi: 10.1007/s10040-022-02564-8 |
[14] | W. H. Hassan, H. H. Hussein, B. K. Nile, The effect of climate change on groundwater recharge in unconfined aquifers in the western desert of Iraq, Groundwater Sustainable Dev., 16 (2022), 100700. https://doi.org/10.1016/j.gsd.2021.100700 doi: 10.1016/j.gsd.2021.100700 |
[15] | W. Tao, F. Shao, L. Su, Q. Wang, B. Zhou, Y. Sun, An analytical model for simulating the rainfall-interception-infiltration-runoff process with non-uniform rainfall, J. Environ. Manage., 344 (2023), 118490. https://doi.org/10.1016/j.jenvman.2023.118490 doi: 10.1016/j.jenvman.2023.118490 |
[16] | P. C. Hsieh, M. C. Wu, Changes in groundwater flow in an unconfined aquifer adjacent to a river under surface recharge, Hydrol. Sci. J., 68 (2023), 920–937. https://doi.org/10.1080/02626667.2023.2193295 doi: 10.1080/02626667.2023.2193295 |
[17] | S. E. Serrano, The Theis solution in heterogeneous aquifers, Groundwater, 35 (1997), 463–467. https://doi.org/10.1111/j.1745-6584.1997.tb00106.x doi: 10.1111/j.1745-6584.1997.tb00106.x |
[18] | T. Scheibe, S. Yabusaki, Scaling of flow and transport behavior in heterogeneous groundwater systems, Adv. Water Resour., 22 (1998), 223–238. https://doi.org/10.1016/S0309-1708(98)00014-1 doi: 10.1016/S0309-1708(98)00014-1 |
[19] | P. M. Meier, J. Carrera, X. Sanchez-Vila, A numerical study on the relationship between transmissivity and specific capacity in heterogeneous aquifers, Groundwater, 37 (1999), 611–617. https://doi.org/10.1111/j.1745-6584.1999.tb01149.x doi: 10.1111/j.1745-6584.1999.tb01149.x |
[20] | C. L. Winter, D. M. Tartakovsky, Groundwater flow in heterogeneous composite aquifers, Water Resour. Res., 38 (2002), 23-1-23-11. https://doi.org/10.1029/2001WR000450 doi: 10.1029/2001WR000450 |
[21] | K. Hemker, M. Bakker, Analytical solutions for whirling groundwater flow in two-dimensional heterogeneous anisotropic aquifers, Water Resour. Res., 42 (2006), W12419. https://doi.org/10.1029/2006WR004901 doi: 10.1029/2006WR004901 |
[22] | X. Sanchez-Vila, A. Guadagnini, J. Carrera, Representative hydraulic conductivities in saturated groundwater flow, Rev. Geophys., 44 (2006), RG3002. https://doi.org/10.1029/2005RG000169 doi: 10.1029/2005RG000169 |
[23] | M. Huysmans, A. Dassargues, Application of multiple-point geostatistics on modelling groundwater flow and transport in a cross-bedded aquifer, In: P. Atkinson, C. D. Lloyd, geoENV VⅡ-Geostatistics for environmental applications, Springer, 16 (2010), 135–190. https://doi.org/10.1007/978-90-481-2322-3_13 |
[24] | M. H. Chuang, C. S. Huang, G. H. Li, H. D. Yeh, Groundwater fluctuations in heterogeneous coastal leaky aquifer systems, Hydrol. Earth Syst. Sci., 14 (2010), 1819–1826. https://doi.org/10.5194/hess-14-1819-2010 doi: 10.5194/hess-14-1819-2010 |
[25] | V. A. Zlotnik, M. B. Cardenas, D. Toundykov, Effects of multiscale anisotropy on basin and hyporheic groundwater flow, Groundwater, 49 (2011), 576–583. https://doi.org/10.1111/j.1745-6584.2010.00775.x doi: 10.1111/j.1745-6584.2010.00775.x |
[26] | X. Liang, Y. K. Zhang, Analytic solutions to transient groundwater flow under time-dependent sources in a heterogeneous aquifer bounded by fluctuating river stage, Adv. Water Resour., 58 (2013), 1–9. https://doi.org/10.1016/j.advwatres.2013.03.010 doi: 10.1016/j.advwatres.2013.03.010 |
[27] | S. K. Das, S. J. Ganesh, T. S. Lundström, Modeling of a groundwater mound in a two-dimensional heterogeneous unconfined aquifer in response to precipitation recharge, J. Hydrol. Eng., 20 (2015), 04014081. https://doi.org/10.1061/(ASCE)HE.1943-5584.0001071 doi: 10.1061/(ASCE)HE.1943-5584.0001071 |
[28] | Q. Wang, H. Zhan, Z. Tang, Two-dimensional flow response to tidal fluctuation in a heterogeneous aquifer-aquitard system, Hydrol. Process., 29 (2015), 927–935. https://doi.org/10.1002/hyp.10207 doi: 10.1002/hyp.10207 |
[29] | P. C. Hsieh, P. C. Lee, Analytical modeling of groundwater flow of vertically multilayered soil stratification in response to temporally varied rainfall recharge, Appl. Math. Modell., 96 (2021), 584–597. https://doi.org/10.1016/j.apm.2021.03.030 doi: 10.1016/j.apm.2021.03.030 |
[30] | A. Hartmann, T. Gleeson, Y. Wada, T. Wagener, Enhanced groundwater recharge rates and altered recharge sensitivity to climate variability through subsurface heterogeneity, Proc. Natl. Acad. Sci., 114 (2017), 2842–2847. https://doi.org/10.1073/pnas.1614941114 doi: 10.1073/pnas.1614941114 |
[31] | S. K. Joshi, S. Gupta, R. Sinha, A. L. Densmore, S. P. Rai, S. Shekhar, et al., Strongly heterogeneous patterns of groundwater depletion in Northwestern India. J. Hydrol., 598 (2021) 126492. https://doi.org/10.1016/j.jhydrol.2021.126492 doi: 10.1016/j.jhydrol.2021.126492 |
[32] | L. Wang, C. Dai, L. Xue, A semianalytical model for pumping tests in finite heterogeneous confined aquifers with arbitrarily shaped boundary, Water Resour. Res., 54 (2018), 3207–3216. https://doi.org/10.1002/2017WR022217 doi: 10.1002/2017WR022217 |
[33] | T. G. Chapman, Modeling groundwater flow over sloping beds, Water Resour. Res., 16 (1980), 1114–1118. https://doi.org/10.1029/WR016i006p01114 doi: 10.1029/WR016i006p01114 |
[34] | R. K. Bansal, Groundwater flow in sloping aquifer under localized transient recharge: analytical study, J. Hydraul. Eng., 139 (2013), 1165–1174. https://doi.org/10.1061/(ASCE)HY.1943-7900.0000784 doi: 10.1061/(ASCE)HY.1943-7900.0000784 |
[35] | M. A. Marino, Water-table fluctuation in semipervious stream-unconfined aquifer systems, J. Hydrol., 19 (1973), 43–52. https://doi.org/10.1016/0022-1694(73)90092-9 doi: 10.1016/0022-1694(73)90092-9 |
[36] | N. E. Verhoest, P. A. Troch, Some analytical solutions of the linearized Boussinesq equation with recharge for a sloping aquifer, Water Resour. Res., 36 (2000), 793–800. https://doi.org/10.1029/1999WR900317 doi: 10.1029/1999WR900317 |
[37] | P. A. Troch, E. van Loon, A. Hilberts, Analytical solutions to a hillslope-storage kinematic wave equation for subsurface flow, Adv. Water Resour., 25 (2002), 637–649. https://doi.org/10.1016/S0309-1708(02)00017-9 doi: 10.1016/S0309-1708(02)00017-9 |
[38] | R. K. Bansal, S. K. Das, Analytical study of water table fluctuation in unconfined aquifers due to varying bed slopes and spatial location of the recharge basin, J. Hydrol. Eng., 15 (2010), 909–917. https://doi.org/10.1061/(ASCE)HE.1943-5584.0000267 doi: 10.1061/(ASCE)HE.1943-5584.0000267 |
[39] | R. K. Bansal, Groundwater fluctuations in sloping aquifers induced by time-varying replenishment and seepage from a uniformly rising stream, Transp. Porous Media, 94 (2012), 817–836. https://doi.org/10.1007/s11242-012-0026-9 doi: 10.1007/s11242-012-0026-9 |