Comparison of LiDAR waveform processing methods for very shallow water bathymetry using Raman,near‐infrared and green signals

Airborne light detection and ranging (LiDAR) bathymetry appears to be a useful technology for bed topography mapping of non‐navigable areas, offering high data density and a high acquisition rate. However, few studies have focused on continental waters, in particular, on very shallow waters (<2 m) where it is difficult to extract the surface and bottom positions that are typically mixed in the green LiDAR signal. This paper proposes two new processing methods for depth extraction based on the use of different LiDAR signals [green, near‐infrared (NIR), Raman] of the SHOALS‐1000T sensor. They have been tested on a very shallow coastal area (Golfe du Morbihan, France) as an analogy to very shallow rivers. The first method is based on a combination of mathematical and heuristic methods using the green and the NIR LiDAR signals to cross validate the information delivered by each signal. The second method extracts water depths from the Raman signal using statistical methods such as principal components analysis (PCA) and classification and regression tree (CART) analysis. The obtained results are then compared to the reference depths, and the performances of the different methods, as well as their advantages/disadvantages are evaluated. The green/NIR method supplies 42% more points compared to the operator process, with an equivalent mean error (?4·2 cm verusu ?4·5 cm) and a smaller standard deviation (25·3 cm verusu 33·5 cm). The Raman processing method provides very scattered results (standard deviation of 40·3 cm) with the lowest mean error (?3·1 cm) and 40% more points. The minimum detectable depth is also improved by the two presented methods, being around 1 m for the green/NIR approach and 0·5 m for the statistical approach, compared to 1·5 m for the data processed by the operator. Despite its ability to measure other parameters like water temperature, the Raman method needed a large amount of reference data to provide reliable depth measurements, as opposed to the green/NIR method. Copyright © 2010 John Wiley & Sons, Ltd.     

Qualitative and quantitative applications of LiDAR imagery in fluvial geomorphology

The potential for geomorphological mapping and quantitative calculations of light detection and ranging (LiDAR) data within fluvial geomorphology was studied for two river catchments within Belgium (Dijle and Amblève), which differ in physical settings and floodplain morphology. Two commercial, of‐the‐shelf LiDAR datasets with different specifications (horizontal resolution and vertical accuracy) were available for parts of the floodplains of both catchments. Real‐time kinematic (RTK) Global Positioning System (GPS) data were used as ground truth for error calculations. Qualitative analysis of LiDAR data allowed the identification of former channel patterns, levees, colluvial hillslope and fan deposits. These results were confirmed by field data, topographic surveys and historical maps. The pixel resolution proved to be an important factor in the identification of small landforms: only features with a width equal to or larger than LiDAR resolution can be detected. This poses limits on the usability of regionally available LiDAR data, which often have a horizontal resolution of several metres. The LiDAR data were also used in a quantitative analysis of channel dynamics. In the study area, the width of the Dijle River channel increased 3 m on average between 1969 and 2003. A sediment budget of channel processes for the period 1969–2003 indicated a total river bank erosion of 16·1 10³ m³ and a total within channel deposition of 7·1 10³ m³, resulting in a net river erosion of 9·0 10³ m³ or c. 0·4 Mg year^?1 per metre river length. Sequential LiDAR data can in theory be used to calculate vertical sedimentation rates, as long as there is control on the error of the reference levels used. Copyright © 2008 John Wiley and Sons, Ltd.     

Improvement of streams hydro‐geomorphological assessment using LiDAR DEMs

Hydro‐geomorphological assessments are an essential component for riverine management plans. They usually require costly and time‐consuming field surveys to characterize the spatial variability of key variables such as flow depth, width, discharge, water surface slope, grain size and unit stream power throughout the river corridor. The objective of this research is to develop automated tools for hydro‐geomorphological assessments using high‐resolution LiDAR digital elevation models (DEMs). More specifically, this paper aims at developing geographic information system (GIS) tools to extract channel slope, width and discharge from 1 m‐resolution LiDAR DEMs to estimate the spatial distribution of unit stream power in two contrasted watersheds in Quebec: a small agricultural stream (Des Fèves River) and a large gravel‐bed river (Matane River). For slope, the centreline extracted from the raw LiDAR DEM was resampled at a coarser resolution using the minimum elevation value. The channel width extraction algorithm progressively increased the centerline from the raw DEM until thresholds of elevation differences and slopes were reached. Based on the comparison with over 4000 differential global positioning system (GPS) measurements of the water surface collected in a 50 km reach of the Matane River, the longitudinal profile and slope estimates extracted from the raw and resampled LiDAR DEMs were in very good agreement with the field measurements (correlation coefficients ranging from 0 · 83 to 0 · 87) and can thus be used to compute stream power. The extracted width also corresponded very well to the channel as seen from ortho‐photos, although the presence of bars in the Matane River increased the level of error in width estimates. The estimated maximum unit stream power spatial patterns corresponded well with field evidence of bank erosion, indicating that LiDAR DEMs can be used with confidence for initial hydro‐geomorphological assessments. Copyright © 2013 John Wiley & Sons, Ltd.     

Quantification of braided river channel change using archival digital image analysis

Historical archives of grey‐scale river channel imagery are extensive. Here, we present and test a methodology to extract detailed quantitative topographic data from such imagery of sand‐bed rivers. Extracting elevation information from rivers is difficult as they are characterized by a low relative relief (<4 m); the area of interest may be spatially extensive (e.g. active channel widths >500 m in large braided rivers); the rate of change of surface elevation is generally low except in the vicinity of individual channel banks where the rate of change is very high; there is the complication that comes from inundation; and there may be an added complication caused by blockage of the field of view by vegetation. Here, we couple archival photogrammetric techniques with image processing methods and test these for quantification of sand‐bed braided river dynamics, illustrated for a 500 m wide, 3 km long reach of the South Saskatchewan River, Canada. Digital photogrammetry was used to quantify dry areas and water edge elevations. A methodology was then used to calibrate the spectral signature of inundated areas by combining established two media digital photogrammetric methods and image matching. This allowed determination of detailed depth maps for inundated areas and, when combined with dry area data, creation of complete digital elevation models. Error propagation methods were used to determine the erosion and deposition depths detectable from sequential digital elevation models. The result was a series of elevation models that demonstrate the potential for acquiring detailed and precise elevation data from any historical aerial imagery of rivers without needing associated calibration data, provided that imagery is of the necessary scale to capture the features of interest. We use these data to highlight several aspects of channel change on the South Saskatchewan River, including bar movement, bank erosion and channel infilling. Copyright © 2010 John Wiley & Sons, Ltd.     

Evaluating the performance of topobathymetric LiDAR to support multi-dimensional flow modelling in a gravel-bed mountain stream

Stream biophysical processes are commonly studied using multi-dimensional numerical modelling that quantifies flow hydraulics from which parameters such as habitat suitability, stream carrying capacity, and bed mobility are derived. These analyses would benefit from accurate high-resolution stream bathymetries spanning tens of kilometres of channel, especially in small streams or where navigation is difficult. Traditional ground-based survey methods are limited by survey time, dense vegetation and stream access, and are usually only feasible for short reaches. Conversely, airborne topobathymetric LiDAR surveys may overcome these limitations, although limited research is available on how errors in LiDAR-derived digital elevation models (DEMs) might propagate through flow models. This study investigated the performance of LiDAR-derived topobathymetry in support of multi-dimensional flow modelling and ecohydraulics calculations in two gravel-bedded reaches (approximately 200 m long), one morphologically complex and one morphologically simple, and at the segment scale (32 km-long stream segment) along a 15 m-wide river in central Idaho, USA. We compared metre and sub-metre-resolution DEMs generated from RTK-GPS ground and Experimental Advanced Airborne Research LiDAR-B (EAARL-B) surveys and water depths, velocities, shear stresses, habitat suitability, and bed mobility modelled with two-dimensional (2D) hydraulic models supported by LiDAR and ground-surveyed DEMs. Residual statistics, bias (B), and standard deviation (SD) of the residuals between depth and velocity predicted from the model supported by LiDAR and ground-survey topobathymetries were up to −0.04 (B) and 0.09 m (SD) for depth and −0.09 (B) and 0.20 m s⁻¹ (SD) for velocity. The accuracy (B = 0.05 m), precision (SD = 0.09 m), and point density (1 point m⁻²) of the LiDAR topobathymetric survey (regardless of reach complexity) were sufficient to support 2D hydrodynamic modelling and derivative stream habitat and process analyses, because these statistics were comparable to those of model calibration with B = 0 m and SD = 0.04 m for water surface elevation and B = 0.05 m s⁻¹ and SD = 0.22 m s⁻¹ for velocity in our investigation. © 2020 John Wiley & Sons, Ltd.     

Quantifying submerged fluvial topography using hyperspatial resolution UAS imagery and structure from motion photogrammetry

Quantifying the topography of rivers and their associated bedforms has been a fundamental concern of fluvial geomorphology for decades. Such data, acquired at high temporal and spatial resolutions, are increasingly in demand for process‐oriented investigations of flow hydraulics, sediment dynamics and in‐stream habitat. In these riverine environments, the most challenging region for topographic measurement is the wetted, submerged channel. Generally, dry bed topography and submerged bathymetry are measured using different methods and technology. This adds to the costs, logistical challenges and data processing requirements of comprehensive river surveys. However, some technologies are capable of measuring the submerged topography. Through‐water photogrammetry and bathymetric LiDAR are capable of reasonably accurate measurements of channel beds in clear water. While the cost of bathymetric LiDAR remains high and its resolution relatively coarse, the recent developments in photogrammetry using Structure from Motion (SfM) algorithms promise a fundamental shift in the accessibility of topographic data for a wide range of settings. Here we present results demonstrating the potential of so called SfM‐photogrammetry for quantifying both exposed and submerged fluvial topography at the mesohabitat scale. We show that imagery acquired from a rotary‐winged Unmanned Aerial System (UAS) can be processed in order to produce digital elevation models (DEMs) with hyperspatial resolutions (c. 0.02 m) for two different river systems over channel lengths of 50–100 m. Errors in submerged areas range from 0.016 m to 0.089 m, which can be reduced to between 0.008 m and 0.053 m with the application of a simple refraction correction. This work therefore demonstrates the potential of UAS platforms and SfM‐photogrammetry as a single technique for surveying fluvial topography at the mesoscale (defined as lengths of channel from c.10 m to a few hundred metres). Copyright © 2014 John Wiley & Sons, Ltd.     

Test of statistical means for the extrapolation of soil depth point information using overlays of spatial environmental data and bootstrapping techniques

Hydrological modelling depends highly on the accuracy and uncertainty of model input parameters such as soil properties. Since most of these data are field surveyed, geostatistical techniques such as kriging, classification and regression trees or more sophisticated soil‐landscape models need to be applied to interpolate point information to the area. Most of the existing interpolation techniques require a random or regular distribution of points within the study area but are not adequate to satisfactorily interpolate soil catena or transect data. The soil landscape model presented in this study is predicting soil information from transect or catena point data using a statistical mean (arithmetic, geometric and harmonic mean) to calculate the soil information based on class means of merged spatial explanatory variables. A data set of 226 soil depth measurements covering a range of 0–6·5 m was used to test the model. The point data were sampled along four transects in the Stubbetorp catchment, SE‐Sweden. We overlaid a geomorphology map (8 classes) with digital elevation model‐derived topographic index maps (2–9 classes) to estimate the range of error the model produces with changing sample size and input maps. The accuracy of the soil depth predictions was estimated with the root mean square error (RMSE) based on a testing and training data set. RMSE ranged generally between 0·73 and 0·83 m ± 0·013 m depending on the amount of classes the merged layers had, but were smallest for a map combination with a low number of classes predicted with the harmonic mean (RMSE = 0·46 m). The results show that the prediction accuracy of this method depends on the number of point values in the sample, the value range of the measured attribute and the initial correlations between point values and explanatory variables, but suggests that the model approach is in general scale invariant. Copyright © 2009 John Wiley & Sons, Ltd.     

Remote measurement of river morphology via fusion of LiDAR topography and spectrally based bathymetry

This study developed and evaluated a hybrid approach to remote measurement of river morphology that combines LiDAR topography with spectrally based bathymetry. Comparison of filtered LiDAR point clouds with surveyed cross‐sections indicated that subtle features on low‐relief floodplains were accurately resolved by LiDAR but that submerged areas could not be detected due to strong absorption of near‐infrared laser pulses by water. The reduced number of returns made the active channel evident in a LiDAR point density map. A second dataset suggested that pulse intensity also could be used to discriminate land from water via a threshold‐based masking procedure. Fusion of LiDAR and optical data required accurate co‐registration of images to the LiDAR, and we developed an object‐oriented procedure for achieving this alignment. Information on flow depths was derived by correlating pixel values with field measurements of depth. Highly turbid conditions dictated a positive relation between green band radiance and flow depth and contributed to under‐prediction of pool depths. Water surface elevations extracted from the LiDAR along the water's edge were used to produce a continuous water surface that preserved along‐channel variations in slope. Subtracting local flow depths from this surface yielded estimates of the bed elevation that were then combined with LiDAR topography for exposed areas to create a composite representation of the riverine terrain. The accuracy of this terrain model was assessed via comparison with detailed field surveys. A map of elevation residuals showed that the greatest errors were associated with underestimation of pool depths and failure to capture cross‐stream differences in water surface elevation. Nevertheless, fusion of LiDAR and passive optical image data provided an efficient means of characterizing river morphology that would not have been possible if either dataset had been used in isolation. Copyright © 2011 John Wiley & Sons, Ltd.     

Development of a river delta: a case study of Cimanuk river mouth,Indonesia

This paper describes delta development processes with particular reference to Cimanuk Delta in Indonesia. Cimanuk river delta, the most rapidly growing river delta in Indonesia, is located on the northern coast of Java Island. The delta is subject to ocean waves of less than 1 m height due to its position in the semi‐enclosed Java Sea in the Indonesian archipelago. The study has been carried out using a hydrodynamic model that accounts for sediment movement through the rivers and estuaries. As an advanced approach to management of river deltas, a numerical model, namely MIKE‐21, is used as a tool in the management of Cimanuk river delta. From calibration and verification of hydrodynamic model, it was found that the best value of bed roughness was 0·1 m. For the sediment‐transport model, the calibration parameters were adjusted to obtain the most satisfactory results of suspended sediment concentration and volume of deposition. By comparing the computed and observed data in the calibration, the best values of critical bed shear stress for deposition, critical bed shear stress for erosion and erosion coefficient were 0·05 N m^?2, 0·15 N m^?2, and 0·00001 kg m^?2 s^?1, respectively. The calibrated model was then used to analyse sensitivity of model parameters and to simulate delta development during the periods 1945–1963 and 1981–1997. It was found that the sensitive model parameters were bed shear stresses for deposition and erosion, while the important model inputs were river suspended sediment concentration, sediment characteristics and hydrodynamic. The model result showed reasonable agreement with the observed data. As evidenced by field data, the mathematical model proves that the Cimanuk river delta is a river‐dominated delta because of its protrusion pattern and very high sediment loads from the Cimanuk river. It was concluded that 86% of sediment load from the Cimanuk river was deposited in the Cimanuk delta. Copyright © 2008 John Wiley & Sons, Ltd.     

The topographic data source of digital terrain models as a key element in the accuracy of hydraulic flood modelling

The effects of the topographic data source and resolution on the hydraulic modelling of floods were analysed. Seven digital terrain models (DTMs) were generated from three different altimetric sources: a global positioning system (GPS) survey and bathymetry; high‐resolution laser altimetry data LiDAR (light detection and ranging); and vectorial cartography (1:5000). Hydraulic results were obtained, using the HEC‐RAS one‐dimensional model, for all seven DTMs. The importance of the DTM's accuracy on the hydraulic modelling results was analysed within three different hydraulic contexts: (1) the discharge and water surface elevation results from the hydraulic model; (2) the delineation of the flooded area; and (3) the relative sensitivity of the hydraulic model to changes in the Manning's n roughness coefficient. The contour‐based DTM was the least accurate with a root mean square error (RMSE) of 4·5 m in the determination of the water level and a variation of up to 50 per cent in the estimation of the inundated area of the floodplain. The GPS‐based DTM produced more realistic water surface elevation results and variations of up to 8 per cent in terms of the flooded area. The laser‐based model's RMSE for water level was 0·3 m, with the flooded area varying by less than 1 per cent. The LiDAR data also showed the greatest sensitivity to changes in the Manning's roughness coefficient. An analysis of the effect of mesh resolution indicated an influence on the delineation of the flooded area with variations of up to 7·3 per cent. In addition to determining the accuracy of the hydraulic modelling results produced from each DTM, an analysis of the time–cost ratio of each topographic data source illustrates that airborne laser scanning is a cost‐effective means of developing a DTM of sufficient accuracy, especially over large areas. Copyright © 2006 John Wiley & Sons, Ltd.     

Evaluating the capabilities of the CASI hyperspectral imaging system and Aquarius bathymetric LiDAR for measuring channel morphology in two distinct river environments

Remote sensing is a powerful tool for examining river morphology. This study used detailed field surveys to assess the capability of the CASI hyperspectral imaging system and Aquarius bathymetric LiDAR to measure bed elevations in rivers with disparate optical characteristics. Field measurements of water column optical properties in the clear Snake River, the more complex Blue and Colorado, and highly turbid Muddy Creek were used to calculate depth retrieval precision and dynamic range. Differences in depth of a few centimeters were detectable via passive optical techniques in the clearest stream, but precision was greatly reduced under turbid conditions. The bathymetric LiDAR evaluated in this study could not detect shallow depths or differences in depth smaller than 11 cm owing to the difficulty of distinguishing water surface and bottom returns in laser waveforms. In clear water and with high radiometric resolution, hyperspectral systems such as CASI could detect depths approaching 10 m, but semi‐empirical analysis of the Aquarius LiDAR indicated that maximum detectable depths were of the order of 2–3 m in the clear‐flowing Snake River, and closer to 1 m in the more turbid streams. Turbidity also constrained spectrally based depth retrieval, and depth estimates from the Blue/Colorado were far less reliable than on the Snake. Both sensors yielded positively biased (0.03 m for CASI, 0.08 m for Aquarius) bed elevations on the Snake, with precisions of 0.16–0.17 m. For the Blue/Colorado, mean errors were of the order of 0.2 m, biased shallow for optical data and biased deep for LiDAR, although no Aquarius laser returns were recorded from the deepest parts of these channels; precisions were reduced to 0.29–0.32 m. Both approaches have advantages and limitations, and prospective users must understand the capabilities and constraints associated with various types of remote sensing to ensure efficient use of these evolving technologies. Copyright © 2015 John Wiley & Sons, Ltd.     

Evaluation of ASTER GDEM and SRTM and their suitability in hydraulic modelling of a glacial lake outburst flood in southeast Tibet

The resolution and accuracy of digital elevation models (DEMs) can affect the hydraulic simulation results for predicting the effects of glacial lake outburst floods (GLOFs). However, for the Tibetan Plateau, high‐quality DEM data are often not available, leaving researchers with near‐global, freely available DEMs, such as the Advanced Spaceborne Thermal Emission and Reflection Radiometer Global Digital Elevation Model (ASTER GDEM) and the Shuttle Radar Topography Mission data (SRTM) for hydraulic modelling. This study explores the suitability of these two freely available DEMs for hydraulic modelling of GLOFs. Our study focused on the flood plain of a potentially dangerous glacial lake in southeastern Tibet, to evaluate the elevation accuracy of ASTER GDEM and SRTM, and their suitability for hydraulic modelling of GLOFs. The elevation accuracies of ASTER GDEM and SRTM were first validated against field global position system (GPS) survey points, and then evaluated with reference to the relatively high precision of 1:50 000 scale DEM (DEM5) constructed from aerial photography. Moreover, the DEM5, ASTER GDEM and SRTM were used as basic topographic data to simulate peak discharge propagation, as well as flood inundation extent and depth in the Hydrologic Engineering Center's River Analysis System one‐dimensional hydraulic model. Results of the three DEM predictions were compared to evaluate the suitability of ASTER GDEM and SRTM for GLOF hydraulic modelling. Comparisons of ASTER GDEM and SRTM each with DEM5 in the flood plain area show root‐mean‐square errors between the former two as ± 15·4 m and between the latter two as ± 13·5 m. Although SRTM overestimates and ASTER GDEM underestimates valley floor elevations, both DEMs can be used to extract the elevations of required geometric data, i.e. stream centre lines, bank lines and cross sections, for flood modelling. However, small errors still exist in the cross sections that may influence the propagation of peak discharge. The flood inundation extent and mean water depths derived from ASTER GDEM predictions are only 2·2% larger and 2·3‐m deeper than that of the DEM5 predictions, whereas the SRTM yields a flood zone extent 6·8% larger than the DEM5 prediction and a mean water depth 2·4‐m shallower than the DEM5 prediction. The modelling shows that, in the absence of high‐precision DEM data, ASTER GDEM or SRTM DEM can be relied on for simulating extreme GLOFs in southeast Tibet. Copyright © 2011 John Wiley & Sons, Ltd.     

Photogrammetric monitoring of small streams under a riparian forest canopy

The recent advent of digital photogrammetry has enabled the modeling and monitoring of river beds at relatively high spatial resolution (0·01 to 1 m) through the extraction of digital elevation models (DEMs). The traditional approach to image capture has been to mount a metric camera to an aircraft, although non‐metric cameras have been mounted to a variety of novel aerial platforms to acquire river‐based imagery (e.g. helicopters, radio‐controlled motorized vehicles, tethered blimps and balloons). However, most of these techniques are designed to acquire imagery at flying heights above the riparian tree canopy. In relatively narrow channels (e.g. <20 m bankfull width), streamside trees can obscure the channel and limit continuous photogrammetric data acquisition of both the channel bed and banks, while still providing useful information regarding the riparian canopy and even spot elevations of the channel. This paper presents a technique for the capture and analysis of close‐range photogrammetric data acquired from a vertically mounted non‐metric camera suspended 10 m above the channel bed by a unipod. The camera is positioned under the riparian forest canopy so that the channel bed can be imaged without obstruction. The system is portable and permits relatively rapid image acquisition over rough terrain and in dense forest. The platform was used to generate DEMs with a nominal ground resolution of 0·03 m. DEMs generated from this platform required post‐possessing to either adjust or eliminate erroneous cells introduced by the extraction process, overhanging branches, and by the effects of refraction at the air–water interface for submerged portions of the channel bed. The vertical precision in the post‐processed surface generally ranged from ± 0·01 to 0·1 m depending on the quality of triangulation and the characteristics of the surface being imaged. Copyright © 2010 John Wiley & Sons, Ltd.     

Annual bed‐elevation regime in the alluvial channel of Squamish River,Southwestern British Columbia,Canada

The aim of this study is to examine the annual regime of channel scour and fill by monitoring bed‐elevation changes in a reach of Squamish River in southwestern British Columbia, Canada. Sonar surveys of 13 river cross‐sections in a sandy gravel‐bed single‐channel study reach were repeated biweekly over a full hydrologic year (1995/6). The survey results show that bedload movement occurs as waves or pulses forming bedwaves that appear to maintain an overall coherence with movement downstream. These bedwaves propagate downstream by a mode here termed pulse scour and pulse fill, a process distinguished from the conventional mode of scour and fill commonly associated with flood events (here termed local scour and local fill). Bedwave celerity was estimated to be about 15·5 m d⁻¹ corresponding to a bedwave residence time in the study reach of almost one hydrologic year. The total amount of local bed‐elevation change ranged between 0·22 m and 2·41 m during the period of study. Analysis of the bed‐elevation and flow data reveals that, because of the bedwave phenomenon, there is no simple relation between the mean bed‐elevation and discharge nor any strong linear correlation among cross‐sectional behaviour. The bed‐elevation data also suggest that complex changes to the bed within a cross‐section are masked when the bed is viewed in one dimension, although no definitive trends in bed behaviour were found in the two‐dimensional analysis. Although a weak seasonal effect is evident in this study, the bed‐elevation regime is dominated by sediment supply‐driven fluctuations in bedload transport occurring at timescales shorter than the seasonal fluctuation in discharge. The study also indicates that bed‐elevation monitoring on Squamish River, and others like it, for purposes of detecting and measuring aggradation/degradation must take into account very considerable and normal channel‐bed variability operating at timescales from hours to months. Copyright © 2000 John Wiley & Sons, Ltd.     

The spatial and temporal patterns of aggradation in a temperate,upland, gravel‐bed river

Intensive field monitoring of a reach of upland gravel‐bed river illustrates the temporal and spatial variability of in‐channel sedimentation. Over the six‐year monitoring period, the mean bed level in the channel has risen by 0·17 m with a maximum bed level rise of 0·5 m noted at one location over a five month winter period. These rapid levels of aggradation have a profound impact on the number and duration of overbank flows with flood frequency increasing on average 2·6 times and overbank flow time increasing by 12·8 hours. This work raises the profile of coarse sediment transfer in the design and operation of river management, specifically engineering schemes. It emphasizes the need for the implementation of strategic monitoring programmes before engineering work occurs to identify zones where aggradation is likely to be problematic. Exploration of the sediment supply and transfer system can explain patterns of channel sedimentation. The complex spatial, seasonal and annual variability in sediment supply and transfer raise uncertainties into the system's response to potential changes in climate and land‐use. Thus, there is a demand for schemes that monitor coarse sediment transfer and channel response. Copyright © 2009 John Wiley & Sons, Ltd.     

In situ quantification of spatial and temporal variability of hyporheic exchange in static and mobile gravel‐bed rivers

Seepage meters modified for use in flowing water were used to directly measure rates of exchange between surface and subsurface water in a gravel‐ and cobble bed river in western Pennsylvania, USA (Allegheny River, Q_mean = 190 m³/s) and a sand‐ and gravel‐bed river in Colorado, USA (South Platte River, Q_mean = 9·7 m³/s). Study reaches at the Allegheny River were located downstream from a dam. The bed was stable with moss, algae, and river grass present in many locations. Median seepage was + 0·28 m/d and seepage was highly variable among measurement locations. Upward and downward seepage greatly exceeded the median seepage rate, ranging from + 2·26 (upward) to ? 3·76 (downward) m/d. At the South Platte River site, substantial local‐scale bed topography as well as mobile bedforms resulted in spatial and temporal variability in seepage greatly in exceedence of the median groundwater discharge rate of 0·24 m/d. Both upward and downward seepage were recorded along every transect across the river with rates ranging from + 2·37 to ? 3·40 m/d. Despite a stable bed, which commonly facilitates clogging by fine‐grained or organic sediments, seepage rates at the Allegheny River were not reduced relative to those at the South Platte River. Seepage rate and direction depended primarily on measurement position relative to local‐ and meso‐scale bed topography at both rivers. Hydraulic gradients were small at nearly all seepage‐measurement locations and commonly were not a good indicator of seepage rate or direction. Therefore, measuring hydraulic gradient and hydraulic conductivity at in‐stream piezometers may be misleading if used to determine seepage flux across the sediment‐water interface. Such a method assumes that flow between the well screen and sediment‐water interface is vertical, which appears to be a poor assumption in coarse‐grained hyporheic settings. Copyright © 2011 John Wiley & Sons, Ltd.     

Stability,morphology and surface grain size patterns of channel bifurcation in gravel–cobble bedded anabranching rivers

This study presents the first detailed field‐based analysis of the morphology of bifurcations within anabranching cobble–gravel rivers. Bifurcations divide the flow of water and sediment into downstream anabranches, thereby influencing the characteristics of the anabranches and the longevity of river islands. The history, morphology, bed grain size, and flow vectors at five bifurcations on the Renous River, New Brunswick, Canada, were studied in detail. The angles of bifurcations within five anabranching rivers in the Miramichi basin were investigated. The average bifurcation angle was 47°, within the range of values cited for braided river bifurcations. Bifurcation angle decreased when anabranches were of similar length. Shields stresses in channels upstream of bifurcations were lower than reported values for braided rivers. Stable bifurcations displayed lower Shields stresses than unstable bifurcations, contrary to experimental results from braided river bifurcations. Bifurcations in anabranching rivers are stabilized by vegetation that slows channel migration and helps to maintain a uniform upstream flow field. The morphology of stable bifurcations enhances their stability. A large bar, shaped like a shallow ramp that increases in elevation to floodplain level, forms at stable bifurcations. Floodplains at stable bifurcations accrete upstream at rates between 0·9 and 2·5 m a^?1. Bars may also form within the entrance of an anabranch downstream of the bifurcation node. These bars are associated with bifurcation instability, forming after a period of stability or an avulsion. Channel abandonment occurs when a bar completely blocks the entrance to one anabranch. The stability of channels upstream of bifurcations and the location of bars at bifurcations influence bifurcation stability and the maintenance of river anabranching in the long term. Copyright © 2006 John Wiley & Sons, Ltd.     

3D float tracking: in situ floodplain roughness estimation

This paper presents a novel technique to quantify in situ hydrodynamic roughness of submerged floodplain vegetation: 3D float tracking. This method uses a custom‐built floating tripod that is released on the inundated floodplain and tracked from shore by a robotic total station. Simultaneously, an acoustic Doppler current profiler (ADCP) collects flow velocity profiles and water depth data. Roughness values are derived from two methods based on (1) run‐averaged values of water depth, slope and flow velocity to compute the roughness based on the Chézy equation, assuming uniform flow, (2) the equation for one‐dimensional free surface flow in a moving window. A sensitivity analysis using synthetic data proved that the median value of the roughness, derived using method 2, is independent of (1) the noise in water levels, up to 9 mm, (2) bottom surface slope, and (3) topographic undulations. The window size should be at least 40 m for a typical lowland river setup. Field measurements were carried out on two floodplain sections with an average vegetation height of 0·030 (Arnhem) and 0·043 m (Dreumel). Method 1 resulted in a Nikuradse roughness length of 0·08 m for both locations. Method 2 gave 0·12 m for Arnhem and 0·19 m for Dreumel. In Arnhem, a spatial pattern of roughness values was present, which might be related to fractional vegetation cover or vegetation density during the flood peak. 3D float tracking proved a flexible and detailed method for roughness determination in the absence of waves, and provided an unrestricted view from shore. Copyright © 2008 John Wiley & Sons, Ltd.     

Spatial prediction of river channel topography by kriging

Topographic information is fundamental to geomorphic inquiry, and spatial prediction of bed elevation from irregular survey data is an important component of many reach‐scale studies. Kriging is a geostatistical technique for obtaining these predictions along with measures of their reliability, and this paper outlines a specialized framework intended for application to river channels. Our modular approach includes an algorithm for transforming the coordinates of data and prediction locations to a channel‐centered coordinate system, several different methods of representing the trend component of topographic variation and search strategies that incorporate geomorphic information to determine which survey data are used to make a prediction at a specific location. For example, a relationship between curvature and the lateral position of maximum depth can be used to include cross‐sectional asymmetry in a two‐dimensional trend surface model, and topographic breaklines can be used to restrict which data are retained in a local neighborhood around each prediction location. Using survey data from a restored gravel‐bed river, we demonstrate how transformation to the channel‐centered coordinate system facilitates interpretation of the variogram, a statistical model of reach‐scale spatial structure used in kriging, and how the choice of a trend model affects the variogram of the residuals from that trend. Similarly, we show how decomposing kriging predictions into their trend and residual components can yield useful information on channel morphology. Cross‐validation analyses involving different data configurations and kriging variants indicate that kriging is quite robust and that survey density is the primary control on the accuracy of bed elevation predictions. The root mean‐square error of these predictions is directly proportional to the spacing between surveyed cross‐sections, even in a reconfigured channel with a relatively simple morphology; sophisticated methods of spatial prediction are no substitute for field data. Copyright © 2007 John Wiley & Sons, Ltd.     

Radio‐tracking gravel particles in a large braided river in New Zealand: a field test of the stochastic theory of bed load transport proposed by Einstein

Hans A. Einstein initiated a probabilistic approach to modelling sediment transport in rivers. His formulae were based on theory and were stimulated by laboratory investigations. The theory assumes that bed load movement occurs in individual steps of rolling, sliding or saltation and rest periods. So far very few attempts have been made to measure stochastic elements in nature. For the first time this paper presents results of radio‐tracing the travel path of individual particles in a large braided gravel bed river: the Waimakariri River of New Zealand. As proposed by Einstein, it was found that rest periods can be modelled by an exponential distribution, but particle step lengths are better represented by a gamma distribution. Einstein assumed an average travel distance of 100 grain‐diameters for any bed load particle between consecutive points of deposition, but larger values of 6·7 m or 150 grain‐diameters and 6·1 m or 120 grain‐diameters were measured for two test particle sizes. Together with other available large scale field data, a dependence of the mean step length on particle diameter relative to the D₅₀ of the bed surface was found. During small floods the time used for movement represents only 2·7% of the total time from erosion to deposition. The increase in percentage of time being used for transport means that it then has to be regarded in stochastic transport models. Tracing the flow path of bed load particles between erosion and deposition sites is a step towards explaining the interactions between sediment transport and river morphology. Copyright © 2001 John Wiley & Sons, Ltd.