Studying the hydrodynamic regime of groundwater, particularly the formation of new freshwater reserves and monitoring their quantitative and qualitative indicators, holds significant scientific and practical importance. In this research, the main causes of groundwater level variations in two-layer media were analyzed, including precipitation and evaporation, water extraction and recharge sources, geological structure and interlayer permeability characteristics, hydraulic gradient, flow directions, water infiltration due to irrigation, filtration coefficient, porosity, layer thickness, drainage, and artesian wells. To accurately model the variations in the water levels of unconfined and confined aquifers, a mathematical model was developed that accounts for the physical-geological and hydrogeological parameters of the study area. The problem was described through mathematical and numerical modeling of geofiltration and geomigration processes. The model is represented by nonlinear differential equations, which do not have an analytical solution due to the presence of free variables. Therefore, a fully stable numerical solution scheme based on high-precision approximation was proposed. The solutions were obtained using iterative calculations and forward-backward substitution methods. The research also considered parameters, such as soil density, effective porosity, and third-order open boundary conditions. As a result, the model provides reliable forecasts for identifying groundwater movement and quality changes. This method offers practical solutions that can be applied in water resource management and planning.
Daliev et al. (Mon,) studied this question.