Search for: All records

The Theis equation is an important mathematical model used for analyzing drawdown data obtained from pumping tests to estimate aquifer parameters. Since the Theis model is a nonlinear equation, a complex graphical procedure is employed for fitting this equation to pump test data. This graphical method was originally proposed by Theis in the late 1930s, and since then, all the groundwater textbooks have included this fitting method. Over the past 90 years, every groundwater hydrologist has been trained to use this tedious procedure for estimating the values of aquifer transmissivity (T) and storage coefficient (S). Unfortunately, this mechanical procedure does not provide any intuition for understanding the inherent limitations in this manual fitting procedure. Furthermore, it does not provide an estimate for the parameter error. In this study, we employ the public domain coding platform Python to develop a script, namely, PyTheis, which can be used to simultaneously evaluate T and S values, and the error associated with these two parameters. We solve nine test problems to demonstrate the robustness of the Python script. The test problems include several published case studies that use real field data. Our tests show that the proposed Python script can efficiently solve a variety of pump test problems. The code can also be easily adapted to solve other hydrological problems that require nonlinear curve fitting routines. 

Groundwater recharge moves downward from the land surface and reaches the groundwater to replenish aquifers. Despite its importance, methods to directly measure recharge remain cost and time‐intensive. Recharge is usually estimated using indirect methods, such as the widely used water‐table fluctuation (WTF) method, which is based on the premise that rises in groundwater levels are due to recharge. In the WTF method, recharge is calculated as the difference between the observed groundwater hydrograph and the hydrograph obtained in the absence of water input. The hydrograph in the absence of rise‐producing input is estimated based on a characteristic master recession curve (MRC), which describes an average behavior for a declining water‐table. Previous studies derive MRC using recession data from all seasons. We hypothesize that for sites where groundwater table is shallow, using recession data from periods with high groundwater‐influenced evapotranspiration (ET) rates versus all periods will yield significantly different MRC, and consequently different estimates of recharge. We test this hypothesis and show that groundwater recession rates are significantly greater in warm months when the groundwater‐influenced ET rates are higher. Since obtaining seasonal recession rates is challenging for locations with a limited amount of data and is prohibitive if it is to be obtained for any given season of a particular year, we propose two novel parsimonious methods to obtain recession time constants for distinct seasons. The proposed methods show the potential to significantly improve the estimates of seasonal recession time constants and provide a better understanding of seasonal variations in recharge estimates.