Defining hydrochemical evolution of streamflow through flowpath dynamics in Kawakami headwater catchment,Central Japan

The hydrochemical behaviour of catchments is often investigated by inferring stream chemistry through identification of source areas involved in hydrograph separation analysis, yet its dynamic evolution of hydrologic pathways has received little attention. Intensive hydrometric and hydrochemical measurements were performed during two different storms on March 29, 2001 and August 21–22, 2001 to define hydrochemical evolution under the dynamic of flow pathways in a 5·2 ha first‐order drainage of the Kawakami experimental basin (KEB), Central Japan, a forested headwater catchment with various soil depths (1·8 to 5 m) overlying late Neogene of volcanic bedrocks. The hydraulic potential distribution and flow lines data showed that the change in flow direction, which was controlled by rainfall amount and antecedent wetness of the soil profile, agreed well with the hydrochemical change across the slope segment during the storm. Hydrograph separation predicted by end‐member mixing analysis (EMMA) using Ca²⁺ and SiO₂ showed that near surface riparian, hillslope soil water and deep riparian groundwater were important in stream flow generation. The evidence of decrease in solutes concentration at a depth of 1 m in the hillslope and 0·6 m in the near surface riparian during peak storm suggested a flushing of high solutes concentration. Most of the solutes accumulated in the deep riparian groundwater zone, which was due to prominent downward flow and agreed well with the residence time. The distinct flow pathways and chemistry between the near surface riparian and deep riparian groundwater zones and the linkage hillslope aquifer and near surface riparian reservoir, which controls rapid flow and solutes flushing during the storm event, are in conflict with the typical assumption that the whole riparian zone resets flow pathways and chemical signature of hillslope soil water, as has been reported in a previous study. Copyright © 2005 John Wiley & Sons, Ltd.     

Flow directions of stream-groundwater exchange in a headwater catchment during the hydrologic year

Understanding near-stream groundwater dynamics and flow directions is important for predicting hillslope-stream connectivity, streamflow generation, and hydrologic controls of streamwater quality. To determine the drivers of groundwater flow in the stream corridor (i.e., the stream channel and the adjacent groundwater in footslopes and riparian areas), we observed the water levels of 36 wells and 7 piezometers along a headwater stream section over a period of 18 months. Groundwater dynamics during events were controlled by the initial position of the groundwater table relative to the subsurface structure. The near-stream groundwater table displayed a fast and pronounced response to precipitation events when lying in fractured bedrock with low storage capacity, and responded less frequently and in a less pronounced way when lying in upper layers with high storage capacity. Precipitation depth, intensity, regolith thickness above the fractured bedrock, and proximity to and elevation above the stream channel also had an effect on the groundwater dynamics, which varied with hydrologic conditions. Our high-frequency and spatially dense measurements highlight the competing influence of groundwater inflow from upslope locations, streamwater level and bedrock properties on the spatiotemporal dynamics of flowpaths in the stream corridor. Near-stream groundwater pointed uniformly towards the stream channel when the stream corridor was hydrologically connected to upslope groundwater. However, local interruptions of the water inflow from upslope locations caused flow reversals towards the footslopes. The direction of near-stream groundwater followed the local fractured bedrock topography during dry hydrologic conditions on a few occasions after events. The outcomes of this research contribute to a better understanding of the drivers controlling spatiotemporal changes in near-stream groundwater dynamics and flow directions in multiple wetness states of the stream corridor.     

Effects of hillslope topography on hydrological responses in a weathered granite mountain,Japan: comparison of the runoff response between the valley‐head and the side slope

To evaluate the effects of hillslope topography on storm runoff in a weathered granite mountain, discharge rate, soil pore water pressures, and water chemistry were observed on two types of hillslope: a valley‐head (a concave hillslope) and a side slope (a planar hillslope). Hydrological responses on the valley‐head and side slope reflected their respective topographic characteristics and varied with the rainfall magnitude. During small rainfall events (<35 mm), runoff from the side slope occurred rapidly relative to the valley‐head. The valley‐head showed little response in storm runoff. As rainfall amounts increased (35–60 mm), the valley‐head yielded a higher flow relative to the side slope. For large rainfall events (>60 mm), runoff from both hillslopes increased with rainfall, although that from the valley‐head was larger than that from the side slope. The differences in the runoff responses were caused by differences in the roles of lower‐slope soils and the convergence of the hillslope. During small rainfall events, the side slope could store little water; in contrast, all rainwater could be stored in the soils at the valley‐head hollow. As the amount of rainfall increased, the subsurface saturated area of the valley‐head extended from the bottom to the upper portion of the slope, with the contributions of transient groundwater via lateral preferential flowpaths due to the high concentration of subsurface water. Conversely, saturated subsurface flow did not contribute to runoff responses, and the subsurface saturated area at the side slope did not extend to the upper slope for the same storm size. During large rainfall events, expansion of the subsurface saturated area was observed in both hillslopes. Thus, differences in the concentration of subsurface water, reflecting hillslope topography, may create differences in the extension of the subsurface saturated area, as well as variability in runoff responses. Copyright © 2007 John Wiley & Sons, Ltd.     

A low‐dimensional physically based model of hydrologic control of shallow landsliding on complex hillslopes

Hillslopes have complex three‐dimensional shapes that are characterized by their plan shape, profile curvature of surface and bedrock, and soil depth. To investigate the stability of complex hillslopes (with different slope curvatures and plan shapes), we combine the hillslope‐storage Boussinesq (HSB) model with the infinite slope stability method. The HSB model is based on the continuity and Darcy equations expressed in terms of storage along the hillslope. Solutions of the HSB equation account explicitly for plan shape by introducing the hillslope width function and for profile curvature through the bedrock slope angle and the hillslope soil depth function. The presented model is composed of three parts: a topography model conceptualizing three‐dimensional soil mantled landscapes, a dynamic hydrology model for shallow subsurface flow and water table depth (HSB model) and an infinite slope stability method based on the Mohr–Coulomb failure law. The resulting hillslope‐storage Boussinesq stability model (HSB‐SM) is able to simulate rain‐induced shallow landsliding on hillslopes with non‐constant bedrock slope and non‐parallel plan shape. We apply the model to nine characteristic hillslope types with three different profile curvatures (concave, straight, convex) and three different plan shapes (convergent, parallel, divergent). In the presented model, the unsaturated storage has been calculated based on the unit head gradient assumption. To relax this assumption and to investigate the effect of neglecting the variations of unsaturated storage on the assessment of slope stability in the transient case, we also combine a coupled model of saturated and unsaturated storage and the infinite slope stability method. The results show that the variations of the unsaturated zone storage do not play a critical role in hillslope stability. Therefore, it can be concluded that the presented dynamic slope stability model (HSB‐SM) can be used safely for slope stability analysis on complex hillslopes. Our results show that after a certain period of rainfall the convergent hillslopes with concave and straight profiles become unstable more quickly than others, whilst divergent convex hillslopes remain stable (even after intense rainfall). In addition, the relation between subsurface flow and hillslope stability has been investigated. Our analyses show that the minimum safety factor (FS) occurs when the rate of subsurface flow is a maximum. In fact, by increasing the subsurface flow, stability decreases for all hillslope shapes. Copyright © 2008 John Wiley & Sons, Ltd.     

A low‐dimensional hillslope‐based catchment model for layered groundwater flow

Despite the strong interaction between surface and subsurface waters, groundwater flow representation is often oversimplified in hydrological models. For instance, the interplay between local or shallow aquifers and deeper regional‐scale aquifers is typically neglected. In this work, a novel hillslope‐based catchment model for the simulation of combined shallow and deep groundwater flow is presented. The model consists of the hillslope‐storage Boussinesq (hsB) model representing shallow groundwater flow and an analytic element (AE) model representing deep regional groundwater flow. The component models are iteratively coupled via a leakage term based on Darcy's law, representing delayed recharge to the regional aquifer through a low conductivity layer. Simulations on synthetic single hillslopes and on a two‐hillslope open‐book catchment are presented, and the results are compared against a benchmark three‐dimensional Richards equation model. The impact of hydraulic conductivity, hillslope plan geometry (uniform, convergent, divergent), and hillslope inclination (0.2%, 5%, and 30%) under drainage and recharge conditions are examined. On the single hillslopes, good matches for heads, hydrographs, and exchange fluxes are generally obtained, with the most significant differences in outflows and heads observed for the 30% slope and for hillslopes with convergent geometry. On the open‐book catchment, cumulative outflows are overestimated by 1–4%. Heads in the confined and unconfined aquifers are adequately reproduced throughout the catchment, whereas exchange fluxes are found to be very sensitive to the hillslope drainable porosity. The new model is highly efficient computationally compared to the benchmark model. The coupled hsB/AE model represents an alternative to commonly used groundwater flow representations in hydrological models, of particular appeal when surface–subsurface exchanges, local aquifer–regional aquifer interactions, and low flows play a key role in a watershed's dynamics. Copyright © 2011 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.     

A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 2. Controls on weathering and permeability profiles,drainage hydraulics,and solute export pathways

The advance of a chemical weathering front into the bedrock of a hillslope is often limited by the rate weathering products that can be carried away, maintaining chemical disequilibrium. If the weathering front is within the saturated zone, groundwater flow downslope may affect the rate of transport and weathering—however, weathering also modifies the rock permeability and the subsurface potential gradient that drives lateral groundwater flow. This feedback may help explain why there tends to be neither “runaway weathering” to great depth nor exposed bedrock covering much of the earth and may provide a mechanism for weathering front advance to keep pace with incision of adjacent streams into bedrock. This is the second of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory. Here, we show how a simplified kinetic model of 1‐D rock weathering can be extended to consider lateral flow in a 2‐D hillslope. Exact and approximate analytical solutions for the location and thickness of weathering within the hillslope are obtained for a number of cases. A location for the weathering front can be found such that lateral flow is able to export weathering products at the rate required to keep pace with stream incision at steady state. Three pathways of solute export are identified: “diffusing up,” where solutes diffuse up and away from the weathering front into the laterally flowing aquifer; “draining down,” where solutes are advected primarily downward into the unweathered bedrock; and “draining along,” where solutes travel laterally within the weathering zone. For each pathway, a different subsurface topography and overall relief of unweathered bedrock within the hillslope is needed to remove solutes at steady state. The relief each pathway requires depends on the rate of stream incision raised to a different power, such that at a given incision rate, one pathway requires minimal relief and, therefore, likely determines the steady‐state hillslope profile.     

Identifying and linking source areas of flow and P transport in dairy‐grazed headwater catchments,North Island,New Zealand

We examined the applicability of the critical‐source area (CSA) concept to the dairy‐grazed 192‐ha Upper Toenepi catchment and its 8·7‐ha Kiwitahi sub‐catchment, New Zealand. We evaluated if phosphorus (P) transport from land into stream is dominated by saturation‐excess (SE) and infiltration‐excess (IE) runoff during stormflow and by sub‐surface (<1·5 m depth) flows during baseflow. We measured stream flow and shallow groundwater levels, collected monthly stream, tile drain (TDA) and groundwater samples, and flow‐proportional stream samples from the Kiwitahi sub‐catchment, and determined their dissolved reactive phosphorus (DRP) and total phosphorus (TP) concentrations. In the Kiwitahi sub‐catchment, during storm events, IE contributions were significant. Contributions from SE appeared significant in the Upper Toenepi catchment. However, in both catchments, sub‐surface contributions dominated stormflow and baseflow periods. Absence of water table at the surface and the water table gradient towards the stream indicated that P transport during events was not limited to surface runoff. The dynamics of the groundwater table and the occurrence of SE areas were influenced by proximity to the stream and hillslope positions. Baseflow accounted for 42% of the annual flow in the Kiwitahi sub‐catchment, and contributed 37 and 52% to the DRP and TP loads, respectively. The P transport during baseflow appeared equally important as P losses from CSAs during stormflow. The close resemblance in P levels between groundwater and stream samples during baseflow demonstrates the importance of shallow groundwater for stream flow. In the Upper Toenepi catchment, contributions from effluent ponds (EFFs) dominated P loads. Management strategies should focus on controlling P release from EFFs, and on decreasing Olsen P concentrations in soil to minimize leaching of P via sub‐surface flow to streams. Research is needed to quantify the role of sub‐surface flow as well as to expand management strategies to minimize P transfers during stormflow and baseflow conditions. Copyright © 2010 John Wiley & Sons, Ltd.     

Contributions of bedrock groundwater to the upscaling of storm‐runoff generation processes in weathered granitic headwater catchments

We examined the contributions of bedrock groundwater to the upscaling of storm‐runoff generation processes in weathered granitic headwater catchments by conducting detailed hydrochemical observations in five catchments that ranged from zero to second order. End‐member mixing analysis (EMMA) was performed to identify the geographical sources of stream water. Throughfall, hillslope groundwater, shallow bedrock groundwater, and deep bedrock groundwater were identified as end members. The contribution of each end member to storm runoff differed among the catchments because of the differing quantities of riparian groundwater, which was recharged by the bedrock groundwater prior to rainfall events. Among the five catchments, the contribution of throughfall was highest during both baseflow and storm flow in a zero‐order catchment with little contribution from the bedrock groundwater to the riparian reservoir. In zero‐order catchments with some contribution from bedrock groundwater, stream water was dominated by shallow bedrock groundwater during baseflow, but it was significantly influenced by hillslope groundwater during storms. In the first‐order catchment, stream water was dominated by shallow bedrock groundwater during storms as well as baseflow periods. In the second‐order catchment, deeper bedrock groundwater than that found in the zero‐order and first‐order catchments contributed to stream water in all periods, except during large storm events. These results suggest that bedrock groundwater influences the upscaling of storm‐runoff generation processes by affecting the linkages of geomorphic units such as hillslopes, riparian zones, and stream channels. Our results highlight the need for a three‐dimensional approach that considers bedrock groundwater flow when studying the upscaling of storm‐runoff generation processes. Copyright © 2014 John Wiley & Sons, Ltd.     

How to model shallow water‐table depth variations: the case of the Kervidy‐Naizin catchment,France

The aim of this work is threefold: (1) to identify the main characteristics of water‐table variations from observations in the Kervidy‐Naizin catchment, a small catchment located in western France; (2) to confront these characteristics with the assumptions of the Topmodel concepts; and (3) to analyse how relaxation of the assumptions could improve the simulation of distributed water‐table depth. A network of piezometers was installed in the Kervidy‐Naizin catchment and the water‐table depth was recorded every 15 min in each piezometer from 1997 to 2000. From these observations, the Kervidy‐Naizin groundwater appears to be characteristic of shallow groundwaters of catchments underlain by crystalline bedrock, in view of the strong relation between water distribution and topography in the bottom land of the hillslopes. However, from midslope to summit, the water table can attain a depth of many metres, it does not parallel the topographic surface and it remains very responsive to rainfall. In particular, hydraulic gradients vary with time and are not equivalent to the soil surface slope. These characteristics call into question some assumptions that are used to model shallow lateral subsurface flow in saturated conditions. We investigate the performance of three models (Topmodel, a kinematic model and a diffusive model) in simulating the hourly distributed water‐table depths along one of the hillslope transects, as well as the hourly stream discharge. For each model, two sets of parameters are identified following a Monte Carlo procedure applied to a simulation period of 2649 h. The performance of each model with each of the two parameter sets is evaluated over a test period of 2158 h. All three models, and hence their underlying assumptions, appear to reproduce adequately the stream discharge variations and water‐table depths in bottom lands at the foot of the hillslope. To simulate the groundwater depth distribution over the whole hillslope, the steady‐state assumption (Topmodel) is quite constraining and leads to unacceptable water‐table depths in midslope and summit areas. Once this assumption is relaxed (kinematic model), the water‐table simulation is improved. A subsequent relaxation of the hydraulic gradient (diffusive model) further improves water‐table simulations in the summit area, while still yielding realistic water‐table depths in the bottom land. Copyright © 2004 John Wiley & Sons, Ltd.     

Experimental evidence of fast groundwater responses in a hillslope/floodplain area in the Black Forest Mountains,Germany

Analysis of water flow pathways from hillslopes to streams is essential for the optimal protection of water resources as well as for ecohydrological studies. This study addresses runoff generation processes at a hillslope and near‐stream shallow groundwater system in the Black Forest Mountains, southwestern Germany. The changing spatial and temporal flow patterns during differing hydrological situations were examined using a combined hydraulic and hydrochemical approach. Groundwater levels at 10 wells, discharge at a near‐stream saturated area, and several natural tracers (deuterium, dissolved silica, and major anions and cations) were observed at different locations during high and low flows. The importance of the groundwater component during flood formation was clearly demonstrated: its contribution was about 80% during a double peak flood event at the saturated area. In addition, a rapid change of the shallow groundwater levels was observed along two transects of groundwater wells in the floodplain. This led to an enhanced groundwater discharge into the saturated area located at the end of one study transect. The amount of groundwater additionally activated during the event was about 30% of total discharge recorded at the outlet of the saturated area. Two alternative hypotheses are discussed to explain this phenomenon: the establishment of locally confined conditions and the development of a pressure wave (hypothesis A), or the significant change of the three‐dimensional groundwater flow lines that caused a large increase of the groundwater catchment at the saturated area during the investigated event (hypothesis B). Even if the exact flow paths and mechanisms could not be clearly identified, the importance of rapid responding hillslope groundwater was undoubtedly demonstrated by a combination of tracer and hydrometric methods. Copyright © 2004 John Wiley & Sons, Ltd.     

Numerical Analysis of Groundwater Ridging Processes Considering Water‐Air Flow in a Hillslope

In this study, a water‐air two‐phase flow model was employed to investigate the formation, extension, and dissipation of groundwater ridging induced by recharge events in a hypothetical hillslope‐riparian zone, considering interactions between the liquid and gas phases in soil voids. The simulation results show that, after a rain begins, the groundwater table near the stream is elevated instantaneously and significantly, thereby generating a pressure gradient driving water toward both the stream (the discharge of groundwater to the stream) and upslope (the extension of groundwater ridging into upslope). Meanwhile, the airflow upslope triggered by the advancing wetting front moves downward gradually. Therefore, the extension of groundwater ridging into upslope and the downward airflow interact within a certain region. After the rain stops, groundwater ridging near the stream declines quickly while the airflow in the lower part of upslope is still moving into the hillslope. Thus, the airflow upslope mitigates the dissipation of groundwater ridging. Additionally, the development of groundwater ridging under different conditions, including rain intensity, intrinsic permeability, capillary fringe height, and initial groundwater table, was analyzed. Changes in intrinsic permeability affect the magnitude of groundwater ridging near the stream, as well as the downward speed of airflow, thereby generating highly complex responses. The capillary fringe is not a controlling factor but an influence factor on the formation of groundwater ridging, which is mainly related to the antecedent moisture. It was demonstrated that groundwater ridging also occurs where an unsaturated zone occurs above the capillary fringe with a subsurface lateral flow.     

Joint interpretation of hydrological and geophysical data: electrical resistivity tomography results from a process hydrological research site in the Black Forest Mountains,Germany

The use of electrical resistivity tomography (ERT; non‐intrusive geophysical technique) was assessed to identify the hydrogeological conditions at a surface water/groundwater test site in the southern Black Forest, Germany. A total of 111 ERT transects were measured, which adopted electrode spacings from 0·5 to 5 m as well as using either Wenner or dipole‐dipole electrode arrays. The resulting two‐dimensional (2D) electrical resistivity distributions are related to the structure and water content of the subsurface. The images were interpreted with respect to previous classical hillslope hydrological investigations within the same research basin using both tracer methods and groundwater level observations. A raster‐grid survey provided a quasi 3D resistivity pattern of the floodplain. Strong structural heterogeneity of the subsurface could be demonstrated, and (non)connectivities between surface and subsurface bodies were mapped. Through the spatial identification of likely flow pathways and source areas of runoff, the deep groundwater within the steeper valley slope seems to be much more connected to runoff generation processes within the valley floodplain than commonly credited in such environmental circumstances. Further, there appears to be no direct link between subsurface water‐bodies adjacent to the stream channel. Deep groundwater sources are also able to contribute towards streamflow from exfiltration at the edge of the floodplain as well as through the saturated areas overlying the floodplain itself. Such exfiltrated water then moves towards the stream as channelized surface flow. These findings support previous tracer investigations which showed that groundwater largely dominates the storm hydrograph of the stream, but the source areas of this component were unclear without geophysical measurements. The work highlighted the importance of using information from previous, complementary hydrochemical and hydrometric research campaigns to better interpret the ERT measurements. On the other hand, the ERT can provide a better spatial understanding of existing hydrochemical and hydrometric data. Copyright © 2009 John Wiley & Sons, Ltd.     

The effect of urbanization on stream hydrology in hillslope watersheds in central Texas

This study examined the effect of urbanization on stream hydrology in hillslope watersheds. Ten streams (seven in hillslope and three in gentle slope watersheds) around Austin, Texas were selected for analysis. For each stream, we compared parameters of transfer function (TF) models estimated from daily rainfall and streamflow data collected in two study periods (October 1988–September 1992 and October 2004–September 2008) representing different degrees of watershed urbanization. As expected, the streams became more intermittent as the watersheds were more urbanized in all the study streams. However, the effect of urbanization on peakflow differs between hillslope and gentle slope watersheds. After watershed urbanization, peakflow increased in gentle slope watersheds, but decreased in hillslope watersheds. Based on the results of the TF models, we found that urbanization made stream not flashier but drier in hillslope watersheds. Overpumpage of aquifer has been recognized as a problem that leads to the stream dryness in the study area. However, the overpumpage alone cannot explain the differences in hydrological changes between the two types of watersheds. We attributed the reduced peakflow and stream dryness in the hillslope watersheds to land grading for construction forming stair‐stepped or terraced landscape. Compared with natural hillslope, a stair‐stepped landscape could infiltrate more stormwater by slowing down surface runoff on tread portions of the stair. Our findings suggest that a watershed management scheme should take into account local hydrogeologic conditions to mitigate the stream dryness resulting from urbanization in hillslope watersheds. Copyright © 2010 John Wiley & Sons, Ltd.     

Spatial and temporal water column and streambed temperature dynamics within an alpine catchment: implications for benthic communities

Despite the significant influence of temperature upon alpine stream benthic communities, thermal regimes of the water column and hyporheic zone of these mountain streams have received limited attention. This paper reports upon a detailed spatio‐temporal study of water column and streambed temperatures undertaken within the Taillon–Gabiétous catchment, French Pyrénées, that aims: (1) to characterize the nature and dynamics of alpine stream water column and streambed thermal patterns; (2) to investigate stream thermal variability under a range of hydroclimatological conditions; and (3) to consider the implications of (1) and (2) for alpine stream benthic communities. The catchment contains four highly dynamic hydrological sources and pathways: (1) two cirque glaciers (Taillon and Gabiétous); (2) seasonal snowpacks; (3) a karst groundwater system; and (4) hillslope aquifers. Water column temperatures were monitored continuously at four sites located along the Taillon glacial stream and at three groundwater springs (two karstic and one hillslope) over the 2002 summer melt season. An eighth site (Tourettes) was established on a predominantly groundwater‐fed stream with limited meltwater input. Bed temperatures (0·05, 0·20 and 0·40 m depth) and river discharge were measured at three sites: (1) the Taillon stream; (2) the Tourettes stream; and (3) below the confluence of (1) and (2). Air temperatures, incoming short‐wave radiation and precipitation were recorded to characterize atmospheric conditions. Glacial stream water column temperatures increased downstream, although groundwater tributaries punctuated longitudinal patterns. Karstic groundwater streams were cooler and more thermally stable than the glacial stream (except at the glacier snout). Hillslope groundwater stream temperatures were most variable and, on average, the warmest of all sites. Streambed temperatures in the glacial stream were coldest and most variable whilst the warmest and least variable streambed temperatures were recorded in an adjacent groundwater tributary. Temperature variability was strongly related to: (1) dynamic water source and pathway contributions; (2) proximity to source; and (3) prevailing hydroclimatological conditions. The high thermal heterogeneity within this catchment may sustain relatively diverse benthic communities, including some endemic Pyrénéan macroinvertebrate taxa. Copyright © 2004 John Wiley & Sons, Ltd.     

Groundwater flushing of solutes at wetland and hillslope positions during storm events in a small glaciated catchment in western New York,USA

While the role of groundwater in flushing of solutes has long been recognized, few studies have explicitly studied the within‐event changes in groundwater chemistry. We compared the changes in groundwater chemistry during storm events for a wetland and hillslope position in a small (1·5 ha) glaciated, forested catchment in western New York. Flushing responses for dissolved organic carbon (DOC) and nitrogen (DON), nitrate (NO₃) and sulfate (SO₄) in wetland and hillslope groundwaters were also compared against the corresponding responses in stream water. Eight storm events with varying intensity, amount, and antecedent moisture conditions were evaluated. Solute flushing patterns for wetland and hillslope groundwaters differed dramatically. While DOC concentrations in wetland groundwater followed a dilution trend, corresponding values for hillslope groundwater showed a slight increase. Concentrations for NO₃ in wetland groundwater were below detection limits, but hillslope groundwaters displayed high NO₃ concentrations with a pronounced increase during storm events. Flushing responses at all positions were also influenced by the size of the event and the time between events. We attributed the differences in flushing to the differences in hydrologic flow paths and biogeochemical conditions. Flushing of the wetland did appear to influence storm‐event stream chemistry but the same could not be said for hillslope groundwaters. This suggests that while a variety of flushing responses may be observed in a catchment, only a subset of these responses affect the discharge chemistry at the catchment outlet. Copyright © 2009 John Wiley & Sons, Ltd.     

Physically based modelling of sheet erosion (detachment and deposition processes) in complex hillslopes

A one‐dimensional uncoupled model governed by this research is a physics‐based modelling of the rainfall‐runoff induced erosion process. The presented model is composed of three parts of a three‐dimensional (3D) hillslope geometry, a nonlinear storage (kinematic wave) model for hillslope hydrological response, and an unsteady physically based surface erosion model. The 3D hillslope geometry model allows describing of the hillslope morphology by defining their plan shape and profile curvature. By changing these two topographic parameters, nine basic hillslope types are derived. The modelling of hillslope hydrological response is based on a flow continuity equation as the relation of discharge and flow depth is passed on kinematic wave approximation. The erosion model is based on a mass conservation equation for unsteady flow. The model assumes that suspended sediment does not affect flow dynamics. The model also accounts for the effect of flow depth plus loose soil depth on soil detachment. The presented model was run for two different precipitations, slope content, and length, and results were plotted for sediment detachment/deposition rate. Based on the obtained results, in hillslopes with convex and straight profile curvatures, sediment detachment only occurred in the whole length of the hillslope. However, in concave ones, sediment detachment and deposition only occurred together in hillslope. The hillslopes with straight profiles and convergent plans have the highest rate of detachment. Also, results show that most detachment rates occur in convex profile curvatures, which are about 15 times more than in straight profiles. Copyright © 2015 John Wiley & Sons, Ltd.     

Bidirectional stream–groundwater flow in response to ephemeral and intermittent streamflow and groundwater seasonality

It is often assumed that the net groundwater flow direction is towards the channel in headwater streams in humid climates, with magnitudes dependent on flow state. However, studies that characterize stream–groundwater interactions in ephemeral and intermittent streams in humid landscapes remain sparse. Here, we examined seasonally driven stream–groundwater interactions in response to temporary streamflow on the basis of field observations of streamflow and groundwater on an adjacent hillslope. The direction of hydraulic head gradients between the stream and groundwater shifted seasonally. The stream gained water (head gradients were towards the stream) when storage state was high. During this period, streamflow was persistent. The stream lost water to the groundwater system (head gradients were away from the stream) when storage state was low. During this period, streamflow only occurred in response to precipitation events, and head gradients remained predominantly away from the stream during events. This suggested that mechanisms other than deep groundwater contributions produced run‐off when storage was low, such as surface and perched subsurface flowpaths above the water table. Analysis of the annual water balance for the study period showed that the residual between precipitation inputs and streamflow and evapotranspiration outputs, which were attributed to the loss of water to the deeper, regional groundwater system, was similar in magnitude to streamflow. This, coupled with results that showed bidirectionality in stream–groundwater head gradients, indicated that headwaters composed of temporary (e.g., ephemeral and intermittent) streams can be important focal areas for regional groundwater recharge, and both contribute to and receive water, solutes, and materials from the groundwater system.     

Dynamics of stormflow generation—A hillslope‐scale field study in east‐central Pennsylvania,USA

A 40 m × 20 m mowed, grass hillslope adjacent to a headwater stream within a 26‐ha watershed in east‐central Pennsylvania, USA, was instrumented to identify and map the extent and dynamics of surface saturation (areas with the water table at the surface) and surface runoff source areas. Rainfall, stream flow and surface runoff from the hillslope were recorded at 5‐min intervals from 11 August to 22 November 1998, and 13 April to 12 November 1999. The dynamics of the water table (0 to 45 cm depth from the soil surface) and the occurrence of surface runoff source areas across the hillslope were recorded using specially designed subsurface saturation and surface runoff sensors, respectively. Detailed data analyses for two rainfall events that occurred in August (57·7 mm in 150 min) and September (83·6 mm in 1265 min) 1999, illustrated the spatial and temporal dynamics of surface saturation and surface runoff source areas. Temporal data analyses showed the necessity to measure the hillslope dynamics at time intervals comparable to that of rainfall measurements. Both infiltration excess surface runoff (runoff caused when rainfall intensity exceeds soil infiltration capacity) and saturation excess surface runoff (runoff caused when soil moisture storage capacity is exceeded) source areas were recorded during these rainfall events. The August rainfall event was primarily an infiltration excess surface runoff event, whereas the September rainfall event produced both infiltration excess and saturation excess surface runoff. Occurrence and disappearance of infiltration excess surface runoff source areas during the rainfall events appeared scattered across the hillslope. Analysis of surface saturation and surface runoff data showed that not all surface saturation areas produced surface runoff that reached the stream. Emergence of subsurface flow to the surface during the post‐rainfall periods appeared to be a major flow process dominating the hillslope after the August rainfall event. Copyright © 2002 John Wiley & Sons, Ltd.