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.     

Simulating saturation‐excess surface run‐off in a semi‐distributed hydrological model

The level of complexity, and the number of parameters, to include in a hydrological model is a relatively contentious issue in hydrological modelling. However, it can be argued that explicitly representing important run‐off generation processes can improve the practical value of a model's outputs. This paper explores the benefits of including a new function into an existing semi‐distributed hydrological model (the Pitman model) that is widely used in the sub‐Saharan Africa region. The new function was designed to represent saturation‐excess surface run‐off processes at subcatchment scales and was motivated by the evidence of dambo (low topography riparian areas) type features in many sub‐Saharan river basins. The results for uncertainty versions of the model, with and without the new function, were compared for 25 catchments, which were divided up into those where evidence of dambos exists and those where there is no such evidence. The results suggest that the new function certainly improves the model results for the catchments where dambos exist, but not in situations where saturation‐excess surface run‐off is not expected to occur. The overall conclusion is therefore that the addition of the new function is justified.     

Run‐off processes from mountains to foothills: The role of soil stratigraphy and structure in influencing run‐off characteristics across high to low relief landscapes

The critical zone features that control run‐off generation, specifically at the regional watershed scale, are not well understood. Here, we addressed this knowledge gap by quantitatively and conceptually linking regional watershed‐scale run‐off regimes with critical zone structure and climate gradients across two physiographic provinces in the Southeastern United States. We characterized long‐term (~20 years) discharge and precipitation regimes for 73 watersheds with United States Geological Survey in‐stream gaging stations across the Appalachian Mountain and Piedmont physiographic provinces of North Carolina. Watersheds included in this analysis had <10% developed land and ranged in size from 14.1–4,390 km². Thirty‐four watersheds were located in the Piedmont physiographic province, which is typically classified as a low relief landscape with deep, highly weathered soils and regolith. Thirty‐nine watersheds were located in the Appalachian Mountain physiographic province, which is typically classified as a steeper landscape with highly weathered, but shallower soils and regolith. From the United States Geological Survey daily mean run‐off time series, we calculated annual and seasonal baseflow indices (BFI), minimum, mean, and maximum daily run‐off, and Pearson's correlation coefficients between precipitation and baseflow. Our results showed that Appalachian Mountain watersheds systematically had higher minimum daily flows and BFI values. Piedmont watersheds displayed much larger deviations from mean annual BFI in response to year‐to‐year variability in precipitation. A series of linear regression models between 21 landscape metrics and annual BFIs showed non‐linear and complex terrestrial–hydrological relationships across the two provinces. From these results, we discuss how distinct features of critical zone architecture, with specific focus on soil depth and stratigraphy, may be dominating the regulation of hydrological processes and run‐off regimes across these provinces.     

Identifying run‐off contributions during melt‐induced run‐off events in a glacierized alpine catchment

We analysed contributions to run‐off using hourly stream water samples from seven individual melt‐induced run‐off events (plus one rainfall event) during 2011, 2012 and 2013 in two nested glacierized catchments in the Eastern Italian Alps. Electrical conductivity and stable isotopes of water were used for mixing analysis and two‐component and three‐component hydrograph separation. High‐elevation snowmelt, glacier melt and autumn groundwater were identified as major end‐members. Discharge and tracers in the stream followed the diurnal variations of air temperature but markedly reacted to rainfall inputs. Hysteresis patterns between discharge and electrical conductivity during the melt‐induced run‐off events revealed contrasting loop directions at the two monitored stream sections. Snowmelt contribution to run‐off was highest in June and July (up to 33%), whereas the maximum contribution of glacier melt was reached in August (up to 65%). The maximum groundwater and rainfall contributions were 62% and 11%, respectively. Run‐off events were generally characterized by decreasing snowmelt and increasing glacier melt fractions from the beginning to the end of the summer 2012, while run‐off events in 2013 showed less variable snowmelt and lower glacier melt contributions than in 2012. The results provided essential insights into the complex dynamics of melt‐induced run‐off events and may be of further use in the context of water resource management in alpine catchments. Copyright © 2015 John Wiley & Sons, Ltd.     

Surface moisture area during rainfall–run‐off events to understand the hydrological dynamics of a basin in a plain region

The understanding of the hydrology of plain basins may be improved by the combined analysis of rainfall–run‐off records and remote sensed surface moisture data. Our work evaluates the surface moisture area (SMA) produced during rainfall–run‐off events in a plain watershed of the Argentine Pampas Region, and studies which hydrological variables are related to the generated SMA. The study area is located in the upper and middle basins of the Del Azul stream, characterized by the presence of small gently hilly areas surrounded by flat landscapes. Data from 9 rainfall–run‐off events were analysed. MODIS surface reflectance data were processed to calculate SMA subsequent to the peak discharge (post‐SMA), and previous to the rainfall events (prev‐SMA), to consider the antecedent wetness. Rainfall–run‐off data included total precipitation depth (P), maximum intensity of rainfall over 6 hr (I6max), surface run‐off registered between the beginning of the event and the day previous to the analysed MODIS scene (R), peak flow (Qp), and flood intensity (IF). In contrast with other works, post‐SMA showed a negative relationship with the R. Three groups of cases were identified: (a) Events of low I6max, high prev‐SMA, and low R were associated with slow and weakly channelized flow over plain areas, leading to saturated overland flow (SOF), with large SMA; (b) events of high I6max, low prev‐SMA, and medium to high R were rapidly transported along the gentle slopes of the basin, related to Hortonian overland flow (HOF) and low post‐SMA; and (c) events of medium to high I6max and prev‐SMA with medium R were related to heterogeneous input‐antecedent‐run‐off conditions combined: Local spatial conditions may have produced HOF or SOF, leading to an averaged response with medium SMA. The interactions between the geomorphology of the basin, the characteristics of the events, and the antecedent conditions may explain the obtained results. This analysis is relevant for the general knowledge of the hydrology of large plains, whose functioning studies are still in their early stages.     

Determination of runoff components using path analysis and isotopic measurements in a glacier‐covered alpine catchment (upper Hailuogou Valley) in southwest China

Monitoring of stable water isotopes (δ¹⁸O and δ²H) at the watershed scales can improve our understanding of complex hydrology and hydroclimatology of the watershed, especially in remote regions. Previous studies that used tracers for hydrograph separation are largely based on end‐member mixing approach (EMMA), but one drawback of this approach is that at least two independent tracers are required for multi‐component separation. Here we introduce a new approach—path analysis, in combination with isotopic measurements to investigate the runoff generation in a glacier‐covered alpine catchment (upper Hailuogou Valley) in southwest China. This newly developed method can not only provide a multi‐component hydrograph separation with the aid of only one tracer but also determine the direct and indirect influence of sources on streamflow. Path analysis show that the majority of streamflow is dominated by ice/snow meltwater that represents about 63–78% of the total discharge, whereas precipitation and groundwater contribute approximately 19–39% and 2–4% of the streamflow discharge, respectively. These results are in good agreement with those derived from EMMA (using ¹⁸O and Cl^? as tracers), corroborating that our proposed approach is successful in hydrograph separation of the catchment. This approach may provide new opportunities for the hydrograph separation of catchment with sparse data and be of interest to catchment hydrologists who seek to understand the behaviour of hydrologic systems. Copyright © 2015 John Wiley & Sons, Ltd.     

Linking snowmelt‐derived fluxes and groundwater flow in a high elevation meadow system,Sierra Nevada Mountains,California

Quantifying snowmelt‐derived fluxes at the watershed scale within hillslope environments is critical for investigating local meadow scale groundwater dynamics in high elevation riparian ecosystems. In this article, we investigate the impact of snowmelt‐derived groundwater flux from the surrounding hillslopes on water table dynamics in Tuolumne Meadows, which is located in the Sierra Nevada Mountains of California, USA. Results show water levels within the meadow are controlled by a combination of fluxes at the hillslope boundaries, snowmelt within the meadow and changes in the stream stage. Observed water level fluctuations at the boundaries of the meadow show the hydrologic connection and subsequent disconnection between the hillslope and meadow aquifers. Timing of groundwater flux entering the meadow as a result of spring snowmelt can vary over 20 days based on the location, aspect, and local geology of the contributing area within the larger watershed. Identifying this temporal and spatial variability in flux entering the meadow is critical for simulating changes in water levels within the meadow. Model results can vary significantly based on the temporal and spatial scales at which watershed processes are linked to local processes within the meadow causing errors when boundary fluxes are lumped in time or space. Without a clear understanding of the surrounding hillslope hydrology, it is difficult to simulate groundwater dynamics within high elevation riparian ecosystems with the accuracy necessary for understanding ecosystem response. Copyright © 2010 John Wiley & Sons, Ltd.     

Modelling soil solute release and transport in run‐off on a loessial slope with and without surface stones

Stone covers on loessial slopes can increase the time of infiltration by slowing the velocity of the overland flow, which reduces the transport of solutes, but few mechanistic models have been tested under water‐scouring conditions. We carried out field experiments to test a previously proposed, physically based model of water and solute transport. The area of soil infiltration was calculated from the uncovered surface area, and Richards' equation and the kinematic wave equation were used to describe water infiltration and flow along slopes with stone covers. The transport of chemicals into the run‐off from the surface soil, presumably by diffusion, and their movement in the soil profile could be described by the convection–diffusion equations of the model. The simulated and measured data correlated well. The stones on the soil surface reduced the area available for infiltration but increased the Manning coefficient, eventually leading to increased water infiltration and decreased solute loss with run‐off. Our results indicated that the traditional model of water movement and solute migration could be used to simulate water transport and solute migration for stone‐covered soil on loessial slopes.     

Projected glacier meltwater and river run‐off changes in the Upper Reach of the Shule River Basin,north‐eastern edge of the Tibetan Plateau

Glacier meltwater change in the north‐eastern edge of the Tibetan Plateau is greatly important for the projection of the impact of future climate change on local water resource management. Although the glaciated area is only approximately 4% of the Upper Reach of the Shule River Basin (URSRB), the average glacier meltwater contribution to river run‐off was approximately 23.6% during the periods 1971/1972 to 2012/2013. A new glacier melting module coupled with the macroscale hydrologic Variable Infiltration Capacity model (VIC‐CAS) was adopted to simulate and project changes in the glacier meltwater and river run‐off of the URSRB forced by downscaled output of the BCC‐CSM1.1(m), CANESM2, GFDL‐CM3, and IPSL‐CM5A‐MR models. Comparisons between the observed and simulated river run‐offs and glacier area changes during the periods 2000/2001, 2004/2006, 2008/2009, and 2012/2013 suggest that the simulation is reasonable. Due to increases in precipitation, the annual total run‐off is projected to increase by approximately 2.58–2.73 × 10⁸ m³ in the 2050s and 0.28–1.87 × 10⁸ m³ in the 2100s compared with run‐off in the 2010s based on the RCP2.6 (low greenhouse gas emission) and RCP4.5 (moderate greenhouse gas emission) scenarios, respectively. The contribution of glacier meltwater to river run‐off will more likely decrease to approximately 10% and less than 5% during the 2050s and 2100s, respectively.     

Stream and bed temperature variability in a coastal headwater catchment: influences of surface‐subsurface interactions and partial‐retention forest harvesting

Stream temperature was recorded between 2002 and 2005 at four sites in a coastal headwater catchment in British Columbia, Canada. Shallow groundwater temperatures, along with bed temperature profiles at depths of 1 to 30 cm, were recorded at 10‐min intervals in two hydrologically distinct reaches beginning in 2003 or 2004, depending on the site. The lower reach had smaller discharge contributions via lateral inflow from the hillslopes and fewer areas with upwelling (UW) and/or neutral flow across the stream bed compared to the middle reach. Bed temperatures were greater than those of shallow groundwater during summer, with higher temperatures in areas of downwelling (DW) flow compared to areas of neutral and UW flow. A paired‐catchment analysis revealed that partial‐retention forest harvesting in autumn 2004 resulted in higher daily maximum stream and bed temperatures but smaller changes in daily minima. Changes in daily maximum stream temperature, averaged over July and August of the post‐harvest year, ranged from 1.6 to 3 °C at different locations within the cut block. Post‐harvest changes in bed temperature in the lower reach were smaller than the changes in stream temperature, greater at sites with DW flow, and decreased with depth at both UW and DW sites, dropping to about 1 °C at a depth of 30 cm. In the middle reach, changes in daily maximum bed temperature, averaged over July and August, were generally about 1 °C and did not vary significantly with depth. The pre‐harvest regression models for shallow groundwater were not suitable for applying the paired‐catchment analysis to estimate the effects of harvesting. However, shallow groundwater was warmer at the lower reach following harvesting, despite generally cooler weather compared to the pre‐harvest year. Copyright © 2012 John Wiley & Sons, Ltd.     

Application of the WEPP model to determine sources of run‐off and sediment in a forested watershed

This study investigates critical run‐off and sediment production sources in a forested Kasilian watershed located in northern Iran. The Water Erosion Prediction Project (WEPP) watershed model was set up to simulate the run‐off and sediment yields. WEPP was calibrated and validated against measured rainfall–run‐off–sediment data. Results showed that simulated run‐off and sediment yields of the watershed were in agreement with the measured data for the calibration and validation periods. While low and medium values of run‐off and sediment yields were adequately simulated by the WEPP model, high run‐off and sediment yield values were underestimated. Performance of the model was evaluated as very good and satisfactory during the calibration and validation stages, respectively. Total soil erosion and sediment load of the study watershed during the study period were determined to be 10 108 t yr^?1 and 8735 t yr^?1, respectively. The northern areas of the watershed with dry farming were identified as the critical erosion prone zones. To prioritize the subwatersheds based on their contribution to the run‐off and sediment production at the watershed's main outlet, unit response approach (URA) was applied. In this regard, subwatersheds close to the main outlet were found to have the highest contribution to sediment yield of the whole watershed. Results indicated that depending on the objective of land and water conservation practices, particularly, for controlling sediment yield at the main outlet, critical areas for implementing the best management practices may be identified through conjunctive application of WEPP and URA. Copyright © 2014 John Wiley & Sons, Ltd.     

Changes in run‐off and sediment load in the three parts of the Yellow River basin,in response to climate change and human activities

Hydrological regimes in the Yellow River have changed significantly because of climate change and intensive human interventions. These changes present severe challenges to water resource utilization and ecological development. Variation of run‐off, suspended sediment load (SSL), and eight precipitation indices (P1: 0–12 mm·day^?1, P12: 12–25 mm·day^?1, P25: 25–50 mm·day^?1, P50: P ≥ 50 mm·day^?1 and corresponding rainfall day: Pd1, Pd12, Pd25, Pd50 day year^?1) in three critical parts of the Yellow River basin (source region: SRYRB, upper reaches: URYRB, middle reaches: MRYRB) were investigated for the period from 1960 to 2015. The results show that run‐off and SSL significantly decreased (P < 0.01) in the URYRB and the MRYRB, whereas their decline in the SRYRB was insignificant (P > 0.05). Moreover, run‐off in the URYRB had one change point in 1987, and SSL in the URYRB as well as run‐off and SSL in the MRYRB had two change points (in the 1970s and the 1990s). Over the same period, only Pd1 and Pd12 in the SRYRB showed significant increasing trends, and an abrupt change appeared in 1981. The optimal precipitation indices for assessing the effects of precipitation on run‐off and SSL in the URYRB and MRYRB were Pd50 and P12, respectively. A double‐mass curve analysis showed that precipitation and human activities contributed to approximately 20% and 80% of the reduction in run‐off, respectively, for both the SRYRB and the MRYRB. However, the contribution rate of precipitation and human activities on SSL reduction was approximately 40% and 60% in the URYRB and 5% and 95% in the MRYRB, respectively. Human activities, primarily soil and water conservation measures and water extraction (diversion), were the main factors (>50%) that reduced the run‐off. However, the dominant driving factors for SSL reduction were soil and water conservation measures and reservoir interception, for which the contribution rate was higher than 70% in the MRYRB. This work strengthens the understanding of hydrological responses to precipitation change and provides a useful reference for regional water resource utilization.     

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.     

Testing the influence of topography and material properties on catchment‐scale soil moisture patterns using remotely sensed vegetation patterns in a humid temperate catchment,northern Britain

In order to evaluate the relationship between the apparent complexity of hillslope soil moisture and the emergent patterns of catchment hydrological behaviour and water quality, we need fine‐resolution catchment‐wide data on soil moisture characteristics. This study proposes a methodology whereby vegetation patterns obtained from high‐resolution orthorectified aerial photographs are used as an indicator of soil moisture characteristics. This enables us to examine a set of hypotheses regarding what drives the spatial patterns of soil moisture at the catchment scale (material properties or topography). We find that the pattern of Juncus effusus vegetation is controlled largely by topography and mediated by the catchment's material properties. Characterizing topography using the topographic index adds value to the soil moisture predictions relative to slope or upslope contributing area (UCA). However, these predictions depart from the observed soil moisture patterns at very steep slopes or low UCAs. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimating groundwater recharge for an arid karst system using a combined approach of time‐lapse camera monitoring and water balance modelling

Groundwater is the principal water resource in semi‐arid and arid environments. Therefore, quantitative estimates of its replenishment rate are important for managing groundwater systems. In dry regions, karst outcrops often show enhanced recharge rates compared with other surface and sub‐surface conditions. Areas with exposed karst features like sinkholes or open shafts allow point recharge, even from single rainfall events. Using the example of the As Sulb plateau in Saudi Arabia, this study introduces a cost‐effective and robust method for recharge monitoring and modelling in karst outcrops. The measurement of discharge of a representative small catchment (4.0 · 10⁴ m²) into a sinkhole, and hence the direct recharge into the aquifer, was carried out with a time‐lapse camera. During the monitoring period of two rainy seasons (autumn 2012 to spring 2014), four recharge events were recorded. Afterwards, recharge data as well as proxy data about the drying of the sediment cover are used to set up a conceptual water balance model. The model was run for 17 years (1971 to 1986 and 2012 to 2014). Simulation results show highly variable seasonal recharge–precipitation ratios between 0 and 0.27. In addition to the amount of seasonal precipitation, this ratio is influenced by the interannual distribution of rainfall events. Overall, an average annual groundwater recharge for the doline (sinkhole) catchment on As Sulb plateau of 5.1 mm has estimated for the simulation period. Copyright © 2015 John Wiley & Sons, Ltd.     

Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

Phosphorus and dissolved silicon dynamics in the River Swale catchment,UK: a mass‐balance approach

The internal riverine processes acting upon phosphorus and dissolved silicon were investigated along a 55 km stretch of the River Swale during four monitoring campaigns. Samples of river water were taken at 3 h intervals at sites on the main river and the three major tributaries. Samples were analysed for soluble reactive phosphorus, total dissolved phosphorus, total phosphorus, dissolved silicon and suspended solid concentration. Mass‐balances for each determinand were calculated by comparing the total load entering the river with the total load measured at the downstream site. The difference, i.e. the residual load, showed that there was a large retention of phosphorus and silicon within the system during the March 1998 flood event, but the other three campaigns produced net‐exports. Cumulative residual loads were calculated for each determinand at 6 h intervals throughout each campaign. This incremental approach showed that the mass‐balance residuals followed relatively consistent patterns under various river discharges. During stable low‐flow, there was a retention of particulate phosphorus within the system and also a retention of total dissolved phosphorus and soluble reactive phosphorus, most likely caused by the sorption of soluble phosphorus by bed‐sediments. In times of high river‐discharge, there was a mobilization and export of stored bed‐sediment phosphorus. During overbank flooding, there was a large retention (58% of total input) of particulate phosphorus within the system, due to the mass deposition of phosphorus‐rich sediment onto the floodplain. Soluble phosphorus was also retained within the system by sequestration from the water column by the high concentration of suspended solids. The dissolved silicon mass‐balance residuals had a less consistent pattern in relation to river discharge. There was a large retention of dissolved silicon during overbank flooding, possibly due to sorption onto floodplain soil, and net‐exports during periods of both stable low‐flow and rising limbs of hydrographs, due to release of dissolved silicon from pore‐waters. Copyright © 2001 John Wiley & Sons, Ltd.     

Temporal trends of bulk precipitation and stream water chemistry (1977–1997) in a small forested area,Krusné hory,northern Bohemia,Czech Republic

The Krusné hory (Erzgebirge or Ore Mountains) has been heavily affected by high atmospheric pollutant deposition caused by fossil fuel combustion in an adjacent Tertiary coal basin. Long‐term routine sampling of bulk precipitation (1977–1996) and stream water (1977–1998) in a forested area on the south‐eastern slope of the mountains were used to evaluate trends and patterns in solute concentration and flux with respect to controlling processes. From 1977 to 1996, the annual volume‐weighted Ca²⁺ and SOconcentrations decreased in bulk precipitation. However, after 1989, when a pronounced and continuous decrease occurred in coal production, annual volume‐weighted concentrations decreased for most solutes, except H⁺. The concentration decreases were marked, with 1996 levels at or below 50% of those in 1989. The lack of a trend in H⁺ is attributed to similar decreases in both acid anions and neutralizing base cations. Stream water concentrations of most solutes, i.e. H⁺, Ca²⁺, Mg²⁺, SONOwere highest at the onset of sampling in 1977, decreased markedly from 1977 to 1983 and decreased more gradually from 1983 to 1998. The spruce forest die‐back and removal reduced dry deposition of these solutes by reducing the filtering action, which was provided by the forest canopy. A notable decrease in stream water Ca²⁺ concentrations occurred after 1995 and may be due to the depletion of Ca²⁺, which was provided by catchment liming in 1986, 1988 and 1989. Solute flux trends in bulk atmospheric deposition and stream water generally were not significant and the lack of trend is attributed to the large interannual variability in precipitation quantity and runoff, respectively. All solutes except Na⁺ varied seasonally. The average seasonal concentrations varied between the solutes, but for most solutes were highest in winter and spring and lowest in summer, correlating with the seasonal trend in runoff. For Ca²⁺, Mg²⁺ and SOthe concentration minimum occurs in September and the maximum occurs in February or March, correlating with the seasonal baseflow. These solutes are primarily controlled by the contribution of soil water and groundwater to stream flow. During snowmelt, the meltwater generally causes concentrations to decrease as soil water and groundwater are diluted. For NO₃ , average minimum concentrations occur in August at the end of the growing season concurrent with the lowest stream flow, and the maximum occurs in February and March with high stream flow during snowmelt. Seasonal stream water NOconcentration variations are large compared with the long‐term decrease. Copyright © 1999 John Wiley & Sons, Ltd.     

Detecting the flow relationships between deep and shallow aquifers in an exploited groundwater system,using long‐term monitoring data and quantitative hydrogeology: the Acque Albule basin case (Rome,Italy)

Five years of hydrogeological monitoring and field activities performed in the complex hydrogeological system of the Acque Albule basin (AAB) were conducted to define the hydrogeological setting, the relationship between deep and shallow aquifers and a conceptual groundwater flow model of this exploited area using conventional quantitative techniques. The basin, which is located close to Rome (Italy) on the west side of the Apennine chain and just north of the Colli Albani volcano, subsided after development of a north–south fault system (about 115 000 y bp). The AAB experiences intense hydrothermal activity, which has produced a large travertine deposit (80‐m thick). The travertine deposit constitutes a fractured aquifer that is the final destination of more than 5 m³ s^‐1 of water and is strongly dewatered by quarry activities. The complex hydrogeology of this basin was investigated, revealing two main hydraulically connected aquifers, one thermalised and partly confined into the limestone bedrock and one unconfined in the travertine. The two aquifers are separated by a non‐continuous clayey aquiclude. The hydrogeological survey and geological characterisation contributed to the development of the groundwater flow conceptual model. Analysis and comparison of the monitored levels highlighted the pattern of flow between the deep and shallow parts of the flow system. Copyright © 2012 John Wiley & Sons, Ltd.     

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.