A comparison of fibre‐optic distributed temperature sensing to traditional methods of evaluating groundwater inflow to streams

There are several methods for determining the spatial distribution and magnitude of groundwater inputs to streams. We compared the results of conventional methods [dye dilution gauging, acoustic Doppler velocimeter (ADV) differential gauging, and geochemical end‐member mixing] to distributed temperature sensing (DTS) using a fibre‐optic cable installed along 900 m of Ninemile Creek in Syracuse, New York, USA, during low‐flow conditions (discharge of 1·4 m³ s^?1). With the exception of differential gauging, all methods identified a focused, contaminated groundwater inflow and produced similar groundwater discharge estimates for that point, with a mean of 66·8 l s^?1 between all methods although the precision of these estimates varied. ADV discharge measurement accuracy was reduced by non‐ideal conditions and failed to identify, much less quantify, the modest groundwater input, which was only 5% of total stream flow. These results indicate ambient tracers, such as heat and geochemical mixing, can yield spatially and quantitatively refined estimates of relatively modest groundwater inflow even in large rivers. DTS heat tracing, in particular, provided the finest spatial characterization of groundwater inflow, and may be more universally applicable than geochemical methods, for which a distinct and consistent groundwater end member may be more difficult to identify. Copyright © 2011 John Wiley & Sons, Ltd.     

Transfer of carbon dioxide and methane through the soil‐water‐atmosphere system at Mer Bleue peatland,Canada

Surface waters associated with peatlands, supersaturated with CO₂ and CH₄ with respect to the atmosphere, act as important pathways linking a large and potentially unstable global repository of C to the atmosphere. Understanding the drivers and mechanisms which control C release from peatland systems to the atmosphere will contribute to better management and modelling of terrestrial C pools. We used non‐dispersive infra‐red (NDIR) CO₂ sensors to continuously measure gas concentrations in a beaver pond at Mer Bleue peatland (Canada); measurements were made between July and August 2007. Concentrations of CO₂ in the surface water (10 cm) reached 13 mg C l^?1 (epCO₂ 72), and 26 mg C l^?1 (epCO₂ 133) at depth (60 cm). The study also showed large diurnal fluctuations in dissolved CO₂ which ranged in amplitude from ～1·6 mg C l^?1 at 10 cm to ～0·2 mg C l^?1 at 60 cm depth. CH₄ concentration and supersaturation (epCH₄) measured using headspace analysis averaged 1·47 mg C l^?1 and 3252, respectively; diurnal cycling was also evident in CH₄ concentrations. Mean estimated evasion rates of CO₂ and CH₄ over the summer period were 44·92 ± 7·86 and 0·44 ± 0·25 µg C m^?2s^?1, respectively. Open water at Mer Bleue is a significant summer hotspot for greenhouse gas emissions within the catchment. Our results suggest that CO₂ concentrations during the summer in beaver ponds at Mer Bleue are strongly influenced by biological processes within the water column involving aquatic plants and algae (in situ photosynthesis and respiration). In terms of carbon cycling, soil‐stream connectivity at this time of year is therefore relatively weak. Copyright © 2008 John Wiley & Sons, Ltd.     

Wind effects on sediment transport by raindrop‐impacted shallow flow: a wind‐tunnel study

In wind‐driven rains, wind velocity and direction are expected to affect not only energy input of rains but also shallow ?ow hydraulics by changing roughness induced by raindrop impacts with an angle on ?ow and the unidirectional splashes in the wind direction. A wind‐tunnel study under wind‐driven rains was conducted to determine the effects of horizontal wind velocity and direction on sediment transport by the raindrop‐impacted shallow ?ow. Windless rains and the rains driven by horizontal wind velocities of 6 m s^?1, 10 m s^?1, and 14 m s^?1 were applied to three agricultural soils packed into a 20 by 55 cm soil pan placed on both windward and leeward slopes of 7 per cent, 15 per cent, and 20 per cent. During each rainfall application, sediment and runoff samples were collected at 5‐min intervals at the bottom edge of the soil pan with wide‐mouth bottles and were determined gravimetrically. Based on the interrill erosion mechanics, kinetic energy ?ux (E_rn) as a rainfall parameter and product of unit discharge and slope in the form of q^bS^c_o as a ?ow parameter were used to explain the interactions between impact and ?ow parameters and sediment transport (q_s). The differential sediment transport rates occurred depending on the variation in raindrop trajectory and rain intensity with the wind velocity and direction. Flux of rain energy computed by combining the effects of wind on the velocity, frequency, and angle of raindrop impact reasonably explained the characteristics of wind‐driven rains and acceptably accounted for the differences in sediment delivery rates to the shallow ?ow transport (R₂ ≥ 0·78). Further analysis of the Pearson correlation coef?cients between E_rn and qS_o and q_s also showed that wind velocity and direction signi?cantly affected the hydraulics of the shallow ?ow. E_rn had a smaller correlation coef?cient with the q_s in windward slopes where not only reverse splashes but also reverse lateral raindrop stress with respect to the shallow ?ow direction occurred. However, E_rn was as much effective as qS_o in the sediment transport in the leeward slopes where advance splashes and advance lateral raindrop stress on the ?ow occurred. Copyright © 2004 John Wiley & Sons, Ltd.     

Factors controlling the chemical evolution of travertine‐depositing rivers of the Barkly karst,northern Australia

Groundwaters feeding travertine‐depositing rivers of the northeastern segment of the Barkly karst (NW Queensland, Australia) are of comparable chemical composition, allowing a detailed investigation of how the rate of downstream chemical evolution varies from river to river. The discharge, pH, temperature, conductivity and major‐ion concentrations of five rivers were determined by standard field and laboratory techniques. The results show that each river experiences similar patterns of downstream chemical evolution, with CO₂ outgassing driving the waters to high levels of calcite supersaturation, which in turn leads to widespread calcium carbonate deposition. However, the rate at which the waters evolve, measured as the loss of CaCO₃ per kilometre, varies from river to river, and depends primarily upon discharge at the time of sampling and stream gradient. For example, Louie Creek (Q = 0·11 m³ s^?1) and Carl Creek (Q = 0·50 m³ s^?1) have identical stream gradients, but the loss of CaCO₃ per kilometre for Louie Creek is twice that of Carl Creek. The Gregory River (Q = 3·07 m³ s^?1), O'Shanassy River (Q = 0·57 m³ s^?1) and Lawn Hill Creek (Q = 0·72 m³ s^?1) have very similar gradients, but the rate of hydrochemical evolution of the Gregory River is significantly less than either of the other two systems. The results have major implications for travertine deposition: the stream reach required for waters to evolve to critical levels of calcite supersaturation will, all others things being equal, increase with increasing discharge, and the length of reach over which travertine is deposited will also increase with increasing discharge. This implies that fossil travertine deposits preserved well downstream of modern deposition limits are likely to have been formed under higher discharge regimes. Copyright © 2002 John Wiley & Sons, Ltd.     

Magnetic enhancement in wildfire‐affected soil and its potential for sediment‐source ascription

Intense rainfall following wildfire can cause substantial soil and sediment redistribution. With concern for the increasing magnitude and frequency of wildfire events, research needs to focus on hydrogeomorphological impacts of fire, particularly downstream fluxes of sediment and nutrients. Here, we investigate variation in magnetic enhancement of soil by fire in burnt eucalypt forest slopes to explore its potential as a post‐fire sediment tracer. Low‐frequency magnetic susceptibility values (χ_lf) of <10 µm material sourced from burnt slopes (c. 8·0–10·4 × 10^?6 m³ kg^?1) are an order of magnitude greater than those of <10 µm material derived from long‐unburnt areas (0·8 × 10^?6 m³ kg^?1). Susceptibility of anhysteretic remanent magnetization (χARM) and saturation isothermal remanent magnetization (SIRM) values are similarly enhanced. Signatures are strongly influenced by soil and sediment particle size and storage of previously burnt material in footslope areas. Whilst observations indicate that signatures based on magnetic enhancement show promise for post‐fire sediment tracing, problems arise with the lack of dimensionality in such data. Magnetic grain size indicators χ_fd%, χARM/SIRM and χ_fd/χARM offer further discrimination of source material but cannot be included in numerical unmixing models owing to non‐linear additivity. This leads to complications in quantitatively ascribing downstream sediment to source areas of contrasting burn severity since sources represent numerical multiples of each other, indicating the need to involve additional indicators, such as geochemical evidence, to allow a more robust discrimination. Copyright © 2005 John Wiley & Sons, Ltd.     

Sediment routing hypothesis for pool‐riffle maintenance

This paper provides comprehensive evidence that sediment routing around pools is a key mechanism for pool‐riffle maintenance in sinuous upland gravel‐bed streams. The findings suggest that pools do not require a reversal in energy for them to scour out any accumulated sediments, if little or no sediments are fed into them. A combination of clast tracing using passive integrated transponder (PIT) tagging and bedload traps (positioned along the thalweg on the upstream riffle, pool entrance, pool exit and downstream riffle) are used to provide information on clast pathways and sediment sorting through a single pool‐riffle unit. Computational fluid dynamics (CFD) is also used to explore hydraulic variability and flow pathways. Clast tracing results provide a strong indication that clasts are not fed through pools, rather they are transported across point bar surfaces, or around bar edges (depending upon previous clast position, clast size, and event magnitude). Spatial variations in bedload transport were found throughout the pool‐riffle unit. The pool entrance bedload trap was often found to be empty, when the others had filled, further supporting the notion that little or no sediment was fed into the pool. The pool exit slope trap would occasionally fill with sediment, thought to be sourced from the eroding outer bank. CFD results demonstrate higher pool shear stresses (τ ≈ 140 N m^–2) in a localized zone adjacent to an eroding outer bank, compared to the upstream and downstream riffles (τ ≈ 60 N m^–2) at flows of 6 · 2 m³ s^–1 (≈ 60% of the bankfull discharge) and above. There was marginal evidence for near‐bed velocity reversal. Near‐bed streamlines, produced from velocity vectors indicate that flow paths are diverted over the bar top rather than being fed through the thalweg. Some streamlines appear to brush the outer edge of the pool for the 4 · 9 m³ s^–1 to 7 · 8 m³ s^–1 (between 50 and 80% of the bankfull discharge) simulations, however complete avoidance was found for discharges greater than this. Copyright © 2013 John Wiley & Sons, Ltd.     

Particle travel distances and bed and sediment compositions associated with rain‐impacted flows

Laboratory experiments were performed with rain of uniform drop size (2·7 mm, 5·1 mm) impacting flows over non‐cohesive beds of uniform sized sand (0·11–0·9 mm) and coal (0·2–0·9 mm) particles with flow velocities (20 mm s^?1, 40 mm s^?1) that were insufficient for the flow to entrain the particles without the aid of raindrop impact. Measurement of particle travel distance under rain made up of 2·7 mm drops confirmed a theoretical relationship between settling velocity and the distance particles travel after being disturbed by drop impact. Although, in theory, a relationship between settling velocity and particle travel distance exists, settling velocity by itself was unable to account for the effect of changes in both particle size and density on sediment discharge from beds of uniform non‐cohesive material. Particle density was also a factor. Further study of how particle characteristics influence sediment discharge will aid modelling of the impact of the soil in process‐based models of erosion by rain‐impacted flow. Copyright © 2001 John Wiley & Sons, Ltd.     

Products and Kinetics of the OH‐radical‐induced Dealkylation of Atrazine

Hydroxyl radicals, generated radiolytically in N₂O/O₂‐saturated solutions, yield in their reaction with atrazine equal amounts of deethylatrazine and acetaldehyde (40% of OH radical yield) and deisopropylatrazine and acetone (16%), respectively. The precursors of deethylatrazine and acetaldehyde is their Schiff base which hydrolyzes slowly (OH^–‐catalyzed: k = 5.2 dm³ mol^–1 s^–1). The hydrolysis of the Schiff base of deisopropylatrazine and acetone is too fast to be detected. In a pulse radiolysis experiment, the intermediate formed upon OH‐radical attack (k = 3·10⁹ dm³ mol^–1 s^–1) has a strong absorption at 440 nm. It decays in the presence of oxygen (k = 1.3·10⁹ dm³ mol^–1 s^–1), and upon deprotonation [pK_a(peroxyl radicals) ≈ 10.5] the peroxyl radicals thus‐formed eliminate superoxide radicals (k = 2.9·10⁵ s^–1). s‐Triazine itself reacts much more slowly with OH radicals (k = 9.7·10⁷ dm³ mol^–1 s^–1). This can explain, why in the case of atrazine in comparison to other aromatic compounds, e.g. toluene, the addition of the OH radical to the ring (estimated at ca. 40%) is of relatively little importance as compared to an H‐abstraction from (activated) positions of the side groups.     

Effects of sediment load on hydraulics of overland flow on steep slopes

Eroded sediment may have significant effects on the hydraulics of overland flow, but few studies have been performed to quantify these effects on steep slopes. This study investigated the potential effects of sediment load on Reynolds number, Froude number, flow depth, mean velocity, Darcy–Weisbach friction coefficient, shear stress, stream power, and unit stream power of overland flow in a sand‐glued hydraulic flume under a wide range of hydraulic conditions and sediment loads. Slope gradients were varied from 8·7 to 34·2%, unit flow rates from 0·66 to 5·26×10^?3 m² s^?1, and sediment loads from 0 to 6·95 kg m^?1 s^?1. Both Reynolds number (Re) and Froude number (Fr) decreased as sediment load increased, implying a decrease in flow turbulence. This inverse relationship should be considered in modeling soil erosion processes. Flow depth increased as sediment load increased with a mean value of 1·227 mm, caused by an increase in volume of sediment‐laden flow (contribution 62·4%) and a decrease in mean flow velocity (contribution 37·6%). The mean flow velocity decreased by up to 0·071 m s^?1 as sediment load increased. The Darcy–Weisbach friction coefficient (f) increased with sediment load, showing that the total energy consumption increased with sediment load. The effects of sediment load on f depended on flow discharge: as flow discharge increased, the influence of sediment load on f decreased due to increased flow depth and reduced relative roughness. Flow shear stress and stream power increased with sediment load, on average, by 80·5% and 60·2%, respectively; however, unit stream power decreased by an average of 11·1% as sediment load increased. Further studies are needed to extend and apply the insights obtained under these controlled conditions to real‐world overland flow conditions. Copyright © 2010 John Wiley & Sons, Ltd.     

Towards a conceptual model of hydrological change on an abandoned cutover bog,Quebec

Cutover bogs do not return to functional peatland ecosystems after abandonment because re‐establishment of peat‐forming mosses is poor. This paper presents a conceptual model of bog disturbance caused by peat harvesting (1942–1972), and the hydrological evolution that occurred after abandonment (1973–1998). Two adjacent bogs of similar size and origin, one harvested and the other essentially undisturbed, provide the basis for understanding what changes occurred. The model is based on historical trends evident from previous surveys of land‐use, bog ecology and resource mapping; and from recent hydrological and ecological data that characterize the current condition. Water balance data and historical information suggest that runoff increased and evapotranspiration decreased following drainage, but tended towards pre‐disturbance levels following abandonment, as vegetation recolonized the surface and drainage became less efficient over time. Dewatering of soil pores after drainage caused shrinkage and oxidation of the peat and surface subsidence of approximately 80 cm over 57 years. Comparisons with a nearby natural bog suggest that bulk density in the upper 50 cm of cutover peat increased from 0·07 to 0·13 g cm^?3, specific yield declined from 0·14 to 0·07, water table fluctuations were 67% greater, and mean saturated hydraulic conductivity declined from 4·1 × 10^?5 to 1·3 × 10^?5 cm s^?1. More than 25 years after abandonment, Sphagnum mosses were distributed over broad areas but covered less than 15% of the surface. Areas with ‘good’ Sphagnum regeneration (>10% cover) were strongly correlated with high water tables (mean ?22 cm), especially in zones of seasonal groundwater discharge, artefacts of the extraction history. Forest cover expanded from 5 to 20% of the study area following abandonment. The effect of forest growth (transpiration and interception) and drainage on lowering water levels eventually will be countered by slower water movement through the increasingly dense soil, and by natural ditch deterioration. However, without management intervention, full re‐establishment of natural hydrological functions will take a very long time. Copyright © 2002 John Wiley & Sons, Ltd.     

Solute movement through intact columns of cryoturbated Upper Chalk

Cryoturbated Upper Chalk is a dichotomous porous medium wherein the intra‐fragment porosity provides water storage and the inter‐fragment porosity provides potential pathways for relatively rapid flow near saturation. Chloride tracer movement through 43 cm long and 45 cm diameter undisturbed chalk columns was studied at water application rates of 0·3, 1·0, and 1·5 cm h^?1. Microscale heterogeneity in effluent was recorded using a grid collection system consisting of 98 funnel‐shaped cells each 3·5 cm in diameter. The total porosity of the columns was 0·47 ± 0·02 m³ m^?3, approximately 13% of pores were ≥ 15 µm diameter, and the saturated hydraulic conductivity was 12·66 ± 1·31 m day^?1. Although the column remained unsaturated during the leaching even at all application rates, proportionate flow through macropores increased as the application rate decreased. The number of dry cells (with 0 ml of effluent) increased as application rate decreased. Half of the leachate was collected from 15, 19 and 22 cells at 0·3, 1·0, 1·5 cm h^?1 application rates respectively. Similar breakthrough curves (BTCs) were obtained at all three application rates when plotted as a function of cumulative drainage, but they were distinctly different when plotted as a function of time. The BTCs indicate that the columns have similar drainage requirement irrespective of application rates, as the rise to the maxima (C/C_o) is almost similar. However, the time required to achieve that leaching requirement varies with application rates, and residence time was less in the case of a higher application rate. A two‐region convection–dispersion model was used to describe the BTCs and fitted well (r² = 0·97–0·99). There was a linear relationship between dispersion coefficient and pore water velocity (correlation coefficient r = 0·95). The results demonstrate the microscale heterogeneity of hydrodynamic properties in the Upper Chalk. Copyright © 2007 John Wiley & Sons, Ltd.     

Heat tracing of heterogeneous hyporheic exchange adjacent to in‐stream geomorphic features

There are many field techniques used to quantify rates of hyporheic exchange, which can vary in magnitude and direction spatially over distances of only a few metres, both within and between morphological features. We used in‐stream mini‐piezometers and heat transport modelling of stream and streambed temperatures to quantify the rates and directions of water flux across the streambed interface upstream and downstream of three types of in‐stream geomorphic features: a permanent dam, a beaver dam remnant and a stream meander. We derived hyporheic flux estimates at three different depths at six different sites for a month and then paired those flux rates with measurements of gradient to derive hydraulic conductivity (K) of the streambed sediments. Heat transport modelling provided consistent daily flux estimates that were in agreement directionally with hydraulic gradient measurements and also identified vertical heterogeneities in hydraulic conductivity that led to variable hyporheic exchange. Streambed K varied over an order of magnitude (1·9 × 10^?6 to 5·7 × 10^?5 m/s). Average rates of hyporheic flux ranged from static (q < ±0·02 m/day) to 0·42 m/day. Heat transport modelling results suggest three kinds of flow around the dams and the meander. Copyright © 2010 John Wiley & Sons, Ltd.     

Effects of eastern redcedar encroachment on soil hydraulic properties along Oklahoma's grassland‐forest ecotone

In north‐central Oklahoma eastern redcedar (Juniperus virginiana), encroachment into grassland is widespread and is suspected of reducing streamflow, but the effects of this encroachment on soil hydraulic properties are unknown. This knowledge gap creates uncertainty in understanding the hydrologic effects of eastern redcedar encroachment and obstructs fact‐based management of encroached systems. The objective of this study was to quantify the effects of eastern redcedar encroachment into tallgrass prairie on soil hydraulic properties. Leaf litter depth, soil organic matter, soil water repellency, soil water content, sorptivity, and unsaturated hydraulic conductivity were measured near Stillwater, OK, along 12 radial transects from eastern redcedar trunks to the center of the grassy intercanopy space. Eastern redcedar encroachment in the second half of the 20th century caused the accumulation of 3 cm of hydrophobic leaf litter near the trunks of eastern redcedar trees. This leaf litter was associated with increased soil organic matter in the upper 6 cm of soil under eastern redcedar trees (5.96% by mass) relative to the grass‐dominated intercanopy area (3.99% by mass). Water repellency was more prevalent under eastern redcedar than under grass, and sorptivity under eastern redcedar was 0.10 mm s^?1/2, one seventh the sorptivity under adjacent prairie grasses (0.68 mm s^?1/2). Median unsaturated hydraulic conductivity under grass was 2.52 cm h^?1, four times greater than under eastern redcedar canopies (0.57 cm h^?1). Lower sorptivity and unsaturated hydraulic conductivity would tend to decrease infiltration and increase runoff, but other factors such as rainfall interception by the eastern redcedar canopy and litter layer, and preferential flow induced by hydrophobicity must be examined before the effects of encroachment on streamflow can be predicted. Copyright © 2011 John Wiley & Sons, Ltd.     

Impacts of sediment load on Manning coefficient in supercritical shallow flow on steep slopes

The Manning equation is one of the most widely used formulae for calculating the velocity of shallow overland flow in hydrological and erosion models. Precise estimation of the Manning's friction coefficient (n) is critical to determining overland flow and soil erosion processes. Few studies have been conducted to quantify the effects of sediment load on Manning's n on steep slopes. This study was conducted to investigate the potential effects of sediment load on Manning's n in a flume with a fixed bed, under wide ranges of hydraulics and sediment loads. Slope gradient varied from 8·7 to 34·2%, unit flow rate from 0·66 to 5·26 × 10^?3 m² s^?1, and sediment load from 0 to 6·95 kg m^?1 s^?1. The Reynolds number ranged from 350 to 5899. Results showed that Manning's n varied in both sediment‐free and sediment‐laden flows ranging from 0·012 to 0·055. The apparent Manning's coefficients of sediment‐laden flow were much greater than those of sediment‐free flow. The mean Manning coefficient of sediment‐laden flow was 51·27% greater than the mean value of sediment‐free flow. For sediment‐laden flow, Manning's n could be estimated with a power function of unit flow discharge and sediment content. Further studies are needed to quantify the potential effects of sediment load on the Manning's n on erodible beds and in fields. Copyright © 2010 John Wiley & Sons, Ltd.     

Heuristic numerical and analytical models of the hydrologic controls over vertical solute transport in a domed peat bog,Jura Mountains,Switzerland

We report the results of numerical and analytical simulations to test the hypothesis that downward vertical flow of porewater from the crests of domed alpine and kettle bogs controls vertical porewater distributions of major solutes such as Ca and Mg. The domed Etang de la Gruère bog (EGr), Switzerland, characterized by a vertical downward gradient of 0·04 and stratified layers of peat, is chosen as a field site for the model calibration and evaluation. The middle 4‐m section of the 6·5 m thick bog peat is heavily humified and has a hydraulic conductivity of ～10^?5·6 cm s^?1. Above and below, peat is less humified with a hydraulic conductivity of ～10^?3 cm s^?1. Heuristic finite difference simulations, using Visual MODFLOW, of the bog hydraulics show that the higher conductivity peat at the bog base is critical to create the observed deep, local flow cells that substantively recharge porewater. Model results and Peclet number calculations show that before ～7000 ¹⁴C yr BP diffusion of solutes from underlying mineral soils controlled the vertical distribution of porewater chemistry. From 7000 to ～1250 ¹⁴C BP the porewater chemistry was probably controlled by both upward diffusion and downward advection, and after ～1250 ¹⁴C yr BP porewater chemistry was probably controlled by downward advection. Concentrations of conservative major solutes in the porewaters of alpine, ombrotrophic bogs are the net effect of both downward vertical porewater movement and upward vertical diffusion, the magnitudes of which are delicately poised to the configuration of the bog water table over time and subsurface peat stratigraphy. Copyright © 2002 John Wiley & Sons, Ltd.     

Controls on effective settling velocities of mud flocs in the Changjiang Estuary,China

Theoretical and experimental studies were undertaken to gain insight into physical parameters controlling the flocculation and settling properties of mud flocs in the Changjiang Estuary, China. The Rouse equation is applied to vertical profiles of suspended sediment concentration to determine the bulk mean settling velocity (w_s) of sediment suspended in the Changjiang Estuary. Both in situ point‐sampled and acoustically measured profiles of suspended mud concentrations were fit selectively. The calculated settling velocities w_s mainly ranged from 0·4 to 4·1 mm s^?1 for the point‐sampled data set, and from 1·0 to 3·0 mm s^?1 for the acoustically measured data set. Furthermore, the settling velocities of mud flocs increased with mean concentration (C?) of mud flocs in suspension and were proportional to increasing bottom shear stress (τ_b) of tidal flow. The best equation for the field settling velocity of mud flocs in the Changjiang Estuary can be expressed by the power law: w_s = mC?ⁿ (m, 1·14–2·37; n, 0·84–1·03). It is suggested that C? and τ_b were the dominant physical parameters controlling the flocculation and w_s of mud flocs in suspension. Copyright © 2004 John Wiley & Sons, Ltd.     

Mechanics of interrill erosion with wind‐driven rain

The vector physics of wind‐driven rain (WDR) differs from that of wind‐free rain, and the interrill soil detachment equations in the Water Erosion Prediction Project (WEPP) model were not originally developed to deal with this phenomenon. This article provides an evaluation of the performance of the interrill component of the WEPP model for WDR events. The interrill delivery rates were measured in the wind tunnel facility of the International Center for Eremology (ICE), Ghent University, Belgium with an experimental setup to study different raindrop impact velocity vectors. Synchronized wind and rain simulations with wind velocities of 6, 10 and 14 m s^–1 were applied to a test surface placed on windward and leeward slopes of 7, 15 and 20%. Since both rainfall intensity and raindrop impact velocity varied greatly depending on differences in the horizontal wind velocity under WDRs, the resultant kinetic energy flux (KE_r, in J m^–2 s^–1) was initially used in place of the WEPP model intensity term in order to incorporate the effect of wind on impact velocity and frequency of raindrops. However, our results showed only minor improvement in the model predictions. For all research data, the model Coefficients of Determination (r²) were 0·63 and 0·71, when using the WEPP and the KE_r approaches, respectively. Alternately, integrating the angle of rain incidence into the model by vectorally partitioning normal kinetic energy flux (KE_rn, in J m^–2 s^–1) from the KE_r greatly improved the model's ability to estimate the interrill sediment delivery rates (r² = 0·91). This finding suggested that along with the fall trajectory of wind‐driven raindrops with a given frequency, raindrop velocity and direction at the point of impact onto the soil surface provided sufficient physical information to improve WEPP sediment delivery rate predictions under WDR. Copyright © 2012 John Wiley & Sons, Ltd.     

Persistence of road runoff generation in a logged catchment in Peninsular Malaysia

Measurements of saturated hydraulic conductivity (K_s) and diagnostic model simulations show that all types of logging road/trail in the 14·4 ha Bukit Tarek Experimental Catchment 3 (BTEC3) generate substantial Horton overland flow (HOF) during most storms, regardless of design and level of trafficking. Near‐surface K_s(0–0·05 m) on the main logging road, skid trails and newly constructed logging terraces was less than 1, 2 and 34 mm h^?1, respectively. Near‐surface K_s on an abandoned skid trail in an adjacent basin was higher (62 mm h^?1), owing to the development of a thin organic‐rich layer on the running surface over the past 40 years. Saturated hydraulic conductivity measured at 0·25 m below the surface of all roads was not different (all <6 mm h^?1) and corresponded to the K_s of the adjacent hillslope subsoil, as most roads were excavated into the regolith more than 0·5–1 m. After 40 years, only limited recovery in near‐surface K_s occurred on the abandoned skid trail. This road generated HOF after the storage capacity of the upper near‐surface layer was exceeded during events larger than about 20 mm. Thus, excavation into low‐K_s substrate had a greater influence on the persistence of surface runoff production than did surface compaction by machinery during construction and subsequent use during logging operations. Overland flow on BTEC3 roads was also augmented by the interception of shallow subsurface flow traveling along the soil–saprolite/bedrock interface and return flow emerging from the cutbank through shallow biogenic pipes. The most feasible strategy for reducing long‐term road‐related impacts in BTEC3 is limiting the depth of excavation and designing a more efficient road network, including minimizing the length and connectivity of roads and skid trails. Copyright © 2007 John Wiley & Sons, Ltd.     

Characterization of rainfall generated by dripper‐type rainfall simulator using piezoelectric transducers and its impact on splash soil erosion

Rainfall erosivity represents the primary driver for particle detachment in splash soil erosion. Several raindrop erosivity indices have been developed in order to quantify the potential of rainfall to cause soil erosion. Different types of rainfall simulators have been used to relate rainfall characteristics to soil detachment. However, rainfall produced by different rainfall simulators has different characteristics, specifically different relationships between rainfall intensity and rainfall erosivity. For this reason, the effect of rainfall characteristics produced by a dripper‐type rainfall simulator on splash soil erosion (D_s) has been investigated. The simulated rainfall kinetic energy (KE) and drop size distribution (DSD) were measured using piezoelectric transducers, modified from the Vaisala RAINCAP^® rain sensor. The soil splash was evaluated under various simulated rainfall intensities ranging from 10 to 100 mm h^?1 using the splash‐cup method. The simulated rainfall intensity (I) and kinetic energy relationship (I–KE) was found to be different from natural rainfall. The simulated rainfall intensity and splash soil erosion relationship (I–D_s) also followed this same trend. The I–KE relationship was found to follow the natural rainfall trend until the rainfall intensity reached 30 mm h^?1 and above this limit the KE started to decrease. This emphasizes the importance of the I–KE relationship in determining the I‐D_s relationship, which can differ from one rainfall simulator to another. D_s was found to be highly correlated with KE (r = 0·85, P < 0·001), when data produced by the rainfall intensity ranged from 10 to 100 mm h^?1. However, when the threshold rainfall intensity (30 mm h^?1) was considered, the correlation coefficient further improved (r = 0·89, P = 0·001). Accordingly, to improve the soil splash estimation of simulated rainfall under various rainfall intensities the I–KE characterization relationship for rainfall simulators has to be taken into account. Copyright © 2010 John Wiley & Sons, Ltd.     

Effects of overlying velocity,particle size,and biofilm growth on stream–subsurface exchange of particles

Stream–subsurface exchange strongly influences the transport of contaminants, fine particles, and other ecologically relevant substances in streams. We used a recirculating laboratory flume (220 cm long and 20 cm wide) to study the effects of particle size, overlying velocity, and biofilm formation on stream–subsurface exchange of particles. Sodium chloride was used as a non‐reactive dissolved tracer and 1‐ and 5‐µm fluorescent microspheres were used as particulate tracers. Surface–subsurface exchange was observed with a clean sand bed and a bed colonized by an autotrophic–heterotrophic biofilm under two different overlying velocities, 0·9 and 5 cm s^?1. Hydrodynamic interactions between the overlying flow and sand bed resulted in a reduction of solute and particle concentrations in the water column, and a corresponding accumulation of particles in both the sediments and in the biofilm. Increasing overlying velocity and particle size resulted in faster removal from the overlying water due to enhanced mass transfer to the bed. The presence of the biofilm did not affect solute exchange under any flow condition tested. The presence of the biofilm significantly increased the deposition of particles under an overlying velocity of 5 cm s^?1, and produced a small but statistically insignificant increase at 0·9 cm^?1. The particles preferentially deposited within the biofilm matrix relative to the underlying sand. These results demonstrate that hydrodynamic transport conditions, particle size, and biofilm formation play a key role in the transport of suspended particles, such as inorganic sediments, particulate organic matter, and pathogenic microorganisms in freshwater ecosystems, and should be taken into consideration when predicting the fate and transport of particles and contaminants in the environment. Copyright © 2009 John Wiley & Sons, Ltd.