The use of caesium‐137 to estimate agricultural erosion on steep slopes in a tropical watershed

The estimation of erosion and sediment delivery rates in tropical mountain watersheds is difficult and most of the methods widely used for estimating soil erosion over large areas have serious limitations. The ¹³⁷Cs approach has potential for quantifying soil erosion because it can provide retrospective estimates of long‐term (since 1963) net sediment redistribution relatively quickly. Despite its great potential, ¹³⁷Cs has not yet been used in an extensive, reconnaissance level survey of erosion in complex tropical mountain environments. The objective of this study was to examine the applicability of the ¹³⁷Cs method to estimate erosion on steep tropical agricultural lands (23 to 80% slopes) in the Nizao watershed, a humid, tropical mountain area of the Dominican Republic. In this study we (i) examine the variation of ¹³⁷Cs in ten reference sites—eight coffee groves and two forested sites—and (ii) estimate erosion from 14 cultivated fields. The soil pool of ¹³⁷Cs ranged from to 150 to 192 mBq cm⁻² on reference sites with minimal erosion. Variability among reference sites was less than expected for such complex mountain terrain. The variability within coffee and forested reference (average CV=28%) sites was similar to the variability found on grassland and forested reference sites in the temperate zone. The estimated annual soil loss from 14 sampled fields ranged from 6 to 61 t ha⁻¹ year⁻¹ with an overall mean of 26 t ha⁻¹ year⁻¹. Overall, the soil erosion estimates found using the ¹³⁷Cs method were much lower than those often assumed for such steep tropical hillsides. These erosion estimates account for soil loss since 1963 only and it seems likely that soil losses may have been much higher in earlier decades immediately after initial forest clearing earlier in the 20th century. Copyright © 2000 John Wiley & Sons, Ltd.     

Hydrograph separation using tracers and digital filters to quantify runoff components in a semi‐arid mesoscale catchment

Chemical hydrograph separation using electrical conductivity and digital filters is applied to quantify runoff components in the 1,640 km² semi‐arid Kaap River catchment and its subcatchments in South Africa. A rich data set of weekly to monthly water quality data ranging from 1978 to 2012 (450 to 940 samples per site) was analysed at 4 sampling locations in the catchment. The data were routinely collected by South Africa's national Department of Water and Sanitation, using standard sampling procedures. Chemical hydrograph separation using electrical conductivity (EC) as a tracer was used as reference and a recursive digital filter was then calibrated for the catchment. Results of the two‐component hydrograph separation indicate the dominance of baseflow in the low flow regime, with a contribution of about 90% of total flow; however, during the wet season, baseflow accounts for 50% of total flow. The digital filter parameters were very sensitive and required calibration, using chemical hydrograph separation as a reference. Calibrated baseflow estimates ranged from 40% of total flow at the catchment outlet to 70% in the tributaries. The study demonstrates that routinely monitored water quality data, especially EC, can be used as a meaningful tracer, which could also aid in the calibration of a digital filter method and reduce uncertainty of estimated flow components. This information enhances our understanding of how baseflow is generated and contributed to streamflow throughout the year, which can aid in quantification of environmental flows, as well as to better parameterize hydrological models used for water resources planning and management. Baseflow estimates can also be useful for groundwater and water quality management.     

Evaluation of the hydrology of the IBIS land surface model in a semi‐arid catchment

This paper evaluates the Integrated BIosphere Simulator (IBIS) land surface model using daily soil moisture data over a 3‐year period (2005–2007) at a semi‐arid site in southeastern Australia, the Stanley catchment, using the Monte Carlo generalized likelihood uncertainty estimation (GLUE) approach. The model was satisfactorily calibrated for both the surface 30 cm and full profile 90 cm. However, full‐profile calibration was not as good as that for the surface, which results from some deficiencies in the evapotranspiration component in IBIS. Relatively small differences in simulated soil moisture were associated with large discrepancies in the predictions of surface runoff, drainage and evapotranspiration. We conclude that while land surface schemes may be effective at simulating heat fluxes, they may be ineffective for prediction of hydrology unless the soil moisture is accurately estimated. Sensitivity analyses indicated that the soil moisture simulations were most sensitive to soil parameters, and the wilting point was the most identifiable parameter. Significant interactions existed between three soils parameters: porosity, saturated hydraulic conductivity and Campbell ‘b’ exponent, so they could not be identified independent of each other. There were no significant differences in parameter sensitivity and interaction for different hydroclimatic years. Even though the data record contained a very dry year and another year with a very large rainfall event, this indicated that the soil model could be calibrated without the data needing to explore the extreme range of dry and wet conditions. IBIS was much less sensitive to vegetation parameters. The leaf area index (LAI) could affect the mean of daily soil moisture time series when LAI < 1, while the variance of the soil moisture time series was sensitive to LAI > 1. IBIS was insensitive to the Jackson rooting parameter, suggesting that the effect of the rooting depth distribution on predictions of hydrology was insignificant. Copyright © 2014 John Wiley & Sons, Ltd.     

Soil moisture states,lateral flow,and streamflow generation in a semi‐arid,snowmelt‐driven catchment

Hydraulic connectivity on hillslopes and the existence of preferred soil moisture states in a catchment have important controls on runoff generation. In this study we investigate the relationships between soil moisture patterns, lateral hillslope flow, and streamflow generation in a semi‐arid, snowmelt‐driven catchment. We identify five soil moisture conditions that occur during a year and present a conceptual model based on field studies and computer simulations of how streamflow is generated with respect to the soil moisture conditions. The five soil moisture conditions are (1) a summer dry period, (2) a transitional fall wetting period, (3) a winter wet, low‐flux period, (4) a spring wet, high‐flux period, and (5) a transitional late‐spring drying period. Transitions between the periods are driven by changes in the water balance between rain, snow, snowmelt and evapotranspiration. Low rates of water input to the soil during the winter allow dry soil regions to persist at the soil–bedrock interface, which act as barriers to lateral flow. Once the dry‐soil flow barriers are wetted, whole‐slope hydraulic connectivity is established, lateral flow can occur, and upland soils are in direct connection with the near‐stream soil moisture. This whole‐slope connectivity can alter near‐stream hydraulics and modify the delivery of water, pressure, and solutes to the stream. Copyright © 2005 John Wiley & Sons, Ltd.     

Tunnel development over a 12 year period in a semi‐arid catchment of the Loess Plateau,China

The tunnel systems in a semi‐arid catchment of the Loess Plateau of China were repeatedly surveyed prior to the rainy seasons of 1989, 1999 and 2001. The surveys aimed to: (1) measure tunnel development over the 12 year period 1989–2001; (2) explore how the physiographical conditions affect the spatio‐temporal variability of tunnel development; and (3) to identify the geomorphic processes associated with tunnel development. The ultimate goal was to quantify the geomorphic significance of tunnel systems in the catchment. Over the 12 year period, the number of tunnel inlets was more than doubled and most of the newly increased inlets were initiated in the few catastrophic storm events. However, tunnel enlargement can occur in storm or inter‐storm periods, mainly through earth falls and slumps in inlets, and water erosion and roof cave‐in collapses in tunnel paths. Tunnel development varied with material properties, land uses and topographic conditions. Net tunnel erosion may contribute at least 25–30% of the catchment sediment yield and was mainly produced by the initiation and enlargement of tunnel inlets rather than tunnel paths. To protect the areas against tunnel erosion, terracing of the upper slopes seems to be more effective than planting vegetation on the lower slopes. Copyright © 2003 John Wiley & Sons, Ltd.     

Use of caesium‐137 data to evaluate SHETRAN simulated long‐term erosion patterns in arable lands

The caesium‐137 method of quantifying soil erosion is used to provide field data for validating the capability of the SHETRAN modelling system for predicting long‐term (30‐year) erosion rates and their spatial variability. Simulations were carried out for two arable farm sites (area 3–5 ha) in central England for which average annual erosion rates of 6·5 and 10·4 t ha^?1 year^?1 had already been determined using caesium‐137 measurements. These rates were compared with a range of simulated values representing the uncertainty in model output derived from uncertainty in the evaluation of model parameters. A successful validation was achieved in that the simulation range contained the measured rate at both sites, whereas the spatial variability was reproduced excellently at one site and partially at the other. The results indicate that, as the caesium‐137 technique measures the erosion caused by all the processes acting at a site, it is relevant to hydrologically based models such as SHETRAN only if erosion by wind, agricultural activities and other processes not represented in the model are insignificant. The results also indicate a need to reduce the uncertainty in model parameter evaluation. More generally, the caesium‐137 technique is shown to provide field data that improve the severity of the validation procedure (accounting for internal as well as outlet conditions) and that add spatial variability to magnitude as a condition for identifying unrealistic parameter sets when seeking to reduce simulation uncertainty. Copyright © 2004 John Wiley & Sons, Ltd.     

Hydrological response of a Brazilian semi‐arid catchment to different land use and climate change scenarios: a modelling study

Brazilian semi‐arid regions are characterized by water scarcity, vulnerability to desertification, and climate variability. The investigation of hydrological processes in this region is of major interest not only for water planning strategies but also to address the possible impact of future climate and land‐use changes on water resources. A hydrological distributed catchment‐scale model (DiCaSM) has been applied to simulate hydrological processes in a small representative catchment of the Brazilian northeast semi‐arid region, and also to investigate the impact of climate and land‐use changes, as well as changes associated with biofuel/energy crops production. The catchment is part of the Brazilian network for semi‐arid hydrology, established by the Brazilian Federal Government. Estimating and modelling streamflow (STF) and recharge in semi‐arid areas is a challenging task, mainly because of limitation in in situ measurements, and also due to the local nature of some processes. Direct recharge measurements are very difficult in semi‐arid catchments and contain a high level of uncertainty. The latter is usually addressed by short‐ and long‐time‐scale calibration and validation at catchment scale, as well as by examining the model sensitivity to the physical parameters responsible for the recharge. The DiCaSM model was run from 2000 to 2008, and streamflow was successfully simulated, with a Nash–Sutcliffe (NS) efficiency coefficient of 0·73, and R² of 0·79. On the basis of a range of climate change scenarios for the region, the DiCaSM model forecasted a reduction by 35%, 68%, and 77%, in groundwater recharge (GWR), and by 34%, 65%, and 72%, in streamflow, for the time spans 2010–2039, 2040–2069, and 2070–2099, respectively, could take place for a dry future climate scenario. These reductions would produce severe impact on water availability in the region. Introducing castor beans to the catchment would increase the GWR and streamflow, mainly if the caatinga areas would be converted into castor beans production. Changing an area of 1000 ha from caatinga to castor beans would increase the GWR by 46% and streamflow by 3%. If the same area of pasture is converted into castor beans, there would be an increase in GWR and streamflow by 24% and 5%, respectively. Such results are expected to contribute towards environmental policies for north‐east Brazil (NEB), and to biofuel production perspectives in the region. Copyright © 2010 John Wiley & Sons, Ltd.     

Spatial patterns and temporal stability of soil moisture across a range of scales in a semi‐arid environment

The antecedent soil moisture status of a catchment is an important factor in hydrological modelling. Traditional Hortonian infiltration models assume that the initial moisture content is constant across the whole catchment, despite the fact that even in small catchments antecedent soil moisture exhibits tremendous spatial heterogeneity. Spatial patterns of soil water distribution across three transects (two in a burnt area and one in an unburnt area) in a semi‐arid area were studied. At the transect scale, when the factors affecting soil moisture were limited to topographical position or local topography, spatial patterns showed time stability, but when other factors, such as vegetation, were taken into account, the spatial patterns became time unstable. At the point scale, and in the same areas, topographical position was the main factor controlling time stability. Scale dependence of time stability was studied and local topography and vegetation presence were observed to play an important role for the correlation between consecutive measures depending on the scale. Copyright © 2000 John Wiley & Sons, Ltd.     

River–aquifer interactions in a semi‐arid environment stressed by groundwater abstraction

Rivers and aquifers are, in many cases, a connected resource and as such the interactions between them need to be understood and quantified for the resource to be managed appropriately. The objective of this paper is to advance the understanding of river–aquifer interactions processes in semi‐arid environments stressed by groundwater abstraction. This is performed using data from a specific catchment where records of precipitation, evapotranspiration, river flow, groundwater levels and groundwater abstraction are analysed using basic statistics, hydrograph analysis and a simple mathematical model to determine the processes causing the spatial and temporal changes in river–aquifer interactions. This combined approach provides a novel but simple methodology to analyse river–aquifer interactions, which can be applied to catchments worldwide. The analysis revealed that the groundwater levels have declined (~ 3 m) since the onset of groundwater abstraction. The decline is predominantly due to the abstraction rather than climatic changes (r = 0.84 for the relationship between groundwater abstraction and groundwater levels; r = 0.92 for the relationship between decline in groundwater levels and magnitude of seasonal drawdown). It is then demonstrated that, since the onset of abstraction, the river has changed from being gaining to losing during low‐flow periods, defined as periods with flow less than 0.5, 1.0 or 1.5 GL/day (1 GL/day = 1 × 10⁶ m³/day). If defined as < 1.0 GL/day, low‐flow periods constitute approximately 65% of the river flows; the periods where the river is losing at low‐flow conditions are thus significant. Importantly, there was a significant delay (> 10 years) between the onset of groundwater abstraction and the changeover from gaining to losing conditions. Finally, a relationship between the groundwater gradient towards the river and the river flow at low‐flow is demonstrated. The results have important implications for water management as well as water ecology and quality. Copyright © 2012 John Wiley & Sons, Ltd.     

Predicting catchment flow in a semi‐arid region via an artificial neural network technique

A model of rainfall–runoff relationships is an essential tool in the process of evaluation of water resources projects. In this paper, we applied an artificial neural network (ANN) based model for flow prediction using the data for a catchment in a semi‐arid region in Morocco. Use of this method for non‐linear modelling has been demonstrated in several scientific fields such as biology, geology, chemistry and physics. The performance of the developed neural network‐based model was compared against multiple linear regression‐based model using the same observed data. It was found that the neural network model consistently gives superior predictions. Based on the results of this study, artificial neural network modelling appears to be a promising technique for the prediction of flow for catchments in semi‐arid regions. Accordingly, the neural network method can be applied to various hydrological systems where other models may be inappropriate. Copyright © 2004 John Wiley & Sons, Ltd.     

Groundwater–surface water interactions in a large semi‐arid floodplain: implications for salinity management

Flow regulation and water diversion for irrigation have considerably impacted the exchange of surface water between the Murray River and its floodplains. However, the way in which river regulation has impacted groundwater–surface water interactions is not completely understood, especially in regards to the salinization and accompanying vegetation dieback currently occurring in many of the floodplains. Groundwater–surface water interactions were studied over a 2 year period in the riparian area of a large floodplain (Hattah–Kulkyne, Victoria) using a combination of piezometric surface monitoring and environmental tracers (Cl⁻, δ²H, and δ¹⁸O). Despite being located in a local and regional groundwater discharge zone, the Murray River is a losing stream under low flow conditions at Hattah–Kulkyne. The discharge zone for local groundwater, regional groundwater and bank recharge is in the floodplain within ∼1 km of the river and is probably driven by high rates of transpiration by the riparian Eucalyptus camaldulensis woodland. Environmental tracers data suggest that the origin of groundwater is principally bank recharge in the riparian zone and a combination of diffuse rainfall recharge and localized floodwater recharge elsewhere in the floodplain. Although the Murray River was losing under low flows, bank discharge occurred during some flood recession periods. The way in which the water table responded to changes in river level was a function of the type of stream bank present, with point bars providing a better connection to the alluvial aquifer than the more common clay‐lined banks. Understanding the spatial variability in the hydraulic connection with the river channel and in vertical recharge following inundations will be critical to design effective salinity remediation strategies for large semi‐arid floodplains. Copyright © 2005 John Wiley & Sons, Ltd.     

Groundwater quality in the semi‐arid region of the Chahardouly basin,West Iran

Chahardouly basin is located in the western part of Iran and is characterized by semi‐arid climatic conditions and scarcity in water resources. The main aquifer systems are developed within alluvial deposits. The availability of groundwater is rather erratic owing to the occurrence of hard rock formation and a saline zone in some parts of the area. The aquifer systems of the area show signs of depletion, which have taken place in recent years due to a decline in water levels. Groundwater samples collected from shallow and deep wells were analysed to examine the quality characteristics of groundwater. The major ion chemistry of groundwater is dominated by Ca²⁺ and HCO₃^?, while higher values of total dissolved solids (TDS) in groundwater are associated with high concentrations of all major ions. An increase in salinity is recorded in the down‐gradient part of the basin. The occurrence of saline groundwater, as witnessed by the high electrical conductivity (EC), may be attributed to the long residence time of water and the dissolution of minerals, as well as evaporation of rainfall and irrigation return flow. Based on SAR values and sodium content (%Na), salinity appears to be responsible for the poor groundwater quality, rendering most of the samples not suitable for irrigation use. Copyright © 2007 John Wiley & Sons, Ltd.     

Modelling surface runoff to evaluate the effects of wildfires in multiple semi‐arid,shrubland‐dominated catchments

Wildfires change the infiltration properties of soil, reduce the amount of interception and result in increased runoff. A wildfire at Northeast Attica, Central Greece, in August 2009, destroyed approximately one third of a study area consisting of a mixture of shrublands, pastures and pines. The present study simultaneously models multiple semi‐arid, shrubland‐dominated Mediterranean catchments and assesses the hydrological response (mean annual and monthly runoff and runoff coefficients) during the first few years following wildfires. A physically based, hydrological model (MIKE SHE) was chosen. Calibration and validation results of mean monthly discharge presented very good agreement with the observed data for the pre‐wildfire and post‐wildfire period for two subcatchments (Nash–Sutcliffe Efficiency coefficient of 79.7%). The model was then used to assess the pre‐wildfire and post‐wildfire runoff responses for each of seven catchments in the study area. Mean annual surface runoff increased for the first year and after the second year following the wildfires increased by 112% and 166%, respectively. These values are within the range observed in similar cases of monitored sites. This modelling approach may provide a way of prioritizing catchment selection with respect to post‐fire remediation activities. Additionally, this modelling assessment methodology would be valuable to other semi‐arid areas because it provides an important means for comprehensively assessing post‐wildfire response over large regions and therefore attempts to address some of the scaled issues in the specific literature field of research. Copyright © 2015 John Wiley & Sons, Ltd.     

Evaporation of perennial semi‐arid woodland in southeastern Australia is adapted for irregular but common dry periods

Measurements of water vapour flux from semi‐arid perennial woodland (mallee) were made for 3 years using eddy covariance instrumentation. There have been no previous long‐term, detailed measures of water use in this ecosystem. Latent energy flux (LE) on a half hourly basis was the measure of the combined soil and plant evaporation, ‘evapotranspiration’ (E_LE) of the site. Aggregation over 3 years of the site measured rain (1136 mm) and the estimated evaporation (794 mm) suggests that 342 mm or 30% of rain had moved into or past the root zone of the vegetation. Above average rainfall during 2011 and the first quarter of 2012 (633 mm, 15 months) would likely have been the period during which significant groundwater recharge occurred. At times immediately after rainfall, E_LE rates were the same or exceeded estimates of potential E calculated from a suitably parameterized Penman–Monteith (E_PMo) equation. Apparent free water E from plant interception and soil evaporation was about 2.3 mm and lasted for 1.3 days following rainfall in summer, while in autumn, E was 5.1 mm that lasted over 5.4 days. The leaf area index (LAI) needed to adjust a wind function calibrated Penman equation (E_PMe) to match the E_LE values could be back calculated to generate seasonal change in LAI from 0.12 to 0.46 and compared well with normalized difference vegetation index; r = 0.38 and p = 0.0213* and LAI calculated from digital cover photography. The apparently conservative response of perennial vegetation evaporation to available water in these semi‐arid environments reinforces the conclusion that these ecosystems use this mechanism to survive the reasonably common dry periods. Plant response to soil water availability is primarily through gradual changes in leaf area. Copyright © 2015 John Wiley & Sons, Ltd.     

Evaluating the use of a gridded climate surface for modelling groundwater recharge in a semi‐arid region (Okanagan Basin,Canada)

Spatially distributed groundwater recharge was simulated for a segment of a semi‐arid valley using three different treatments of meteorological input data and potential evapotranspiration (PET). For the same area, timeframe, land cover characteristics and soil properties, groundwater recharge was estimate using (i) single‐station climate data with monthly PET calculated by the Thornthwaite method; (ii) single‐station climate data with daily PET calculated by the Penman–Monteith method; and (iii) daily gridded climate data with spatially distributed PET calculated using the Penman–Monteith method. For each treatment, the magnitude and distribution of actual evapotranspiration (AET) for summer months compared well with those estimated for a 5‐year crop study, suggesting that the near‐surface hydrological processes were replicated and that subsequent groundwater recharge rates are realistic. However, for winter months, calculated AET was near zero when using the Thornthwaite PET method. Mean annual groundwater recharge varied from ～3·2 to 10·0 mm when PET was calculated by the Thornthwaite method, and from ～1·8 to 7·5 mm when PET was calculated by the Penman–Monteith method. Comparisons of bivariate plots of seasonal recharge rates estimated from single‐station versus gridded surface climate reveal that there is greater variability between the different methods for spring months, which is the season of greatest recharge. Furthermore, these seasonal differences are shown to provide different results when compared to the depth to water table, which could lead to different results of evaporative extinction depth. These findings illustrate potential consequences of using different approaches for representing spatial meteorological input data, which could provide conflicting predictions when modelling the influence of climate change on groundwater recharge. Copyright © 2010 John Wiley & Sons, Ltd.     

Gully‐head erosion processes on a semi‐arid valley floor in Kenya: a case study into temporal variation and sediment budgeting

A three year monitoring programme of gully‐head retreat was established to assess the significance of sediment production in a drainage network that expanded rapidly by gully‐head erosion on the low‐angled alluvio‐lacustrine Njemps Flats in semi‐arid Baringo District, Kenya. This paper discusses the factors controlling the large observed spatial and temporal variation in gully‐head retreat rates, ranging from 0 to 15 m a^?1. The selected gullies differed in planform and in runoff‐contributing catchment area but soil material and land use were similar. The data were analysed at event and annual timescales. The results show that at annual timescale rainfall amount appears to be a good indicator of gully‐head retreat, while at storm‐event timescale rainfall distribution has to be taken into account. A model is proposed, including only rainfall (P) and the number of dry days (DD) between storms: which explains 56 per cent of the variation in retreat rate of the single‐headed gully of Lam1. A detailed sediment budget has been established for Lam1 and its runoff‐contributing area (RCA). By measuring sediment input from the RCA, the sediment output by channelized flow and linear retreat of the gully head for nine storms, it can be seen that erosion shifts between different components of the budget depending on the duration of the dry period (DD) between storms. Sediment input from the RCA was usually the largest component for the smaller storms. The erosion of the gully head occurred as a direct effect of runoff falling over the edge (GH_waterfall) and of the indirect destabilization of the adjacent walls by the waterfall erosion and by saturation (GH_mass/storage). The latter component (GH_mass/storage) was usually much larger that the former (GH_waterfall). The sediment output from the gully was strongly related to the runoff volume while the linear retreat, because of its complex behaviour, was not. Overall, the results show that the annual retreat is the optimal timescale to predict retreat patterns. More detailed knowledge about relevant processes and interactions is necessary if gully‐head erosion is to be included in event‐based soil erosion models. Copyright © 2001 John Wiley & Sons, Ltd.     

Use of magnetic and radiometric measurements to investigate erosion and sedimentation in a British upland catchment

In response to the potential shortcomings of single-technique strategies in the investigation of erosion and sedimentation, a combined magnetic and radiometric (¹³Cs and ²¹⁰Pb) approach has been undertaken in the upland, watershed-lake system of Howden Reservoir, Derbyshire. By combining these techniques, some assessment of sediment sources and the erosion status of the catchment has been achieved. Alone, each approach would have been unable to determine unequivocally sediment provenance. Furthermore, the parallel use of these measurements has highlighted limitations and/or uncertainties in both the magnetic and ¹³⁷Cs techniques. These problems reflect the particular soil characteristics and drainage conditions of this upland catchment. Despite the documented severity of peat erosion in the region, Howden Reservoir has a mixture of sediment sources and a relatively moderate rate of sedimentation. Sediment yields (total 127·7 t km^?2 yr^?1 including organic fraction 31·3 t km^?2 yr^?1) are, however, higher than in other British upland areas.     

The influence of preferential flow on hillslope hydrology in a semi‐arid watershed (in the Spanish Dehesas)

Preferential flow is known to influence hillslope hydrology in many areas around the world. Most research on preferential flow has been performed in temperate regions. Preferential infiltration has also been found in semi‐arid regions, but its impact on the hydrology of these regions is poorly known. The aim of this study is to describe and quantify the influence of preferential flow on the hillslope hydrology from small scale (infiltration) to large scale (subsurface stormflow) in a semi‐arid Dehesa landscape. Precipitation, soil moisture content, piezometric water level and discharge data were used to analyse the hydrological functioning of a catchment in Spain. Variability of soil moisture content during the transition from dry to wet season (September to November) within horizontal soil layers leads to the conclusion that there is preferential infiltration into the soils. When the rainfall intensity is high, a water level rapidly builds up in the piezometer pipes in the area, sometimes even reaching soil surface. This water level also drops back to bedrock within a few hours (under dry catchment conditions) to days (under wet catchment conditions). As the soil matrix is not necessarily wet while this water layer is built up, it is thought to be a transient water table in large connected pores which drain partly to the matrix, partly fill up bedrock irregularities and partly drain through subsurface flow to the channels. When the soil matrix becomes wetter the loss of water from macropores to the matrix and bedrock decreases and subsurface stormflow increases. It may be concluded that the hillslope hydrological system consists of a fine matrix domain and a macropore domain, which have their own flow characteristics but which also interact, depending on the soil matrix and macropore moisture contents. The macropore flow can result in subsurface flow, ranging from 13% contribution to total discharge for a large event of high intensity rainfall or high discharge to 80% of total discharge for a small event with low intensity rainfall or low discharge. During large events the fraction of subsurface stormflow in the discharge is suppressed by the large amount of surface runoff. Copyright © 2008 John Wiley & Sons, Ltd.     

The importance of surface controls on overland flow connectivity in semi‐arid environments: results from a numerical experimental approach

In semi‐arid environments, the characteristics of the land surface determine how rainfall is transformed into surface runoff and influences how this runoff moves from the hillslopes into river channels. Whether or not water reaches the river channel is determined by the hydrological connectivity. This paper uses a numerical experiment‐based approach to systematically assess the effects of slope length, gradient, flow path convergence, infiltration rates and vegetation patterns on the generation and connectivity of runoff. The experiments were performed with the Connectivity of Runoff Model, 2D version distributed, physically based, hydrological model. The experiments presented are set within a semi‐arid environment, characteristic of south‐eastern Spain, which is subject to low frequency high rainfall intensity storm events. As a result, the dominant hydrological processes are infiltration excess runoff generation and surface flow dynamics. The results from the modelling experiments demonstrate that three surface factors are important in determining the form of the discharge hydrograph: the slope length, the slope gradient and the infiltration characteristics at the hillslope‐channel connection. These factors are all related to the time required for generated runoff to reach an efficient flow channel, because once in this channel, the transmission losses significantly decrease. Because these factors are distributed across the landscape, they have a fundamental role in controlling the landscape hydrological response to storm events. Copyright © 2013 John Wiley & Sons, Ltd.