CC BY
49
УДК 631.43
Ahmed Y.Mady13 and Evgeny V. Shein2
1Lomonosov Moscow State University, Soil Science Faculty, Soil Physics and Reclamation Department. leninskie Gory. 119991, Moscow, russia 2Dokuchaev Soil Science institute, Pyzhevskiy per. 119017, Moscow, russia;
3Ain Shams University, faculty of Agriculture, Soil Science Department, Hadayek Shoubra 11241, cairo, Egypt E-mail: evgeny.shein@gmail.com
modeling soil thermal diffusivity as a function
of soil moisture
Soil thermal parameters are playing an important role in soil heat transfer. Mathematical models and Pedotransfer functions (PTFs) are considered the modern soil physics. Thermal diffusivity is an important parameter of the soil thermal properties, used to study the soil heat flow. Recently, PTFs are widely used to forecast thermal properties. The aim of the work is to ascertain the relation between thermal diffusivity and soil moisture by a model. The relationship between thermal diffusivity and soil moisture by is described by a quadratic equation, and then determined the parameters of this equation using PTFs. The determination coefficient (R2), Root Mean Square Error (RMSE) and Willams-Kloot test exhibited that the proposed PTF of the quadratic equation is better than Arkhangelskaya (2004) model for predicting soil thermal diffusivity as a function of soil physical properties.
Key words: soil physical properties, thermal parameters, Soil moisture, pedotransfer functions (PTFs), mathematical models.
Introduction
Soil thermal parameters are mainly input for modeling of soil heat flux. Mathematical models and PTFs are important tools for predicting soil heat transfer. PTFs are depended on soil physical characteristics (particle size distribution, bulk density and organic matter). Thermal diffusivity is one of the important soil thermal parameters. Thermal diffusivity
K =
x
c • v
is defined as the ratio of the thermal conductivity (k) to the volumetric heat capacity (Cv) [5]. It is a parameter that quantifies the ability of soil to store thermal energy during soil heat transfer processes. Thermal diffusivity is the controlling in soil heat flux during transient conductive heating processes. The literature contains many studies the influence of soil moisture content, bulk density, tillage system, salt concentration, and organic matter on soil thermal parameters [1]-[3]. There are many methods used to measure thermal dif-fusivity, which may be divided into two types: (1) direct methods have been measured in the laboratory according to Kondratieff [5] are more accrue however, costly and time-consuming and (2) indirect methods based on mathematical models and Pedotransfer functions (PTFs), it is depending
on describe the relationship between soil thermal diffusivity and soil moisture by a model, then estimate the parameters of this model by PTFs. PTFs estimate soil thermal parameters based on various physical and chemical soil characteristics. Particle size distribution, soil bulk density and soil organic matter are commonly used predictor variables in PTFs.
There are many mathematical models study soil thermal parameters by defining equation, such as a model of Chung and Horton [3] describe the relation between soil thermal conductivity and soil moisture
x -(e )= b + b2e + b3e05.
The empirical parameters bx, b2 and b3are depended on soil physical properties. HYDRUS1-D program can be estimated those parameters depend on soil texture.
On the other hand, Tikhonravova and Khitrov [6] estimated soil thermal diffusivity as a function of soil moisture by polynomial equation,
K = K0 + aft + a2e2 + a3e5.
Where, K0, ap a2 and a3 are the parameters of the equation. Moreover, Arhangelskaya [1] suggested another kind of model: a lognormal equation dependence on thermal diffusivity from soil moisture. She described the relationship
56 ВЕСТНИК Оренбургского государственного университета 2016 № 12 (200)
between thermal diffusivity and water content by S shape curve.
K = K0 + a exp
f ( W
w
- 0.5
V wo J
Where, w is the water content, K is the corresponding thermal diffusivity; K0, a, w0, and b are parameters ofthe curve: K0 is thermal diffusivity of dry soil, w0 is the water content corresponding to the maximum thermal diffusivity, and (K0+a) is the maximum thermal diffusivity at w=w0.The parameter b characterizes the peak width of the curve. However, this model includes many (four) parameters. So, it is important to calculate soil heat transfer, and thermal diffusivity without data, is to estimate a model describes the relation between thermal diffusivity and soil physical parameters using PTFs. The objective of this work is to propose a mathematical model that will allow to estimate the thermal diffusivity from soil water content based on soil physical properties, and the tasks were several: (1) to measure soil thermal diffusivity under different values of soil moisture by the direct method; (2) to describe a relation between thermal diffusivity and soil moisture by a model; (3) to estimate parameters of this model based on soil physical characteristics using PTFs and (4) compare computed K that obtained by the proposed PTF with the Arkhangel'skaya PTF [1] to estimate a validation of this model.
Material and methods
Soil samples
The silt loam and silt clay loam soil representing major Agriculture in sod-podzolic soils of the Moscow region, Russia, Zelenograd field laboratory of Soil Science Institute. Four soil profiles were dug and, undisturbed soil samples were collected according to the difference in the depth. Soil samples (13 samples) were collected and soil thermal diffusivity was measured under different values of soil moisture (7 levels for each sample), N=91 repeats. PTF was estimated based on soil physical properties. Three soil profiles were dug down, to determine and propose PTFs. Moreover, one soil profile was dug down to test and evaluates the efficiency of this PTF.
Soil Physical Properties
Particle Size Distribution was measured using the pipette Method according to [5]. Soil bulk density (B.D.) was determined by core method, organic matter (OM), which was measured using high temperature using (Express analyzer AN-7529), soil particle density was determined by pycnometer according to [5].
Experimental thermal diffusivity
One method for measuring thermal diffusivity K directly in the laboratory is depended on placing a heat source, having a constant temperature in contact with the surface of a soil column having constant cross-sectional area and insulated sides. Soil thermal diffusivity (Kexp) was determined in the laboratory by the direct method using Kondratieff method according to [5]. Soil cores were initially saturated with water and then subjected to drying, thermal diffusivity was measured under different values of soil moisture. The levels of soil water content were 0.05, 0.15, 0.2, 0.25, 0.28, 0.35, and 0.4 g/g.
Mathematical models
The model of Arkhangel'skaya [1] (KPTFs-2), and the proposed quadratic equation (KPTFs-1) were used to estimate thermal diffusivity as shown in the table (1).
1-The parameters ofa model ofArkhangel'skaya, K0, a, w0, and b were determined from PTFs.
( w ^
K = K0 + a exp
- 0.5
V wo J
K0 =-1.06 + 1.98.B.D + 0.20.C;
a = -0.58 + 2.43.B.D - 0.08.C;
b = 0.12 + 0.12.B.D + 0.12.C .
Where, B.D is the bulk density (Mg/m3) and C (%) is the organic Carbone concentration.
2-The quadratic interpolation equation was to estimate thermal diffusivity K (9 ) as a function of soil moisture.
K(9)= b + b29 - b392.
b
b
BECTHUK Opeheypeckoeo eocyoapctbehhoeo yhubepcuteta 2016 № 12 (200) 57
The parameters of equation bp b2 and b3 were estimated by PTFs based on soil physical properties
[4].
Statistical analysis
The efficiency of the Pedotransfer functions (PTFs) was determined using the determination coefficient (R2), Root Mean Square Error (RMSE) and Williams-Kloot test.
RMSE =
ZYm - Yc )
N
1
ZYm - Yc 1
N
N
Where, Ym is the measured value of K, Yc is the corresponding calculated.
Williams-Kloot test:
(Dl + D2)
2
= a • (D2 - Dl),
where, D1 and D2 denote the errors of the model calculated and experimental values, a - is a regression coefficient, if t value for ta >0, the first model is better, if t <0, the second model is better.
Software
Software tools were Microsoft excel, and SPSS program for the Statistical Analysis (Table.1).
Results and Discussion
1-Modeling thermal diffusivity as function of soil moisture.
Variation of thermal diffusivity of silt loam and silt clay loam of soil samples as a function of soil moisture is shown in Fig. 1. The thermal diffusivity values of soil samples were varied in the range
from 1.21x10-7 to 2.417x10-7 m2/s, the range of soil moisture was 0.05 to 40 g/g. Thermal diffusivity at first increased rapidly with increasing water content to reach the maximum, then decreased at a slower rate. The reason of that water content increase thermal contact between soil particles and replaces the air, which has lower thermal conductivity than water and increases volumetric heat capacity between soil partials sequence increasing soil thermal diffusivity. But the thermal diffusivity is increasing more rapidly than the volumetric heat capacity, as a result of that decrease thermal diffusivity. While soil thermal diffusivity increased liner with increasing soil bulk density. This result was agreed with [2]-[3].
We could describe the represented data in Fig.1: by H shaped curve using a quadratic interpolation equation.
k(e)= b1 + b2e - b3e2.
Where, bp b2 and b3 are the experimental parameters depended on soil physical properties and e is a fraction of mass soil moisture [4].
2-Parameterization and Pedotransfer functions (PTFs)
The parameters of a quadratic equation were estimated by the fitting curve using SPSS program. Take into account the relations between parameters of the quadratic equation and soil physical properties which more effect on soil thermal diffusivity. Table.2: shows that Correlations between soil physical properties (bulk density, organic matter and clay content) and parameters of the quadratic equation.
We could determine PTFs to calculate parameters of the quadratic equation b1, b2 and b3 as a function of soil physical properties using statistical regression as flowing equations:
Table.1: Different methods used to estimate thermal diffusivity
N
Thermal diffusivity Model Method Parameter
(K ) v exp7 In Lab (Kondratieff Methods) Kondratieff Methods
(KPTFs-1) K(e)=b+b2e-b3e2 Quadratic equation b1, b2 and b3
(KPTFs-2) i- i w T 2' Arkhangel'skaya (2009) Lognormal equation K0, a, w, w0 and b
K = K0 + a exp - 0.5 I w0)
b \
58 BECTHUK OPEHEyPECKOEO EOCy/JAPCTBEHHOEO yHUBEPCUTETA 2016 № 12 (200)
к (e )=ь1 + ь2е - b3e2;
b1 = 2.00 • 1 0-7 exp • (- 0.09.B.D - 0.035.0.М - 0.067.CL); b2 = 2.46 • 1 0-7 exp • (0.57.B.D - 0.079.0.М + 0.011.CL); b3 = 2.53 •l 0-7 exp • (0.79.B.D - 0.107.0.M + 0.015.CL).
Where, (B.D) is bulk density (Mg/m3), (O.M) is organic matter content, and (CL%) is clay content (Table.2).
Table.2: Correlations (R) between Bulk Density (B.D), Organic matter (O.M) and Clay content (CL) with parameters of a quadratic equation b1, b2, and b3
Parameters R
B.D O.M CL
b -0.509 0.684 -0.630
Ь2 0.879 -0.850 0.883
b3 0.799 -0.743 0.785
The efficiency of PTFs for predicating thermal diffusivity.
Table.4: shows that the calculated values of K obtained by two PTFs compared to the corresponding measured value of (Kexp) by the direct method. The results observed that the best PTF could be used to calculate soil thermal diffusivity as a function of soil moisture was the quadratic equation (KPTFs-1), for silt loam and silt clay loam sod-podzolic soil. R2 for (KPTFs-1) was higher than (KPTFs-2) were 0.873 and 0.730. Moreover RMSE
2.8
10
| 2.6
&
> 2.4
(Л ¡fc
T3 2.2
ro
^ 2 <D 2 .C
1.8
0.1 0.2 0.3 soil moisture g/g
0.4
0.1 0.2 0.3 soil moisture g/g
0.4
x 10
-7
30-40 cm
со
I 2.8 I 2.6
(Л ¡fc
t5 2.4 ro E
<D .C
2.2 2
0 0.1 0.2 0.3 soil moisture g/g
0.4
x 10
-7
40-50 cm
со
™ 2.8
1 2.6
w n
тз 2.4
га
^22
Ф 22
*
* *
* * ■
*
*
0.1 0.2 0.3 soil moisture g/g
0.4
Fig.1: Thermal diffusivity (m2/s)at different values of soil moisture (g/g) for different soil depths.
ВЕСТНИК Оренбургского государственного университета 2016 № 12 (200)
59
0
0
3
3
2
0
for (KPTFs-1) was lower than (KPTFs-2) were 8.8 *10-8 and 1.4*10-7. The reason of that the model of Arkhangel'skaya [1] described the relationship between thermal diffusivity and water content by S -shaped curve. However, this description differs from ours as we described thermal diffusivity by H shaped curve that described represented data from experimental thermal diffusivity. Williams-Kloot criterion, demonstrated that Student test was 13.33 and was significant <0.005. It means and ensures that the first PTFs of the quadratic equation is better than PTFs of Arkhangel'skaya [1] for predicting thermal diffusivity under silt loam and silt clay loam soil.
Conclusion
Soil moisture and soil bulk density are soil physical properties have great effect on thermal diffusivity. The effect of soil bulk density and soil moisture are more than clay content on thermal diffusivity. Thermal diffusivity at first increased rapidly with increasing water content then decreased at a slower rate. We described the relationship be-
tween thermal diffusivity and soil moisture by H shaped curve using a quadratic equation (KPTFs-1). Moreover, estimate the parameters of this equation by PTFs. PTFs based on soil physical properties (particle size distribution, bulk density and organic matter). There are strong correlations between parameters of the quadratic equation (b1, b2 and b3) and soil physical properties (bulk density, organic matter content and clay content). As well as there are correlations between parameters of the quadratic equation that measured by fitting curve and by PTF: KPTFs-1 by quadratic equation and KPTFs-2 by the lognormal model of Arkhangel'skaya. R2, RMSE, and Williams-Kloot criterion exhibited that (KPTFs-1) is better than (KPTFs-2). The PTF of the quadratic equation is the accurate equation for forecasting soil thermal diffusivity, for silt loam and silt clay loam sod-podzolic type of soils genesis.
Acknowledgment
This work was supported by the Russian Scientific Fund (the project № 14-16-00065).
10.11.2017
Reference:
1. Arkhangel'skaya T.A. Parametrization and mathematical modeling dependence of termodiffusivity from soil moisture. Pochvove-denie. 2009. № 2. C. 178-188. [in Russian].
2. Bouma, J. Using soil survey data for quantitative land evaluation. Advances in Soil Science.1989. Vol 9. Pp. 177-213.
3. Chung J., R.Horton . Soil heat and water flow with a partial surface mulch. Water Resour. Res. 1987. Vol.23. Pp. 2175-2186.
4. Shein, E.V., A.Y. Mady. Soil Thermal Parameters Assessment by Direct Method and Mathematical Models. Journal of Soil Science and Environmental Management. 2016.7 (10):166-172.
5. Theories and Methods of Soil Physics. 2007. Shein, E. V. and L.O. Karpachevskyi (Eds.) Izd. Grif i K, Moscow [in Russian].
6. Tikhonravova, P.I. and N. B. Khitrov (2003). Estimation of Thermal Conductivity in Vertisols of the Central Ciscaucasus Region." Pochvovedenie, 2003. No. 3, 342-351.
"Article is published for the first time"
Evgeny V. Shein, professor at Soil Physics and Reclamation Department, Soil Science Faculty Moscow State University named M.V.Lomonosov Leninskie Gory 119991 Moscow, Russia. Professor at Dokuchaev Soil Science Institute, Pyzhevskiy per.,7. 119017, Moscow, Russia. Tel: +7 903-539-88-19 E-mail: evgeny.shein@gmail.com
Ahmed Y. Mady, PhD. Student at Lomonosov Moscow State University, Soil Science Faculty, Soil Physicis and Reclamation Department,119991, Moscow Russia. Assistant Lecture at Ain shams university, Faculty of Agriculture, Soil science Department,
11241, Cairo, Egypt. Tel: +7916-990-90-50 E-mail: ahmed_mady@agr.asu.edu.eg
60 BECTHUKopehsyprckoro eocymapctbehhoeo yhubepchteta 2016 № 12(200)
