The performance of time series models is assessed using synthetic head series simulated with a numerical model that solves Richards' equation for variably saturated flow. Heads were simulated in a homogeneous unconfined aquifer between two parallel canals; measured daily precipitation and potential evaporation are specified at the land surface and root water uptake is simulated. The head response to a precipitation event is nonlinear and depends on the saturation degree and rainfall before and after the precipitation event while evaporation reduction occurs during summers. Synthetic series were generated for 27 years and three different soil types; the unsaturated zone thickness varies between 0 and >5 m. The synthetic head series were simulated with a linear and nonlinear time series model. Performance of a linear time series model with four parameters, using a scaled Gamma response, gave R² values ranging from 0.67 to 0.96. The nonlinear time series model with five parameters simulates recharge using a root zone reservoir after which the head response to recharge is simulated with a scaled Gamma response function. The nonlinear time series model was able to simulate all synthetic head series very well with R² values above 0.9 for almost all models. The head response of the nonlinear model to a single precipitation event compares well to the response of the variably saturated groundwater model. The provided scripts may be used to simulate synthetic head series for other climates or for systems with additional complexity to assess the performance of other data-driven models.

Time series analysis is an increasingly popular method to analyze heads measured in an observation well. Common applications include the quantification of the effect of different stresses (rainfall, pumping, etc.), and the detection of trends and outliers. Pastas is a new and open source Python package for the analysis of hydrogeological time series. The objective of Pastas is twofold: to provide a scientific framework to develop and test new methods, and to provide a reliable ready-to-use software tool for groundwater practitioners. Transfer function noise modeling is applied using predefined response functions. For example, the head response to rainfall is simulated through the convolution of measured rainfall with a Gamma response function. Pastas models are created and analyzed through scripts, ensuring reproducibility and providing a transparent report of the entire modeling process. A Pastas model can be constructed in seven simple steps: import Pastas, read the time series, create a model, specify the stresses and the types of response functions, estimate the model parameters, visualize output, and analyze the results. These seven steps, including the corresponding Python code, are applied to investigate how rainfall and reference evaporation can explain measured heads in an observation well in Kingstown, Rhode Island, USA. The second example demonstrates the use of scripts to analyze a large number of observation wells in batch to estimate the extent of the drawdown caused by a well field in the Netherlands. Pastas is free and open source software available under the MIT-license at http://github.com/pastas/pastas.

An approach is presented to determine groundwater flow in unconsolidated aquifers with a heat pulse response test using a heating cable and a fiber-optic cable. The cables are installed together using direct push so that the cables are in direct contact with the aquifer. The temperature response is measured for multiple days along the fiber-optic cable with Distributed Temperature Sensing (DTS). The new approach fits a two-dimensional analytical solution to the temperature measurements, so that the specific discharge can be estimated without knowledge of the position of the fiber-optic cable relative to the heating cable. Two case studies are presented. The first case study is at a managed aquifer recharge system where fiber-optic cables are inserted 15 m deep at various locations to test the fitting procedure. Similar and relatively large specific discharges are found at the different locations with little vertical variation (0.4–0.6 m/day). The second case study is at a polder, where the water level is maintained 2 m below the surrounding lakes, resulting in significant groundwater flow. The heating and fiber-optic cables are inserted to a depth of 45 m. The specific discharge varies 0.07–0.1 m/day and is significantly larger in a thin layer at 30-m depth. It is shown with numerical experiments that the estimated specific discharge is smoother than in reality due to vertical conduction, but the peak specific discharge is estimated correctly for layers thicker than ∼1.5 m.

An approach is presented to determine the seasonal variations in travel time in a bank filtration system using a passive heat tracer test. The temperature in the aquifer varies seasonally because of temperature variations of the infiltrating surface water and at the soil surface. Temperature was measured with distributed temperature sensing along fiber optic cables that were inserted vertically into the aquifer with direct push equipment. The approach was applied to a bank filtration system consisting of a sequence of alternating, elongated recharge basins and rows of recovery wells. A SEAWAT model was developed to simulate coupled flow and heat transport. The model of a two-dimensional vertical cross section is able to simulate the temperature of the water at the well and the measured vertical temperature profiles reasonably well. MODPATH was used to compute flowpaths and the travel time distribution. At the study site, temporal variation of the pumping discharge was the dominant factor influencing the travel time distribution. For an equivalent system with a constant pumping rate, variations in the travel time distribution are caused by variations in the temperature-dependent viscosity. As a result, travel times increase in the winter, when a larger fraction of the water travels through the warmer, lower part of the aquifer, and decrease in the summer, when the upper part of the aquifer is warmer.

An analytic approach is presented for the simulation of variations in the groundwater level due to temporal variations of recharge in surficial aquifers. Such variations, called groundwater dynamics, are computed through convolution of the response function due to an impulse of recharge with a measured time series of recharge. It is proposed to approximate the impulse response function with an exponential function of time which has two parameters that are functions of space only. These parameters are computed by setting the zeroth and first temporal moments of the approximate impulse response function equal to the corresponding moments of the true impulse response function. The zeroth and first moments are modeled with the analytic element method. The zeroth moment may be modeled with existing analytic elements, while new analytic elements are derived for the modeling of the first moment. Moment matching may be applied in the same fashion with other approximate impulse response functions. It is shown that the proposed approach gives accurate results for a circular island through comparison with an exact solution; both a step recharge function and a measured series of 10 years of recharge were used. The presented approach is specifically useful for modeling groundwater dynamics in aquifers with shallow groundwater tables as is demonstrated in a practical application. The analytic element method is a gridless method that allows for the precise placement of ditches and streams that regulate groundwater levels in such aquifers; heads may be computed analytically at any point and at any time. The presented approach may be extended to simulate the effect of other transient stresses (such as fluctuating surface water levels or pumping rates), and to simulate transient effects in multi-aquifer systems.