Simulating short‐circuiting flow in a constructed wetland: the implications of bathymetry and vegetation effects

Short‐circuiting flow, commonly experienced in many constructed wetlands, reduces hydraulic retention times in unit wetland cells and decreases the treatment efficiency. A two‐dimensional (2‐D), physically based, distributed modelling approach was used to systematically address the effects of bathymetry and vegetation on short‐circuiting flow, which previously have been neglected or lumped in one‐dimensional wetland flow models. In this study, a 2‐D transient hydrodynamics with advection‐dispersion model was developed using MIKE 21 and calibrated with bromide tracer data collected at the Orlando Easterly Wetland Cell 7. The estimated topographic difference between short‐circuiting flow zone and adjacent area ranged from 0·3 to 0·8 m. A range of the Manning roughness coefficient at the short‐circuiting flow zone was estimated (0·022–0·045 s m^?1/3). Sensitivity analysis of topographical and vegetative heterogeneity deduced during model calibration shows that relic ditches or other ditch‐shaped landforms and the associated sparse vegetation along the main flow direction intensify the short‐circuiting pattern, considerably affecting 2‐D solute transport simulation. In terms of hydraulic efficiency, this study indicates that the bathymetry effect on short‐circuiting flow is more important than the vegetation effect. Copyright © 2009 John Wiley & Sons, Ltd.     

Modelling the bulk flow of a bedrock‐constrained,multi‐channel reach of the Mekong River,Siphandone, southern Laos

The general nature of bulk flow within bedrock single‐channel reaches has been considered by several studies recently. However, the flow structure of a bedrock‐constrained, large river with a multiple channel network has not been investigated previously. The multiple channel network of the Siphandone wetlands in Laos, a section of the Mekong River, was modelled using a steady one‐dimensional hydraulic model. The river network is characterized by a spatially‐varying channel‐form leading to significant changes in the bulk flow properties between and along the channels. The challenge to model the bulk flow in such a remote region was the lack of ideal boundary conditions. The flow models considered both low flow, high inbank and overbank flows and were calibrated using SPOT satellite sensor imagery and limited field data concerning water levels. The application of the model highlighted flow characteristics of a large multi‐channel network and also further indicated the field data that would be required to properly characterize the flow field empirically. Important results included the observation that adjacent channels within the network had different water surface slopes for the same moments in time; thus calibration data for modelling similar systems needs to account for these significant local differences. Further, the in‐channel hydraulic roughness coefficient strongly varied from one cross‐section to the next (Manning's ‘n’ range: 0·01 to 0·10). These differences were amplified during low flow but persisted in muted form during high discharges. Copyright © 2012 John Wiley & Sons, Ltd.     

Comparison of AnnAGNPS and SWAT model simulation results in USDA‐CEAP agricultural watersheds in south‐central Kansas

This study was conducted under the USDA‐Conservation Effects Assessment Project (CEAP) in the Cheney Lake watershed in south‐central Kansas. The Cheney Lake watershed has been identified as ‘impaired waters’ under Section 303(d) of the Federal Clean Water Act for sediments and total phosphorus. The USDA‐CEAP seeks to quantify environmental benefits of conservation programmes on water quality by monitoring and modelling. Two of the most widely used USDA watershed‐scale models are Annualized AGricultural Non‐Point Source (AnnAGNPS) and Soil and Water Assessment Tool (SWAT). The objectives of this study were to compare hydrology, sediment, and total phosphorus simulation results from AnnAGNPS and SWAT in separate calibration and validation watersheds. Models were calibrated in Red Rock Creek watershed and validated in Goose Creek watershed, both sub‐watersheds of the Cheney Lake watershed. Forty‐five months (January 1997 to September 2000) of monthly measured flow and water quality data were used to evaluate the two models. Both models generally provided from fair to very good correlation and model efficiency for simulating surface runoff and sediment yield during calibration and validation (correlation coefficient; R², from 0·50 to 0·89, Nash Sutcliffe efficiency index, E, from 0·47 to 0·73, root mean square error, RMSE, from 0·25 to 0·45 m³ s^?1 for flow, from 158 to 312 Mg for sediment yield). Total phosphorus predictions from calibration and validation of SWAT indicated good correlation and model efficiency (R² from 0·60 to 0·70, E from 0·63 to 0·68) while total phosphorus predictions from validation of AnnAGNPS were from unsatisfactory to very good (R² from 0·60 to 0·77, E from ? 2·38 to 0·32). The root mean square error–observations standard deviation ratio (RSR) was estimated as excellent (from 0·08 to 0·25) for the all model simulated parameters during the calibration and validation study. The percentage bias (PBIAS) of the model simulated parameters varied from unsatisfactory to excellent (from 128 to 3). This study determined SWAT to be the most appropriate model for this watershed based on calibration and validation results. Copyright © 2008 John Wiley & Sons, Ltd.     

Hydrology of hillslope‐wetland streams,Polar Bear Pass,Nunavut, Canada

Polar Bear Pass is a large High Arctic low‐gradient wetland (100 km²) bordered by low‐lying hills which are notched by a series of v‐shaped valleys. The spring and summer hydrology of two High Arctic hillslope‐wetland catchments, a first‐order stream, 0·2 km² Landing Strip Creek (LSC) and a larger second‐order basin, 4·2 km² Windy Creek (WC), is described here. A water balance framework was employed in 2008 to examine the movement of water from upland reaches into the low‐lying wetland. Snowcover was low in both basins (<50 mm in water equivalent units), but they both exhibited nival‐type regimes. After the main snowmelt season ended, runoff ceased in the smaller catchment (LSC), but not at the larger basin (WC) which continued to flow throughout the summer. Both basins responded to summer rains in different ways. At LSC, late‐summer continuous streamflow occurred only when rainfall satisfied the large soil moisture deficit in the upper bowl‐shaped zone of the basin. At WC, the presence of thinly thawed, ice‐rich polygonal terrain within the stream channel and in the upper reaches of the catchment likely limited infiltration in these near‐stream zones and enhanced runoff in response to both moderate and high rainfall. Subsequently, seasonal runoff ratios differed between the two sites (0·19 vs 0·68) as did the seasonal storage + residual (+16 vs ?50 mm). This suggests that the post‐snowmelt season runoff response to summer precipitation is very much modified by the unique basin characteristics (soil‐type, vegetation, ground ice) and their location within each stream order type. Copyright © 2010 John Wiley & Sons, Ltd.     

Heat tracing of heterogeneous hyporheic exchange adjacent to in‐stream geomorphic features

There are many field techniques used to quantify rates of hyporheic exchange, which can vary in magnitude and direction spatially over distances of only a few metres, both within and between morphological features. We used in‐stream mini‐piezometers and heat transport modelling of stream and streambed temperatures to quantify the rates and directions of water flux across the streambed interface upstream and downstream of three types of in‐stream geomorphic features: a permanent dam, a beaver dam remnant and a stream meander. We derived hyporheic flux estimates at three different depths at six different sites for a month and then paired those flux rates with measurements of gradient to derive hydraulic conductivity (K) of the streambed sediments. Heat transport modelling provided consistent daily flux estimates that were in agreement directionally with hydraulic gradient measurements and also identified vertical heterogeneities in hydraulic conductivity that led to variable hyporheic exchange. Streambed K varied over an order of magnitude (1·9 × 10^?6 to 5·7 × 10^?5 m/s). Average rates of hyporheic flux ranged from static (q < ±0·02 m/day) to 0·42 m/day. Heat transport modelling results suggest three kinds of flow around the dams and the meander. Copyright © 2010 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.     

Using time‐domain reflectometry to characterize shallow solute transport in an oak woodland hillslope in northern California,USA

The natural heterogeneity of water and solute movement in hillslope soils makes it difficult to accurately characterize the transport of surface‐applied pollutants without first gathering spatially distributed hydrological data. This study examined the application of time‐domain reflectometry (TDR) to measure solute transport in hillslopes. Three different plot designs were used to examine the transport of a conservative tracer in the first 50 cm of a moderately sloping soil. In the first plot, which was designed to examine spatial variability in vertical transport in a 1·2 m² plot, a single probe per meter was found to adequately characterize vertical solute travel times. In addition, a dye and excavation study in this plot revealed lateral preferential flow in small macropores and a transport pattern where solute is focused vertically into preferential flow pathways. The bypass flow delivers solute deeper in the soil, where lateral flow occurs. The second plot, designed to capture both vertical and lateral flow, provided additional evidence confirming the flow patterns identified in the excavation of the first plot. The third plot was designed to examine lateral flow and once again preferential flow of the tracer was observed. In one instance rapid solute transport in this plot was estimated to occur in as little as 3% of the available pore space. Finally, it was demonstrated that the soil anisotropy, although partially responsible for lateral subsurface transport, may also homogenize the transport response across the hillslope by decreasing vertical solute spreading. Copyright © 2002 John Wiley & Sons, Ltd.     

Infiltration and solute transport under a seasonal wetland: bromide tracer experiments in Saskatoon,Canada

In the northern glaciated plain of North America, the duration of surface water in seasonal wetlands is strongly influenced by the rate of infiltration and evaporation. Infiltration also plays important roles in nutrient exchange at the sediment–water interface and groundwater recharge under wetlands. A whole‐wetland bromide tracer experiment was conducted in Saskatchewan, Canada to evaluate infiltration and solute transport processes. Bromide concentrations of surface water, groundwater, sediment pore water and plant tissues were monitored as the pond water‐level gradually dropped until there was no surface water. Hydraulic head gradients showed strong lateral flow from under the wetland to the treed riparian zone during the growing season. The bromide mass balance analysis showed that in early spring, almost 50% of water loss from the wetland was by infiltration, and it increased to about 70% in summer as plants in and around the wetland started to transpire more actively. The infiltration contributed to recharging the shallow, local groundwater under the wetland, but much of it was taken up by trees without recharging the deeper groundwater system. Emergent plants growing in the wetlands incorporated some bromide, but overall uptake of bromide by vegetation was less than 10% of the amount initially released. After one summer, most of the subsurface bromide was found within 40–80 cm of the soil surface. However, some bromide penetrated as deep as 2–3 m, presumably owing to preferential flow pathways provided by root holes or fractures. Copyright © 2004 Crown in the Right of Canada. Published by John Wiley & Sons, Ltd.     

Understanding changes in annual runoff following land use changes: a systematic data‐based approach

A simple modelling framework for assessing the response of ungauged catchments to land use change in South‐Western Australia is presented. The framework uses knowledge of transpiration losses from native vegetation and pasture and then partitions the ‘excess’ water (resulting from reduced transpiration after land use change) between runoff and deep storage. The simple partitioning is achieved by using soft information (satellite imagery, previous mapping and field assessment) to delimit the spread of the permanently saturated area close to the stream. Runoff is then assumed to increase in proportion to the saturated area, with the residual difference going to deep storage. The model parameters to describe the annual water yield are obtained a priori and no calibration to streamflow is required. We tested the model using gauged records over 25 years from paired catchment experiments in South‐Western Australia. Very good estimates of runoff were obtained from high rainfall (>1100 mm yr⁻¹) catchments (R² > 0·9) and for low rainfall (<900 mm yr⁻¹) catchments after clearing (R² = 0·96) but results were poorer (R² = 0·55) for an uncleared low rainfall catchment, although overall balances were reasonable. In the drier uncleared catchments, the within‐year distributions of rainfall may exert a substantial influence on runoff response that is not completely captured by the presented model. Copyright © 2005 John Wiley & Sons, Ltd.     

Wetland nitrogen dynamics in an Adirondack forested watershed

Wetlands often form the transition zone between upland soils and watershed streams, however, stream–wetland interactions and hydrobiogeochemical processes are poorly understood. We measured changes in stream nitrogen (N) through one riparian wetland and one beaver meadow in the Archer Creek watershed in the Adirondack Mountains of New York State, USA from 1 March to 31 July 1996. In the riparian wetland we also measured changes in groundwater N. Groundwater N changed significantly from tension lysimeters at the edge of the peatland to piezometer nests within the peatland. Mean N concentrations at the peatland perimeter were 1·5, 0·5 and 18·6 µmol L^?1 for NH₄⁺, NO₃^? and DON (dissolved organic nitrogen), respectively, whereas peatland groundwater N concentration was 56·9, 1·5 and 31·6 µmol L^?1 for NH₄⁺, NO₃^? and DON, respectively. The mean concentrations of stream water N species at the inlet to the wetlands were 1·5, 10·1 and 16·9 µmol L^?1 for NH₄⁺, NO₃^? and DON, respectively and 1·6, 28·1 and 8·4 µmol L^?1 at the wetland outlet. Although groundwater total dissolved N (TDN) concentrations changed more than stream water TDN through the wetlands, hydrological cross‐sections for the peatland showed that wetland groundwater contributed minimally to stream flow during the study period. Therefore, surface water N chemistry was affected more by in‐stream N transformations than by groundwater N transformations because the in‐stream changes, although small, affected a much greater volume of water. Stream water N input–output budgets indicated that the riparian peatland retained 0·16 mol N ha^?1 day^?1 of total dissolved N and the beaver meadow retained 0·26 mol N ha^?1 day^?1 during the study period. Nitrate dominated surface water TDN flux from the wetlands during the spring whereas DON dominated during the summer. This study demonstrates that although groundwater N changed significantly in the riparian peatland, those changes were not reflected in the stream. Consequently, although in‐stream changes of N concentrations were less marked than those in groundwater, they had a greater effect on stream water chemistry—because wetland groundwater contributed minimally to stream flow. Copyright © 2004 John Wiley & Sons, Ltd.     

Modelling Canadian prairie wetland hydrology using a semi‐distributed streamflow model

A hydrological model (SLURP) that was designed for simulating hydrological processes taking place in large river basins was, with minimal modification, used successfully to simulate water level variations over a 28‐year period (1969–1996) for a 3‐ha prairie wetland in Saskatchewan. The model calculates a water balance based on precipitation, snowmelt, evaporation, surface runoff and subsurface flow on a daily time‐step. The model was first calibrated for two periods (1969–1973 for cropland and 1987–1990 for grassland), then it was applied to records outside the calibration periods. The model reproduced the wetland water level variations during a 28‐year period with good accuracy. The wetland water levels were most sensitive to the infiltration coefficient of surface soil under frozen conditions and to maximum soil moisture storage. The applicability of the model and the calibrated parameters to a smaller wetland, with an area of 0·24 ha, was examined. This simulation indicated that scale effects are important, probably largely in relation to snow redistribution by wind. Copyright © 2000 John Wiley & Sons, Ltd.     

Soil partitioning and surface store controls on spring runoff from a boreal forest peatland basin in north‐central Manitoba,Canada

The mid‐ to high‐boreal forest in Canada occupies the discontinuous permafrost zone, and is often underlain by glaciolacustrine sediments mantled by a highly porous organic mat. The result is a poorly drained landscape dominated by wetlands. Frost‐table dynamics and surface storage conditions help to control runoff contributions from various landscape elements, hydrological linkages between these elements, and basin streamflow during spring snowmelt. Runoff components and pathways in a forested peatland basin were assessed during two spring snowmelts with contrasting input and basin conditions. Runoff from relatively intense melt (up to 16 mm day^?1) on slopes with limited soil thawing combined with large pre‐melt storage in surface depressions to produce high flows composed primarily of meltwater (78% of the 0·29 m³ s^?1 peak discharge) routed over wetland surfaces and through permeable upper peat layers. Melt intensity was less in the subsequent year (maximum of 10 mm day^?1) and active layer development was relatively greater (0·2 m deeper at the end of spring melt), resulting in less slope runoff. Coupling of reduced slope contributions with lower storage levels in basin wetlands led to relatively subdued streamflows dominated by older water (73% of the 0·09 m³ s^?1 peak discharge) routed through less‐permeable deeper peat layers and mineral soil. Interannual differences in runoff conditions provide important insight for the development of distributed hydrological models for boreal forest basins and into potential influences on biogeochemical cycling in this landscape under a warming climate. Copyright © 2001 John Wiley & Sons, Ltd.     

Distributed and integrated response of a geographic information system‐based hydrologic model in the eastern Palouse region,Idaho

Verification of distributed hydrologic models is rare owing to the lack of spatially detailed field measurements and a common mismatch between the scale at which soil hydraulic properties are measured and the scale of a single modelling unit. In this study, two of the most commonly calibrated parameters, i.e. soil depth and the vertical distribution of lateral saturated hydraulic conductivity K_s, were eliminated by a spatially detailed soil characterization and results of a hillslope‐scale field experiment. The soil moisture routing (SMR) model, a geographic information system‐based hydrologic model, was modified to represent the dominant hydrologic processes for the Palouse region of northern Idaho. Modifications included K_s as a double exponential function of depth in a single soil layer, a snow accumulation and melt algorithm, and a simple relationship between storage and perched water depth (PWD) using the drainable porosity. The model was applied to a 2 ha catchment without calibration to measured data. Distributed responses were compared with observed PWD over a 3‐year period on a 10 m × 15 m grid. Integrated responses were compared with observed surface runoff at the catchment outlet. The modified SMR model simulated the PWD fluctuations remarkably well, especially considering the shallow soils in this catchment: a 0·20 m error in PWD is equivalent to only a 1·6% error in predicted soil moisture content. Simulations also captured PWD fluctuations during a year with high spatial variability of snow accumulation and snowmelt rates at upslope, mid‐slope, and toe slope positions with errors as low as 0·09 m, 0·12 m, and 0·12 m respectively. Errors in distributed and integrated model simulations were attributed mostly to misrepresentation of rain events and snowmelt timing problems. In one location in the catchment, simulated PWD was consistently greater than observed PWD, indicating a localized recharge zone, which was not identified by the soil morphological survey. Copyright © 2006 John Wiley & Sons, Ltd.     

Runoff and sediment generation on bench‐terraced hillsides: measurements and up‐scaling of a field‐based model

Despite widespread bench‐terracing, stream sediment yields from agricultural hillsides in upland West Java remain high. We studied the causes of this lack of effect by combining measurements at different spatial scales using an erosion process model. Event runoff and sediment yield from two 4‐ha terraced hillside subcatchments were measured and field surveys of land use, bench‐terrace geometry and storage of sediment in the drainage network were conducted for two consecutive years. Runoff was 3·0–3·9% of rainfall and sediment yield was 11–30 t ha⁻¹ yr⁻¹ for different years, subcatchments and calculation techniques. Sediment storage changes in the subcatchment drainage network were less than 2 t ha⁻¹, whereas an additional 0·3–1·5 t ha⁻¹ was stored in the gully between the subcatchment flumes and the main stream. This suggests mean annual sediment delivery ratios of 86–125%, or 80–104% if this additional storage is included. The Terrace Erosion and Sediment Transport (TEST) model developed and validated for the studied environment was parameterized using erosion plot studies, land use surveys and digital terrain analysis to simulate runoff and sediment generation on the terraced hillsides. This resulted in over‐estimates of runoff and under‐estimates of runoff sediment concentration. Relatively poor model performance was attributed to sample bias in the six erosion plots used for model calibration and unaccounted covariance between important terrain attributes such as slope, infiltration capacity, soil conservation works and vegetation cover. Copyright © 2005 John Wiley & Sons, Ltd.     

Evaluation of a dynamic multi‐class sediment transport model in a catchment under soil‐conservation agriculture

Soil erosion models are essential tools for the successful implementation of effective and adapted soil conservation measures on agricultural land. Therefore, models are needed that predict sediment delivery and quality, give a good spatial representation of erosion and deposition and allow us to account for various soil conservation measures. Here, we evaluate how well a modified version of the spatially distributed multi‐class sediment transport model (MCST) simulates the effectiveness of control measures for different event sizes. We use 8 year runoff and sediment delivery data from two small agricultural watersheds (0·7 and 3·7 ha) under optimized soil conservation. The modified MCST model successfully simulates surface runoff and sediment delivery from both watersheds; one of which was dominated by sheet and the other was partly affected by rill erosion. Moreover, first results of modelling enrichment of clay in sediment delivery are promising, showing the potential of MCST to model sediment enrichment and nutrient transport. In general, our results and those of an earlier modelling exercise in the Belgian Loess Belt indicate the potential of the MCST model to evaluate soil erosion and deposition under different agricultural land uses. As the model explicitly takes into account the dominant effects of soil‐conservation agriculture, it should be successfully applicable for soil‐conservation planning/evaluation in other environments. Copyright © 2008 John Wiley & Sons, Ltd.     

Climatic control on river discharge simulations,Zackenberg River drainage basin,northeast Greenland

A spatially distributed, physically based, hydrologic modeling system (MIKE SHE) was applied to quantify intra‐ and inter‐annual discharge from the snow and glacierized Zackenberg River drainage basin (512 km²; 20% glacier cover) in northeast Greenland. Evolution of snow accumulation, distribution by wind‐blown snow, blowing‐snow sublimation, and snow and ice surface melt were simulated by a spatially distributed, physically based, snow‐evolution modelling system (SnowModel) and used as input to MIKE SHE. Discharge simulations were performed for three periods 1997–2001 (calibration period), 2001–2005 (validation period), and 2071–2100 (scenario period). The combination of SnowModel and MIKE SHE shows promising results; the timing and magnitude of simulated discharge were generally in accordance with observations (R² = 0·58); however, discrepancies between simulated and observed discharge hydrographs do occur (maximum daily difference up to 44·6 m³ s^?1 and up to 9% difference between observed and simulated cumulative discharge). The model does not perform well when a sudden outburst of glacial dammed water occurs, like the 2005 extreme flood event. The modelling study showed that soil processes related to yearly change in active layer depth and glacial processes (such as changes in yearly glacier area, seasonal changes in the internal glacier drainage system, and the sudden release of glacial bulk water storage) need to be determined, for example, from field studies and incorporated in the models before basin runoff can be quantified more precisely. The SnowModel and MIKE SHE model only include first‐order effects of climate change. For the period 2071–2100, future IPCC A2 and B2 climate scenarios based on the HIRHAM regional climate model and HadCM3 atmosphere–ocean general circulation model simulations indicated a mean annual Zackenberg runoff about 1·5 orders of magnitude greater (around 650 mmWE year^?1) than from today 1997–2005 (around 430 mmWE year^?1), mainly based on changes in negative glacier net mass balance. Copyright © 2007 John Wiley & Sons, Ltd.     

A temperature‐precipitation‐based model of thirty‐year mean snowpack accumulation and melt in Oregon,USA

High‐resolution, spatially extensive climate grids can be useful in regional hydrologic applications. However, in regions where precipitation is dominated by snow, snowmelt models are often used to account for timing and magnitude of water delivery. We developed an empirical, nonlinear model to estimate 30‐year means of monthly snowpack and snowmelt throughout Oregon. Precipitation and temperature for the period 1971–2000, derived from 400‐m resolution PRISM data, and potential evapotranspiration (estimated from temperature and day length) drive the model. The model was calibrated using mean monthly data from 45 SNOTEL sites and accurately estimated snowpack at 25 validation sites: R² = 0·76, Nash‐Sutcliffe Efficiency (NSE) = 0·80. Calibrating it with data from all 70 SNOTEL sites gave somewhat better results (R² = 0·84, NSE = 0·85). We separately applied the model to SNOTEL stations located < 200 and ≥ 200 km from the Oregon coast, since they have different climatic conditions. The model performed equally well for both areas. We used the model to modify moisture surplus (precipitation minus potential evapotranspiration) to account for snowpack accumulation and snowmelt. The resulting values accurately reflect the shape and magnitude of runoff at a snow‐dominated basin, with low winter values and a June peak. Our findings suggest that the model is robust with respect to different climatic conditions, and that it can be used to estimate potential runoff in snow‐dominated basins. The model may allow high‐resolution, regional hydrologic comparisons to be made across basins that are differentially affected by snowpack, and may prove useful for investigating regional hydrologic response to climate change. Published in 2011 by John Wiley & Sons, Ltd.     

Parameter estimation in semi‐distributed hydrological catchment modelling using a multi‐criteria objective function

Output generated by hydrologic simulation models is traditionally calibrated and validated using split‐samples of observed time series of total water flow, measured at the drainage outlet of the river basin. Although this approach might yield an optimal set of model parameters, capable of reproducing the total flow, it has been observed that the flow components making up the total flow are often poorly reproduced. Previous research suggests that notwithstanding the underlying physical processes are often poorly mimicked through calibration of a set of parameters hydrologic models most of the time acceptably estimates the total flow. The objective of this study was to calibrate and validate a computer‐based hydrologic model with respect to the total and slow flow. The quick flow component used in this study was taken as the difference between the total and slow flow. Model calibrations were pursued on the basis of comparing the simulated output with the observed total and slow flow using qualitative (graphical) assessments and quantitative (statistical) indicators. The study was conducted using the Soil and Water Assessment Tool (SWAT) model and a 10‐year historical record (1986–1995) of the daily flow components of the Grote Nete River basin (Belgium). The data of the period 1986–1989 were used for model calibration and data of the period 1990–1995 for model validation. The predicted daily average total flow matched the observed values with a Nash–Sutcliff coefficient of 0·67 during calibration and 0·66 during validation. The Nash–Sutcliff coefficient for slow flow was 0·72 during calibration and 0·61 during validation. Analysis of high and low flows indicated that the model is unbiased. A sensitivity analysis revealed that for the modelling of the daily total flow, accurate estimation of all 10 calibration parameters in the SWAT model is justified, while for the slow flow processes only 4 out of the set of 10 parameters were identified as most sensitive. Copyright © 2007 John Wiley & Sons, Ltd.     

SWAT‐simulated hydrological impact of land‐use change in the Zanjanrood basin,Northwest Iran

Understanding the impacts of land‐use changes on hydrology at the watershed scale can facilitate development of sustainable water resource strategies. This paper investigates the hydrological effects of land‐use change in Zanjanrood basin, Iran. The water balance was simulated using the Soil and Water Assessment Tool (AVSWAT2000). Model calibration and uncertainty analysis were performed with sequential uncertainty fitting (SUFI‐2). Simulation results from January 1998 to December 2002 were used for parameter calibration, and then the model was validated for the period of January 2003 to December 2004. The predicted monthly streamflow matched the observed values: during calibration the correlation coefficient was 0·86 and the Nash–Sutcliffe coefficient 0·79, compared with 0·80 and 0·79, respectively, during validation. The model was used to simulate the main components of the hydrological cycle, in order to study the effects of land‐use changes in 1967, 1994 and 2007. The study reveals that during 1967 a 34·5% decrease of grassland with concurrent increases of shrubland (13·9%), rain‐fed agriculture (12·1%), bare ground (5·5%) irrigated agriculture (2·2%), and urban area (0·7%) led to a 33% increase in the amount of surface runoff and a 22% decrease in the groundwater recharge. Furthermore, the area of sub‐basins that was influenced by high runoff (14–28 mm) increased. The results indicate that the hydrological response to overgrazing and the replacing of rangelands (grassland and shrubland) with rain‐fed agriculture and bare ground (badlands) is nonlinear and exhibits a threshold effect. The runoff rises dramatically when more than 60% of the rangeland is removed. For groundwater this threshold lies at an 80% decrease in rangeland. Copyright © 2009 John Wiley & Sons, Ltd.     

Wetlands with controlled drainage and sub‐irrigation systems—modelling of the water balance

Over the past centuries, the agricultural use of wetlands in Central Europe has required interference with the natural wetland water balance. Often this has consisted of drainage measures alone. In low‐precipitation areas, it has also involved the operation of combined drainage and sub‐irrigation systems. Model studies conducted as part of planning processes, or with a view to finding out the impact of changing climate conditions on the water balance of wetlands, must take these facts into account. For this reason, a water balance model has been devised for wetlands whose water balance is governed by water resources management systems. It is based on the WBalMo model system. Special modules were integrated into WBalMo to calculate the water balance of wetland areas (WABI module) and to regulate inflow partitioning within the wetland (REGINF module). When calculating the water balance, the WABI module takes into account precipitation and potential evapotranspiration, groundwater levels below surface, soil types, land‐use classes, inflows via the running water system, and data for target water levels. It provides actual evapotranspiration, discharge into the running water system, and groundwater levels in the area. The example of the Spreewald, a major wetland area in north‐eastern Germany, was used to design and test the WBalMo Spreewald model. The comparison of measured and calculated water balance parameters of the wetland area confirms the suitability of the model for water balance studies in wetlands with complex water resources management systems. The results reveal the strong influence of water management on the water balance of such areas. The model system has proved to be excellently suited for planning and carrying out water management measures aimed at the sustainable development of wetlands. Furthermore, scenario analyses can be used to assess the impact of global change on the water balance of wetlands. Copyright © 2006 John Wiley & Sons, Ltd.