Frankenberger - Tile Drainage in SWAT
Frankenberger - Tile Drainage in SWAT
Frankenberger - Tile Drainage in SWAT
Jane Frankenberger
Dept. of Agricultural & Biological Engineering, Purdue University
With contributions from Chelsie Boles and Daniel Moriasi
Outline
Hydrologic effects of tile drainage in
Midwest landscapes
Factors controlling tile flow
SWAT drainage routines
Remaining needs
X
Surface
Runoff?
Surface runoff lasts a few hours; Tile flow lasts
a few days. Both are hydrograph peaks
2.0
1.8
Surface runoff Tile Flow
1.6
1.4
Flow, cms
1.2
1.0
0.8
0.6
0.4
0.2
0.0
Apr-06 May-06 Jun-06 Jul-06 Aug-06 Sep-06 Oct-06 Nov-06
Tile flow is water that has infiltrated;
Curve number needs to reflect that.
SCS Curve number was developed to
reproduce the hydrograph
It is empirical with no physical basis
Logically, if it is used to separate surface
runoff and infiltration, it needs to be greatly
reduced in tile-drained landscapes, often
by 30% or more
4. 1.
FactorsSaturated
Impermeable
2. Surface
hydraulic
controlling
layer
roughnesstile
– drain
depth
3. Drainage coefficient
conductivity (size
of the soilof tiles)
but flow
especially permeability
SWAT tile drainage routines
Older version (since at least SWAT2005)
based on drawdown time
– Defined TDRAIN as the time to drain soils to
field capacity, set by the user as a static
parameter.
– GDRAIN, a lag coefficient
– Large storm or small, the time of drainage is
the same.
SWAT tile drainage routines
Moriasi et al. (2007a) developed a new
drainage simulation method using the
Houghoudt and Kirkham drainage
equations and a drainage coefficient
– Some call this the “DRAINMOD routines”
– I think H-K-DC would be a good term
Hooghoudt - when water table
below the surface (no surface flow)
Kirkham – when water table
above the surface
1.2
1.0
0.8
0.6
0.4
0.2
0.0
Apr-06 May-06 Jun-06 Jul-06 Aug-06 Sep-06 Oct-06 Nov-06
Remaining Needs
The impermeable layer – Depth
and permeability need to be separated
“DEP_IMP” described as a depth, but actually
controls permeability (seepage through the layer)
Seepage factor
should generally
be < 0.01.
Drainage is highly
sensitive to this
factor
Simulating phosphorus in tile
drains
Need to simulate
through macropores
50
Discharge (MGD)
40
30
20
10
0
1 2 3 4 5 6 7 8 9 10 11 12
Month
EPA-ECHO point source –
Phosphorus by month
0.9
2008 2009 2010 2011 ave
0.8
0.7
0.6
TP (mg/L)
0.5
0.4
0.3
0.2
0.1
0.0
1 2 3 4 5 6 7 8 9 10 11 12
Month
Effect of conservation practices if
implemented on all appropriate land