Variable Source Area Research Articles

Catchment Coevolution and the Geomorphic Origins of Variable Source Area Hydrology

AbstractFeatures of landscape morphology—including slope, curvature, and drainage dissection—are important controls on runoff generation in upland landscapes. Over long timescales, runoff plays an essential role in shaping these same features through surface erosion. This feedback between erosion and runoff generation suggests that modeling long‐term landscape evolution together with dynamic runoff generation could provide insight into hydrological function. Here we examine the emergence of variable source area runoff generation in a new coupled hydro‐geomorphic model that accounts for water balance partitioning between surface flow, subsurface flow, and evapotranspiration as landscapes evolve over millions of years. We derive a minimal set of dimensionless numbers that provide insight into how hydrologic and geomorphic parameters together affect landscapes. Across the parameter space we investigated, model results collapsed to a single inverse relationship between the dimensionless relief and the ratio of catchment quickflow to discharge. Furthermore, we found an inverse relationship between the Hillslope number, which describes topographic relief relative to aquifer thickness, and the proportion of the landscape that was variably saturated. While the model generally produces fluvial topography visually similar to simpler landscape evolution models, certain parameter combinations produce wide valley bottom wetlands and non‐dendritic, trellis‐like drainage networks, which may reflect real conditions in some landscapes where aquifer gradients become decoupled from topography. With these results, we demonstrate the power of hydro‐geomorphic models for generating new insights into hydrological processes, and also suggest that subsurface hydrology may be integral for modeling aspects of long‐term landscape evolution.

Watershed scale modeling of Dissolved organic carbon export from variable source areas

Dissolved organic carbon (DOC) in surface water influences the global carbon cycle, ecosystem productivity, and water quality in potable water supply systems. Few physically based watershed models can simulate carbon cycling and predict DOC in surface waters under the influence of natural and anthropogenic drivers. In this work we transform the SWAT-Carbon (C) model to simulate DOC from variable source runoff areas in a humid forested watershed in the northeastern United States. Remotely sensed data were used to parameterize and simulate forest growth and evapotranspiration. The calibrated model accurately simulated streamflow and DOC flux at six sites across the watershed compared to measured data. The DOC predictions across sites showed model sensitivity to soil properties, particularly soil depth and available water capacity. The spatial distribution of DOC export across the watershed followed the pattern of surface runoff from variable source areas. Model sensitivity of DOC flux showed greater sensitivity to changes in precipitation than temperature. The overall good performance of the model makes it a valuable tool for watershed-scale modeling of DOC to understand the influence of climate and watershed management towards developing mitigation strategies. The methods presented in this study can be used in forested watersheds in regions where runoff from variable source areas is important for water quality predictions.

Riparian zone as a variable source area for the estimation of evapotranspiration through the analysis of daily fluctuations in streamflow

AbstractEvapotranspiration is a critical component of the water balance of a catchment, and riparian zones play a crucial role in the hydrological process. However, the influence of riparian zones on evapotranspiration, especially at a catchment scale, needs more investigation. Since evapotranspiration is negligible at night, the difference in recharge between day and night generates diurnal fluctuations in streamflows. A new method is developed to estimate riparian evapotranspiration by analysing the streamflow response to diurnal evapotranspiration. The potential evapotranspiration (PET) calculated by FAO‐56 is used as a true reference for validation. The estimated evapotranspiration shows acceptable goodness of fit with the hourly PET calculated from FAO 56. However, the new approach estimates average evapotranspiration across the entire riparian area than a point observation and provides a better understanding of the impact of the riparian zone on streamflow patterns.

Targeting the critical source areas of phosphorus based on hydrological sensitive area delineation to control nonpoint source pollution in watersheds

To pursue higher agricultural output, farmers have invested heavily in fertilizers, pesticides and other factors of production, which inevitably causes agricultural nonpoint source pollution (ANPSP) and seriously restricted the sustainable development of global agriculture. Accurate and efficient integrated prevention and control of ANPSP requires scientific identification of critical source areas and influencing factors of nonpoint source pollution, which is particularly important in geographic conditions where steep-slope agriculture dominates and landscape features are complex. In this study, a typical small watershed in the mountainous area of the Three Gorges Reservoir area was selected, and the risk of nonpoint source phosphorus loss in the watershed was assessed using the Phosphorus Index (PI) on the basis of Soil Topographic Index (STI) prediction of runoff production area in the watershed. Finally, the critical source areas (CSAs) for its loss were identified. The results showed that the area of hydrologically sensitive areas in the watershed was approximately 22.87% of the total study area, with an area of 37.05 hm2. The hydrology of variable source areas within hydrologically sensitive areas plays a dominant role and is more likely to produce runoff and increase the risk of agricultural nonpoint source phosphorus pollution. The results of PI analysis showed that the risk of phosphorus loss was low in most areas of the watershed (45.77%), the areas at high risk of phosphorus loss accounted for 29.33% of the entire watershed with an area of 47.52 hm2, and tea plantation was the most dominant land use type. The critical source area of nonpoint source phosphorus output accounts for 5.47% of the watershed area and is mainly distributed along the watershed system, 82.55% of which is tea plantation. Compared with the traditional phosphorus index evaluation system, this method is important for the accurate identification of critical source areas of nonpoint source pollutants, for strengthening the efficient prevention, control and management of ANPSP, and for improving the utilization rate of arable land as well as land management and sustainable development of watersheds.

Development and evaluation of temperature-induced variable source area runoff generation model

In permafrost regions, temperature and precipitation play decisive roles in hydrological processes. Soil freeze–thaw cycles in the active layer cause the runoff generation process to show multiple mechanisms and seasonal alternating patterns. In this study, based on the precipitation-runoff modeling system implemented in the Java modeling framework Object Modeling System (PRMS-OMS), we developed a runoff generation model with temperature-induced variable source area (TVSA) for permafrost regions by introducing an active layer parameterization scheme, an active layer freeze–thaw module, a glacier module, and the sub-permafrost groundwater module. The TVSA model was calibrated and validated using field observations in Fenghuoshan (FHS), a typical permafrost watershed, and Tuotuohe (TTH), a typical permafrost-glacier catchment, on the Qinghai-Tibet Plateau (QTP). During the calibration period, parameters related to the thermal conditions of the active layer and runoff processes were calibrated using the observed freeze and thaw depth data and the discharge data. For the validation period, the model successfully reproduced the freeze–thaw processes (average root mean squared error = 0.205 m) and discharges (FHS, Nash–Sutcliffe efficiency coefficient (NSE) = 0.99; TTH, NSE = 0.85). The TVSA model is a powerful tool for the system identification of variable source runoff generation processes and the associated physical mechanism under temperature control in permafrost basins; additionally, the model accurately reproduces hydrological processes and the associated response to climate change.

A Modified Distributed CN-VSA Method for Mapping of the Seasonally Variable Source Areas

Many watershed models employ the Soil Conservation Service Curve Number (SCS-CN) approach for runoff simulation based on soil and land use information. These models implicitly assume that runoff is generated by the Hortonian process and; therefore, cannot correctly account for the effects of topography, variable source area (VSA) and/or soil moisture distribution in a watershed. This paper presents a new distributed CN-VSA method that is based on the SCS-CN approach to estimate runoff amount and uses the topographic wetness index (TWI) to distribute the runoff-generating areas within the watershed spatially. The size of the saturated-watershed areas and their spatial locations are simulated by assuming an average annual value of potential maximum retention. However, the literature indicates significant seasonal variation in potential maximum retention which can considerably effect water balance and amount of nonpoint source pollution. This paper focuses on developing a modified distributed CN-VSA method that accounts for the seasonal changes in the potential maximum retention. The results indicate that the modified distributed CN-VSA approach is better than distributed CN-VSA to simulate runoff amount and spatial distribution of runoff-generating areas. Overall, the study results are significant for improved understanding of hydrological response of watershed where seasonal factors describe the potential maximum retention, and, thus, saturation excess runoff generation in the watershed.

An insight into the hydrological aspects of landslides of 2018 in Kodagu, South India

The goal of the work has been to assess the factors that led to the series of landslides, which ravaged the District of Kodagu situated in the mountainous Western Ghats, in 2018. The work has been accomplished through field investigations and analysis of field data. While commonly, studies in India lay emphasis on the geological aspects and anthropogenic factors associated with landslides, the present work intended to study in detail the hydrological aspects of the disaster, which wreaked havoc in the region characterised by deep soils, very heavy rainfall and a dense vegetation cover. All the slides of significant magnitude along a 25-km stretch of the Payaswini valley, draining the western slopes of the Madikeri hill range, were identified, and their land characteristics delineated using field surveys and maps. Soil characteristics of the region were determined through tests on samples collected from ten slide sites. Results from tests on texture, density and permeability were used to learn about the part played by subsurface water. Characteristics of rainfall during the year were investigated using field data from the worst affected region. Studies on progress of the Monsoon rain, an analysis of probability of extreme values and establishment of intensity–duration relationships were accomplished. A variable source area watershed model incorporating MODFLOW, the popular ground water flow simulation model, has been applied on a typical micro-catchment in order to study the growth of the saturated zone during rains. It has been inferred from the investigations that an unprecedented combination of extremely wet catchment and a very heavy spell of rainfall, rather than any commonly attributed factor, has been the primary cause of the disaster in the district. It is concluded that it would be a futile exercise working on this natural phenomenon devoid of extensive field studies, because of the uniqueness of the region.

Integrating Daily CO2 Concentrations in SWAT-VSA to Examine Climate Change Impacts on Hydrology in a Karst Watershed

HighlightsWe used SWAT-VSA to assess the effects of climate change with rising CO2 on the water balance of a karst basin.For future climate, SWAT-VSA with rising CO2 yielded 7.1% less ET and 6.3% more runoff than standard SWAT-VSA.Rising CO2 also affected variable source areas, with greater ET declines and runoff increases in the wettest soils.Findings suggest CO2 effects on water balance should be included in future climate change studies with SWAT-VSA.Abstract. Characterizing the effects of climate change on hydrology is important to watershed management. In this study, we used SWAT-VSA to examine the effects of climate change and increasing atmospheric CO2 (CO2) on the water balance of Spring Creek watershed, a mixed land-use karst basin in the Upper Chesapeake Bay watershed. First, we modified the stomatal conductance and leaf area index (LAI) routines of SWAT-VSA’s Penman-Monteith evapotranspiration (ET) procedure and enabled the model to accept daily CO2 data. Using downscaled climate projections from nine global climate models (GCMs), we then compared water balance estimations from baseline SWAT-VSA against two modified versions of SWAT-VSA. One SWAT-VSA version integrated daily CO2 levels (SWAT-VSA_CO2), while another version added flexible stomatal conductance and LAI routines (SWAT-VSA_CO2+Plant) to the dynamic CO2 capacity. Under current climate (1985-2015), the three SWAT-VSA models produced generally similar water balance estimations, with 51% of precipitation lost to ET and the remainder converted to runoff (10%), lateral flow (9%), and percolate (30%). For future climate (2020-2065), water balance simulations diverged between baseline SWAT-VSA and the two modified SWAT-VSA models with CO2. Notably, variable stomatal conductance and LAI routines produced no detectable effects beyond that of CO2. For the 2020-2065 period, baseline SWAT-VSA projected ET increases of 0.7 mm year-1, while SWAT-VSA models with CO2 suggested that annual ET could decline by approximately -0.4 mm year-1 over the same period. As a result, the two CO2-based SWAT-VSA models predicted streamflow increases of almost 1.6 mm year-1 over the 2020-2065 period, which were roughly double the streamflow increases projected by baseline SWAT-VSA. In general, SWAT-VSA models with CO2 effects produced 22.4% more streamflow in 2045-2065 than the SWAT-VSA model without CO2. Results also showed that adding daily CO2 to SWAT-VSA reduced ET in wetter parts of Spring Creek watershed, leading to greater runoff losses from variable source areas compared to baseline SWAT-VSA. Findings from the study highlight the importance of considering increasing atmospheric CO2 concentrations in water balance simulations with SWAT-VSA in order to gain a fuller appreciation of the hydrologic uncertainties with climate change. Keywords: Carbon dioxide, Climate change, Hydrologic model, Water balance, Watershed.

A comparative evaluation of the continuous and event-based modelling approaches for identifying critical source areas for sediment and phosphorus losses

Proper identification of critical source areas (CSAs) is important for economic viability of any best management practices (BMPs) aimed at reducing sediment and phosphorus loads to receiving water bodies. Both continuous and event-based hydrologic and water quality models are widely used to identify and assess CSAs, however, their comparative assessment is lacking. In this study, we have used continuous Soil and Water Assessment Tool (SWAT) and event-based Agriculture Non-Point Source (AGNPS) pollution models to identify CSAs for sediment and phosphorus in a watershed in Ontario, Canada. Along with their original version, both models were re-conceptualized to incorporate saturation excess mechanism of runoff generation, which is also refereed as variable source area (VSA) integration. The models were set-up using high resolution spatial, crop- and land-management, and meteorological dataset; and calibrated with reasonable accuracy against streamflow, sediment and phosphorus concentration data at multiple locations. Threshold value (t-value) approach was used to identify CSA areas in the watershed. Results showed that both models were in agreement (up to 96% of fields) that summer season did not constitute hot-moments (<6% of the watershed area as CSAs) for both sediment and phosphorus. SWAT models identified winter (~50% of watershed area as CSA) and AGNPS models identified early spring (~50% of watershed areas as CSAs) season as the hot-moment for both sediment and phosphorus. Contrasting result, as indicated by low (1%) matching in field CSA potential, was observed in autumn season. In the same season, VSA integrated SWAT and AGNPS models showed better matching (43% for sediment and 31% for phosphorus), highlighting the importance of VSA integration in the models. Qualitative validation of model-based CSA potential with oblique aerial-photograph-based CSA potential in two soil moisture conditions (wetter and drier) indicated slightly better performance of the SWAT models, and over-prediction of the AGNPS models. However, a more comprehensive analysis based on more detailed field observations is needed to further confirm the results.

On the contribution of groundwater to streamflow in laterite catchments of the Darling Range, south‐western Australia

AbstractIn deeply weathered laterite catchments of the Darling Range in south‐western Australia, the direct contribution (i.e., discharge) of permanent groundwater to streamflow has long been considered as minor. Instead, downslope shallow throughflow was thought to dominate, generating more than 90% of streamflow. We used a chemical hydrograph separation approach to estimate annual groundwater discharge for three catchments over periods of up to 39 years, and found that direct groundwater contributions to streamflow were far more variable across catchments and through time than has previously been acknowledged. The estimated proportion of annual streamflow sourced directly from groundwater ranged from 0 to 93% and was related linearly to the size of the groundwater discharge area in the catchment valley floor. In contrast, contributions from shallow sources including shallow throughflow varied primarily and linearly with annual rainfall. However, the response to rainfall was “amplified” in a predictable way by the size of the groundwater discharge area, consistent with the variable source area concept. We derived a functional relationship between catchment annual rainfall‐runoff ratio and groundwater discharge area and successfully applied this to a further four catchments, inferring that the results were broadly applicable across the Darling Range. The implications for an improved understanding of streamflow generating processes in the study region, and for laterite catchments generally, are discussed.

Modeling of phosphorus loss from field to watershed: A review.

Phosphorus (P) losses from nonpoint sources into surface water resources through surface runoff and tile drainage play a significant role in eutrophication. Accordingly, the number of studies involving the modeling of agricultural P losses, the uncertainties of such models, and the best management practices (BMPs) supported by the modeling of hypothetical P loss reduction scenarios has increased significantly around the world. Many improvements have been made to these models: separate manure P pools, variable source areas allowing the determination of critical source areas of P loss, analyses of modeling uncertainties, and understanding of legacy P. However, several elements are still missing or have yet to be sufficiently addressed: the incorporation of preferential flow into models, the modification of P sorption-desorption processes considering recent research data (e.g., pedotransfer functions for labile, active, or stable P, along with P sorption coefficients), BMP parameterization, and scale-up issues, as well as stakeholder-scientist and experimentalist-modeler interactions. The accuracy of P loss modeling can be improved by (a) incorporating dynamic P sorption-desorption processes and new P subroutines for direct P loss from manure, fertilizer, and dung, (b) modeling preferential flow, connectivity between field and adjacent water bodies, and P in-stream processes, (c) including an assessment of model uncertainty, (d) integrating field and watershed models for BMP calibration and scaling field results up to larger areas, and (e) building a holistic interaction between stakeholders, experimentalists, and modelers.

Pairing soil sampling with very-high resolution UAV imagery: An examination of drivers of soil and nutrient movement and agricultural productivity in southern Ontario

Soil erosion from agricultural lands continues to be a global societal problem. The movement of soils is often accompanied by nitrogen and phosphorus that are crucial to crop growth, but their redistribution from farm fields to waterways can reduce crop yields and degrade water quality. While within-field sediment and nutrient movement has been quantified using small plots and edge-of-field monitoring, these approaches fail to capture their spatial distribution. The pairing of soil sampling with unmanned aerial vehicle (UAV) data offers a novel and low-cost approach to map the spatial distribution of soil characteristics and nutrient concentrations within a farm field. UAV data are used to generate a digital terrain model and subsequently map within-field topographic variation and erosional flow pathways. Topographic variation is discretized into landform elements (flat, shoulder, backslope, footslope) that capture within-field heterogeneity and have potential for scaling out soil sampling to larger spatial extents. Our results show the controlling factor of water content and organic matter on crop yield, as represented by normalized difference vegetation index (NDVI). Significant differences in water content and organic matter were found across landform elements with increases in both parameters downslope. Upslope landform elements contained more sand content (9–20%) and had lower NDVI values than downslope elements. Complementing these findings, significant differences in organic matter, soluble nitrogen, and soluble reactive phosphorus occurred along erosional flow pathways. Our within-field results have implications for farmers, as our analysis of soil characteristics indicated that NDVI was positively correlated with water content (0.05), organic matter (0.15), silt (0.36), and clay (0.17) content and negatively correlated with soluble nitrogen (−0.47) and phosphorus (−0.30) concentrations. In addition to discussing the challenges and opportunities for expanding upon the presented research, we use a simple proof-of-concept hydrological model to demonstrate the potential role of hydrological connectivity and variable source area as a driver of within-field nutrient movement. The combination of our empirical results showing water content and organic matter as controlling factors on agricultural yield, the role of hydrological connectivity, and climate predictions of increased storm intensity suggest that additional research into the generation of novel time-series soil sampling and UAV-generated erosion and flow data could advance our understanding of the variations in soil characteristics and nutrient concentrations within individual farm fields.

Mapping runoff generating areas using AGNPS-VSA model

ABSTRACT In humid regions, surface runoff is often generated by saturation-excess runoff mechanisms from relatively small variable source areas (VSAs). However, the majority of the current hydrologic models are based on infiltration-excess mechanisms. In this study, the AGricultural Non-Point Source Pollution (AGNPS) model was used to integrate the VSA concept using topographic wetness index (TWI). Both the original and AGNPS-VSA models were evaluated for a small agricultural field in Ontario, Canada. The results indicate that the AGNPS-VSA model performed better than original model. The AGNPS-VSA model predicted that only the saturated portion of the field with higher TWI values produced runoff, whereas the original AGNPS model showed uniform hydrologic response from the entire field. The results of this study are important for accurately mapping the locations of VSAs. This new model could be a powerful tool in identifying critical source areas for applying targeted best management practices to minimize pollutant loads to receiving waters.

Sensitivity of nitrate concentration‐discharge patterns to soil nitrate distribution and drainage properties in the vertical dimension

We propose a conceptual model that examines the ‘variable source area’ (VSA) and ‘nitrate flushing’ hypothesis in the vertical direction, and use this approach to explain nitrate concentration–discharge relationships in a semi-arid watershed. We use an eco-hydrology simulation model (RHESSys) to show that small changes in the vertical distribution of nitrate mass and their interaction with soil hydraulic conductivity can result in abrupt changes in the nitrate concentration–discharge relationship. We show that the estimated concentration–discharge relationship is sensitive to the parameters governing soil vertical nitrate distribution and soil hydraulic conductivity, at both patch scale and watershed scale, where lateral redistribution of water and nitrate is also accounted for. Given heterogeneity in nitrogen inputs, uptake processes, soil drainage and storage processes, substantial variation in parameters that describe rate of changes in vertical distribution of soil nitrate and hydraulic properties is likely both within and between watersheds. Thus, we argue that vertical ‘variable source area’ processes may be as important as lateral VSA in determining concentration discharge relationships.

Forest management influences the effects of streamside wet areas on stream ecosystems

Riparian zones contain areas of strong hydrological connectivity between land and stream, referred to as variable source areas (VSAs), and are considered biogeochemical control points. However, little is known about whether VSAs influence stream communities and whether this connectivity is affected by forest management. To address this, we used multiple biotic and abiotic indicators to (1) examine the influence of VSAs on riparian vegetation and stream ecosystems by comparing VSA and non-VSA reaches and (2) explore how forest management may affect the influence of VSAs on stream ecosystems. We detected some significant differences between VSA and non-VSA reaches in the riparian vegetation (greater understory and lower tree density) and stream ecosystem indicators (greater dissolved organic matter aromaticity, microbial biomass, peroxidase activity and collector-gatherer density, and lower dissolved organic carbon concentrations, algal biomass, and predatory macroinvertebrate density), which suggests that VSAs may create a more heterotrophic ecosystem locally. However, we show some evidence that forest management activities (specifically, road density) can alter the influence of VSAs and eliminate the differences observed at lower forest management intensities, and that the most hydrologically connected areas seem more sensitive to disturbance. Therefore, we suggest that the heterogeneity in hydrological connectivity along riparian zones should be considered when planning forest harvesting operations and road building (e.g., wider riparian buffers around VSAs).

Simulation and statistical modelling approaches to investigate hydrologic regime transformations following Eastern hemlock decline

AbstractEcohydrological processes occurring at or near the Earth's surface are strongly influenced by Eastern hemlock [EH; Tsuga canadensis (L.) Carrière], a foundation tree species of eastern North American forests. EH populations are currently threatened by the invasive hemlock woolly adelgid (HWA; Adelges tsugae Annand). HWA HWA populations have been expanding rapidly throughout the EH's range. Catchment‐scale research examining the hydrological consequences of HWA infestation is lacking, and plot‐scale studies remain conflicted in their findings. Given the complex relationships between canopy interception, unsaturated and saturated groundwater storage, and root water uptake, it is not immediately clear how EH loss will affect the hydrologic cycle. We investigated the impact of EH mortality on stream discharge characteristics across a regional sample of catchments utilizing both simulation and statistical modelling approaches. We first examined the relationship between various catchment characteristics, including EH health, and three hydrological variables through regression analysis. We then employed a non‐parametric statistical test to evaluate differences in hydrologic regime trends between non‐infested and infested catchments. Finally, we calibrated a physically based hydrologic model and considered differences in optimal model parameter values and simulated overland runoff between non‐infested and infested catchments. HWA presence modified several ecohydrological characteristics and precipitation partitioning between groundwater flows and surface runoff, potentially driving higher stream flashiness and overland flow, lower baseflow contributions and catchment storage, shorter flow‐path lengths, and variable source area dilation at infested sites. Our results suggest that EH decline will augment flooding potential associated with the increasing frequency and intensity of Atlantic Basin tropical cyclone events. Further, our physically based simulation provides more determinate results than regression analysis, indicating that a purely statistical methodology, commonly utilized in studying the relationship between landcover characteristics and hydrologic regime, neglects dynamic physical ecohydrologic relationships.

From Points to Patterns: Using Groundwater Time Series Clustering to Investigate Subsurface Hydrological Connectivity and Runoff Source Area Dynamics

AbstractGroundwater levels are typically measured at only a limited number of points in a catchment. Thus, upscaling these point measurements to the catchment scale is necessary to determine subsurface flow paths and runoff source areas. Here we present a data‐driven approach composed of time series clustering and topography‐based upscaling of shallow, perched groundwater dynamics using groundwater data from 51 monitoring sites in a 20‐ha prealpine headwater catchment in Switzerland. The agreement between the upscaled (modeled) and measured groundwater dynamics was strong for most of the 19‐month study period for the upslope and footslope locations but weaker at the beginning of events and for the midslope locations. However, these differences between measured and modeled groundwater levels did not significantly affect modeled groundwater activation, that is, the time when groundwater levels were within the more transmissive soil layers near the soil surface. The resulting groundwater activation maps represent the groundwater response across the catchment and highlight the dynamic expansion and contraction of the subsurface runoff source areas, particularly along the channel network. This is in agreement with the variable source area concept. However, there were also isolated active zones that did not get connected to the stream during rainfall events, highlighting the need to distinguish between variable active and variable stream‐connected runoff source areas. Our data‐driven approach to upscale point measurements of shallow groundwater levels appears useful for studying catchment‐scale variations in groundwater storage and connectivity and thus may help to better understand runoff generation in mountain catchments.

Assessing runoff generation in riparian wetlands: monitoring groundwater-surface water dynamics at the micro-catchment scale.

Riparian wetlands (RW) are important variable source areas for runoff generation. They are usually characterised by a combination of groundwater exfiltration-which maintains saturated conditions in low-lying organic-rich soils-and direct precipitation. Both processes interact to generate overland flow as a dominant runoff process. The small-scale details of groundwater-surface water (GW-SW) interactions are usually not well understood in RW. Here, we report the results of a study from an experimental catchment in the Scottish Highlands where spatio-temporal runoff processes in RW were investigated using isotopes, alkalinity and hydrometric measurements. We focused on perennial micro-catchments within the RW and ephemeral zero-order channels draining peatland hollows and hummocks to better understand the heterogeneity in GW-SW interactions. The 12-month study period was dominated by the wettest winter (December/January) period on record. Runoff generation in the RW was strongly controlled by the local groundwater response to direct rainfall, but also the exfiltration of groundwater from upslope. This groundwater drainage is focused in the hollows in ephemeral and perennial drainage channels, but in wet conditions, as exfiltration rates increase, can affect hummocks as well. The hollows provide the dominant areas for mixing groundwater, soil water and direct rainfall to deliver water to the stream network as hollows "fill and spill" to increase connectivity. They also provide wet areas for evaporation which is evident in enriched isotope signatures in summer. Although there is some degree of heterogeneity in the extent to which groundwater influences specific micro-catchments, particularly under low flows, the overall isotopic response is quite similar, especially when the catchment is wet and this responses can explain the isotope signatures observed in the stream. In the future, more longitudinal studies of micro-catchments are needed to better explain the heterogeneity observed.

The worthiness of using information on land-use–land-cover in watershed models for Western Ghats: A case study

The variable source area (VSA) theory of runoff generation mechanisms has been proved to hold good in many wet mountainous areas, decades ago. According to this theory, infiltration-excess overland flow is limited to very small areas in mountainous and forested catchments. But, the perception that the land surface characteristics, including land-use–land-cover (LULC), form the major factors influencing the response of the catchment to rainfall has dominated the thought in hydrology to such an extent that models based on the overland flow theory continue to be used even in such areas. The present study was taken up in order to understand the worthiness of using parameters, including the curve number (CN), that are based on the physiographic characteristics of the catchment in a watershed model designed to estimate runoff in the wet mountainous areas of the Western Ghats in southern India, where the VSA theory has been proved to hold good. The study has been accomplished by applying the NITK model developed for estimating runoff using daily rainfall data. This model is believed to estimate reliably the streamflow in the region using parameter values that can be computed from catchment characteristics. In the present study, it is applied on three gauged streams in the region of Western Ghats in Karnataka. Initially, the performance of the model has been studied with the parameters fixed using the catchment characteristics. Later, the model has been used as a tool to test hypotheses concerning the catchment response, by varying the parameter values, adopting a trial and error procedure. Initial results showed that the model performance is poor as the coefficients of efficiency vary between –66.9 and 82%. The sensitivity analysis carried out subsequently showed that the model parameters are required to be altered greatly for good performance and that the model simulations are not sensitive to the parameter CN. Further, the performance of this model was compared with that of a VSA model, known to suit the region well. This showed that even after all the changes in the model parameters, the model results are not highly reliable. Hence, in order to understand the reasons for the poor performance of the model, a technique was developed to compute the CN values that would be actually necessary to simulate daily direct runoff (DRO) reliably in this method, the daily values of CN are computed by applying backwards the expression for runoff on the DRO estimated by the VSA model. The variations in the values of CN computed using this method are then studied. It is found that the variations in daily CN are high and highly random too, whereas the NITK model uses only three fixed values of CN. It is thus concluded that factors other than those on which the CN is popularly believed to depend control the runoff generation in the region and that influence of LULC on runoff is not discernible at all from the kind of data that is commonly available.

The Use of River Flow Discharge and Sediment Load for Multi-Objective Calibration of SWAT Based on the Bayesian Inference

The soil and water assessment tool (SWAT) is widely used to quantify the spatial and temporal patterns of sediment loads for watershed-scale management of sediment and nonpoint-source pollutants. However few studies considered the trade-off between flow and sediment objectives during model calibration processes. This study proposes a new multi-objective calibration method that incorporates both flow and sediment observed information into a likelihood function based on the Bayesian inference. For comparison, two likelihood functions, i.e., the Nash–Sutcliffe efficiency coefficient (NSE) approach that assumes model residuals follow the Gaussian distribution, and the BC-GED approach that assumes model residuals after Box–Cox transformation (BC) follow the generalized error distribution (GED), are applied for calibrating the flow and sediment parameters of SWAT with the water balance model and the variable source area concept (SWAT-WB-VSA) in the Baocun watershed, Eastern China. Compared with the single-objective method, the multi-objective approach improves the performance of sediment simulations without significantly impairing the performance of flow simulations, and reduces the uncertainty of flow parameters, especially flow concentration parameters. With the NSE approach, SWAT-WB-VSA captures extreme flood events well, but fails to mimic low values of river discharge and sediment load, possibly because the NSE approach is an informal likelihood function, and puts greater emphasis on high values. By contrast, the BC-GED approach approximates a formal likelihood function, and balances consideration of the high- and low- values. As a result, inferred results of the BC-GED method are more reasonable and consistent with the field survey results and previous related-studies. This method even discriminates the nonerodible characteristic of main channels.