**Abstract**

Electromagnetic methods have been widely used in the measurement of the water content of the soil. These methods utilize the permittivity as electrical properties of the soil, to determine the moisture content of the soil. Since the measurements are carried out indirectly, a calibration between permittivity and the water content of the soil is needed. Generally, the calibration method is generated by using an empirical and mixing model. This study presents a proposed model of calibration by using a normalization approach to calibrate the value of the permittivity of the water content of the soil. Then the model was applied using electrical capacitance volume tomography (ECVT) to image soil water content during infiltration of water in a soil column. Granular and silty sand were used as soil material in the experiments. The result showed that the model for measuring moisture water content can be seen in each layer during soil water infiltration in the soil column.

**Keywords:** dielectric analysis model, soil moisture water content, ECVT

#### **1. Introduction**

Measurement of soil moisture water content has become an important part of the analysis of various fields of study, especially those involving irrigation in agriculture, forestry, hydrology, and land activity. For example, in agriculture study, it is required to determine water source and ensure the quality of crop [1]. Moreover, soil moisture water content is useful to analyze soil water contamination by observing changes in water content during the addition of substance [2]. Soil water content also plays an important role in slope stability analysis [3–6].

Various techniques in the measurement of contamination and water content of soil have been discussed in literatures (e.g., [7–12]). Based on previous studies, soil water content measurement techniques are widely used as an electromagnetic method [13] such as time domain reflectometry [14], ground-penetrating radar (GPR) [8], and electrical capacitance [15]. These all methods measure the value of relative permittivity of soil to find soil moisture water content.

Tomography is a promising technique for measurement of water content in soil, especially for capacitance-based tomography. It is because this technique is not only capable of measuring water content in the soil but also is capable of imaging the

distribution of water in the soil. Tomography technique is also preferred because it is nondestructive and noninvasive. As a tomography technique, ECVT is a system used to view enclosed objects by measuring changes in capacitance and then compute relative permittivity distribution to create three-dimensional images in real time [16]. The shape of geometry sensor on ECVT is not confined to one form; it can be in the form of arbitrary shape of geometries [17]. This possibility gives an extra advantage of ECVT in measuring soil water content.

The previous study has successfully monitored the propagation of distribution of water in soil column [18]. In this study, the equation of normalized volumetric water content from ECVT method is proposed and compared with other equations. The proposed model is then used to analyze the volumetric water content during soil water infiltration in a vessel.

### **2. ECVT principle**

The basic measurement of ECVT is derived from Poisson's equation:

$$\nabla \cdot \left( \varepsilon(\mathbf{x}, \mathbf{y}, \mathbf{z}) \nabla \rho(\mathbf{x}, \mathbf{y}, \mathbf{z}) = -\rho(\mathbf{x}, \mathbf{y}, \mathbf{z}) \right) \tag{1}$$

where *ε* is relative permittivity distribution, *ϕ* is electric potential, and *ρ* is charge distribution. From Eq. (1), capacitance value can be obtained by using the equation below:

$$C = \frac{1}{\Delta V} \oint\_{\Gamma} \varepsilon(\mathbf{x}, \mathbf{y}, z) \nabla \rho(\mathbf{x}, \mathbf{y}, z) \cdot \hat{n} \, dl \tag{2}$$

where Δ*V* is potential difference and *C* is capacitance. By using matrix expression, Eq. (2) can be written like the following equation:

$$\mathbf{C} = \mathbf{S} \cdot \mathbf{G} \tag{3}$$

model across a range of soil water content. This model can be divided into two categories, which are model with one parameter and model with two or more

*Dielectric Analysis Model for Measurement of Soil Moisture Water Content Using Electrical…*

There are some models that proposed the ε-θ relationship. Topp et al. [20] introduced successfully the ε-θ relationship that is commonly used in geotechnical

where *ε* is the relative permittivity or dielectric constant and *θ* is the volumetric

Equation (1) is derived empirically through experiments of various of mineral soil using a time-domain reflectometer (TDR) at a frequency between 1 MHz and 1 GHz, with an estimated error value of 0.013 [20]. In another form, Topp's equa-

In addition, Topp et al. [20] also proposed equations of ε-θ relationship for

Then Roth et al. [21] proposed another empirical equation by an experiment using miniprobe TDR which has been used previously by [22]. The ε-θ relationship

Calibration of these equations has a volumetric water content error value of

Simple equation of ε-θ relationship was proposed by [23]. This equation resulted

Schaap et al. [24] also introduce a simple equation by performing experiments of

The next equation comes from [25], which is using coaxial transmission system

0.015 cm<sup>3</sup> cm�<sup>3</sup> for mineral soil and 0.035 cm3 cm�<sup>3</sup> for organic soil [22].

*<sup>θ</sup>* <sup>¼</sup> <sup>0</sup>*:*<sup>1181</sup> ffiffi

505 measurements of organic forest floor sample by using TDR where the ε-θ

*<sup>θ</sup>* <sup>¼</sup> <sup>0</sup>*:*<sup>136</sup> ffiffi

and using soil samples with wide range of soil textures:

from the principal of dielectric mixing models and analyzed the TDR without coatings that are considered potential sources of error in measurement:

while the ε-θ relationship for organic soil and material is

*<sup>ε</sup>* <sup>¼</sup> <sup>1</sup>*:*<sup>74</sup> � <sup>0</sup>*:*34*<sup>θ</sup>* <sup>þ</sup> <sup>135</sup>*θ*<sup>2</sup> � <sup>55</sup>*:*3*θ*<sup>3</sup> Organic soil (7)

*<sup>ε</sup>* <sup>¼</sup> <sup>3</sup>*:*<sup>57</sup> <sup>þ</sup> <sup>31</sup>*:*7*<sup>θ</sup>* <sup>þ</sup> <sup>114</sup>*θ*<sup>2</sup> � <sup>68</sup>*:*2*θ*<sup>3</sup> <sup>450</sup> � *<sup>μ</sup>*m glass beads (8)

*<sup>θ</sup>* ¼ �0*:*<sup>0728</sup> <sup>þ</sup> <sup>0</sup>*:*0448*<sup>ε</sup>* � <sup>0</sup>*:*<sup>00195</sup> *<sup>ε</sup>*<sup>2</sup> <sup>þ</sup> <sup>0</sup>*:*000036 1*ε*<sup>3</sup> (9)

*<sup>θ</sup>* <sup>¼</sup> ‐0*:*<sup>0233</sup> <sup>þ</sup> <sup>0</sup>*:*0285*<sup>ε</sup>* � <sup>0</sup>*:*<sup>000431</sup> *<sup>ε</sup>*<sup>2</sup> <sup>þ</sup> <sup>0</sup>*:*00000304*ε*<sup>3</sup> (10)

*<sup>ε</sup>* <sup>¼</sup> <sup>3</sup>*:*<sup>03</sup> <sup>þ</sup> <sup>9</sup>*:*3*<sup>θ</sup>* <sup>þ</sup> <sup>146</sup>*:*0*θ*<sup>2</sup> � <sup>76</sup>*:*7*θ*<sup>3</sup> (5)

*<sup>ε</sup>* � <sup>5</sup>*:*<sup>5</sup> � <sup>10</sup>�<sup>4</sup>*ε*<sup>2</sup> <sup>þ</sup> <sup>4</sup>*:*<sup>3</sup> � <sup>10</sup>�<sup>6</sup>*ε*<sup>3</sup> (6)

*<sup>ε</sup>* <sup>p</sup> � <sup>0</sup>*:*<sup>1841</sup> (11)

*<sup>ε</sup>* <sup>p</sup> � <sup>0</sup>*:*<sup>119</sup> (12)

parameters.

**3.1 Model with one parameter**

water content of soil.

area. The relationship is as shown below:

*DOI: http://dx.doi.org/10.5772/intechopen.89057*

tion can also be written as follows:

organic soil and 450 μm glass beads.

for mineral soil proposed by [21] is

relationship is

**75**

*<sup>θ</sup>* ¼ �5*:*<sup>3</sup> � <sup>10</sup>�<sup>2</sup> <sup>þ</sup> <sup>2</sup>*:*<sup>92</sup> � <sup>10</sup>�<sup>2</sup>

where **C** is capacitance matrix, **G** is distribution of relative permittivity matrix, and **S** is sensitivity matrix. The sensitivity matrix is generated from the sensor and geometry design and number of sensors. In matrix operation, the value of **G** can be obtained by inversing matrix **S** and multiplying it with matrix **C**. For non-square matrix, matrix inversion is very difficult to solve, so the approximation could be attempted by using transpose matrix. The equation to calculate matrix **G** becomes

$$\mathbf{G} = \mathbf{S}^{\mathrm{T}} \; \mathbf{C} \tag{4}$$

Equations (3) and (4) are known as the forward and inverse problem, respectively. Inverse problem is used to reconstruct the capacitance measurements to become relative permittivity distribution. The simple method for reconstruction is using linear back projection (LBP) [19].

#### **3. Models of soil moisture water content and relative permittivity relationship**

The relationship between relative permittivity and volumetric water content has been used by previous researchers to determine the volumetric water content. Many numbers of functions have been proposed to describe the ε-θ relationship

*Dielectric Analysis Model for Measurement of Soil Moisture Water Content Using Electrical… DOI: http://dx.doi.org/10.5772/intechopen.89057*

model across a range of soil water content. This model can be divided into two categories, which are model with one parameter and model with two or more parameters.

#### **3.1 Model with one parameter**

distribution of water in the soil. Tomography technique is also preferred because it is nondestructive and noninvasive. As a tomography technique, ECVT is a system used to view enclosed objects by measuring changes in capacitance and then compute relative permittivity distribution to create three-dimensional images in real time [16]. The shape of geometry sensor on ECVT is not confined to one form; it can be in the form of arbitrary shape of geometries [17]. This possibility gives an

The previous study has successfully monitored the propagation of distribution of

∇ � ð Þ *ε*ð Þ *x, y, z* ∇*φ*ð Þ¼� *x, y, z ρ*ð Þ *x, y, z* (1)

*ε*ð Þ *x, y, z* ∇*φ*ð Þ� *x, y, z n dl* ^ (2)

**C** ¼ **S G** (3)

**<sup>G</sup>** <sup>¼</sup> **ST <sup>C</sup>** (4)

water in soil column [18]. In this study, the equation of normalized volumetric water content from ECVT method is proposed and compared with other equations. The proposed model is then used to analyze the volumetric water content during

The basic measurement of ECVT is derived from Poisson's equation:

where *ε* is relative permittivity distribution, *ϕ* is electric potential, and *ρ* is charge distribution. From Eq. (1), capacitance value can be obtained by using the

where Δ*V* is potential difference and *C* is capacitance. By using matrix expres-

where **C** is capacitance matrix, **G** is distribution of relative permittivity matrix, and **S** is sensitivity matrix. The sensitivity matrix is generated from the sensor and geometry design and number of sensors. In matrix operation, the value of **G** can be obtained by inversing matrix **S** and multiplying it with matrix **C**. For non-square matrix, matrix inversion is very difficult to solve, so the approximation could be attempted by using transpose matrix. The equation to calculate matrix **G** becomes

Equations (3) and (4) are known as the forward and inverse problem, respectively. Inverse problem is used to reconstruct the capacitance measurements to become relative permittivity distribution. The simple method for reconstruction is

The relationship between relative permittivity and volumetric water content has

**3. Models of soil moisture water content and relative permittivity**

been used by previous researchers to determine the volumetric water content. Many numbers of functions have been proposed to describe the ε-θ relationship

extra advantage of ECVT in measuring soil water content.

*Modern Applications of Electrostatics and Dielectrics*

*<sup>C</sup>* <sup>¼</sup> <sup>1</sup> Δ*V* þ Γ

sion, Eq. (2) can be written like the following equation:

using linear back projection (LBP) [19].

**relationship**

**74**

soil water infiltration in a vessel.

**2. ECVT principle**

equation below:

There are some models that proposed the ε-θ relationship. Topp et al. [20] introduced successfully the ε-θ relationship that is commonly used in geotechnical area. The relationship is as shown below:

$$
\epsilon = \mathbf{3.03} + \mathbf{9.39} + \mathbf{146.09}^2 - \mathbf{76.70}^3 \tag{5}
$$

where *ε* is the relative permittivity or dielectric constant and *θ* is the volumetric water content of soil.

Equation (1) is derived empirically through experiments of various of mineral soil using a time-domain reflectometer (TDR) at a frequency between 1 MHz and 1 GHz, with an estimated error value of 0.013 [20]. In another form, Topp's equation can also be written as follows:

$$\theta = -5.3 \times 10^{-2} + 2.92 \times 10^{-2} \text{e} - 5.5 \times 10^{-4} \text{e}^2 + 4.3 \times 10^{-6} \text{e}^3 \tag{6}$$

In addition, Topp et al. [20] also proposed equations of ε-θ relationship for organic soil and 450 μm glass beads.

$$\varepsilon = 1.74 - 0.34\theta + 135\theta^2 - 55.3\theta^3 \quad \text{Organic soil} \tag{7}$$

$$\boldsymbol{\alpha} = \mathbf{3.57} + \mathbf{31.70} + \mathbf{1140}^2 - \mathbf{68.2} \boldsymbol{\theta}^3 \quad \mathbf{450} - \mu \mathbf{m} \text{ glass leads} \tag{8}$$

Then Roth et al. [21] proposed another empirical equation by an experiment using miniprobe TDR which has been used previously by [22]. The ε-θ relationship for mineral soil proposed by [21] is

$$\theta = -0.0728 + 0.0448\varepsilon - 0.00195\,\varepsilon^2 + 0.000036\,\text{1}\varepsilon^3 \tag{9}$$

while the ε-θ relationship for organic soil and material is

$$\theta = \text{-0.0233} + \text{0.0285}\epsilon - \text{0.000431}\,\epsilon^2 + \text{0.000000304}\epsilon^3 \tag{10}$$

Calibration of these equations has a volumetric water content error value of 0.015 cm<sup>3</sup> cm�<sup>3</sup> for mineral soil and 0.035 cm3 cm�<sup>3</sup> for organic soil [22].

Simple equation of ε-θ relationship was proposed by [23]. This equation resulted from the principal of dielectric mixing models and analyzed the TDR without coatings that are considered potential sources of error in measurement:

$$
\theta = 0.1181 \sqrt{\varepsilon} - 0.1841 \tag{11}
$$

Schaap et al. [24] also introduce a simple equation by performing experiments of 505 measurements of organic forest floor sample by using TDR where the ε-θ relationship is

$$
\theta = 0.136\sqrt{\varepsilon} - 0.119 \tag{12}
$$

The next equation comes from [25], which is using coaxial transmission system and using soil samples with wide range of soil textures:

$$\theta = -0.0286 + 0.02435e - 0.0003421e^2 + 0.00000237e^3 \tag{13}$$

Equations (5)–(8) are used for measurement at a frequency of 100 MHz.

[13] proposed an empirical model where permittivity measurement was measured based on capacitance. This experiment used a type of quartz sand with a range of particle sizes between 0.15 and 0.9 mm:

$$\varepsilon = A \left( \frac{1}{1 + (a(1 - \theta))^n} \right)^{1 - \frac{1}{n}} + B \tag{14}$$

using multiple linear regression analysis to best fit the measurement data, resulting

*Dielectric Analysis Model for Measurement of Soil Moisture Water Content Using Electrical…*

*<sup>ε</sup>* <sup>p</sup> <sup>þ</sup> <sup>1</sup>*:*<sup>208</sup> � <sup>2</sup>*:*454*<sup>ρ</sup>*

Robinson et al. [14] give the equation used for coarse-textured, layered soils by

!

ffiffi *<sup>ε</sup>* <sup>p</sup> � ffiffiffiffiffiffiffi *εdry* p

ffiffiffiffiffiffi *<sup>ε</sup>sat* <sup>p</sup> � ffiffiffiffiffiffiffi *εdry* p

Data from previous research (e.g., [20, 21, 25–27, 29, 31–35]) are used to compare the patterns of the equations which are discussed in this study. This data

**Figure 1** shows several curves representing Eqs. (5)–(14). All equations look occupied by all the available data. However, each equation appears to have certain characteristics to the data. Eqs. (8) and (12) only cover the boundary area of the data, while the other equations lie mostly in the central part of the data. These would seem to depend on the properties of the soil types used. One of the charac-

**Figure 2** shows the influence of porosity (η = 0.3 to η = 0.7), on the suitability of the equations (Eqs. (15), (16), (18) and (19)–(21)) with data, and also displayed some of the data with a value of porosity (0.33, 0.44, and 0.62) in order to see the effect of porosity on predictions of water content of the equation. From the image it can be seen that the different porosity values of the data will result in different patterns. In **Figure 2a**, the equations already have the same pattern with the data, but a change of porosity in the equation does not give a significant effect on the pattern of the line, so it is only fit in certain small area of data though with different porosity. **Figure 2b** shows that the equation is such as the linear equation that has not affected on the changes in porosity. These equations also appear not to follow the pattern of distribution of data. The same thing happened in **Figure 2c** with a shift in

An overestimated result is produced in **Figure 2e**, where the equation is not able to cover all areas of data; whereas in **Figure 2f**, it has been seen covering almost all areas except the data on water content values smaller than 0.2, but this equation has not been able to adjust to the data that have a value of porosity. The only equation that gives the better approach is **Figure 2d**. This equation produces a pattern in accordance with the existing data. The equation is also seen fit to data that has a

In this study the relative permittivity was analyzed by the ECVT system generated in the form of normalization. Normalized volumetric water content can be

> <sup>Θ</sup> <sup>¼</sup> *<sup>θ</sup>* � *<sup>θ</sup><sup>r</sup> θ<sup>s</sup>* � *θ<sup>r</sup>*

where *εdry* and *εsat* are the permittivity measured at oven dry soil and saturation

<sup>9</sup>*:*<sup>93</sup> (20)

(21)

(22)

*θ* ¼

using TDR and coarse-grained, glass bead, and quartz grains:

consists of various soil types with different properties.

teristics analyzed in this study is the porosity of the soil.

values on the x-axis which is more to the left.

value of porosity.

defined as

**77**

*θ* ¼ *η*

ffiffi

in the following equation:

soil (**Table 1**).

where *ρ* is dry bulk density.

*DOI: http://dx.doi.org/10.5772/intechopen.89057*

**3.3 Comparison using existing data**

where A = 33, B = 2, α = 1.5, and n = 14.

#### **3.2 Model with two or more parameters**

Some relationship equations between 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 others study [26–28], [29] proposed the following equations:

$$
\varepsilon = \theta \left( \varepsilon\_i + (\varepsilon\_w - \varepsilon\_i) \frac{\theta}{\theta\_t} \gamma \right) + (\eta - \theta) \varepsilon\_d + (\mathbf{1} - \eta) \varepsilon\_r \tag{15}
$$

Equation (11) is used for *θ* ≤ *θt*, while *θ* > *θ<sup>t</sup>* used the following equation:

$$\varepsilon = \theta\_l(\varepsilon\_i + (\varepsilon\_w - \varepsilon\_i)\eta) + (\theta - \theta\_l)\varepsilon\_w + (\eta - \theta)\varepsilon\_d + (1 - \eta)\varepsilon\_r \tag{16}$$

where *εi*, *εw*, *εa*, and *ε<sup>r</sup>* are the permittivity of ice, water, air and rock, respectively (i.e., *ε<sup>i</sup>* = 3.2, *ε<sup>w</sup>* = 80, and *ε<sup>a</sup>* = 1), while *θ<sup>t</sup>* is transition moisture (0.16–0.33), *η* is the porosity of soil (0.5), and γ is the fitting parameter (0.3–0.5) [29].

In [30], the equation based on dielectric mixing model, which has been described by [31], is proposed. Experiments carried out by measuring a wide range of soil types using TDR with the error value of soil water content is not more than 0.013 cm3 cm�<sup>3</sup> [30], with forms of the equation below:

$$\theta = \frac{\varepsilon^{\mathbb{Y}} - (\mathbf{1} - \boldsymbol{\eta})\varepsilon\_{\mathfrak{s}}^{\mathbb{Y}} - \eta \varepsilon\_{\mathfrak{a}}^{\mathbb{Y}}}{\varepsilon\_{\mathfrak{w}}^{\mathbb{Y}} - \varepsilon\_{\mathfrak{a}}^{\mathbb{Y}}} \quad ; \boldsymbol{\gamma} = -\mathbf{1} \tag{17}$$

$$\theta = \frac{\varepsilon^{\mathbb{Y}} - (\mathbb{1} - \eta)\varepsilon\_{\varepsilon}^{\mathbb{Y}} - \eta \varepsilon\_{a}^{\mathbb{Y}}}{\varepsilon\_{w}^{\mathbb{Y}} - \varepsilon\_{a}^{\mathbb{Y}}} \quad ; \gamma = \mathbb{1} \tag{18}$$

where γ = �1 for three phases in series and γ = 1 for three phases in parallel.

Another model is proposed by [32]. They conducted the experiments by using TDR and 62 kinds of soil sample that consist of mineral soils, organic soil, standard pot soils, artificial peat loess and peat sand, sea and river sand, forest litter, etc. which differ in terms of texture and bulk density [32].

$$\theta = \frac{\sqrt{\varepsilon} - 3.47 + 6.22\eta - 3.82\eta^2}{7.01 + 6.89\eta - 7.83\eta^2} \tag{19}$$

Equation (13) gives the uncertainty of soil water content value of 0.03.

Gardner et al. [15] used capacitance measurement methods to obtain soil water content with the soil dry bulk density values ranging between 1.08 and 1.49, then

*Dielectric Analysis Model for Measurement of Soil Moisture Water Content Using Electrical… DOI: http://dx.doi.org/10.5772/intechopen.89057*

using multiple linear regression analysis to best fit the measurement data, resulting in the following equation:

$$\theta = \frac{\sqrt{\varepsilon} + 1.208 - 2.454\rho}{9.93} \tag{20}$$

where *ρ* is dry bulk density.

*<sup>θ</sup>* ¼ �0*:*<sup>0286</sup> <sup>þ</sup> <sup>0</sup>*:*02435*<sup>ε</sup>* � <sup>0</sup>*:*0003421*ε*<sup>2</sup> <sup>þ</sup> <sup>0</sup>*:*00000237*ε*<sup>3</sup> (13)

*n*

þ *B* (14)

þ ð Þ *η* � *θ ε<sup>a</sup>* þ ð Þ 1 � *η ε<sup>r</sup>* (15)

*; γ* ¼ �1 (17)

*; γ* ¼ 1 (18)

Equations (5)–(8) are used for measurement at a frequency of 100 MHz. [13] proposed an empirical model where permittivity measurement was measured based on capacitance. This experiment used a type of quartz sand with a

> <sup>1</sup> <sup>þ</sup> ð Þ *<sup>α</sup>*ð Þ <sup>1</sup> � *<sup>θ</sup> <sup>n</sup>* � �1�<sup>1</sup>

Some relationship equations between 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 others study [26–28], [29] proposed

> *θ θt γ*

Equation (11) is used for *θ* ≤ *θt*, while *θ* > *θ<sup>t</sup>* used the following equation:

is the porosity of soil (0.5), and γ is the fitting parameter (0.3–0.5) [29]. In [30], the equation based on dielectric mixing model, which has been described by [31], is proposed. Experiments carried out by measuring a wide range of soil types using TDR with the error value of soil water content is not more than

*<sup>θ</sup>* <sup>¼</sup> *εγ* � ð Þ <sup>1</sup> � *<sup>η</sup> <sup>ε</sup><sup>γ</sup>*

*ε γ <sup>w</sup>* � *ε γ a*

*<sup>θ</sup>* <sup>¼</sup> *<sup>ε</sup><sup>γ</sup>* � ð Þ <sup>1</sup> � *<sup>η</sup> εγ*

*ε γ <sup>w</sup>* � *ε γ a*

*ε* ¼ *θt*ð*ε<sup>i</sup>* þ ð Þ *ε<sup>w</sup>* � *ε<sup>i</sup> γ*Þ þ ð Þ *θ* � *θ<sup>t</sup> ε<sup>w</sup>* þ ð Þ *η* � *θ ε<sup>a</sup>* þ ð Þ 1 � *η ε<sup>r</sup>* (16)

*<sup>s</sup>* � *ηεγ a*

> *<sup>s</sup>* � *ηεγ a*

where γ = �1 for three phases in series and γ = 1 for three phases in parallel. Another model is proposed by [32]. They conducted the experiments by using TDR and 62 kinds of soil sample that consist of mineral soils, organic soil, standard pot soils, artificial peat loess and peat sand, sea and river sand, forest litter, etc.

*<sup>ε</sup>* <sup>p</sup> � <sup>3</sup>*:*<sup>47</sup> <sup>þ</sup> <sup>6</sup>*:*22*<sup>η</sup>* � <sup>3</sup>*:*82*η*<sup>2</sup>

Gardner et al. [15] used capacitance measurement methods to obtain soil water content with the soil dry bulk density values ranging between 1.08 and 1.49, then

Equation (13) gives the uncertainty of soil water content value of 0.03.

<sup>7</sup>*:*<sup>01</sup> <sup>þ</sup> <sup>6</sup>*:*89*<sup>η</sup>* � <sup>7</sup>*:*83*η*<sup>2</sup> (19)

where *εi*, *εw*, *εa*, and *ε<sup>r</sup>* are the permittivity of ice, water, air and rock, respectively (i.e., *ε<sup>i</sup>* = 3.2, *ε<sup>w</sup>* = 80, and *ε<sup>a</sup>* = 1), while *θ<sup>t</sup>* is transition moisture (0.16–0.33), *η*

*<sup>ε</sup>* <sup>¼</sup> *<sup>A</sup>* <sup>1</sup>

range of particle sizes between 0.15 and 0.9 mm:

*Modern Applications of Electrostatics and Dielectrics*

where A = 33, B = 2, α = 1.5, and n = 14.

*ε* ¼ *θ ε<sup>i</sup>* þ ð Þ *ε<sup>w</sup>* � *ε<sup>i</sup>*

0.013 cm3 cm�<sup>3</sup> [30], with forms of the equation below:

which differ in terms of texture and bulk density [32].

*θ* ¼

**76**

ffiffi

� �

**3.2 Model with two or more parameters**

the following equations:

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

$$\theta = \eta \left( \frac{\sqrt{\varepsilon} - \sqrt{\varepsilon\_{dry}}}{\sqrt{\varepsilon\_{sat}} - \sqrt{\varepsilon\_{dry}}} \right) \tag{21}$$

where *εdry* and *εsat* are the permittivity measured at oven dry soil and saturation soil (**Table 1**).

#### **3.3 Comparison using existing data**

Data from previous research (e.g., [20, 21, 25–27, 29, 31–35]) are used to compare the patterns of the equations which are discussed in this study. This data consists of various soil types with different properties.

**Figure 1** shows several curves representing Eqs. (5)–(14). All equations look occupied by all the available data. However, each equation appears to have certain characteristics to the data. Eqs. (8) and (12) only cover the boundary area of the data, while the other equations lie mostly in the central part of the data. These would seem to depend on the properties of the soil types used. One of the characteristics analyzed in this study is the porosity of the soil.

**Figure 2** shows the influence of porosity (η = 0.3 to η = 0.7), on the suitability of the equations (Eqs. (15), (16), (18) and (19)–(21)) with data, and also displayed some of the data with a value of porosity (0.33, 0.44, and 0.62) in order to see the effect of porosity on predictions of water content of the equation. From the image it can be seen that the different porosity values of the data will result in different patterns.

In **Figure 2a**, the equations already have the same pattern with the data, but a change of porosity in the equation does not give a significant effect on the pattern of the line, so it is only fit in certain small area of data though with different porosity.

**Figure 2b** shows that the equation is such as the linear equation that has not affected on the changes in porosity. These equations also appear not to follow the pattern of distribution of data. The same thing happened in **Figure 2c** with a shift in values on the x-axis which is more to the left.

An overestimated result is produced in **Figure 2e**, where the equation is not able to cover all areas of data; whereas in **Figure 2f**, it has been seen covering almost all areas except the data on water content values smaller than 0.2, but this equation has not been able to adjust to the data that have a value of porosity. The only equation that gives the better approach is **Figure 2d**. This equation produces a pattern in accordance with the existing data. The equation is also seen fit to data that has a value of porosity.

In this study the relative permittivity was analyzed by the ECVT system generated in the form of normalization. Normalized volumetric water content can be defined as

$$
\Theta = \frac{\theta - \theta\_r}{\theta\_s - \theta\_r} \tag{22}
$$


**Figure 1.**

**Figure 2.**

**79**

*e) Eq. 20; f) Eq. 21.*

*Comparison all equations with different porosity, by using: a) Eq. 15; b) Eq. 16; c) Eq. 18; d) Eq. 19;*

*Dielectric Analysis Model for Measurement of Soil Moisture Water Content Using Electrical…*

*Comparisons using all data for Eqs. (5)–(14).*

*DOI: http://dx.doi.org/10.5772/intechopen.89057*

**Table 1.**

*Summary of all equations of ε-θ relationship.*

*Dielectric Analysis Model for Measurement of Soil Moisture Water Content Using Electrical… DOI: http://dx.doi.org/10.5772/intechopen.89057*

**Figure 1.** *Comparisons using all data for Eqs. (5)–(14).*

**Eq. Source Experimental method Soil type Properties of soil**

• Mineral soil • Organic soil • Vermiculite • Glass beads

(7) Organic soil — 0.422 —

beads

9 mineral soils

(10) 7 organic soils 0.527–0.785 0.2–0.77 0.70–

25 samples of forest floors

22 different samples

different field

soil samples

• Brown earths • Silica materials

grained, quartz grain, sandy soil

Model with one parameter

(6)

(5) [20] ε: using TDR Tektronix

(9) [21] ε: TDR miniprobe 250 ps

(11) [23] Using model of inverse

(12) [24] ε: TDR Tektronix 1502B

(13) [25] Coaxial transmission/

(14) [13] ε: based on capacitance measurement

Model with two or more parameters (15) [29] Modeling by using

[26–28]

(20) [15] Capacitance probe 80– 150 MHz

*Summary of all equations of ε-θ relationship.*

(16)

**Table 1.**

**78**

7S12 model to perform 18 experiments with different treatments θ: using gravimetric technique

*Modern Applications of Electrostatics and Dielectrics*

(8) 450 μm glass

rise-time needle pulse θ: gravimetric technique

averaging for TDR method by analyzing the

θ: gravimetric technique

θ: gravimetric technique

data from other studies

(17) [30] TDR From 11

(19) [32] TDR CAMI 62 kinds of

(21) [14] TDR Tektronix 1502B Coarse-

sites (18)

reflection apparatus controlled by a Hewlett-Packard 8510C Vector Network Analyzer system 45 MHz to 26.5 GHz

mixing model

**Porosity (cm<sup>3</sup> cm<sup>3</sup>**

0.418– 0.482

— —— —

— —— —

Quartz sand —— —

**)**

— • 1.04–1.44 • 0.422 • 1.08 • 1.49–1.61

— 1.60–1.61 —

— 0.086–0.263 1.3

0.4–0.6 1.1–1.7 2.6–2.75

0.33–0.95 0.13–1.66 1.06–2.7

• 1.24–1.63

—— —

—

— • 1.08–1.49

**Bulk density (g cm<sup>3</sup> )**

**Particle density (g cm<sup>3</sup> )**

—

1.26–1.55 2.28–2.67

1.63

#### **Figure 2.**

*Comparison all equations with different porosity, by using: a) Eq. 15; b) Eq. 16; c) Eq. 18; d) Eq. 19; e) Eq. 20; f) Eq. 21.*


**Table 2.**

*Material properties.*

where *θ* is the volumetric water content, *θ<sup>r</sup>* is the residual volumetric water content, and *θ<sup>s</sup>* is the saturated volumetric water content.

The normalization of relative permittivity gives privilege to define normalized volumetric water content, in this study it is assumed in three models:

$$\begin{aligned} \boldsymbol{\Theta} &= \boldsymbol{\varepsilon}\_{\text{N}} \, ^{0.5} \qquad & \qquad (\mathbf{a})\\ \boldsymbol{\Theta} &= \boldsymbol{\varepsilon}\_{\text{N}} \qquad & \qquad (\mathbf{b})\\ \boldsymbol{\Theta} &= \boldsymbol{\varepsilon}\_{\text{N}} \, ^{2} \qquad & \qquad (\mathbf{c}) \end{aligned} \tag{23}$$

**4. Experimental setup**

*DOI: http://dx.doi.org/10.5772/intechopen.89057*

**5. Result and discussion**

1.55 g cm<sup>3</sup>

**Figure 4.**

**81**

*compared with established models.*

as the top layer and 32nd layer as the bottom layer.

ECVT system consists of three parts: (i) sensors, (ii) data acquisition, and (iii) system reconstruction and visualization as shown in **Figure 3**. In this experiment we used 32 channel hexagonal sensors for the first experiment and 24 channel hexagonal sensors for the second experiment, as a soil column with diameter and height 11.5 and 27 cm, respectively. This column is divided into 32 layers with the 1st layer

*Dielectric Analysis Model for Measurement of Soil Moisture Water Content Using Electrical…*

In the first experiment, 3 l by volume of soil in the column was supplied. The soil material that was used in this study was sand collected from the Cisadane River in Tangerang, Indonesia. The soil contained 17% fine sand and 83% medium sand with

, respectively. In this experiment, the soil in the vessel was supplied with

porosity of soil 41.79%. The specific gravity and soil density were 2.663 and

water flow with a discharge of 7.2 ml/s until ponded condition, and the discharge was stopped when the pond of water level was at 2 cm above the surface soil. During ponded condition, the data capacitances were measured iteratively and sent to the computer. The data acquisition frequency was set to one frame per second. In the second experiment, 3 kg of silty sand was supplied into the column. After that, 1.4 liter of water was filled into the soil column using constant head method, in which height of the water is maintained constant by a distance of 4 cm from the surface of soil.

Normalized volumetric water content and relative permittivity relationship of several equations is shown in **Figure 4**. The first model of the proposed model seems to have a similar pattern with Topp et al. [20] and Malicki et al. [32] models. By adding some constants value to the first model, this model will be quite fit with

*Normalized volumetric water content and relative permittivity relationship of the three proposed models*

where *ε<sup>N</sup>* is normalized permittivity which can be calculated as

$$
\varepsilon\_N = \frac{\varepsilon - \varepsilon\_{air}}{\varepsilon\_{water} - \varepsilon\_{air}} \tag{24}
$$

where *ε*, *εair*, and εwater are actual relative permittivity measurement, relative permittivity of air, and relative permittivity of water, respectively (**Table 2**).

**Figure 3.** *Experimental setup [15].*

*Dielectric Analysis Model for Measurement of Soil Moisture Water Content Using Electrical… DOI: http://dx.doi.org/10.5772/intechopen.89057*

#### **4. Experimental setup**

where *θ* is the volumetric water content, *θ<sup>r</sup>* is the residual volumetric water

Dry soil 3–10 N, P, K, Ca, Mg, S, Cu, Zn, Fe, Mn, B, Cl, Na, H

The normalization of relative permittivity gives privilege to define normalized

Θ ¼ *ε<sup>N</sup>* ð Þ b

*<sup>ε</sup><sup>N</sup>* <sup>¼</sup> *<sup>ε</sup>* � *<sup>ε</sup>air εwater* � *εair*

where *ε*, *εair*, and εwater are actual relative permittivity measurement, relative permittivity of air, and relative permittivity of water, respectively (**Table 2**).

<sup>0</sup>*:*<sup>5</sup> ð Þ<sup>a</sup>

(23)

(24)

<sup>2</sup> ð Þ<sup>c</sup>

content, and *θ<sup>s</sup>* is the saturated volumetric water content.

*Modern Applications of Electrostatics and Dielectrics*

**Table 2.** *Material properties.*

**Figure 3.**

**80**

*Experimental setup [15].*

volumetric water content, in this study it is assumed in three models:

**Material Relative permittivity Elements** Air 1 N2, O2 Water 80 H2O

Θ ¼ *ε<sup>N</sup>*

Θ ¼ *ε*<sup>N</sup>

where *ε<sup>N</sup>* is normalized permittivity which can be calculated as

ECVT system consists of three parts: (i) sensors, (ii) data acquisition, and (iii) system reconstruction and visualization as shown in **Figure 3**. In this experiment we used 32 channel hexagonal sensors for the first experiment and 24 channel hexagonal sensors for the second experiment, as a soil column with diameter and height 11.5 and 27 cm, respectively. This column is divided into 32 layers with the 1st layer as the top layer and 32nd layer as the bottom layer.

In the first experiment, 3 l by volume of soil in the column was supplied. The soil material that was used in this study was sand collected from the Cisadane River in Tangerang, Indonesia. The soil contained 17% fine sand and 83% medium sand with porosity of soil 41.79%. The specific gravity and soil density were 2.663 and 1.55 g cm<sup>3</sup> , respectively. In this experiment, the soil in the vessel was supplied with water flow with a discharge of 7.2 ml/s until ponded condition, and the discharge was stopped when the pond of water level was at 2 cm above the surface soil. During ponded condition, the data capacitances were measured iteratively and sent to the computer. The data acquisition frequency was set to one frame per second.

In the second experiment, 3 kg of silty sand was supplied into the column. After that, 1.4 liter of water was filled into the soil column using constant head method, in which height of the water is maintained constant by a distance of 4 cm from the surface of soil.

### **5. Result and discussion**

Normalized volumetric water content and relative permittivity relationship of several equations is shown in **Figure 4**. The first model of the proposed model seems to have a similar pattern with Topp et al. [20] and Malicki et al. [32] models. By adding some constants value to the first model, this model will be quite fit with

**Figure 4.**

*Normalized volumetric water content and relative permittivity relationship of the three proposed models compared with established models.*

Modification to the first model (Eq. (23a)) to fit with previous models can be done by adding the constants by trial and error. The modification is shown below:

*Dielectric Analysis Model for Measurement of Soil Moisture Water Content Using Electrical…*

**Figure 5** shows the first model with modification (Eq. (25)) quite fit with the

Based on this result, Eq. (25) was used to analyze the normalized volumetric water content during water infiltration in soil. The results of water infiltration in soil can be seen in **Figure 6**. **Figure 6a** shows the images of normalized volumetric water content from red (i.e., dry condition, *ε<sup>N</sup>* = 0) to blue (i.e., saturated condition, *ε<sup>N</sup>* = 1) colors. The scale of the color means normalization value of relative

**Figure 6a** shows the image sequencing of water infiltration methods from 1, 50, 100, 250, 400, 600, 700, and 1000 s, respectively. In this figure, the position and movement of water per second can be seen clearly. **Figure 6b** shows normalized

*Result from experiment 2, soil water infiltration (a) and normalized volumetric water content of soil (blue line, 1st layer; red line, 9th layer; green line, 17th layer; cyan line, 25th layer; and yellow line, 32nd layer) (b).*

Topp and Malicki models.

**Figure 7.**

**83**

permittivity distribution in image.

*DOI: http://dx.doi.org/10.5772/intechopen.89057*

<sup>Θ</sup> <sup>¼</sup> <sup>0</sup>*:*<sup>9</sup> � ð Þ *<sup>ε</sup>N*‐0*:*<sup>015</sup> <sup>0</sup>*:*<sup>65</sup> (25)

*Normalized volumetric water content and relative permittivity relationship of the first model modification compared with established models.*

#### **Figure 6.**

*Result from experiment 1, soil water infiltration (a) and normalized volumetric water content of soil (blue line, 1st layer; red line, 9th layer; green line, 17th layer; cyan line, 25th layer; and yellow line, 32nd layer) (b).*

Topp and Malicki models. However, the second model looks similar to Roth et al. [21] model, which gives a linear relationship between normalized volumetric water content and relative permittivity. In contrast, the third model is out of fit with other models.

*Dielectric Analysis Model for Measurement of Soil Moisture Water Content Using Electrical… DOI: http://dx.doi.org/10.5772/intechopen.89057*

Modification to the first model (Eq. (23a)) to fit with previous models can be done by adding the constants by trial and error. The modification is shown below:

$$
\Theta = \mathbf{0}.\Theta \times \left(\varepsilon\_{N^-} \mathbf{0}.\mathbf{0}15\right)^{0.65} \tag{25}
$$

**Figure 5** shows the first model with modification (Eq. (25)) quite fit with the Topp and Malicki models.

Based on this result, Eq. (25) was used to analyze the normalized volumetric water content during water infiltration in soil. The results of water infiltration in soil can be seen in **Figure 6**. **Figure 6a** shows the images of normalized volumetric water content from red (i.e., dry condition, *ε<sup>N</sup>* = 0) to blue (i.e., saturated condition, *ε<sup>N</sup>* = 1) colors. The scale of the color means normalization value of relative permittivity distribution in image.

**Figure 6a** shows the image sequencing of water infiltration methods from 1, 50, 100, 250, 400, 600, 700, and 1000 s, respectively. In this figure, the position and movement of water per second can be seen clearly. **Figure 6b** shows normalized

#### **Figure 7.**

*Result from experiment 2, soil water infiltration (a) and normalized volumetric water content of soil (blue line, 1st layer; red line, 9th layer; green line, 17th layer; cyan line, 25th layer; and yellow line, 32nd layer) (b).*

Topp and Malicki models. However, the second model looks similar to Roth et al. [21] model, which gives a linear relationship between normalized volumetric water content and relative permittivity. In contrast, the third model is out of fit with other

*Result from experiment 1, soil water infiltration (a) and normalized volumetric water content of soil (blue line, 1st layer; red line, 9th layer; green line, 17th layer; cyan line, 25th layer; and yellow line, 32nd layer) (b).*

*Normalized volumetric water content and relative permittivity relationship of the first model modification*

models.

**82**

**Figure 6.**

**Figure 5.**

*compared with established models.*

*Modern Applications of Electrostatics and Dielectrics*

volumetric water content of soil for 32 layers during water infiltration in the soil. Blue, green, red, cyan, and yellow lines indicated 1, 9, 17, 25, and 32 soil layers, respectively. In the first layer, the normalized volumetric water content increases very fast and reaches a saturated condition at around 100–300 s. The degree of saturation for 1st layer starts to decrease after supplying of water was stopped at 300 s and relative constant was at around 0.3 from 700 to 1000 s. For 9th and 17th layers, the degree of saturation increases sequencing and reaches a saturated condition at a similar time (i.e., at around 450 s) and starts to decrease at 500 s, while for the 25th layer, the degree of saturation increases at around 300 s and reaches a stable condition at 700 s with a degree of saturation at around 0.85. Moreover at the 32nd layer, the degree of saturation increases at 630 s, and the degree of saturation reaches 0.2 at the end of the experiment (i.e., 1000 s). This experiment showed clearly the availability of air trapped at the bottom of soil in the vessel (see **Figure 6a**).

**Figure 7** also shows the image sequencing of water infiltration (**Figure 7a**) and normalized volumetric water content of each layer of the soil column (**Figure 7b**). In **Figure 7a**, it can be seen that the silty sand has a normalized relative permittivity value around 0.3 before water infiltrate to the soil column. It can be caused by two possibilities, (i) because of the moisture content stored in the soil and (ii) because the soil particles are very small and compact, so the porosity of soil is also small which causes no air cavity in the soil. From **Figure 7b** we can see how the water infiltrates into the each layer clearly. The first layer has increased drastically around the first 30 s. It is easily understood that the top layers will reach the maximum value of normalized volumetric water content first, because these layers get the first supply of water.

#### **6. Conclusion**

Some equations between the relative permittivity and volumetric water content have been described in this study. From this study, there is an equation that has demonstrated efficacy in conformity with the data of experimental results, that is, equation proposed by [32]. This equation uses the porosity factor as a parameter in the relationship between volumetric water content and relative permittivity.

Normalized volumetric water content of soil has been analyzed in this study using ECVT system. Normalized volumetric water content can be shown and analyzed layer per layer of soil column for every second. We found that the ECVT system has advantages in measuring soil water content which are nondestructive and noninvasive to the sample object, 3D image, and real-time monitoring for water infiltration.

**Author details**

**85**

Mukhlisin Muhammad<sup>1</sup>

\* and Saputra Almushfi<sup>2</sup>

*Dielectric Analysis Model for Measurement of Soil Moisture Water Content Using Electrical…*

*DOI: http://dx.doi.org/10.5772/intechopen.89057*

© 2019 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/ by/3.0), which permits unrestricted use, distribution, and reproduction in any medium,

1 Politeknik Negeri Semarang, Semarang, Indonesia

2 CTECH Labs Edwar Technology, Tangerang, Indonesia

\*Address all correspondence to: mmukhlis2@gmail.com

provided the original work is properly cited.

#### **Acknowledgements**

The authors wish to thank the Ministry of Research, Technology and Higher Education of Republic of Indonesia, for providing the research grand.

#### **Conflict of interest**

The authors declare that there is no conflict of interest.

*Dielectric Analysis Model for Measurement of Soil Moisture Water Content Using Electrical… DOI: http://dx.doi.org/10.5772/intechopen.89057*
