Sources and sinks of nutrients in a New Zealand hill pasture catchment I. Stormflow generation

The processes of stormflow generation were studied in a hill pasture catchment near Hamilton, New Zealand. Although rainfall was relatively evenly distributed throughout the year, stormflow was highly seasonal and over 65 per cent occurred during the winter. Three main processes contributing to stormflow were identified which could be related to soil type and physiographic position. On gleyed soils derived from rhyolitic colluvium, saturation overland flow was the dominant process. Hydrographs from ‘Whipkey’ throughflow troughs also indicated that there was a subsurface response (saturated wedge) from this soil type. On steeper convex slopes, more permeable soils were derived from weathered greywacke. The presence of ephemeral springs on the hillslopes and direct observation during storm events indicated that storm runoff was generated as return flow from this soil. It was noted that nitrate concentrations from subsurface sources were 5–10 times higher than surface runoff. This difference in concentration was utilized in a chemical mixing equation which partitioned stormflow sources. This was compared with the stormflow predicted from rain falling on to saturated areas. There was good agreement between the two models for winter-spring events with respect to the volumes of surface runoff predicted, however the saturated areas model underestimated total stormflow. The results of the study are briefly discussed in terms of the potential for water quality management.     

Parametric study of a physically-based,plot-scale hillslope hydrological model through virtual experiments

Abstract

A physically-based hillslope hydrological model with shallow overland flow and rapid subsurface stormflow components was developed and calibrated using field experiments conducted on a preferential path nested hillslope in northeast India. Virtual experiments were carried out to perform sensitivity analysis of the model using the automated parameter estimation (PEST) algorithm. Different physical parameters of the model were varied to study the resulting effects on overland flow and subsurface stormflow responses from the theoretical hillslopes. It was observed that topographical shapes had significant effects on overland flow hydrographs. The slope profiles, surface storage, relief, rainfall intensity and infiltration rates primarily controlled the overland flow response of the hillslopes. Prompt subsurface stormflow responses were mainly dominated by lateral preferential flow, as soil matrix flow rates were very slow. Rainfall intensity and soil macropore structures were the most influential parameters on subsurface stormflow. The number of connected soil macropores was a more sensitive parameter than the size of macropores. In hillslopes with highly active vertical and lateral preferential pathways, saturation excess overland flow was not evident. However, saturation excess overland flow was generated if the lateral macropores were disconnected. Under such conditions, rainfall intensity, duration and preferential flow rate governed the process of saturation excess overland flow generation from hillslopes.

Editor D. Koutsoyiannis; Associate editor C. Perrin     

Stormflow generation in a small rainforest catchment in the Luquillo Experimental Forest,Puerto Rico

Various complementary techniques were used to investigate the stormflow generating processes in a small headwater catchment in northeastern Puerto Rico. Over 100 samples were taken of soil matrix water, macropore flow, streamflow and precipitation, mainly during two storms of contrasting magnitude, for the analysis of calcium, magnesium, silicon, potassium, sodium and chloride. These were combined with hydrometric information on streamflow, return flow, precipitation, throughfall and soil moisture to distinguish water following different flow paths. Geo‐electric sounding was used to survey the subsurface structure of the catchment, revealing a weathering front that coincided with the elevation of the stream channel instead of running parallel to surface topography. The hydrometric data were used in combination with soil physical data, a one‐dimensional soil water model (VAMPS ) and a three‐component chemical mass‐balance mixing model to describe the stormflow response of the catchment. It is inferred that most stormflow travelled through macropores in the top 20 cm of the soil profile. During a large event, saturation overland flow also accounted for a considerable portion of the stormflow, although it was not possible to quantify the associated volume fully. Although the mass‐balance mixing model approach gave valuable information about the various flow paths within the catchment, it was not possible to distill the full picture from the model alone; additional hydrometric and soil physical evidence was needed to aid in the interpretation of the model results. Copyright © 2004 John Wiley & Sons, Ltd.     

Soil pipe contribution to steady subsurface stormflow

Flow through a saturated idealized hillslope with a single soil pipe was simulated using a finite difference solution to the equation for three-dimensional Darcian flow in saturated heterogeneous media. The proportions of hillslope discharge originating from flow through the soil matrix and from flow through a soil pipe were determined, considering such factors as the radius, depth and length of the pipe, pipe spacing, and the length and slope of the hillslope. Results demonstrate that soil pipes can contribute a significant amount, and in many instances, the majority of total subsurface stormflow.     

A hydrochemical approach to quantify the role of return flow in a surface flow‐dominated catchment

Stormflow generation in headwater catchments dominated by subsurface flow has been studied extensively, yet catchments dominated by surface flow have received less attention. We addressed this by testing whether stormflow chemistry is controlled by either (a) the event‐water signature of overland flow, or (b) the pre‐event water signature of return flow. We used a high‐resolution hydrochemical data set of stormflow and end‐members of multiple storms in an end‐member mixing analysis to determine the number of end‐members needed to explain stormflow, characterize and identify potential end‐members, calculate their contributions to stormflow, and develop a conceptual model of stormflow. The arrangement and relative positioning of end‐members in stormflow mixing space suggest that saturation excess overland flow (26–48%) and return flow from two different subsurface storage pools (17–53%) are both similarly important for stormflow. These results suggest that pipes and fractures are important flow paths to rapidly release stored water and highlight the value of within‐event resolution hydrochemical data to assess the full range and dynamics of flow paths.     

Streamflow generation in a headwater basin on the precambrian shield

Current conceptual runoff models hypothesize that stormflow generation on the Canadian Shield is a combination of subsurface stormflow and saturation overland flow. This concept was tested during spring runoff in a small (3.3 ha) headwater basin using: (1) isotopic and chemical hydrograph separation and (2) field mapping and direct tracing of saturated areas. Isotopic and chemical hydrograph separation indicated three runoff components: (1) pre-melt subsurface flow; (2) subsurface flow of new (event) water; and (3) direct precipitation on to saturated areas (DPS). During early thaw-freeze cycles, their relative contributions to total flow remained constant (65 per cent, 30 per cent, and 5 per cent respectively). It is hypothesized that lateral flow along the bedrock/mineral soil interface, possibly through macropores, supplied large volumes of subsurface flow (of both old and new water) rapidly to the stream channel. Much higher contributions of DPS were observed during an intensive rain-on-snow event (15 per cent of total flow). Mapping and direct tracing of saturated areas using lithium bromide, suggested that saturated area size was positively correlated to stream discharge but its response lagged behind that of discharge. These observations suggest that the runoff mechanisms, and hence the sources of stream flow, will vary depending on storm characteristics.     

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.     

Evaluating the spatial distribution of water balance in a small watershed,Pennsylvania

A conceptual water‐balance model was modified from a point application to be distributed for evaluating the spatial distribution of watershed water balance based on daily precipitation, temperature and other hydrological parameters. The model was calibrated by comparing simulated daily variation in soil moisture with field observed data and results of another model that simulates the vertical soil moisture flow by numerically solving Richards' equation. The impacts of soil and land use on the hydrological components of the water balance, such as evapotranspiration, soil moisture deficit, runoff and subsurface drainage, were evaluated with the calibrated model in this study. Given the same meteorological conditions and land use, the soil moisture deficit, evapotranspiration and surface runoff increase, and subsurface drainage decreases, as the available water capacity of soil increases. Among various land uses, alfalfa produced high soil moisture deficit and evapotranspiration and lower surface runoff and subsurface drainage, whereas soybeans produced an opposite trend. The simulated distribution of various hydrological components shows the combined effect of soil and land use. Simulated hydrological components compare well with observed data. The study demonstrated that the distributed water balance approach is efficient and has advantages over the use of single average value of hydrological variables and the application at a single point in the traditional practice. Copyright © 2000 John Wiley & Sons, Ltd.     

Evaluating a new landform evolution model: A case study using a proposed mine rehabilitation landform

Landform evolution models are powerful tools for determining long-term erosional stability and denudation rates spanning geological timescales. SIBERIA, CAESAR and CHILD are examples of these model. The newly developed State Space Soil Production and Assessment Model (SSSPAM) coupled soilscape-landform evolution model has the ability to assess overall erosion rates of catchment scale landforms either using short-term precipitation events, variable precipitation or time-averaged precipitation (annual average). In addition, SSSPAM has the capability of developing the subsurface soil profile through weathering and armouring. In SSSPAM, physical processes of pedogenesis such as erosion and armouring, diffusion, sediment deposition and weathering are modelled using a state space matrix approach. In this article we simulate the short-term evolution (100 years) of a proposed post-mining landform using both SIBERIA and SSSPAM and compare the erosion and sediment output results. For the short-term simulations SSSPAM's armouring capability was disabled. The models were then used to simulate the evolution of the catchment for 10,000 years. Results demonstrate that the short-term SSSPAM simulation results compare well with the results from the established landform evolution model SIBERIA. The long-term armouring disabled SSSPAM simulations produces simulated erosion rates comparable with SIBERIA simulations both of which are similar to upper limit of field measured denudation rates. The SSSPAM simulation using armouring demonstrated that armouring reduced the erosion rate of the catchment by a factor of 4 which is comparable with the lower limit of field measured denudation rates. This observation emphasizes the importance of armouring in long-term evolution of landforms. Soil profile cross-sections developed from the same results show that SSSPAM can also reproduce subsurface soil evolution and stratification and spatial variability of soil profile characteristics typically observed in the field.     

The influence of preferential flow on hillslope hydrology in a semi‐arid watershed (in the Spanish Dehesas)

Preferential flow is known to influence hillslope hydrology in many areas around the world. Most research on preferential flow has been performed in temperate regions. Preferential infiltration has also been found in semi‐arid regions, but its impact on the hydrology of these regions is poorly known. The aim of this study is to describe and quantify the influence of preferential flow on the hillslope hydrology from small scale (infiltration) to large scale (subsurface stormflow) in a semi‐arid Dehesa landscape. Precipitation, soil moisture content, piezometric water level and discharge data were used to analyse the hydrological functioning of a catchment in Spain. Variability of soil moisture content during the transition from dry to wet season (September to November) within horizontal soil layers leads to the conclusion that there is preferential infiltration into the soils. When the rainfall intensity is high, a water level rapidly builds up in the piezometer pipes in the area, sometimes even reaching soil surface. This water level also drops back to bedrock within a few hours (under dry catchment conditions) to days (under wet catchment conditions). As the soil matrix is not necessarily wet while this water layer is built up, it is thought to be a transient water table in large connected pores which drain partly to the matrix, partly fill up bedrock irregularities and partly drain through subsurface flow to the channels. When the soil matrix becomes wetter the loss of water from macropores to the matrix and bedrock decreases and subsurface stormflow increases. It may be concluded that the hillslope hydrological system consists of a fine matrix domain and a macropore domain, which have their own flow characteristics but which also interact, depending on the soil matrix and macropore moisture contents. The macropore flow can result in subsurface flow, ranging from 13% contribution to total discharge for a large event of high intensity rainfall or high discharge to 80% of total discharge for a small event with low intensity rainfall or low discharge. During large events the fraction of subsurface stormflow in the discharge is suppressed by the large amount of surface runoff. Copyright © 2008 John Wiley & Sons, Ltd.     

Development and application of a spatially‐distributed Arctic hydrological and thermal process model (ARHYTHM)

A process‐based, spatially distributed hydrological model was developed to quantitatively simulate the energy and mass transfer processes and their interactions within arctic regions (arctic hydrological and thermal model, ARHYTHM). The model first determines the flow direction in each element, the channel drainage network and the drainage area based upon the digital elevation data. Then it simulates various physical processes: including snow ablation, subsurface flow, overland flow and channel flow routing, soil thawing and evapotranspiration. The kinematic wave method is used for conducting overland flow and channel flow routing. The subsurface flow is simulated using the Darcian approach. The energy balance scheme was the primary approach used in energy‐related process simulations (snowmelt and evapotranspiration), although there are options to model snowmelt by the degree‐day method and evapotranspiration by the Priestley–Taylor equation. This hydrological model simulates the dynamic interactions of each of these processes and can predict spatially distributed snowmelt, soil moisture and evapotranspiration over a watershed at each time step as well as discharge in any specified channel(s). The model was applied to Imnavait watershed (about 2·2 km²) and the Upper Kuparuk River basin (about 146 km²) in northern Alaska. Simulated results of spatially distributed soil moisture content, discharge at gauging stations, snowpack ablations curves and other results yield reasonable agreement, both spatially and temporally, with available data sets such as SAR imagery‐generated soil moisture data and field measurements of snowpack ablation, and discharge data at selected points. The initial timing of simulated discharge does not compare well with the measured data during snowmelt periods mainly because the effect of snow damming on runoff was not considered in the model. Results from the application of this model demonstrate that spatially distributed models have the potential for improving our understanding of hydrology for certain settings. Finally, a critical component that led to the performance of this modelling is the coupling of the mass and energy processes. Copyright © 2000 John Wiley & Sons, Ltd.     

Development and testing of a physically based,three-dimensional model of surface and subsurface hydrology

We present a numerical, catchment-scale model that solves flow equations of surface and subsurface flow in a three-dimensional domain. Surface flow is described by the two-dimensional parabolic approximation of the St. Venant equation, using Manning’s equation of motion; subsurface flow is described by the three-dimensional Richards’ equation for the unsaturated zone and by three-dimensional Darcy’s law for the saturated zone, using an integrated finite difference formulation. The hydrological component is a dynamic link library implemented within a comprehensive model which simulates surface energy, radiation budget, snow melt, potential evapotranspiration, plant development and plant water uptake. We tested the model by comparing distributed and integrated three-dimensional simulated and observed perched water depth (PWD), stream flow data, and soil water contents for a small catchment. Additional tests were performed for the snow melting algorithm as well as the different hydrological processes involved. The model successfully described the water balance and its components as evidenced by good agreement between measured and modelled data.     

Physically‐based modelling of double‐peak discharge responses at Slapton Wood catchment

Heavy winter rainfall produces double‐peak hydrographs at the Slapton Wood catchment, Devon, UK. The first peak is saturation‐excess overland flow in the hillslope hollows and the second (i.e. the delayed peak) is subsurface stormflow. The physically‐based spatially‐distributed model SHETRAN is used to try to improve the understanding of the processes that cause the double peaks. A three‐stage (multi‐scale) approach to calibration is used: (1) water balance validation for vertical one‐dimensional flow at arable, grassland and woodland plots; (2) two‐dimensional flow for cross‐sections cutting across the stream valley; and (3) three‐dimensional flow in the full catchment. The main data are for rainfall, stream discharge, evaporation, soil water potential and phreatic surface level. At each scale there was successful comparison with measured responses, using as far as possible parameter values from measurements. There was some calibration but all calibrated values at one scale were used at a larger scale. A large proportion of the subsurface runoff enters the stream from three dry valleys (hillslope hollows), and previous studies have suggested convergence of the water in the three large hollows as being the major mechanism for the production of the delayed peaks. The SHETRAN modelling suggests that the hillslopes that drain directly into the stream are also involved in producing the delayed discharges. The model shows how in the summer most of the catchment is hydraulically disconnected from the stream. In the autumn the catchment eventually ‘wets up’ and shallow subsurface flows are produced, with water deflected laterally along the soil‐bedrock interface producing the delayed peak in the stream hydrograph. Copyright © 2007 John Wiley & Sons, Ltd.     

INFERRING HYDROLOGICAL PROCESSES IN A TEMPERATE BASIN USING ISOTOPIC AND GEOCHEMICAL HYDROGRAPH SEPARATION: A RE-EVALUATION

Simultaneous monitoring of conservative and non-conservative tracers in streamflow offers a valuable means of obtaining information on the age and flow paths of water reaching the basin outlet. Previous studies of stormflow generation in a small forested basin on the Canadian Shield used isotopic (IHS) and geochemical hydrograph separations (GHS) to infer that some event water during snowmelt reaches the stream via subsurface pathways, and that surface water runoff is generated by direct precipitation on to saturated areas (DPSA) in the stream valley. These hypotheses were tested for rainfall inputs using simultaneous IHS (¹⁸O) and GHS (dissolved silica) of basin stormflow, supplemented by hydrochemical and hydrometric data from throughflow troughs installed on basin slopes. Comparison of pre-event and subsurface water hydrographs did not provide conclusive evidence for subsurface movement of event water to the stream, owing to the appreciable uncertainty associated with the hydrograph separations. However, IHSs of runoff at the soil–bedrock interface on basin slopes indicated that event water comprised 25–50% of total runoff from areas with deep soil cover, and that these contributions supplied event water flux from the basin in excess of that attributable to DPSA. The surface water component of stormflow estimated from the GHS was also largely the result of DPSA. GHS assumes that dissolved silica is rapidly and uniformly taken up by water infiltrating the soil and that water moving via surface pathways retains the low dissolved silica level of rainfall; however, neither assumption was supported by the hillslope results. Instead, results suggest that the observed depression of silica levels in basin stormflow previously attributed to dilution by DPSA was partly a function of transport of dilute event water to the channel via preferential pathways. Implications of these results for the general use of simultaneous IHS and GHS to infer hydrological processes are discussed. © 1997 by John Wiley & Sons, Ltd.     

Lumped hysteretic model for subsurface stormflow developed using downward approach

The terms ‘downward’ and ‘upward’ (synonymous with ‘top‐down’ and ‘bottom‐up’ respectively) are sometimes used when describing methods for developing hydrological models. A downward approach is used here to develop a lumped catchment‐scale model for subsurface stormflow at the 0·94 km² Slapton Wood catchment. During the development, as few assumptions as possible are made about the behaviour of subsurface stormflow at the catchment scale, and no assumptions are made about its behaviour at smaller scales. (In an upward approach, in contrast, the modelling would be based on assumptions about, and data for, the behaviour at smaller scales, such as the hillslope, plot, and point scales.) The model has a single store with a relatively simple relationship between discharge and storage, based on equations describing hysteretic patterns seen in a graph of discharge against storage. Double‐peaked hydrographs have been observed at the catchment outlet. Rainfall on the channel and infiltration‐excess and saturation‐excess runoff give a rapid response, and shallow subsurface stormflow gives a delayed response. Hydrographs are successfully simulated for the large delayed responses observed in 1971–1980 and 1989–1991, then a lumped model for the rapid response is coupled to the lumped hysteretic model and some double‐peaked hydrographs simulated. A physical interpretation is developed for the lumped hysteretic model, making use of information on patterns of perched saturation observed in 1982 on a hillslope at the Slapton Wood catchment. Downward and upward approaches are complementary, and the most robust way to develop and improve lumped catchment models is to iterate between downward and upward steps. Possible next steps are described. Copyright © 2006 John Wiley & Sons, Ltd.     

Stormflow generation involving pipe flow in a zero‐order basin of Peninsular Malaysia

Hydrological responses in a zero‐order basin (ZOB), a portion of whose discharge emerged via preferential flow through soil pipes, were examined over a 2‐year period in Peninsular Malaysia to elucidate primary stormflow generation processes. Silicon (Si) and specific conductance (EC) in various runoff components were also measured to identify their sources. ZOB flow response was dependent on antecedent precipitation amount; runoff increased linearly with precipitation during events >20 mm in relatively wet antecedent moisture conditions. Runoff derived from direct precipitation falling onto saturated areas accounted for <0·2% of total ZOB flow volume during the study period, indicating the predominance of subsurface pathways in ZOB flow. ZOB flow (high EC and low Si) was distinct from perennial baseflow via bedrock seepage (low EC and high Si) 5 m downstream of the ZOB outlet. Pipe flow responded quickly to ZOB flow rate and was characterized by a threshold flow capacity unique to each pipe. Piezometric data and pipe flow records demonstrated that pipes located deeper in the soil initiated first, followed by those at shallower depths; initiation of pipe flow corresponded to shallow groundwater rise above the saprolite‐soil interface. Chemical signatures of pipe flow were similar to each other and to the ZOB flow, suggesting that the sources were well‐mixed soil‐derived shallow groundwater. Based upon the volume of pipe flow during storms, the combined contribution of the pipes monitored accounted for 48% of total ZOB flow during the study period. Our results suggest that shallow groundwater, possibly facilitated by preferential flow accreted above the saprolite–soil interface, provides dominant stormflow, and that soil pipes play an important role in the rapid delivery of solute‐rich water to the stream system. Copyright © 2006 John Wiley & Sons, Ltd.     

Effect of bedrock permeability on subsurface stormflow and the water balance of a trenched hillslope at the Panola Mountain Research Watershed,Georgia, USA

The effect of bedrock permeability on subsurface stormflow initiation and the hillslope water balance is poorly understood. Previous hillslope hydrological studies at the Panola Mountain Research Watershed (PMRW), Georgia, USA, have assumed that the bedrock underlying the trenched hillslope is effectively impermeable. This paper presents a series of sprinkling experiments where we test the bedrock impermeability hypothesis at the PMRW. Specifically, we quantify the bedrock permeability effects on hillslope subsurface stormflow generation and the hillslope water balance at the PMRW. Five sprinkling experiments were performed by applying 882–1676 mm of rainfall over a ～5·5 m × 12 m area on the lower hillslope during ～8 days. In addition to water input and output captured at the trench, we measured transpiration in 14 trees on the slope to close the water balance. Of the 193 mm day^?1 applied during the later part of the sprinkling experiments when soil moisture changes were small, <14 mm day^?1 was collected at the trench and <4 mm day^?1 was transpired by the trees, with residual bedrock leakage of >175 mm day^?1 (91%). Bedrock moisture was measured at three locations downslope of the water collection system in the trench. Bedrock moisture responded quickly to precipitation in early spring. Peak tracer breakthrough in response to natural precipitation in the bedrock downslope from the trench was delayed only 2 days relative to peak tracer arrival in subsurface stormflow at the trench. Leakage to bedrock influences subsurface stormflow at the storm time‐scale and also the water balance of the hillslope. This has important implications for the age and geochemistry of the water and thus how one models this hillslope and watershed. Copyright © 2006 John Wiley & Sons, Ltd.     

Hydrological behaviour and modelling of a volcanic tropical cultivated catchment

The hydrological behaviour of the cultivated Féfé catchment (17·8 ha) on the tropical volcanic island of Guadeloupe was studied to identify flow paths, to quantify water fluxes, and finally, to build a lumped model to simulate discharge and piezometer levels. The approach combined two steps, an experimental step and a modelling step, which covered two time scales, the annual and the storm event scale. The hydrological measurements were conducted over 2 years. The Féfé catchment is characterized by heavy rainfall (4229 mm year^?1) on permeable Andosols; the results showed that underground flow paths involved two overlapping aquifers, and that the annual water balance in 2003 was shared among outflows of the deep aquifer (42%), evapotranspiration (31%), and streamflow (27%). On the event scale, the surface runoff coefficient ranges between 6·2% and 24·4% depending on antecedent dry or wet moisture conditions. Hortonian overland flow predominated over subsurface and saturation overland flow processes. Recharge of the shallow aquifer is mainly governed by a constant infiltration capacity of the Andosols with depth in the vadose zone. Outflows of this shallow aquifer were the baseflow of the main stream and the recharge of the deep aquifer. Volcanic deposits at Féfé promoted the underground flow path, and cultivated areas seemed to explain the high stormflow values relative to other tropical small catchments under rain forest. A conceptual lumped model integrating runoff, infiltration, evapotranspiration, and fluctuations of the two overlapping aquifers was developed. The model has six parameters and was calibrated and validated on the hydrograph at the outlet and on the two piezometers of the shallow and the deep aquifers. The results show fair to good agreement between measured and simulated variables, and consequently, the model was consistent with the main hydrological processes observed from experimental results in wet conditions. Copyright © 2008 John Wiley & Sons, Ltd.