Hydrological Applications on
AWCI Reference Basins

Oliver SAAVEDRA
Earth Observation Data Integration and Fusion Research Initiative(EDITORIA)
Department of Civil Engineering

Update on AWCI Reference basins

Ciliwung..to last one

Upper Tone River Basin
Yagisawa
Fujiwara

Shimagawa

Naramata

AimataSonohara

Yamba

Multipurpose Multi-Reservoir
Error Forecast

Short-term
Weather forecast Output

Spatially distributed data

Observed rainfall
Optimization of existing dam network

Coupled DHM & Dam Operation

Adjusted
OperationalRule

Efficient water resourcesmanagement

Optimization
schemeOptimal dam water
release scheduleSpatial input data

Weather forecast data
Grid Point Value (GPV)- Meso-scale (2002/07~)- 0.125º spatial resolution- issued every 6 hours-18 hours lead-time (hourly)

000612

1
8000
612Day 1

at 00

at 06

1
800

Day
2

at 12

at 18

Flood reduction with GPV 7~12

Optimized release
rules
Optimized

Flood peakreduction

Outflow eq. 0
2000

560

800

550

700

540

600

Sim outflow
Sim inflow
Sim water levelObs water level500
400300

530
520
510
500

200

490

100

480

0
7/8.1z470
7/9.1z7/10.1z

1500

7/11.1z

7/12.1z

2002

1000

Water level increase due to storage

500

Fujiwara dam

0

7/11.1z

7/12.1z

2002

Peak created due towater release from dams..awci_basins_map

37/10.1zdischarge [m /s]

7/9.1z

1000
900
800
700
600
500
400
300
200
100
0
7/8.1z7/9.1z7/10.1z

7/11.1z
2002

650
645
640
635
630
Sim outflow
625
Sim inflow
620
Sim water level 615
Obs water level 610
605
600
7/12.1z

waterlevel [m]3

Outflow eq. inflow

7/8.1z570

Sonohara dam

3

3000

discharge [m /s]

900
discharge [m /s]

Iwamoto gauge

25001000waterlevel [m]Water is stored until
max capacity is reachedApplication on Pampanga River••••

Simulated area: 10,000 km2
Computing Grid: 1800 m
Target event:Oct-Dec 2004
Input Rainfall:– Observed• TRRM, 3hr, 0.25°– Forecast• Global JMA/GPV,
24 h lead every 12 h

Angat damAngat DAM

Bypass #2
Auxilliary Units
Bypass #1

POWER TUNNELPENSTOCKSan Jose Substation

IPO DAM

LOW LEVEL OUTLET

BUSTOS DAM

Click here for:

Main Units

AHEPP OPERATING
CONCEPTANGAT HYDROELECTRIC POWER PLANT

Customize of Cost Function
Not only to minimize flood damage
but also to consider future water use and power generation

Flood Control
Angat River Basin has a long
history of floods.
serving as a flood control facilityof about 63MCM of floodwaters

Water Use-To supply the domestic and industrial
water requirements of residents in
Metro Manila (97%)-To provide irrigation to about 31,000
hectares of land in 20 municipalities andtowns in Pampanga and Bulacan

Maximize {the benefit fromwater use and power generation}

Power Generation

Minimize {the risk of flooddamage}To generate hydroelectric power tofeed the Luzon Grid

TRMM rainfall distribution

Simulated Hourly Q at Manila Bay
30000
2004.11.29

25000
20000
15000
2004.12.02

10000
5000
0

GLOBAL DATASETSRTM90 m DEM

FAO
Soil Map

DERIVEDINPUT DATA & PARAMETERS•Elevation & slope•River channel network•Geomorphology properties•Sub- basins delineation•Soil unit
•Sat. Hydraulic Conductivity•Sat. Soil Moisture•Alpha & water N

LOCALDATA SET

50 m DEM

DigitalSoil Map

Satellite
Land use•Max. canopy storage•Max. surface storage•Vegetation coverage

DigitalLand use

NDVI orLAIVegetation indexes
e.g. Max. Leaf Area Index

Field work
Inspection

Satellite
RainfallSpatial distributed rainfall

Observed rainfall•Rain gauge•RadarSatellite
Imagery•River morphology•River cross- sections

Measured Riverchannels

Globaldischarge

River stream flow

Observed
River discharge

Satellite
Imagery..awci_basins_map

Estimation of overflow

Flood Inundationmaps

Huong River, Vietnam

Optimal flood warning
Error Forecast

Short-term
Weather forecast Output

Spatially distributed data

Observed rainfall

Distributed Hydrological Model
Adjusted
Flood warning
Level

Efficientflood warning

Optimizationscheme

Optimal flood warning

Facts about Huong Model• Working window:• Simulated area:• Computing Grid:• Target event:• Input Rainfall:

107.5-107.75E 16.0-16.5 N
1500 km2 down to Kim Long
500 m ( from 50 m DEM)
22-26 Nov 2004

– Observed• Rain gauge network (daily)• Global Satellite TRRM, 3hr, 0.25°– Forecast• Meso- scale HRM: at hydro-met. stations (24 hr lead time)• Global JMA/GPV, 24 hr lead time issued every 12 hours• NWPO: UKMO, NCEP

– Assimilated• Clouds microphysics ARPS, IMDASDischarge [m3/s]

24-30 November Simulation
7000
6000
5000
4000
3000
2000
1000
0

Difference due tooverflow during
Inundation

23Gauge

24Obs

25

26
27Nov, 2004

TRMM

28

HRM-6hr

2930
GPV

Inundation areas in Hue City..awci_basins_map

Application on Meghna River

SylhetDhaka

GLOBAL DATA
SET

SRTM
90 m DEM

FAO
Soil Map

DERIVEDINPUT DATA & PARAMETERS•Elevation & slope•River channel network•Geomorphology properties•Sub- basins delineation•Soil unit
•Sat. Hydraulic Conductivity•Sat. Soil Moisture•Alpha & water N

LOCALDATA SET

50 m DEM

DigitalSoil MapSatellite
Land use•Max. canopy storage•Max. surface storage•Vegetation coverage

DigitalLand use

NDVI orLAI

Vegetation indexes
e.g. Max. Leaf Area Index

Field work
Inspection

Satellite
RainfallSpatial distributed rainfall

Observed rainfall•Rain gauge•RadarSatellite
Imagery•River morphology•River cross- sections

Measured Riverchannels

Globaldischarge

River stream flow

Observed
River discharge

Satellite
Imagery

Estimation of overflow

Flood Inundationmaps

GIS thematic Maps

Thousands

Hourly Simulation Comparison
100

gpv

80

trmm

60
40
20
0
1/1z

3/13z

6/1z

8/13z11/1z 13/13z 16/1z

DAERAHDAERAH GENANGANGENANGAN DI
DI DASDAS CILIWUNGCILIWUNG

Daerah Genangan..awci_basins_map

12 1113 1014 98PA Karet

63247K.Melayu

2Kali Bata

8

Bukit Duri

3Pengadegan

9Jati Pinggir

4Gang Arus

10Grogol5

Kebon Baru

11Teluk Gong

6

Bidara Cina

13Poglar

12

Kapuk Muara

14

Taman Ratu/ Greenville

1PA Manggarai

Pos Bendung Katulampa

Rawajati

754

Pos Depok

1Prasarana yang perlu
pengendalian
5Pintu Air Manggarai
Pintu Air Kapitol
Pintu Air Jemb. Merah
Pintu Air Tangki
Pintu Air Pasar Ikan
Pompa PluitOptimal gates operation
Error Forecast

Short-term
Weather forecast Output

Spatially distributed data

Observed rainfall

Distributed Hydrological Model
Adjusted
Gates spill

Efficientflood reduction

Optimization
schemeOptimal gates operation..awci_basins_map

Thanks for your kind hosting!

Error forecast Integration to System[mm/h][h]x (t) = Forecast_rain(t) - Obs_rain(t)weightlead _ time =1

σ lead _ time

σ lead _ time

∑ (x − x)=2n∑ (Qsim _ gauge − Qopt _ gauge) }MIN ......{Error = weight *2

N..back to strategy

OBJECTIVE FUNCTION: Flood reduction & future water use
QQ optRun DHM without dam effect

QsimTimeSystem is not
activated

noQsim > Qopt
yesQ

Run DHM with dam effectQsim without dams

Q opt

flood_vol

Qsim without dams
QQ optTimeQuantify flood_volat protecting point
Time

Qsim with dams

Identify critical hour

Qsim with dams critical hour
Release water up to critical hour, but
without exceeding the total inflow

The total amount of water
released, released_vol, is expected
to be similar to flood_vol

adjusted release
Close dams’ gates after critical hour

MIN…..ErrorÆ abs[ flood_vol – released_vol ]..application