Hyperscale terrain modelling of braided rivers: fusing mobile terrestrial laser scanning and optical bathymetric mapping

Quantifying the morphology of braided rivers is a key task for understanding braided river behaviour. In the last decade, developments in geomatics technologies and associated data processing methods have transformed the production of precise, reach‐scale topographic datasets. Nevertheless, generating accurate Digital Elevation Models (DEMs) remains a demanding task, particularly in fluvial systems. This paper identifies a threefold set of challenges associated with surveying these dynamic landforms: complex relief, inundated shallow channels and high rates of sediment transport, and terms these challenges the ‘morphological’, ‘wetted channel’ and ‘mobility’ problems, respectively. In an attempt to confront these issues directly, this paper presents a novel survey methodology that combines mobile terrestrial laser scanning and non‐metric aerial photography with data reduction and surface modelling techniques to render DEMs from the resulting very high resolution datasets. The approach is used to generate and model a precise, dense topographic dataset for a 2.5 km reach of the braided Rees River, New Zealand. Data were acquired rapidly between high flow events and incorporate over 5 x 10⁹ raw survey observations with point densities of 1600 pts m^‐2 on exposed bar and channel surfaces. A detailed error analysis of the resulting sub‐metre resolution is described to quantify DEM quality across the entire surface model. This reveals unparalleled low vertical errors for such a large and complex surface model; between 0.03 and 0.12 m in exposed and inundated areas of the model, respectively. Copyright © 2013 John Wiley & Sons, Ltd.     

Geostatistical analysis of ground‐survey elevation data to elucidate spatial and temporal river channel change

A digital elevation model (DEM) of a fluvial environment represented landform surface variability well and provided a medium for monitoring morphological change over time. Elevation was measured above an arbitrary datum using a ground‐based three‐dimensional tacheometric survey in two reaches of the River Nent, UK, in July 1998, October 1998 (after flood conditions) and June 1999. A detailed geostatistical analysis of the elevation data was used to model the spatial variation of elevation and to produce DEMs in each reach and for each survey period. Maps of the difference in elevation were produced and volumetric change was calculated for each reach and each survey period. The parameters of variogram models were used to describe the morphological character of each reach and to elucidate the linkages between process and the form of channel change operating at different spatial and temporal scales. The analysis of channel change on the River Nent shows the potential of geostatistics for investigating the magnitude and frequency of geomorphic work in other rivers. A flood modified the channel features, but low magnitude and high frequency flows rationalized the morphology. In spite of relatively small amounts of net flux the channel features changed as a consequence of the reworking of existing material. The blocking of chute entrances and redirection of the channel had a considerable effect on the behaviour of the channel. Such small changes suggested that the distributary system was sensitive to variation in sediment regime. Plots of the kriging variances against sampling intervals were used to quantify the temporal variation in sampling redundancy (ranging between ?11 per cent and +93 per cent). These curves illustrated the importance of bespoke sampling designs to reduce sampling effort by incorporating anisotropic variation in space and geomorphic information on flow regime. Variation in the nugget parameter of the variogram models was interpreted as sampling inaccuracy caused by variability in particle size and is believed to be important for future work on surface roughness. Copyright © 2003 John Wiley & Sons, Ltd.     

Modelling the effect of vegetation on channel pattern in bedload rivers

We modify a simple numerical stream‐pattern model to examine the effect of sediment stabilization by roots on the channel pattern of bedload rivers. In the model, vegetation enhances bank resistance to erosion, causing the development of a single channel instead of a rapidly changing, multiple channel (braided) pattern. Net aggradation resulting from a high sediment supply, however, causes frequent avulsions that destroy vegetation locally, leading to the development of a multiple‐channel pattern. A stability diagram representing multiple model runs predicts whether a river will exhibit single or multiple channels, based on plant‐enhanced bank strength, and on the time scale of plant development relative to a time scale for change in unvegetated channels. A second stability diagram predicts the way in which the amplitude and period of a fluctuating imposed sediment load influence whether a single or multiple‐channel pattern develops. Copyright © 2003 John Wiley & Sons, Ltd.     

Comparison of hydraulic geometry between sand‐ and gravel‐bed rivers in relation to channel pattern discrimination

A comparison has been made between the hydraulic geometry of sand‐ and gravel‐bed rivers, based on data from alluvial rivers around the world. The results indicate a signi?cant difference in hydraulic geometry among sand‐ and gravel‐bed rivers with different channel patterns. On this basis, some diagrams for discrimination of meandering and braided channel patterns have been established. The relationships between channel width and water discharge, between channel depth and water discharge, between width–depth ratio and water discharge and between channel slope and water discharge can all be used for channel pattern discrimination. The relationship between channel width and channel depth can also be used for channel pattern discrimination. However, the accuracy of these relationships for channel pattern discrimination varies, and the depth–discharge relationship is a better discriminator of pattern type than the classic slope–discharge function. The cause for this difference has been explained qualitatively. To predict the development of channel patterns under different natural conditions, the pattern discriminator should be searched on the basis of independent or at least semi‐independent variables. The relationship between stream power and bed material grain size can be used to discriminate channel patterns, which shows a better result than the discriminator using the slope–discharge relationship. Copyright © 2004 John Wiley & Sons, Ltd.     

Geostatistical interpolation of space–time rainfall on Tamshui River basin,Taiwan

Taiwan suffers from heavy storm rainfall during the typhoon season. This usually causes large river runoff, overland flow, erosion, landslides, debris flows, loss of power, etc. In order to evaluate storm impacts on the downstream basin, a real‐time hydrological modelling is used to estimate potential hazard areas. This can be used as a decision‐support system for the Emergency Response Center, National Fire Agency Ministry, to make ‘real‐time’ responses and minimize possible damage to human life and property. This study used 34 observed events from 14 telemetered rain‐gauges in the Tamshui River basin, Taiwan, to study the spatial–temporal characteristics of typhoon rainfall. In the study, regionalized theory and cross‐semi‐variograms were used to identify the spatial‐temporal structure of typhoon rainfall. The power form and parameters of the cross‐semi‐variogram were derived through analysis of the observed data. In the end, cross‐validation was used to evaluate the performance of the interpolated rainfall on the river basin. The results show the derived rainfall interpolator represents the observed events well, which indicates the rainfall interpolator can be used as a spatial‐temporal rainfall input for real‐time hydrological modelling. Copyright © 2007 John Wiley & Sons, Ltd.     

Computational fluid dynamics modelling of boundary roughness in gravel‐bed rivers: an investigation of the effects of random variability in bed elevation

Results from a series of numerical simulations of two‐dimensional open‐channel flow, conducted using the computational fluid dynamics (CFD) code FLUENT, are compared with data quantifying the mean and turbulent characteristics of open‐channel flow over two contrasting gravel beds. Boundary roughness effects are represented using both the conventional wall function approach and a random elevation model that simulates the effects of supra‐grid‐scale roughness elements (e.g. particle clusters and small bedforms). Results obtained using the random elevation model are characterized by a peak in turbulent kinetic energy located well above the bed (typically at y/h = 0·1–0·3). This is consistent with the field data and in contrast to the results obtained using the wall function approach for which maximum turbulent kinetic energy levels occur at the bed. Use of the random elevation model to represent supra‐grid‐scale roughness also allows a reduction in the height of the near‐bed mesh cell and therefore offers some potential to overcome problems experienced by the wall function approach in flows characterized by high relative roughness. Despite these benefits, the results of simulations conducted using the random elevation model are sensitive to the horizontal and vertical mesh resolution. Increasing the horizontal mesh resolution results in an increase in the near‐bed velocity gradient and turbulent kinetic energy, effectively roughening the bed. Varying the vertical resolution of the mesh has little effect on simulated mean velocity profiles, but results in substantial changes to the shape of the turbulent kinetic energy profile. These findings have significant implications for the application of CFD within natural gravel‐bed channels, particularly with regard to issues of topographic data collection, roughness parameterization and the derivation of mesh‐independent solutions. Copyright © 2001 John Wiley & Sons, Ltd.     

Random controls on semi‐rhythmic spacing of pools and riffles in constriction‐dominated rivers

Average pool spacing between five and seven bankfull widths has been documented in environments throughout the world, but has limited theoretical justification in coarse‐bedded and bedrock environments. Pool formation in coarse‐bedded and bedrock channels has been attributed to bedrock and boulder constrictions. Because the spacing of these constrictions may be irregular in nature, it is difficult to reconcile pool‐formation processes with the supposedly rhythmic spacing of pools and riffles. To address these issues, a simulation model for pool and riffle formation is used to demonstrate that semi‐rhythmic spacing of pools with an approximate spacing of five to seven bankfull widths can be recreated from a random distribution of obstructions and minimum pool‐ and riffle‐length criteria. It is assumed that a pool–riffle couplet will achieve a minimum length based on dominant‐discharge conditions. Values for the minimum‐length assumption are based on field data collected in New England and California, while the theoretical basis relies on the demonstrated hydraulic response of individual pools to elongation. Results from the simulations show that the location of pools can be primarily random in character, but still assume an average spacing between four and eight bankfull widths for a variety of conditions. Field verification data generally support the model but highlight a highly skewed distribution of pool‐forming elements and pool spacing. The relation between pool spacing and bankfull widths is attributed to the common geometric response of these features to dominant‐discharge conditions. Copyright © 2001 John Wiley & Sons, Ltd.     

Effect of growing watershed imperviousness on hydrograph parameters and peak discharge

An increasing impervious area is quickly extending over the Wu‐Tu watershed due to the endless demands of the people. Generally, impervious paving is a major result of urbanization and more recently has had the potential to produce more enormous flood disasters than those of the past. In this study, 40 available rainfall–runoff events were chosen to calibrate the applicable parameters of the models and to determine the relationships between the impervious surfaces and the calibrated parameters. Model inputs came from the outcomes of the block kriging method and the non‐linear programming method. In the optimal process, the shuffled complex evolution method and three criteria were applied to compare the observed and simulated hydrographs. The tendencies of the variations of the parameters with their corresponding imperviousness were established through regression analysis. Ten cases were used to examine the established equations of the parameters and impervious covers. Finally, the design flood routines of various return periods were furnished through use of approaches containing a design storm, block kriging, the SCS model, and a rainfall‐runoff model with established functional relationships. These simulated flood hydrographs were used to compare and understand the past, present, and future hydrological conditions of the watershed studied. In the research results, the time to peak of flood hydrographs for various storms was diminished approximately from 11 h to 6 h in different decrements, whereas peak flow increased respectively from 127 m³ s^?1 to 629 m³ s^?1 for different storm intensities. In addition, this study provides a design diagram for the peak flow ratio to help engineers and designers to construct hydraulic structures efficiently and prevent possible damage to human life and property. Copyright © 2007 John Wiley & Sons, Ltd.     

The role of beach morphology on coastal cliff erosion under extreme waves

Erosion of hard‐rock coastal cliffs is understood to be caused by a combination of both marine and sub‐aerial processes. Beach morphology, tidal elevation and significant wave heights, especially under extreme storm conditions, can lead to variability in wave energy flux to the cliff‐toe. Wave and water level measurements in the nearshore under energetic conditions are difficult to obtain and in situ observations are rare. Here we use monthly cliff‐face volume changes detected using terrestrial laser scanning alongside beach morphological changes and modelled nearshore hydrodynamics to examine how exposed cliffs respond to changes in extreme wave conditions and beach morphology. The measurements cover the North Atlantic storms of 2013 to 2014 and consider two exposed stretches of coastline (Porthleven and Godrevy, UK) with contrasting beach morphology fronting the cliffs; a flat dissipative sandy beach at Godrevy and a steep reflective gravel beach at Porthleven. Beach slope and the elevation of the beach–cliff junction were found to influence the frequency of cliff inundation and the power of wave–cliff impacts. Numerical modelling (XBeach‐G) showed that under highly energetic wave conditions, i.e. those that occurred in the North Atlantic during winter 2013–2014, with H_s = 5.5 m (dissipative site) and 8 m (reflective site), the combination of greater wave height and steeper beach at the reflective site led to amplified wave run‐up, subjecting these cliffs to waves over four times as powerful as those impacting the cliffs at the dissipative site (39 kWm^‐1 compared with 9 kWm^‐1). This study highlighted the sensitivity of cliff erosion to extreme wave conditions, where the majority (over 90% of the annual value) of cliff‐face erosion ensued during the winter. The significance of these short‐term erosion rates in the context of long‐term retreat illustrates the importance of incorporating short‐term beach and wave dynamics into geomorphological studies of coastal cliff change. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Morphological changes on meander point bars associated with flow structure at different discharges

This study investigates the impact of flow structure of different discharges on meander point bar morphology. We carried out mobile and terrestrial laser scanning campaigns before and after a flood on two sandy‐bed point bars. Between the scans, the flow structure was examined using an Acoustic Doppler Current Profiler at three flow stages. The results indicated that a meander point bar both affects and in turn, is itself modified by the flow at different discharges. The lower flow stages also have a significant effect on point bar morphology, especially on deposition over the bar head. Secondary circulation is responsible for scroll bar formation on the point bar margin beyond the apex. Flow separation at the inner bank, by contrast, does not require secondary circulation, but is dependent on flow depth over the point bar. A sudden increase in depth beyond the point bar top causes decreased stream power over the bar tail. The flow separation and decreased stream power cause a slow flow zone and net deposition over point bar tail. The backwater effect, if evident, may strengthen the process. Thus, filling over the bar tail seems generic for point bars and independent on secondary flow. Chutes and chute bars, scroll bars, bar head filling and bar platform filling, by contrast, require special fluvio‐morphological circumstances discussed in this paper. Whilst this paper confirms that the three‐dimensional flow structure has a major effect on point bar morphology, the flow structure seems to depend on how the point bar affects the flow trajectory which, in turn, depends upon the flow stage. Finally, the shape of the bend and the grain size distribution control the impacts of the flow structure, leading to divergent morphologies of point bars with certain generic features. Copyright © 2012 John Wiley & Sons, Ltd.     

南海西南部深水水道的多波束地形与多道反射地震研究

基于高分辨率多波束地形与多道地震资料精细解释,发现了南海西南部陆坡区5条大型深水水道.这5条水道可分成两类：1、5号水道延伸短,横剖面呈V形,弯曲度低,属顺直水道,侵蚀作用主导;2-4号水道延伸长,横剖面呈U形,弯曲度高,属弯曲水道,侵蚀-沉积作用主导.水道走向NEE-NE,长度大于350 km的有2条,宽度最大达2 km,切割深度最大达150 m.推测这些水道是陆源物质重力流作用形成的,是该区陆源物质通过陆架输送到南海西南次海盆的重要通道.陆坡区这些水道基本发育天然堤.在下游区域第二长的2号水道谷底明显高于周缘海底,水道-天然堤沉积已逐渐演变为朵叶沉积.这些水道的深部对应基底隆起区,说明深部构造的晚期活动控制了水道的平面位置,重力流与深部构造晚期活动共同控制了水道的形成演化.     

Effects of urbanization variables on model parameters for watershed divisions

Abstract

A semi-distributed model with a parallel connection was applied to examine the effects of urbanization variables. Data were obtained from watershed divisions that were characterized by various degrees of urbanization. The mean rainfall was calculated using the kriging method. The model inputs were obtained by subtracting mean rainfall from Ф-index values, based on the spatially uniform loss assumption. Regression analysis was applied to determine the relationship between the parameters of 64 calibrations and urbanization variables among the divisions. The results showed that overland parameters produced more consistent change in response to imperviousness than to population. Conversely, the channel parameter was unaffected by changes in urbanization. The verification results of 46 cases showed that power linkage was a potential option for linking division parameters with the corresponding imperviousness based on four evaluation criteria. The changes in imperviousness on overland parameters show the hydrological effects of division urbanizations.

Editor D. Koutsoyiannis; Associate editor T. Wagener

Citation Chen, R., Chuang, W.-N., and Cheng, S., 2014. Effects of urbanization variables on model parameters for watershed divisions. Hydrological Sciences Journal, 59 (6), 1167–1183. http://dx.doi.org/10.1080/02626667.2014.910305     

Effect of spatial variation characteristics on contouring of design storm depth

The ordinary kriging method, a geostatistical interpolation technique, was applied for developing contour maps of design storm depth in northern Taiwan using intensity–duration–frequency (IDF) data. Results of variogram modelling on design storm depths indicate that the design storms can be categorized into two distinct storm types: (i) storms of short duration and high spatial variation and (ii) storms of long duration and less spatial variation. For storms of the first category, the influence range of rainfall depth decreases when the recurrence interval increases, owing to the increasing degree of their spatial independence. However, for storms of the second category, the influence range of rainfall depth does not change significantly and has an average of approximately 72 km. For very extreme events, such as events of short duration and long recurrence interval, we do not recommend usage of the established design storm contours, because most of the interstation distances exceed the influence ranges. Our study concludes that the influence range of the design storm depth is dependent on the design duration and recurrence interval and is a key factor in developing design storm contours. Copyright © 2003 John Wiley & Sons, Ltd.     

Assessment of geomorphological bank evolution of the alluvial threshold rivers based on entropy concept parameters

The complex stream bank profiles in alluvial channels and rivers that are formed after reaching equilibrium has been a popular topic of research for many geomorphologists and river engineers. The entropy theory has recently been successfully applied to this problem. However, the existing methods restrict the further application of the entropy parameter to determine the cross-section slope of the river banks. To solve this limitation, we introduce a novel approach in the extraction of the equation based on the calculation of the entropy parameter (λ) and the transverse slope of the bank profile at threshold channel conditions. The effects of different hydraulic and geometric parameters are evaluated on a variation of the entropy parameter. Sensitivity analysis on the parameters affecting the entropy parameter shows that the most effective parameter on the λ-slope multiplier is the maximum slope of the bank profile and the dimensionless lateral distance of the river banks.     

Field evidence for the influence of weathering on rock erodibility and channel form in bedrock rivers

Erosion processes in bedrock‐floored rivers shape channel cross‐sectional geometry and the broader landscape. However, the influence of weathering on channel slope and geometry is not well understood. Weathering can produce variation in rock erodibility within channel cross‐sections. Recent numerical modeling results suggest that weathering may preferentially weaken rock on channel banks relative to the thalweg, strongly influencing channel form. Here, we present the first quantitative field study of differential weathering across channel cross‐sections. We hypothesize that average cross‐section erosion rate controls the magnitude of this contrast in weathering between the banks and the thalweg. Erosion rate, in turn, is moderated by the extent to which weathering processes increase bedrock erodibility. We test these hypotheses on tributaries to the Potomac River, Virginia, with inferred erosion rates from ~0.1 m/kyr to >0.8 m/kyr, with higher rates in knickpoints spawned by the migratory Great Falls knickzone. We selected nine channel cross‐sections on three tributaries spanning the full range of erosion rates, and at multiple flow heights we measured (1) rock compressive strength using a Schmidt hammer, (2) rock surface roughness using a contour gage combined with automated photograph analysis, and (3) crack density (crack length/area) at three cross‐sections on one channel. All cross‐sections showed significant (p < 0.01 for strength, p < 0.05 for roughness) increases in weathering by at least one metric with height above the thalweg. These results, assuming that the weathered state of rock is a proxy for erodibility, indicate that rock erodibility varies inversely with bedrock inundation frequency. Differences in weathering between the thalweg and the channel margins tend to decrease as inferred erosion rates increase, leading to variations in channel form related to the interplay of weathering and erosion rate. This observation is consistent with numerical modeling that predicts a strong influence of weathering‐related erodibility on channel morphology. Copyright © 2017 John Wiley & Sons, Ltd.     

Gully morphological characteristics in the loess hilly‐gully region based on 3D laser scanning technique

Gully morphology characteristics can be used effectively to describe the status of gully development. The Chabagou watershed, located in the hilly‐gully region of the Loess Plateau in China, was selected to investigate gully morphological characteristics using a 3D laser scanning technique (LIDAR). Thirty‐one representative gullies located at different watershed locations and gully orders were chosen to quantitatively describe gully morphology and establish empirical equations for estimating gully volume based on gully length and gully surface area. Images and point cloud data for the 31 gullies were collected, and digital elevation models (DEMs) with 10‐cm resolution were generated. ArcGIS 10.1 was then used to extract fundamental gully morphological parameters covering gully length (L), gully width (W_T) and gully depth (D), and some derivative morphological parameters, including gully head curvature (C), gully width–depth ratio (w/d), gully bottom‐to‐top width ratio (W_B/W_T), gully surface area (A_g) and gully volume (V_g). The results indicated that gullies in the upper watershed and the second order were more developed based on their high values of gully head curvature. The potential for gully development increased from the second order to the fourth order. Within the same gully orders, gullies in the lower watershed were more active with more development potential. A method for differentiating between gully head and gully sidewalls based on the gully head curvature value was proposed with a mean relative error of 8.77%. U‐shaped cross‐sections were widely distributed in the upper watershed and upper positions of a gully, while V‐shaped cross‐sections were widely distributed in the lower watershed and lower positions of a gully. V–L and V–A_g empirical equations with acceptable accuracy were established and can be used to estimate gully erosion in the Loess hilly‐gully region. Copyright © 2017 John Wiley & Sons, Ltd.     

Use of terrestrial photosieving and airborne topographic LiDAR to assess bed grain size in large rivers: a study on the Rhine River

Most grain size monitoring is still being conducted by manual sampling in the field, which is time consuming and has low spatial representation. Due to new remote sensing methods, some limitations have been partly overcome, but methodological progress is still needed for large rivers as well as in underwater conditions. In this article, we tested the reliability of two methods along the Old Rhine River (France/Germany) to estimate the grain size distribution (GSD) in above-water conditions: (i) a low-cost terrestrial photosieving method based on an automatic procedure using Digital Grain Size (DGS) software and (ii) an airborne LiDAR topo-bathymetric survey. We also tested the ability of terrestrial photosieving to estimate the GSD in underwater conditions. Field pebble counts were performed to compare and calibrate both methods. The results showed that the automatic procedure of terrestrial photosieving is a reliable method to estimate the GSD of sediment patches in both above-water and underwater conditions with clean substrates. Sensitivity analyses showed that environmental conditions, including solar lighting conditions and petrographic variability, significantly influence the GSD from the automatic procedure in above-water conditions. The presence of biofilm in underwater conditions significantly altered the GSD estimation using the automatic procedure, but the proposed manual procedure overcame this problem. The airborne LiDAR topographic survey is an accurate method to estimate the GSD of above-water bedforms and is able to generate grain size maps. The combination of terrestrial photosieving and airborne topographic LiDAR methods is adapted to assess the GSD over several kilometers long reaches of large rivers. © 2020 John Wiley & Sons, Ltd.     

The impact of stress history on bed structure

Recent research has started to focus on how prolonged periods of sub‐threshold flows may be capable of imparting structural changes that contribute to increased bed stability. To date, this effect (termed ‘stress history’) has been found to be significant in acting to increase a bed's critical shear stress at entrainment threshold. However, it is supported by only limited, qualitative and often speculative information on the mechanisms of this stabilization process in grade‐specific studies. As such, this paper uses high resolution laser scanning to quantitatively ascertain the granular mechanics underpinning the relationship between stress history and entrainment threshold for beds of a range of grain size distributions. Employing a bed slope of 1/200, three grain size distributions with median grain sizes (D₅₀) of 4·8 mm [uniform (σ_g = (D₈₄/D₁₆)^0.5 = 1·13; bimodal (σ_g = 2·08); and, unimodal (σ_g = 1·63)] were exposed to antecedent stress histories of 60 and 960 minutes duration. Antecedent shear stress magnitude was set at 50% of the critical shear stress for the D₅₀ when no stress history period was employed. Two laser displacement scans of the bed surface (approximate area 100 mm × 117 mm) were taken, one prior to the antecedent period and one after this period, so that changes to surface topography could be quantified (resolution of x = 0·10 mm, y = 0·13 mm and z = 0·24 mm). Rearrangement of bed surface structure is described using statistical analysis and two‐dimensional (2D) semi‐variograms to analyse scaling behaviour. Results reveal vertical settlement, changes to bed roughness and particle repositioning. However, the bed grain size distribution influences the relative importance of each mechanism in determining stress history induced bed stability; this is the focus of discussion in this paper. Copyright © 2012 John Wiley & Sons, Ltd.     

基于无损多源异构数据融合的水利工程隐蔽结构隐患探测技术

本文介绍了水利工程外露结构和隐蔽结构的不同检测方法和特点以及相关研究现状,总结了之前水利工程隐蔽结构隐患检测探索经验,结合横琴新区澳门大学中心沟东闸翼墙断裂探测等工程应用为例,通过三维激光扫描技术、探地雷达技术相结合的多源无损探测技术为例,对水利工程隐蔽结构的多源异构探测数据的融合应用做了相应探索研究。