A coupled vertically integrated model to describe lateral exchanges between surface and subsurface in large alluvial floodplains with a fully penetrating river

This paper presents a vertically averaged model for studying water and solute exchanges between a large river and its adjacent alluvial aquifer. The hydraulic model couples horizontal 2D Saint Venant equations for river flow and a 2D Dupuit equation for aquifer flow. The dynamic coupling between river and aquifer is provided by continuity of fluxes and water level elevation between the two domains. Equations are solved simultaneously by linking the two hydrological system matrices in a single global matrix in order to ensure the continuity conditions between river and aquifer and to accurately model two‐way coupling between these two domains. The model is applied to a large reach (about 36 km²) of the Garonne River (south‐western France) and its floodplain, including an instrumented site in a meander. Simulated hydraulic heads are compared with experimental measurements on the Garonne River and aquifer in the floodplain. Model verification includes comparisons for one point sampling date (27 piezometers, 30 March 2000) and for hydraulic heads variations measured continuously over 5 months (5 piezometers, 1 January to 1 June 2000). The model accurately reproduces the strong hydraulic connections between the Garonne River and its aquifer, which are confirmed by the simultaneous variation of the water level in the river and in piezometers located near the river bank. The simulations also confirmed that the model is able to reproduce groundwater flow dynamics during flood events. Given these results, the hydraulic model was coupled with a solute‐transport component, based on advection‐dispersion equations, to investigate the theoretical dynamics of a conservative tracer over 5 years throughout the 36 km² reach studied. Meanders were shown to favour exchanges between river and aquifer, and although the tracer was diluted in the river, the contamination moved downstream from the injection plots and affected both river banks. Copyright © 2008 John Wiley & Sons, Ltd.     

Transfer function‐noise modelling of an aquifer system in NE Italy

In this article, we use a transfer function‐noise (TFN) modelling strategy with single output and multiple/single inputs to study the relationships among a large unconfined aquifer in the upper Venetian plain (Northeast Italy), its recharge components (rainfalls and losing river) and a multi‐layered confined system located in the middle Venetian plain. Model identification from the data raises a range of difficulties when seeking models with consistent physical behaviour, but no information related to the transfer function order and the lags with no zero weights is available. Therefore we use an automatic identification procedure for TFN models. The obtained results suggest that the rainfall component is more important than the river discharge in the unconfined aquifer, and the behaviours of the deep‐confined aquifers are synchronous with that observed in the Badoere area. Copyright © 2010 John Wiley & Sons, Ltd.     

Assessment of the quantitative and qualitative buffer function of an alluvial wetland: hydrological modelling of a large floodplain (Garonne River,France)

For 2 years, water flow‐patterns in the Garonne floodplain of south‐western France were studied in the field and through hydrodynamic modelling (MARTHE Hydrodynamic Software developed by BRGM). Water flow‐paths and the transport of dissolved elements between river and aquifer have been investigated and modelled. In order to quantify the buffer function of the alluvial floodplain, we focused our work on the effect of a major flood on the water flow‐direction, and on nitrate transport. Thus, we showed that the effect of a large flood in the river was rapidly lost with increasing distance from the river. During the observation period, a hydrologically active strip only 300 m wide on either side of the riverbed played a buffering role in absorbing the flood crest. It was also found that meanders favour the exchange between river and alluvial aquifer, shown by the creation of bypasses between the upstream and downstream parts of meanders. This, in turn, contributes to a dilution of nitrates in the phreatic aquifer, which here has higher nitrate content than the surface water; such dilution may result in an overestimation of the denitrification process in the wooded riverbanks. The coupling of chemical measurements—especially of chlorides and nitrate—with modelling of the dissolved‐element transport allows us to establish the water balance for the riparian wetland, and to separate the effect of dilution and denitrification on nitrate concentration. This indicated the existence of areas in the riparian wetlands where denitrification is particularly strong, leading to reductions in nitrate concentrations of 10 to 30 mg/l NO₃^? during the flood. Copyright © 2003 John Wiley & Sons, Ltd.     

Reducing monitoring gaps at the aquifer–river interface by modelling groundwater–surface water exchange flow patterns

This study investigates spatial patterns and temporal dynamics of aquifer–river exchange flow at a reach of the River Leith, UK. Observations of sub‐channel vertical hydraulic gradients at the field site indicate the dominance of groundwater up‐welling into the river and the absence of groundwater recharge from surface water. However, observed hydraulic heads do not provide information on potential surface water infiltration into the top 0–15 cm of the streambed as these depths are not covered by the existing experimental infrastructure. In order to evaluate whether surface water infiltration is likely to occur outside the ‘window of detection’, i.e. the shallow streambed, a numerical groundwater model is used to simulate hydrological exchanges between the aquifer and the river. Transient simulations of the successfully validated model (Nash and Sutcliff efficiency of 0·91) suggest that surface water infiltration is marginal and that the possibility of significant volumes of surface water infiltrating into non‐monitored shallow streambed sediments can be excluded for the simulation period. Furthermore, the simulation results show that with increasing head differences between river and aquifer towards the end of the simulation period, the impact of streambed topography and hydraulic conductivity on spatial patterns of exchange flow rates decreases. A set of peak flow scenarios with altered groundwater‐surface water head gradients is simulated in order to quantify the potential for surface water infiltration during characteristic winter flow conditions following the observation period. The results indicate that, particularly at the beginning of peak flow conditions, head gradients are likely to cause substantial increase in surface water infiltration into the streambed. The study highlights the potential for the improvement of process understanding of hyporheic exchange flow patterns at the stream reach scale by simulating aquifer‐river exchange fluxes with a standard numerical groundwater model and a simple but robust model structure and parameterization. Copyright © 2011 John Wiley & Sons, Ltd.     

A geomorphology‐based integrated stream–aquifer interaction model for semi‐gauged catchments

Many of the existing stream–aquifer interaction models available in the literature are very complex with limited applicability in semi‐gauged and ungauged catchments. In this study, to estimate the influent and effluent subsurface water fluxes under limited geo‐hydrometeorological data availability conditions, a simple stream–aquifer interaction model, namely, the variable parameter McCarthy–Muskingum (VPMM) hillslope‐storage Boussinesq (hsB) model, has been developed. This novel model couples the VPMM streamflow transport with the hsB groundwater flow transport modules in online mode. In this integrated model, the surface water–groundwater flux exchange process is modelled by the Darcian approach with the variable hydraulic heads between the river stage and groundwater table accounting for the rainfall forcing. Considering the exchange fluxes in the hyporheic zone and lateral overland flow contribution, this approach is field tested in a typical 48‐km stretch of the Brahmani River in eastern India to simulate the streamflow and its depth with the minimum Nash–Sutcliffe efficiency of 94% and 88%; the maximum root mean square error of 134 m³/s and 0.35 m; and the minimum index of agreement of 98% and 97%, respectively. This modelling approach could be very well utilized in data‐scarce world‐river basins to estimate the stream–aquifer exchange flux due to rainfall forcings.     

Linking hydrology and sediment dynamics of large alluvial rivers to landscape diversity in the Ganga dispersal system,India

Discharge and sediment load data for several stations along the Ganga River and its major tributaries in the western Ganga plains (WGP) for a period of ~30 years have been analysed to understand the hydrological characteristics and sediment dynamics. In terms of hydrology, the rivers are less flood‐prone than believed, exceeding bankfull discharges less frequently than the expected 1.5 year return interval. This has been attributed to the rivers of this region occupying incised valleys formed in the Late Quaternary period. Rivers draining the WGP are supply‐limited systems compared to those draining the eastern Ganga plains (EGP) which have been characterized as transport‐limited systems. We suggest that such geomorphic diversity as a function of spatial variability in precipitation regime and hinterland geology has existed for at least the Late Quaternary period and they in turn influence the modern day hydrology of the river systems in a significant way. Copyright © 2016 John Wiley & Sons, Ltd.     

Groundwater quality in the semi‐arid region of the Chahardouly basin,West Iran

Chahardouly basin is located in the western part of Iran and is characterized by semi‐arid climatic conditions and scarcity in water resources. The main aquifer systems are developed within alluvial deposits. The availability of groundwater is rather erratic owing to the occurrence of hard rock formation and a saline zone in some parts of the area. The aquifer systems of the area show signs of depletion, which have taken place in recent years due to a decline in water levels. Groundwater samples collected from shallow and deep wells were analysed to examine the quality characteristics of groundwater. The major ion chemistry of groundwater is dominated by Ca²⁺ and HCO₃^?, while higher values of total dissolved solids (TDS) in groundwater are associated with high concentrations of all major ions. An increase in salinity is recorded in the down‐gradient part of the basin. The occurrence of saline groundwater, as witnessed by the high electrical conductivity (EC), may be attributed to the long residence time of water and the dissolution of minerals, as well as evaporation of rainfall and irrigation return flow. Based on SAR values and sodium content (%Na), salinity appears to be responsible for the poor groundwater quality, rendering most of the samples not suitable for irrigation use. Copyright © 2007 John Wiley & Sons, Ltd.     

Modelling the water level of the alluvial aquifer of an ephemeral river in south-western Zimbabwe

ABSTRACT

Water from the alluvium of ephemeral rivers in Zimbabwe is increasingly being used. These alluvial aquifers are recharged annually from infiltrating floodwater. Nonetheless, the size of this water resource is not without limit and an understanding of the hydrological processes of an alluvial aquifer is required for its sustainable management. This paper presents the development of a water balance model, which estimates the water level in an alluvial aquifer recharged by surface flow and rainfall, while allowing for abstraction, evaporation and other losses. The model is coupled with a watershed model, which generates inflows from upland catchment areas and tributaries. Climate, hydrological, land cover and geomorphological data were collected as inputs to both models as well as observed flow and water levels for model calibration and validation. The sand river model was found to be good at simulating the observed water level and was most sensitive to porosity and seepage.     

Hydrogeological interactions between fault zones and alluvial aquifers in regional flow systems

The hydrological influence of fault zones in tectonic areas is usually difficult to depict from field data. Numerical simulation allows representation of such flow systems and an estimation of flow lines and rates. This paper reports on simulations of the groundwater flow in a range‐and‐basin area affected by a regional fault zone, which may drain or recharge an overlaying alluvial aquifer. Different hydraulic conductivity values for the range rocks, the fault‐zone, and the sedimentary infill of the basin are considered, as well as different fault‐zone widths and boundary conditions. Results show that upward and downward fluxes develop in the upper part of the fault zone, controlled by the action of the alluvial aquifer, influencing the recharge of the sedimentary basin. This paper shows the hydrological efficiency of fault zones as preferential flow; it also analyses the constraints that determine groundwater recharge to the surrounding basins. These results contribute to the understanding of hydrogeological dynamics in tectonic areas. Copyright © 2008 John Wiley & Sons, Ltd.     

Dissolution in a variably confined carbonate platform: effects of allogenic runoff,hydraulic damming of groundwater inputs,and surface–groundwater exchange at the basin scale

In variably confined carbonate platforms, impermeable confining units collect rainfall over large areas and deliver runoff to rivers or conduits in unconfined portions of platforms. Runoff can increase river stage or conduit heads in unconfined portions of platforms faster than local infiltration of rainfall can increase groundwater heads, causing hydraulic gradients between rivers, conduits and the aquifer to reverse. Gradient reversals cause flood waters to flow from rivers and conduits into the aquifer where they can dissolve limestone. Previous work on impacts of gradient reversals on dissolution has primarily emphasized individual caves and little research has been conducted at basin scales. To address this gap in knowledge, we used legacy data to assess how a gradient of aquifer confinement across the Suwannee River Basin, north‐central Florida affected locations, magnitudes and processes of dissolution during 2005–2007, a period with extreme ranges of discharge. During intense rain events, runoff from the confining unit increased river stage above groundwater heads in unconfined portions of the platform, hydraulically damming inputs of groundwater along a 200 km reach of river. Hydraulic damming allowed allogenic runoff with SI_CAL < ?4 to fill the entire river channel and flow into the aquifer via reversing springs. Storage of runoff in the aquifer decreased peak river discharges downstream and contributed to dissolution within the aquifer. Temporary storage of allogenic runoff in karst aquifers represents hyporheic exchange at a scale that is larger than found in streams flowing over non‐karst aquifers because conduits in karst aquifers extend the area available for exchange beyond river beds deep into aquifers. Post‐depositional porosity in variably confined carbonate platforms should thus be enhanced along rivers that originate on confining units. This distribution should be considered in models of porosity distribution used to manage water and hydrocarbon resources in carbonate rocks. Copyright © 2013 John Wiley & Sons, Ltd.     

Geomorphology,hydrology, and aquatic vegetation drive seasonal hyporheic flow patterns across a gravel‐dominated floodplain

Across 1·7 km² of the Umatilla River floodplain (Oregon, USA), we investigated the influences of an ephemeral tributary and perennial ‘spring channel’ (fed only by upwelling groundwater) on hyporheic hydrology. We derived maps of winter and summer water‐table elevations from data collected at 46 monitoring wells and 19 stage gauges and used resulting maps to infer groundwater flow direction. Groundwater flow direction varied seasonally across the floodplain and was influenced by main channel stage, flooding, the tributary creek, and the location and direction of hyporheic exchange in the spring channel. Hyporheic exchange in the spring channel was evaluated with a geochemical mixing model, which confirmed patterns of floodplain groundwater movement inferred from water‐table maps and showed that the spring channel was fed predominantly by hyporheic water from the floodplain aquifer (87% during winter, 80% during summer), with its remaining flow supplied by upslope groundwater from the adjacent catchment aquifer. Summertime growth of aquatic macrophytes in the spring channel also influenced patterns of hyporheic exchange and groundwater flow direction in the alluvial aquifer by increasing flow resistance in the spring channel, locally raising surface water stage and adjacent water‐table elevation, and thereby altering the slope of the water‐table in the hyporheic zone. The Umatilla River floodplain is larger than most sites where hyporheic hydrology has been investigated in detail. Yet, our results corroborate other research that has identified off‐channel geomorphic features as important drivers of hyporheic hydrology, including previously published modeling efforts from a similar river and field observations from smaller streams. Copyright © 2007 John Wiley & Sons, Ltd.     

Radial groundwater flow to a finite diameter well in a leaky confined aquifer with a finite‐thickness skin

A mathematical model that describes the drawdown due to constant pumpage from a finite radius well in a two‐zone leaky confined aquifer system is presented. The aquifer system is overlain by an aquitard and underlain by an impermeable formation. A skin zone of constant thickness exists around the wellbore. A general solution to a two‐zone leaky confined aquifer system in Laplace domain is developed and inverted numerically to the time‐domain solution using the modified Crump (1976) algorithm. The results show that the drawdown distribution is significantly influenced by the properties and thickness of the skin zone and aquitard. The sensitivity analyses of parameters of the aquifer and aquitard are performed to illustrate their effects on drawdowns in a two‐zone leaky confined aquifer system. For the negative‐skin case, the drawdown is very sensitive to the relative change in the formation transmissivity. For the positive‐skin case, the drawdown is also sensitive to the relative changes in the skin thickness, and both the skin and formation transmissivities over the entire pumping period and the well radius and formation storage coefficient at early pumping time. Copyright © 2009 John Wiley & Sons, Ltd.     

A patch hierarchy approach to modeling surface and subsurface hydrology in complex flood‐plain environments

Geomorphology interacts with surface‐ and ground‐water hydrology across multiple spatial scales. Nonetheless, hydrologic and hydrogeologic models are most commonly implemented at a single spatial scale. Using an existing hydrogeologic computer model, we implemented a simple hierarchical approach to modeling surface‐ and ground‐water hydrology in a complex geomorphic setting. We parameterized the model to simulate ground‐ and surface‐water ?ow patterns through a hierarchical, three‐dimensional, quantitative representation of an anabranched montane alluvial ?ood plain (the Nyack Flood Plain, Middle Fork Flathead River, Montana, USA). Comparison of model results to ?eld data showed that the model provided reasonable representations of spatial patterns of aquifer recharge and discharge, temporal patterns of ?ood‐water storage on the ?ood plain, and rates of ground‐water movement from the main river channel into a large lateral spring channel on the ?ood plain, and water table elevation in the alluvial aquifer. These results suggest that a hierarchical approach to modeling ground‐ and surface‐water hydrology can reproduce realistic patterns of surface‐ and ground‐water ?ux on alluvial ?ood plains, and therefore should provide an excellent ‘quantitative laboratory’ for studying complex interactions between geomorphology and hydrology at and across multiple spatial scales. Copyright © 2004 John Wiley & Sons, Ltd.     

Assimilation of historical head data to estimate spatial distributions of stream bed and aquifer hydraulic conductivity fields

Management of water resources in alluvial aquifers relies mainly on understanding interactions between hydraulically connected streams and aquifers. Numerical models that simulate this interaction often are used as decision support tools for water resource management. However, the accuracy of numerical predictions relies heavily on unknown system parameters (e.g., streambed conductivity and aquifer hydraulic conductivity), which are spatially heterogeneous and difficult to measure directly. This paper employs an ensemble smoother to invert groundwater level measurements to jointly estimate spatially varying streambed and alluvial aquifer hydraulic conductivity along a 35.6‐km segment of the South Platte River in Northeastern Colorado. The accuracy of the inversion procedure is evaluated using a synthetic experiment and historical groundwater level measurements, with the latter constituting the novelty of this study in the inversion and validation of high‐resolution fields of streambed and aquifer conductivities. Results show that the estimated streambed conductivity field and aquifer conductivity field produce an acceptable agreement between observed and simulated groundwater levels and stream flow rates. The estimated parameter fields are also used to simulate the spatially varying flow exchange between the alluvial aquifer and the stream, which exhibits high spatial variability along the river reach with a maximum average monthly aquifer gain of about 2.3 m³/day and a maximum average monthly aquifer loss of 2.8 m³/day, per unit area of streambed (m²). These results demonstrate that data assimilation inversion provides a reliable and computationally affordable tool to estimate the spatial variability of streambed and aquifer conductivities at high resolution in real‐world systems.     

Using chlorofluorocarbons (CFCs) and tritium (3H) to estimate groundwater age and flow velocity in Hohhot Basin,China

The concentrations of chlorofluorocarbons (CFC‐11, CFC‐12 and CFC‐113) and tritium (³H) content in groundwater were used to date groundwater age, delineate groundwater flow systems and estimate flow velocity in the Hohhot basin. The estimated young groundwater age is fallen in the bracket of 21 ~ 50 a and indicates the presence of two different age profiles and flow systems in the shallow groundwater system. Older age waters occur under the topographically low areas, where the aquifer is double‐layer aquifer system consisting of shallow unconfined‐semi‐confined aquifer and deep confined aquifer. This reflects long flow paths associated with regional flow. Groundwater (range from 21 to 34 years) in the north piedmont and east hilly areas, where the aquifer is a single‐layer aquifer consisting of alluvial fans, are typically younger than those in the low areas. The combination of CFCs dating with hydrogeological information indicates that both local and regional flow systems are present at the basin. The regional groundwater flow mainly flows from the north and east to the southwest, the local groundwater flow system occurs nearby the Hohhot city. The mean regional groundwater flow velocity of the shallow groundwater is estimated about 0.73 km/a. These findings can aid in refining hydrogeological conceptual model of the study area. Copyright © 2012 John Wiley & Sons, Ltd.     

Influence of partial confinement and Holocene river formation on groundwater flow and dissolution in the Florida carbonate platform

Much of what is known about groundwater circulation and geochemical evolution in carbonate platforms is based on platforms that are fully confined or unconfined. Much less is known about groundwater flow paths and geochemical evolution in partially confined platforms, particularly those supporting surface water. In north‐central Florida, sea level rise and a transition to a wetter climate during the Holocene formed rivers in unconfined portions of the Florida carbonate platform. Focusing on data from the Santa Fe River basin, we show river formation has led to important differences in the hydrological and geochemical evolution of the Santa Fe River basin relative to fully confined or unconfined platforms. Runoff from the siliciclastic confining layer drove river incision and created topographic relief, reorienting the termination of local and regional groundwater flow paths from the coast to the rivers in unconfined portions of the platform. The most chemically evolved groundwater occurs at the end of the longest and deepest flow paths, which discharge near the center of the platform because of incision of the Santa Fe River at the edge of the confining unit. This pattern of discharge of mineralized water differs from fully confined or unconfined platforms where discharge of the most mineralized water occurs at the coast. Mineralized water flowing into the Santa Fe River is diluted by less evolved water derived from shorter, shallower flow paths that discharge to the river downstream. Formation of rivers shortens flow path lengths, thereby decreasing groundwater residence times and allowing freshwater to discharge more quickly to the oceans in the newly formed rivers than in platforms that lack rivers. Similar dynamic changes to groundwater systems should be expected to occur in the future as climate change and sea level rise develop surface water on other carbonate platforms and low lying coastal aquifer systems. Copyright © 2012 John Wiley & Sons, Ltd.     

Application of an integrated surface water‐groundwater model to multi‐aquifers modeling in Choushui River alluvial fan,Taiwan

This research proposes a combination of SWAT and MODFLOW, MD‐SWAT‐MODFLOW, to address the multi‐aquifers condition in Choushui River alluvial fan, Taiwan. The natural recharge and unidentified pumping/recharge are separately estimated. The model identifies the monthly pumping/recharge rates in multi‐aquifers so that the daily streamflow can be simulated correctly. A multi‐aquifers condition means a subsurface formation composed of at least the unconfined aquifer, the confined aquifer, and an in‐between aquitard. In such a case, the variation of groundwater level is related to pumping/recharge activities in vertically adjacent aquifer and the river‐aquifer interaction. Both factors in turn affect the streamflow performance. Results show that MD‐SWAT‐MODFLOW performs better than SWAT alone in terms of simulated streamflow, especially during low flow period, when pumping/recharge rates are properly estimated. A sensitivity analysis of individual parameter suggests that the vertical leakance may be the most sensitive among all investigated MODFLOW parameters in terms of the estimated pumping/recharge among aquifers, and the Latin‐Hypercube‐One‐factor‐At‐a‐Time sensitivity analysis indicates that the hydraulic conductivity of channel is the most sensitive to the model performance. It also points out the necessity to simultaneously estimate pumping/recharge rates in multi‐aquifers. The estimated net pumping rate can be treated as a lower bound of the actual local pumping rate. As a whole, the model provides the spatio‐temporal groundwater use, which gives the authorities insights to manage groundwater resources. Copyright © 2012 John Wiley & Sons, Ltd.     

A test of equilibrium theory and a demonstration of its practical application for predicting the morphodynamics of the Yangtze River

Taking the width/depth ratio of a river channel as an independent variable, a variational analysis of basic flow relationships shows that alluvial‐channel flow adjusts channel geometry to achieve stationary equilibrium when the condition of maximum flow efficiency (MFE) is satisfied. As a test of the veracity of MFE and to examine if this theory of self‐adjusting channel morphodynamics can be practically applied to large river systems, this study examines the degree of correspondence between theoretically determined equilibrium channel geometries and actual measurements along the middle and lower Yangtze River. Using four different forms of the Meyer‐Peter and Müller bedload relation and relations of flow continuity and resistance we show that the Meyer‐Peter and Müller bedload relation modified on the basis of MFE theory predicts channel dimensions most accurately when applied to the middle and lower Yangtze River. This provides convincing evidence supporting MFE equilibrium theory. Copyright © 2014 John Wiley & Sons, Ltd.     

A two‐dimensional analytical solution for groundwater flow in a leaky confined aquifer system near open tidal water

Groundwater in coastal areas is commonly disturbed by tidal fluctuations. A two‐dimensional analytical solution is derived to describe the groundwater fluctuation in a leaky confined aquifer system near open tidal water under the assumption that the groundwater head in the confined aquifer fluctuates in response to sea tide whereas that of the overlying unconfined aquifer remains constant. The analytical solution presented here is an extension of the solution by Sun for two‐dimensional groundwater flow in a confined aquifer and the solution by Jiao and Tang for one‐dimensional groundwater flow in a leaky confined aquifer. The analytical solution is compared with a two‐dimensional finite difference solution. On the basis of the analytical solution, the groundwater head distribution in a leaky confined aquifer in response to tidal boundaries is examined and the influence of leakage on groundwater fluctuation is discussed. Copyright © 2001 John Wiley & Sons, Ltd.     

Groundwater response to tidal fluctuations in a leaky confined coastal aquifer with a finite length

The study on the hydraulic properties of coastal aquifers has significant implications both in hydrological sciences and environmental engineering. Although many analytical solutions are available, most of them are based on the same basic assumption that assumes aquifers extend landward semi‐infinitely, which does not necessarily reflect the reality. In this study, the general solutions for a leaky confined coastal aquifer have been developed that consider both finitely landward constant‐head and no‐flow boundaries. The newly developed solutions were then used to examine theoretically the joint effects of leakage and aquifer length on hydraulic head fluctuations within the leaky confined aquifer, and the validity of using the simplified solution, which assumes the aquifer is semi‐infinite. The results illustrated that the use of the simplified solution may cause significant errors, depending on joint effects of leakage and aquifer length. A dimensionless characteristic parameter was then proposed as an index for judging the applicability of the simplified solution. In addition, practical application of the general solution for the constant‐head inland boundary was used to characterize the hydraulic properties of a leaky confined aquifer using the data collected from a field site at the Seine River estuary, France, and the versatility of the general solution was further justified.