WISTOO en
WISTOO en
WISTOO en
RAINFALL-RUNOFF
TRANSFORMATION - WISTOO
Ozga-Zielińska Maria
Institute of Environmental Engineering Systems
Warsaw University of Technology, Warsaw, Poland
WISTOO model is adopted for simulation with one hour step for
mountainous and hilly watersheds. Long-term simulations can be
conducted – up to 185 days (summer half-year). Model is
implemented for personal computers with Windows operating
system, and uses digital thematic layers (Fig. 1):
• digital elevation model,
• stream network,
• soil structure,
• land-use structure.
In the algorithm, mainly raster square model of spatial phenomena
description is used. Raster size depends on scale, but computer
program can handle maximum 4 000 000 raster elements (cells).
Computer implementation of algorithm enables tracing of:
• hydrographs in ten arbitrary river cross-sections,
• spatial distribution of: net precipitation, velocity and
depth of surface and subsurface flow, soil moisture,
• influence of water reservoir.
1. INTRODUCTION
EVAPOTRANSPIRATION
PRECIPITATION
INTERCEPTION
NET
PRECIPITATION
EXCESS SURFACE
INFILTRATION
RAINFALL RUNOFF
TOTAL RUNOFF
j
1
å(d n * Pk )
k =1
Pi = j
k
(1)
1
å n
k =1 d k
P2
P1
P4
x1 x2
x4
P5 P3
x 5 Px x3
R i - R i+1
A= (3)
(H i - H i+1 ) R i
where: Ri, Ri+1 - recorded precipitation depth for the lowest and the
highest station in the watershed, [mm],
Hi, Hi+1 - elevation above sea level for the highest and the
lowest station in the watershed, [m].
2.3. Interception
dc c
= (1 - p ) R - E p dla c < s
dt s
(4)
dc
= (1 - p ) R - E p - D o e b ( c - s) dla c ≥ s
dt
Rnett = p R dla c ≤ s
(5)
Rnett = R dla c > s
co
c t = co + [ (1 - p) R - E p ] ∆t (6)
s
where: ct - current water depth on plants at time t, [mm],
co - water depth on plants at time t-1, [mm].
2.4. Evapotranspiration
soil
moisture
change
subsoil
∆ (Rn - G ) + ρ cp (ε - e e ) / ra
Ea = (7)
λ [∆ + γ (1 + rc /ra )]
where: Ea - current evapotranspiration, [mm],
∆ - inclination of water vapor pressure curve vs.
temperature, [hPa/K],
Rn - net radiation, [W/m2],
G - soil heat flux, [W/m2],
ρ cp - volumetric dry air heat capacity, [J/m3K],
ε, ee - maximal and current water vapor pressure, [hPa],
ra - aerodynamic resistance function, [s/m],
rc - superficial resistance function, [s/m],
λ - latent heat of evaporation, [J/kg],
γ - psychrometric constant, [hPa/K].
ba
rc = a + (8)
R n + ca
2
é æ z v - d h öù
êlnçç ÷÷ú
è z øû
ra = 8.06 ë
o
(9)
1 + 0.864 v w
2.5. Infiltration
subsoil supply
soil supply
SUBSOIL
subsurface
INFILTRATION supply
gruntwater
supply
æ h + hk ö
v = K Sh = K ç 1 + w ÷ (11)
è z ø
Hydraulic conductivity.
Hydraulic conductivity is estimated from formula:
α
æ Θ - Θs ö
K(Θ) = K o çç ÷÷ (12)
è Θo - Θs ø
æ ψ g -ψ d h ö
v = K çç + 1 + w ÷÷ (14)
è z z ø
3
2
1
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7
humidity
V = R ∆t (15)
3 3
V = å vi ∆t ∆Θi = å z i ∆Θi (16)
i =1 i =1
where: ∆Θi – width of humidity front, [-].
K (Θ i ) - K (Θ i-1 )
vi = (17)
Θ i - Θ i-1
Outputs from the model are: exceed rainfall, input values for
subsurface runoff and groundwater runoff.
Initially conditions.
Key element is determination of initially conditions of the process,
i.e. initial humidity. This problem is solved with assumption that
volume of water supplying groundwater reservoir is equal
watershed runoff:
nz nz
Q k = å Fi Rg i = å Fi Rp i (18)
i =1 i =1
parcolation
supply from
subsurface runoff
surface runoff
supply
infiltration into
subsoil
∂ qx ∂ q y ∂ h
+ + = RE
∂x ∂ y ∂t (19)
q = αk h N
In used model parameters are set like for turbulent flow, according
Manning theory:
Sz
αk =
nh
(20)
5
N=
3
∂ qx q i+1 - q i2+1 qi - qi
≈β 4 + (1 - β ) 4 2 (21)
∂ x ∆x ∆x
∂ qy q i+1 - q i3+1 qi - qi
≈β 1 + (1 - β ) 1 3 (22)
∂ y ∆y ∆y
∂ h h iśr+1 - h iśr
≈ (23)
∂t ∆t
RE = β R iE+1 + (1 - β ) R iE (24)
q3
y k+1
q4 q2
hśr
X
yk
y
q1
xi x i+1
x
αβ
where: a = -
∆x
1
b =
∆t
β q id+1 (1 - β ) (q id - q io ) h iśr
c= + - - β R iE+1 - (1 - β ) R iE
∆x ∆x ∆t
qd - inflow sum for cell,
qo - outflow sum for cell.
h iśr+1 = h io + ∆h (26)
f(h śr )
where: ∆h = - ,
f(h o )
ho – water depth prom previous time step.
R gr =
( )
(q grj - q dj ) β + q oj-1 - q dj-1 (1 - β )
+
Hg µ - h śrj-1 (1 - β ) R Ej-1
- (27)
∆x β ∆t β β
where: qgr – maximal outflow from cell calculated from soil layer
depth Hg and filtering coefficient µ.
Rw =
( )
- q dj β + q oj-1 - q dj-1 (1 - β ) h śrj-1 (1 - β ) R Ej-1
- - (28)
∆x β ∆t β β
q g = K f h g Sgr (29)
qg
qi = (30)
Nz
L r K f Sgr ∆t j R g ∆t
h gj+1 = 1 - hg + (31)
µ Fg µ
where: Lr - total streams length, [m],
Fg - groundwater reservoir supply area, [m2],
Rg - average groundwater supply from aeration zone, [m/s].
∂Q ∂F
+ = qb
∂x ∂t (32)
Q = α k R hN-1F
groundwater
runoff
Figure13. Transformation process in streams calculation scheme
DIGITAL ELEVATION
MODEL PROCEDURES
TOOL WATERSHED PARAMETERS
PROCEDURES PROCEDURES
METEOROLOGICAL DATA
PROCEDURES
CALCULATION PROCEDURES
RAINFALL-RUNOFF
TRANSFORMATION WITH
HYDROGRAPH VISUALIZATION
RAINFALL-RUNOFF
TEMPORAL TRANSFORMATION WITH DATA
VISUALIZATION
SPACIAL VISUALIZATION OF STORAGE
OF PROCESSES SELECTED PROCESSES
RUNOFF TRANSFORMATION
THROUGH
WATER RESERVOIR
0 0 0
0 1 1 1 0
2 3 2 0
Figure 17. Slopes in Łososina watershed – lighter areas correspond to flat terrain
(minimal slope)
There are five types of soil in the watershed: clayey, skeletal,
fen, dust, sandy. Skeletal soils are in initial phase of evolution.
Most of the area is covered by clay and silt soils with rock waste,
stones, fine stone bits, sand.
Land-use (Fig.18) was divided into five classes: fields and
croplands (1), forest (2), sparse rural areas (3), dense rural areas
(4), main asphalt roads (width above 15m) (5).