The Scientific World Journal

Volume 2013 (2013), Article ID 421762, 7 pages

http://dx.doi.org/10.1155/2013/421762

Research Article

Performance Evaluation of Volumetric Water Content and Relative Permittivity Models

Muhammad Mukhlisin1,2 and Almushfi Saputra¹

¹Department of Civil and Structural Engineering, Faculty of Engineering and Built Environment, Universiti Kebangsaan Malaysia, 43600 Bangi, Selangor, Malaysia

²Department of Civil Engineering, Polytechnic Negeri Semarang, Jl. Prof. Soedarto, SH, Tembalang, Semarang 50275, Indonesia

Received 25 July 2013; Accepted 12 September 2013

Academic Editors: G. Ganjegunte and J. J. Wang

Copyright © 2013 Muhammad Mukhlisin and Almushfi Saputra. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

In recent years many models have been proposed for measuring soil water content () based on the permittivity () value. Permittivity is one of the properties used to determine in measurements using the electromagnetic method. This method is widely used due to quite substantial differences in values of for air, soil, and water, as it allows the value to be measured accurately. The performance of six proposed models with one parameter (i.e., permittivity) and five proposed models with two or more parameters (i.e., permittivity, porosity, and dry bulk density of soil) is discussed and evaluated. Secondary data obtained from previous studies are used for comparison to calibrate and evaluate the models. The results show that the models with one parameter proposed by Roth et al. (1992) and Topp et al. (1980) have the greatest data errors, while for the model with two parameters, the model proposed by Malicki et al. (1996) agrees very well with the data compared with other models.

1. Introduction

Measurement of water content () in soil has become a major component of the various fields of geotechnical analysis. Measurement of is needed to support many bodies of research related to the soil [1]. For example, in agriculture, is an important factor for irrigation and crop quality maintenance. On the other hand, in hydrology, determining the rate and quantity of water movement in soil requires the measurement. Meanwhile, in forestry, is required for information on the water storage capacity of soil. Besides, also affects the stability of the slope of soil due to its relationship with the soil strength [2–9].

Measurements of can be categorised as direct and indirect measurements. Gravimetric measurement is a direct measurement which is categorised as a conventional method. In this method, the value of is determined by subtracting dry from wet soil sample weights. This method is very accurate but it is not practical due to the long time it takes to get the result. However this method is used as a calibration for other techniques.

For indirect measurement, electrical methods for measuring have primarily been subjected to extensive study due to their ease and practicality of use. These methods have been widely used and discussed in many previous studies (e.g., [10–14]). Moreover, issues of the instruments used for measurement of water content, from small-scale (<1 m²) to large-scale (100 m²), and suitable methods for measurement at those various scales have also been discussed in detail by [15]. In their study, Robinson et al. [15] concluded that the method for measuring water content requires improvement. In some other studies (e.g., [16–18]) electrical properties are measured to get the characteristic of . In their study it can be seen that permittivity () measurement can be used to predict . Permittivity () is the most common of electrical properties that are used to measure . Although there are also some techniques by measuring the capacitance of soil, then it converted into (e.g., [19, 20]).

To represent the relationship, there are several models that have been proposed in the last few decades (e.g., [14, 18–30]). To produce the relationship model, most of them used gravimetric measurement for data calibrations. However, data from previous studies were also used (e.g., [25]).

All models can be categorised as having one parameter or two or more parameters. Models with one parameter only involve the relationship between the permittivity and water content, whereas models with two or more parameters include other parameters such as porosity or dry bulk density. This classification is used to analyse the influence of parameters other than the permittivity that affect the value of water content, and no previous study has tried to analyse this.

In this study, the models proposed by [14, 18–23, 25, 28–30] were reviewed and compared using secondary data from previous studies (e.g., [22, 27, 31, 32]). The data are then used to determine which model has a significant relationship.

2. Theory and Method

Many equations have been proposed for calibration of and . These models can be divided into two categories: models with one parameter and models with two or more parameters.

2.1. Model with One Parameter

There are some proposed models that show the relationship. Topp et al. [22] successfully introduced the relationship that is commonly used in the geotechnical area for the first time. The relationship is

where is the relative permittivity or dielectric constant and is the volumetric water content of the soil. The experiment uses Time Domain Reflectometry (TDR) at a frequency between 1 MHz and 1 GHz to measure of several mineral soils. Then a polynomial fitting is used empirically to obtain the relationship model. The estimated error in this model is 0.013 [22]. In their study, Topp et al. [22] also provide another form of (1a) as follows:

The relationship models for organic soil and 450 m glass beads are also shown as follows:

Roth et al. [21] used miniprobe TDR for their experiment to propose another empirical equation for the relationship, which had been used previously by [33]. The relationship for mineral soil proposed by [21] is

while the relationship for organic soil and material is

The error estimations of these equations for mineral soil and organic soil are 0.015 and 0.035 cm³ cm⁻³, respectively [21].

Ferré et al. [25] proposed a simple equation for the relationship. This equation was generated from the principle of dielectric mixing models and using TDR without coatings:

A simple equation was also introduced by [29]. They used 505 measurements from organic forest floor sample experiments using TDR, where the relationship is

A coaxial transmission system at a frequency of 100 MHz was used by [18] to produce another model of the relationship, which used a wide range of soil textures samples:

Permittivity based on capacitance measurement was investigated by [20]. They proposed an empirical model from experiment using a type of quartz sand with particle sizes in the range 0.15–0.9 mm:

where , , , and .

2.2. Model with Two or More Parameters

Some relationship equations for permittivity and soil water content were also influenced by other parameters such as porosity and bulk density. By using the concept of mixing models and using data from other studies (e.g., [34–36]), Wang and Schmugge [30] proposed the following equations:

Equation (7a) is used for , while for the following equation is used:

where , , , and are the permittivity of ice, water, air, and rock, respectively (i.e., , , and ), while is the transition moisture (0.16–0.33), is the porosity of soil (0.5), and is the fitting parameter (0.3–0.5) [30].

Roth et al. [28] proposed the equation based on the dielectric mixing model which has been described by [24]. The experiments were carried out by measuring a wide range of soil types using TDR with the error value of soil water content, no more than 0.013 cm³ cm⁻³ [28], with forms of the following equation:

where for three phases in series and for three phases in parallel.

Another model was proposed by [23]. They conducted experiments using TDR and 62 kinds of soil samples consisting of mineral soils, organic soil, standard pot soils, artificial peat-loess and peat-sand, sea and river sand, forest litter, and so forth, which differ in terms of texture and bulk density, which gives an uncertainty of soil water content of 0.03 [23]:

Gardner et al. [19] used capacitance measurement methods to obtain soil water content with soil dry bulk density values ranging rom 1.08 to 1.49 and then used multiple linear regression analysis to best fit the measurement data, resulting in the following equation:

where is dry bulk density.

Robinson et al. [14] developed an equation for coarse textured, layered soils by using TDR and coarse-grained, glass bead, and quartz grains:

where and are the permittivity values for dry and saturated soil, respectively.

Table 1 shows the 11 proposed equations of the relationship for one and two or more parameters. It provides a brief explanation, including the experimental method, soil type, properties of the soil, and the sources information for each proposed equation.

Table 1: Summary of all equations of relationship.

Equations
Source
Experimental method
Soil type
Properties of soil

Porosity

(cm³cm⁻³)
Bulk density

(g cm⁻³)
Particle density

(g cm⁻³)

Model with one parameter

(1a)

(i) Mineral soil

(ii) Organic soil

(iii) Vermiculite

(iv) Glass beads
—
(i) 1.04–1.44

(ii) 0.422

(iii) 1.08

(iv) 1.49–1.61
—

(1b)
Topp et al. [22]
ε: using TDR Tektronik Model 7S12 to perform 18 experiments with different treatments

θ: using gravimetric technique

(1c)
Organic soil
—
0.422
—

(1d)

450 μm glass beads
—
1.60–1.61
—

(2a)
Roth et al. [21]
ε: TDR miniprobe 250 ps rise time needle pulse

θ: gravimetric technique
9 Mineral soils
0.418–0.482
1.26–1.55
2.28–2.67

(2b)
7 Organic soils
0.527–0.785
0.2–0.77
0.70–1.63

(3)
Ferré et al. [25]
Using model of inverse averaging for TDR method by analysing the mixing model
—
—
—
—

(4)
Schaap et al. [29]
ε: TDR Tektronix 1502B

θ: gravimetric technique
25 samples of forest floors
—
0.086–0.263
1.3

(5)
Curtis [18]
Coaxial Transmission/reflection apparatus controlled by a Hewlett-Packard 8510C Vector Network Analyzer system 45 MHz to 26.5 GHz
—
—
—
—

(6)
Wu et al. [20]
ε: based on capacitance measurement

θ: gravimetric technique
Quartz sand
—
—
—

Model with two parameters

(7a)
Wang and Schmugge [30]
Modelling using data from other studies [34–36]
22 different samples
0.4–0.6
1.1–1.7
2.6–2.75

(7b)

(8a)
Roth et al. [21]
TDR
From 11 different field sites

(8b)

(9)
Malicki et al. [23]
TDR CAMI
62 kinds of soil samples
0.33–0.95
0.13–1.66
1.06–2.7

(10)
Gardner et al. [19]
Capacitance probe 80–150 MHz
(i) Brown earths

(ii) Silica materials
—
(i) 1.08–1.49

(ii) 1.24–1.63
—

(11)
Robinson et al. [14]
TDR Tektronix 1502B
Coarse grained, quartz grain, sandy soil
—
—
—

Models with one parameter use only a single parameter or variable to calculate soil moisture content. This parameter is permittivity. Most of these models are defined using empirical methods. Models with two or more parameters have other parameters besides permittivity, such as porosity and bulk density.

2.3. Secondary Data

In this study, reference data are needed to test the ability of all these models. For this purpose, secondary data obtained from previous experimental studies that showed the relationship between soil water content and permittivity were used. Overall, there are 44 secondary data, and the data sources and soil porosity used in this study can be seen in Table 2. These data were generated from experiments using different methods such as TDR (e.g., [14, 21, 22, 26, 28, 31, 32]), capacitance probe [20], and frequency domain [27] and also from varying types of soil.

Table 2: Source of secondary data and porosity.

Porosity
Data source

0.30–0.38
Friedman [26]; Sabouroux and Ba [31]

0.40–0.48
Friedman [26]; Hilhorst et al. [27]

0.5
Skierucha et al. [32]

0.62–0.66
Roth et al. [21]; Friedman [26]

—
Topp et al. [22]; Malicki et al. [23]; Robinson et al. [14]; Curtis [18]; Gardner et al. [19]; Dobson et al. [24]; Wang and Schmugge [30]

The value of relative permittivity air, water, dry soil, and saturated soil can be seen in Table 3. Relative permittivity of dry soil and saturated soil is obtained from the study of the previous researches [14]. Relative permittivity of material is affected by the chemical components of its constituent and can be calculated by using the mixture model [37].

Table 3: Relative permittivity of material properties.

Material
Relative permittivity
Chemical elements

Air (ε air)
1
N₂, O₂

Water (ε water)
80
H₂O

Dry soil (ε dry soil)
2–4
N, P, K, Ca, Mg, S, Cu, Zn, Fe, Mn, B, Cl, Na, H

Saturated soil (ε saturated soil)
23–28
H₂O, N, P, K, Ca, Mg, S, Cu, Zn, Fe, Mn, B, Cl, Na, H

3. Results and Discussion

The effects of porosity on the relationship are shown in Figure 1. Four secondary data samples which have porosity ranging from 0.30 to 0.66 are highlighted.

Figure 1: Secondary data of volumetric water content and permittivity as a function of porosity.

Figure 1 shows that the smaller the soil porosity, the greater the value of permittivity for a given value of volumetric water content. In this condition, the pores in the soil will be filled by water and air. Therefore when the porosity is large, then the rest of the pores are filled by air. This corresponds to the concept of dielectric mixing used in models by [24]. When most of the volume fraction of soil pores is filled by air, it donates a small value of the total permittivity of the soil. This figure also shows that the spread of data does not occur significantly for small water content (0–0.1). In this condition, the value of permittivity is in the range of 1–5. Otherwise, when the water content begins to increase, it produces scattered data values.

3.1. Model with One Parameter

Figure 2 shows curves for (1a) to (6), which have one parameter. All equations appeared to cover all of the available data. However, each equation appears in a certain position within the data. Equations (1a), (1b), (1d), (3), and (5) are quite close to each other and tend towards the upper part of the data which have relatively small porosity (<0.5). On the other hand, the larger porosity (>0.6) is occupied by (4).

Figure 2: Comparisons using all data for (1a) to (6).

Almost all of the curves show a similar trend, except for (2a) and (6). The curve for (2a) indicates a lower increase in permittivity when the water content is greater than 0.4, while the curve for (6) gives a constant permittivity value when the water content is greater than 0.4. Wu et al. [20] explain that this reduction in the increase of permittivity is due to the effect of saturation in the soil.

Equations (1c) and (2b) occupy the central part of data distribution. These curves provide a reasonably safe prediction of the relationship for equations with one parameter or without any parameters of porosity. Table 4 shows the -square of each curve to data and also the Root Mean Square Error (RMSE) of each equation to data. Equation (2b) gives a better result for -square and RMSE compared with other equations with one parameter.

Table 4: -square and root mean square error (RMSE) of the equations to data.

Equations
(1a)
(1b)
(1c)
(1d)
(2a)
(2b)
(3)
(4)
(5)
(6)

0.749
0.8189
0.824
0.480
0.807
0.877
0.834
0.612
0.834
0.705

RMSE
6.131
0.078
5.133
8.824
0.071
0.058
0.075
0.114
0.075
6.536

3.2. Model with Two Parameters

Figure 3 shows the effect of porosity ( to ) on the suitability of (7a) to (11) with data and also displays some of the data with a value of porosity (0.33, 0.44, and 0.62) in order to see the fit between data with a curve based on the value of porosity. In Figure 3(a), it can be seen that (7a) only fits in a certain small area of the data, though with different porosity. In this equation, the effect of changing porosity is not significant. Figures 3(b) and 3(c), with all the possible values of porosity, show that neither of these equations is quite enough to follow the pattern of the data. In these equations, the trend is linear for both of the graphs.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 3: Comparison of (7a) to (11) with all data and different porosity.

A better plot is shown in Figure 3(d), where the equation occupies all of the data well. This figure shows that (9), which was proposed by [23], has a significant effect on changing porosity. The figure also shows that curves merge very well in the range of the secondary data.

Overestimated results are produced in Figure 3(e). Curves with small porosity parameters can not even cover the area of data. There are only two curves passing through the area of data. Nevertheless, these curves are inconsistent with the position of the porosity of the data.

Figure 3(f) shows a wide spread of curves for changes of porosity as the water content increases. Almost all areas are covered except data for water content values​ ​ smaller than 0.2. However, when viewed in terms of the porosity data, the curves in this figure do not look quite as good because they spread without following the porosity data.

4. Conclusion

A comparison of some equations for the relationship was performed to provide an overview of the efficacy and ability of each equation in describing the relationship and its correlation to the porosity of soil. In this study, secondary data were used as a reference to compare the equations. Secondary data with porosity values were plotted in one graph to show the effect of soil porosity on the relationship between water content and permittivity. For the same water content, the permittivity of soil decreases with increasing porosity. In this case, porosity should be taken into account when considering the relationship.

There are some models of the relationship that do not fit with data in a certain range. For equations with one parameter, the model of [21] for mineral soil and the model of [20] indicate curves which tend to slow down the increase in permittivity at near saturation conditions of soil. Furthermore, equations with one parameter are also not able to cover all areas of the data.

However, of the models that included a porosity parameter, apparently not all of them could explain the effect of porosity on the relationship very well. From this study, only the equation proposed by [23] gives a fairly good conformity of data for different porosity values of soil.

Acknowledgment

The authors would like to express their appreciation to the financial support provided by the Ministry of Higher Education (MOHE) of Malaysia.

References

N. Osman and S. S. Barakbah, “Parameters to predict slope stability-Soil water and root profiles,” Ecological Engineering, vol. 28, no. 1, pp. 90–95, 2006.

M. Mukhlisin, K. Kosugi, Y. Satofuka, and T. Mizuyama, “Effects of soil porosity on slope stability and debris flow runout at a weathered granitic hillslope,” Vadose Zone Journal, vol. 5, no. 1, pp. 283–295, 2006.

M. Mukhlisin, M. R. Taha, and K. Kosugi, “Numerical analysis of effective soil porosity and soil thickness effects on slope stability at a hillslope of weathered granitic soil formation,” Geosciences Journal, vol. 12, no. 4, pp. 401–410, 2008.

M. Mukhlisin and M. R. Taha, “Slope stability analysis of a weathered granitic hillslope as effects of soil thickness,” European Journal of Scientific Research, vol. 30, no. 1, pp. 36–44, 2009.

M. Mukhlisin, M. R. Baidillah, M. R. Taha, and A. El-Shafie, “Effect of soil water retention model on slope stability analysis,” International Journal of Physical Sciences, vol. 6, no. 19, pp. 4629–4635, 2011.

M. Mukhlisin, I. Idris, W. Z. Wan Yaacob, A. ElShafie, and M. R. Taha, “Soil slope deformation behavior in relation to soil water interaction based on centrifuge physical modeling,” International Journal of Physical Sciences, vol. 6, no. 13, pp. 3126–3133, 2011.

M. Mukhlisin, A. El-Shafie, and M. R. Taha, “Regularized versus non-regularized neural network model for prediction of saturated soil-water content on weathered granite soil formation,” Neural Computing and Applications, vol. 21, no. 3, pp. 543–553, 2012.

M. Mukhlisin and M. R. Taha, “Numerical model of antecedent rainfall effect on slope stability at a hillslope of weathered granitic soil formation,” Journal Geological Society of India, vol. 79, no. 5, pp. 525–531, 2012.

M. Mukhlisin, A. Saputra, A. El-Shafie, and M. R. Taha, “Measurement of dynamic soil water content based on electrochemical capacitance tomography,” International Journal of Electrochemical Science, vol. 7, pp. 5467–5483, 2012.

R. Černý, “Time-domain reflectometry method and its application for measuring moisture content in porous materials: a review,” Measurement, vol. 42, no. 3, pp. 329–336, 2009.

S. R. Evett and G. W. Parkin, “Advances in soil water content sensing: the continuing maturation of technology and theory,” Vadose Zone Journal, vol. 4, no. 4, pp. 986–991, 2005.

J. A. Huisman, S. S. Hubbard, J. D. Redman, and A. P. Annan, “Measuring soil water content with ground penetrating radar: a review,” Vadose Zone Journal, vol. 2, no. 4, pp. 476–491, 2003.

K. Noborio, “Measurement of soil water content and electrical conductivity by time domain reflectometry: a review,” Computers and Electronics in Agriculture, vol. 31, no. 3, pp. 213–237, 2001.

D. A. Robinson, S. B. Jones, J. M. Blonquist Jr., and S. P. Friedman, “A physically derived water content/permittivity calibration model for coarse-textured, layered soils,” Soil Science Society of America Journal, vol. 69, no. 5, pp. 1372–1378, 2005.

D. A. Robinson, C. S. Campbell, J. W. Hopmans et al., “Soil moisture measurement for ecological and hydrological watershed-scale observatories: a review,” Vadose Zone Journal, vol. 7, no. 1, pp. 358–389, 2008.

J. Curtis and R. Narayanan, “Effects of laboratory procedures on soil electrical property measurements,” IEEE Transactions on Instrumentation and Measurement, vol. 47, no. 6, pp. 1474–1480, 1998.

T. Saarenketo, “Electrical properties of water in clay and silty soils,” Journal of Applied Geophysics, vol. 40, no. 1–3, pp. 73–88, 1998.

J. O. Curtis, “Moisture effects on the dielectric properties of soils,” IEEE Transactions on Geoscience and Remote Sensing, vol. 39, no. 1, pp. 125–128, 2001.

C. M. K. Gardner, T. J. Dean, and J. D. Cooper, “Soil water content measurement with a high-frequency capacitance sensor,” Journal of Agricultural Engineering Research, vol. 71, no. 4, pp. 395–403, 1998.

S. Y. Wu, Q. Y. Zhou, G. Wang, L. Yang, and C. P. Ling, “The relationship between electrical capacitance-based dielectric constant and soil water content,” Environmental Earth Sciences, vol. 62, no. 5, pp. 999–1011, 2011.

C. H. Roth, M. A. Malicki, and R. Plagge, “Empirical evaluation of the relationship between soil dielectric constant and volumetric water content as the basis for calibrating soil moisture measurements by TDR,” Journal of Soil Science, vol. 43, no. 1, pp. 1–13, 1992.

G. C. Topp, J. L. Davis, and A. P. Annan, “Electromagnetic determination of soil water content: measurements in coaxial transmission lines,” Water Resources Research, vol. 16, no. 3, pp. 574–582, 1980.

M. A. Malicki, R. Plagge, and C. H. Roth, “Improving the calibration of dielectric TDR soil moisture determination taking into account the solid soil,” European Journal of Soil Science, vol. 47, no. 3, pp. 357–366, 1996.

M. C. Dobson, F. T. Ulaby, M. T. Hallikainen, and M. A. El-Rayes, “Microwave dielectric behavior of wet soil—part II: dielectric mixing models,” IEEE Transactions on Geoscience and Remote Sensing, vol. 23, no. 1, pp. 35–46, 1985.

P. A. Ferré, D. L. Rudolph, and R. G. Kachanoski, “Spatial averaging of water content by time domain reflectometry: implications for twin rod probes with and without dielectric coatings,” Water Resources Research, vol. 32, no. 2, pp. 271–279, 1996.

S. P. Friedman, “A saturation degree-dependent composite spheres model for describing the effective dielectric constant of unsaturated porous media,” Water Resources Research, vol. 34, no. 11, pp. 2949–2961, 1998.

M. A. Hilhorst, C. Dirksen, F. W. H. Kampers, and R. A. Feddes, “New dielectric mixture equation for porous materials based on depolarization factors,” Soil Science Society of America Journal, vol. 64, no. 5, pp. 1581–1587, 2000.

K. Roth, R. Schulin, H. Fluhler, and W. Attinger, “Calibration of time domain reflectometry for water content measurement using a composite dielectric approach,” Water Resources Research, vol. 26, no. 10, pp. 2267–2273, 1990.

M. G. Schaap, L. De Lange, and T. J. Heimovaara, “TDR calibration of organic forest floor media,” Soil Technology, vol. 11, no. 2, pp. 205–217, 1997.

J. R. Wang and T. J. Schmugge, “An empirical model for the complex dielectric permittivity of soils as a function of water content,” IEEE Transactions on Geoscience and Remote Sensing, vol. 18, no. 4, pp. 288–295, 1980.

P. Sabouroux and D. Ba, “Epsimu, a tool for dielectric properties measurement of porous media: application in wet granular materials characterization,” Progress In Electromagnetics Research B, no. 29, pp. 191–207, 2011.

W. Skierucha, A. Wilczek, and O. Alokhina, “Calibration of a TDR probe for low soil water content measurements,” Sensors and Actuators A, vol. 147, no. 2, pp. 544–552, 2008.

M. A. Malicki, R. Plagge, M. Renger, and R. T. Walczak, “Application of time-domain reflectometry (TDR) soil moisture miniprobe for the determination of unsaturated soil water characteristics from undisturbed soil cores,” Irrigation Science, vol. 13, no. 2, pp. 65–72, 1992.

J. R. Lundien, “Terrain analysis by electromagnetic means,” U.S. Army Engineer Waterways Experiment Station, 1971.

R. W. Newton, Terrain Analysis by Electromagnetic Means, Texas A&M University, College Station, Tex, USA, 1977.

J. Wang, T. Schmugge, and D. Williams, “Dielectric constants of soils at microwave frequencies—II,” National Aeronautics and Space Administration, 1978.

S. Sharif, “Chemical and mineral composition of dust and its effect on the dielectric constant,” IEEE Transactions on Geoscience and Remote Sensing, vol. 33, no. 2, pp. 353–359, 1995.