Transcript Slide 1

SABAE-HW3D
A Hydro-meteorological
Model Coupling
the Land Surface and Groundwater Flow
Youssef Loukili and Allan D. Woodbury
[email protected]
[email protected]
Department of Civil Engineering – University of Manitoba
Winnipeg , Manitoba , Canada
SABAE-HW
Soil Atmosphere Boundary,
Accurate Evaluations of Heat and Water
1. Free choice of soil column depth
and layers in SABAE-HW
2. More accuracy in the water mass
balance terms
3. Efficient coupling with groundwater
models through water table b.c.
4. Complete update of CLASS v 2.6
to Fortran 90/95
5. Extensive validation and
inter-comparisons to other models
(HELP, SHAW, HYDRUS)
Meteorological inputs required by CLASS and SABAE-HW :
1.
2.
4.
6.
Incoming long wave and short wave solar radiation
Precipitation
3. Wind speed
Air temperature
5. Atmospheric pressure
Specific humidity , or relative humidity
Soil characteristics ( Clapp and Hornberger model )
Vegetation parameters
Output :
1.
2.
3.
4.
5.
6.
7.
Surface, canopy, snow and soil layer temperatures
Liquid and frozen water contents
Surface heat flux
Net absorbed short and long wave radiations
Sensible and latent heat fluxes
Water evaporation, pond and bottom drainage
Snow accumulation and snow melt
Comparisons
• Loukili, Y. , Snelgrove, K.R., and A.D. Woodbury, SABAE-HW – An
enhancement of the water balance prediction in the Canadian Land
Surface Scheme, Vadose Zone J., (accepted), 2008.
• Loukili, Y. and D. D. Woodbury, ABAE-HW3D: a MeteorHydrological Model Coupling the Land Surface to Groundwater Flow,
Eos Trans. AGU, 88(52), Fall Meet. Suppl., 2007
• Loukili, Y., Woodbury, A.D. and K.R. Snelgrove, AccuCLASS - an
Enhancement of the Canadian Land Surface Scheme for Climate
Assessment Over the Prairies, Eos Trans. AGU, 87(52), Fall Meet.
Suppl., 2006
• Field Testing Assiniboine Delta Aquifer (ADA)
Qu’A
ppel
le
Rive
r
Lake
Winnipeg
Assiniboine
Lake
Delta Aquifer Manitoba
Assiniboine
River
.
Winnipeg
LL-98
Sour
is
Rive
r
Canada
U.S.A.
Red River
Pine Creek North of Assiniboine Delta Aquifer of Manitoba
n
sa
d 0.50 m
LL-98
y
am 3.25 m
o
L nd
sa
sa
nd
o
Winnipeg
0.35 m
NARR meteorological data (1999-2003) – Resolution 32 x 32 km2
Depth = 20 cm
0.4
Total moisture (cm/cm)
0.35
SABAE
CLASS
0.3
0.25
0.2
0.15
0.1
0.05
1-Jan-99
2-Jan-00
2-Jan-01
3-Jan-02
4-Jan-03
5-Jan-04
4-Jan-03
5-Jan-04
Depth = 20 cm
0.4
Total moisture (cm/cm)
0.35
SABAE
SHAW
0.3
0.25
0.2
0.15
0.1
0.05
1-Jan-99
2-Jan-00
2-Jan-01
3-Jan-02
40
SABAE
SHAW
35
25
20
15
10
5
0
1-Jan-99
2-Jan-00
2-Jan-01
3-Jan-02
4-Jan-03
5-Jan-04
70
SABAE
SHAW
60
50
Daily pond (mm)
Snow depth (cm)
30
40
30
20
10
0
1-Jan-99
2-Jan-00
2-Jan-01
3-Jan-02
4-Jan-03
5-Jan-04
31-Jan-03
0.4
SABAE (liquid)
SHAW (liquid)
CLASS (liquid)
0.35
Moisture (cm/cm)
0.3
SABAE (frozen)
SHAW (frozen)
CLASS (frozen)
0.25
0.2
0.15
0.1
0.05
0
0
0.5
1
1.5
2
2.5
Depth (m)
3
3.5
4
4.5
The meteor-hydrological coupled model SABAE-HW3D
• SABAE-HW operating in each soil column
Time step = 30 min
• 2D horizontal saturated flow using Finite Volume Method
Time step = 5 days
( conservative fluxes )
• Soil columns and grids extend to water table
Updated each 5 days
( grid regeneration allowed )
• Bottom drainage summed up and provided as cell (groundwater) recharge
80
40
40
60
40
40
120
80
140
12
0
120
40
1 60
1 40
100
180
100
200
20
40
20
20
ADA Transmissivity x 103 U.S. Gals/Ft/Day
Z
Y
X
base
390
370
350
330
310
290
270
ADA’s Bedrock elevation
y (m)
5.56E+06
March 25, 2000
top-h
13
9
6
2
5.55E+06
5.54E+06
5.53E+06
5.52E+06
5.51E+06
5.5E+06
5.49E+06
460000
480000
500000
520000
x (m)
Depth to the water table in the ADA (m)
y (m)
5.56E+06
March 25, 2000
h-base
28
21
14
7
5.55E+06
5.54E+06
5.53E+06
5.52E+06
5.51E+06
5.5E+06
5.49E+06
460000
480000
500000
520000
x (m)
Saturated thickness in the ADA (m)
y (m)
5.56E+06
March 25, 2000
h
380.423
376.971
373.519
370.067
366.615
363.162
359.71
356.258
352.806
349.354
345.902
342.45
338.998
335.546
332.094
328.641
325.189
321.737
318.285
314.833
311.381
307.929
304.477
301.025
297.573
5.55E+06
5.54E+06
5.53E+06
5.52E+06
5.51E+06
5.5E+06
5.49E+06
460000
480000
500000
520000
x (m)
Piezometric head in the ADA (m)
y (m)
5.56E+06
March 25, 2000
dlydrain
0.611269
0.577801
0.544334
0.510867
0.477399
0.443932
0.410464
0.376997
0.34353
0.310062
0.276595
0.243127
0.20966
0.176193
0.142725
0.109258
0.0757905
0.0423231
0.00885575
-0.0246116
-0.058079
-0.0915464
-0.125014
-0.158481
-0.191949
5.55E+06
5.54E+06
5.53E+06
5.52E+06
5.51E+06
5.5E+06
5.49E+06
460000
480000
500000
520000
x (m)
Recharge to Groundwater in the ADA (mm)
y (m)
5.56E+06
February 4, 2000
Mice1
0.358116
0.353091
0.348066
0.343041
0.338016
0.332991
0.327966
0.322941
0.317916
0.312891
0.307866
0.302841
0.297816
0.292791
0.287766
0.282741
0.277716
0.272691
0.267666
0.262641
0.257616
0.252591
0.247566
0.242541
0.237516
5.55E+06
5.54E+06
5.53E+06
5.52E+06
5.51E+06
5.5E+06
5.49E+06
460000
480000
500000
520000
x (m)
Distributed frozen moisture in the top 10 cm of soil in the ADA
y (m)
5.56E+06
March 25, 2000
5.55E+06
5.54E+06
5.53E+06
5.52E+06
5.51E+06
5.5E+06
5.49E+06
460000
480000
500000
520000
x (m)
Flow directions and the water divides in the ADA
Frozen moisture – Cross section West-East ADA
Soil Temperatures – Cross section West-East ADA
Cumulative Recharge to WT – Entire ADA
Summary and Conclusions
1. SABAE-HW is inter-compared with SHAW and CLASS models
and proved to be effective for a cold region hydrology.
2. SABAE-HW3D efficiently couples the unsaturated and saturated
flows using weather forcing data.
3. SABAE-HW3D is benchmarked against Seep/w (Geo-Slope
International, 2002).
4. The coupled model is also successfully applied
to the Assiniboine Delta Aquifer of Manitoba.
5. Further simulations are planed using downscaled weather data
and actual vegetation cover.
Liquid moisture – Cross section West-East ADA
Team and Milestones
• Dr. Youssef Loukili, Research Associate
• Alireza Hejazi, Ph.D. student
Team and Milestones
• Field experiments & hypothesis testing:
September 2008 (Hayashi, Loukili, Hejazi)
• Further code development, Fall 2008
• GMS-MODFLOW of ADA (3-D aquifer): May
2008 (Hejazi)
• Refinement of stand-alone SABAE runs
(Hanesiak, NARCCAP of ADA), Fall 2008
(Loukili)
• Full linkage SABAE-HW3D and GCM (single
column), Fall 2008
• Detailed runs over ADA, Fall 2008 (Hejazi)