Performance Evaluation of Volumetric Water Content and Relative Permittivity Models

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 R 2 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.


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][3][4][5][6][7][8][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][11][12][13][14]). Moreover, issues of the instruments used for measurement of water content, from small-scale (<1 m 2 ) to large-scale (100 m 2 ), 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][17][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]).
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.

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.

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:  Roth et al. [21] used miniprobe TDR for their experiment to propose another empirical equation for therelationship, which had been used previously by [33]. Therelationship for mineral soil proposed by [21] is = −0.0728 + 0.0448 − 0.00195 2 + 0.0000361 3 , (2a) 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 3 cm −3 , respectively [21].
Ferré et al. [25] proposed a simple equation for therelationship. 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 therelationship, 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 = 33, = 2, = 1.5, and = 14.

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][35][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., = 3.2, = 80, and = 1), 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 3 cm −3 [28], with forms of the following equation: The Scientific World Journal 3 where = −1 for three phases in series and = 1 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 [ 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 dry and sat are the permittivity values for dry and saturated soil, respectively.

4
The Scientific World Journal  Table 1 shows the 11 proposed equations of therelationship 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.
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.

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.
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].

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 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. 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),   Equation (11) (f) Figure 3: Comparison of (7a) to (11) with all data and different porosity.

6
The Scientific World Journal porosity (<0.5). On the other hand, the larger porosity (>0.6) is occupied by (4). 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 thesquare 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. Figure 3 shows the effect of porosity ( = 0.3 to = 0.7) 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.

Model with Two Parameters.
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.

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.