Improved simulation of river water and groundwater exchange in an alluvial plain using the SWAT model

Hydrological interaction between surface and subsurface water systems has a significant impact on water quality, ecosystems and biogeochemistry cycling of both systems. Distributed models have been developed to simulate this function, but they require detailed spatial inputs and extensive computation time. The soil and water assessment tool (SWAT) model is a semi‐distributed model that has been successfully applied around the world. However, it has not been able to simulate the two‐way exchanges between surface water and groundwater. In this study, the SWAT‐landscape unit (LU) model – based on a catena method that routes flow across three LUs (the divide, the hillslope and the valley) – was modified and applied in the floodplain of the Garonne River. The modified model was called SWAT‐LUD. Darcy's equation was applied to simulate groundwater flow. The algorithm for surface water‐level simulation during flooding periods was modified, and the influence of flooding on groundwater levels was added to the model. Chloride was chosen as a conservative tracer to test simulated water exchanges. The simulated water exchange quantity from SWAT‐LUD was compared with the output of a two‐dimensional distributed model, surface–subsurface water exchange model. The results showed that simulated groundwater levels in the LU adjoining the river matched the observed data very well. Additionally, SWAT‐LUD model was able to reflect the actual water exchange between the river and the aquifer. It showed that river water discharge has a significant influence on the surface–groundwater exchanges. The main water flow direction in the river/groundwater interface was from groundwater to river; water that flowed in this direction accounted for 65% of the total exchanged water volume. The water mixing occurs mainly during high hydraulic periods. Flooded water was important for the surface–subsurface water exchange process; it accounted for 69% of total water that flowed from the river to the aquifer. The new module also provides the option of simulating pollution transfer occurring at the river/groundwater interface at the catchment scale. Copyright © 2015 John Wiley & Sons, Ltd.     

Simulated groundwater interaction with rivers and springs in the Heihe river basin

Surface water and groundwater in the Heihe river basin of China are interconnected and the pattern of water resources exploitation has a direct effect on the interaction of groundwater and surface water, especially on a downstream oasis. A three‐dimensional groundwater flow simulation model with eight model layers was established to simulate the regional groundwater flow in the multilayered aquifer system and the interaction among the rivers, springs, and groundwater. The model was calibrated not only with historical water levels but also with the investigated baseflow and spring flux. The simulation results of the numerical model match reasonably well with the observed groundwater levels, baseflow to rivers, and spring flux. The numerical simulation also demonstrates that the hydraulic connection between the river and the aquifers has transferred from the coupling to decoupling at some reaches. It is suggested that there is a vital need to reduce groundwater withdrawal and to rationalize the use of both groundwater and surface water in order to maintain sustainable development in the study area. Copyright © 2007 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.     

Is there a superior conceptual groundwater model structure for baseflow simulation?

Previous work has shown that streamflow response during baseflow conditions is a function of storage, but also that this functional relationship varies among seasons and catchments. Traditionally, hydrological models incorporate conceptual groundwater models consisting of linear or non‐linear storage–outflow functions. Identification of the right model structure and model parameterization however is challenging. The aim of this paper is to systematically test different model structures in a set of catchments where different aquifer types govern baseflow generation processes. Nine different two‐parameter conceptual groundwater models are applied with multi‐objective calibration to transform two different groundwater recharge series derived from a soil‐atmosphere‐vegetation transfer model into baseflow separated from streamflow data. The relative performance differences of the model structures allow to systematically improve the understanding of baseflow generation processes and to identify most appropriate model structures for different aquifer types. We found more versatile and more aquifer‐specific optimal model structures and elucidate the role of interflow, flow paths, recharge regimes and partially contributing storages. Aquifer‐specific recommendations of storage models were found for fractured and karstic aquifers, whereas large storage capacities blur the identification of superior model structures for complex and porous aquifers. A model performance matrix is presented, which highlights the joint effects of different recharge inputs, calibration criteria, model structures and aquifer types. The matrix is a guidance to improve groundwater model structures towards their representation of the dominant baseflow generation processes of specific aquifer types. Copyright © 2014 John Wiley & Sons, Ltd.     

Implementing a nonlinear groundwater module in the soil and water assessment tool (SWAT)

The Soil and Water Assessment Tool (SWAT) is widely used in modeling water quantity and quality. In the original SWAT, groundwater flow is calculated using a linear‐reservoir model, with outflow proportional to storage. However, observations show that this assumption is not always applicable; for example, macropores in Karst formations would seriously affect the groundwater behavior. A nonlinear groundwater algorithm was introduced in a new version of the SWAT model, called ISWAT. The Shenandoah Valley area in the Eastern U.S., which includes a number of geologic formations including Karst, was selected to test the modified ISWAT model. Parameter ESTimation (PEST) was coupled with ISWAT to auto‐calibrate the nonlinear parameter values. Ten years of record at 15 stream gauges were used to calibrate the model. The nonlinear ISWAT, statistically and visually, performed better in stream discharge estimation especially during baseflow recession and low‐flow periods. This indicated that the nonlinear algorithm can better represent groundwater behavior. The coupled ISWAT‐PEST approach can be used in future stream discharge simulation. Copyright © 2013 John Wiley & Sons, Ltd.     

A multi‐storage groundwater concept for the SWAT model to emphasize nonlinear groundwater dynamics in lowland catchments

Hydrological models are useful tools to analyze present and future conditions of water quantity and quality. The integrated modelling of water and nutrients needs an adequate representation of the different discharge components. In common with many lowlands, groundwater contribution to the discharge in the North German lowlands is a key factor for a reasonable representation of the water balance, especially in low flow periods. Several studies revealed that the widely used Soil and Water Assessment Tool (SWAT) model performs poorly for low flow periods. This paper deals with the extension of the groundwater module of the SWAT model to enhance low flow representation. The current two‐storage concept of SWAT was further developed to a three‐storage concept. This was realized due to modification of the groundwater module by splitting the active groundwater storage into a fast and a slow contributing aquifer. The results of this study show that the groundwater module with three storages leads to a good prediction of the overall discharge especially for the recession limbs and the low flow periods. The improved performance is reflected in the signature measures for the mid‐segment (percent bias ?2.4% vs ?15.9%) and the low segment (percent bias 14.8% vs 46.8%) of the flow duration curve. The three‐storage groundwater module is more process oriented than the original version due to the introduction of a fast and a slow groundwater flow component. The three‐storage version includes a modular approach, because groundwater storages can be activated or deactivated independently for subbasin and hydrological response unit level. Copyright © 2013 John Wiley & Sons, Ltd.     

Predicting saturation‐excess runoff distribution with a lumped hillslope model: SWAT‐HS

Saturation‐excess runoff is the major runoff mechanism in humid well‐vegetated areas where infiltration rates often exceed rainfall intensity. Although the Soil and Water Assessment Tool (SWAT) is one of the most widely used models, it predicts runoff based mainly on soil and land use characteristics, and is implicitly an infiltration‐excess runoff type of model. Previous attempts to incorporate the saturation‐excess runoff mechanism in SWAT fell short due to the inability to distribute water from one hydrological response unit to another. This paper introduces a modified version of SWAT, referred to as SWAT‐Hillslope (SWAT‐HS). This modification improves the simulation of saturation‐excess runoff by redefining hydrological response units based on wetness classes and by introducing a surface aquifer with the ability to route interflow from “drier” to “wetter” wetness classes. Mathematically, the surface aquifer is a nonlinear reservoir that generates rapid subsurface stormflow as the water table in the surface aquifer rises. The SWAT‐HS model was tested in the Town Brook watershed in the upper reaches of the West Branch Delaware River in the Catskill region of New York, USA. SWAT‐HS predicted discharge well with a Nash‐Sutcliffe Efficiency of 0.68 and 0.87 for daily and monthly time steps. Compared to the original SWAT model, SWAT‐HS predicted less surface runoff and groundwater flow and more lateral flow. The saturated areas predicted by SWAT‐HS were concentrated in locations with a high topographic index and were in agreement with field observations. With the incorporation of topographic characteristics and the addition of the surface aquifer, SWAT‐HS improved streamflow simulation and gave a good representation of saturated areas on the dates that measurements were available. SWAT‐HS is expected to improve water quality model predictions where the location of the surface runoff matters.     

The Effect of Modeled Recharge Distribution on Simulated Groundwater Availability and Capture

Simulating groundwater flow in basin‐fill aquifers of the semiarid southwestern United States commonly requires decisions about how to distribute aquifer recharge. Precipitation can recharge basin‐fill aquifers by direct infiltration and transport through faults and fractures in the high‐elevation areas, by flowing overland through high‐elevation areas to infiltrate at basin‐fill margins along mountain fronts, by flowing overland to infiltrate along ephemeral channels that often traverse basins in the area, or by some combination of these processes. The importance of accurately simulating recharge distributions is a current topic of discussion among hydrologists and water managers in the region, but no comparative study has been performed to analyze the effects of different recharge distributions on groundwater simulations. This study investigates the importance of the distribution of aquifer recharge in simulating regional groundwater flow in basin‐fill aquifers by calibrating a groundwater‐flow model to four different recharge distributions, all with the same total amount of recharge. Similarities are seen in results from steady‐state models for optimized hydraulic conductivity values, fit of simulated to observed hydraulic heads, and composite scaled sensitivities of conductivity parameter zones. Transient simulations with hypothetical storage properties and pumping rates produce similar capture rates and storage change results, but differences are noted in the rate of drawdown at some well locations owing to the differences in optimized hydraulic conductivity. Depending on whether the purpose of the groundwater model is to simulate changes in groundwater levels or changes in storage and capture, the distribution of aquifer recharge may or may not be of primary importance.     

How to improve the representation of hydrological processes in SWAT for a lowland catchment – temporal analysis of parameter sensitivity and model performance

Model diagnostic analyses help to improve the understanding of hydrological processes and their representation in hydrological models. A detailed temporal analysis detects periods of poor model performance and model components with potential for model improvements, which cannot be found by analysing the whole discharge time series. In this study, we aim to improve the understanding of hydrological processes by investigating the temporal dynamics of parameter sensitivity and of model performance for the Soil and Water Assessment Tool model applied to the Treene lowland catchment in Northern Germany. The temporal analysis shows that the parameter sensitivity varies temporally with high sensitivity for three groundwater parameters (groundwater time delay, baseflow recession constant and aquifer fraction coefficient) and one evaporation parameter (soil evaporation compensation factor). Whereas the soil evaporation compensation factor dominates in baseflow and resaturation periods, groundwater time delay, baseflow recession constant and aquifer fraction coefficient are dominant in the peak and recession phases. The temporal analysis of model performance identifies three clusters with different model performances, which can be related to different phases of the hydrograph. The lowest performance, when comparing six performance measures, is detected for the baseflow cluster. A spatially distributed analysis for six hydrological stations within the Treene catchment shows similar results for all stations. The linkage of periods with poor model performance to the dominant model components in these phases and with the related hydrological processes shows that the groundwater module has the highest potential for improvement. This temporal diagnostic analysis enhances the understanding of the Soil and Water Assessment Tool model and of the dominant hydrological processes in the lowland catchment. Copyright © 2013 John Wiley & Sons, Ltd.     

The effects of climate changes on aquifer storage and river baseflow

Abstract

The effects of changes in climate on aquifer storage and groundwater flow to rivers have been investigated using an idealized representation of the aquifer/river system. The generalized aquifer/river model can incorporate spatial variability in aquifer transmissivity and is applied with parameters characteristic of Chalk and Triassic sandstone aquifers in the United Kingdom, and is also applicable to other aquifers elsewhere. The model is run using historical time series of recharge, estimated from observed rainfall and potential evaporation data, and with climate inputs perturbed according to a number of climate change scenarios. Simulations of baseflow suggest large proportional reductions at low flows from Chalk under high evaporation change scenarios. Simulated baseflow from the slower responding Triassic sandstone aquifer shows more uniform and less severe reductions. The change in hydrological regime is less extreme for the low evaporation change scenario, but remains significant for the Chalk aquifer.     

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.     

Modeling Steady Sea Water Intrusion with Single‐Density Groundwater Codes

Steady interface flow in heterogeneous aquifer systems is simulated with single‐density groundwater codes by using transformed values for the hydraulic conductivity and thickness of the aquifers and aquitards. For example, unconfined interface flow may be simulated with a transformed model by setting the base of the aquifer to sea level and by multiplying the hydraulic conductivity with 41 (for sea water density of 1025 kg/m³). Similar transformations are derived for unconfined interface flow with a finite aquifer base and for confined multi‐aquifer interface flow. The head and flow distribution are identical in the transformed and original model domains. The location of the interface is obtained through application of the Ghyben‐Herzberg formula. The transformed problem may be solved with a single‐density code that is able to simulate unconfined flow where the saturated thickness is a linear function of the head and, depending on the boundary conditions, the code needs to be able to simulate dry cells where the saturated thickness is zero. For multi‐aquifer interface flow, an additional requirement is that the code must be able to handle vertical leakage in situations where flow in an aquifer is unconfined while there is also flow in the aquifer directly above it. Specific examples and limitations are discussed for the application of the approach with MODFLOW. Comparisons between exact interface flow solutions and MODFLOW solutions of the transformed model domain show good agreement. The presented approach is an efficient alternative to running transient sea water intrusion models until steady state is reached.     

Focused Groundwater Controlled Feedbacks into the Hyporheic Zone During Baseflow Recession

Groundwater surface water interaction in the hyporheic zone remains an important challenge for water resources management and ecosystem restoration. In heterogeneous stratified glacial sediments, reach‐scale environments contain an uneven distribution of focused groundwater flow occurring simultaneously with diffusely discharging groundwater. This results in a variation of stream‐aquifer interactions, where focused flow systems are able to temporally dominate exchange processes. The research presented here investigates the direct and indirect influences focused groundwater discharge exerts on the hyporheic zone during baseflow recession. Field results demonstrate that as diffuse sources of groundwater deplete during baseflow recession, focused groundwater discharge remains constant. During baseflow recession the hyporheic zone is unable to expand, while the high nitrate concentration from focused discharge changes the chemistry of the stream. The final result is a higher concentration of nitrate in the hyporheic zone as this altered surface water infiltrates into the subsurface. This indirect coupling of focused groundwater discharge and the hyporheic zone is unaccounted for in hyporheic studies at this time. Results indicate important implications for the potential reduction of agricultural degradation of water quality.     

A low‐dimensional hillslope‐based catchment model for layered groundwater flow

Despite the strong interaction between surface and subsurface waters, groundwater flow representation is often oversimplified in hydrological models. For instance, the interplay between local or shallow aquifers and deeper regional‐scale aquifers is typically neglected. In this work, a novel hillslope‐based catchment model for the simulation of combined shallow and deep groundwater flow is presented. The model consists of the hillslope‐storage Boussinesq (hsB) model representing shallow groundwater flow and an analytic element (AE) model representing deep regional groundwater flow. The component models are iteratively coupled via a leakage term based on Darcy's law, representing delayed recharge to the regional aquifer through a low conductivity layer. Simulations on synthetic single hillslopes and on a two‐hillslope open‐book catchment are presented, and the results are compared against a benchmark three‐dimensional Richards equation model. The impact of hydraulic conductivity, hillslope plan geometry (uniform, convergent, divergent), and hillslope inclination (0.2%, 5%, and 30%) under drainage and recharge conditions are examined. On the single hillslopes, good matches for heads, hydrographs, and exchange fluxes are generally obtained, with the most significant differences in outflows and heads observed for the 30% slope and for hillslopes with convergent geometry. On the open‐book catchment, cumulative outflows are overestimated by 1–4%. Heads in the confined and unconfined aquifers are adequately reproduced throughout the catchment, whereas exchange fluxes are found to be very sensitive to the hillslope drainable porosity. The new model is highly efficient computationally compared to the benchmark model. The coupled hsB/AE model represents an alternative to commonly used groundwater flow representations in hydrological models, of particular appeal when surface–subsurface exchanges, local aquifer–regional aquifer interactions, and low flows play a key role in a watershed's dynamics. Copyright © 2011 John Wiley & Sons, Ltd.     

Semi‐Arid Aquifer Responses to Forest Restoration Treatments and Climate Change

The purpose of this study was to develop an interpretive groundwater‐flow model to assess the impacts that planned forest restoration treatments and anticipated climate change will have on large regional, deep (>400 m), semi‐arid aquifers. Simulations were conducted to examine how tree basal area reductions impact groundwater recharge from historic conditions to 2099. Novel spatial analyses were conducted to determine areas and rates of potential increases in groundwater recharge. Changes in recharge were applied to the model by identifying zones of basal area reduction from planned forest restoration treatments and applying recharge‐change factors to these zones. Over a 10‐year period of forest restoration treatment, a 2.8% increase in recharge to one adjacent groundwater basin (the Verde Valley sub‐basin) was estimated, compared to conditions that existed from 2000 to 2005. However, this increase in recharge was assumed to quickly decline after treatment due to regrowth of vegetation and forest underbrush and their associated increased evapotranspiration. Furthermore, simulated increases in groundwater recharge were masked by decreases in water levels, stream baseflow, and groundwater storage resulting from surface water diversions and groundwater pumping. These results indicate that there is an imbalance between water supply and demand in this regional, semi‐arid aquifer. Current water management practices may not be sustainable into the far future and comprehensive action should be taken to minimize this water budget imbalance.     

Predicting Pesticide Attenuation in a Fractured Aquifer Using Lumped‐Parameter Models

Neighboring springs draining fractured‐rock aquifers can display large differences in water quality and flow regime, depending on local variations of the connectivity and the aperture size distribution of the fracture network. Consequently, because homogeneous equivalent parameters cannot be assumed a priori for the entire regional aquifer, the vulnerability to pollution of such springs has to be studied on a case by case basis. In this paper, a simple lumped‐parameter model usually applied to estimate the mean transit time of water (or tracer) is presented. The original exponential piston‐flow model was modified to take land‐use distribution into account and applied to predict the evolution of atrazine concentration in a series of springs draining a fractured sandstone aquifer in Luxembourg, where despite a nationwide ban in 2005, atrazine concentrations still had not begun to decrease in 2009. This persistence could be explained by exponentially distributed residence times in the aquifer, demonstrating that in some real world cases, models based on the groundwater residence time distribution can be a powerful tool for trend reversal assessments as recommended for instance by current European Union guidelines.     

Analytic element modelling for strip aquifers

The Analytic Element Method (AEM) provides a convenient tool for groundwater flow analysis in unbounded continuous domains. The AEM is based on the superposition of analytic functions, known as elements, useful at both regional and local scales. In this study, analytic elements for strip aquifers are presented. Such aquifers occur in riverine or coastal deposits and in outcrop zones of confined aquifers. Local flow field is modelled indirectly, using a reference plane related to the aquifer domain through the Schwarz‐Christoffel transform. The regional flow is obtained as a solution of the one‐dimensional flow equation. The proposed methodology was tested by modelling two hypothetical situations, which were compared to exact solutions. It is shown that regional boundaries can be reproduced exactly while local fields are adequately reproduced with analytic elements. The developed elements are applied to simulate a real flow field in northeastern Brazil showing good agreement with measured water levels. Copyright © 2011 John Wiley & Sons, Ltd.     

Hydrology of a headwater basin wetland: Groundwater discharge and wetland maintenance

The link between groundwater and surface hydrology in a small headwater drainage basin in the zone of glacial deposition of southern Ontario south of the Precambrian Shield was examined for two years. The basin is situated in a discharge zone of a regional aquifer and contains a small treed spring-fed swamp. The swamp exists because of the groundwater and has little effect on the maintenance of streamflow. Groundwater input to the swamp is an order of magnitude larger than precipitation. Groundwater of local and regional origin passes through the swamp by two routes: surface streamlets, where groundwater that emerges at specific seepage points in the swamp is conveyed over the ground surface with little interaction with the swamp itself, and by diffuse seepage in the swamp and through the bed of the stream. While the diffuse seepage input is the smaller component of groundwater it maintains the swamp's saturation. Groundwater input to the swamp from the specific seepage points and diffuse flow varies little over a year; therefore the saturation of the swamp and baseflow from the basin display little seasonal variation compared to other wetland types. The existence of the valley bottom in the headwater basin alters the seasonal and storm hydrology and is important to biogeochemical transformation of emerging groundwater.     

Dynamics of a headwater system and peatland under current conditions and with climate change

Interactions between headwater aquifers and peatlands have received limited scientific attention. Hydrological stresses, including those related to climate change, may adversely impact these interactions. In this study, the dynamics of a southern Québec headwater system where a peatland is present is simulated under current conditions and with climate change. The model is calibrated in steady state on field‐measured data and provides satisfactory results for transient‐state conditions. Under current conditions, simulations confirm that the peatland is fed by the fractured bedrock aquifer year‐round and provides continuous baseflow to its outlets. Climate change is simulated through its impact on groundwater recharge. Predicted precipitation and temperature data from a suite of regional climate model scenarios provide a net precipitation variation range from +10% to ?30% for the 2041–2070 horizon. Calibrated recharge is modified within this range to perform a sensitivity analysis of the headwater model to recharge variations (+10%, ?15% and ?30%). Total contribution from the aquifer to rivers and streams varies from +14% to ?44% of the baseline for +10% to ?30% recharge changes from spring 2010 data, for example. With higher recharge, the peatland receives more groundwater, which could significantly change its vegetation pattern and eventually ecosystem functions. For a ?30% recharge, the peatland becomes perched above the aquifer during the summer, fall and winter. Recharge reductions also induce sharp declines in groundwater levels and drying streams. Copyright © 2013 John Wiley & Sons, Ltd.     

Normal-mode analysis of the structure of baseflow-recession curves

It has been observed for several years that semi-logarithmic plots of the baseflow recession for many streams which partially penetrate aquifers in the United Kingdom are not single straight-line plots. Simple modelling studies indicate that this behaviour, for typical values of aquifer parameters, occurs for even the simplest of groundwater systems. In general such studies show that the baseflow component arising from even a structurally uncomplicated aquifer should be viewed as a superposition of many, distinct exponential terms. The origin of these terms does not necessarily lie in a varied hydrogeological structure, but can arise purely from the dynamics of groundwater flow. Such a form for the baseflow recession can be used to perform empirical fits to observed data. Moreover, where more than one exponential term is contributing to the total baseflow, great care must be taken in the interpretation of associated data.