Top Banner
Sub-hourly Modeling of Sheet Erosion and Sediment Transport using the SWAT Model José María Bodoque; Ana Lucia; José Francisco Martín Duque; Juan Antonio Ballesteros; Bouchra Haddad
18

José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Jul 09, 2020

Download

Documents

dariahiddleston
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Sub-hourly Modeling of Sheet Erosion and Sediment Transport using the SWAT Model

José María Bodoque; Ana Lucia; José Francisco Martín Duque;

Juan Antonio Ballesteros; Bouchra Haddad

Page 2: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Why ?

Gullies developed on sands:

Rare worldwide

Few studies on such gullies

Apparent high activity and hillslope-channel coupling

No previous hydrogeomorphic studies in this area

Page 3: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Objectives

Understand the geomorphic processes operating within these gullies.

Demonstrate sheet erosion is the predominant erosive process

Estimate sediment yield

Page 4: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Study Area

Continental mediterranean climate – Precipitation 650 mm/year

Cretaceous bedrock sediments

Mesas and cuestas Sandy cambisols Holm oak and savin

forests

Zona de estudio

PROVINCIA DE SEGOVIA

Segovia

Sierra

de

Guadarra

ma

Rio Duero

Rio CegaRio M

oros

Page 5: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Experimental catchment: Barranca de los Pinos

Landform map

– High gradient slopes

Gravitational processes

– Low gradient slopes

Erosive processes due to overland flow and splash

– Channels

Fluvial processes

Page 6: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Experimental catchment: Barranca de los Pinos

Page 7: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Methods

Estimate sheet erosion rates by using a dendrogeomorphic aproach

Page 8: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Methods

Estimate sheet erosion rates by using a dendrogeomorphic aproach

CES: Sediment delivery ratio

A: Catchment area

P: Main channel slope

BR: Bifurcation ratio (f – stream order)

Page 9: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Field measurement of total soil erosion and runoff

Methods

Page 10: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Hudson, 1997

Methods

Parshall flume water discharge

Field measurement of total soil erosion and runoff

Page 11: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Siphon samplers

suspended sediments and solutes

Hudson, 1997

Methods

6 ½-liter bottles

Field measurement of total soil erosion and runoff

Page 12: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Methods

Reid slot sampler automatic and

continuous bedload sampling

Field measurement of total soil erosion and runoff

Page 13: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

SolutesSuspended SedimentBedload

Methods Channel: Basin outlet

Total sediment yield : 3.3 Mg

Effective rain: 9.2 mm Runoff coefficient 5.4 %

0

20

40

60

80

100

120

140

0.0

1.0

2.0

3.0

4.0

5.0

6.0

7.0

8.0

0 10 20 30 40 50 60

Cum

ula

tive

be

dlo

ad

ma

ss (

kg)

Wate

r d

ep

th (

cm

) W

ate

r d

ischa

rge

(l/s)

Be

dlo

ad

flu

x (

kg/s

m)

Time (min)

water depth

bedload flux

water discharge

bedload mass

Page 14: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Methods Subhourly SWAT model -Motivation

The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended for daily upland erosion and sediment transport modelling in overland flow

The MUSLE is an empirical model developed for predicting long term average soil loss and is NOT adequate for subdaily continous simulation

Subdaily erosion and sediment transport is not available in SWAT 2005/2009

Jeong et al., 2010

Page 15: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Methods Subhourly SWAT model -Motivation

1

)()(

tFKtf

Where K is the hydraulic conductivity, y is the wetting front suction head, is the change in moisture content and F(t) is the cumulative infiltration depth.

Estimation of infiltration and excess rainfall

Page 16: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Methods Subhourly SWAT model -Motivation

Surface runoff lag

Page 17: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Methods Subhourly SWAT model -Motivation

Unit hydrograph

tp

Tr/2

Tr

tb

Tp

Lluvia neta

Escorrentía

directa

qp

Page 18: José María Bodoque; Ana Lucia; José Francisco Martín Duque ... · Methods Subhourly SWAT model -Motivation The Universal Soil Loss Equation (MUSLE) in SWAT 2005/2009 is intended

Methods Subhourly SWAT model -Motivation

Channel impoundment routing

Advancing

Flood

Wave

I > Q

IQ

QQ

QI

II

IQ

I Q

Receding

Flood

Wave

Q > I

KQS Prism

)(Wedge QIKXS

K = travel time of peak through the reach

X = weight on inflow versus outflow (0 ≤ X ≤ 0.5)

X = 0 Reservoir, storage depends on outflow, no

wedge

X = 0.0 - 0.3 Natural stream

)( QIKXKQS

])1([ QXXIKS