On the Panday and Huyakorn surface–subsurface hydrology test case: analysis of internal flow dynamics

A number of previous studies using models of integrated surface‐subsurface hydrology have adopted the Panday and Huyakorn (P&H) tilted V‐catchment test case (Panday S, Huyakorn PS. 2004. A fully coupled spatially distributed model for evaluating surface/subsurface flow. Advances in Water Resources 27: 361–382) to show inter‐code comparability. The P&H test case is used to evaluate models that simulate a broad range of hydrological processes, and yet only the catchment outflow hydrograph has been presented as verification of the consistency between codes. Therefore, a more comprehensive evaluation of the surface‐subsurface hydrology of the P&H case is needed. This study explores the internal catchment functioning of the P&H case, using the popular catchment simulator MODHMS. The processes leading to streamflow generation in the model are illustrated, including separation of overland flow (OLF) and groundwater discharge to the stream. The results identify non‐physical flow processes due to the problem set‐up, and modifications to the P&H case are suggested that include changes to stream roughness and incision of the stream channel to overcome these shortcomings. A modified P&H case produced more plausible transfers between OLF and the stream, and an increased groundwater discharge to the stream (6·5% of streamflow in the modified case compared to 0·5% in the original case). Despite changes to internal flow processes, near‐identical outflow hydrographs were obtained, showing the importance of considering and comparing internal flow processes when using surface‐subsurface hydrology test cases to evaluate integrated hydrological simulators. Copyright © 2011 John Wiley & Sons, Ltd.     

Water and salt balance modelling to predict the effects of land-use changes in forested catchments. 2. coupled model of water and salt balances

A long-term salt balance model is coupled with the small catchment water balance model presented in Part 1 of this series of papers. The salt balance model was designed as a simple robust, conceptually based model of the fundamental salt fluxes and stores in forested and cleared catchments. The model has four interdependent stores representing salt storage in the unsaturated zone, the deep permanent saturated groundwater system, the near-stream perched groundwater system and in a ‘salt bulge’ just above the permanent water-table. The model has performed well in simulations carried out on Salmon and Wights, two small experimental catchments in south-west Western Australia. When applied to Wights catchment the salt balance model was able to predict the stream salinities prior to clearing of native forests, and the increased salinities after the clearing.     

A comprehensive assessment framework for attributing trends in streamflow and groundwater storage to climatic and anthropogenic changes: A case study in the typical semi-arid catchments of southern India

The clearest signs of hydrologic change can be observed from the trends in streamflow and groundwater levels in a catchment. During 1980–2007, significant declines in streamflow (−3.03 mm/year) and groundwater levels (−0.22 m/year) were observed in Himayat Sagar (HS) catchment, India. We examined the degree to which hydrologic changes observed in the HS catchment can be attributed to various internal and external drivers of change (climatic and anthropogenic changes). This study used an investigative approach to attribute hydrologic changes. First, it involves to develop a model and test its ability to predict hydrologic trends in a catchment that has undergone significant changes. Second, it examines the relative importance of different causes of change on the hydrologic response. The analysis was carried out using Modified Soil and Water Assessment Tool (SWAT), a semi-distributed rainfall-runoff model coupled with a lumped groundwater model for each sub- catchment. The model results indicated that the decline in potential evapotranspiration (PET) appears to be partially offset by a significant response to changes in rainfall. Measures that enhance recharge, such as watershed hydrological structures, have had limited success in terms of reducing impacts on the catchment-scale water balance. Groundwater storage has declined at a rate of 5 mm/y due to impact of land use changes and this was replaced by a net addition of 2 mm/y by hydrological structures. The impact of land use change on streamflow is an order of magnitude larger than the impact of hydrological structures and about is 2.5 times higher in terms of groundwater impact. Model results indicate that both exogenous and endogenous changes can have large impacts on catchment hydrology and should be considered together. The proposed comprehensive framework and approach demonstrated here is valuable in attributing trends in streamflow and groundwater levels to catchment climatic and anthropogenic changes.     

Variable water‐saturated areas and streamflow generation in the small Ringelbach catchment (Vosges Mountains,France): the master recession curve as an equilibrium curve for interactions between atmosphere,surface and ground waters

In the small Ringelbach research catchment, where studies on the water cycle components in a granitic mountainous environment have been conducted since 1976, the water‐saturated areas that are hydraulically connected to the outlet play a major role in the streamflow generation, as it is here that complex interactions between atmosphere, surface and ground waters take place. During baseflow recession periods, which may last several months between two groundwater recharge events, the atmospheric inputs of water and energy on these contributing areas only explain the streamflow fluctuations observed around the master recession curve, which defines the groundwater contribution: fluctuating above it in the case of precipitation input on these areas, below it in the case of evaporation output from these areas. Streamflow may therefore largely deviate from the master recession curve in the case of long, hot, dry spells. Detailed mapping has shown that their variable extent is well related to baseflow by a loglinear curve. On the other hand, a synthetic master recession curve, well fitted by a second‐order hyperbolic function, has been obtained from numerous pure recession periods. Both based on these two curves, a simple procedure and a simple model have been used to (i) validate the hypothesis that the connected saturated areas are the only permanent variable contributing areas and (ii) simulate the daily streamflow volumes over long baseflow recession periods by a water balance of the aquifer below these areas only. The storm runoff ratio for small to moderate rainfall events is indeed corresponding to the catchment saturated fraction at that time. The volume of daily streamflow oscillations is indeed corresponding to the evaporation at the potential rate from the saturated areas only. In both cases, streamflow naturally tends towards the master recession curve after the end of any atmospheric perturbation. Introducing these findings into TOPMODEL led to significantly improved simulation results during baseflow recession periods. The master recession curve may therefore be considered as a dynamic equilibrium curve. Together with the relationship between saturated extent and baseflow, it provides the main characteristics necessary to understand and model the interactions at this complex interface and the resulting daily streamflow variations during baseflow recession periods in this type of catchment. Copyright © 2016 John Wiley & Sons, Ltd.     

Field investigation and modelling of coupled stream discharge and shallow water‐table dynamics in a small Mediterranean catchment (Sardinia)

In the semi‐arid Mediterranean environment, the rainfall–runoff relationships are complex because of the markedly irregular patterns in rainfall, the seasonal mismatch between evaporation and rainfall, and the spatial heterogeneity in landscape properties. Watersheds often display considerable non‐linear threshold behavior, which still make runoff generation an open research question. Our objectives in this context were: to identify the primary processes of runoff generation in a small natural catchment; to test whether a physically based model, which takes into consideration only the primary processes, is able to predict spatially distributed water‐table and stream discharge dynamics; and to use the hydrological model to increase our understanding of runoff generation mechanisms. The observed seasonal dynamics of soil moisture, water‐table depth, and stream discharge indicated that Hortonian overland‐flow was negligible and the main mechanism of runoff generation was saturated subsurface‐flow. This gives rise to base‐flow, controls the formation of the saturated areas, and contributes to storm‐flow together with saturation overland‐flow. The distributed model, with a 1D scheme for the kinematic surface‐flow, a 2D sub‐horizontal scheme for the saturated subsurface‐flow, and ignoring the unsaturated flow, performed efficiently in years when runoff volume was high and medium, although there was a smoothing effect on the observed water‐table. In dry years, small errors greatly reduced the efficiency of the model. The hydrological model has allowed to relate the runoff generation mechanisms with the land‐use. The forested hillslopes, where the calibrated soil conductivity was high, were never saturated, except at the foot of the slopes, where exfiltration of saturated subsurface‐flow contributed to storm‐flow. Saturation overland‐flow was only found near the streams, except when there were storm‐flow peaks, when it also occurred on hillslopes used for pasture, where soil conductivity was low. The bedrock–soil percolation, simulated by a threshold mechanism, further increased the non‐linearity of the rainfall–runoff processes. Copyright © 2013 John Wiley & Sons, Ltd.     

Multi-Constrained Catchment Scale Optimization of Groundwater Abstraction Using Linear Programming

Due to increasing water demands globally, freshwater ecosystems are under constant pressure. Groundwater resources, as the main source of accessible freshwater, are crucially important for irrigation worldwide. Over-abstraction of groundwater leads to declines in groundwater levels; consequently, the groundwater inflow to streams decreases. The reduction in baseflow and alteration of the streamflow regime can potentially have an adverse effect on groundwater-dependent ecosystems. A spatially distributed, coupled groundwater–surface water model can simulate the impacts of groundwater abstraction on aquatic ecosystems. A constrained optimization algorithm and a simulation model in combination can provide an objective tool for the water practitioner to evaluate the interplay between economic benefits of groundwater abstractions and requirements to environmental flow. In this study, a holistic catchment-scale groundwater abstraction optimization framework has been developed that allows for a spatially explicit optimization of groundwater abstraction, while fulfilling a predefined maximum allowed reduction of streamflow (baseflow [Q95] or median flow [Q50]) as constraint criteria for 1484 stream locations across the catchment. A balanced K-Means clustering method was implemented to reduce the computational burden of the optimization. The model parameters and observation uncertainties calculated based on Bayesian linear theory allow for a risk assessment on the optimized groundwater abstraction values. The results from different optimization scenarios indicated that using the linear programming optimization algorithm in conjunction with integrated models provides valuable information for guiding the water practitioners in designing an effective groundwater abstraction plan with the consideration of environmental flow criteria important for the ecological status of the entire system.     

Implementation of Solute Transport in the Vadose Zone into the “HYDRUS Package for MODFLOW”

The “HYDRUS package for MODFLOW” is an existing MODFLOW package that allows MODFLOW to simultaneously evaluate transient water flow in both unsaturated and saturated zones. The package is based on incorporating parts of the HYDRUS-1D model (to simulate unsaturated water flow in the vadose zone) into MODFLOW (to simulate saturated groundwater flow). The coupled model is effective in addressing spatially variable saturated-unsaturated hydrological processes at the regional scale. However, one of the major limitations of this coupled model is that it does not have the capability to simulate solute transport along with water flow and therefore, the model cannot be employed for evaluating groundwater contamination. In this work, a modified unsaturated flow and transport package (modified HYDRUS package for MODFLOW and MT3DMS) has been developed and linked to the three-dimensional (3D) groundwater flow model MODFLOW and the 3D groundwater solute transport model MT3DMS. The new package can simulate, in addition to water flow in the vadose zone, also solute transport involving many biogeochemical processes and reactions, including first-order degradation, volatilization, linear or nonlinear sorption, one-site kinetic sorption, two-site sorption, and two-kinetic sites sorption. Due to complex interactions at the groundwater table, certain modifications of the pressure head (compared to the original coupling) and solute concentration profiles were incorporated into the modified HYDRUS package. The performance of the newly developed model is evaluated using HYDRUS (2D/3D), and the results indicate that the new model is effective in simulating the movement of water and contaminants in the saturated-unsaturated flow domains.     

Modeling Variably Saturated Subsurface Solute Transport with MODFLOW‐UZF and MT3DMS

The MT3DMS groundwater solute transport model was modified to simulate solute transport in the unsaturated zone by incorporating the unsaturated‐zone flow (UZF1) package developed for MODFLOW. The modified MT3DMS code uses a volume‐averaged approach in which Lagrangian‐based UZF1 fluid fluxes and storage changes are mapped onto a fixed grid. Referred to as UZF‐MT3DMS, the linked model was tested against published benchmarks solved analytically as well as against other published codes, most frequently the U.S. Geological Survey's Variably‐Saturated Two‐Dimensional Flow and Transport Model. Results from a suite of test cases demonstrate that the modified code accurately simulates solute advection, dispersion, and reaction in the unsaturated zone. Two‐ and three‐dimensional simulations also were investigated to ensure unsaturated‐saturated zone interaction was simulated correctly. Because the UZF1 solution is analytical, large‐scale flow and transport investigations can be performed free from the computational and data burdens required by numerical solutions to Richards' equation. Results demonstrate that significant simulation runtime savings can be achieved with UZF‐MT3DMS, an important development when hundreds or thousands of model runs are required during parameter estimation and uncertainty analysis. Three‐dimensional variably saturated flow and transport simulations revealed UZF‐MT3DMS to have runtimes that are less than one tenth of the time required by models that rely on Richards' equation. Given its accuracy and efficiency, and the wide‐spread use of both MODFLOW and MT3DMS, the added capability of unsaturated‐zone transport in this familiar modeling framework stands to benefit a broad user‐ship.     

Water and salt balance modelling to predict the effects of land-use changes in forested catchments. 1. Small catchment water balance model

A long-term water balance model has been developed to predict the hydrological effects of land-use change (especially forest clearing) in small experimental catchments in the south-west of Western Australia. This small catchment model has been used as the building block for the development of a large catchment-scale model, and has also formed the basis for a coupled water and salt balance model, developed to predict the changes in stream salinity resulting from land-use and climate change. The application of the coupled salt and water balance model to predict stream salinities in two small experimental catchments, and the application of the large catchment-scale model to predict changes in water yield in a medium-sized catchment that is being mined for bauxite, are presented in Parts 2 and 3, respectively, of this series of papers. The small catchment model has been designed as a simple, robust, conceptually based model of the basic daily water balance fluxes in forested catchments. The responses of the catchment to rainfall and pan evaporation are conceptualized in terms of three interdependent subsurface stores A, B and F. Store A depicts a near-stream perched aquifer system; B represents a deeper, permanent groundwater system; and F is an intermediate, unsaturated infiltration store. The responses of these stores are characterized by a set of constitutive relations which involves a number of conceptual parameters. These parameters are estimated by calibration by comparing observed and predicted runoff. The model has performed very well in simulations carried out on Salmon and Wights, two small experimental catchments in the Collie River basin in south-west Western Australia. The results from the application of the model to these small catchments are presented in this paper.     

Hydrological response of a Brazilian semi‐arid catchment to different land use and climate change scenarios: a modelling study

Brazilian semi‐arid regions are characterized by water scarcity, vulnerability to desertification, and climate variability. The investigation of hydrological processes in this region is of major interest not only for water planning strategies but also to address the possible impact of future climate and land‐use changes on water resources. A hydrological distributed catchment‐scale model (DiCaSM) has been applied to simulate hydrological processes in a small representative catchment of the Brazilian northeast semi‐arid region, and also to investigate the impact of climate and land‐use changes, as well as changes associated with biofuel/energy crops production. The catchment is part of the Brazilian network for semi‐arid hydrology, established by the Brazilian Federal Government. Estimating and modelling streamflow (STF) and recharge in semi‐arid areas is a challenging task, mainly because of limitation in in situ measurements, and also due to the local nature of some processes. Direct recharge measurements are very difficult in semi‐arid catchments and contain a high level of uncertainty. The latter is usually addressed by short‐ and long‐time‐scale calibration and validation at catchment scale, as well as by examining the model sensitivity to the physical parameters responsible for the recharge. The DiCaSM model was run from 2000 to 2008, and streamflow was successfully simulated, with a Nash–Sutcliffe (NS) efficiency coefficient of 0·73, and R² of 0·79. On the basis of a range of climate change scenarios for the region, the DiCaSM model forecasted a reduction by 35%, 68%, and 77%, in groundwater recharge (GWR), and by 34%, 65%, and 72%, in streamflow, for the time spans 2010–2039, 2040–2069, and 2070–2099, respectively, could take place for a dry future climate scenario. These reductions would produce severe impact on water availability in the region. Introducing castor beans to the catchment would increase the GWR and streamflow, mainly if the caatinga areas would be converted into castor beans production. Changing an area of 1000 ha from caatinga to castor beans would increase the GWR by 46% and streamflow by 3%. If the same area of pasture is converted into castor beans, there would be an increase in GWR and streamflow by 24% and 5%, respectively. Such results are expected to contribute towards environmental policies for north‐east Brazil (NEB), and to biofuel production perspectives in the region. Copyright © 2010 John Wiley & Sons, Ltd.     

A quasi three-dimensional model for flow and transport in unsaturated and saturated zones: 1. Implementation of the quasi two-dimensional case

A quasi three-dimensional (QUASI 3-D) model is presented for simulating the subsurface water flow and solute transport in the unsaturated and in the saturated zones of soil. The model is based on the assumptions of vertical flow in the unsaturated zone and essentially horizontal groundwater flow. The 1-D Richards equation for the unsaturated zone is coupled at the phreatic surface with the 2-D flow equation for the saturated zone. The latter was obtained by averaging 3-D flow equation in the saturated zone over the aquifer thickness. Unlike the Boussinesq equation for a leaky-phreatic aquifer, the developed model does not contain a storage term with specific yield and a source term for natural replenishment. Instead it includes a water flux term at the phreatic surface through which the Richards equation is linked with the groundwater flow equation. The vertical water flux in the saturated zone is evaluated on the basis of the fluid mass balance equation while the horizontal fluxes, in that equation, are prescribed by Darcy law. A 3-D transport equation is used to simulate the solute migration. A numerical algorithm to solve the problem for the general quasi 3-D case was developed. The developed methodology was exemplified for the quasi 2-D cross-sectional case (QUASI2D). Simulations for three synthetic problems demonstrate good agreement between the results obtained by QUASI2D and two fully 2-D flow and transport codes (SUTRA and 2DSOIL). Yet, simulations with the QUASI2D code were several times faster than those by the SUTRA and the 2DSOIL codes.     

Response time and water origin in a steep nested catchment in the Italian Dolomites

In this study, we investigate the surface flow time of rise in response to rainfall and snowmelt events at different spatial scales and the main sources originating channel runoff and spring water in a steep nested headwater catchment (Rio Vauz, Italian Dolomites), characterized by a marked elevation gradient. We monitored precipitation at different elevations and measured water stage/streamflow at the outlet of two rocky subcatchments of the same size, representative of the upper part of the catchment dominated by outcropping bedrock, at the outlet of a soil‐mantled and vegetated subcatchment of similar size but different morphology, and at the outlet of the main catchment. Hydrometric data are coupled with stable isotopes and electrical conductivity sampled from different water sources during five years, and used as tracers in end‐member mixing analysis, application of two component mixing models and analysis of the slope of the dual‐isotope regression line. Results reveal that times of rise are slightly shorter for the two rocky subcatchments, particularly for snowmelt and mixed rainfall/snowmelt events, compared to the soil‐mantled catchment and the entire Rio Vauz Catchment. The highly‐variable tracer signature of the different water sources reflects the geomorphological and geological complexity of the study area. The principal end‐members for channel runoff and spring water are identified in rainfall and snowmelt, which are the dominant water sources in the rocky upper part of the study catchment, and soil water and shallow groundwater, which play a relevant role in originating baseflow and spring water in the soil‐mantled and vegetated lower part of the catchment. Particularly, snowmelt contributes up to 64 ± 8% to spring water in the concave upper parts of the catchment and up to 62 ± 11% to channel runoff in the lower part of the catchment. These results offer new experimental evidences on how Dolomitic catchments capture and store rain water and meltwater, releasing it through a complex network of surface and subsurface flow pathways, and allow for the construction of a preliminary conceptual model on water transmission in snowmelt‐dominated catchments featuring marked elevation gradients.     

MODFLOW‐Based Coupled Surface Water Routing and Groundwater‐Flow Simulation

In this paper, we present a flexible approach for simulating one‐ and two‐dimensional routing of surface water using a numerical surface water routing (SWR) code implicitly coupled to the groundwater‐flow process in MODFLOW. Surface water routing in SWR can be simulated using a diffusive‐wave approximation of the Saint‐Venant equations and/or a simplified level‐pool approach. SWR can account for surface water flow controlled by backwater conditions caused by small water‐surface gradients or surface water control structures. A number of typical surface water control structures, such as culverts, weirs, and gates, can be represented, and it is possible to implement operational rules to manage surface water stages and streamflow. The nonlinear system of surface water flow equations formulated in SWR is solved by using Newton methods and direct or iterative solvers. SWR was tested by simulating the (1) Lal axisymmetric overland flow, (2) V‐catchment, and (3) modified Pinder‐Sauer problems. Simulated results for these problems compare well with other published results and indicate that SWR provides accurate results for surface water‐only and coupled surface water/groundwater problems. Results for an application of SWR and MODFLOW to the Snapper Creek area of Miami‐Dade County, Florida, USA are also presented and demonstrate the value of coupled surface water and groundwater simulation in managed, low‐relief coastal settings.     

Influence of bedrock groundwater on streamflow characteristics in a volcanic catchment

Groundwater movements in volcanic mountains and their effects on streamflow discharge and representative elementary area (REA) have remained largely unclear. We surveyed the discharge and chemical composition of spring and stream water in two catchments: the Hontani river (NR) catchment (6.6 km²) and the Hosotani river (SR) catchment (4.0 km²) at the southern part of Daisen volcano, Japan. Daisen volcano is a young volcano (17 × 10³ years) at an early stage of erosion. Our study indicated that deep groundwater that moved through thick lava and pyroclastic flows and that could not be explained by shallow movements controlled by surface topography contributed dominantly to streamflow at larger catchment areas. At the NR catchment, the deep groundwater contribution clearly increased at a catchment boundary defined by an area of 3.0 km² and an elevation of 800 m. At the SR catchment, the contribution deep groundwater to the stream also increased suddenly at a boundary threshold of 2.0 and 700 m. Beyond these thresholds, the contributions of deep bedrock groundwater remained constant, indicating that the REA is between 2 and 3 km² at the observed area. These results indicate that the hydrological conditions of base flow were controlled mainly by the deep bedrock groundwater that moved through thick lava and pyroclastic flows in the undissected volcanic body of the upper part of the catchment. Our study demonstrates that deep and long groundwater movements via a deep bedrock layer including thick deposits of volcanic materials at the two catchments on Daisen volcano strongly determined streamflow discharge instead of the mixing of small‐scale hydrological conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

Multi‐variable and multi‐site calibration and validation of SWAT in a large mountainous catchment with high spatial variability

Many methods developed for calibration and validation of physically based distributed hydrological models are time consuming and computationally intensive. Only a small set of input parameters can be optimized, and the optimization often results in unrealistic values. In this study we adopted a multi‐variable and multi‐site approach to calibration and validation of the Soil Water Assessment Tool (SWAT) model for the Motueka catchment, making use of extensive field measurements. Not only were a number of hydrological processes (model components) in a catchment evaluated, but also a number of subcatchments were used in the calibration. The internal variables used were PET, annual water yield, daily streamflow, baseflow, and soil moisture. The study was conducted using an 11‐year historical flow record (1990–2000); 1990–94 was used for calibration and 1995–2000 for validation. SWAT generally predicted well the PET, water yield and daily streamflow. The predicted daily streamflow matched the observed values, with a Nash–Sutcliffe coefficient of 0·78 during calibration and 0·72 during validation. However, values for subcatchments ranged from 0·31 to 0·67 during calibration, and 0·36 to 0·52 during validation. The predicted soil moisture remained wet compared with the measurement. About 50% of the extra soil water storage predicted by the model can be ascribed to overprediction of precipitation; the remaining 50% discrepancy was likely to be a result of poor representation of soil properties. Hydrological compensations in the modelling results are derived from water balances in the various pathways and storage (evaporation, streamflow, surface runoff, soil moisture and groundwater) and the contributions to streamflow from different geographic areas (hill slopes, variable source areas, sub‐basins, and subcatchments). The use of an integrated multi‐variable and multi‐site method improved the model calibration and validation and highlighted the areas and hydrological processes requiring greater calibration effort. Copyright © 2005 John Wiley & Sons, Ltd.     

Modelling the potential impacts of groundwater hydrology on long‐term drainage basin evolution

Fluvial erosion processes are driven by water discharge on the land surface, which is produced by surface runoff and groundwater discharge. Although groundwater is often neglected in long‐term landscape evolution problems, water table levels control patterns of Dunne runoff production, and groundwater discharge can contribute significantly to storm flows. In this analysis, we investigate the role that groundwater movement plays in long‐term drainage basin evolution by modifying a widely used landscape evolution model to include a more detailed representation of basin hydrology. Precipitation is generated by a stochastic process, and the precipitation is partitioned between surface runoff and groundwater recharge using a specified infiltration capacity. Groundwater flow is simulated by a dynamic two‐dimensional Dupuit equation for an unconfined aquifer with an irregular underlying impervious layer. The model is applied to the WE‐38 basin, an experimental catchment in Pennsylvania, because 60–80 per cent of the discharge is derived from groundwater and substantial hydrologic and geomorphic information is available. The hydrologic model is first calibrated to match the observed streamflows, and then the combined hydrologic/geomorphic model is used to simulate scenarios with different infiltration capacities. The results of this modelling exercise indicate that the basin can be divided into three zones with distinct streamflow‐generating characteristics, and different parts of the basin can have different geomorphic effective events. Over long periods of time, scenarios in which groundwater discharge is large tend to modify the topography in a way that promotes groundwater discharge and inhibits Dunne runoff. Copyright © 2006 John Wiley & Sons, Ltd.     

Separation of base flow from streamflow using groundwater levels—illustrated for the Pang catchment (UK)

A new filter to separate base flow from streamflow has developed that uses observed groundwater levels. To relate the base flow to the observed groundwater levels, a non‐linear relation was used. This relation is suitable for unconfined aquifers with deep groundwater levels that do not respond to individual rainfall event. Because the filter was calibrated using total streamflow, an estimate of the direct runoff was also needed. The direct runoff was estimated from precipitation and potential evapotranspiration using a water balance model. The parameters for the base flow and direct runoff were estimated simultaneously using a Monte Carlo approach. Instead of one best solution, a range of satisfactory solutions was accepted. The filter was applied to data from two nested gauging stations in the Pang catchment (UK). Streamflow at the upstream station (Frilsham) is strongly dominated by base flow from the main aquifer, whereas at the downstream station (Pangbourne) a significant component of direct runoff also occurs. The filter appeared to provide satisfactory estimates at both stations. For Pangbourne, the rise of the base flow was strongly delayed compared with the rise of the streamflow. However, base flow exceeded streamflow on several occasions, especially during summer and autumn, which might be explained by evapotranspiration from riparian vegetation. To evaluate the results, the base flow was also estimated using three existing base‐flow separation filters: an arithmetic filter (BFI), a digital filter (Boughton) and another filter based on groundwater levels (Kliner and Knĕz̆ek). Both the BFI and Boughton filters showed a much smaller difference in base flow between the two stations. The Kliner and Knĕz̆ek filter gave consistently lower estimates of the base flow. Differences and lack of clarity in the definition of base flow complicated the comparison between the filters. An advantage of the method introduced in this paper is the clear interpretation of the separated components. A disadvantage is the high data requirement. Copyright © 2004 John Wiley & Sons, Ltd.     

Simulated response of an intermittent stream to rainfall frequency patterns

The response of intermittent catchments to rainfall is complex and difficult to model. This study uses the spatially distributed CATchment HYdrology (CATHY) model to explore how the frequency of daily rainfall (λ) can affect the hydrologic regime of intermittent catchments. After a multi-objective calibration and validation of CATHY against experimental measurements of streamflow and groundwater levels in a catchment used as a pasture, the role of λ in affecting streamflow characteristics was explored using different scenarios. With different values of λ for the dry and wet periods of the year, CATHY showed that a series of frequent rainfall events was often associated with incipient streamflow, independent of the season. Activation of streamflow during the wet season was related to multiple factors and was not often associated with the shallow groundwater levels near the outlet of the catchment. The interplay between rainfall depth and intensity acted as the most important factor for the generation of streamflow. Using the difference between accumulated rainfall and evapotranspiration as a measure of wetness, saturated subsurface flow mechanism generated streamflow in simulations with wetness at least three times larger than mean wetness of other simulations. Although groundwater uprise near the outlet did not effectively contribute to streamflow in the initial days of flow, it strongly correlated with the magnitude of the runoff coefficient. Values of λ close or equal to the maximum value in the wet season can sustain the connectivity between groundwater and streamflow in the riparian zone. This connectivity increases the catchment wetness, which consequently results in an increase of the generated streamflow. Our study showed that rainfall regimes characterized by different λ were able to identify distinct flow regimes typical of either intermittent, ephemeral, or nonflowing catchments. Decrease of λ in the wet season is likely associated with a reduction of streamflow, with a shift of flow regime from intermittent to ephemeral or no-flow.     

Modification of the SWAT model to simulate regional groundwater flow using a multicell aquifer

The soil and water assessment tool (SWAT) has been widely used and thoroughly tested in many places in the world. The application of the SWAT model has pointed out that 2 of the major weaknesses of SWAT are related to the nonspatial reference of the hydrologic response unit concept and to the simplified groundwater concept, which contribute to its low performance in baseflow simulation and its inability to simulate regional groundwater flow. This study modified the groundwater module of SWAT to overcome the above limitations. The modified groundwater module has 2 aquifers. The local aquifer, which is the shallow aquifer in the original SWAT, represents a local groundwater flow system. The regional aquifer, which replaces the deep aquifer of the original SWAT, represents intermediate and regional groundwater flow systems. Groundwater recharge is partitioned into local and regional aquifer recharges. The regional aquifer is represented by a multicell aquifer (MCA) model. The regional aquifer is discretized into cells using the Thiessen polygon method, where centres of the cells are locations of groundwater observation wells. Groundwater flow between cells is modelled using Darcy's law. Return flow from cell to stream is conceptualized using a non‐linear storage–discharge relationship. The SWAT model with the modified aquifer module, the so‐called SWAT‐MCA, was tested in 2 basins (Wipperau and Neetze) with porous aquifers in a lowland area in Lower Saxony, Germany. Results from the Wipperau basin show that the SWAT‐MCA model is able (a) to simulate baseflow in a lowland area (where baseflow is a dominant source of streamflow) better than the original model and (b) to simulate regional groundwater flow, shown by the simulated groundwater levels in cells, quite well.