![Page 1: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/1.jpg)
Numerical Simulations of Variably Saturated Flow with Energy and Water Phase Change in
Northern Latitude PeatlandCollin A. Macheel
Dept of Mining and Geological EngineeringUniversity of Alaska, Fairbanks
![Page 2: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/2.jpg)
Additional support and contributors Daanen, P. Ronald, Geophysical Institute, University
of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and
Geological Engineering, University of Alaska Fairbanks, P.O. Box 755800, Fairbanks, AK 99775
McGuire, A. David, Institute of Arctic Biology, University of Alaska Fairbanks, Fairbanks, AK 99775
Turetsky, Merritt, Department of Integrative Biology, University of Guelph, Guelph, Ontario N1G2W1
Waddington, Mike, School of Geography and Earth Sciences, McMaster University, Hamilton,Ontario L8S4L8
Kane, Evan, Department of Plant Biology, Michigan State University, East Lansing, MI 48824
![Page 3: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/3.jpg)
•25-30% of the world’s soil C
•Predominant in N Latitudes
•Up to 21% global CH4emissions
•Net sink of CO2
![Page 4: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/4.jpg)
Ecological Research Site Subsurface temperature Water fluxes and climatic data Water level manipulations Microbial populations Vegetative response to WT and
climate manipulations
![Page 5: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/5.jpg)
![Page 6: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/6.jpg)
![Page 7: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/7.jpg)
Hydraulic ConductivityDepth 10-30cm 30-50cm 70-90cm 200-250cm
cm/s cm/s cm/s cm/s
0.03 0.02 2.2x10-3 6x10-7
Abvove MineralBulk
Density Porosity
cm g/ccm %
74.0 0.07 95
69.0 0.09 94
64.0 0.10 94
59.0 0.11 93
54.0 0.14 91
49.0 0.14 91
44.0 0.15 91
39.0 0.14 91
33.5 0.16 90
29.5 0.17 89
24.5 0.18 88
19.0 0.17 89
16.5 0.23 85
14.5 0.33 79
10.0 0.12 92
7.0 0.14 91
4.0 0.25 84
0.0 0.23 85
-4.0 0.44 73
![Page 8: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/8.jpg)
00.20.40.60.8
11.21.41.61.8
0 5 10 15 20
Time (Hr)
WT
(cm
)
0
0.01
0.02
0.03
0.04
0.05
0.06
Prec
ip (c
m)
-0.2
-0.1
0
0.1
0.2
0.3
0.4
0.5
0.6
0 5 10 15 20 25 30
Time (Hr)
WT
(cm
)
00.0050.010.0150.020.0250.030.0350.040.0450.05
Prec
ip (c
m)
Diurnal Temp Fluctuations
-5
0
5
10
15
20
25
30
35
0 100 200 300 400 500 600 700 800
Hours
Deg
ree
C
Air-2cm-10cm-25 cm-50 cm
![Page 9: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/9.jpg)
ρB increases w
/ depth
λ,С
Porosity decreases w/ depth
Thermal Conductivity Models
0.0000
0.1000
0.2000
0.3000
0.4000
0.5000
0.6000
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
VWC (%)
Ther
mal
Con
duct
ivit
y (W
/mK
Kettridge,2008
Kellner, 2009
Letts, 1999
Hayashi, 2008
Weiss, 2006
McKenzie, 2007a
McKenzie, 2007b
![Page 10: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/10.jpg)
Thermal conductivity functions
Weighted average of soil constituents(De Vries, 1963)
Similar to an Arithmetic mean (Forouki, 1986)
Summation of Heat Capacities
)()()()( θθ
θθθ ++++−= oo
ofwf
f
ooa kkfknk
aaoow
aaaoooww
fyfyfkfykfykfk
++++
=
i
n
iiCfC ∑
=
=
=3
1
![Page 11: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/11.jpg)
Soil Moisture Retention
0 500 1000 1500 2000Pressure (cm)
Vol.
Moi
stur
e C
onte
nt
SFCC
-13-11-9-7-5-3-1
Temp (C)
Vo
l. M
ois
ture
Co
nte
nt
Fre
eze
Fu
nct
ion
(1/
C)
SFCFrz Fuct
![Page 12: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/12.jpg)
Applied in HYDRUS (2D/3D)
Kettridge vs Observed
1011
1213
141516
1718
1920
0 50 100 150 200 250 300
t (hours)
T (d
egre
e C
)
ObservedSimulation
R2=0.54, RMSE=7.5C
![Page 13: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/13.jpg)
Applied in HYDRUS (2D/3D)
Weiss vs Observed
10
11
12
13
14
15
16
17
18
19
20
0 50 100 150 200 250 300
t (hours)
T (d
egre
e C
)
ObservedSimulated
R2=0.75 RMSE=1.7C
![Page 14: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/14.jpg)
Inverse Solution Tools
Observed vs Simulated
10
11
12
13
14
15
16
17
18
0 50 100 150 200 250 300
time (hours)
T (d
egre
e C
)
ObservedSimulated
R2=0.90, RMSE=1.2C
![Page 15: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/15.jpg)
Inverse Solution Tools
Observed vs Simulated
10
11
12
13
14
15
16
17
18
0.000 50.000 100.000 150.000 200.000 250.000 300.000
t (hours)
T (d
egre
e C
)
Observed
Simulated
R2=0.989, RMSE=0.37C
![Page 16: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/16.jpg)
R=0.83, R2=0.74, SE=1.78 n=113
•Assuming Soil Surface Temp dependent on three variables
•Air T, WT and time/season
•Regression characterized by multi-variable paraboloid
![Page 17: Numerical Simulations of Variably Saturated Flow with ... · of Alaska Fairbanks, Fairbanks AK 99775-7320 Misra, Debasmita, Department of Mining and Geological Engineering, University](https://reader034.vdocument.in/reader034/viewer/2022042708/5faf5f70177fc808d9730f2f/html5/thumbnails/17.jpg)
Acknowledgements Dr. David McGuire for making travel to the
annual AEG Conference at Lake Tahoe possible. Support from the National Science Foundation grant DEB-0425328, the Bonanza Creek Long-Term Ecological Research program (funded jointly by NSF grant DEB-0423442 and USDA Forest Service, Pacific Northwest Research grant PNW01-JV11261952-231. Also legal advice from the attorneys at Dewey, Cheatum and Howe.