Simulated soil water storage effects on streamflow generation in a mountainous snowmelt environment,Idaho, USA

Although soil processes affect the timing and amount of streamflow generated from snowmelt, they are often overlooked in estimations of snowmelt‐generated streamflow in the western USA. The use of a soil water balance modelling approach to incorporate the effects of soil processes, in particular soil water storage, on the timing and amount of snowmelt generated streamflow, was investigated. The study was conducted in the Reynolds Mountain East (RME) watershed, a 38 ha, snowmelt‐dominated watershed in southwest Idaho. Snowmelt or rainfall inputs to the soil were determined using a well established snow accumulation and melt model (Isnobal). The soil water balance model was first evaluated at a point scale, using periodic soil water content measurements made over two years at 14 sites. In general, the simulated soil water profiles were in agreement with measurements (P < 0·05) as further indicated by high R² values (mostly > 0·85), y‐intercept values near 0, slopes near 1 and low average differences between measured and modelled values. In addition, observed soil water dynamics were generally consistent with critical model assumptions. Spatially distributed simulations over the watershed for the same two years indicate that streamflow initiation and cessation are closely linked to the overall watershed soil water storage capacity, which acts as a threshold. When soil water storage was below the threshold, streamflow was insensitive to snowmelt inputs, but once the threshold was crossed, the streamflow response was very rapid. At these times there was a relatively high degree of spatial continuity of satiated soils within the watershed. Incorporation of soil water storage effects may improve estimation of the timing and amount of streamflow generated from mountainous watersheds dominated by snowmelt. Copyright © 2008 John Wiley & Sons, Ltd.     

A simple concept for calibrating runoff thresholds in quasi‐distributed variable source area watershed models

Most semi‐distributed watershed water quality models divide the watershed into hydrologic response units (HRU) with no flow among them. This is problematic when watersheds are delineated to include variable source areas (VSAs) because it is the lateral flows from upslope areas to downslope areas that generate VSAs. Although hydrologic modellers have often successfully calibrated these types of models, there can still be considerable uncertainty in model results. In this paper, a topographic‐index‐based method is described and tested to distribute effective soil water holding capacity among HRUs, which can be subsequently adjusted using the watershed baseflow coefficient. The method is tested using a version of the Soil and Water Assessment Tool (SWAT) model that simulates VSA runoff and is applied to two watersheds: a New York State (NYS) watershed, and one in the head waters of the Blue Nile Basin (BNB) in Ethiopia. Daily streamflow predicted using effective soil water storage capacities based only on the topographic index were reassuringly accurate in both the NYS watershed (daily Nash Sutcliffe (E) = 0·73) and in the BNB (E = 0·70). Using the baseflow coefficient to adjust the effective soil water storage capacity only slightly improved streamflow predictions in NYS (E = 0·75) but substantially improved the BNB predictions (E = 0·80). By comparison, the standard SWAT model, which uses the traditional look‐up tables to determine a runoff curve number, performed considerably less accurately in un‐calibrated form (E = 0·51 for NYS and E = 0·45 for BNB), but improved substantially when explicitly calibrated to streamflow measurements (E = 0·76 for NYS and E = 0·67 for the BNB). The calibration method presented here provides a parsimonious, systematic approach to using established models in VSA watersheds that reduces the ambiguity inherent in model calibration. Copyright © 2011 John Wiley & Sons, Ltd.     

Modelling potential hydrological impact of abandoned underground mines in the Monday Creek Watershed,Ohio

Abandoned underground mines (AUM) have caused dramatic environmental effects that are closely linked to regional sustainability. This paper explores the potential hydrological impact of AUM in the Monday Creek Watershed, a typically mined area in Appalachian region, using the Soil and Water Assessment Tool (SWAT 2005) model and Sequential Uncertainty Fitting (SUFI‐2), calibrated at both the global and local scales. The locally calibrated model better incorporates those key parameters relevant to AUM for specific sub‐basins and hydrologic response units. Data from the years 2003–2004 were used for calibration and 2005–2006 for validation. The results were quite satisfactory; both the coefficient of determination (R²) and the Nash–Sutcliffe efficiency statistic were over 0.80. The potential influences of AUM were assessed by modelling an alternative scenario assuming no AUM for the period 2003–2009. Results show that the hydrological process of lateral subsurface flow plays a dominant role in linking AUM to overall watershed hydrology. The potential hydrological impact of AUM is an increased annual lateral flow of 82.1%, and a decrease in annual surface flow by 15%, leading to an increase of 16.9% in annual water yield for the Monday Creek Watershed. The seasonal fluctuation of water yield has a similar trend to lateral flow, decreasing from March to August and increasing from August to January. Higher volume, higher flow peaks and higher recession constants characterized the hydrograph of daily streamflow from AUM. The results indicate that more emphasis should be put on lateral flow for further study of acid mine drainage and flooding control in those watersheds with AUM. Copyright © 2012 John Wiley & Sons, Ltd.     

Spatial patterns of soil nitrate in Japanese forested watersheds: importance of the near-stream zone as a source of nitrate in stream water

Spatial patterns of N dynamics in soil were evaluated within two small forested watersheds in Japan. These two watersheds were characterized by steep slopes (>30°) and high stream NO⁻₃ drainage rates (8·4 to 25·1 kg N ha⁻¹ yr⁻¹) that were greater than bulk precipitation N input rates (7·5 to 13·5 kg N ha⁻¹ yr⁻¹). Higher rates of nitrification potential at near-stream zones were reflected in greater NO⁻₃ contents for soil at the near-stream zones compared with ridge zones. Both stream discharge rates and NO⁻₃ concentrations in deep unsaturated soil at the near-stream zones were positively correlated to NO⁻₃ concentrations in stream water. These relationships, together with high soil NO⁻₃ contents at the near-stream zones, suggest that the near-stream zone was an important source of NO⁻₃ to stream water. Nitrate flux from these near-stream zones was also related to the drainage of cations (K⁺, Ca²⁺ and Mg²⁺). The steep slope of the watersheds resulted in small saturated areas that contributed to the high NO⁻₃ production (high nitrification rates) in the near-stream zone. © 1998 John Wiley & Sons, Ltd.     

Controls of surface soil moisture spatial patterns and their temporal stability in a semi‐arid steppe

Temporal stability of soil moisture spatial patterns has important implications for optimal soil and water management and effective field monitoring. The aim of this study was to investigate the temporal stability of soil moisture spatial patterns over four plots of 105 m × 135 m in grid size with different grazing intensities in a semi‐arid steppe in China. We also examined whether a time‐stable location can be identified from causative factors (i.e. soil, vegetation, and topography). At each plot, surface soil moisture (0–6 cm) was measured about biweekly from 2004 to 2006 using 100 points in each grid. Possible controls of soil moisture, including soil texture, organic carbon, bulk density, vegetation coverage, and topographic indices, were determined at the same grid points. The results showed that the spatial patterns of soil moisture were considerably stable over the 3‐y monitoring period. Soil moisture under wet conditions (averaged volumetric moisture contents > 20%) was more stable than that under dry ( ) or moist ( ) conditions. The best representative point for the whole field identified in each plot was accurate in representing the field mean moisture over time (R² ≥ 0·97; p < 0·0001). The degree of temporal persistence varied with grazing intensity, which was partly related to grazing‐induced differences in soil and vegetation properties. The correlation analysis showed that soil properties, and to a lesser extent vegetation and topographic properties, were important in controlling the temporal stability of soil moisture spatial patterns in this relatively flat grassland. Response surface regression analysis was used to quantitatively identify representative monitoring locations a priori from available soil‐plant parameters. This allows appropriate selection of monitoring locations and enhances efficiency in managing soil and water resources in semi‐arid environments. Copyright © 2010 John Wiley & Sons, Ltd.     

Using GIS and a digital elevation model to assess the effectiveness of variable grade flow diversion terraces in reducing soil erosion in northwestern New Brunswick,Canada

Flow diversion terraces (FDT) are commonly used beneficial management practice (BMP) for soil conservation on sloped terrain susceptible to water erosion. A simple GIS‐based soil erosion model was designed to assess the effectiveness of the FDT system under different climatic, topographic, and soil conditions at a sub‐basin level. The model was used to estimate the soil conservation support practice factor (P‐factor), which inherently considered two major outcomes with its implementation, namely (1) reduced slope length, and (2) sediment deposition in terraced channels. A benchmark site, the agriculture‐dominated watershed in northwestern New Brunswick (NB), was selected to test the performance of the model and estimated P‐factors. The estimated P‐factors ranged from 0·38–1·0 for soil conservation planning objectives and ranged from 0·001 to 0·45 in sediment yield calculations for water‐quality assessment. The model estimated that the average annual sediment yield was 773 kg ha^?1 yr ^?1 compared with a measured value of 641 kg ha^?1 yr^?1. The P‐factors estimated in this study were comparable with predicted values obtained with the revised universal soil loss equation (RUSLE2). The P‐factors from this study have the potential to be directly used as input in hydrological models, such as the soil and water assessment tool (SWAT), or in soil conservation planning where only conventional digital elevation models (DEMs) are available. Copyright © 2009 John Wiley & Sons, Ltd.     

Estimation of tree water stress from stem and soil water monitoring with time‐domain reflectometry in two small forested basins in Spain

Soil‐tree water relationships were studied using time domain reflectometry (TDR) in two small forested basins in Spain. The stem water content of two Mediterranean Quercus species (Quercus pyrenaica and Quercus rotundifolia) was measured using previously constructed species‐specific equations. To monitor soil moisture, a TDR station network was used in both cases. Sixteen Q. pyrenaica and six Q. rotundifolia individuals were selected to install two TDR probes in their trunks (at 20 and 120 cm above the ground) to monitor stem water content. Stem and soil water contents were measured fortnightly. The stem water content of both species showed a similar temporal trend for the period studied. A spring maximum (0·654 cm³ cm^?3 for Q. pyrenaica and 0·568 cm³ cm^?3 for Q. rotundifolia) was found to be associated with high transpiration and no soil moisture deficit, and a late‐summer minimum (0·520 cm³ cm^?3 for Q. pyrenaica and 0·426 cm³ cm^?3 for Q. rotundifolia) was associated with the end of the dry season. This drop in stem water content occurs when the available water in the soil decreases. This seasonal difference presumably reflects water withdrawn from stem storage to support the transpirational demands of the tree. Since plant water stress results in reduced stem water content and since this drop can be measured by TDR, it may be concluded that this technology offers a suitable tool for detecting plant water stress. Copyright © 2007 John Wiley & Sons, Ltd.     

TDR pressure cell for monitoring water content retention and bulk electrical conductivity curves in undisturbed soil samples

The water retention curve (θ(ψ)), which defines the relationship between soil volumetric water content (θ) and matric potential (ψ), is of paramount importance in characterizing the hydraulic behaviour of soils. However, few methods are so far available for estimating θ(ψ) in undisturbed soil samples. We present a new design of TDR‐pressure cell (TDR‐Cell) for estimating θ(ψ) in undisturbed soil samples. The TDR‐Cell consists of a 50‐mm‐long and 50‐mm internal diameter stainless steel cylinder (which constitutes the outer frame of a coaxial line) attached to a porous ceramic disc and closed at the ends with two aluminium lids. A 49‐mm‐long and 3‐mm‐diameter stainless steel rod, which runs longitudinally through the centre of the cylinder, constitutes the inner rod of a coaxial TDR probe. The TDR‐Cell was used to determine the θ(ψ) curves of a packed sand and seven undisturbed soil samples from three profiles of agricultural soils. These θ(ψ) curves were subsequently compared to those obtained from the corresponding 2‐mm sieved soils using the pressure plate method. Measurements of bulk electrical conductivity, σ_a, as a function of the water content, σ_a(θ), of the undisturbed soil samples were also performed. An excellent correlation (R² = 0·988) was found between the θ values measured by TDR on the different undisturbed soils and the corresponding θ obtained from the soil gravimetric water content. A typical bimodal θ(ψ) function was found for most of the undisturbed soil samples. Comparison between the θ(ψ) curves measured with the TDR‐Cell and those obtained from the 2‐mm sieved soils showed that the pressure plate method overestimates θ at low ψ values. The σ_a(θ) relationship was well described by a simple power expression (R² > 0·95), in which the power factor, defined as tortuosity, ranged between 1·18 and 3·75. Copyright © 2011 John Wiley & Sons, Ltd.     

Development and application of a physically based landscape water balance in the SWAT model

Watershed scale hydrological and biogeochemical models rely on the correct spatial‐temporal prediction of processes governing water and contaminant movement. The Soil and Water Assessment Tool (SWAT) model, one of the most commonly used watershed scale models, uses the popular curve number (CN) method to determine the respective amounts of infiltration and surface runoff. Although appropriate for flood forecasting in temperate climates, the CN method has been shown to be less than ideal in many situations (e.g. monsoonal climates and areas dominated by variable source area hydrology). The CN model is based on the assumption that there is a unique relationship between the average moisture content and the CN for all hydrologic response units (HRUs), and that the moisture content distribution is similar for each runoff event, which is not the case in many regions. Presented here is a physically based water balance that was coded in the SWAT model to replace the CN method of runoff generation. To compare this new water balance SWAT (SWAT‐WB) to the original CN‐based SWAT (SWAT‐CN), two watersheds were initialized; one in the headwaters of the Blue Nile in Ethiopia and one in the Catskill Mountains of New York. In the Ethiopian watershed, streamflow predictions were better using SWAT‐WB than SWAT‐CN [Nash–Sutcliffe efficiencies (NSE) of 0·79 and 0·67, respectively]. In the temperate Catskills, SWAT‐WB and SWAT‐CN predictions were approximately equivalent (NSE > 0·70). The spatial distribution of runoff‐generating areas differed greatly between the two models, with SWAT‐WB reflecting the topographical controls imposed on the model. Results show that a water balance provides results equal to or better than the CN, but with a more physically based approach. Copyright © 2010 John Wiley & Sons, Ltd.     

Vulnerability of water resources,vegetation productivity and soil erosion to climate change in Mediterranean watersheds

Climate change is expected to increase temperatures and lower rainfall in Mediterranean regions; however, there is a great degree of uncertainty as to the amount of change. This limits the prediction capacity of models to quantify impacts on water resources, vegetation productivity and erosion. This work circumvents this problem by analysing the sensitivity of these variables to varying degrees of temperature change (increased by up to 6·4 °C), rainfall (reduced by up to 40%) and atmospheric CO₂ concentrations (increased by up to 100%). The SWAT watershed model was applied to 18 large watersheds in two contrasting regions of Portugal, one humid and one semi‐arid; incremental changes to climate variables were simulated using a stochastic weather generator. The main results indicate that water runoff, particularly subsurface runoff, is highly sensitive to these climate change trends (down by 80%). The biomass growth of most species showed a declining trend (wheat down by 40%), due to the negative impacts of increasing temperatures, dampened by higher CO₂ concentrations. Mediterranean species, however, showed a positive response to milder degrees of climate change. Changes to erosion depended on the interactions between the decline in surface runoff (driving erosion rates downward) and biomass growth (driving erosion rates upward). For the milder rainfall changes, soil erosion showed a significant increasing trend in wheat fields (up to 150% in the humid watersheds), well above the recovery capacity of the soil. Overall, the results indicate a shift of the humid watersheds to acquire semi‐arid characteristics, such as more irregular river flows and increasingly marginal conditions for agricultural production. Copyright © 2007 John Wiley & Sons, Ltd.     

Impact of soil and water conservation measures on catchment hydrological response—a case in north Ethiopia

Impact studies of catchment management in the developing world rarely include detailed hydrological components. Here, changes in the hydrological response of a 200‐ha catchment in north Ethiopia are investigated. The management included various soil and water conservation measures such as the construction of dry masonry stone bunds and check dams, the abandonment of post‐harvest grazing, and the establishment of woody vegetation. Measurements at the catchment outlet indicated a runoff depth of 5 mm or a runoff coefficient (RC) of 1·6% in the rainy season of 2006. Combined with runoff measurements at plot scale, this allowed calculating the runoff curve number (CN) for various land uses and land management techniques. The pre‐implementation runoff depth was then predicted using the CN values and a ponding adjustment factor, representing the abstraction of runoff induced by the 242 check dams in gullies. Using the 2006 rainfall depths, the runoff depth for the 2000 land management situation was predicted to be 26·5 mm (RC = 8%), in line with current RCs of nearby catchments. Monitoring of the ground water level indicated a rise after catchment management. The yearly rise in water table after the onset of the rains (ΔT) relative to the water surplus (WS) over the same period increased between 2002–2003 (ΔT/WS = 3·4) and 2006 (ΔT/WS >11·1). Emerging wells and irrigation are other indicators for improved water supply in the managed catchment. Cropped fields in the gullies indicate that farmers are less frightened for the destructive effects of flash floods. Due to increased soil water content, the crop growing period is prolonged. It can be concluded that this catchment management has resulted in a higher infiltration rate and a reduction of direct runoff volume by 81% which has had a positive influence on the catchment water balance. Copyright © 2010 John Wiley & Sons, Ltd.     

Synergistic effects of water repellency and macropore flow on the hydraulic conductivity of a burned forest soil,south‐east Australia

Research shows that water repellency is a key hydraulic property that results in reduced infiltration rates in burned soils. However, more work is required in order to link the hydrological behaviour of water repellent soils to observed runoff responses at the plot and hillslope scale. This study used 5 M ethanol and water in disc infiltrometers to quantify the role of macropore flow and water repellency on spatial and temporal infiltration patterns in a burned soil at plot (<10 m²) scale in a wet eucalypt forest in south‐east Australia. In the first summer and winter after wildfire, an average of 70% and 60%, respectively, of the plot area was water repellent and did not contribute to infiltration. Macropores (r > 0·5 mm), comprising just 5·5% of the soil volume, contributed to 70% and 95%, respectively, of the field‐saturated and ponded hydraulic conductivity (K_p). Because flow occurred almost entirely via macropores in non‐repellent areas, this meant that less than 2·5% of the soil surface effectively contributed to infiltration. The hydraulic conductivity increased by a factor of up to 2·5 as the hydraulic head increased from 0 to 5 mm. Due to the synergistic effect of macropore flow and water repellency, the coefficient of variation (CV) in K_p was three times higher in the water‐repellent soil (CV = 175%) than under the simulated non‐repellent conditions (CV = 66%). The high spatial variability in K_p would act to reduce the effective infiltration rate during runoff generation at plot scale. Ponding, which tend to increase with increasing scale, activates flow through macropores and would raise the effective infiltration rates at larger scales. Field experiments designed to provide representative measurements of infiltration after fire in these systems must therefore consider both the inherent variability in hydraulic conductivity and the variability in infiltration caused by interactions between surface runoff and hydraulic conductivity. Copyright © 2010 John Wiley & Sons, Ltd.     

Daily streamflow modelling and assessment based on the curve‐number technique

A cell‐based long‐term hydrological model (CELTHYM) that can be integrated with a geographical information system (GIS) was developed to predict continuous stream flow from small agricultural watersheds. The CELTHYM uses a cell‐by‐cell soil moisture balance approach. For surface runoff estimation, the curve number technique considering soil moisture on a daily basis was used, and release rate was used to estimate baseflow. Evapotranspiration was computed using the FAO modified Penman equation that considered land‐use‐based crop coefficients, soil moisture and the influence of topography on radiation. A rice paddy field water budget model was also adapted for the specific application of the model to East Asia. Model sensitivity analysis was conducted to obtain operational information about the model calibration parameters. The CELTHYM was calibrated and verified with measured runoff data from the WS#1 and WS#3 watersheds of the Seoul National University, Department of Agricultural Engineering, in Hwaseong County, Kyounggi Province, South Korea. The WS#1 watershed is comprised of about 35·4% rice paddy fields and 42·3% forest, whereas the WS#3 watershed is about 85·0% forest and 11·5% rice paddy fields. The CELTHYM was calibrated for the parameter release rate, K, and soil moisture storage coefficient, STC, and results were compared with the measured runoff data for 1986. The validation results for WS#1 considering all daily stream flow were poor with R², E² and RMSE having values of 0·40, ?6·63 and 9·69 (mm), respectively, but validation results for days without rainfall were statistically significant (R² = 0·66). Results for WS#3 showed good agreement with observed data for all days, and R², E² and RMSE were 0·92, 0·91 and 2·23 (mm), respectively, suggesting potential for CELTHYM application to other watersheds. The direct runoff and water balance components for watershed WS#1 with significant areas of paddy fields did not perform well, suggesting that additional study of these components is needed. Copyright © 2002 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.     

Potentials and limitations of modelling vertical distributions of root water uptake of an Austrian pine forest on a sandy soil

Root water uptake patterns are often studied with simulation models of the unsaturated soil water flow, as they are difficult to measure directly. Calibration of these models is not straightforward and causes uncertainties in simulated uptake distributions. In this paper we study how uncertainties in the calibration of the SWIF model affect uncertainty intervals in simulated uptake patterns of an Austrian pine stand (Pinus nigra var. nigra) on a sandy soil. After calibrating and validating SWIF with a large data set of more than 125 000 measured soil water contents over a three year period, uncertainty ranges in simulated soil water dynamics and root water uptake distributions were estimated with a Monte Carlo analysis. In general, uncertainties in root uptake patterns were small (typically <2 10⁻⁴ m³ m⁻³ day⁻¹) and were higher for trees with a shallow rooting system (0·8 m) than for trees with a deep rooting system (2·5 m). Uncertainties arose mainly from uncertainties in simulated soil water fluxes and from variations in the reduction of uptake during periods of drought. Uncertainties in soil water contents were far higher (typically 0·01 m³ m⁻³) than uncertainties in uptake, illustrating that uncertainties in uptake parameters and those in the distribution of water uptake hardly affect the modelling of soil water dynamics. Root water uptake models should therefore be validated against measured uptake distributions, which can be determined on sandy soils during dry periods with a high water use when soil fluxes are negligible to uptake. Copyright © 2000 John Wiley & Sons, Ltd.     

Integration of remotely sensed C factor into SWAT for modelling sediment yield

The C factor, representing the impact of plant and ground cover on soil loss, is one of the important factors of the Modified Universal Soil Loss Equation (MUSLE) in the Soil and Water Assessment Tool (SWAT) to model sediment yield. The daily update of C factors in SWAT was originally determined by land use types and plant growth cycles. This does not reflect the spatial variation of C values that exists within a large land use area. We present a new approach to integrate remotely sensed C factors into SWAT for highlighting the effect of detailed vegetative cover data on soil erosion and sediment yield. First, the C factor was estimated using the abundance of ground components extracted from remote sensing images. Then, the gridding data of the C factor were aggregated to hydrological response units (HRUs), instead of to land use units of SWAT. In the end, the C factor values in HRUs were integrated into SWAT to predict sediment yield by modifying the ysed subroutine. This substitution work not only increases the spatial variation of the C factor in SWAT, but also makes it possible to utilize other sources of C databases rather than those from the United States. The demonstration in the Dage basin shows that the modified SWAT produces reasonable results in water flow simulation and sediment yield prediction using remotely sensed C values. The Nash–Sutcliffe efficiency coefficient (E_NS) and R² for surface runoff range from 0·69 to 0·77 and 0·73 to 0·87, respectively. The coefficients E_NS and R² for sediment yield were generally above 0·70 and 0·60, respectively. The soil erosion risk map based on sediment yield prediction at the HRU level illustrates instructive details on spatial distribution of soil loss. Copyright © 2011 John Wiley & Sons, Ltd.     

Effect of watershed subdivision on simulation of water balance components

The present effect of watershed subdivision on simulated water balance components using the thoroughly tested Soil and Water Assessment Tool (SWAT) model has been evaluated for the Nagwan watershed in eastern India. Observed meteorological and hydrological data (daily rainfall, temperature, relative humidity and runoff) for the years 1995 to 1998 were collected and used. The watershed and sub‐watershed boundaries, slope and soil texture maps were generated using a geographical information system. A supervised classification method was used for land‐use/cover classification from satellite imagery of 1996. In order to study the effect of watershed subdivision, the watershed was spatially defined into three decomposition schemes, namely a single watershed, and 12 and 22 sub‐watersheds. The simulation using the SWAT model was done for a period of 4 years (1995 to 1998). Results of the study showed a perfect water balance for the Nagwan watershed under all of the decomposition schemes. Results also revealed that the number and size of sub‐watersheds do not appreciably affect surface runoff. Except for runoff, there was a marked variation in the individual components of the water balance under the three decomposition schemes. Though the runoff component of the water balance showed negligible variation among the three cases, variations were noticed in the other components: evapotranspiration (5 to 48%), percolation (2 to 26%) and soil water content (0·30 to 22%). Thus, based on this study, it is concluded that watershed subdivision has a significant effect on the water balance components. Copyright © 2005 John Wiley & Sons, Ltd.     

A bimodal four‐parameter lognormal linear model of soil water repellency persistence

Soil water repellency may be characterized in terms of the delayed infiltration time of a water droplet resting on the soil surface, which is, water drop penetration time (WDPT), or repellency persistence. Such repellency persistence varies nonlinearly with soil water content (θ_g), although no models have been proposed to reproduce the variation of WDPT with θ_g in soils. Dynamic factor analysis (DFA) is used to identify two common patterns of unexplained variability in a scattered dataset of WDPT versus θ_g measurements. A four‐parameter lognormal distribution was fitted to both common patterns obtained by DFA, and these were combined additively in a weighted multiple linear bimodal model. We show how such an empirical model is capable of reproducing a large variety of WDPT versus θ_g curve shapes (N = 80) both within a wide range of measured WDPTs (0–17 000 s) and for samples with organic matter content ranging from 21·7 to 80·6 g (100 g)^?1. Copyright © 2009 John Wiley & Sons, Ltd.     

Near‐saturated surface soil hydraulic properties under different land uses in the St Denis National Wildlife Area,Saskatchewan, Canada

Surface soil hydraulic properties are key factors controlling the partition of rainfall and snowmelt into runoff and soil water storage, and their knowledge is needed for sound land management. The objective of this study was to evaluate the effects of three land uses (native grass, brome grass and cultivated) on surface soil hydraulic properties under near‐saturated conditions at the St Denis National Wildlife Area, Saskatchewan, Canada. For each land use, water infiltration rates were measured using double‐ring and tension infiltrometers at ?0·3, ?0·7, ?1·5 and ?2·2 kPa pressure heads. Macroporosity and unsaturated hydraulic properties of the surface soil were estimated. Mean field‐saturated hydraulic conductivity (K_fs), unsaturated hydraulic conductivity at ?0·3 kPa pressure head, inverse capillary length scale (α) and water‐conducting macroporosity were compared for different land uses. These parameters of the native grass and brome grass sites were significantly (p < 0·1) higher than that of the cultivated sites. At the ?0·3 kPa pressure head, hydraulic conductivity of grasslands was two to three times greater than that of cultivated lands. Values of α were about two times and values of K_fs about four times greater in grasslands than in cultivated fields. Water‐conducting macroporosity of grasslands and cultivated fields were 0·04% and 0·01% of the total soil volume, respectively. Over 90% of the total water flux at ?0·06 kPa pressure head was transmitted through pores > 1·36 × 10^?4 m in diameter in the three land uses. Land use modified near‐saturated hydraulic properties of surface soil and consequently may alter the water balance of the area by changing the amount of surface runoff and soil water storage. Copyright © 2004 John Wiley & Sons, Ltd.     

Estimates of spatial variation in evaporation using satellite‐derived surface temperature and a water balance model

Evaporation dominates the water balance in arid and semi‐arid areas. The estimation of evaporation by land‐cover type is important for proper management of scarce water resources. Here, we present a method to assess spatial and temporal patterns of actual evaporation by relating water balance evaporation estimates to satellite‐derived radiometric surface temperature. The method is applied to a heterogeneous landscape in the Krishna River basin in south India using 10‐day composites of NOAA advanced very high‐resolution radiometer satellite imagery. The surface temperature predicts the difference between reference evaporation and modelled actual evaporation well in the four catchments (r² = 0·85 to r² = 0·88). Spatial and temporal variations in evaporation are linked to vegetation type and irrigation. During the monsoon season (June–September), evaporation occurs quite uniformly over the case‐study area (1·7–2·1 mm day^?1), since precipitation is in excess of soil moisture holding capacity, but it is higher in irrigated areas (2·2–2·7 mm day^?1). In the post‐monsoon season (December–March) evaporation is highest in irrigated areas (2·4 mm day^?1). A seemingly reasonable estimate of temporal and spatial patterns of evaporation can be made without the use of more complex and data‐intensive methods; the method also constrains satellite estimates of evaporation by the annual water balance, thereby assuring accuracy at the seasonal and annual time‐scales. Copyright © 2007 John Wiley & Sons, Ltd.