Parameter identification for leaky aquifers using global optimization methods

In the past, graphical or computer methods were usually employed to determine the aquifer parameters of the observed data obtained from field pumping tests. Since we employed the computer methods to determine the aquifer parameters, an analytical aquifer model was required to estimate the predicted drawdown. Following this, the gradient‐type approach was used to solve the nonlinear least‐squares equations to obtain the aquifer parameters. This paper proposes a novel approach based on a drawdown model and a global optimization method of simulated annealing (SA) or a genetic algorithm (GA) to determine the best‐fit aquifer parameters for leaky aquifer systems. The aquifer parameters obtained from SA and the GA almost agree with those obtained from the extended Kalman filter and gradient‐type method. Moreover, all results indicate that the SA and GA are robust and yield consistent results when dealing with the parameter identification problems. Copyright © 2006 John Wiley & Sons, Ltd.     

An automated disc infiltrometer for infiltration rate measurements using a microflowmeter

This work presents a new design of disc infiltrometer, which, associated with a microflowmeter (MF) and a solenoid valve set, makes it possible to automate the infiltration rate (Q) measurements at different soil pressure heads (ψ). The MF consists of a 13·8‐cm long and 1·5 mm i.d. pipe, with a pressure transducer connecting the two ends of the MF, inserted in a water‐flow pipe that connects the Mariotte tube and the water‐supply reservoir of the disc infiltrometer. Water flow is calculated from the head losses in the MF. Changes in ψ in the bubble tower, automatically affected when the infiltration rate reaches steady state, are controlled by a datalogger connected to four solenoid valves. The new design was tested in laboratory and field conditions, and the results showed that the MF allows the soil water infiltration rates to be correctly estimated for different soil characteristics. The solenoid valve set plus datalogger system satisfactorily monitored the changes in ψ and allowed the measurement time to be optimized. Copyright © 2011 John Wiley & Sons, Ltd.     

Parameter Estimation for Groundwater Models under Uncertain Irrigation Data

The success of modeling groundwater is strongly influenced by the accuracy of the model parameters that are used to characterize the subsurface system. However, the presence of uncertainty and possibly bias in groundwater model source/sink terms may lead to biased estimates of model parameters and model predictions when the standard regression‐based inverse modeling techniques are used. This study first quantifies the levels of bias in groundwater model parameters and predictions due to the presence of errors in irrigation data. Then, a new inverse modeling technique called input uncertainty weighted least‐squares (IUWLS) is presented for unbiased estimation of the parameters when pumping and other source/sink data are uncertain. The approach uses the concept of generalized least‐squares method with the weight of the objective function depending on the level of pumping uncertainty and iteratively adjusted during the parameter optimization process. We have conducted both analytical and numerical experiments, using irrigation pumping data from the Republican River Basin in Nebraska, to evaluate the performance of ordinary least‐squares (OLS) and IUWLS calibration methods under different levels of uncertainty of irrigation data and calibration conditions. The result from the OLS method shows the presence of statistically significant (p < 0.05) bias in estimated parameters and model predictions that persist despite calibrating the models to different calibration data and sample sizes. However, by directly accounting for the irrigation pumping uncertainties during the calibration procedures, the proposed IUWLS is able to minimize the bias effectively without adding significant computational burden to the calibration processes.     

Optimal pumping locations of skimming wells

Abstract

A real-life problem involving pumping of groundwater from a series of existing wells along a river flood plain underlain with geologically saline water is examined within a conceptual framework. Unplanned pumping results in upconing of saline water. Therefore, it is necessary to determine optimal locations of fixed capacity pumping wells in space and time from a set of pre-selected candidate wells that minimize total salinity concentration in space and time. The nonlinear, non-convex, combinatorial problem involving zero—one decision variables is solved in a simulation—optimization (S/O) framework. Optimization is accomplished by using simulated annealing (SA)—a search algorithm. The computational burden is primarily managed by replacing the numerical model with a surrogate simulator—artificial neural network (ANN). The computational burden is further reduced through intuitive algorithmic guidance. The model results suggest that the skimming wells must be operated from optimal locations such that they are staggered in space and time to obtain least saline water.     

Identification of Groundwater Parameters Using an Adaptative Multiscale Method

The identification of groundwater parameters in heterogeneous systems is a major challenge in groundwater modeling. Flexible parameterization methods are needed to assess the complexity of the spatial distributions of these parameters in real aquifers. In this article, we introduce an adaptative parameterization to identify the distribution of hydraulic conductivity within the large‐scale (4400 km²) Upper Rhine aquifer. The method is based on adaptative multiscale triangulation (AMT) coupled with an inverse problem procedure that identifies the parameters' distributions by reducing the error between measured and simulated heads. The AMT method has the advantage of combining both zonation and interpolation approaches. The AMT method uses area‐based interpolation rather than an interpolation based on stochastic features. The method is applied to a standard 2D groundwater model that takes into account the interactions between the aquifer and surface water bodies, groundwater recharge, and pumping wells. The simulation period covers 204 months, from January 1986 to December 2002. Recordings at 109 piezometers are used for model calibration. The simulated heads are globally quite accurate and reproduce the main dynamics of the system. The local hydraulic conductivities resulting from the AMT method agree qualitatively with existing local experimental observations across the Rhine aquifer.     

The variable projection method for waveform inversion with an unknown source function

This paper compares three alternative algorithms for simultaneously estimating a source wavelet at the same time as an earth model in full‐waveform inversion: (i) simultaneous descent, (ii) alternating descent and (iii) descent with the variable projection method. The latter is a technique for solving separable least‐squares problems that is well‐known in the applied mathematics literature. When applied to full‐waveform inversion, it involves making the source wavelet an implicit function of the earth model via a least‐squares filter‐estimation process. Since the source wavelet becomes purely a function of medium parameters, it no longer needs to be treated as a separate unknown in the inversion. Essentially, the predicted data are projected onto the measured data in a least‐squares sense at every function evaluation, making use of the fact that the filter estimation problem is trivial when compared to the full‐waveform inversion problem. Numerical tests on a simple 1D model indicate that the variable projection method gives the best result; actually producing results in quality that are very similar to control experiments with a known, correct wavelet.     

Capture Versus Capture Zones: Clarifying Terminology Related to Sources of Water to Wells

The term capture, related to the source of water derived from wells, has been used in two distinct yet related contexts by the hydrologic community. The first is a water‐budget context, in which capture refers to decreases in the rates of groundwater outflow and (or) increases in the rates of recharge along head‐dependent boundaries of an aquifer in response to pumping. The second is a transport context, in which capture zone refers to the specific flowpaths that define the three‐dimensional, volumetric portion of a groundwater flow field that discharges to a well. A closely related issue that has become associated with the source of water to wells is streamflow depletion, which refers to the reduction in streamflow caused by pumping, and is a type of capture. Rates of capture and streamflow depletion are calculated by use of water‐budget analyses, most often with groundwater‐flow models. Transport models, particularly particle‐tracking methods, are used to determine capture zones to wells. In general, however, transport methods are not useful for quantifying actual or potential streamflow depletion or other types of capture along aquifer boundaries. To clarify the sometimes subtle differences among these terms, we describe the processes and relations among capture, capture zones, and streamflow depletion, and provide proposed terminology to distinguish among them.     

An inverse procedure to estimate transmissivity from heads and SP signals

Experimental hydraulic heads and electrical (self-potential) signals associated with a pumping test were used in an inverse model to estimate the transmissivity distribution of a real aquifer. Several works reported in the literature show that there is a relatively good linear relationship between the hydraulic heads in the aquifer and electrical signals measured at the ground surface. In this experimental test field, first, the current coupling coefficient was determined by the best fit between experimental and modeled self-potential signals at the end of the pumping phase. Soon afterward, with the hydraulic heads obtained from the self-potential signals, the transmissivity distribution of the aquifer was conditioned by means an inverse model based on the successive linear estimator (SLE). To further substantiate the estimated T field from the SLE analysis, we analyzed the drawdown rate, the derivative of the drawdown with respect to the ln(t), because the drawdown rate is highly sensitive to the variability in the transmissivity field. In our opinion, these results show that self-potential signals allow the monitoring of subsurface flow in the course of pumping experiments, and that electrical potentials serve as a good complement to piezometric observations to condition and characterize the transmissivity distribution of an aquifer.     

INFLUENCE OF PUMPING ON FOUR-COMPONENT BOREHOLE STRAIN AT TAI'AN SEISMIC STATION AND ITS INTERFERENCE MECHANISM

Four-component borehole strainmeter (FCBS) is one kind of high-precision borehole strain observation instruments invented in China. As a kind of near-surface deformation observation instrument, FCBS is also easily disturbed by the external environment factors. As a common factor, pumping has significant influence on FCBS observation. Existing studies mostly identify the pumping interference from the perspective of observation curve morphology, relatively few studies focus on its interference mechanism. In order to truly capture earthquake precursor information, it is necessary to study the interference mechanism. In recent years, RZB-3 type FCBS at Tai'an seismic station has been seriously affected by pumping, so it is necessary and also feasible to study the interference mechanism of pumping. Since the influence of pumping interference on borehole strainmeter is common, this work would be very practical and be used for reference by other borehole strain observation stations. We find that the original observation curves and observed surface strain, shear strain from RZB-3 type FCBS at Tai'an seismic station have the characteristics of synchronous change with the borehole water level, in which the linear correlation coefficient between the two observed shear strain curves and borehole water level reached 0.70 and 0.82 respectively. We further find that the principal strain direction of borehole and borehole water level after normalization meet the nonlinear function as y=1.217arctan(x)^0.224-0.284. The above phenomenon indicates that the observation of RZB-3 type FCBS at Tai'an seismic station is significantly affected by the borehole water level, and the influence is more obvious and the gradient is larger at the stage of low water level. Pumping interference often appears in low water level stage and changes the rock pore pressure state. Statistics show that pumping interference affects the borehole strain state. To investigate the interference mechanism of pumping to RZB-3 type FCBS at Tai'an seismic station, we take a known pumping as an example, in which we study the principal strain state of the borehole in three periods of normal pumping, interruption of pumping and resuming pumping respectively. During each period, we solve 3 parameters of the principal plane strain state, i.e. the maximum principal strain rate, the minimum principal strain rate and the maximum principal strain direction from four observation equations of FCBS by nonlinear iterative least squares algorithm. On the other hand, concentrated load model (CLM) is used to simulate the mechanical mechanism of pumping. Firstly, the depth of FCBS relative to pumping source and the concentrated load at pumping source are inversed, then, the strain state surrounding the pumping well, including the state at RZB-3 borehole, is simulated by forward modeling. By comparing these results, we find that: (1)The concentrated load at pumping source inversed by CLM during periods of normal pumping and resuming pumping are both located at or near the bottom of the pumping well, which is consistent with the actual situation, indicating that mechanism and degree of the influence of pumping on borehole strain are well simulated by CLM. (2)The observed strain state is consistent with the simulation result of pumping interference by forward modeling, indicating that the principal strain state of borehole calculated based on observation of FCBS reflects the true strain state of borehole under different pumping states.(3)The inversed concentrated load at pumping source during pumping periods is significant greater than the load of the pumped water, indicating that the pumping process has more significant influence on the pore pressure of rocks than the load of the pumped water. Even though CLM is an approximate simulation since it's based on some elastic assumptions, the interference mechanism of pumping on RZB-3 type FCBS at Tai'an seismic station is well explained, which is maybe very helpful for studying the influence of pumping interference on other deformation instruments, locating the unknown pumping source and studying the characteristics of pore pressure of rocks.     

Fecal Indicator Bacteria Variability in Samples Pumped from Monitoring Wells

The detection of microbiological contamination in drinking water from groundwater wells is often made with a limited number of samples that are collected using traditional geochemical sampling protocols. The objective of this study is to examine the variability of fecal indicator bacteria, as observed using discrete samples, due to pumping. Two wells were instrumented as multilevel piezometers in a bedrock aquifer, and bacterial enumeration was conducted on a total of 166 samples (for total coliform, fecal coliform, Escherichia coli, and fecal streptococci) using standard membrane filtration methods. Five tests were conducted using pumping rates ranging from 0.3 to 17 L/min in a variety of purging scenarios, which included constant and variable (incremental increase and decrease) flow. The results clearly show a rapid and reproducible, 1 to 2 log‐unit decrease in fecal indicator bacteria at the onset of pumping to stabilized, low‐level concentrations prior to the removal of three to five well volumes. The pumping rate was not found to be correlated with the magnitude of observed bacterial counts. Based on the results, we suggest sampling protocols for fecal indicator bacteria that include multiple collections during the course of pumping, including early‐time samples, and consider other techniques such as microscopic enumeration when assessing the source of bacteria from the well‐aquifer system.     

Field Study of Subsurface Heterogeneity with Steady‐State Hydraulic Tomography

Remediation of subsurface contamination requires an understanding of the contaminant (history, source location, plume extent and concentration, etc.), and, knowledge of the spatial distribution of hydraulic conductivity (K) that governs groundwater flow and solute transport. Many methods exist for characterizing K heterogeneity, but most if not all methods require the collection of a large number of small‐scale data and its interpolation. In this study, we conduct a hydraulic tomography survey at a highly heterogeneous glaciofluvial deposit at the North Campus Research Site (NCRS) located at the University of Waterloo, Waterloo, Ontario, Canada to sequentially interpret four pumping tests using the steady‐state form of the Sequential Successive Linear Estimator (SSLE) ( Yeh and Liu 2000 ). The resulting three‐dimensional (3D) K distribution (or K‐tomogram) is compared against: ( 1 ) K distributions obtained through the inverse modeling of individual pumping tests using SSLE, and ( 2 ) effective hydraulic conductivity (K_eff) estimates obtained by automatically calibrating a groundwater flow model while treating the medium to be homogeneous. Such a K_eff is often used for designing remediation operations, and thus is used as the basis for comparison with the K‐tomogram. Our results clearly show that hydraulic tomography is superior to the inversions of single pumping tests or K_eff estimates. This is particularly significant for contaminated sites where an accurate representation of the flow field is critical for simulating contaminant transport and injection of chemical and biological agents used for active remediation of contaminant source zones and plumes.     

Generalized Radial Flow in Synthetic Flow Systems

Traditional analysis methods used to determine hydraulic properties from pumping tests work well in many porous media aquifers, but they often do not work in heterogeneous and fractured‐rock aquifers, producing non‐plausible and erroneous results. The generalized radial flow model developed by Barker (1988) can reveal information about heterogeneity characteristics and aquifer geometry from pumping test data by way of a flow dimension parameter. The physical meaning of non‐integer flow dimensions has long been a subject of debate and research. We focus on understanding and interpreting non‐radial flow through high permeability conduits within fractured aquifers. We develop and simulate flow within idealized non‐radial flow conduits and expand on this concept by simulating pumping in non‐fractal random fields with specific properties that mimic persistent sub‐radial flow responses. Our results demonstrate that non‐integer flow dimensions can arise from non‐fractal geometries within aquifers. We expand on these geometric concepts and successfully simulate pumping in random fields that mimic well‐test responses seen in the Culebra Dolomite above the Waste Isolation Pilot Plant.     

Redundant and Nonredundant Information for Model Calibration or Hydraulic Tomography

Drawdown data from independent pumping tests have widely been used to validate the estimated hydraulic parameters from inverse modeling or hydraulic tomography (HT). Yet, the independent pumping test has not been clearly defined. Therefore, the goal of this paper is to define this independent pumping test concept, based on the redundant or nonredundant information about aquifer heterogeneity embedded in the observed heads during cross-hole pumping tests. The definition of complete, moderate redundancy and high nonredundancy of information are stipulated using cross-correlation analysis of the relationship between the head and heterogeneity. Afterward, data from numerical experiments and field sequential pumping test campaigns reinforce the concept and the definition.     

Solving Groundwater Flow Problems with Time Series Analysis: You May Not Even Need Another Model

Time series analysis is a data-driven approach to analyze time series of heads measured in an observation well. Time series models are commonly much simpler and give much better fits than regular groundwater models. Time series analysis with response functions gives insight into why heads vary, while such insight is difficult to gain with black box models out of the artificial intelligence world. An important application is to quantify the contributions to the head variation of different stresses on the aquifer, such as rainfall and evaporation, pumping, and surface water levels. Time series analysis may be applied to answer many groundwater questions without the need for a regular groundwater model, such as what is the drawdown caused by a pumping station? Or, how long will it take before groundwater levels recover after a period of drought? Even when a regular groundwater model is needed to solve a groundwater problem, time series analysis can be of great value. It can be used to clean up the data, identify the major stresses on the aquifer, determine the most important processes that affect flow in the aquifer, and give an indication of the fit that can be expected. In addition, it can be used to determine calibration targets for steady-state models, and it can provide several alternative calibration methods for transient models. In summary, the overarching message of this paper is that it would be wise to do time series analysis for any application that uses measured groundwater heads.     

Prospective modelling of 3D hyporheic exchange based on high‐resolution topography and stream elevation

River managers and scientists interested in hyporheic processes need adequate tools for characterizing hyporheic exchange flow (HEF) at local sites where only poor information on subsurface properties are available. This study evaluates a three‐dimensional modelling approach, on the basis of detailed surface parameterization and a simplified subsurface structure, for comparison of potential HEF characteristics at three experimental reaches at the channel‐unit scale. First, calibration is conducted to determine the best fit‐of‐heads given the model simplification, then the structure of residuals are used to evaluate the origin of the misfit, and finally, a sensitivity analysis is conducted to identify inter‐site differences in HEF. Results show that such an approach can highlight potential magnitude differences in HEF characteristics between reaches. The sensitivity analysis is successful in delineating the small area of exchange that remains under conditions of high groundwater discharge. In this case, however, the calibrated model performs poorly in representing the exchange pattern at the sediment–water interface, thus suggesting that the approach is less adequate for a deterministic simulation of observed heads. The summary statistics are in the range of similar published models, for which the reported indicator is the root mean square error on heads normalized by the head drop over the reach. We recommend, however, that modellers use a more comparable indicator, such as a measure of the residuals normalized by a measure of observed vertical head differences. Overall, when subsurface data are unavailable or sparse, a three‐dimensional groundwater model based on high‐resolution topographic data combined with a sensitivity analysis appears as a useful tool for a preliminary characterization of HEF. Copyright © 2013 John Wiley & Sons, Ltd.     

Validation of a soil water balance model using soil water content and pressure head data

The validation of soil water balance models and the evaluation of the quality of the model predictions at field‐scale require time‐series of in situ measured model outputs. In our study, we have validated such a model using a 6‐year period with time‐series of automatically recorded, daily volumetric soil water contents measured with the time‐domain reflectometry with intelligent microelements (TRIME) method and daily pressure heads measured with tensiometers. The comparisons of simulated with measured soil water contents and pressure heads were analysed using the modelling efficiency index (IA) and the square root of the mean square error (RMSE) in order to evaluate the prediction quality of the model. In our study, IA and RMSE, obtained either from the comparison of simulated with measured soil water contents or the comparison of calculated with observed pressure heads, in some cases lead to different results regarding the evaluation of the simulation quality of the soil water balance model. For example, a good fit between simulated and observed soil water contents does not necessarily result in a comparably good fit between the corresponding calculated and measured pressure heads. Therefore, a combined use of both measurement techniques, which takes into account their respective advantages and disadvantages, gives a more complete overview on the simulation quality of the soil water balance model than the single use of one of those techniques. Copyright © 2004 John Wiley & Sons, Ltd.     

Coupled Surface and Groundwater Hydrological Modeling in a Changing Climate

Many current watershed modeling efforts now incorporate surface water and groundwater for managing water resources since the exchanges between groundwater and surface water need a special focus considering the changing climate. The influence of groundwater dynamics on water and energy balance components is investigated in the Snake River Basin (SRB) by coupling the Variable Infiltration Capacity (VIC) and MODFLOW models (VIC‐MF) for the period of 1986 through 2042. A 4.4% increase in base flows and a 10.3% decrease in peak flows are estimated by VIC‐MF compared to the VIC model in SRB. The VIC‐MF model shows significant improvement in the streamflow simulation (Nash‐Sutcliffe efficiency [NSE] of 0.84) at King Hill, where the VIC model could not capture the effect of spring discharge in the streamflow simulation (NSE of ?0.30); however, the streamflow estimates show an overall decreasing trend. Two climate scenarios representing median and high radiative‐forcings such as representative concentration pathways 4.5 and 8.5 show an average increase in the water table elevations between 2.1 and 2.6 m (6.9 and 8.5 feet) through the year 2042. The spatial patterns of these exchanges show a higher groundwater elevation of 15 m (50 feet) in the downstream area and a lower elevation of up to 3 m (10 feet) in the upstream area. Broadly, this study supports results of previous work demonstrating that integrated assessment of groundwater‐surface water enables stakeholders to balance pumping, recharge and base flow needs and to manage the watersheds that are subjected to human pressures more sustainably.     

Influence of a thin veneer of low‐hydraulic‐conductivity sediment on modelled exchange between river water and groundwater in response to induced infiltration

A thin layer of fine‐grained sediment commonly is deposited at the sediment–water interface of streams and rivers during low‐flow conditions, and may hinder exchange at the sediment–water interface similar to that observed at many riverbank‐filtration (RBF) sites. Results from a numerical groundwater‐flow model indicate that a low‐permeability veneer reduces the contribution of river water to a pumping well in a riparian aquifer to various degrees, depending on simulated hydraulic gradients, hydrogeological properties, and pumping conditions. Seepage of river water is reduced by 5–10% when a 2‐cm thick, low‐permeability veneer is present on the bed surface. Increasing thickness of the low‐permeability layer to 0·1 m has little effect on distribution of seepage or percentage contribution from the river to the pumping well. A three‐orders‐of‐magnitude reduction in hydraulic conductivity of the veneer is required to reduce seepage from the river to the extent typically associated with clogging at RBF sites. This degree of reduction is much larger than field‐measured values that were on the order of a factor of 20–25. Over 90% of seepage occurs within 12 m of the shoreline closest to the pumping well for most simulations. Virtually no seepage occurs through the thalweg near the shoreline opposite the pumping well, although no low‐permeability sediment was simulated for the thalweg. These results are relevant to natural settings that favour formation of a substantial, low‐permeability sediment veneer, as well as central‐pivot irrigation systems, and municipal water supplies where river seepage is induced via pumping wells. Published in 2011 by John Wiley & Sons, Ltd.     

Two-stage aquifer pumping subject to slow desorption and persistent sources

This work focuses on improving pump-and-treat remediation by optimizing a two-stage operational scheme to reduce volumes extracted when confronted with nonequilibrium desorption, low-permeability units, and continuous contaminant sources such as non-aqueous phase liquids (NAPL). Q1 and Q2 are the initial short-term high pumping rate and later long-term low pumping rate, respectively. A two-dimensional ground water flow and transport management model was used to test the proposed strategy for plumes developed from finite (NAPL-free) and continuous (NAPL-driven) contaminant sources in homogeneous and nonhomogeneous (zoned) aquifers. Remediation scenarios were simulated over durations of 2000, 6000, and 15,000 d to determine (1) the optimal time to switch from a preset Q1 to Q2 and (2) the value of Q2. The problem was constrained by mass removal requirements, maximum allowable downgradient concentrations, and practical bounds on Q2. Q1 was fixed at preset values 50% to 200% higher than the single-stage pumping rates (i.e., steady pumping rates during entire remediation period) necessary to achieve a desired cleanup level and capture the plume. Results for the NAPL-free homogeneous case under nonequilibrium desorption conditions achieved the same level of cleanup as single-stage pumping, while reducing extracted volumes by up to 36%. Comparable savings were obtained with NAPL-driven sources only when the source concentration was reduced by at least 2 orders of magnitude. For the zoned aquifer, the proposed strategy provided volume savings of up to 24% under NAPL-free and reduced source conditions.     

Boundary depletion rate and drawdown in leaky wedge‐shaped aquifers

This study presents analytical solutions of the three‐dimensional groundwater flow to a well in leaky confined and leaky water table wedge‐shaped aquifers. Leaky wedge‐shaped aquifers with and without storage in the aquitard are considered, and both transient and steady‐state drawdown solutions are derived. Unlike the previous solutions of the wedge‐shaped aquifers, the leakages from aquitard are considered in these solutions and unlike similar previous work for leaky aquifers, leakage from aquitards and from the water table are treated as the lower and upper boundary conditions. A special form of finite Fourier transforms is used to transform the z‐coordinate in deriving the solutions. The leakage induced by a partially penetrating pumping well in a wedge‐shaped aquifer depends on aquitard hydraulic parameters, the wedge‐shaped aquifer parameters, as well as the pumping well parameters. We calculate lateral boundary dimensionless flux at a representative line and investigate its sensitivity to the aquitard hydraulic parameters. We also investigate the effects of wedge angle, partial penetration, screen location and piezometer location on the steady‐state dimensionless drawdown for different leakage parameters. Results of our study are presented in the form of dimensionless flux‐dimensionless time and dimensionless drawdown‐leakage parameter type curves. The results are useful for evaluating the relative role of lateral wedge boundaries and leakage source on flow in wedge‐shaped aquifers. This is very useful for water management problems and for assessing groundwater pollution. The presented analytical solutions can also be used in parameter identification and in calculating stream depletion rate and volume. Copyright © 2011 John Wiley & Sons, Ltd.