The nahal yatir bedload database: Sediment dynamics in a gravel-bed ephemeral stream

The Nahal Yatir Bedload Transport Database arises from the Northern Negev River Sediment Monitoring Programme and represents the first body of bedload information to be collected during flash floods in desert gravel-bed streams. Bedload flux was established automatically with three slot-samplers of the Birkbeck type. This was complemented by hydraulic measurements that allow sediment transport to be rated against channel average shear stress and specific stream power. Bed material and bedload grain-size distributions are also provided. The database is given in both printed and electronic tabular format.     

Using a landform evolution model to study ephemeral gullying in agricultural fields: the effects of rainfall patterns on ephemeral gully dynamics

Water driven soil erosion is a major cause of land degradation worldwide. Ephemeral gullies (EGs) are considered key contributors to agricultural catchment soil loss. Despite their importance, the parameters and drivers controlling EG dynamics have not been adequately quantified. Here we investigate the effects of rainfall characteristics on EGs, using the physically based landform evolution model (LEM) CAESAR‐Lisflood. An initial goal of this study was to test the feasibility of using a LEM to estimate EG dynamics based on an easily obtainable and moderate spatial resolution (2 × 2 m) Digital Elevation Model (DEM). EG evolution was simulated for two rainfall seasons in a 0.37 km² agricultural plot situated in a semiarid catchment in central Israel. The 2014 rainfall season was used to calibrate the model and the 2015 season was used for validation. The model overall well predicted the EG network structure and average depth but tended to underestimate the EG length. The effects of rainfall characteristics on EG dynamics were investigated by comparing simulations employing seven rainfall scenarios. Four of these scenarios differ in their overall rainfall volume relative to observed precipitation (+20%, +10%, ?10%, ?20%). The remaining three scenarios vary in the temporal distribution of rainfall during each storm, allowing us to isolate the effect of rainfall intensity on EG evolution. The results show that: (1) EG dynamics strongly correlated with changes in rainfall volume; (2) small‐scale morphological behavior varies between rainfall scenarios, resulting in different meandering and connectivity variability; (3) EG evolution is divided into two main stages, an initial rapid development occurring after the first two weeks of the rainy season, followed by a stable development period; (4) a 12 mm h^?1 intensity threshold was observed to initiate and, later, modify EGs; and (5) inner storm rainfall variability can have a considerable effect on EG evolution. Copyright © 2016 John Wiley & Sons, Ltd.     

Transmission losses to alluvium and associated moisture dynamics in a semiarid ephemeral channel system in Southern Africa

Transmission losses through the bed of ephemeral rivers in arid and semiarid regions can account for a large proportion of the total amount of runoff generated upstream. Losses have typically been estimated by measuring discharge at two points in the channel system. This paper presents some results from a semiarid catchment in South Africa based on moisture observations of the alluvial material using neutron probe access tubes. Channel flow in this catchment is a very infrequent occurrence and consequently presents few opportunities to observe the processes of transmission loss. Only two events have occurred since the beginning of observations in this area (3 October, and 14 November, 1989). the available information is far from complete but suggests that some 75 per cent of the upstream flow in the first event and 22 per cent in the second event was lost to the alluvial material.     

Spatial heterogeneity and temporal dynamics of soil water tension in a mature Norway spruce stand

In ecosystem research great effort is made in measuring soil water tension, because this is a critical calibration variable for modelling soil water fluxes. In this paper the spatial heterogeneity and temporal dynamics of soil tensions and their consequences for the determination of water fluxes are investigated. Studies were carried out at a Norway spruce stand in the Fichtelgebirge (NE Bavaria). Standard tensiometers were installed at three soil depths (20 each) on the whole experimental plot, as well as 45 microtensiometers as a dense grid in a small soil pit. Microtensiometry at the centimetre scale showed that, depending on rain intensity and initial soil water tension, even a soil without discernible macrostructure may show preferential water infiltration. At the stand scale the variability of soil hydraulic properties and tree root distribution causes substantial heterogeneity of soil water tension, as observed by standard tensiometers. A functional relationship between increasing spatial heterogeneity of tensiometer readings and increasing soil water tension was found, which was particularly pronounced after longer dry periods. Also at low soil water tension, where spatial heterogeneity was low, the calculation of water fluxes from tensiometer values was critical, owing to the fact that small differences in measuring soil water tension resulted in big differences in calculated water fluxes. At high soil water tension in summer the spatial heterogeneity of tensiometer readings was extremely high. At our experimental site, since 30% of the total rain in summer falls in events having a precipitation rate greater than 5 mm h⁻¹, preferential water and solute flow was an important phenomenon. We conclude that the validation of calculated water fluxes using measured soil water tension at the stand scale is not an appropriate tool, because of measurement difficulties, considerable spatial heterogeneity, especially in dry periods, and the great variability of soil hydraulic properties. © 1998 John Wiley & Sons, Ltd.     

Aquatic vegetation dynamics in two pit lakes related to interannual water level fluctuation

Water regime characteristics have been recognized as critical factors for aquatic vegetation. In this study, we examined changes in aquatic vegetation coverage area in two shallow sub-lakes of Poyang Lake (Bang Lake and Cuoji Lake) during the dry season from 1987 to 2017. The relationships between eight water regime components (annual average water level, annual maximum water level, annual minimum water level, and flooded days at five water levels [11, 13, 15, 17, and 19 m]) and aquatic vegetation coverage area were determined. The most critical water regimes were identified and results demonstrated that aquatic vegetation coverage area in Bang Lake and Cuoji Lake peaked in drier years (2005 and 2009, respectively) with no obvious up or down trend. Water regimes indicating high flow events such as annual maximum water level, flooded days at water level 19 m, and annual average water level were found to be more important for predicting aquatic vegetation. High-flow events appear to be essential for understanding aquatic vegetation dynamics in pit lakes, yet overall the influences of water level fluctuation on aquatic vegetation varied among wetland units of Poyang Lake. This study helps to understand the hydroecological dynamics in connected lakes further and provide a reference for the lake management and protection.     

Population dynamics of autotrophic picoplankton in relation to environmental factors in a productive lake

Spatial and seasonal fluctuations in autotrophic picoplankton (APP) abundance in a eutrophic, dimictic lake (Lake Aydat, France) were measured concurrently with a variety of environmental variables. Cell number ranged from 0.03 to 2.36×10⁶ cells·ml^–1 (highest concentrations were >5-fold higher than in oligotrophic lakes) and averaged 24 ± 7% of total picoplankton abundance (APP + heterotrophic bacteria). APP abundance (1) peaked in spring simultaneously with heterotrophic flagellate and ciliate densities, (2) decreased during the nitrogen-limited and summer stratification period, and (3) increased with fall turnover. In summer-autumn, the contribution of single-cell eukaryotic (up to 66%) and colonial prokaryotic (18%) forms to total abundance peaked in the bottom waters. Multivariate regression analyses suggest that >40% variance in APP number changes may be explained by ciliate abundance (at 0–4 m depth-range), heterotrophic flagellate number and oxygen concentration (5–9 m), and ciliate carbon biomass (10–14 m). The model accounting for changes in heterotrophic bacterial abundance (5–9 m) indicates chlorophylla concentration (r ²=58%) and ciliate abundance (r ²=34%) as dominant covariates. The data presented here suggest that micrograzers control APP abundance in Lake Aydat.     

Nutrient dynamics in an alpine headwater stream: use of continuous water quality sensors to examine responses to wildfire and precipitation events

Stream water quality can change substantively during diurnal cycles, discrete flow events, and seasonal time scales. In this study, we assessed event responses in surface water nutrient concentrations and biogeochemical parameters through the deployment of continuous water quality sensors from March to October 2011 in the East Fork Jemez River, located in northern New Mexico, USA. Events included two pre‐fire non‐monsoonal precipitation events in April, four post‐fire precipitation events in August and September (associated with monsoonal thunderstorms), and two post‐fire non‐monsoonal precipitation events in October. The six post‐fire events occurred after the Las Conchas wildfire burned a significant portion of the contributing watershed (36%) beginning in June 2011. Surface water nitrate (NO₃? N) concentrations increased by an average of 50% after pre‐fire and post‐fire non‐monsoonal precipitation events and were associated with small increases in turbidity (up to 15 NTU). Beginning 1 month after the start of the large regional wildfire, monsoonal precipitation events resulted in large multi‐day increases in dissolved NO₃? N (6 × background levels), dissolved phosphate (100 × background levels), specific conductance (5 × background levels), and turbidity (>100 × background levels). These periods also corresponded with substantial sags in dissolved oxygen (<4 mg l^?1) and pH (<6.5). The short duration and rapid rates of change during many of these flow events, particularly following wildfire, highlight the importance of continuous water quality monitoring to quantify the timing and magnitude of event responses in streams and to examine large water quality excursions linked to catchment disturbance. Copyright © 2015 John Wiley & Sons, Ltd.     

Rill and ephemeral gully erosion in a small olive grove catchment in Spain: Interactions between management and conservation measures

Concentrated flow erosion in Mediterranean cultivated areas is considered a major process of land degradation. Rills and ephemeral gullies in a 6.4 ha olive orchard catchment located in an intensive commercial farm, which could represent a common degradation scenario in hilly areas, were measured and compared with the sediment loads at the outlet. Four GPS survey campaigns were carried out between March 2009 and March 2014, for periods with cumulative precipitation ranging between 728 and 121 mm. Chopped pruning residues were kept on the lanes during campaigns 2–4, whereas a grass cover crop was seeded in campaign 4. Rainfall, runoff and sediment loads were measured in a flume gauge station at the catchment outlet. The ratio of concentrated flow erosion to catchment sediment load varied between 1.0 and 35.0. Total concentrated flow erosion ranged between 25.0 and 0.1 t ha⁻¹. Rill erosion was the dominant process on ephemeral gullies for three campaigns, with a mean contribution to the total concentrated flow erosion of 55%. Rills clearly followed tractor tracks along the most parallel lanes to the maximum slope. Therefore, a change in traffic direction would be helpful to reduce the connectivity, controlled by the spatial distribution of rows and lanes in the farm. Olive plantation distributions should prioritize hydrological criteria to adapt traffic patterns in olive farms, as long as the tractor manoeuvrability and risk of overturning are not adversely affected. Although the farmer only kept the cover crop for one campaign, it proved to be an efficient measure for interrupting rills along the lanes, whereas pruning residues were effective for the control of interrill erosion. On-site application of pruning residues reduced their handling and transport costs, while increasing the soil fertility and soil cover. However, guidelines for the efficient application of pruning residues are still necessary in terms of residue rates and orientation.     

The influence of climate on water chemistry states and dynamics in rivers across Australia

For effective water quality management and policy development, spatial variability in the mean concentrations and dynamics of riverine water quality needs to be understood. Using water chemistry (calcium, electrical conductivity, nitrate-nitrite, soluble reactive phosphorus, total nitrogen, total phosphorus and total suspended solids) data for up to 578 locations across the Australian continent, we assessed the impact of climate zones (arid, Mediterranean, temperate, subtropical, tropical) on (i) inter-annual mean concentration and (ii) water chemistry dynamics as represented by constituent export regimes (ratio of the coefficients of variation of concentration and discharge) and export patterns (slope of the concentration-discharge relationship). We found that inter-annual mean concentrations vary significantly by climate zones and that spatial variability in water chemistry generally exceeds temporal variability. However, export regimes and patterns are generally consistent across climate zones. This suggests that intrinsic properties of individual constituents rather than catchment properties determine export regimes and patterns. The spatially consistent water chemistry dynamics highlights the potential to predict riverine water quality across the Australian continent, which can support national riverine water quality management and policy development.     

Groundwater dynamics in a till hillslope: flow directions,gradients and delay

Knowledge of groundwater dynamics is important for the understanding of hydrological controls on chemical processes along the water flow pathways. To increase our knowledge of groundwater dynamics in areas with shallow groundwater, the groundwater dynamics along a hillslope were studied in a boreal catchment in Southern Sweden. The forested hillslope had a 1‐ to 2‐m deep layer of sandy till above bedrock. The groundwater flow direction and slope were calculated under the assumption that the flow followed the slope of the groundwater table, which was computed for different triangles, each defined by three groundwater wells. The flow direction showed considerable variations over time, with a maximum variation of 75°. During periods of high groundwater levels the flow was almost perpendicular to the stream, but as the groundwater level fell, the flow direction became gradually more parallel to the stream, directed in the downstream direction. These findings are of importance for the interpretation of results from hillslope transects, where the flow direction usually is assumed to be invariable and always in the direction of the hillslope. The variations in the groundwater flow direction may also cause an apparent dispersion for groundwater‐based transport. In contrast to findings in several other studies, the groundwater level was most responsive to rainfall and snowmelt in the upper part of the hillslope, while the lower parts of the slope reached their highest groundwater level up to 40 h after the upper parts. This can be explained by the topography with a wetter hollow area in the upper part. Copyright © 2011 John Wiley & Sons, Ltd.     

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.     

Life cycle phenology, secondary production, and trophic guilds of caddisfly species in a lake-outlet stream of Patagonia

Lake-outlets are transitional areas recognized as highly productive ecosystems in terms of density and biomass of aquatic insects. Life cycle, secondary production, trophic guilds and environmental constraints of caddisfly assemblages were investigated on a natural lake-outlet stream (Nant y Fall) in Patagonia, Argentina. We investigated the site monthly from May 2007 to April 2008 by recording environmental data and sampling caddis larvae using a Surber sampler at riffle areas (n = 36). Mastigoptila longicornuta, Smicridea annulicornis, Smicridea frequens, Brachysetodes quadrifidus and Parasericostoma ovale displayed well synchronized univoltine life cycles, while Neoatopsyche brevispina, Neopsilochorema tricarinatum, showed an asynchronous development pattern, although most of them had an extended recruitment, similar to those reported for non lake-outlet streams in the area. Annual secondary production per species varied from 8.22 (B. quadrifidus) to 3568.83 mg m⁻² y⁻¹ (P. ovale), with overall caddisfly production amounting to 4.8 g m⁻² y⁻¹. Shredder/collector-filterer ratio was 3.3/1 suggesting that the system was detrital based. Redundancy analysis indicated that seasonally dynamic variables such as discharge, benthic particulate organic matter and temperature were the main predictors of seasonal caddisfly assemblage variation. We propose that the variety of food resources and the significant spatial heterogeneity at lake-outlet streams contribute to sustain a rich caddisfly community.     

Water storage dynamics in a till hillslope: the foundation for modeling flows and turnover times

Studies on hydrology, biogeochemistry, or mineral weathering often rely on assumptions about flow paths, water storage dynamics, and transit times. Testing these assumptions requires detailed hydrometric data that are usually unavailable at the catchment scale. Hillslope studies provide an alternative for obtaining a better understanding, but even on such well‐defined and delimited scales, it is rare to have a comprehensive set of hydrometric observations from the water divide down to the stream that can constrain efforts to quantify water storage, movement, and turnover time. Here, we quantified water storage with daily resolution in a hillslope during the course of almost an entire year using hydrological measurements at the study site and an extended version of the vertical equilibrium model. We used an exponential function to simulate the relationship between hillslope discharge and water table; this was used to derive transmissivity profiles along the hillslope and map mean pore water velocities in the saturated zone. Based on the transmissivity profiles, the soil layer transmitting 99% of lateral flow to the stream had a depth that ranged from 8.9 m at the water divide to under 1 m closer to the stream. During the study period, the total storage of this layer varied from 1189 to 1485 mm, resulting in a turnover time of 2172 days. From the pore water velocities, we mapped the time it would take a water particle situated at any point of the saturated zone anywhere along the hillslope to exit as runoff. Our calculations point to the strengths as well as limitations of simple hydrometric data for inferring hydrological properties and water travel times in the subsurface.     

The dynamics of sediment‐associated contaminants over a transition from drought to multiple flood events in a lowland UK catchment

Fine sediment in suspended form, recently deposited overbank and in temporary storage on or in channel beds, was collected in the Nene basin during a period of drought through to a period of four high flows. The sediment was analysed for arsenic, copper, lead, phosphorus and zinc concentrations with the aim of investigating their sources, movement, temporary storage and potential for environmental harm. Copper, lead and zinc probably originated from urban street dusts, phosphorus (originally in dissolved form) from sewage effluent and arsenic from natural soils developed over ironstone. There was little difference in the metal or arsenic concentrations in the sediment under different flow conditions; instead, proximity to pollutant sources appeared to control their concentrations. Phosphorus in tributary sub‐catchments probably adsorbed to sediment during periods of low flow but these sediments were flushed away during high flows and replaced by sediment with lower concentrations. However, concentrations of all pollutants in overbank sediments along the Nene's main channel were not reduced in successive flood events, suggesting no first flush effect. Only phosphorus accumulated on sediment at concentrations exceeding those of its catchment‐based sources (e.g. street dusts, channel banks and catchment soils). This scavenging of aqueous phosphate by sediment explained the difference in behaviour between phosphorus, arsenic and heavy metals. The surface area and organic matter content were shown to have a small effect on contaminant concentrations. Street dust contaminants only exceeded predicted effect levels in close proximity to urban areas, suggesting a small potential for harm to the aquatic environment. Arsenic concentrations exceeded predicted effect levels in most sediment samples. However, it has been shown to be largely non‐bioavailable in previously published research on the Nene, limiting its potential for significant environmental harm. Phosphorus concentrations in river sediments are high in comparison to the soils in the catchment and in comparison with sediment–P concentrations in other published lowland catchment studies, indicating a large potential for eutrophication should the Phosphorus be, or become, bioavailable. Copyright © 2015 John Wiley & Sons, Ltd.     

Modification of water masses in the Barents Sea and its coupling to ice dynamics: a model study

The physical and biological environment of the Barents Sea is characterised by large variability on a wide range of scales. Results from a numerical ocean model, SINMOD, are presented showing that the physical variability is partly forced by changes in annual net ice import. The mean contribution from ice import in the simulation period (1979–2007) is about 40% of the total amount of ice melted each year. The annual ice import into the Barents Sea varies between 143 and 1,236 km³, and this causes a substantial variability in the amount of annual ice melt in the Barents Sea. This in turn impacts the freshwater content. The simulated freshwater contribution from ice is 0.02 Sv on average and 0.04 Sv at maximum. When mixed into a mean net Atlantic Water (AW) inflow of 1.1 Sv with a salinity of 35.1, this freshwater addition decreases the salinity of the modified AW to 34.4 and 33.9 for the mean and maximum freshwater fluxes, respectively. Ice import may thus be important for the Barents Sea production of Arctic Ocean halocline water which has salinity of about 34.5. The changes in the ice melt the following summer due to ice import also affect the formation of dense water in the Barents Sea by changing stratification, altering the vertical mixing rates and affecting heat loss from the warm AW. The model results thus indicate that ice import from the Arctic has a great impact on water mass modification in the Barents Sea which in turn impacts the ventilation of the Arctic Ocean.     

Estimation of discharges of water flows and debris floods in a small watershed

Floods in small mountainous watersheds cover a wide spectrum of flow. They can range from clear water flows and hyperconcentrated flows to debris floods and debris flows, and calculation of the peak discharge is crucial for predicting and mitigating such hazards. To determine the optimal approach for discharge estimation, this study compared water flow monitoring hydrographs to investigate the performance of five hydrological models that incorporate different runoff yields and influx calculation methods. Two of the models performed well in simulating the peak discharge, peak time, and total flow volume of the water flood. The ratio (γ) of the monitored debris flood discharge (Q_d) to the simulated water flow discharge (Q_w) was investigated. Qualitatively, γ initially increased with Q_w but then decreased when Q_w exceeded a certain threshold, which corresponded to rainfall of 95 and 120 mm in a 6- and 24-h event with a normal distribution of precipitation, respectively. The decrease might be attributable to a threshold of sediment availability being reached, beyond which increased flow rate is not matched by increased sediment input in the large watershed. Uncertainty of hydrological calculation was evaluated by dividing the catchment into sub-basins and adopting different rainfall time steps as input. The efficiency of using a distributed simulation exhibited marginal improvement potential compared with a lumped simulation. Conversely, the rainfall time step input significantly affected the simulation results by delaying the peak time and decreasing the peak discharge. This research demonstrates the applicability of a discharge estimation method that combines a hydrological water flow simulation and an estimation of γ. The results were verified on the basis of monitored flow densities and videos obtained in two watersheds with areas of 2.34 and 32.4 km².     

Lateral dynamics and associated transport of sediment in the upper reaches of a partially mixed estuary,Chesapeake Bay,USA

Data from time series of transects made over a tidal period across a section of the upper Chesapeake Bay, USA, reveal the influence of lateral dynamics on sediment transport in an area with a deep channel and broad extents of shallower flanks. Contributions to lateral momentum by rotation (Coriolis plus channel curvature), cross channel density gradients and cross channel surface slope were estimated, and the friction and acceleration terms needed to complete the balance were compared to patterns of observed lateral circulation. During ebb, net rotation effects were larger because of river velocity and reinforcement of Coriolis by curvature. During flood, stratification was greater because of landward advection of strong vertical density gradients. Together, the ebb intensified lateral circulation and flood intensified stratification focused sediment and sediment transport along the left side of the estuary (looking seaward). The tendency for greater stratification on flood and net sediment flux toward the left-hand shoal are contrary to more common models which, in the northern hemisphere, predict greater resuspension on flood and move sediment toward the right-hand shoal. These tidal asymmetries interact with the lateral circulation to focus net sediment flux on the left side of the estuary, and to produce net ebb directed sediment transport at the surface of the same order of magnitude as net flood directed sediment transport at the bottom.     

Spatial and temporal dynamics of soil moisture after rainfall events along a slope in Regosols of southwest China

High frequency seasonal drought in purple soils (Regosols in FAO taxonomy) of the hilly upland areas of Sichuan basin, China, is one of the key restrictive factors for crop production. In order to manage irrigation and fertilizer application in these soils effectively, the soil water content in a sloped plot with 60 cm soil depth was measured by neutron probe devices to investigate the soil moisture regime during the 1998 rainy season after various amounts of rainfall events. The results showed that variation of soil moisture along the slope positions was highest in the top soil layer during the period of sporadic rainfall that did not induce any runoff. The coefficients of variation of soil moisture at various slope positions (upper, middle, and lower) are 17·36%, 8·95%, 10·25%, 8·58%, 8·05% and 9·21% at the 10 cm, 20 cm, 30 cm, 40 cm, 50 cm and 60 cm soil depths respectively. When surface runoff occurred, the soil moisture dynamics at various positions on the plot were then very different. Soil water content decreased more rapidly on the upper slope than on the middle and lower slope positions. When both surface runoff and throughflow occurred, the soil moisture dynamics in the various layers showed a stable period (soil water content is near constant as time elapses) that lasted about 1 week. Also, the pattern of moisture dynamics is ‘decreasing–stabilization–decreasing’. Thus, irrigation and fertilization management according to the spatial and temporal features of soil moisture dynamics on sloped land can increase the water and fertilizer utilization efficacy by reducing their losses during the stable period. Copyright © 2006 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