Modelling water-harvesting systems in the arid south of Tunisia using SWAT

In many arid countries, runoff water-harvesting systems support the livelihood of the rural population. Little is known, however, about the effect of these systems on the water balance components of arid watersheds. The objective of this study was to adapt and evaluate the GIS-based watershed model SWAT (Soil Water Assessment Tool) for simulating the main hydrologic processes in arid environments. The model was applied to the 270-km(2) watershed of wadi Koutine in southeast Tunisia, which receives about 200 mm annual rain. The main adjustment for adapting the model to this dry Mediterranean environment was the inclusion of water-harvesting systems, which capture and use surface runoff for crop production in upstream subbasins, and a modification of the crop growth processes. The adjusted version of the model was named SWAT-WH. Model evaluation was performed based on 38 runoff events recorded at the Koutine station between 1973 and 1985. The model predicted that the average annual watershed rainfall of the 12-year evaluation period (209 mm) was split into ET (72%), groundwater recharge (22%) and outflow (6%). The evaluation coefficients for calibration and validation were, respectively, R-2 (coefficient of determination) 0.77 and 0.44; E (Nash-Sutcliffe coefficient) 0.73 and 0.43; and MAE (Mean Absolute Error) 2.6 mm and 3.0 mm, indicating that the model could reproduce the observed events reasonably well. However, the runoff record was dominated by two extreme events, which had a strong effect on the evaluation criteria. Discrepancies remained mainly due to uncertainties in the observed daily rainfall and runoff data. Recommendations for future research include the installation of additional rainfall and runoff gauges with continuous data logging and the collection of more field data to represent the soils and land use. In addition, crop growth and yield monitoring is needed for a proper evaluation of crop production, to allow an economic assessment of the different water uses in the watershed

Effects of watershed and riparian zone characteristics on nutrient concentrations in the River Scheldt Basin

International audienceThe relative influence of a set of watershed characteristics on surface water nutrient concentrations was examined in 173 watersheds within two subcatchments (Upper-Scheldt and Nete) of the River Scheldt Basin (Flanders, Belgium). Each watershed was described by seasonal rainfall, discharge loading of point sources, morphological characteristics (area, average slope, drainage density, elongation), land use and soil properties (soil texture and drainage). Partial regression analysis revealed that soil drainage variables had the strongest influence on nutrient concentrations. Additional influence was exerted by land use and point source loading variables. Nitrate concentrations were positively correlated with effluent loadings coming from wastewater treatment plants and with the area of agricultural land. Phosphate concentrations were best explained by effluent loadings of industrial point sources and by the area of urban land. Land use close to the river was not a better predictor of nitrate and phosphate concentrations than land use away from the river. This suggests that the mediating impact of riparian zones is rather explained by the hydrologic pathways within the buffer strip

Evaluating parameterizations of aerodynamic resistance to heat transfer using field measurements

Parameterizations of aerodynamic resistance to heat and water transfer have a significant impact on the accuracy of models of land – atmosphere interactions and of estimated surface fluxes using spectro-radiometric data collected from aircrafts and satellites. We have used measurements from an eddy correlation system to derive the aerodynamic resistance to heat transfer over a bare soil surface as well as over a maize canopy. Diurnal variations of aerodynamic resistance have been analyzed. The results showed that the diurnal variation of aerodynamic resistance during daytime (07:00 h–18:00 h) was significant for both the bare soil surface and the maize canopy although the range of variation was limited. Based on the measurements made by the eddy correlation system, a comprehensive evaluation of eight popularly used parameterization schemes of aerodynamic resistance was carried out. The roughness length for heat transfer is a crucial parameter in the estimation of aerodynamic resistance to heat transfer and can neither be taken as a constant nor be neglected. Comparing with the measurements, the parameterizations by Choudhury et al. (1986), Viney (1991), Yang et al. (2001) and the modified forms of Verma et al. (1976) and Mahrt and Ek (1984) by inclusion of roughness length for heat transfer gave good agreements with the measurements, while the parameterizations by Hatfield et al. (1983) and Xie (1988) showed larger errors even though the roughness length for heat transfer has been taken into account

Modelling direct episodic recharge in the Western Australian wheatbelt

In agricultural regions of Western Australia, salinity is spreading. This is because the area taken up by groundwater discharge is increasing as a result of increased groundwater recharge following the replacement of native vegetation systems by annual crops and pasture species. Attempts to reduce groundwater recharge are now being made as it is hoped that this will decrease the rate of land salinisation. At several sites, average recharge rates have been estimated to be from 2% to 13% of the average annu ter hydrographs from the agricultural regions indicate that at some sites recharge does not occur as small amounts every year in a regular manner, but as infrequent, unpredictable, relatively large ©·pisodic\u27 events. At such sites, relatively small incr ses in the water use of annual crops and patures are unlikely to have a significant effect on the magnitude of large episodic pulses of recharge

Runoff thresholds in derived flood frequency distributions

In general, different mechanisms may be identified as responsible of runoff generation during ordinary events or extraordinary events at the basin scale. In a simplified scheme these mechanisms may be represented by different runoff thresholds. In this context, the derived flood frequency model, based on the effect of partial contributing areas on peak flow, proposed by Iacobellis and Fiorentino (2000), was generalized by providing a new formulation of the derived distribution where two runoff components are explicitly considered. The model was tested on a group of basins in Southern Italy characterized by annual maximum flood distributions highly skewed. The application of the proposed model provided good results in terms of descriptive ability. Model parameters were also found to be well correlated with geomorphological basin descriptors. Two different threshold mechanisms, associated respectively to ordinary and extraordinary events, were identified. In fact, we found that ordinary floods are mostly due to rainfall events exceeding a threshold infiltration rate in a small source area, while the so-called outlier events, responsible of the high skewness of flood distributions, are triggered when severe rainfalls exceed a threshold storage in a large portion of the basin

Runoff regime estimation at high-elevation sites: a parsimonious water balance approach

Abstract. We develop a water balance model, parsimonious both in terms of parameterization and of required input data, to characterize the average runoff regime of high-elevation and scarcely monitored basins. The model uses a temperature threshold to partition precipitation into rainfall and snowfall, and to estimate evapotranspiration volumes. The role of snow in the transformation of precipitation into runoff is investigated at the monthly time scale through a specific snowmelt module that estimates melted quantities by a non-linear function of temperature. A probabilistic representation of temperature is also introduced, in order to mimic its sub-monthly variability. To account for the commonly reported rainfall underestimation at high elevations, a two-step precipitation adjustment procedure is implemented to guarantee the closure of the water balance. The model is applied to a group of catchments in the North-Western Italian Alps, and its performances are assessed by comparing measured and simulated runoff regimes both in terms of total bias and anomalies, by means of a new metric, specifically conceived to compare the shape of the two curves. The obtained results indicates that the model is able to predict the observed runoff seasonality satisfactorily, notwithstanding its parsimony (the model has only two parameters to be estimated). In particular, when the parameter calibration is performed separately for each basin, the model proves to be able to reproduce the runoff seasonality. At the regional scale (i.e., with uniform parameters for the whole region), the performance is less positive, but the model is still able to discern among different mechanisms of runoff formation that depend on the role of the snow storage. Because of its parsimony and the robustness in the approach, the model is suitable for application in ungauged basins and for large scale investigations of the role of climatic variables on water availability and runoff timing in mountainous regions

Uncertainty in geological and hydrogeological data

Uncertainty in conceptual model structure and in environmental data is of essential interest when dealing with uncertainty in water resources management. To make quantification of uncertainty possible is it necessary to identify and characterise the uncertainty in geological and hydrogeological data. This paper discusses a range of available techniques to describe the uncertainty related to geological model structure and scale of support. Literature examples on uncertainty in hydrogeological variables such as saturated hydraulic conductivity, specific yield, specific storage, effective porosity and dispersivity are given. Field data usually have a spatial and temporal scale of support that is different from the one on which numerical models for water resources management operate. Uncertainty in hydrogeological data variables is characterised and assessed within the methodological framework of the HarmoniRiB classification

Recent trends in groundwater levels in a highly seasonal hydrological system: the Ganges-Brahmaputra-Meghna Delta

Groundwater levels in shallow aquifers underlying Asian mega-deltas are characterized by strong seasonal variations associated with monsoon rainfall. To resolve trend and seasonal components in weekly groundwater levels in the Ganges-Brahmaputra-Meghna (GBM) Delta, we apply a nonparametric seasonal-trend decomposition procedure (STL) to observations compiled from 1985-2005 in Bangladesh. Seasonality dominates observed variance in groundwater levels but declining groundwater levels (> 1 m/yr) are detected in urban and peri-urban areas around Dhaka as well as in north-central, northwestern, and southwestern parts of the country (0.1-0.5 m/yr) where intensive abstraction of groundwater is conducted for dry-season rice cultivation. Rising groundwater levels (0.5-2.5 cm/yr) are observed in the estuarine and southern coastal regions. This novel application of the STL procedure reveals, for the first time, the unsustainability of irrigation supplied by shallow aquifers in some areas (e.g., High Barind Tract) of the GBM Delta and the hydrological impact of potential seawater intrusion of coastal aquifers associated with sea-level rise. Our findings provide important insight into the hydrological impacts of groundwater-fed irrigation and sea-level rise in other Asian mega-deltas where monitoring data are limited

An approach to propagate streamflow statistics along the river network

Streamflow at ungauged sites is often predicted by means of regional statistical procedures. The standard regional approaches do not preserve the information related to the hierarchy among gauged stations deriving from their location along the river network. However, this information is important when estimating runoff at a site located immediately upstream or downstream of a gauging station. We propose here a novel approach, referred to as the Along-Stream Estimation (ASE) method, to improve runoff estimation at ungauged sites. The ASE approach starts from the regional estimate at an ungauged (target) site, and corrects it based on regional and sample estimates of the same variable at a donor site, where sample data are available. A criterion to define the domain of application around each donor site of the ASE approach is proposed, and the uncertainty inherent in the estimates obtained is evaluated. This allows one to compare the variance of the along-stream estimates to that of other models that eventually become available for application (e.g. regional models), and thus to choose the most accurate method (or to combine different estimates). The ASE model was applied in the northwest of Italy in connection with an existing regional model for flood frequency analysis. The analysed variables are the first L-moments of the annual discharge maxima. The application demonstrates that the ASE approach can be used effectively to improve the regional estimates for the L-moment of order one (the index flood), particularly when the area ratio of a pair of donor-target basins is less than or equal to ten. However, in this case study, the method does not provide significant improvements to the estimation of higher-order L-moment