Landscape Evolution Modelling of naturally dammed rivers

Natural damming of upland river systems, such as landslide or lava damming, occurs worldwide. Many dams fail shortly after their creation, while other dams are long‐lived and therefore have a long‐term impact on fluvial and landscape evolution. This long‐term impact is still poorly understood and landscape evolution modelling (LEM) can increase our understanding of different aspects of this response. Our objective was to simulate fluvial response to damming, by monitoring sediment redistribution and river profile evolution for a range of geomorphic settings. We used LEM LAPSUS, which calculates runoff erosion and deposition and can deal with non‐spurious sinks, such as dam‐impounded areas. Because fluvial dynamics under detachment‐limited and transport‐limited conditions are different, we mimicked these conditions using low and high erodibility settings, respectively. To compare the relative impact of different dam types, we evaluated five scenarios for each landscape condition: one scenario without a dam and four scenarios with dams of increasing erodibility. Results showed that dam‐related sediment storage persisted at least until 15 000 years for all dam scenarios. Incision and knickpoint retreat occurred faster in the detachment‐limited landscape than in the transport‐limited landscape. Furthermore, in the transport‐limited landscape, knickpoint persistence decreased with increasing dam erodibility. Stream capture occurred only in the transport‐limited landscape due to a persisting floodplain behind the dam and headward erosion of adjacent channels. Changes in sediment yield variation due to stream captures did occur but cannot be distinguished from other changes in variation of sediment yield. Comparison of the model results with field examples indicates that the model reproduces several key phenomena of damming response in both transport‐limited and detachment‐limited landscapes. We conclude that a damming event which occurred 15 000 years ago can influence present‐day sediment yield, profile evolution and stream patterns. Copyright © 2014 John Wiley & Sons, Ltd.     

How do streamflow generation mechanisms affect watershed hypsometry?

Several studies have shown that the dominant streamflow generation mechanism in a river basin can leave distinct geomorphological signatures in basin topography. In particular, it has been suggested previously that basins generated by groundwater discharge tend to have a larger hypsometric integral than surface runoff basins because fluvial erosion is more focused in the valleys where groundwater discharge tends to occur. In this analysis, we aim to clarify this relationship by developing an alternative method to quantify the effects of streamflow generation mechanisms on basin hypsometry and by using a numerical model that can generate streamflow by different processes to evaluate the sensitivity of the results to the hydrological and geomorphological properties of the basin. The model results suggest that the hypsometric characteristics that are usually associated with groundwater discharge basins, such as a larger hypsometric integral, occur primarily when drainage networks are still advancing in the watershed. During later stages of development, an additional factor such as lithological controls or a distinct geomorphological process would be needed to preserve these features. The model results also show that the hypsometric effects are stronger when the parameters of the fluvial erosion process promote the influence of small discharge rates. Copyright © 2007 John Wiley & Sons, Ltd.     

Tipping points in Texas rivers

Anticipating geomorphic tipping points requires that we learn from the past. Major geomorphic changes in coastal plain rivers of Texas resulting in river metamorphosis or regime shifts were identified and the major driving factors determined. Eleven such transformations – possible tipping points – were identified from contemporary observations, historical records, and Quaternary reconstructions. Two of the tipping points (between general aggrading and degrading valley states) are associated with reversals in a fundamental system control (sea‐level). One (stable or aggrading versus degrading channels) is associated with an abrupt change in sediment supply due to dam construction, and two others (changes from meandering to anastomosing channel patterns, and different anastomosis styles) are similarly related to changes in sediment supply and/or transport capacity, but with additional elements of historical contingency. Three tipping points are related to avulsions. One, from a regime dominated by re‐occupation of former channels to one dominated by progradation into flood basins, is driven by progressive long‐term filling of incised valleys. Another, nodal avulsions, is triggered by disturbances associated with tectonic uplift or listric faults. The third, avulsions and related valley metamorphosis in unfilled incised valleys, is due to fundamental dynamical instabilities within the fluvial system. This synthesis and analysis suggests that geomorphic tipping points are sometimes associated with general extrinsic or intrinsic (to the fluvial system) environmental change, independent of any disturbances or instabilities. Others are associated with natural (e.g. tectonic) or human (dams) disturbances, and still others with intrinsic geomorphic instabilities. This suggests future tipping points will be equally diverse with respect to their drivers and dynamics. Copyright © 2018 John Wiley & Sons, Ltd.     

Minimum froude number and the equilibrium of alluvial sand rivers

Rivers adjust towards an equilibrium condition, the stability of which depends upon a set of controlling factors expressed by the Froude number. As alluvial river channels approach stable conditions, the Froude number of the channel flow will tend to attain a minimum value which reflects minimum bed material motion and maximum channel stability, under the constraints imposed by water discharge, sediment load, and particle size. Computer simulations for sand bed rivers show that the Froude number of the flow tends to a minimum value when the equilibrium river tends to a certain hydraulic geometry. Evidence from 57 alluvial sand material rivers and stable canals shows that this simulated hydraulic geometry with minimum Froude number corresponds to the natural equilibrium state.     

The influence of sediment size,relative grain size and channel slope on initiation of sediment motion in boulder bed rivers. A lichenometric study

Sediment transport of four boulder bed rivers is studied using lichenometry. The presence of lichens on boulders in the river channel is used to date the last mobilization of the blocks. Using size frequency diagrams and regional growth curves calibrated with dated reference points it is possible to determine the flood event responsible for the last mobilization of each boulder with lichens present. The specific stream power of flood events over the last 60 years is then calculated, and thresholds of sediment transport based on the sediment size are calculated. The results from the four studied rivers are compared to similar relationships in the literature. Sediment motion thresholds appear to be very variable within the same type of river (mountainous boulder bed rivers). The critical specific stream power necessary to mobilize a particle of a given diameter may vary by up to 10 times from one river to the next. Bed sediment size and river slope may explain this large range of stream powers. Calculation of the relative size of the transported particles (D_i/D₅₀) also shows that both hiding and protrusion effects, as well as channels slope, are important factors in sediment transport. Copyright © 2010 John Wiley & Sons, Ltd.     

Channel pattern of proglacial rivers: topographic forcing due to glacier retreat

Glacier retreat leads to changes in channel pattern during deglaciation, in response to changing water, sediment and base level controls. Recent ongoing retreat at Skaftafellsjökull, Iceland (c. 50 m per year since 1998) has resulted in the formation of a sequence of river terraces, and several changes in river channel pattern. This paper compares widely used models of river channel pattern against the changes observed at Skaftafellsjökull. Doing this reveals the role of topographic forcing in determining proglacial channel pattern, whilst examining the predictive power and limitations of the various approaches to classifying river channels. Topography was found to play a large role in determining channel pattern in proglacial environments for two reasons: firstly, glacier retreat forces rivers to flow through confined moraine reaches. In these reaches, channels which theory predicts should be braided are forced to adopt a single channel. Secondly, proximal incision of proglacial rivers, accompanied by downstream aggradation, leads to changes in slope which force the river to cross channel pattern thresholds. The findings of this work indicate that in the short term, the majority of channel pattern change in proglacial rivers is due to topographic forcing, and that changes due to changing hydrology and sediment supply are initially relatively minor, although likely to increase in significance as deglaciation progresses. These findings have implications for palaeohydraulic studies, where changes in proglacial channel pattern are frequently interpreted as being due to changes in water or sediment supply. This paper shows that channel pattern can change at timescales faster than hydrological or sediment budget changes usually occur, in association with relatively minor changes in glacier mass balance. Copyright © 2014 John Wiley & Sons, Ltd.     

Longitudinal distributions of river flood power: the combined automated flood,elevation and stream power (CAFES) methodology

Stream power can be an extremely useful index of fluvial sediment transport, channel pattern, river channel erosion and riparian habitat development. However, most previous studies of downstream changes in stream power have relied on field measurements at selected cross‐sections, which are time consuming, and typically based on limited data, which cannot fully represent important spatial variations in stream power. We present here, therefore, a novel methodology we call CAFES (combined automated flood, elevation and stream power), to quantify downstream change in river flood power, based on integrating in a GIS framework Flood Estimation Handbook systems with the 5 m grid NEXTMap Britain digital elevation model derived from IFSAR (interferometric synthetic aperture radar). This provides a useful modelling platform to quantify at unprecedented resolution longitudinal distributions of flood discharge, elevation, floodplain slope and flood power at reach and basin scales. Values can be resolved to a 50 m grid. CAFES approaches have distinct advantages over current methodologies for reach‐ and basin‐scale stream power assessments and therefore for the interpretation and prediction of fluvial processes. The methodology has significant international applicability for understanding basin‐scale hydraulics, sediment transport, erosion and sedimentation processes and river basin management. Copyright © 2008 John Wiley & Sons, Ltd.     

Semi-analytical evaluation of integral forms associated with Lamb's problem

The present work gives solutions of integrals resulting from Lamb's problem solution by an accurate and efficient approximation procedure. The formulation of the solutions requires the study of Green's function of Lamb problem solution. Unfortunately, this problem is associated to great mathematical difficulties and needs the evaluation of some complex integrals. We handled these integrals by an appropriate strategy based on decoupling their effects mathematically and solved by using the least squares method. The present solutions can be applied to study the dynamical interaction between soil and structures, to be used to determine the transient wave fields caused by a seismic source, and also to assess numerical computations with a different numerical methods program.     

Fluvial hydraulics of hyperconcentrated floods in Chinese rivers

Hyperconcentrated floods, with sediment concentrations higher than 200 kg/m³, occur frequently in the Yellow River and its tributaries on the Loess Plateau. This paper studies the fluvial hydraulics of hyperconcentrated floods by statistical analysis and comparison with low sediment concentration floods. The fluvial process induced by hyperconcentrated floods is extremely rapid. The river morphology may be altered more at a faster rate by one hyperconcentrated flood than by low sediment concentration floods over a decade. The vertical sediment concentration distribution in hyperconcentrated floods is homogeneous. The Darcy–Weisbach coefficient of hyperconcentrated floods varies with the Reynolds number in the same way as normal open channel flows but a representative viscosity is used to replace the viscosity, η. If the concentration is not extremely high and the Reynolds number is larger than 2000, the flow is turbulent and the Darcy–Weisbach coefficient for the hyperconcentrated floods is almost the same as low sediment concentration floods. Serious channel erosion, which is referred to as ‘ripping up the bottom’ in Chinese, occurs in narrow‐deep channels during hyperconcentrated floods. However, in wide‐shallow channels, hyperconcentrated floods may result in serious sedimentation. Moreover, a hyperconcentrated flood may cause the channel to become narrower and deeper, thus, reducing the flood stage by more than 1 m if the flood event lasts longer than one day. The fluvial process during hyperconcentrated floods also changes the propagation of flood waves. Successive waves may catch up with and overlap the first wave, thus, increasing the peak discharge of the flood wave during flood propagation along the river course. Copyright © 2009 John Wiley & Sons, Ltd.     

Knickpoint recession rate and catchment area: the case of uplifted rivers in Eastern Scotland

Knickpoint behaviour is a key to understanding both the landscape responses to a base‐level fall and the corresponding sediment fluxes from rejuvenated catchments, and must be accommodated in numerical models of large‐scale landscape evolution. Knickpoint recession in streams draining to glacio‐isostatically uplifted shorelines in eastern Scotland is used to assess whether knickpoint recession is a function of discharge (here represented by its surrogate, catchment area). Knickpoints are identified using DS plots (log slope versus log downstream distance). A statistically significant power relationship is found between distance of headward recession and catchment area. Such knickpoint recession data may be used to determine the values of m and n in the stream power law, E = KA^mSⁿ. The data have too many uncertainties, however, to judge definitively whether they are consistent with m = n = 1 (bedrock erosion is proportional to stream power and KPs should be maintained and propagate headwards) or m = 0·3, n = 0·7 (bedrock incision is proportional to shear stress and KPs do not propagate but degrade in place by rotation or replacement). Nonetheless, the E Scotland m and n values point to the dominance of catchment area (discharge) in determining knickpoint retreat rates and are therefore more consistent with the stream power law formulation in which bedrock erosion is proportional to stream power. Copyright © 2005 John Wiley & Sons, Ltd.     

Remotely sensed rivers in the Anthropocene: state of the art and prospects

The rivers of the world are undergoing accelerated change in the Anthropocene, and need to be managed at much broader spatial and temporal scales than before. Fluvial remote sensing now offers a technical and methodological framework that can be deployed to monitor the processes at work and to assess the trajectories of rivers in the Anthropocene. In this paper, we review research investigating past, present and future fluvial corridor conditions and processes using remote sensing and we consider emerging challenges facing fluvial and riparian research. We introduce a suite of remote sensing methods designed to diagnose river changes at reach to regional scales. We then focus on identification of channel patterns and acting processes from satellite, airborne or ground acquisitions. These techniques range from grain scales to landform scales, and from real time scales to inter-annual scales. We discuss how remote sensing data can now be coupled to catchment scale models that simulate sediment transfer within connected river networks. We also consider future opportunities in terms of datasets and other resources which are likely to impact river management and monitoring at the global scale. We conclude with a summary of challenges and prospects for remotely sensed rivers in the Anthropocene. © 2019 John Wiley & Sons, Ltd.     

Optical remote mapping of rivers at sub‐meter resolutions and watershed extents

At watershed extents, our understanding of river form, process and function is largely based on locally intensive mapping of river reaches, or on spatially extensive but low density data scattered throughout a watershed (e.g. cross sections). The net effect has been to characterize streams as discontinuous systems. Recent advances in optical remote sensing of rivers indicate that it should now be possible to generate accurate and continuous maps of in‐stream habitats, depths, algae, wood, stream power and other features at sub‐meter resolutions across entire watersheds so long as the water is clear and the aerial view is unobstructed. Such maps would transform river science and management by providing improved data, better models and explanation, and enhanced applications. Obstacles to achieving this vision include variations in optics associated with shadows, water clarity, variable substrates and target–sun angle geometry. Logistical obstacles are primarily due to the reliance of existing ground validation procedures on time‐of‐flight field measurements, which are impossible to accomplish at watershed extents, particularly in large and difficult to access river basins. Philosophical issues must also be addressed that relate to the expectations around accuracy assessment, the need for and utility of physically based models to evaluate remote sensing results and the ethics of revealing information about river resources at fine spatial resolutions. Despite these obstacles and issues, catchment extent remote river mapping is now feasible, as is demonstrated by a proof‐of‐concept example for the Nueces River, Texas, and examples of how different image types (radar, lidar, thermal) could be merged with optical imagery. The greatest obstacle to development and implementation of more remote sensing, catchment scale ‘river observatories’ is the absence of broadly based funding initiatives to support collaborative research by multiple investigators in different river settings. Copyright © 2007 John Wiley & Sons, Ltd.     

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.     

Changes in stream flow regime in headwater catchments of the Yellow River basin since the 1950s

The headwater catchments of the Yellow River basin generate over 35% of the basin's total stream flow and play a vital role in meeting downstream water resources requirements. In recent years the Yellow River has experienced significant changes in its hydrological regime, including an increased number of zero‐flow days. These changes have serious implications for water security and basin management. We investigated changes in stream flow regime of four headwater catchments since the 1950s. The rank‐based non‐parametric Mann–Kendall test was used to detect trends in annual stream flow. The results showed no significant trend for the period 1956 to 2000. However, change‐point analysis showed that a significant change in annual stream flow occurred around 1990, and hence the stream‐flow data can be divided into two periods: 1956–1990 and 1991–2000. There was a considerable difference in average annual stream flow between the two periods, with a maximum reduction of 51%. Wet‐season rainfall appears to be the main factor responsible for the decreasing trend in annual stream flow. Reductions in annual stream flow were associated with decreased interannual variability in stream flow. Seasonal stream flow distribution changed from bimodal to unimodal between the two periods, with winter stream flow showing a greater reduction than other seasons. Daily stream flow regime represented by flow duration curves showed that all percentile flows were decreased in the second period. The high flow index (Q₅/Q₅₀) reduced by up to 28%, whereas the reduction in the low flow index (Q₉₅/Q₅₀) is more dramatic, with up to 100% reduction. Copyright © 2006 John Wiley & Sons, Ltd.     

基于面平均雨量误差修正的实时洪水预报修正方法

空间集总式水文模型的洪水预报精度会受到面平均雨量估计误差的严重影响.点雨量监测值的误差类型、误差大小以及流域的雨量站点密度和站点的空间分布都会影响到面平均雨量的计算.为提高实时洪水预报精度,本文提出了一种基于降雨系统响应曲线洪水预报误差修正方法.通过此方法估计降雨输入项的误差,从而提高洪水预报精度.此方法将水文模型做为输入和输出之间的响应系统,用实测流量和计算流量之间的差值做为信息,通过降雨系统响应曲线,使用最小二乘估计原理,对面平均雨量进行修正,再用修正后的面平均雨量重新计算出流过程.将此修正方法结合新安江模型使用理想案例进行检验,并应用于王家坝流域的16场历史洪水以及此流域不同雨量站密度的情况下,结果证明均有明显修正效果,且在雨量站密度较低时修正效果更加明显.该方法是一种结构简单且不增加模型参数和复杂度的实时洪水修正的新方法.     

中小学校疏散演练人员疏散模型研究

本文以元胞自动机模型理论为指导,以人员疏散时间最短为目标,构建中小学校疏散演练人员疏散模型。本文从开展地震应急疏散模型研究出发,为地震应急演练疏散训练软件提供核心算法。通过对空间基础构建元素及搭建方法的分析研究,设计完成疏散虚拟场景的建模方法。最后开展人员疏散算法构建方式的研究,并结合应急疏散演练人员行为规范的研究成果,完成疏散模型的总体设计,达到为应急演练疏散软件提供人员疏散演化模型以及为人员疏散过程合理性评判提供依据的目的。     

Modelling temperature change downstream of forest harvest using Newton's law of cooling

We adapted Newton's law of cooling to model downstream water temperature change in response to stream‐adjacent forest harvest on small and medium streams (average 327 ha in size) throughout the Oregon Coast Range, USA. The model requires measured stream gradient, width, depth and upstream control reach temperatures as inputs and contains two free parameters, which were determined by fitting the model to measured stream temperature data. This model reproduces the measured downstream temperature responses to within 0.4 °C for 15 of the 16 streams studied and provides insight into the physical sources of site‐to‐site variation among those responses. We also use the model to examine how the pre‐harvest to post‐harvest change in daily maximum stream temperature depends on distance from the harvest reach. The model suggests that the pre‐harvest to post‐harvest temperature change approximately 300 m downstream of the harvest will range from roughly 82% to less than 1% of that temperature change that occurred within the harvest reach, depending primarily on the downstream width, depth and gradient. Using study‐averaged values for these channel characteristics, the model suggests that for a stream representative of those in the study, the temperature change approximately 300 m downstream of the harvest will be 56% of the temperature change that occurred within the harvest reach. This adapted Newton's law of cooling procedure represents a highly practical means for predicting stream temperature behaviour downstream of timber harvests relative to conventional heat budget approaches and is informative of the dominant processes affecting stream temperature. Copyright © 2015 John Wiley & Sons, Ltd.     

太湖流域印染业企业环境行为分析

在分析企业环境行为类型及决策机制的基础上,认为内外环境因素的综合作用促成了企业环境行为的发生.据此,通过相关性和主成分分析对太湖流域印染企业环境行为影响因素进行选择,并构建了环境行为决策数量模型,结果表明:私营合资企业环保投入水平高于国有集体企业;中型企业环保投入高于小型企业,此外,企业环境行为演变也受政府和经济环境制约.最后提出了优化企业环境行为的对策,即加快企业改制,强化用地供应管理,改革排污费制度以及规划产业布局.     

On the inadequacies of Long's model for steady two-dimensional stratified flows

Nonlinear analysis of two-dimensional steady flows with density stratification in the presence of gravity is considered. Inadequacies of Long's model for steady stratified flow over topography are explored. These include occurrence of closed streamline regions and waves propagating upstream. The usual requirements in Long's model of constant dynamic pressure and constant vertical density gradient in the upstream condition are believed to be the cause of these inadequacies. In this article, we consider a relaxation of these requirements, and also provide a systematic framework to accomplish this. As illustrations of this generalized formulation, exact solutions are given for the following two special flow configurations: the stratified flow over a barrier in an infinite channel; the stratified flow due to a line sink in an infinite channel. These solutions exhibit again closed-streamline regions as well as waves propagating upstream. The persistence of these inadequacies in the generalized Long's model appears to indicate that they are not quite consequences of the assumptions of constant dynamic pressure and constant vertical density gradient in Long's model, contrary to previous belief.

On the other hand, solutions admitted by the generalized Long's model show that departures from Long's model become small as the flow becomes more and more supercritical. They provide a nonlinear mechanism for the generation of columnar disturbances upstream of the obstacle and lead in subcritical flows to qualitatively different streamline topological patterns involving saddle points, which may describe the lee-wave-breaking process in subcritical flows and could serve as seats of turbulence in real flows. The occurrences of upstream disturbances in the presence of lee-wave-breaking activity described by the present solution are in accord with the experiments of Long (Long, R.R., “Some aspects of the flow of stratified fluids, Part 3. Continuous density gradients”, Tellus 7, 341--357 (1955)) and Davis (Davis, R.E., “The two-dimensional flow of a stratified fluid over an obstacle”, J. Fluid Mech. 36, 127–143 ()).