The goal … To collect channel velocity and cross-section data on Monday and operationally transmit real-time stage, velocity, area, and discharge on Tuesday regardless of:
The United States Geological Survey (USGS) and Environment Canada (EC) are exploring the use of velocity radars to measure surface-water velocities and compute real time mean-channel velocity and discharge in streams and rivers. Velocity radars can (1) deliver real-time discharge at new stations where stage-discharge, index-velocity, or slope-discharge ratings are not available; (2) extend ratings; (3) corroborate indirect measurements; and (4) provide an alternative for measuring discharge at sites with complex ratings . The ultimate goal of this effort is to transition the use of surface-water velocity radars from a proof-of-concept to an operational mode and to more clearly determine operational limitations.
These guidelines are an extension of the non-contact methods initiated by HYDRO 21 (Costa et al., 2006). The USGS Hydrologic Instrumentation Facility (HIF), Project Chiefs from 6 USGS Water Science Centers (WSCs), and two vendors (Stalker and Hydrological Services of America; HSA) participated in proof-of-concept testing. Fourteen (14) sites, which exhibited different hydraulic flow regimes, were identified and collocated with existing USGS streamgages. Laboratory testing was conducted in parallel with field deployments by the HIF (Fulford, 2015) and the Switzerland Federal Institute of Metrology (METAS). Testing at the HIF and METAS was conducted in the spring 2015 using standard methods including flumes, carriage-tow tanks, and tuning forks with known frequencies. Stalker and HSA offered their units gratis in exchange for an assessment on their performance. Similarly, OTT HydroMet (OTT) offered their radar for testing; however, the unit is still under development. Fulton et al. (2020) compiled a comprehensive radar analysis conducted at 10 collocated USGS streamgages.
You’ll need to acquire a small amount of velocity and channel data beyond what a normal field trip demands. Data should be collected at the cross section-of-interest and in the vicinity of the radar footprint. Keep in mind that the radar’s ability to return a surface-water velocity is influenced by (1) the quality of scatterers or waveforms on the water surface, (2) the air gap or the distance between the bridge deck and the water surface, and (3) the potential noise imposed by wind drift, eddies, secondary flows, and macro turbulence.
This document was designed to (1) evaluate whether radars can accurately measure surface-water velocities, (2) compute the mean velocity and discharge at a channel cross-section, (3) identify the environmental and hydraulic factors that influence surface-water velocity measurements, and (4) establish a protocol for transitioning the this proof-of-concept to an operational streamgage platform.
Velocity radars can be deployed by hand (figure 1) or fixed on bridges (figure 2), light cableways or cable stays (figure 3).
Channel, velocity and discharge data should be co-collected using methods consistent with Turnipseed and Sauer (2010) or Office of Surface Water Hydroacoustics webpage (http://hydroacoustics.usgs.gov/index.shtml).
Channel, velocity and discharge data measured at USGS streamgages will be used to validate the radar-derived data; however, when siting a surface-water velocity radar, the same hydrodynamic conditions used to site a conventional streamgage (Site Selection, p. 9; Turnipseed and Sauer, 2010) should be followed, which include:
It is important to note the stationing or location of the vertical (“termed the y-axis”) in a cross section where the maximum-instream or maximum-surface water velocity occurs. Generally, the maximum-surface-water velocity occurs at the same vertical as the maximum-instream velocity. The y-axis is also where all velocity and depth data should be collected to translate surface-water velocities into a mean-channel velocity. Velocity data can be collected using current meters, ADCPs, and ADVs. The following data will be recorded at each measured cross section:
The y-axis, which contains the maximum in-stream or surface-water velocity, should be recorded at the cross section-of-interest. Research indicates the y-axis is quite stable and does not vary with changing hydraulic conditions (Fulton, 2011; Fulton and Ostrowski, 2008; Chiu et al., 2001; Chiu and Chen, 1999; Fulton, 1999). It should be noted that y-axis rarely coincides with the channel thalweg in open or engineered channels. Generally, the location of the y-axis and the parameters used to compute discharge (see Surf Board link, https://my.usgs.gov/confluence/pages/viewpage.action?pageId=552933693) remain constant for a cross-section regardless of variations in stage, velocity, flow, flow, channel geometry, bed form and material, slope, or alignment; however, field verification of these parameters must be conducted periodically and a stage-area rating must be maintained. Although used for verification and confirmation, historical data such as stage, velocity, and discharge, which are commonly collected for a range of low and high flow events, is not needed, but it is advised. By following the procedures listed below, real-time mean-channel velocity and discharge can be computed straight away.
If the channel has not been surveyed, it is recommended that the cross section be surveyed and include the wetted perimeter and above-water points-of-interest in the floodplain. The stage-area rating is used to compute area, which is required to compute discharge. Horizontal and vertical control should be surveyed relative to the gage datum using a total station survey, GPS receiver, or an equivalent. The program AreaComp (Lant and Mueller, 2012) can be used to generate a synthetic stage-area rating, when estimating areas above the water surface during the day of the siting.
Current meters, ADVs, ADCPs, and surface-velocity radars can be used to site a radar. The Radarsubmittal_form. xlsx (https://drive.google.com/drive/folders/0BzD_MCxxmLFyRDZhUzV1a1dMd28) summarizes the channel and hydraulic data needs. It’s important to avoid wind-dominated reaches, eddies, secondary flows, and macro turbulence (figures 4 through 9) caused by piers when collecting velocity data. Generally, data collection and radar deployments point should be upstream of bridges or structures.
When wading is possible, measure velocities and compute discharge in accordance with Turnipseed and Sauer (2010). Select the y-axis from the 25 - 30 verticals comprising the measurement that exhibits the greatest velocity value based on either the maximum-surface velocity; 0.2D and 0.8D velocities; or 0.6D velocity. At the y-axis, measure the surface velocity using handheld or portable velocity radar concurrently with point velocities immediately below the water surface to the channel bottom at an interval that can be used to adequately define the velocity distribution along the selected vertical. Depending on water depth, this should include a minimum of 6 point velocities (near the channel bottom, 0.2D, 0.4D, 0.6D, 0.8D, close to the water surface while minimizing air entrainment); however, it is preferred that the surface-water velocity, point velocities near the water surface and close to the channel bottom, 0.2D, 0.3D, 0.4D, 0.5D, 0.6D, 0.7D, 0.8D, 0.9D be collected. Repeat this procedure to the left and right of the y-axis.
When wading is not possible, velocity data should be collected at the y-axis using a Stationary Moving Bed Analysis (SMBA) either by boat, light cableway, tethered from a bridge, or river banks. Surface-velocities should be collected concurrently with the ADCP measurement using handheld or portable velocity radar. Repeat this procedure to the left and right of the y-axis. Discharge should be computed using QRev (https://hydroacoustics.usgs.gov/movingboat/QRev.shtml) or an equivalent. When coupled with a GPS receiver, the lat/long of the y-axis should be recorded or established using VMT (Parsons, 2012). Process the velocity distribution in WinRiver II by choosing:
“Reprocess Selected Transect”
“Averaging Data” to identify the number of ensembles to average to reduce noise
“Tabular, Earth Velocity Magnitude and Direction” and “Graphs, Profile, Velocity”
The depth and velocity data can be copied to a text file.
Velocity radars are used to measure surface velocities and do not penetrate the water surface. Typically, the vertical containing the maximum-surface velocity will contain the maximum-instream velocity. Velocities should be measured relative to bridge stationing or geo-referenced using a GPS receiver. 20 to 25 surface-water velocities are needed to adequately identify the maximum-surface water velocity and y-axis. The velocity radar can be pointed upstream (preferred) or downstream from a bridge or walkway. It should be oriented parallel to flow lines and tilted (from horizontal) at a nominal 45-degree incidence angle. It should be noted that different radar units operate at different incidence angles. It’s important to note when collecting velocity data to avoid wind-dominated reaches, eddies, secondary flows, and macro turbulence.
Conventional methods including stage-discharge and index-velocity ratings, current meters/ADVs/mid-section method, and ADCPs are widely accepted as industry-standards (table 1). As a result, they will be used as a periodic benchmarks and should be conducted in conjunction with the radar-derived discharges for comparison purposes. Real-time radar results from the Tanana River at Nenana, AK and Rio Grande at Embudo, NM are presented in figures 10 and 11, respectively.
Table 1 Method summary.
QRev, WinRiver II or RiverSurveyor
Surface-water velocity radar
Probability Concept (see https://my.usgs.gov/confluence/pages/viewpage.action?pageId=552933693
Surface- velocity radars will not work at every site, particularly where wind drift is dominant and surface-water velocities are less than 0.5 fps and the surface-water scatterers (small wave forms) are ill-defined.
The scatterer size and quality is a function of the frequency of the transmitted radar and environmental factors such as wind drift, turbulence, and rain. The length of these small-scale surface waves or Bragg waves is computed using equation 1:
λb = λ/(2 sin θ) Eq. 1
λb= wavelength of the water wave or Bragg wave, Gigahertz; GHz
λ = wavelength of the transmitted radar signal; GHz
θ is the incidence angle
λ = c / v
c = speed of light = 3E8 meters per second, m/s
v = frequency in hertz; Hz
K-band radar transmission frequency = 24 GHz = 24E9 Hz
Incidence angle = 45 degrees
λb = 3E8 m/s x 100 cm/m / 24E9 Hz / (2 x sin 45 degrees) = 0.88 cm, which represents the small-scale surface waves that serve as scatterers. A wind-drift correction algorithm was developed and applied to the Red River of the North at Grand Forks, ND. The results are illustrated in figures 10 (raw discharge) and 11 (corrected discharge for wind drift).
The duration of sampling has an impact on radar returns and should be optimized in the field based on the spectrum returned from the radar (figure 12). Typically, sample durations range from 40 seconds to 2 minutes. A number of filtering schemes to reduce measurement noise such as high- and low-pass filters, moving averages, LOESS, Savitsky-Golay filter, and Kalman filter are being investigated.
Regardless of the equipment selected, velocity variations associated with turbulence create noise in the data, which complicates the analysis. The objective is to minimize scatter in the velocity data collected at the y-axis by reducing the standard deviation of the velocity data.
When using current meters or ADVs, turbulent flows that can dominate in natural and artificial channels are accompanied by local eddies (Rantz et al., 1982) that can result in variations in the velocity in any direction. It changes rapidly in time and space and is scale-dependent. Pierce (1941) reported that the velocities recorded by current meters in laboratory flumes were greater in magnitude (relative to the mean velocity) for lower velocities when compared to higher velocities. At high velocities, variations have a minor effect on current-meter observations. It’s customary to observe velocity at a point by current meter for a period that ranges from 40 to 70 seconds (Rantz et al., 1982).
When using ADCPs it is important to minimize noisy data related to ADCPs, this is achieved by collecting an even number of (reciprocal) transects with a minimum of 720 seconds total exposure time be made during steady-flow conditions. The measured discharge will be the average of the discharges from all valid reciprocal transects. Reciprocal transects should always be made to reduce potential directional biases. For policy detail, see OSW Technical Memorandum 2011.08. If using a TRDI product, the Correlation Profile should be reviewed to ensure the quality of the measurement is sufficient (reviewing the signal-to-noise ratio, SNR) to ensure the SNR is greater than 128. After locating the y-axis and to provide for redundancy, velocity data should be collected at the y-axis using a an SMBA either by boat or tethered from a bridge or river banks.
When using surface velocity radars and depending on variations in velocities with time, each vertical should be sampled for 40 seconds to 2 minutes. Please note if measurements are collected downstream of a bridge and depending on stage, piering can create secondary flow patterns that can influence the velocity distribution at the y-axis and ultimately the parameters used to compute discharge. It’s preferred that all radar and hydroacoustic measurements be collected upstream of bridges. Using ADCPs downstream of a bridge to compare discharge rates is acceptable. Diagnostic tests, which are run for current meters or acoustic instruments, are not available for velocity radars. However, tuning forks can be used to validate the velocities reported by a radar (equations 2 and 3). By striking the tuning fork and placing the tuning fork in front of the radar antenna, the recorded velocity is measured using the Doppler shift.
v = (f_Doppler Shift × c)/(2 x f_radar ) Eq. 2
v = velocity recorded by the radar, feet per second; fps
f Doppler shift = Doppler shift associated with tuning fork frequency, MHz
c = speed of light, miles per second, mi/s
f radar = radar frequency, gigahertz; GHz
f Doppler Shift = 1055 Hz tuning fork frequency, MHz
c = 186,000 mi/s
f radar = 24.150 GHz
v = (1055 MHz x 186,000 mi/s x 5,280 ft/mi)/(2 x 24.150 GHz x (10^9 Hz)/(1 GHz)) = 21.5 fps Eq. 3
Some portable and fixed-mount radars (Sommer Messtechnik) produce spectra, which offer a quantitative tool that serves as a “spin test” for electromagnetic instruments and can be used to qualify the value of surface-water velocity estimate. Figures 13 through 16 offer examples of good, fair, and poor spectra that can be used to assess the quality of the measurement.
Handheld radar (Stalker Pro II SVR) deployed from a bridge used to measure surface-water velocities concurrently while measuring point velocities with a FlowTracker at the y-axis in a cross section.
Fixed-mount radar (Sommer RQ30 and Stalker Surface Velocity Sensor) deployed from a bridge near upstream of USGS streamgage 05082500 Red River of the North at Grand Forks, ND and used to measure surface-water velocities and compute discharge. Courtesy of Chris Laveau.
Cable stay deployment (Sommer RQ30), where little or no infrastructure exists at USGS streamgages 385309104561101 Middle Waldo and 385254104560401 Lower Waldo. The radars are used to measure surface-water velocities and compute discharge. Data is transmitted via 3G and Iridium modems.
Good surface scatterers (Green for go!) and minimal wind drift. Poor surface scatterers (Red for stop!) downstream of piering. Fixed-mount radars (Sommer RQ30 and Stalker Surface Velocity Sensor) deployed from a bridge upstream of the USGS streamgage 08279500 Rio Grande at Embudo, NM. The radars are installed at the y-axis of the cross section and are used to measure surface-water velocities and compute discharge. Courtesy Jay Cederberg.
Good surface scatterers (Green for go!) and minimal wind drift. Poor surface scatterers (Red for stop!) caused by wind drift. Wind drift creates noise in the radar returns and must be corrected prior to transmitting surface-water velocities. Upstream of the USGS streamgage 15515500 Tanana River at Nenana, AK. Courtesy Heather Best.
Good surface scatterers (Green for go!) and minimal wind drift upstream of the pier, where the radar is pointing. Poor surface scatterers (Red for stop!) located downstream and adjacent to the pier caused by secondary flows and eddies. Fixed-mount radar (Sommer RQ30) deployed from a bridge upstream of the USGS streamgage 06192500 Yellowstone River near Livingston, MT. The radars are installed at the y-axis of the cross section and are used to measure surface-water velocities and compute discharge. Courtesy Steve Holnbeck.
Good surface scatterers (Green for go!) associated with pancake ice. Fixed-mount heterodyne radar (APL-UW RiverScat) is pointing upstream near USGS streamgage 01538700 Susquehanna River at Bloomsburg, PA.
Good surface scatterers (Green for go!) near the USGS streamgage 05056678 Tolna Coulee near Tolna, ND. Courtesy Chris Laveau.
Poor surface scatterers (Red for stop!) and significant wind drift near USGS 06751490 North Fork Cache La Poudre River at Livermore, CO.
Raw discharge data, Red River of the North, Grand Forks, ND (PROVISIONAL).
Wind-drift corrected discharge data, Red River of the North, Grand Forks, ND (PROVISIONAL).
Stage-discharge rating and instantaneous radar-derived discharge, Tanana River at Nenana, AK (PROVISIONAL). Noise is caused by wind drift.
Stage-discharge rating, radar-derived discharge and measured discharge, Rio Grande at Embudo, NM (PROVISIONAL).
Good spectra with a sharp peak and no velocities from opposite directions. Courtesy Wolfram Sommer (Sommer Messtechnik).
Fair spectra with small velocities from opposing directions, yet a significant peak in the direction of flow. Courtesy Wolfram Sommer (Sommer Messtechnik).
Poor spectra with velocities from a variety of directions and multiple peak velocities. The radar will recognize a velocity, but it will not be very accurate. Courtesy Wolfram Sommer (Sommer Messtechnik).
Chiu, C.-L., Tung, N.C., Hsu, S.M., and Fulton, J.W., 2001, Comparison and assessment of methods of measuring discharge in rivers and streams, Research Report No. CEEWR-4, Dept. of Civil & Environmental Engineering, University of Pittsburgh, Pittsburgh, PA.
Chiu, C.-L. and Chen, Y.-C., 1999, “Efficient methods of measuring discharge and reservoir-sediment inflow.” Risk Analysis in Dam Safety Assessment, ed. by J.-T. Kuo and B.-C. Yen, Water Resources Pub., Inc., Highlands Ranch, Colorado, 97-116.
Costa, J.E., Cheng, R.T., Haeni, F.P., Melcher, N., Spicer, K.R., Hayes, E., Plant, W., Hayes, K., Teague, C., and Barrick, D., 2006, Use of radars to monitor stream discharge by noncontact methods, Water Resources Research, 42, W07422, doi:10.1029/2005WR004430.
Fulford, J., draft 2revised 11 Nov 2015, Laboratory Testing of the Sommer RG-30. Internal USGS document.
Fulton, J.W.; Mason, C.A.; Eggleston, J.R.; Nicotra, M.J.; Chiu, C.-L.; Henneberg, M.F.; Best, H.R.; Cederberg, J.R.; Holnbeck, S.R.; Lotspeich, R.R.; Laveau, C.D.; Moramarco, T.; Jones, M.E.; Gourley, J.J.; Wasielewski, D. Near-Field Remote Sensing of Surface Velocity and River Discharge Using Radars and the Probability Concept at 10 U.S. Geological Survey Streamgages. Remote Sens. 2020, 12, 1296..
Fulton, J.W, 2011, Non-contact methods for computing streamflow as part of the 70th Anniversary of the Ministry of Water Resources of the People’s Republic of China, November 2011.
Fulton, J.W. and Ostrowski, J., 2008, Measuring real-time streamflow using emerging technologies: Radar, hydroacoustics, and the probability concept, Journal of Hydrology 357, 1– 10.
Fulton, J.W., 1999, Comparison of conventional and probability-based modeling of open-channel flow in the Allegheny River, Pennsylvania, USA, Unpublished Thesis, University of Pittsburgh, Department of Civil and Environmental Engineering.
Lant, J. and Mueller, D.S., May 2012, Stage Area Rating Application – AreaComp2, USGS, Office of Surface Water, accessed April 10, 2013 http://hydroacoustics.usgs.gov/software/AreaComp2_Users_Guide.pdf
Parsons, D.R., Jackson, P.R., Czuba, J.A., Engel, F.L., Rhoads, B.L., Oberg, K.A., Best, J.L., Mueller, D.S., Johnson, K.K., and Riley, J.D., 2012, Velocity Mapping Toolbox (VMT): a processing and visualization suite for moving-vessel ADCP measurements, Earth Surface Processes and Landforms, Published online in Wiley Online Library (wileyonlinelibrary.com) DOI: 10.1002/esp.3367.
Pierce, C. H., 1941, Investigations of methods and equipment used in stream gaging; Part 1, Performance of current meters in water of shallow depth: U.S. Geol. Survey Water-Supply Paper 868-A, 35 p.
Turnipseed, D.P., and Sauer, V.B., 2010, Discharge measurements at gaging stations: U.S. Geological Survey Techniques and Methods book 3, chap. A8, 87 p. (Also available at http://pubs.usgs.gov/tm/tm3-a8/.)
U.S. Geological Survey
Denver Federal Center
Box 25046, MS 415
Denver, CO 80225
U.S. Geological Survey
3400 Shell St.
Fairbanks, AK 99701
U.S. Geological Survey
Denver Federal Center
Box 25046, MS 415
Denver, CO 80225
|You may also want to use visual panels to communicate related information, tips or things users need to be aware of.|
Related articles appear here based on the labels you select. Click to edit the macro and add or change labels.