DOCSLIB.ORG

A Kolmogoroff-Type Scaling for the Fine Structure of Drainage Basins

Total Page：16

File Type：pdf, Size：1020Kb

Download full-text PDF Read full-text
A Kolmogoroff-Type Scaling for the Fine Structure of Drainage Basins A Kolmogoroff-type Scaling for the Fine Structure of Drainage Basins Gary Parker, University of Minnesota Peter K. Haff and A. Brad Murray, Duke University Kolmogoroff scaling in turbulent flows: Energy Cascade •Turbulent energy is produced in eddies scaling with the size of the "box" (e.g. river depth). •The nonlinear terms in the equations of momentum balance grind larger eddies into ever smaller eddies. •The grinding continues until the eddies are so small that the turbulent energy can be effectively dissipated into heat. •There are no smaller eddies. production of nonlinear transfer of dissipation of turbulent energy energy from large to turbulent energy to at large scales small scales heat at small scales Balance of energy in a steady turbulent flow ∂ (mean kinetic energy of turbulence) = ∂t production rate − dissipation rate + spatial transfer rate Let u = characteristic turbulent velocity (L/T) L = size of the container (e.g. river depth) (L) ν = viscosity of flow (L2/T) P = production rate/mass of turbulent energy (L2/T3) ε = dissipation rate/mass of turbulent energy (L2/T3) uk = Kolmogoroff velocity scale (L/T) ηk = Kolmogoroff length scale ∂u u3 ∂u′ ∂u′ u 2 P = −u′u′ i ~ ε = ν i i ~ ν k i j ∂x ∂x η2 ∂x j L j j k Make Kolmogoroff scales from the viscosity ν (L2/T) and the dissipation rate ε (L2/T3) 1/ 4 ⎛ ν3 ⎞ u η ⎜ ⎟ 1/ 4 k k ηk = ⎜ ⎟ uk = (νε) = 1 ⎝ ε ⎠ ν Recalling that u3 P ≅ ε ε ~ L it is found that 3 / 4 η ⎛ ν ⎞ k ~ ⎜ ⎟ L ⎝ uL ⎠ That is, the higher the Reynolds number of the flow the finer is the turbulence From Tennekes and Lumley Does the idea of a Kolmogoroff scale have any application to drainage basins? How dense does a drainage network have to be in order to "cover" the catchment? Well, how dense? Consider a (statistically) steady state landscape undergoing uplift at constant rate vu (L/T). The channels are undergoing incision at rate vI (L/T). The hillslopes contain a well-developed regolith which moves downslope with a kinematic diffusivity k (L2/T) HYPOTHESIS 1: Channel incision creates elevation fluctuations. HYPOTHESIS 2: Hillslope diffusion obliterates elevation fluctuations. HYPOTHESIS 3: The drainage network must be just sufficiently fine so that rate of creation of elevation fluctuations balances rate of obliteration. Channel incision in the near-absence of hillslope diffusion: the SLOT CANYON: creates amplitude fluctuation Hillslope diffusion with weak channel incision: MELTED ICE-CREAM MORPHOLOGY: destroys amplitude fluctuation Image courtesy Bill Dietrich Some scales Lu = length of headwater tributary 2 Au = αuLu = area of headwater catchment: αu ~o(1) Sh = characteristic slope of hillslope in headwater catchment ridge headwater tributary Some parameters η = bed elevation t = time (x1, x2) = spatial coordinate ζ = incision rate (q1, q2) = vector of volume hillslope transport/width k = hillslope diffusivity r qr = −k∇η Exner equation of sediment balance (neglecting porosity) ridge ∂η = −ς + k∇2η + v headwater ∂t u tributary Now for simplicity we assume that •All the hillslope transport occurs on the hillslope and •All the incision occurs in the channel Integrate Exner on hillslope ∂η r dA = v A − ∇ ⋅ qrdA = v A − q ds ∫∫∂t u u ∫∫u u ∫ n AAuu gain loss due to steady due to transport from state uplift hillslope to channel ∂η q = −k ~ kS ridge n ∂n h Here n is normal headwater to red perimeter of tributary path integral: Continue integration on hillslope: path integral is around red line ∂η q = −k ~ kS ridge n ∂n h qnds ~ 2Luqn ~ 2LukSh ∫ headwater tributary The following scale relation results: 2 vuAu = αuvuLu = 2αhLukSh where αh is an o(1) parameter That is, the rate of hillslope denudation must just balance with rate of uplift Incision rate: we assume ⎧vI in channel ς = ⎨ ⎩0 on hillslope Thus within the channel of the headwater tributary ∂η = −v + v or thus vI = vu ∂t I u steady Channel incision must just keep pace state with uplift Mean and fluctuating bed elevation η′ = η − η η η Video clip Hasbargen and Paola Equation of evolution of amplitude of elevation fluctuation Decompose into mean and fluctuating parts: η = η + η′ ς = ς + ς′ where the overbar denotes an average over an appropriate spatial window and the prime denotes a fluctuation about that average. Local Exner: ∂η = −ς + k∇2 η + v ∂t u ∂η 2 Local Exner: = −ς + k∇ η + v ∂t u ∂η 2 Spatially averaged Exner: = −ς + k∇ η + v ∂t u ∂ ⎛ 1 2 ⎞ 2 Multiply by η ⎜ η ⎟ = −ςη + kη∇ η + vu η = ∂t 2 and reduce: ⎝ ⎠ r r r r − ςη + k∇ ⋅ ()η∇η − k()∇η ⋅ ()∇ ⋅ η + vu η A ∂ ⎛ 1 1 ⎞ Multiply local ⎜ η2 + η′2 ⎟ = −ςη − ς′η′ Exner by η, ∂t ⎝ 2 2 ⎠ average and r r r r + k∇ ⋅ ()η∇η − k()∇η ⋅ ()∇ ⋅ η + vu η reduce: r r r r + k∇ ⋅ η′∇η′ − k ∇η′ ⋅ ∇η′ B ()()() Subtract B from A to get equation of evolution of amplitude of elevation fluctuation: ∂ ⎛ 1 ⎞ r r r r ⎜ η′2 ⎟ = −ς′η′ + k∇ ⋅ (η′∇η′)− k(∇η′)⋅ (∇η′) ∂t ⎝ 2 ⎠ IIIIII IV I. Time rate of change of amplitude of elevation fluctuation II. Rate of creation of amplitude fluctuation due to incision III. Rate of transport of amplitude fluctuation due to diffusion IV. Rate of destruction of amplitude fluctuation by hillslope diffusion Steady state: approximate balance between creation and destruction of amplitude fluctuation r r − ς′η′ ≅ k(∇η′)(⋅ ∇η′ ) Now if most of the destruction is biased toward the finest scale of the drainage basin, i.e. the headwater catchment, then r r 2 k(∇η′)⋅ (∇η′) = αdkSh where αd is an o(1) constant Approximate form for the rate of creation of amplitude fluctuation by incision ridge ηp vI∆t headwater ηc(t) ηc(t+∆t) tributary In the headwater catchment, appoximate the channel as incising into a plain with constant elevation ηp. Channel elevation ηc decreases with speed vI. η A + η BL Thus: η = p u c u where B denotes channel width Au + BLu Now ⎧ 1 2 (ηp − η) on plain 1 2 ⎪ ()η′ = 2 ⎨1 2 ⎪ (η − η)2 in channel ⎩2 c Since η& c = −vI = −vu it follows that ∂ 1 2 ∂ 1 2 ()η′ A + ()η′ BL ∂t 2 u ∂t 2 u ∂ 1 2 − η′ς′ = ()η′ = plain channel ≅ ∂t 2 A + BL incision u u BLu (η − ηc )vu Au Now scaling η − ηc = α cH where H is an "effective" flow depth and αc is an o(1) constant. Thus the rate of creation of fluctuation amplitude by incision is expressed as BLu − η′ς′ = αcHvu Au Balancing this against the rate of destruction of fluctuation amplitude, BLu 2 αcHvu = αdkSh Au Scale relations for headwater catchment: 2 Geometric scaling: Au = αuLu Rate of hillslope denudation must α v L2 = 2α L kS balance with rate of uplift: u u u h u h Rates of creation and destruction of BL α Hv u = α kS2 elevation fluctuation amplitude must c u d h Au balance: FROM THESE WE OBTAIN 1/ 3 2 / 3 L ⎛ k ⎞ ⎛ v A1/ 2 ⎞ u = α ⎜ ⎟ S = α ⎜ u e ⎟ 1/ 2 1⎜ 1/ 2 ⎟ h 2 ⎜ ⎟ A e ⎝ vuA e ⎠ ⎝ k ⎠ 1/ 3 1/ 3 1 ⎛ α ⎞ 1 ⎛ α ⎞ where A = BH = effective ⎜ c 2 ⎟ ⎜ c ⎟ e α1 = ⎜4 αh ⎟ α2 = ⎜4 ⎟ channel area and αu ⎝ αd ⎠ 2 ⎝ αdαh ⎠ Let AT = the total area of the drainage basin and LT = 1/2 (AT) denote a length scale for the total basin. It then follows that 1/ 3 1/ 3 L 1 ⎛ α ⎞ ⎛ k ⎞ ⎛ A1/ 2 ⎞ u = ⎜4 c α2 ⎟ ⎜ ⎟ ⎜ e ⎟ ⎜ h ⎟ ⎜ 1/ 2 ⎟ ⎜ ⎟ LT αu ⎝ αd ⎠ ⎝ vuA e ⎠ ⎝ LT ⎠ Now suppose that Ae is set. Then the ratio Lu/LT becomes smaller (the drainage basin becomes more intricate) as a) the size of the basin L increases, b) the rate of uplift vu increases and c) the hillslope diffusivity k decreases. In addition, Sh becomes larger (headwater hillslopes become steeper) as vu increases and k decreases. Sample calculation (numbers for k, vu courtesy Bill Dietrich, Oregon Pacific Coast Range) "Kolmogoroff" Scaling for Drainage Basins α1 10.86 αc 5 channel incision parameter α2 1.357 αd 0.5 amplitude dissipation parameter B5m αh 2 hillslope diffusion parameter H2m αu 0.5 geometric parameter for headwater basin vu 0.1 cm/yr k300cm^2/yr α1 10.9 α2 1.36 1/ 3 2 / 3 Ae 3.162 1/ 2 Lu ⎛ k ⎞ ⎛ vuA e ⎞ Rek 0.105 m = α ⎜ ⎟ S = α ⎜ ⎟ 1/ 2 1⎜ 1/ 2 ⎟ h 2 ⎜ ⎟ A e ⎝ vuA e ⎠ ⎝ k ⎠ 1/ 3 1/ 3 Lu 72.68 m 1 ⎛ α ⎞ 1 ⎛ α ⎞ α = ⎜4 c α2 ⎟ α = ⎜4 c ⎟ Sh 0.303 1 ⎜ h ⎟ 2 ⎜ ⎟ αu ⎝ αd ⎠ 2 ⎝ αdαh ⎠ θ 16.85 deg Sample calculation (numbers for k, vu courtesy Bill Dietrich, Oregon Pacific Coast Range) "Kolmogoroff" Scaling for Drainage Basins α1 10.86 αc 5 channel incision parameter α2 1.357 αd 0.5 amplitude dissipation parameter B5m αh 2 hillslope diffusion parameter H2m αu 0.5 geometric parameter for headwater basin vu 1cm/yr k300cm^2/yr α1 10.9 α2 1.36 1/ 3 2 / 3 Ae 3.162 1/ 2 Lu ⎛ k ⎞ ⎛ vuA e ⎞ Rek 1.054 m = α ⎜ ⎟ S = α ⎜ ⎟ 1/ 2 1⎜ 1/ 2 ⎟ h 2 ⎜ ⎟ A e ⎝ vuA e ⎠ ⎝ k ⎠ 1/ 3 1/ 3 Lu 33.74 m 1 ⎛ α ⎞ 1 ⎛ α ⎞ α = ⎜4 c α2 ⎟ α = ⎜4 c ⎟ Sh 1.406 1 ⎜ h ⎟ 2 ⎜ ⎟ αu ⎝ αd ⎠ 2 ⎝ αdαh ⎠ θ 54.57 deg THE END!! Or maybe not! But wait! It's not over yet! Can we explain how submarine drainage basins form? Trinity Canyon and associated drainage network, Eel Margin, Northern California And what about the fine scale of tidal drainage networks? Barnstable Salt Marsh, Cape Cod, Massachusetts Image courtesy Tao Sun, Sergio Fagherazzi and David Furbish Sample calculation Consider two rivers: a "prototype" with mean velocity U = 4 m/s and depth H = 2 m, and a "model" with mean velocity U = 1 m/s and depth H = 0.125 m.
Load more

Recommended publications
  • Palu Colostate 0053A 15347.Pdf
    DISSERTATION FLOODWAVE AND SEDIMENT TRANSPORT ASSESSMENT ALONG THE DOCE RIVER AFTER THE FUNDÃO TAILINGS DAM COLLAPSE (BRAZIL) Submitted by Marcos Cristiano Palu Department of Civil and Environmental Engineering In partial fulfillment of the requirements For the Degree of Doctor of Philosophy Colorado State University Fort Collins, Colorado Spring 2019 Doctoral Committee: Advisor: Pierre Julien Christopher Thornton Robert Ettema Sara Rathburn Copyright by Marcos Cristiano Palu 2019 All Rights Reserved ABSTRACT FLOODWAVE AND SEDIMENT TRANSPORT ASSESSMENT ALONG THE DOCE RIVER AFTER THE FUNDÃO TAILINGS DAM COLLAPSE (BRAZIL) The collapse of the Fundão Tailings Dam in November 2015 spilled 32 Mm3 of mine waste, causing a substantial socio-economic and environmental damage within the Doce River basin in Brazil. Approximately 90% of the spilled volume deposited over 118 km downstream of Fundão Dam on floodplains. Nevertheless, high concentration of suspended sediment (≈ 400,000 mg/l) reached the Doce River, where the floodwave and sediment wave traveled at different velocities over 550 km to the Atlantic Ocean. The one-dimensional advection- dispersion equation with sediment settling was solved to determine, for tailing sediment, the longitudinal dispersion coefficient and the settling rate along the river and in the reservoirs (Baguari, Aimorés and Mascarenhas). The values found for the longitudinal dispersion coefficient ranged from 30 to 120 m2/s, which are consistent with those in the literature. Moreover, the sediment settling rate along the whole extension of the river corresponds to the deposition of finer material stored in Fundão Dam, which particle size ranged from 1.1 to 2 . The simulation of the flashy hydrographs on the Doce River after the dam collapse was initially carried out with several widespread one-dimensional flood routing methods, including the Modified Puls, Muskingum-Cunge, Preissmann, Crank Nicolson and QUICKEST.
    [Show full text]
  • Fluvial Sedimentary Patterns
    ANRV400-FL42-03 ARI 13 November 2009 11:49 Fluvial Sedimentary Patterns G. Seminara Department of Civil, Environmental, and Architectural Engineering, University of Genova, 16145 Genova, Italy; email: [email protected] Annu. Rev. Fluid Mech. 2010. 42:43–66 Key Words First published online as a Review in Advance on sediment transport, morphodynamics, stability, meander, dunes, bars August 17, 2009 The Annual Review of Fluid Mechanics is online at Abstract fluid.annualreviews.org Geomorphology is concerned with the shaping of Earth’s surface. A major by University of California - Berkeley on 02/08/12. For personal use only. This article’s doi: contributing mechanism is the interaction of natural fluids with the erodible 10.1146/annurev-fluid-121108-145612 Annu. Rev. Fluid Mech. 2010.42:43-66. Downloaded from www.annualreviews.org surface of Earth, which is ultimately responsible for the variety of sedi- Copyright c 2010 by Annual Reviews. mentary patterns observed in rivers, estuaries, coasts, deserts, and the deep All rights reserved submarine environment. This review focuses on fluvial patterns, both free 0066-4189/10/0115-0043$20.00 and forced. Free patterns arise spontaneously from instabilities of the liquid- solid interface in the form of interfacial waves affecting either bed elevation or channel alignment: Their peculiar feature is that they express instabilities of the boundary itself rather than flow instabilities capable of destabilizing the boundary. Forced patterns arise from external hydrologic forcing affect- ing the boundary conditions of the system. After reviewing the formulation of the problem of morphodynamics, which turns out to have the nature of a free boundary problem, I discuss systematically the hierarchy of patterns observed in river basins at different scales.
    [Show full text]
  • Manuscript, Helped to Motivate the Work
    Long-Profile Evolution of Transport-Limited Gravel-Bed Rivers Andrew D. Wickert1 and Taylor F. Schildgen2,3 1Department of Earth Sciences and Saint Anthony Falls Laboratory, University of Minnesota, Minneapolis, Minnesota, USA. 2Institut für Erd- und Umweltwissenschaften, Universität Potsdam, 14476 Potsdam, Germany. 3Helmholtz Zentrum Potsdam, GeoForschungsZentrum (GFZ) Potsdam, 14473 Potsdam, Germany. Correspondence to: A. D. Wickert ([email protected]) Abstract. Alluvial and transport-limited bedrock rivers constitute the majority of fluvial systems on Earth. Their long profiles hold clues to their present state and past evolution. We currently possess first-principles-based governing equations for flow, sediment transport, and channel morphodynamics in these systems, which we lack for detachment-limited bedrock rivers. Here we formally couple these equations for transport-limited gravel-bed river long-profile evolution. The result is a new predictive 5 relationship whose functional form and parameters are grounded in theory and defined through experimental data. From this, we produce a power-law analytical solution and a finite-difference numerical solution to long-profile evolution. Steady-state channel concavity and steepness are diagnostic of external drivers: concavity decreases with increasing uplift, and steepness increases with increasing sediment-to-water supply ratio. Constraining free parameters explains common observations of river form: To match observed channel concavities, gravel-sized sediments must weather and fine – typically rapidly – and valleys 10 must widen gradually. To match the empirical square-root width–discharge scaling in equilibrium-width gravel-bed rivers, downstream fining must occur. The ability to assign a cause to such observations is the direct result of a deductive approach to developing equations for landscape evolution.
    [Show full text]
  • Introduction to Morphodynamics of Sedimentary Patterns
    MORPHODYNAMICS OF SEDIMENTARY PATTERNS Paolo Blondeaux, Marco Colombini, Giovanni Seminara, Giovanna Vittori Introduction to Morphodynamics of Sedimentary Patterns DIDATTICARICERCA RICERCA Genova University Press Monograph Series Morphodynamics of Sedimentary Patterns Editorial Board: Paolo Blondeaux, Marco Colombini Giovanni Seminara, Giovanna Vittori Advisory Board: Sivaramakrishnan Balachandar (University of Florida U.S.A.) Maurizio Brocchini (Università Politecnica delle Marche, Italy) François Charru (Université Paul Sabatier, France) Giovanni Coco (University of Auckland, New Zealand) Enrico Foti (Università di Catania, Italy) Marcelo H. Garcia (University of Illinois, U.S.A.) Suzanne J.M.H. Hulscher (University of Twente, NL) Stefano Lanzoni (Università di Padova, Italy) Miguel A. Losada (University of Granada, Spain) Chris Paola (University of Minnesota, U.S.A.) Gary Parker (University of Illinois U.S.A.) Luca Ridolfi(Politecnico di Torino, Italy) Andrea Rinaldo (École Polytechnique Fédérale de Lausanne, Switzerland) Yasuyuki Shimizu (Hokkaido University, Japan) Marco Tubino (Università di Trento, Italy) Markus Uhlmann (Karlsruhe Institute of Technology, Germany) Introduction to Morphodynamics of Sedimentary Patterns Paolo Blondeaux, Marco Colombini, Giovanni Seminara, Giovanna Vittori Introduction to Morphodynamics of Sedimentary Patterns is the bookmark of the University of Genoa On the cover: The delta of Lena river (Russia). The image was taken on July 27, 2000 by the Landsat 7 satellite operated by the U.S. Geological Survey and NASA (false-color composite image made using shortwave infrared, infrared, and red wavelengths). Image credit: NASA This book has been object of a double peer-review according with UPI rules. Publisher GENOVA UNIVERSITY PRESS Piazza della Nunziata, 6 16124 Genova Tel. 010 20951558 Fax 010 20951552 e-mail: [email protected] e-mail: [email protected] http://gup.unige.it/ The authors are at disposal for any eventual rights about published images.
    [Show full text]
  • Macro-Roughness Model of Bedrock–Alluvial River Morphodynamics
    Earth Surf. Dynam., 3, 113–138, 2015 www.earth-surf-dynam.net/3/113/2015/ doi:10.5194/esurf-3-113-2015 © Author(s) 2015. CC Attribution 3.0 License. Macro-roughness model of bedrock–alluvial river morphodynamics L. Zhang1, G. Parker2, C. P. Stark3, T. Inoue4, E. Viparelli5, X. Fu1, and N. Izumi6 1State Key Laboratory of Hydroscience and Engineering, Tsinghua University, Beijing, China 2Department of Civil & Environmental Engineering and Department of Geology, Hydrosystems Laboratory, University of Illinois, Urbana, IL, USA 3Lamont-Doherty Earth Observatory, Columbia University, Palisades, NY, USA 4Civil Engineering Research Institute for Cold Regions, Hiragishi Sapporo, Japan 5Dept. of Civil & Environmental Engineering, University of South Carolina, Columbia, SC, USA 6Faculty of Engineering, Hokkaido University, Sapporo, Japan Correspondence to: L. Zhang ([email protected]) Received: 14 April 2014 – Published in Earth Surf. Dynam. Discuss.: 27 May 2014 Revised: 6 January 2015 – Accepted: 6 January 2015 – Published: 11 February 2015 Abstract. The 1-D saltation–abrasion model of channel bedrock incision of Sklar and Dietrich (2004), in which the erosion rate is buffered by the surface area fraction of bedrock covered by alluvium, was a major advance over models that treat river erosion as a function of bed slope and drainage area. Their model is, however, limited because it calculates bed cover in terms of bedload sediment supply rather than local bedload transport. It implicitly assumes that as sediment supply from upstream changes, the transport rate adjusts instantaneously everywhere downstream to match. This assumption is not valid in general, and thus can give rise to unphysical consequences.
    [Show full text]
  • Comparing 2-D and 1-D/2-D Modelling of Agly River Bed Change During a Flood S
    Comparing 2-D and 1-D/2-D modelling of Agly River bed change during a flood S. Mezbache, André Paquier, M. Hasbaia To cite this version: S. Mezbache, André Paquier, M. Hasbaia. Comparing 2-D and 1-D/2-D modelling of Agly River bed change during a flood. River Flow 2020, Jul 2020, Delft, Netherlands. pp.704-710, 10.1201/b22619. hal-03107142 HAL Id: hal-03107142 https://hal.archives-ouvertes.fr/hal-03107142 Submitted on 12 Jan 2021 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Comparing 2-D and 1-D/2-D modelling of Agly River bed change during a flood S. Mezbache M’sila University, M’sila, Algeria A. Paquier INRAE, U.R. RiverLy, Villeurbanne, France M. Hasbaia M’sila University, M’sila, Algeria ABSTRACT: The authors discuss the use of a 1-D model or a 2-D model for the evolution of the bed between the levees while the spreading of water and sediment over the floodplains is simu- lated using a 2-D model. A 1-D model cannot represent the asymmetrical behaviour of the flow in the bends but makes easier the estimate of the variations of the sediment transport capacity along the river if the bed geometry is complex.
    [Show full text]
  • Three Dimensional Mobile Bed Dynamics for Sediment Transport Modeling
    THREE DIMENSIONAL MOBILE BED DYNAMICS FOR SEDIMENT TRANSPORT MODELING DISSERTATION Presented in Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy in the Graduate School of The Ohio State University By Sean O’Neil, B.S., M.S. ***** The Ohio State University 2002 Dissertation Committee: Approved by Professor Keith W. Bedford, Adviser Professor Carolyn J. Merry Adviser Professor Diane L. Foster Civil Engineering Graduate Program c Copyright by Sean O’Neil 2002 ABSTRACT The transport and fate of suspended sediments continues to be critical to the understand- ing of environmental water quality issues within surface waters. Many contaminants of environmental concern within marine and freshwater systems are hydrophobic, thus read- ily adsorbed to bed material or suspended particles. Additionally, management strategies for evaluating and remediating the effects of dredging operations or marine construction, as well as legacy pollution from military and industrial processes requires knowledge of sediment-water interactions. The dynamic properties within the bed, the bed-water column inter-exchange and the transport properties of the flowing water is a multi-scale nonlin- ear problem for which the mobile bed dynamics with consolidation (MBDC) model was formulated. A new continuum-based consolidation model for a saturated sediment bed has been developed and verified on a stand-alone basis. The model solves the one-dimensional, vertical, nonlinear Gibson equation describing finite-strain, primary consolidation for satu- rated fine sediments. The consolidation problem is a moving boundary value problem, and has been coupled with a mobile bed model that solves for bed level variations and grain size fraction(s) in time within a thin layer at the bed surface.
    [Show full text]
  • Autogenic Erosional Surfaces in Fluvio-Deltaic Stratigraphy from Floods, Avulsions, and Backwater Hydrodynamics
    Journal of Sedimentary Research, 2019, v. 89, 815–832 Research Article DOI: http://dx.doi.org/10.2110/jsr.2019.40 AUTOGENIC EROSIONAL SURFACES IN FLUVIO-DELTAIC STRATIGRAPHY FROM FLOODS, AVULSIONS, AND BACKWATER HYDRODYNAMICS 1,2 3 3 VAMSI GANTI, MICHAEL P. LAMB, AND AUSTIN J. CHADWICK 1Department of Geography, University of California Santa Barbara, Santa Barbara, California, U.S.A. 2Department of Earth Science, University of California Santa Barbara, Santa Barbara, California, U.S.A. 3Division of Geological and Planetary Sciences, California Institute of Technology, Pasadena, California, U.S.A. e-mail: [email protected] ABSTRACT: Erosional surfaces set the architecture of fluvio-deltaic stratigraphy, and they have classically been interpreted in terms of changes in boundary conditions such as climate, tectonics, and base level (allogenic forces). Intrinsic dynamics of sedimentary systems (autogenic dynamics) can also create a rich stratigraphic architecture, and a major knowledge gap exists in parsing the relative roles of autogenic versus allogenic processes. Emerging theoretical and experimental work suggests that backwater hydrodynamics play an important role in driving transient channel incision in river deltas, even those experiencing net aggradation. Here, we identify and quantify two autogenic mechanisms that produce broad erosional surfaces in fluvio-deltaic stratigraphy, namely, floods and avulsions. Using a simple mass-balance model for single-threaded delta channel systems, we show that flood-induced scours begin near the shoreline, and avulsion-induced scours begin at the avulsion site, and both propagate upstream over a distance that scales with the backwater length, bed slope, and bed grain size. We also develop scaling relationships for the maximum scour depths arising from these mechanisms, which are functions of characteristic flow depth and formative flood variability.
    [Show full text]
  • Advanced Numerical Modeling of Sediment Transport in Gravel-Bed Rivers
    water Article Advanced Numerical Modeling of Sediment Transport in Gravel-Bed Rivers Van Hieu Bui 1,2,* , Minh Duc Bui 1 and Peter Rutschmann 1 1 Institute of Hydraulic and Water Resources Engineering, Technische Universität München, Arcisstrasse 21, D-80333 München, Germany; [email protected] (M.D.B.); [email protected] (P.R.) 2 Faculty of Mechanical Engineering, Thuyloi University, 175 Tay Son, Dong Da, Hanoi 100000, Vietnam * Correspondence: [email protected]; Tel.: +49-152-134-15987 Received: 15 February 2019; Accepted: 13 March 2019; Published: 17 March 2019 Abstract: Understanding the alterations of gravel bed structures, sediment transport, and the effects on aquatic habitat play an essential role in eco-hydraulic and sediment transport management. In recent years, the evaluation of changes of void in bed materials has attracted more concern. However, analyzing the morphological changes and grain size distribution that are associated with the porosity variations in gravel-bed rivers are still challenging. This study develops a new model using a multi-layer’s concept to simulate morphological changes and grain size distribution, taking into account the porosity variabilities in a gravel-bed river based on the mass conservation for each size fraction and the exchange of fine sediments between the surface and subsurface layers. The Discrete Element Method (DEM) is applied to model infiltration processes and to confirm the effects of the relative size of fine sediment to gravel on the infiltration depth. Further, the exchange rate and the bed porosity are estimated while using empirical formulae. The new model was tested on three straight channels.
    [Show full text]
  • Chapter 4 Alternate Bars and Alluvial Fans 4.1
    CHAPTER 4 ALTERNATE BARS AND ALLUVIAL FANS 4.1 INTRODUCTION At least as interesting as the mechanism of sediment transport itself is the ability of rivers to construct morphologies. Primary among these morphologies is the river itself, with a definable self-formed channel maintaining consistent bankfull geometry, as well as a floodplain. Within the river there can be bedforms such as ripples, dunes and antidunes, and bar forms such as alternate bars and the point bars of bends. The river can meander or braid, and is the agent of construction of larger morphologies such as drainage basins and alluvial fans. In the interest of space, only alternate bars and alluvial fans are treated here. 4.2 ALTERNATE BARS In streams of sufficient width B the bed of the channel does not remain plane, but rather is subject to bar instability. The governing parameter here is the ratio of width to depth B/H. When this number exceeds a value between 12 and 20, the flow becomes unstable and creates the single-row alternate bars shown in Figure 4.1 below. Single-row alternate bars can be considered to be precursors to meandering. If the banks are inerodible, the bars grow to a finite-amplitude equilibrium and then migrate downstream. If the banks are erodible, the bars induce enough sinuosity to allow the primary mechanism for meandering to take hold. At larger values of B/H, multiple-row alternate bars can form. These are precursors to braiding. GP Lecture Notes: Sediment Transport 30 Figure 4.1 Alternate bars. 4.3 FORMULATION FOR ALTERNATE BARS As in most sediment transport problems, the quasi-steady approximation can be employed.
    [Show full text]
  • Morphodynamics of River Bed Variation with Variable Bedload Step Length
    Earth Surf. Dynam., 2, 243–253, 2014 Open Access www.earth-surf-dynam.net/2/243/2014/ Earth Surface doi:10.5194/esurf-2-243-2014 © Author(s) 2014. CC Attribution 3.0 License. Dynamics Morphodynamics of river bed variation with variable bedload step length A. Pelosi1 and G. Parker2 1Department of Civil Engineering, Università degli Studi di Salerno (UNISA), Via Giovanni Paolo II – 84084 Fisciano (Salerno), Italy 2Department of Civil & Environmental Engineering and Department of Geology, Hydrosystems Laboratory, University of Illinois, 301 N. Mathews Ave., Urbana, IL 61801, USA Correspondence to: A. Pelosi ([email protected]) Received: 25 October 2013 – Published in Earth Surf. Dynam. Discuss.: 22 November 2013 Revised: 24 March 2014 – Accepted: 7 April 2014 – Published: 23 April 2014 Abstract. Here we consider the 1-D morphodynamics of an erodible bed subject to bedload transport. Fluvial bed elevation variation is typically modeled by the Exner equation, which, in its classical form, expresses mass conservation in terms of the divergence of the bedload sediment flux. An entrainment form of the Exner equa- tion can be written as an alternative description of the same bedload processes, by introducing the notions of an entrainment rate into bedload and of a particle step length, and assuming a certain probability distribution for the step length. This entrainment form implies some degree of nonlocality, which is absent from the standard flux form, so that these two expressions, which are different ways to look at same conservation principle (i.e., sediment continuity), may no longer become equivalent in cases when channel complexity and flow conditions allow for long particle saltation steps (including, but not limited to the case where particle step length has a heavy tailed distribution) or when the domain of interest is not long compared to the step length (e.g., labora- tory scales, or saltation over relatively smooth surfaces).
    [Show full text]
  • A Two Dimensional Finite Element Hydrodynamic River Morphology and Gravel Transport Model
    2nd Joint Federal Interagency Conference, Las Vegas, NV, June 27 - July 1, 2010 A TWO DIMENSIONAL FINITE ELEMENT HYDRODYNAMIC RIVER MORPHOLOGY AND GRAVEL TRANSPORT MODEL Stephen Kwan, Hydraulic Engineering Concepts, Bellevue, WA, [email protected]; Jose A. Vasquez, Northwest Hydraulic Consultants, Vancouver, BC, Canada, [email protected]; Robert G. Millar, University of British Columbia, Vancouver, BC, Canada, [email protected]; Peter M. Steffler, University of Alberta, Edmonton, AB, Canada, [email protected]. Abstract A depth-averaged finite element two-dimensional hydrodynamic-morphological and gravel transport model is presented. Based on River2D by Steffler et al. (2002) of the University of Alberta, River2D-Morphology (R2DM) uses an unstructured mesh to simulate flow hydraulics, sediment transport for uni-size and mixed size sediment, and morphological changes of a river over time. Changes in bed elevation are calculated by solving Exner’s equation for sediment mass conservation. For uni-size sediment transport, R2DM uses transport functions of Meyer-Peter-Muller, Engelund-Hansen, Van Rijn or an empirical formula. For mixed size sediment transport, the function of Wilcock and Crowe (2003) is used. Secondary flow effects are modeled using the technique of Struiksma et al. (1985) with a unique algorithm that calculates the radius of curvature from the streamlines. R2DM is capable of exploring the dynamics of grain size distributions such as vertical sorting and armoring and can calculate the fraction of sand on the bed surfaces of rivers and streams. Results have been verified with experimental data for aggradation and downstream fining (Paolo, 1992; Seal, 1992; Toro-Escobar, 2000); for degradation (Ashida and Michiue, 1971); and for dispersion (Cui and Parker, 2003).
    [Show full text]