Accuracy of Flowmeters Measuring Horizontal Groundwater Flow in an Unconsolidated Aquifer Simulator

Borehole flowmeters that measure horizontal flow velocity and direction of groundwater flow are being increasingly applied to a wide variety of environmental problems. This study was carried out to evaluate the measurement accuracy of several types of flowmeters in an unconsolidated aquifer simulator. Flowmeter response to hydraulic gradient, aquifer properties, and well‐screen construction was measured during 2003 and 2005 at the U.S. Geological Survey Hydrologic Instrumentation Facility in Bay St. Louis, Mississippi. The flowmeters tested included a commercially available heat‐pulse flowmeter, an acoustic Doppler flowmeter, a scanning colloidal borescope flowmeter, and a fluid‐conductivity logging system. Results of the study indicated that at least one flowmeter was capable of measuring borehole flow velocity and direction in most simulated conditions. The mean error in direction measurements ranged from 15.1° to 23.5° and the directional accuracy of all tested flowmeters improved with increasing hydraulic gradient. The range of Darcy velocities examined in this study ranged 4.3 to 155 ft/d. For many plots comparing the simulated and measured Darcy velocity, the squared correlation coefficient (r²) exceeded 0.92. The accuracy of velocity measurements varied with well construction and velocity magnitude. The use of horizontal flowmeters in environmental studies appears promising but applications may require more than one type of flowmeter to span the range of conditions encountered in the field. Interpreting flowmeter data from field settings may be complicated by geologic heterogeneity, preferential flow, vertical flow, constricted screen openings, and nonoptimal screen orientation.     

The Electromagnetic Borehole Flowmeter: Description and Application

Borehole flowmeters are downhole tools that measure axial flow in a well or borehole. Desirable flowmeter characteristics include low detection limit, a wide range of operation, accuracy, durability, reliable performance, and a small diameter and length. The recently developed electromagnetic (F.M) flowmeter has these trails. The first portion of this paper presents the MM flowmeter design, provides laboratory calibration data, and compares the performance characteristics of MM flowmeters 10 those of impeller and thermal pulse flowmeters. The second portion of the paper discusses applications of the MM flowmeter.     

Evaluating the Use of In‐Well Heat Tracer Tests to Measure Borehole Flow Rates

Recent research has demonstrated the use of in‐well heat tracer tests monitored by a fiber optic distributed temperature sensing (DTS) system to characterize borehole flow conditions in open bedrock boreholes. However, the accuracy of borehole flow rates determined from in‐well heat tracer tests has not been evaluated. The purpose of the research presented here is to determine whether borehole flow rates obtained using DTS‐monitored in‐well heat tracer tests are reasonable, and to evaluate the range of flow rates measureable with this method. To accomplish this, borehole flow rates measured using in‐well heat tracer tests are compared to borehole flow rates measured in the same boreholes using an impeller or heat pulse flowmeter. A comparison of flow rates measured using in‐well heat tracer tests to flow rates measured with an impeller flowmeter under the same conditions showed good agreement. A comparison of in‐well heat tracer test flow rate measurements to previously‐collected heat pulse flowmeter measurements indicates that the heat tracer test results produced borehole flow rates and flow profiles similar to those measured with the heat pulse flowmeter. The results of this study indicate that borehole flow rates determined from DTS‐monitored in‐well heat tracer tests are reasonable estimates of actual borehole flow rates. In addition, the range of borehole flow rates measurable by in‐well heat tracer tests spans from less than 10^?1 m/min to approximately 10¹ m/min, overlapping the ranges typically measurable with an impeller flowmeter or a heat pulse flowmeter, making in‐well heat tracer testing a versatile borehole flow logging tool.     

Accuracy of Flowmeters Measuring Horizontal Flow in Fractured-Rock Simulators

Laboratory evaluations of flowmeter response to flow in fractured-rock simulators are needed to improve understanding of data collected in field settings. The ability of flowmeters to accurately measure the velocity and direction of water flowing between parallel plates was used as a surrogate for instrument response in fractured-rock aquifers. A colloidal borescope flowmeter and a heat-pulse flowmeter were deployed in a fractured rock simulator with 4-inch and 6-inch inner-diameter, uncased wells with 0.39- and 1.0-inch fracture apertures and groundwater velocities from 35 to 975 ft/d. The colloidal borescope measurements and applied velocities were positively correlated in all wells and apertures (the coefficient of determination [r²] = 0.61–0.89) and most accurately measured direction at higher velocities. The mean directional error in colloidal borescope measurements was less than 17° in 6-inch wells and 31° in the 4-inch wells at velocities between 92 and 958 ft/d. Heat-pulse flowmeter measurements were 0.001 to 0.004 times less than applied rates and may indicate that water was moving around rather than through the instrument's integrated packer. The mean directional error of heat-pulse flowmeter measurements were about 18 and 42° in the 0.39- and 1.0-inch fractures, respectively, for groundwater velocities within the manufacturer's suggested range of application (0.5–100 ft/d). Measurements made at vertical increments and fracture positions in the well using the colloidal borescope indicate that laminar flow occurs within the central 50% of the fracture but measurements above or below are likely affected by eddy currents.     

Characterization of a multilayer aquifer using open well dilution tests

An approach to characterization of multilayer aquifer systems using open well borehole dilution is described. The approach involves measuring observation well flow velocities while a nearby extraction well is pumped by introducing a saline tracer into observation wells and collecting dilution vs. depth profiles. Inspection of tracer profile evolution allows discrete permeable layers within the aquifer to be identified. Dilution profiles for well sections between permeable layers are then converted into vertical borehole flow velocities and their evolution, using an analytic solution to the advection-dispersion equation applied to borehole flow. The dilution approach is potentially able to measure much smaller flow velocities that would be detectable using flowmeters. Vertical flow velocity data from the observation wells are then matched to those generated using a hydraulic model of the aquifer system, "shorted" by the observation wells, to yield the hydraulic properties of the constituent layers. Observation well flow monitoring of pumping tests represents a cost-effective alternative or preliminary approach to pump testing each layer of a multilayer aquifer system separately using straddle packers or screened wells and requires no prior knowledge of permeable layer depths and thicknesses. The modification described here, of using tracer dilution rather than flowmeter logging to obtain well flow velocities, allows the approach to be extended to greater well separations, thus characterizing a larger volume of the aquifer. An example of the application of this approach to a multilayer Chalk Aquifer in Yorkshire, Northeast England, is presented.     

River-induced flow dynamics in long-screen wells and impact on aqueous samples

Previously published field investigations and modeling studies have demonstrated the potential for sample bias associated with vertical wellbore flow in conventional monitoring wells constructed with long-screened intervals. This article builds on the existing body of literature by (1) demonstrating the utility of continuous (i.e., hourly measurements for ～1 month) ambient wellbore flow monitoring and (2) presenting results from a field experiment where relatively large wellbore flows (up to 4 L/min) were induced by aquifer hydrodynamics associated with a fluctuating river boundary located approximately 250 m from the test well. The observed vertical wellbore flows were strongly correlated with fluctuations in river stage, alternating between upward and downward flow throughout the monitoring period in response to changes in river stage. Continuous monitoring of ambient wellbore flows using an electromagnetic borehole flowmeter allowed these effects to be evaluated in concert with continuously monitored river-stage elevations (hourly) and aqueous uranium concentrations (daily) in a long-screen well and an adjacent multilevel well cluster. This study demonstrates that when contaminant concentrations within the aquifer vary significantly over the depth interval interrogated, river-induced vertical wellbore flow can result in variations in measured concentration that nearly encompass the full range of variation in aquifer contaminant concentration with depth.     

Estimation of water velocities in boreholes using a new passive flowmeter

Abstract

Preferential flow pathways in a fractured aquifer may yield abrupt reductions of the water velocity in a well. We propose a new device for measuring low (5–13 cm d^-1) velocities in wells originating from fractures at different depths. The presented flowmeter has been applied in a well in the Bari (southern Italy) fractured aquifer. In the same well, the horizontal flowmeter velocity (9.6 cm d^-1) at 0.5 m depth was compared with velocity (8 cm d^-1) derived from a field tracer test, providing a value 16.5% higher. Moreover, the flowmeter measurements at 1.5 m depth gave a horizontal velocity of 7.2 cm d^-1, which is 11% less than water flow velocity estimated from the field test. The new flowmeter implements the tracer point-dilution method in a plastic (PVC) pipe by causing the water flow to pass through an artificial filter. Laboratory calibration tests have confirmed the good performance of the proposed flowmeter technique, even for water flow up to 300 cm d^-1. The flowmeter was sensitive to 0.1 cm d^-1, with a detection limit of 1.5 cm d^-1, i.e. half the measurable flow velocity of existing flowmeters in wells.

Editor D. Koutsoyiannis; Associate editor S. Grimaldi     

Increased Thermal-Pulse Flowmeter Resolution by Adding a Packer and Computer

Measurement accuracy was increased by nearly one order of magnitude by outfitting the thermal-pulse flowmeter (TFM) with an inflatable packer. To accurately measure slow water velocities in boreholes greater than 15 cm diameter, it is necessary to divert borehole fluids through the TFM by inflating a packer. During calibration it was noted that the TFM's accuracy decreased as the borehole diameter increased. With Lhe packer inflated the TFM has a useful flow measurement range of 0.08 to 15 L/min (with flow velocities of 0.24 ± 0.012 cm/inin to 45.7 ± 0.61 cm/min, respectively, in 20-cm-diameter pipe), compared to 0.8 to 57 L/min for a packcrless TFM. A computer interlace was added to the TFM to provide a real-time graphical display of the differential voltage output from the TFM, a running mean and standard deviation of the pulse-response time, and a mean flow rate and velocity based on calibration curve fits.     

Investigation into the application of cross-correlation analysis on acoustic backscattered signals from suspended sediments to measure nearbed current profiles

In the present work, the potential of correlation methods applied to acoustic signals backscattered from suspended sediments, is examined with regard to the measurement of near-bed vertical profiles of the horizontal current velocity. The technique uses the cross-correlation of the backscattered signal between pairs of downward looking horizontally separated transducers, to measure the time taken for suspension structures to advect from one transducer to another. This approach requires a degree of spatial coherence in the suspension field, and advantage is taken of this coherence to measure the flow. To investigate the technique, backscattered data collected using an array of transducers has been examined. The array was deployed along stream in an estuarine environment, which was subject to strong rectilinear currents, and high levels of suspended sediment concentration. Using pairs of transducers with separations between 0.25–5.42 m, the structure of the suspension field has been assessed, the implication of these observations for measuring flow discussed, and nearbed measurements of current profiles obtained. The results show that in the estuarine environment investigated, if the transducer pair separation was of the order of a metre or less, and currents were sufficient to retain sediments in suspension, vertical profiles of the horizontal current could be measured. To assess the results conventional electromagnetic and impeller current meter measurements were collected at four heights above the bed. Comparison of the conventional reference measurements with the cross-correlation velocities show that the correlation method has the potential to provide reliable measurements of near-bed current profiles.     

Vertical Wellbore Flow Monitoring for Assessing Spatial and Temporal Flow Relationships with a Dynamic River Boundary

A useful tool for identifying the temporal and spatial ambient wellbore flow relationships near a dynamic river boundary is to monitor ambient vertical wellbore flow with an electromagnetic borehole flowmeter. This is important because the presence of the wellbore can result in significant mixing or exchange of groundwater vertically across the aquifer. Mixing or exchanging groundwater within the well-screen section can have significant impacts on the distribution of contaminants within the aquifer and adverse effects on the representativeness of groundwater samples collected from the monitoring well. Ambient monitoring data, collected from long screened wells at Hanford’s 300-Area Integrated Field Research Challenge site, located approximately 260 m from the Columbia River, demonstrate that vertical wellbore flow exhibits both a positive and inverse temporal relationship with periodic river-stage fluctuations that can change over short distances between wells. The spatial distribution of these vertical flows across the well field indicates two general regions of ambient wellbore flow behavior. The western region of the site is characterized by vertical flows that are positively related to river-stage fluctuations. In contrast, the eastern region of the site exhibits vertical flows that are inversely related to river-stage fluctuations. The cause of this opposite relationship is not completely understood; however, the positive relationships appear to be associated with high-energy Hanford formation flood deposits. These flood deposits have a well-defined northwest-southeast trend and are believed to coincide with a local paleochannel. The inverse relationships are attributed to an erosional, subsurface high in the Hanford/Ringold Formation contact between the site and the Columbia River. Under these complex hydrogeologic and hydrodynamic conditions, the behavior of ambient vertical wellbore flow in monitoring wells near a dynamic river boundary can have important implications for collecting groundwater-quality samples, for contributing to contaminant distribution within an aquifer system, and for implementing effective remediation strategies.     

The relationship between VHF radar auroral backscatter amplitude and Doppler velocity: a statistical study

A statistical investigation of the relationship between VHF radar auroral backscatter intensity and Doppler velocity has been undertaken with data collected from 8 years operation of the Wick site of the Sweden And Britain Radar-auroral Experiment (SABRE). The results indicate three different regimes within the statistical data set; firstly, for Doppler velocities <200 m s⁻¹, the backscatter intensity (measured in decibels) remains relatively constant. Secondly, a linear relationship is observed between the backscatter intensity (in decibels) and Doppler velocity for velocities between 200 m s⁻¹ and 700 m s⁻¹. At velocities greater than 700 m s⁻¹ the backscatter intensity saturates at a maximum value as the Doppler velocity increases. There are three possible geophysical mechanisms for the saturation in the backscatter intensity at high phase speeds: a saturation in the irregularity turbulence level, a maximisation of the scattering volume, and a modification of the local ambient electron density. There is also a difference in the dependence of the backscatter intensity on Doppler velocity for the flow towards and away from the radar. The results for flow towards the radar exhibit a consistent relationship between backscatter intensity and measured velocities throughout the solar cycle. For flow away from the radar, however, the relationship between backscatter intensity and Doppler velocity varies during the solar cycle. The geometry of the SABRE system ensures that flow towards the radar is predominantly associated with the eastward electrojet, and flow away is associated with the westward electrojet. The difference in the backscatter intensity variation as a function of Doppler velocity is attributed to asymmetries between the eastward and westward electrojets and the geophysical parameters controlling the backscatter amplitude.     

Tidal effects on variations of fresh–saltwater interface and groundwater flow in a multilayered coastal aquifer on a volcanic island (Jeju Island, Korea)

We conducted various field studies at the seawater intrusion monitoring wells located in the eastern part of Jeju Island, Korea, to observe the tidal effect on groundwater–seawater flow in the coastal aquifer. Studies included monitoring the fluctuations of groundwater and tide levels, electrical and temperature logging, and 2-D heat-pulse flowmeter tests. According to time-series analysis, tidal effects on groundwater level reached up to 3 km inland from the coastline. Water-level variation was more sensitive to tidal fluctuations near the coast, and more related to rainfall toward inland areas. Temporal and spatial variations in the shape and location of the freshwater–saltwater interface were analyzed using data from nine monitoring wells. The results indicated that the interface toe is located at a distance of 6–8 km from the coastline and its location was related to geological layers present. Long-term seasonal variations revealed no major changes in the interface; minor variations were due to moving boundary conditions induced by tidal fluctuations. Using the two-dimensional heat-pulse flowmeter, groundwater flow directions and velocities at four tidal stages were measured on three monitoring wells drilled into the multilayered aquifers. This direct measurement enabled us to relate the differences of flow velocities and directions with geology and tidal fluctuations. Combining the results of EC logging and flowmeter tests, we found a zone where freshwater and saltwater moved alternately in opposite directions, as influenced by the tidal fluctuations. Integrating various physical logging and flowmeter data with water-level fluctuations improved our understanding of the behavior of fresh and seawater flow in the coastal aquifers.     

The Influence of Wellbore Inflow on Electromagnetic Borehole Flowmeter Measurements

This paper describes a combined field, laboratory, and numerical study of electromagnetic borehole flowmeter measurements acquired without the use of a packer or skirt to block bypass flow around the flowmeter. The most significant finding is that inflow through the wellbore screen changes the ratio of flow through the flowmeter to wellbore flow. Experiments reveal up to a factor of two differences in this ratio for conditions with and without inflow through the wellbore screen. Standard practice is to assume the ratio is constant. A numerical model has been developed to simulate the effect of inflow on the flowmeter. The model is formulated using momentum conservation within the borehole and around the flowmeter. The model is embedded in the MODFLOW-2000 ground water flow code.     

Temporal changes in the vertical distribution of flow and chloride in deep wells

The combination of flowmeter and depth-dependent water-quality data was used to evaluate the quantity and source of high-chloride water yielded from different depths to eight production wells in the Pleasant Valley area of southern California. The wells were screened from 117 to 437 m below land surface, and in most cases, flow from the aquifer into the wells was not uniformly distributed throughout the well screen. Wells having as little as 6 m of screen in the overlying upper aquifer system yielded as much as 50% of their water from the upper system during drought periods, while the deeper parts of the well screens yielded 15% or less of the total yield of the wells. Mixing of water within wells during pumping degraded higher-quality water with poorer-quality water from deeper depths, and in some cases with poorer-quality water from the overlying upper aquifer system. Changes in the mixture of water within a well, resulting from changes in the distribution of flow into the well, changed the quality of water from the surface discharge of wells over time. The combination of flowmeter and depth-dependent water quality data yielded information about sources of high-chloride water to wells that was not available on the basis of samples collected from nearby observation wells. Changing well design to eliminate small quantities of poor-quality water from deeper parts of the well may improve the quality of water from some wells without greatly reducing well yield.     

An Evaluation of Interpolation Methodologies for Generating Three-Dimensional Hydraulic Property Distributions from Measured Data

The process of attempting to model ground-water systems requires a good understanding of the spatial variation of aquifer hydraulic properties. The capabilities of the more recent innovative flowmeters such as the electromagnetic and heat pulse flowmeters provide the sensitivity to measure ambient flows and pump-induced flows. These flowmeters provide the measurements of pump-induced vertical flows which are analyzed to obtain vertical variations in horizontal hydraulic conductivity, K(z). With discrete areal K-values, K(x, y), and vertical profiles of K, provided by multiwell testing, the essential elements are present to produce a three-dimensional hydraulic conductivity field. The advent of these new flow measuring devices has contributed much to the motivation behind this paper. This paper presents the results of applying deterministic and stochastic methodology to the three-dimensional interpolation of hydraulic properties, specifically, hydraulic conductivity, K. Three of the approaches applied in this paper are deterministic in nature, inverse-distance weighting, inverse-distance-squared weighting, and ordinary kriging, while the fourth is a stochastic approach based on self-affine fractals. All of the methods are applied to measured data collected from 14 wells at a site in the United States near Mobile, Alabama. The three-dimensional K-distributions generated by each of the methods are used as inputs to an advective based transport model with the resulting model output compared to a two-well tracer study run previously at the same site.     

Measurement of ice flow velocities from GPS positions logged by short-period seismographs in East Antarctica

The ice flow velocity is a basic feature of glaciers and ice sheets. Measuring ice flow velocities is very important for estimating the mass balance of ice sheets in the Arctic and Antarctic. Traditional methods for measuring ice flow velocity include the use of stakes, snow pits and on-site geodetic GPS and remote sensing measurement methods. Geodetic GPS measurements have high accuracy, but geodetic GPS monitoring points only sparsely cover the Antarctic ice sheets. Moreover, the resolution and accuracy of ice flow velocities based on remote sensing measurements are low. Although the accuracy of the location data recorded by the navigation-grade GPS receivers embedded in short-period seismographs is not as good as that of geodetic GPS,the ice flow velocity can be accurately measured by these navigation-grade GPS data collected over a sufficiently long period. In this paper, navigation-grade GPS location data obtained by passive seismic observations during the 36 th Chinese National Antarctic Research Expedition were used to accurately track the movement characteristics of the ice sheet in the Larsemann Hills of East Antarctica and the Taishan Station area. The results showed that the ice sheet in the two study areas is basically moving northwestward with an average ice flow velocity of approximately 1 m mon-1. The results in the Taishan Station area are basically consistent with the geodetic GPS results, indicating that it is feasible to use the embedded GPS location data from shortperiod seismographs to track the movement characteristics of ice sheets. The ice flow characteristics in the Larsemann Hills are more complex. The measured ice flow velocities in the Larsemann Hills with a resolution of 200 m help to understand its characteristics. In summary, the ice flow velocities derived from GPS location data are of great significance for studying ice sheet dynamics and glacier mass balance and for evaluating the systematic errors caused by ice sheet movements in seismic imaging.     

Numerical simulation of three‐dimensional flow hydraulics in a braided channel

Results are presented from a numerical simulation of three‐dimensional flow hydraulics around a mid‐channel bar carried out using the FLUENT/UNS computational fluid dynamics (CFD) software package. FLUENT/UNS solves the three‐dimensional Reynolds‐averaged form of the Navier–Stokes equations. Turbulence closure is achieved using a RNG k–ϵ model. Simulated flow velocities are compared with measured two‐dimensional velocities (downstream and cross‐stream) obtained using an electromagnetic current meter (ECM). The results of the simulation are qualitatively consistent with the flow structures observed in the field. Quantitative comparison of the simulated and measured velocity magnitudes indicates a strong positive correlation between the two (r=0·88) and a mean difference of 0·09 m s⁻¹. Deviations between simulated and measured velocities may be identified that are both random and systematic. The former may reflect a number of factors including subgrid‐scale natural spatial variability in flow velocities associated with local bed structures and measurement uncertainty resulting from problems of ECM orientation. Model mesh configuration, roughness parameterization and inlet boundary condition uncertainty may each contribute to systematic differences between simulated and measured flow velocities. These results illustrate the potential for using CFD software to simulate flow hydraulics in natural channels with complex configurations. They also highlight the need for detailed spatially distributed datasets of three‐dimensional flow variables to establish the accuracy and applicability of CFD software. Copyright © 1999 John Wiley & Sons, Ltd.     

Use of an electromagnetic seepage meter to investigate temporal variability in lake seepage

A commercially available electromagnetic flowmeter is attached to a seepage cylinder to create an electromagnetic seepage meter (ESM) for automating measurement of fluxes across the sediment/water interface between ground water and surface water. The ESM is evaluated through its application at two lakes in New England, one where water seeps into the lake and one where water seeps out of the lake. The electromagnetic flowmeter replaces the seepage-meter bag and provides a continuous series of measurements from which temporal seepage processes can be investigated. It provides flow measurements over a range of three orders of magnitude, and contains no protruding components or moving parts. The ESM was used to evaluate duration of seepage disturbance following meter installation and indicated natural seepage rates resumed approximately one hour following meter insertion in a sandy lakebed. Lakebed seepage also varied considerably in response to lakebed disturbances, near-shore waves, and rainfalls, indicating hydrologic processes are occurring in shallow lakebed settings at time scales that have largely gone unobserved.     

Soil erosion by surface water flow on a stony,semiarid hillslope

Soil erosion on hillslopes occurs by processes of soil splash from raindrop impacts and sediment entrainment by surface water flows. This study investigates the process of soil erosion by surface water flow on a stony soil in a semiarid environment. A field experimental method was developed whereby erosion by concentrated flow could be measured in predefined flow areas without disturbing the soil surface. The method allowed for measurements in this study of flow erosion at a much wider range of slopes (2·6 to 30·1 per cent) and unit discharge rates (0·0007 to 0·007 m2 s−1) than have been previously feasible. Flow velocities were correlated to discharge and hydraulic radius, but not to slope. The lack of correlation between velocity and slope might have been due to the greater rock cover on the steeper slopes which caused the surface to be hydraulically rougher and thus counteract the expected effect of slope on flow velocity. The detachment data illustrated limitations in applying a linear hydraulic shear stress model over the entire range of the data collected. Flow detachment rates were better correlated to a power function of either shear stress (r2 = 0·51) or stream power (r2 = 0·59). Published in 1999 by John Wiley & Sons, Ltd.     

Seasonal vegetation and management influence overland flow velocity and roughness in upland grasslands

There is considerable interest in how headwater management may influence downstream flood peaks in temperate humid regions. However, there is a dearth of data on flow velocities across headwater hillslopes and limited understanding of whether surface flow velocity is influenced by seasonal changes in roughness through vegetation cycles or management. A portable hillslope flume was used to investigate overland flow velocities for four common headwater grassland habitats in northern England: Low-density Grazing, Hay Meadow, Rank Grassland and Juncus effusus Rush pasture. Overland flow velocity was measured in replicate plots for each habitat, in response to three applied flow rates, with the experiments repeated during five different periods of the annual grassland cycle. Mean annual overland flow velocity was significantly lower for the Rank Grassland habitat (0.026 m/s) followed by Low-density Grazing and Rushes (0.032 and 0.029 m/s), then Hay Meadows (0.041 m/s), which had the greatest mean annual velocity (examples from 12 L/min flow rate). Applying our mean overland flow velocities to a theoretical 100 m hillslope suggests overland flow is delayed by >1 hr on Rank Grassland when compared to Hay Meadows in an 18 mm storm. Thus grassland management is important for slowing overland flow and delaying peak flows across upland headwaters. Surface roughness was also strongly controlled by annual cycles of vegetation growth, decay, grazing and cutting. Winter overland flow velocities were significantly higher than in summer, varying between 0.004 m/s (Rushes, November) and 0.034 m/s (Rushes, June); and velocities significantly increased after cutting varying between 0.006 m/s (Hay meadows, July) and 0.054 m/s (Hay meadows, September). These results show that seasonal vegetation change should be incorporated into flood modelling, as cycles of surface roughness in grasslands strongly modify overland flow, potentially having a large impact on downstream flood peak and timing. Our data also showed that Darcy-Weisbach roughness approximations greatly over-estimated measured flow velocities.