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

Past terrestrial water storage (1980–2008) in the Amazon Basin reconstructed from GRACE and in situ river gauging data

Terrestrial water storage (TWS) composed of surface waters, soil moisture, groundwater and snow where appropriate, is a key element of global and continental water cycle. Since 2002, the Gravity Recovery and Climate Experiment (GRACE) space gravimetry mission provides a new tool to measure large-scale TWS variations. However, for the past few decades, direct estimate of TWS variability is accessible from hydrological modeling only. Here we propose a novel approach that combines GRACE-based TWS spatial patterns with multi-decadal-long in situ river level records, to reconstruct past 2-D TWS over a river basin. Results are presented for the Amazon Basin for the period 1980–2008, focusing on the interannual time scale. Results are compared with past TWS estimated by the global hydrological model ISBA-TRIP. Correlations between reconstructed past interannual TWS variability and known climate forcing modes over the region (e.g., El Niño-Southern Oscillation and Pacific Decadal Oscillation) are also estimated. This method offers new perspective for improving our knowledge of past interannual TWS in world river basins where natural climate variability (as opposed to direct anthropogenic forcing) drives TWS variations