avances en la modelacion de avenidas subitas, transporte ...pierre/ce_old... · avances en la...

1

Pierre Y. JulienColorado State University

Jiutepec, Mexico

November 20, 2008

Avances en la modelacion de avenidas subitas, transporte

de sedimentos y contaminantes con TREX

PROBLEM: Upland ErosionIntroduct ion

Erosion Picture

2

PROBLEM: … and DepositionDeposition pictureI

ntroduct ion

CASC2D- Julien et al. (1995)CASC2D-SED – Johnson et al. (2000), Rojas (2002)

Rainfall

RetentionInfiltration

CASC2D-SED

3

Goodwin Creek WatershedApplication

• Location: Panola County (MS)

• Area: 21Km2

• Monitored by ARS-NSL (Oxford, MS)

• 37 rain gages

• 14 stream gages (water and sediment)

• Channel surveys

• GIS data

INPUT DATA (DEM)

Digital Elevation Model

Derive:1. Channel Network2. Slope distribution

Application

4

INPUT DATA (soil type)Application

INPUT DATA (land use)Application

5

INPUT DATA (rainfall)

10/17/81 event:

Duration: 3.5 hr.

Depth: 73 mm. 0.0

0.5

1.0

1.5

2.0

2.5

3.0

3.5

4.0

0 50 100 150 200 250 300 350Time (min)

Inte

nsity

(in/

h)

1 2

5 6

7 10

12 14

34 35

41 42

52 61

64 65

Rain Gage Number

Raingage location

Application

Receiving Cell

Receiving Cell

Parent material

Deposition

Suspension

Deposition

Suspension

qsy

qsxOutgoing Cell

SEDIMENT ROUTING

Availablematerial

CASC2D-SED

AdvectionCapacity

vs. supply

6

UPLAND EROSION (2-D)Modified Kilinc and Richardson equation for the overland:

PC15.0

KWQS23210)s*m/tons(q

035.266.1

ot ⎟⎠⎞

⎜⎝⎛=

CASC2D-SED

DEM

Hydraulics

SoilsLand use

GEOVISUALIZATION rainfall-runoff

10/17/81: calibration event:

Rainfall rates, infiltration, water depth and sediment concentration

Results

7

GEOVISUALIZATION rainfall-runoff

10/17/81: calibration event:

Rainfall rates, infiltration, water depth and sediment concentration

Results

012345678

0 200 400 600 800Time [min]

Run

off [

mm

/h]

MODEL CALIBRATION

02468

1012141618

0 200 400 600 800Time [min]

Run

off [

mm

/h]

02468

1012141618

0 200 400 600 800Time [min]

Run

off [

mm

/h]

0

2

4

6

8

10

12

0 200 400 600 800Time [min]

Run

off [

mm

/h]

0

2

4

6

8

10

12

14

0 200 400 600 800Time [min]

Run

off [

mm

/h]

0

2

4

6

8

10

12

14

0 200 400 600 800Time [min]

Run

off [

mm

/h]

ObservedSimulated

Hydrographs

INTERNAL VALIDATION

8

02468

101214161820

0 200 400 600 800Time [min]

Qs [

tons

/ ha

/ da

y] 02468

10121416

0 200 400 600 800Time [min]

Qs [

tons

/ ha

/ da

y]

02468

101214161820

0 200 400 600 800Time [min]

Qs [

tons

/ ha

/ da

y]

0

10

20

30

40

50

60

0 200 400 600 800Time [min]

Qs [

tons

/ ha

/ da

y]

05

101520253035404550

0 200 400 600 800Time [min]

Qs [

tons

/ ha

/ da

y]

MODEL CALIBRATION

ObservedSimulated

012345678

0 200 400 600 800Time [min]

Qs [

tons

/ ha

/ da

y]

Sedigraphs

INTERNAL VALIDATION

SEDIMENT RATING CURVES

9

NET EROSION (Uniform)

Grid

Size

Net Erosion and Deposition*

*Net difference between erosion and deposition.

10

California Gulch Watershed• EPA Superfund Site

• Location: Lake County (CO)

100-year flood: 2-h: 1.73 in

Leadville

11

Model Application to California Gulch

Elevation Range:2909 - 3654 m

Mean Slope: 12.6 %

Topography:Basin is 30.6 km2

Transport Issue:Contaminant sources distributed across site (~2000 sources)

Metals: Cd, Cu, Zn

Toxicity to fish and benthos downstream

Contaminant Characterization: Soils

12

Halgren, Julien CG/AR Research Progress Meeting June 2008

California Gulch, Colorado

24

Real-Time Precip. and FluvialReal-Time FluvialOther Airport Precip.

Data Locations:

CG/AR Progress Review: June 2008

13

Erosion and Sediment Transport

Spatially Distributed Model Inputs

Elevation (m)High : 3650

Low : 2900

Digital ElevationModel

Soils

Soil TypeWet Alluvial Land (Wa)

Gravel Pit (GP)

Rosane Loam (Rtc)

Troutville Gravelly Sandy Loam (TrE)

Perian Soi ls (PIF)

Water (Ponds) (W)

Newfork Gravelly Sandy Loam (NfB)

Perian Gravelly Sandy Loam (PgD)

Leadville Sandy Loam (LeE)

Mine Pits and Dumps (MP)

Slickens (mill taili ngs) (Sw)

Placer Diggings and Tailings (Pn)

Bross Gravelly Sandy Loam (BrF)

Tomichi Sandy Loam (ToE)

Low Densi ty Urban (Leadville Sandy Loam)

High Density Urban (Leadville Sandy Loam)

Commerical (Leadville Sandy Loam)

Land Use

Land UseOpen Water

Perennial Snow/Ice

Low Density Urban

High Density Urban

Commercial

Bare Rock/Sand

Deciduous Forest

Evergreen Forest

Mixed Forest

Shrubland

Grasslands/Herbaceous

Pasture/Hay

Row Crops

14

TREX Schematic

–

Water depths from a rainfall event*

*1-in-100 year intensity, 2 hour duration uniform rainfall event.

15

Event transport of sediment (TSS)*

*Transport is computed by grain size. Total solids shown.

CG1

0

50

100

150

200

250

300

0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08Flow (m3/s)

TSS

(mg/

L)

Simulated (June 12-13, 2003)

Simulated (Sept 5-8, 2003)

Observed (1984-2004)

Simulated vs. Observed TSS

CG4

0

200

400

600

800

1000

1200

0 0.05 0.1 0.15 0.2 0.25 0.3 0.35

Flow (m3/s)

TSS

(mg/

L)

Simulated (June 12-13, 2003)Simulated (Sept 5-8, 2003)Observed (1984-2004)

CG6

0

50

100

150

200

250

300

350

400

450

500

0 0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16Flow (m3/s)

TSS

(mg/

L)

Simulated (June 12-13, 2003)

Simulated (Sept 5-8, 2003)

Observed (1984-2004)

16

California Gulch Watershed

Total Zinc Transport: June 12-13, 2005

17

Simulated vs. Observed ZincCG1

0

5

10

15

20

25

30

35

0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08Flow (m3/s)

Tota

l Zn

(mg/

L)

Simulated (June 12-13, 2003)

Simulated (Sept 5-8, 2003)

Observed (1984-2004)

CG4

0

10

20

30

40

50

60

0.00 0.05 0.10 0.15 0.20 0.25 0.30 0.35Flow (m3/s)

Tota

l Zn

(mg/

L)

Simulated (June 12-13, 2003)

Simulated (Sept 5-8, 2003)

Observed (1984-2004)

CG6

0

10

20

30

40

50

60

0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16Flow (m3/s)

Tota

l Zn

(mg/

L)

Simulated (June 12-13, 2003)

Simulated (Sept 5-8, 2003)

Observed (1984-2004)

Total Zinc Transport: 1-in-100-yr Event

18

Zinc Accumulation: 1-in-100-yr Event

CSU Watershed Model TREX

Flow Depth [ft] 100 yr StormSurface Water Depth [ft]

19

Total rainfall of Typhoon Maemi from 9/8/2003 to 9/13/2003)

Mangun mountain, South Korea

20

Annual average: 3,450 tons/km2/year.

Soil loss mapsTyphoon Maemi: 2,920 tons/km2/year.

Imha reservoir

21

CASC2D-SED Modeling 2004

CONCLUSIONS

• TREX simulates the hydrologic response of watershed at resolutions of 30m.

• The model has been calibrated on several watersheds for water, sediment and metal fluxes.

• MPEG movies enhance the visualization of hydrologic processes including sediment and metals.

• The model is available on the web … for free!

TREX and CASC2D-SED Web Pages

At Colorado State UniversityUnder direction of Dr. Pierre Julien

pierre@engr.colostate.edu

Current manual, source code, example, MPEG movies

22

Acknowledgments

Dr. Mark Velleux, CSU, now HydroqualDr. John England, CSU, also USBRDr. Rosalia Rojas, CSUHyeon Sik Kim, CSU, also KOWACOJames Halgren, PhD candidate at CSUSeema Shah-Fairbank, PhD candidate at CSU

CASC2D-SED Modeling 2004

pierre@engr.colostate.edu

MuchasGracias!