No abstract

In the process of groundwater recharge, surface water usually enters an aquifer by passing an overlying unsaturated zone. Little attention has been given to the development of analytical solutions to a coupled unsaturatedsaturated flow model due to localized recharge up to now. This paper develops a mathematical model to depict threedimensional transient unsaturated-saturated flow in an unconfined aquifer with localized recharge on the ground surface. The model contains Richards' equation for unsaturated flow, a flow equation for saturated formation, and the Gardner constitutive model describing the behavior of unsaturated soil properties. Both flow equations are coupled through the continuity conditions of the head and flux at the water table. The semi-analytical solution to the coupled flow model is derived by the methods of Laplace transform and Fourier cosine transform. A sensitivity analysis is performed to explore the head response to the change in each of the aquifer parameters. A quantitative tool is presented to assess the recharge efficiency signifying the percentage of the water from the recharge to the aquifer. We found that the effect of unsaturated flow on the saturated hydraulic head is negligible if two criteria associated with the unsaturated soil properties and initial aquifer thickness are satisfied. The head distributions predicted from the present solution match well with those from finite-difference simulations. The predictions of the present solution also agree well with the observed data from a field experiment at an artificial recharge pond in Fresno County, California.

A mathematical model describing the transient hydraulic head distribution induced by constant-head pumping/injection at a partially penetrating well in a radial two-zone confined aquifer is a mixed-type boundary value problem. The analytical solution of the model is in terms of an improper integral with an integrand having a singularity at the origin. The solution should rely on numerical methods to evaluate the integral and handle the problems of convergence and singularity. This study aims at developing a new approximate solution describing the transient hydraulic head distribution for a constant-head test (CHT) at a partially penetrating well in the aquifer. This approximate solution is acquired based on a time-dependent diffusion layer approximation proposed in the field of electrochemistry. The diffusion layer can be analogous to the radius of influence in the area of well hydraulics. The approximate solution is in terms of modified Bessel functions for aquifers with a partially penetrating well and can reduce to a simpler form in terms of a natural logarithmic function for the case of well full penetration. The predicted hydraulic heads from the present approximate solution are compared with those estimated by the Laplace-domain solution of the model. The result shows that the predicted spatial head distributions are accurate in the formation zone and fairly good in the skin zone. In addition, the present solution gives an accurate temporal head distribution at a specific location when the radius of influence is far away from the observation wells. This newly developed approximate solution has advantages of easy computing and good accuracy from practical viewpoint, and thus is a handy tool to evaluate temporal and spatial hydraulic head distributions for the CHT.

The drawdown solution has been widely used to analyze pumping test data for the determination of aquifer parameters when coupled with an optimization scheme. The solution can also be used to predict the drawdown due to pumping and design the dewatering system. The drawdown solution for flow toward a finite-radius well with a skin zone in a confined aquifer of infinite extent in radial direction had been developed before. To our best knowledge, the drawdown solution in confined aquifers of finite extent so far has never before been presented in the groundwater literature. This article presents a mathematical model for describing the drawdown distribution due to a constant-flux pumping from a finite-radius well with a skin zone in confined aquifers of finite extent. The analytical solution of the model is developed by applying the methods of Laplace transforms and Bromwich contour integral. This solution can be used to investigate the effects of finite boundary and conductivity ratio on the drawdown distribution. In addition, the inverse relationship between Laplace- and time-domain variables is used to develop the large time solution which can reduce to the Thiem solution if there is no skin zone

Abstract. The drawdown solution has been widely used toanalyze pumping test data for the determination of aquifer parameters when coupled with an optimization scheme. The solution can also be used to predict the drawdown due to pumping and design the dewatering system. The drawdown solution for flow toward a finite-radius well with a skin zone in a confined aquifer of infinite extent in radial direction had been developed before. To our best knowledge, the drawdown solution in confined aquifers of finite extent with a skin zone so far has never before been presented in the groundwater literature. This article presents a mathematical model for describing the drawdown distribution due to a constantflux pumping from a finite-radius well with a skin zone in confined aquifers of finite extent. The analytical solution of the model is developed by applying the methods of Laplace transforms, Bromwich contour integral, and residue theorem. This solution can be used to investigate the effects of finite boundary and conductivity ratio on the drawdown distribution. In addition, the inverse relationship between Laplaceand time-domain variables is used to develop the large time solution which can reduce to the Thiem solution if there is no skin zone.