Spatial and temporal variability of the hydrological response in a small Mediterranean research catchment (Vallcebre,Eastern Pyrenees)

This paper analyses the spatial and temporal variability of the hydrological response in a small Mediterranean catchment (Cal Rodó). The first part of the analysis focuses on the rainfall–runoff relationship at seasonal and monthly scale, using an 8‐year data set. Then, using storm‐flow volume and coefficient, the temporal variability of the rainfall–runoff relationship and its relationship with several hydrological variables are analysed at the event scale from hydrographs observed over a 3‐year period. Finally, the spatial non‐linearity of the hydrological response is examined by comparing the Cal Rodó hydrological response with the Can Vila sub‐catchment response at the event scale. Results show that, on a seasonal and monthly scale, there is no simple relationship between rainfall and runoff depths, and that evapotranspiration is a factor that introduced some non‐linearity in the rainfall–runoff relationship. The analysis of monthly values also reveals the existence of a threshold in the relationship between rainfall and runoff depths, denoting a more contrasted hydrological response than the one usually observed in humid catchments. At the event scale, the storm‐flow coefficient has a clear seasonal pattern with an alternance between a wet period, when the catchment is hydrologically responsive, and a dry summer period, when the catchment is much less reactive to any rainfall. The relationship between the storm‐flow coefficient and rainfall depth, rainfall maximum intensity and base‐flow shows that observed correlations are the same as those observed for humid conditions, even if correlation coefficients are notably lower. Comparison with the Can Vila sub‐catchment highlights the spatial heterogeneity of the rainfall‐runoff relationship at the small catchment scale. Although interpretation in terms of runoff processes remains delicate, heterogeneities between the two catchments seem to be related to changes in the ratio between infiltration excess and saturation processes in runoff formation. Copyright © 2007 John Wiley & Sons, Ltd.     

Spatial and temporal variability of water table dynamics in an afforested catchment of the Central Spanish Pyrenees

Human-induced afforestation has been one of the main policies for environmental management of farmland abandonment in Mediterranean areas. Over the last decades, several studies have reviewed the impact of afforestation activities on geomorphological and hydrological responses and soil properties, although few studies have evaluated the effects on water table dynamics. In parallel to human-induced afforestation activities, natural revegetation occurred in abandoned fields and in fields where the intensity of human activity declined, driving the expansion of shrubs. This research addresses the spatial and temporal variability of water table dynamics in a small afforested sub-catchment located in the Central Spanish Pyrenees. Differences between afforestation (Pinus nigra and Pinus sylvestris) and natural plant colonization (shrubs, mainly Genista scorpius, Buxus sempervirens, and Juniperus communis) and early abandoned meadows (G. scorpius), are analysed in terms of runoff generation and seasonal water table depth dynamics. Precipitation, runoff and water table datasets recorded for the 2014–2019 period are used. Results show a high temporal and spatial variability with large fluctuations in discharge and water table. Groundwater dynamics varied markedly over the year, identifying a wet and dry period with different responses suggesting different runoff generation processes (Hortonian flow during dry and wet periods, and saturation excess runoff during wet conditions). Furthermore, important differences are noted among the various land cover types: (i) in the natural revegetation area (shrubland and meadows) a marked seasonal cycle was observed with short saturation periods during winter and spring; and (ii) in the afforestation areas, the water table dynamics showed a seasonal cycle with a high variability, with fast responses and rapid oscillations. Likewise, the relationship between the depth of water table and hydrological variables was not straightforward, suggesting complex hydrological behaviour.     

Field investigation and modelling of coupled stream discharge and shallow water‐table dynamics in a small Mediterranean catchment (Sardinia)

In the semi‐arid Mediterranean environment, the rainfall–runoff relationships are complex because of the markedly irregular patterns in rainfall, the seasonal mismatch between evaporation and rainfall, and the spatial heterogeneity in landscape properties. Watersheds often display considerable non‐linear threshold behavior, which still make runoff generation an open research question. Our objectives in this context were: to identify the primary processes of runoff generation in a small natural catchment; to test whether a physically based model, which takes into consideration only the primary processes, is able to predict spatially distributed water‐table and stream discharge dynamics; and to use the hydrological model to increase our understanding of runoff generation mechanisms. The observed seasonal dynamics of soil moisture, water‐table depth, and stream discharge indicated that Hortonian overland‐flow was negligible and the main mechanism of runoff generation was saturated subsurface‐flow. This gives rise to base‐flow, controls the formation of the saturated areas, and contributes to storm‐flow together with saturation overland‐flow. The distributed model, with a 1D scheme for the kinematic surface‐flow, a 2D sub‐horizontal scheme for the saturated subsurface‐flow, and ignoring the unsaturated flow, performed efficiently in years when runoff volume was high and medium, although there was a smoothing effect on the observed water‐table. In dry years, small errors greatly reduced the efficiency of the model. The hydrological model has allowed to relate the runoff generation mechanisms with the land‐use. The forested hillslopes, where the calibrated soil conductivity was high, were never saturated, except at the foot of the slopes, where exfiltration of saturated subsurface‐flow contributed to storm‐flow. Saturation overland‐flow was only found near the streams, except when there were storm‐flow peaks, when it also occurred on hillslopes used for pasture, where soil conductivity was low. The bedrock–soil percolation, simulated by a threshold mechanism, further increased the non‐linearity of the rainfall–runoff processes. Copyright © 2013 John Wiley & Sons, Ltd.     

Spatial and temporal variability of event runoff characteristics in a small agricultural catchment

ABSTRACT

The objective of this study is to investigate the factors that control event runoff characteristics at the small catchment scale. The study area is the Hydrological Open Air Laboratory, Lower Austria. Event runoff coefficient (Rc), recession time constant (Tc) and peak discharge (Qp) are estimated from hourly discharge and precipitation data for 298 events in the period 2013–2015. The results show that the Rc and their variability tend to be largest for the tile drainages (mean Rc = 0.09) and the main outlet (mean Rc = 0.08) showing larger Rc in January/February and smaller Rc in July/August. Tc does not vary much between the systems and tends to be largest at the main outlet (mean Tc = 6.5 h) and smallest for the tile drainages (mean Tc = 4.5 h). Groundwater levels explain the temporal variability of Rc and Tc more than soil moisture or precipitation, suggesting a role of shallow flow paths.     

Spatio‐temporal controls of snowmelt and runoff generation during rain‐on‐snow events in a mid‐latitude mountain catchment

A network of 30 standalone snow monitoring stations was used to investigate the snow cover distribution, snowmelt dynamics, and runoff generation during two rain‐on‐snow (ROS) events in a 40 km² montane catchment in the Black Forest region of southwestern Germany. A multiple linear regression analysis using elevation, aspect, and land cover as predictors for the snow water equivalent (SWE) distribution within the catchment was applied on an hourly basis for two significant ROS flood events that occurred in December 2012. The available snowmelt water, liquid precipitation, as well as the total retention storage of the snow cover were considered in order to estimate the amount of water potentially available for the runoff generation. The study provides a spatially and temporally distributed picture of how the two observed ROS floods developed in the catchment. It became evident that the retention capacity of the snow cover is a crucial mechanism during ROS. It took several hours before water was released from the snowpack during the first ROS event, while retention storage was exceeded within 1 h from the start of the second event. Elevation was the most important terrain feature. South‐facing terrain contributed more water for runoff than north‐facing slopes, and only slightly more runoff was generated at open compared to forested areas. The results highlight the importance of snowmelt together with liquid precipitation for the generation of flood runoff during ROS and the large temporal and spatial variability of the relevant processes. Copyright © 2015 John Wiley & Sons, Ltd.     

Real‐time monitoring and estimation of the discharge of flash floods in a steep mountain catchment

Synchronously and accurately estimating the flood discharges and dynamic changes in the fluid density is essential for hydraulic analysis and forecasting of flash floods, as well as for risk assessment. However, such information is rare for steep mountain catchments, especially in regions that are hotspots for earthquakes. Therefore, six hydrological monitoring sites were established in the main stream and tributaries of the 78.3‐km² Longxi River catchment, an affected region of the Wenchuan earthquake region in China. Direct real‐time monitoring equipment was installed to measure the flow depths, velocities, and fluid total pressures of the flood hydrographs. On the basis of field measurements, real‐time mean cross‐sectional velocities during the flood hydrographs could be derived from easily obtainable parameters: cross‐sectional maximum velocities and the calibrated dimensionless parameter K_h . Real‐time discharges were determined on the basis of a noncontact method to establish the effective rating curves of this mountainous stream, ranging from 1.46 to 386.34 m³/s with the root mean square errors of ≤10.22 m³/s. Compared with the traditional point‐velocity method and empirical Manning's formula, the proposed noncontact method was reliable and safe for monitoring whole flood hydrographs. Additionally, the real‐time fluid density during the flood hydrographs was calculated on the basis of the direct monitoring parameters for fluid total pressures and water depths. During the flood hydrograph, transient flow behaviour with higher fluid density generally occurred downstream during the flood peak periods when the flow was in the supercritical flow regime. The observed behaviour greatly increased the threat of damage to infrastructure and human life near the river. Thus, it is important to accurately estimate flood discharge and identify for fluid densities so that people at risk from an impending flash flood are given reliable, advanced warning.     

Seasonal changes in runoff generation in a small forested mountain catchment

This study aimed to investigate the seasonal variability of runoff generation processes, the sources of stream water, and the controls on the contribution of event water to streamflow for a small forested catchment in the Italian pre‐Alps. Hydrometric, isotopic, and electrical conductivity data collected between August 2012 and August 2013 revealed a marked seasonal variability in runoff responses. Noticeable differences in runoff coefficients and hydrological dynamics between summer and fall/spring rainfall events were related to antecedent moisture conditions and event size. Two‐component and three‐component hydrograph separation and end‐member mixing analysis showed an increase in event water contributions to streamflow with event size and average rainfall intensity. Event water fractions were larger during dry conditions in the summer, suggesting that stormflow generation in the summer consisted predominantly of direct channel precipitation and some saturated overland flow from the riparian zone. On the contrary, groundwater and hillslope soil water contributions dominated the streamflow response during wet conditions in fall. Seasonal differences were also noted between event water fractions computed based on isotopic and electrical conductivity data, likely because of the dilution effect during the wetter months. Copyright © 2014 John Wiley & Sons, Ltd.     

Assessing the long-term impact of climatic variability and human activities on the water resources of a meso-scale Mediterranean catchment

Abstract

This article addresses the critical need for a better quantitative understanding of how water resources from the Hérault River catchment in France have been influenced by climate variability and the increasing pressure of human activity over the last 50 years. A method is proposed for assessing the relative impacts of climate and growing water demand on the decrease in discharge observed at various gauging stations in the periods 1961–1980 and 1981–2010. An annual water balance at the basin scale was calculated first, taking into account precipitation, actual evapotranspiration, water withdrawals and water discharge. Next, the evolution of the seasonal variability in hydroclimatic conditions and water withdrawals was studied. The catchment was then divided into zones according to the main geographical characteristics to investigate the heterogeneity of the climatic and human dynamics. This delimitation took into account the distribution of climate, topography, lithology, land cover and water uses, as well as the availability of discharge series. At the area scale, annual water balances were calculated to understand the internal changes that occurred in the catchment between both past periods. The decrease in runoff can be explained by the decrease in winter precipitation in the upstream areas and by the increase during summer in both water withdrawals and evapotranspiration in the downstream areas, mainly due to the increase in temperature. Thus, water stress increased in summer by 35%. This work is the first step of a larger research project to assess possible future changes in the capacity to satisfy water demand in the Hérault River catchment, using a model that combines hydrological processes and water demand.

Editor Z.W. Kundzewicz     

Snow model sensitivity analysis to understand spatial and temporal snow dynamics in a high‐elevation catchment

In this paper, we addressed a sensitivity analysis of the snow module of the GEOtop2.0 model at point and catchment scale in a small high‐elevation catchment in the Eastern Italian Alps (catchment size: 61 km²). Simulated snow depth and snow water equivalent at the point scale were compared with measured data at four locations from 2009 to 2013. At the catchment scale, simulated snow‐covered area (SCA) was compared with binary snow cover maps derived from moderate‐resolution imaging spectroradiometer (MODIS) and Landsat satellite imagery. Sensitivity analyses were used to assess the effect of different model parameterizations on model performance at both scales and the effect of different thresholds of simulated snow depth on the agreement with MODIS data. Our results at point scale indicated that modifying only the “snow correction factor” resulted in substantial improvements of the snow model and effectively compensated inaccurate winter precipitation by enhancing snow accumulation. SCA inaccuracies at catchment scale during accumulation and melt period were affected little by different snow depth thresholds when using calibrated winter precipitation from point scale. However, inaccuracies were strongly controlled by topographic characteristics and model parameterizations driving snow albedo (“snow ageing coefficient” and “extinction of snow albedo”) during accumulation and melt period. Although highest accuracies (overall accuracy = 1 in 86% of the catchment area) were observed during winter, lower accuracies (overall accuracy < 0.7) occurred during the early accumulation and melt period (in 29% and 23%, respectively), mostly present in areas with grassland and forest, slopes of 20–40°, areas exposed NW or areas with a topographic roughness index of ?0.25 to 0 m. These findings may give recommendations for defining more effective model parameterization strategies and guide future work, in which simulated and MODIS SCA may be combined to generate improved products for SCA monitoring in Alpine catchments.     

Spatial variability of rainfall on a sub‐kilometre scale

The variability of rainfall in space and time is an essential driver of many processes in nature but little is known about its extent on the sub‐kilometre scale, despite many agricultural and environmental experiments on this scale. A network of 13 tipping‐bucket rain gauges was operated on a 1·4 km² test site in southern Germany for four years to quantify spatial trends in rainfall depth, intensity, erosivity, and predicted runoff. The random measuring error ranged from 10% to 0·1% in case of 1 mm and 100 mm rainfall, respectively. The wind effects could be well described by the mean slope of the horizon at the stations. Except for one station, which was excluded from further analysis, the relative differences due to wind were in maximum ±5%. Gradients in rainfall depth representing the 1‐km² scale derived by linear regressions were much larger and ranged from 1·0 to 15·7 mm km^?1 with a mean of 4·2 mm km^?1 (median 3·3 mm km^?1). They mainly developed during short bursts of rain and thus gradients were even larger for rain intensities and caused a variation in rain erosivity of up to 255% for an individual event. The trends did not have a single primary direction and thus level out on the long term, but for short‐time periods or for single events the assumption of spatially uniform rainfall is invalid on the sub‐kilometre scale. The strength of the spatial trend increased with rain intensity. This has important implications for any hydrological or geomorphologic process sensitive to maximum rain intensities, especially when focusing on large, rare events. These sub‐kilometre scale differences are hence highly relevant for environmental processes acting on short‐time scales like flooding or erosion. They should be considered during establishing, validating and application of any event‐based runoff or erosion model. Copyright © 2009 John Wiley & Sons, Ltd.     

Runoff and sediment yield from land under various uses in a Mediterranean mountain area: long‐term results from an experimental station

In this study we analyzed runoff and sediment yield from land under various traditional and current land uses in Mediterranean mountain areas, using long‐term data from an experimental station in the Aísa Valley, Central Spanish Pyrenees. Monitoring at this station has provided 20 years of data that can help explain the hydrological and geomorphological changes that have been observed at larger spatial scales, and also the changes that have occurred to some of the most characteristic landscapes of the Mediterranean middle mountains. In spite of the problems associated with the use of small experimental plots, the results obtained are consistent with other studies in the Mediterranean region, and confirm the strong influence of land use changes on runoff generation and sediment yield. The results indicate that: (i) cereal cultivation on steep slopes (both alternating cereal cultivation and fallow on sloping fields and shifting agriculture on the steepest slopes) represents a major problem for soil conservation. This explains the occurrence throughout the Mediterranean mountains of many degraded hillslopes, which show evidence of sheet wash erosion, rilling, gullying and shallow landsliding; (ii) farmland abandonment has led to a marked reduction in runoff and sediment yield as a consequence of rapid plant recolonization, particularly by dense shrubs; (iii) the natural transformation of abandoned fields into grazing meadows has reduced runoff and sediment yield. Land use trends in the Mediterranean mountains are mainly characterized by generalized farmland abandonment and a decrease in livestock pressure. From a hydrological and geomorphological point of view the main consequences have been a reduction in overland flow from the hillslopes, and a reduction in sediment sources, with differences up to one order of magnitude in sediment yield from dense shrub cover and grazing meadow areas compared with areas under shifting agriculture. Copyright © 2012 John Wiley & Sons, Ltd.     

Hydrological processes of a closed catchment‐lake system in a semi‐arid Mediterranean environment

Data collected in 4 years of field observations were used in conjunction with continuous simulation models to study, at the small‐basin scale, the water balance of a closed catchment‐lake system in a semi‐arid Mediterranean environment. The open water evaporation was computed with the Penman equation, using the data set collected in the middle of the lake. The surface runoff was partly measured at the main tributary and partly simulated using a distributed, catchment, hydrological model, calibrated with the observed discharge. The simplified structure of the developed modelling mainly concerns soil moisture dynamics and bedrock hydraulics, whereas the flow components are physically based. The calibration produced high efficiency coefficients and showed that surface runoff is greatly affected by soil water percolation into fractured bedrock. The bedrock reduces the storm‐flow peaks and the interflow and has important multi‐year effects on the annual runoff coefficients. The net subsurface outflow from the lake was calculated as the residual of the lake water balance. It was almost constant in the dry seasons and increased in the wet seasons, because of the moistening of the unsaturated soil. During the years of observation, rainfall 30% higher than average caused abundant runoff and a continuous rise in the lake water levels. The analysis allows to predict that, in years with lower than the average rainfall, runoff will be drastically reduced and will not be able to compensate for negative balance between precipitation and lake evaporation. Such highly unsteady situations, with great fluctuations in lake levels, are typical of closed catchment‐lake systems in the semi‐arid Mediterranean environment. Copyright © 2012 John Wiley & Sons, Ltd.     

Quantifying aggregation and change in runoff source in accordance with catchment area increase in a forested headwater catchment

There has been a great deal of research interest regarding changes in flow path/runoff source with increases in catchment area. However, there have been very few quantitative studies taking subscale variability and convergence of flow path/runoff source into account, especially in relation to headwater catchments. This study was performed to elucidate how the contributions and discharge rates of subsurface water (water in the soil layer) and groundwater (water in fractured bedrock) aggregate and change with catchment area increase, and to elucidate whether the spatial variability of the discharge rate of groundwater determines the spatial variability of stream discharge or groundwater contribution. The study area was a 5‐km² forested headwater catchment in Japan. We measured stream discharge at 113 points and water chemistry at 159 points under base flow conditions. End‐member mixing analysis was used to separate stream water into subsurface water and groundwater. The contributions of both subsurface water and groundwater had large variability below 1 km². The contribution of subsurface water decreased markedly, while that of groundwater increased markedly, with increases in catchment area. The specific discharge of subsurface water showed a large degree of variability and decreased with catchment area below 0.1 km², becoming almost constant above 0.1 km². The specific discharge of groundwater showed large variability below 1 km² and increased with catchment area. These results indicated that the variabilities of stream discharge and groundwater contribution corresponded well with the variability of the discharge rate of groundwater. However, below 0.1 km², it was necessary to consider variations in the discharge rates of both subsurface water and groundwater. Copyright © 2016 John Wiley & Sons, Ltd.     

Seasonal and storm flow dynamics of dissolved organic carbon in a Mediterranean mountain catchment (Vallcebre,eastern Pyrenees)

To improve understanding of DOC dynamics in seasonal Mediterranean environments, rainfall, soil water, groundwater and stream water samples were taken during a 27-month period in the Can Vila catchment (northeast Spain). Using these data, we characterized DOC dynamics in the different hydrological compartments and analysed the factors affecting them. We also analysed DOC dynamics during storm events and the factors that control DOC delivery to the stream. The results show some seasonality in rainwater and soil water DOC concentrations, while no clear seasonality was observed in stream water and groundwater, where DOC dynamics were strongly related to discharge and water table variations. For storm events with several discharge peaks, the slope of the discharge–DOC concentration relationship was higher for the first peak. The rather similar dynamics of stream water DOC concentration in all floods contrast with the observed diversity of hydrological processes. This raises the question of the origin of the observed rapid DOC increase.

EDITOR M.C. Acreman

ASSOCIATE EDITOR K. Heal     

Spatio‐temporal variability of the isotopic input signal in a partly forested catchment: Implications for hydrograph separation

The isotopic composition of precipitation (D and ¹⁸O) has been widely used as an input signal in water tracer studies. Whereas much recent effort has been put into developing methodologies to improve our understanding and modelling of hydrological processes (e.g., transit‐time distributions or young water fractions), less attention has been paid to the spatio‐temporal variability of the isotopic composition of precipitation, used as input signal in these studies. Here, we investigated the uncertainty in isotope‐based hydrograph separation due to the spatio‐temporal variability of the isotopic composition of precipitation. The study was carried out in a Mediterranean headwater catchment (0.56 km²). Rainfall and throughfall samples were collected at three locations across this relatively small catchment, and stream water samples were collected at the outlet. Results showed that throughout an event, the spatial variability of the input signal had a higher impact on hydrograph separation results than its temporal variability. However, differences in isotope‐based hydrograph separation determined preevent water due to the spatio‐temporal variability were different between events and ranged between 1 and 14%. Based on catchment‐scale isoscapes, the most representative sampling location could also be identified. This study confirms that even in small headwater catchments, spatio‐temporal variability can be significant. Therefore, it is important to characterize this variability and identify the best sampling strategy to reduce the uncertainty in our understanding of catchment hydrological processes.     

Formulation and application of the unsaturated/saturated catchment models SUSCAT and WEC‐C

This paper describes the formulation and application of a coupled unsaturated/saturated model framework developed to investigate the impact of mining on catchment water yield and groundwater dynamics. The model conceptualization was implemented in both a finite‐element (SUSCAT) and finite‐difference (WEC‐C) solution scheme and found to give similar results. The model framework simulates a coupled surface‐water and groundwater system in which a physically based solution scheme was used to simulate one‐dimensional movement through the unsaturated zone, and a distributed model was used to simulate two‐dimensional saturated groundwater flow. Each soil column comprises a series of layers, each layer being connected to adjacent cells. Subsurface lateral flow is considered when any cell within a layer develops a saturated thickness. Simulation results presented are based on a catchment in the Darling Range, Western Australia that was progressively mined and subsequently rehabilitated. The results predicted the groundwater system beneath the mine areas to have a peak rise owing to mining of between 2 and 4 m. Six years after mining, and following vegetation rehabilitation, the groundwater rise had reduced to 1 m above simulated unmined levels. The corresponding streamflow increase as a result of mining was estimated to peak at 21 mm/year and declined to 7·4 mm/year eight years after revegetation of the mined areas. The simulated groundwater response and streamflow results derived from both models were found to be consistent with observed data. Copyright © 2002 John Wiley & Sons, Ltd.     

Flow behaviour,suspended sediment transport and transmission losses in a small (sub‐bank‐full) flow event in an Australian desert stream

The behaviour of a discrete sub‐bank‐full flow event in a small desert stream in western NSW, Australia, is analysed from direct observation and sediment sampling during the flow event and from later channel surveys. The flow event, the result of an isolated afternoon thunderstorm, had a peak discharge of 9 m³/s at an upstream station. Transmission loss totally consumed the flow over the following 7·6 km. Suspended sediment concentration was highest at the flow front (not the discharge peak) and declined linearly with the log of time since passage of the flow front, regardless of discharge variation. The transmission loss responsible for the waning and eventual cessation of flow occurred at a mean rate of 13.2% per km. This is quite rapid, and is more than twice the corresponding figure for bank‐full flows estimated by Dunkerley (1992) on the same stream system. It is proposed that transmission losses in ephemeral streams of the kind studied may be minimized in flows near bank‐full stage, and be higher in both sub‐bank‐full and overbank flows. Factors contributing to enhanced flow loss in the sub‐bank‐full flow studied included abstractions of flow to pools, scour holes and other low points along the channel, and overflow abstractions into channel filaments that did not rejoin the main flow. On the other hand, losses were curtailed by the shallow depth of banks wetted and by extensive mud drapes that were set down over sand bars and other porous channel materials during the flow. Thus, in contrast with the relatively regular pattern of transmission loss inferred from large floods, losses from low flows exhibit marked spatial variability and depend to a considerable extent on streamwise variations in channel geometry, in addition to the depth and porosity of channel perimeter sediments. Copyright © 1999 John Wiley & Sons, Ltd.     

Suspended sediment transport in a small Mediterranean agricultural catchment

The aim of this study is to analyze suspended sediment transport in a Mediterranean agricultural catchment under traditional soil and water conservation practices. Field measurements were conducted in Can Revull, a small ephemeral catchment (1.03 km²) on the island of Mallorca. This study uses continuous turbidity records to analyse suspended sediment transport regimes, construct and interpret multiple regression models of total suspended sediment concentration (SSC) and of SSC related to stormflow discharge, and assess the sediment loads and yields of three hydrological years (2004–2005 to 2006–2007). An annual average SSC of 17.3 mg l^?1, with a maximum of 2270 mg l^?1, was recorded in the middle of the winter period when rainfall intensities are high and headwater slopes are ploughed and thus bare. Strong seasonal contrasts of baseflow dynamics associated with different degrees of dilution provide a large scatter in SSC and in the derived rating curves, reflecting that other factors control the supply of suspended sediment. Multiple regression models identify rainfall intensity as the most significant variable in sediment supply. However, under baseflow conditions, physical and biological processes generate sediment in the channel that is subsequently removed during high flow. In contrast, when baseflow is not present, rainfall intensity is the only process that supplies sediment to the channel, mostly from hillslopes. Considering the study period as average in terms of total annual rainfall and intensities, suspended sediment yields were an order of magnitude lower than those obtained in other Mediterranean catchments, a factor that can be related to the historical use of soil conservation practices. Copyright © 2009 John Wiley & Sons, Ltd.     

Rainfall,run‐off,and sediment transport dynamics in a humid mountain badland area: Long‐term results from a small catchment

In this study, we investigated rainfall, run‐off, and sediment transport dynamics (414 run‐off events and 231 events with sediment information) of a humid mountain badland area—the Araguás catchment (Central Pyrenees, Spain)—from October 2005 to September 2016. Use of this long‐term database allows characterization of the hydrological response, which consist of low‐magnitude/high‐frequency events and high‐magnitude/low‐frequency events, and identification of seasonal dynamics and rainfall‐run‐off thresholds. Our results indicate that the Araguás catchment, similarly to other humid badlands, had high hydrological responsiveness (mean annual run‐off coefficient: 0.52), a non‐linear relationship of rainfall with run‐off (common in Mediterranean environments), and seasonal hydrological and sedimentological dynamics. We created and validated a multivariate regression model to characterize the hydrological variables (stormflow and peak discharge) and sedimentological variables (mean and maximum suspended sediment concentrations and total suspended sediment load). In summer and at the beginning of autumn, the response was mainly related to rainfall intensity, suggesting a predomination of Hortonian flows. In contrast, in spring and winter, the responses were mainly related to the antecedent conditions (previous rainfall and baseflow), suggesting the occurrence of saturated excess flow processes, and the contribution of neighbouring vegetated areas. The multivariate analysis also showed that total sediment load is better predicted by a multivariate regression model that integrates pre‐event, rainfall, and run‐off variables. In general, our models provided more accurate predictions of small‐magnitude/high‐frequency events than high‐magnitude/low‐frequency events. This study highlights the high inter‐ and intra‐annual variability response in humid badland areas and that long‐term records are needed to reduce the uncertainty of hydrological and sedimentological responses in Mediterranean badland areas.