# Analysis of the Accuracy of an FDR Sensor in Soil Moisture Measurement under Laboratory and Field Conditions

Authors: Anzhen Qin; Dongfeng Ning; Zhandong Liu; Aiwang Duan
Journal: Journal of Sensors
 (2021)

---

## Abstract

Soil water content (SWC, % vol) is a key factor affecting plant growth and development. SWC measurement is vital to rational use of water resources for irrigation, and the accuracy of sensors in SWC measurement is of significant importance to smart data-driven irrigation. Here, a laboratory experiment and a field lysimetric experiment were conducted to evaluate the accuracy of Insentek sensors under various soil conditions (1.1 to 1.5 bulk densities and sand to clay soil textures) and irrigation levels (30, 45, and 60 mm), in 2018 and 2019. A microweighing lysimeter and oven-drying method were used as standard methods to compare the Insentek method. The root mean square error (RMSE, % vol) and relative prediction deviation (RPD) between the Insentek and microlysimetric SWC values were 0.89–1.04% vol and 5.6–6.8, respectively, under laboratory condition. The RPD value is larger than the threshold value of 4.0, indicating the accuracy of the Insentek sensors is reliable under laboratory condition. Except for 60 mm irrigation treatment, the RMSE between Insentek and the oven-drying method under field condition was 1.44–1.93% vol, and the RPD value was 1.56–1.93, lower than the threshold value of 4.0. The tiny gap between the Insentek sensor and soil may accelerate water infiltration along the probe 0-3 d after irrigation while increase air filling 5–7 d after irrigation, causing greater RMSE and lower RPD values. The dissatisfied performance in field condition may also be associated with the obvious drawbacks of oven-drying method, such as disturbance in soil sampling. When using oven-drying method to analyze the accuracy of the Insentek sensors in field condition, the concerns should be well addressed.


---

## Body

## 1. Introduction

Global demands for cereal grains are projected to increase by 50% by 2050. To produce more food, agriculture has become the largest water consumer accounting for more than 60% of water withdrawn from resources worldwide [1]. With the ever-increasing demands for food and water, modern agriculture should be developed toward the precision agriculture with efficient use of water and nutrient resources [2]. Nevertheless, the present efficiency of irrigation is low. In China, the mean utilization coefficient of irrigation water was 0.45 in 2019 [3]. Precision irrigation according to soil moisture dynamics is the key to achieving high water use efficiency (WUE, kg m–3). However, the main difficulty for precision irrigation lies in real-time and rapid monitoring of soil water content (SWC, % vol) [4]. Over the past decades, scientists and engineers have developed various sensors to monitor the dynamics of soil moisture, including capacitive sensors (e.g., time-domain reflectometers (TDR) [5] and frequency domain reflectometers (FDR) [6]), neutron probes [7], and tensiometers [8]. TDR is a less-disturbing methodology, but the price of TDR sensor is relatively high [9]. As for neutron probes, it is difficult for the sensors to separate the signal of soil moisture from roots, giving rise to higher SWC estimates in upper soil layer [7], whereas tensiometers do not perform well in dry soil environment [8].The Insentek sensor (Zhejiang Oriental Insentek Technology Ltd., Co, Hangzhou, China) is an FDR-type sensor that can automatically monitor SWC data. The sensor is based on the difference of relative dielectric permittivity (DP) in soils [4]. In general, DP of saturated soil at 20°C is 80, while the DP of dry soil is 2–5 [10–12]. Through detecting the DP of soil, SWC can be determined. Compared with TDR sensors, the Insentek sensor has the advantages of rapid speed of measurement, easy access to data, and low price in market [6]. The sensor is powered by a built-in lithium-polymer rechargeable battery charged by a solar panel (see Figure 1(a)). Also, a built-in GPS module with wireless communication function allows the traceability of collected data [13]. Users are able to check the real-time SWC data using their PC or cellphones. This provides convenience for users without frequent visits to experimental sites [14]. So far, the latest figure has shown that there have been more than 35,000 Insentek sensors already installed across China, forming a nationwide soil moisture monitoring network across the country. When the sensors’ year-round data were checked, they showed good continuity and stability [4]. However, in practice, we did find that there was an extremely tiny gap between the probe and soil, especially when the soil was dry. The gap allows the infiltration of water immediately after irrigation or heavy precipitation, which may overestimate SWC (see Figure 1(b)). However, when soils become drier, the gap will be filled with air, causing underestimates of SWC values since the DP of air is very low [15]. The problem of the gap exists in every tube-based SWC measurement methods, and it may inevitably cause errors in data.Figure 1
(a) Insentek sensor installed in maize plots at the Qitai experimental station in 2018. (b) Thermal image (24.8–34.8°C) of soil water infiltration area along an Insentek sensor in a lysimeter experiment 2 hr after a 10 mm irrigation at Xinxiang experiment station, in June 2018.
(a)(b)In addition, the measurement accuracy of FDR sensors can be affected by soil temperature [16, 17]. A temperature-dependent nonlinear model needs to be developed to compensate for the effects of temperature [18, 19]. Generally, the DP of soils decreases with increasing soil temperature, resulting in lower SWC estimates under higher temperature [20]. To minimize the negative effects, a temperature sensor and a moisture sensor were simultaneously installed in Insentek sensors. Engineers of Insentek sensors have performed thousands of experiments under various temperature and moisture levels in order to improve the accuracy of the sensors. Nevertheless, in published literature, the Insentek sensor studied in the present study has been characterized only once before [4]. In their published paper, the authors found that the Insentek sensor was an efficient tool in estimating maize evapotranspiration in the North China Plain with acceptable accuracy. However, the degree to which SWC values from Insentek sensors represent these of oven-drying method has not been well tested under both laboratory and field conditions.In this study, the accuracy of Insentek sensors was tested both in laboratory and field conditions. We filled the gap with diluted mud in the installation of the sensors trying to exclude its adverse influences on SWC measurement. To start the experiment, a microlysimeter system was designed to test the sensors in laboratory, while in the field, a standard oven-drying method was adopted to evaluate the data of Insentek sensors. The experiment may help us validate the accuracy and reliability of Insentek sensors and better understand the uncertainty and variability in actual measurement. We hypothesized that Insentek sensors were able to monitor real-time SWC accurately both in laboratory and field conditions. The objectives of this study are (i) to validate the accuracy of Insentek sensors using microweighing lysimeters in laboratory condition and oven-drying method in field condition and (ii) to find out what factors may contribute to the variability in data if any.


## 2. Materials and Methods

### 2.1. Site Description

The experiment was carried out in June to August 2018 and 2019, at the Xinxiang experimental station of the Institute of Farmland Irrigation, Chinese Academy of Agricultural Sciences (CAAS) (35° 19′ N, 113° 53′ E, 73.2 m a.s.l.). The place has a continent temperate monsoon climate. Long-term (2009–2019) annual mean temperature is 14.1°C, annual precipitation 582 mm, annual sunshine hour 2497 hr, and frost-free period 220 d. The soil is a silt loam soil [21]. Field holding capacity and permanent wilting point in 0–100 cm soil layers was 24% and 9% (by weight), average bulk density across the 0–100 cm soil depth was 1.42 g cm–3, and water table was detected more than 5 m below the soil surface. Soil available N, P, and K contents were 72, 17.8, and 100 mg kg–1, respectively. Soil organic carbon content was 6.54 g kg–1, and soil pH was around 7.5.


### 2.2. Experimental Design

The experiment consisted of a microweighing lysimeter experiment in the Key Laboratory of Crop Water Use and Regulation, CAAS, from June 1st to August 30th, 2018 (see Figure2(a)), and a nonweighing lysimeter experiment under an automatic rainproof shelter at the Xinxiang experimental station from June 1st to August 30th, 2019 (see Figure 2(b)). In the microlysimeter experiment, three microlysimetric systems were designed to evaluate the accuracy of Insentek sensors across various soil bulk densities from 1.1 to 1.5 g cm–3 and soil textures from sand to clay. Soils with different textures were sampled at various locations during April and May 2018 (see Table 1). When soils were collected, all stones and root debris were removed by riddling. The soils were then repacked carefully into a transparent acrylic (PMMA) tube. Different soil bulk densities were achieved by manually compacting. The Insentek sensor was installed vertically through the central circle of the soil column.Figure 2
Schematic diagram of (a) a microweighing lysimetric system and (b) a nonweighing lysimeter system.
(a)(b)Table 1
Soil particle size distribution for soils sampled at different locations, in April and May 2018.
LocationsLongitudeLatitudeAltitudeSand(>0.05 mm)Silt(0.002–0.05 mm)Clay(<0.002 mm)Soil textureHuixian county35° 30′113° 41′94 m95.2%4.8%0SandYuanyang county34° 59′114° 03′86 m53.4%41.2%5.4%Silt sandYanjin county35° 14′114° 15′69 m45.7%28.6%25.7%Sandy loamMuye district35° 20′113° 46′73 m41.1%34.8%24.1%Silt loamHongqi district35° 07′114° 10′71 m37.3%35.9%26.8%Silt loamChangyuan county35° 10′114° 41′62 m17.7%43.8%38.5%Loam clayQinyang county35° 04′112° 55′126 m8.8%35.1%56.1%ClaySoils collected at different locations in 2018 were classified according to their soil particle size distribution. Sand and silt sand were integrated as a sand soil, since its sand particle content was larger than 50% and clay content less than 10%. Similarly, sandy loam and silt loam were integrated as a loam soil, since its sand content was 25–50%, and silt content 25–40%; loam clay and clay were integrated as a clay soil, with its clay content larger than 35% and silt content larger than 30% [21]. The microweighing lysimetric system consisted of an electronic weighing balance, an Insentek sensor, and two soil columns. The precision of the balance was ±10–4 g. Since the maximum measuring radius of an Insentek sensor was 15 cm in soils, and the diameter of the sensor was 6.5 cm, the total diameter of the soil column was designed 36.5 cm and the height 20 cm, which was considered to fully cover the measuring range of the sensor. Each type of tested soil was packed in a transparent acrylic (PMMA) tube. The junction between the sensor and the PMMA tube was sealed with glue to prevent water drainage. During the measuring period, soils were first added distilled water to holding capacity and then were air dried to wilting point under natural environment. The measurements were performed in laboratory condition at a constant temperature of 26±0.5°C. After soils were mixed with water, soil columns were placed for a period of 12 hr to obtain a uniform moisture distribution. The weight changes of soils were automatically recorded by the microweighing lysimeters. In the meantime, SWC (% vol) was monitored by Insentek sensor. The battery of the Insentek sensor was fully charged in advance. The charged battery can power the sensor to work as long as two months without recharging.In large-scale lysimeter experiment, it consisted of nine lysimeters (3.3mlong×2.0mwide×2.3mdeep) packed uniformly with silt loam soil in the 0–100 cm soil layers and sandy loam soil in the 100–200 cm soil layers. The bottom 30 cm was filled with coarse sands and <3 cm gravels to permit drainage toward lysimeter outlet. Collecting buckets were placed at the bottom side of lysimeters to hold gravity drainage effluent. These lysimeters were arranged in two rows. Between the rows, a 2.5 m wide concrete platform was constructed level with the top of the lysimeters. Nine Insentek sensors were installed at the center of each lysimeter to compare SWC measured by the standard oven-drying method. Electrically operated rainproof shelters were installed for the lysimeters. The shelters were moved over all lysimeters before precipitation happened. Three irrigation levels were arranged in a complete random design, including 30 mm, 45 mm, and 60 mm per irrigation, respectively. Each plot was irrigated separately, and different amounts of water were applied using surface irrigation, and a precision flow meter (Shanghai Water Meter Manufacturing Ltd., Co, Shanghai, China) was installed to measure and control the irrigation amounts. Preirrigation SWC was measured both by Insentek sensor and oven-drying method. After irrigation, SWC was continuously monitored by Insentek sensors at 2 hr interval. Due to the labor cost of soil-core sampling, SWC was measured by an oven-drying method at 2 hr interval for the first 8 hr after irrigation, at 4 hr interval for the next 24 hr, and 6 hr interval for the remaining days. The measurement usually lasted 7 d and repeated for three times.


### 2.3. Soil Water Content and Infiltration Measurements

100 cm long Insentek sensors were used to monitor SWC (% vol) at 10 cm increment to a depth of 100 cm in field condition. 40 cm long Insentek sensors were adopted to monitor SWC at 20 cm and 40 cm depth in soil columns in laboratory condition. The sensors measured SWC every 10 min, and the average of two hours was stored in dataloggers and was transmitted to a local server via the 4G wireless network.Microlysimeters were designed to measure the weight changes of soil column every two hours, and the data were stored in dataloggers. We used the average of a whole day to calculate daily changes of soil column weight. SWC on theith day was derived from initial SWC and the difference in SWC (ΔSWC) using Equation (1):
(1)SWCi=SWCini−∑i=1ΔSWC,where SWCi is soil water content (% vol) on the ith day, SWCini is the initial SWC at the beginning of the experiment, and ΔSWC is the daily difference between SWC on the previous and the current measurement date.In this study,ΔSWC is derived from the changes of soil column weight (ΔWt, g) using Equation (2):
(2)ΔSWC=ΔWtH⋅π⋅R2−r2⋅ρ,where ΔSWC is the daily difference in SWC (% vol), ΔWt the daily change in soil column weight (g), H is the height of soil column (cm), R and r are the radius (cm) of the soil column and of the Insentek sensor, and ρ is the density of water (g cm–3).Soil gravimetric water content (%w/w) was measured at 10 cm increment to a depth of 100 cm using a soil auger (3 cm in diameter). Soil cores were randomly sampled about 20–30 cm away from Insentek sensor, then were oven-dried at 105°C for 8 hr after the fresh weight weighed. SWC was measured on a dry weight basis. Soil bulk density (g cm–3) at the same depth and interval was determined before and after the experiment using stainless steel ring method [22]. SWC (% vol) was calculated by multiplying soil gravimetric water content with soil bulk density (g cm–3) from the same layer.In this study, infiltration area was defined as a product of infiltration depth (cm) and duration of infiltration after irrigation (d), which indicated how deep and how long irrigation amount affected SWC. The unit of infiltration area is cm d. Infiltration area was derived from pixel area estimated using the software ImageJ (NIH, Maryland, USA). Since the maximum infiltration depth was 50 cm in the present study, the area was calculated within 0–50 cm soil depth using Equation (3):
(3)Sin=SmaxSpiSmax−pi,where Sin is infiltration area (cm d), Smax is the maximum infiltration area (350 cm d), which is a product of maximum infiltration depth (50 cm) and maximum duration after irrigation (7 d), Spi is pixel area, and Smax–pi is the total pixel area.


### 2.4. Goodness-of-Fit Test

The evaluation factors include slope and coefficient of determination (R2) for linear regression between the microlysimetric and Insentek data, and between the oven-drying and Insentek data. The slope of a line characterizes the direction of a line. The ideal slope is 1.0, indicating SWC data by the Insentek sensor are equal to those of microlysimeter or oven-drying method. The R2 describes the proportion of variance in SWC data using microlysimeter or oven-drying method explained by the Insentek data. In addition, the root mean square error (RMSE) and relative prediction deviation (RPD) were adopted to evaluate the accuracy of Insentek sensors. The RMSE was used to investigate the differences between the microlysimeter and Insentek data, and the between oven-drying and Insentek data. The RMSE was calculated using Equation (4):
(4)RMSE=∑i=1nxinsen−xj2/n,where RMSE is root mean square error, xInsen and xj are corresponding SWC values estimated by the Insentek and microlysimetric or oven-drying method, respectively, and n is the number of values evaluated. The smaller the RMSE values are, the more accurate the Insentek sensor is.RPD is used to indicate the reliability of Insentek data. RPD was calculated using Equation (5):
(5)RPD=STDEVxj/RMSE,where RPD is relative prediction deviation and STDEVxj is standard deviation of microlysimetric or oven-drying SWC values. RPD≥4.0 indicates Insentek data are reliable; 1.4<RPD<4.0 means the data are feasible but needs to be improved; RPD≤1.4 indicates the data are unreliable [4].


### 2.5. Statistical Analysis

Data were analyzed using an analysis of variance with the Statistical Analysis Software (version 19.0, SPSS Inc., Chicago, IL, USA). Significance was declared at the probability level of 0.05, unless otherwise stated. Relationships among SWC by the Insentek, microlysimeter, and oven-drying methods were analyzed by means of the Levenberg–Marquardt Algorithm. Figures were plotted using Original Pro 9.1 (Origin Lab Corporation, Northampton, MA, USA).


## 2.1. Site Description

The experiment was carried out in June to August 2018 and 2019, at the Xinxiang experimental station of the Institute of Farmland Irrigation, Chinese Academy of Agricultural Sciences (CAAS) (35° 19′ N, 113° 53′ E, 73.2 m a.s.l.). The place has a continent temperate monsoon climate. Long-term (2009–2019) annual mean temperature is 14.1°C, annual precipitation 582 mm, annual sunshine hour 2497 hr, and frost-free period 220 d. The soil is a silt loam soil [21]. Field holding capacity and permanent wilting point in 0–100 cm soil layers was 24% and 9% (by weight), average bulk density across the 0–100 cm soil depth was 1.42 g cm–3, and water table was detected more than 5 m below the soil surface. Soil available N, P, and K contents were 72, 17.8, and 100 mg kg–1, respectively. Soil organic carbon content was 6.54 g kg–1, and soil pH was around 7.5.


## 2.2. Experimental Design

The experiment consisted of a microweighing lysimeter experiment in the Key Laboratory of Crop Water Use and Regulation, CAAS, from June 1st to August 30th, 2018 (see Figure2(a)), and a nonweighing lysimeter experiment under an automatic rainproof shelter at the Xinxiang experimental station from June 1st to August 30th, 2019 (see Figure 2(b)). In the microlysimeter experiment, three microlysimetric systems were designed to evaluate the accuracy of Insentek sensors across various soil bulk densities from 1.1 to 1.5 g cm–3 and soil textures from sand to clay. Soils with different textures were sampled at various locations during April and May 2018 (see Table 1). When soils were collected, all stones and root debris were removed by riddling. The soils were then repacked carefully into a transparent acrylic (PMMA) tube. Different soil bulk densities were achieved by manually compacting. The Insentek sensor was installed vertically through the central circle of the soil column.Figure 2
Schematic diagram of (a) a microweighing lysimetric system and (b) a nonweighing lysimeter system.
(a)(b)Table 1
Soil particle size distribution for soils sampled at different locations, in April and May 2018.
LocationsLongitudeLatitudeAltitudeSand(>0.05 mm)Silt(0.002–0.05 mm)Clay(<0.002 mm)Soil textureHuixian county35° 30′113° 41′94 m95.2%4.8%0SandYuanyang county34° 59′114° 03′86 m53.4%41.2%5.4%Silt sandYanjin county35° 14′114° 15′69 m45.7%28.6%25.7%Sandy loamMuye district35° 20′113° 46′73 m41.1%34.8%24.1%Silt loamHongqi district35° 07′114° 10′71 m37.3%35.9%26.8%Silt loamChangyuan county35° 10′114° 41′62 m17.7%43.8%38.5%Loam clayQinyang county35° 04′112° 55′126 m8.8%35.1%56.1%ClaySoils collected at different locations in 2018 were classified according to their soil particle size distribution. Sand and silt sand were integrated as a sand soil, since its sand particle content was larger than 50% and clay content less than 10%. Similarly, sandy loam and silt loam were integrated as a loam soil, since its sand content was 25–50%, and silt content 25–40%; loam clay and clay were integrated as a clay soil, with its clay content larger than 35% and silt content larger than 30% [21]. The microweighing lysimetric system consisted of an electronic weighing balance, an Insentek sensor, and two soil columns. The precision of the balance was ±10–4 g. Since the maximum measuring radius of an Insentek sensor was 15 cm in soils, and the diameter of the sensor was 6.5 cm, the total diameter of the soil column was designed 36.5 cm and the height 20 cm, which was considered to fully cover the measuring range of the sensor. Each type of tested soil was packed in a transparent acrylic (PMMA) tube. The junction between the sensor and the PMMA tube was sealed with glue to prevent water drainage. During the measuring period, soils were first added distilled water to holding capacity and then were air dried to wilting point under natural environment. The measurements were performed in laboratory condition at a constant temperature of 26±0.5°C. After soils were mixed with water, soil columns were placed for a period of 12 hr to obtain a uniform moisture distribution. The weight changes of soils were automatically recorded by the microweighing lysimeters. In the meantime, SWC (% vol) was monitored by Insentek sensor. The battery of the Insentek sensor was fully charged in advance. The charged battery can power the sensor to work as long as two months without recharging.In large-scale lysimeter experiment, it consisted of nine lysimeters (3.3mlong×2.0mwide×2.3mdeep) packed uniformly with silt loam soil in the 0–100 cm soil layers and sandy loam soil in the 100–200 cm soil layers. The bottom 30 cm was filled with coarse sands and <3 cm gravels to permit drainage toward lysimeter outlet. Collecting buckets were placed at the bottom side of lysimeters to hold gravity drainage effluent. These lysimeters were arranged in two rows. Between the rows, a 2.5 m wide concrete platform was constructed level with the top of the lysimeters. Nine Insentek sensors were installed at the center of each lysimeter to compare SWC measured by the standard oven-drying method. Electrically operated rainproof shelters were installed for the lysimeters. The shelters were moved over all lysimeters before precipitation happened. Three irrigation levels were arranged in a complete random design, including 30 mm, 45 mm, and 60 mm per irrigation, respectively. Each plot was irrigated separately, and different amounts of water were applied using surface irrigation, and a precision flow meter (Shanghai Water Meter Manufacturing Ltd., Co, Shanghai, China) was installed to measure and control the irrigation amounts. Preirrigation SWC was measured both by Insentek sensor and oven-drying method. After irrigation, SWC was continuously monitored by Insentek sensors at 2 hr interval. Due to the labor cost of soil-core sampling, SWC was measured by an oven-drying method at 2 hr interval for the first 8 hr after irrigation, at 4 hr interval for the next 24 hr, and 6 hr interval for the remaining days. The measurement usually lasted 7 d and repeated for three times.


## 2.3. Soil Water Content and Infiltration Measurements

100 cm long Insentek sensors were used to monitor SWC (% vol) at 10 cm increment to a depth of 100 cm in field condition. 40 cm long Insentek sensors were adopted to monitor SWC at 20 cm and 40 cm depth in soil columns in laboratory condition. The sensors measured SWC every 10 min, and the average of two hours was stored in dataloggers and was transmitted to a local server via the 4G wireless network.Microlysimeters were designed to measure the weight changes of soil column every two hours, and the data were stored in dataloggers. We used the average of a whole day to calculate daily changes of soil column weight. SWC on theith day was derived from initial SWC and the difference in SWC (ΔSWC) using Equation (1):
(1)SWCi=SWCini−∑i=1ΔSWC,where SWCi is soil water content (% vol) on the ith day, SWCini is the initial SWC at the beginning of the experiment, and ΔSWC is the daily difference between SWC on the previous and the current measurement date.In this study,ΔSWC is derived from the changes of soil column weight (ΔWt, g) using Equation (2):
(2)ΔSWC=ΔWtH⋅π⋅R2−r2⋅ρ,where ΔSWC is the daily difference in SWC (% vol), ΔWt the daily change in soil column weight (g), H is the height of soil column (cm), R and r are the radius (cm) of the soil column and of the Insentek sensor, and ρ is the density of water (g cm–3).Soil gravimetric water content (%w/w) was measured at 10 cm increment to a depth of 100 cm using a soil auger (3 cm in diameter). Soil cores were randomly sampled about 20–30 cm away from Insentek sensor, then were oven-dried at 105°C for 8 hr after the fresh weight weighed. SWC was measured on a dry weight basis. Soil bulk density (g cm–3) at the same depth and interval was determined before and after the experiment using stainless steel ring method [22]. SWC (% vol) was calculated by multiplying soil gravimetric water content with soil bulk density (g cm–3) from the same layer.In this study, infiltration area was defined as a product of infiltration depth (cm) and duration of infiltration after irrigation (d), which indicated how deep and how long irrigation amount affected SWC. The unit of infiltration area is cm d. Infiltration area was derived from pixel area estimated using the software ImageJ (NIH, Maryland, USA). Since the maximum infiltration depth was 50 cm in the present study, the area was calculated within 0–50 cm soil depth using Equation (3):
(3)Sin=SmaxSpiSmax−pi,where Sin is infiltration area (cm d), Smax is the maximum infiltration area (350 cm d), which is a product of maximum infiltration depth (50 cm) and maximum duration after irrigation (7 d), Spi is pixel area, and Smax–pi is the total pixel area.


## 2.4. Goodness-of-Fit Test

The evaluation factors include slope and coefficient of determination (R2) for linear regression between the microlysimetric and Insentek data, and between the oven-drying and Insentek data. The slope of a line characterizes the direction of a line. The ideal slope is 1.0, indicating SWC data by the Insentek sensor are equal to those of microlysimeter or oven-drying method. The R2 describes the proportion of variance in SWC data using microlysimeter or oven-drying method explained by the Insentek data. In addition, the root mean square error (RMSE) and relative prediction deviation (RPD) were adopted to evaluate the accuracy of Insentek sensors. The RMSE was used to investigate the differences between the microlysimeter and Insentek data, and the between oven-drying and Insentek data. The RMSE was calculated using Equation (4):
(4)RMSE=∑i=1nxinsen−xj2/n,where RMSE is root mean square error, xInsen and xj are corresponding SWC values estimated by the Insentek and microlysimetric or oven-drying method, respectively, and n is the number of values evaluated. The smaller the RMSE values are, the more accurate the Insentek sensor is.RPD is used to indicate the reliability of Insentek data. RPD was calculated using Equation (5):
(5)RPD=STDEVxj/RMSE,where RPD is relative prediction deviation and STDEVxj is standard deviation of microlysimetric or oven-drying SWC values. RPD≥4.0 indicates Insentek data are reliable; 1.4<RPD<4.0 means the data are feasible but needs to be improved; RPD≤1.4 indicates the data are unreliable [4].


## 2.5. Statistical Analysis

Data were analyzed using an analysis of variance with the Statistical Analysis Software (version 19.0, SPSS Inc., Chicago, IL, USA). Significance was declared at the probability level of 0.05, unless otherwise stated. Relationships among SWC by the Insentek, microlysimeter, and oven-drying methods were analyzed by means of the Levenberg–Marquardt Algorithm. Figures were plotted using Original Pro 9.1 (Origin Lab Corporation, Northampton, MA, USA).


## 3. Results and Discussion

### 3.1. Comparison of Soil Water Content between Insentek and Microlysimetric Methods

SWC measured by Insentek sensor had a significant positive linear correlation with the data calculated by microlysimetric system (see Figure3). Coefficient of determination (R2) of the linear functions for sand, loam, and clay soils was between 0.995 and 0.997, and that for different soil bulk densities (1.1–1.5 g cm–3) was between 0.996 and 1.000. The correlation between microlysimetric and Insentek data was extremely significant (p<0.001). The soil physical property test showed that the wilting point of sandy and loam soil was significantly lower than that of clay, while the field holding capacity of loam and clay soil was higher than that of sandy soil (see Table 2). Consequently, the available SWC of loam soil is significantly greater than that of sandy and clay soils. According to the goodness-of-fit test, the RMSE between Insentek and microlysimetric SWC values is between 0.89 and 1.04% vol, and the RPD is between 5.6 and 6.8, larger than the threshold value (4.0), which indicates that the accuracy of Insentek sensor is high in laboratory condition. The best range of bulk density for the Insentek sensor in SWC measurement is 1.1–1.4 g cm–3. Under high bulk density of 1.5 g cm–3, the RMSE remains 1.48% vol, and RPD is 3.7, less than the threshold value.Figure 3
Linear correlation between soil water content (% vol) measured using the Insentek and microweighing lysimetric methods across (a) different soil textures (sand, loam, and clay) and (b) soil bulk densities (1.1 to 1.5 g cm–3).
(a)(b)Table 2
Goodness-of-fit test between soil water content (SWC, % vol) calculated based on Insentek sensor and oven-drying method across different soil textures and soil bulk densities.
TreatmentsMicroweighing lysimetric values (% vol)Goodness-of-fitR2pWPFCAverage SWCRMSERPDSoil textureSand3.46 c28.6 b18.09 c1.0446.780.9970.0001Loam7.26 b35.1 a22.71 b0.9277.990.9950.0001Clay17.7 a33.7 a25.35 a0.8975.610.9970.0001Soil bulk density (g cm–3)1.17.47 b25.1 c15.6 c0.9215.930.9970.00011.217.7 a35.1 a26.4 a1.1644.680.9980.00011.37.14 b33.7 a23.8 b0.8637.880.9990.00011.43.46 c28.8 b17.9 c0.32022.81.0000.00011.59.22 b28.2 b21.1 b1.4823.710.9960.0001Note: WP is wilt point; FC is field holding capacity;R2 is coefficient of determination; RMSE is root mean square error; RPD is relative prediction deviation, and values measured by Insentek sensor are reliable with RPD≥4.0; and p is probability. Different letters in each column indicate significant differences at p<0.05.So far, the standard oven-drying method is considered the most accurate method of SWC measurement [22]. Similarly, the weighing lysimetric method is regarded the most reliable approach to estimate crop evapotranspiration [6]. However, the oven-drying method is destructive, laborious, and time-consuming. That promotes the development of nondestructive, low-cost technique for SWC measurement. The Insentek sensor is an inexpensive and easy-to-use FDR probe, compared with TDR sensors. Before being applied to precision irrigation, Insentek sensors should be proven accurate and reliable in SWC measurement under a variety of field conditions. In laboratory condition, Insentek sensors have better performance in clay and loam soils compared to sandy soils, which is consistent with the results of Adla et al. [23]. The reason may be related to more reasonable range of available SWC for loam soils than sandy soils.


### 3.2. Comparison of Soil Water Content and Infiltration between Insentek and Oven-Drying Methods

SWC prior to irrigation was similar between Insentek and oven-drying methods (see Figure4). However, postirrigation SWC measured by Insentek sensors was larger than that of oven-drying method. Specifically, SWC in 0–10 cm soil layer decreased noticeably with time under 60 mm treatment. For 30 mm treatment, SWC in 10–20 cm layer remained stable within 1.0 d after irrigation and then gradually declined. As for SWC in 20–30 cm soil layer, it increased gradually and reached the maximum SWC values 2 d after irrigation under 30 mm treatment, while for SWC in 0–30 cm soil layer, it reached the maximum values within 1 d after irrigation under 45 and 60 mm treatments and remained stable for the remaining 6 d. Under 45 mm and 60 cm treatments, the mean SWC in 30–50 cm soil layers increased gradually with time and eventually reached the maximum values 3.0 d after irrigation. Averaged across the three irrigation treatments, postirrigation SWC of 10–20 cm layer measured by Insentek sensors is 7.0% higher than that of oven-drying method (33.3% vs. 31.1%). Similarly, SWC of 20–30 cm soil layer measured by Insentek sensors is 5.6% higher, and that of 30–50 cm layer is 4.8% higher (p<0.05), respectively, than that of the oven-drying method. Irrigation treatments showed nonsignificant effect (p<0.05) on SWC below 50 cm soil layer. Besides, SWC below 50 cm layer was similar between the Insentek and oven-drying methods.Figure 4
Dynamics of soil volumetric water content (% vol) for each 10 cm soil layer to a depth of 0–100 cm measured using the Insentek and oven-drying methods before and after irrigation events applied at the amounts of 30 mm (a, d), 45 mm (b, e), and 60 mm (c, f), respectively.Figure5 shows the differences in infiltration area between the Insentek and oven-drying methods. The maximum infiltration depth was detected 50 cm under the three irrigation treatments. Although the infiltration depth measured by the two methods was generally the same, the duration of infiltration (d) after irrigation measured by the oven-drying method was much shorter than that of the Insentek method, especially for the 30 mm treatment. Under the 30 mm treatment, the duration of infiltration for SWC ranging from 32 to 36% vol was 4.3 d by the Insentek method, 3.9 d longer than that of the oven-drying method (see Table 3).Figure 5
Contour plot of soil volumetric water content (% vol) in 0–100 cm soil layer measured using the Insentek and oven-drying methods before and after different irrigation events applied at the amounts of 30 mm (a, d), 45 mm (b, e), and 60 mm (c, f), respectively.Table 3
Infiltration area (cm d) within the 0–50 cm soil depth converted using pixel area.
TreatmentsIrrigation quota (mm)Range of SWC (% vol)Pixel area within 0–50 cm depthPercent of total pixel area (%)Infiltration depth (cm)Duration after irrigation (d)Infiltration area within 0–50 cm depth (cm d)Maximum value0–3610100507.0350Insentek method3032–360.686.8184.32428–321.6316.3215.85724–284.6146.1437.01614532–360.656.5311.72328–323.3133.1385.811624–283.6936.9487.01296032–360.717.1361.82528–324.4744.7427.015624–282.0620.6507.072Oven-drying method3032–360.0191.9120.46.728–321.0410.4192.73624–284.2842.8397.01504532–360.414.1251.11428–322.0520.5333.37224–284.1741.7456.91466032–360.232.3260.6828–323.9839.8415.613924–283.8138.1507.0133Note: maximum infiltration depth is 50 cm under the three irrigation treatments. Pixel area is calculated using the software ImageJ (NIH, Maryland, USA). Infiltration area is the product of infiltration depth (cm) and duration after irrigation (d).In this study, the difference in infiltration area between the Insentek and oven-drying methods was larger in the upper 0–10 cm depth, but it became smaller with soil depth. The result was consistent with the findings of Domínguez–Niño et al. [24]. In addition, other factors, such as soil temperature, were considered attributable to influence sensor accuracy [25]. Diurnal variations in soil temperature are heterogeneous across different soil depths [26] and are generally larger in the upper soil layer, resulting in disturbance in sensor measurement. To solve the problem, a temperature sensor and a moisture sensor have been installed in the Insentek sensors, and then, thousands of experiments were performed to develop proper temperature compensation models to improve the accuracy of the sensors [27]. In laboratory with a constant temperature, our results proved that soil texture and bulk density had little impact on the accuracy of Insentek sensors, except sandy soil with high bulk density. The result was consistent with Domsch et al. [28]. In general, sandy soil texture and high bulk density are associated with high thermal conductivity of soil [29], which negatively affects the soil dielectric constant [30], and thus SWC data.


### 3.3. Goodness-of-Fit Test between the Insentek and Oven-Drying Methods

Slopes for equations between Insentek and oven-drying methods were <1.0, indicating SWC measured by the oven-drying method was constantly lower, which was probably due to the destruction of soils causing increases in soil evaporation in upper soils and thus leading to lower SWC (see Table4). In soil sampling, great efforts were made to reduce destructive effects on soils, leaving soils within the measuring radius (15 cm) of the Insentek sensor unchanged. During the experiment, sampling locations were kept a distance of 20–30 cm away from the Insentek sensor. After soil sampling, holes were filled with original soils. However, some artificial factors such as compactness of soil and destruction of macro soil pores inevitably contributed to the variability of experimental results. We suggest that disturbance to the soils should be minimized when choosing the oven-drying method to calibrate SWC sensors. The overestimates of SWC and infiltration area may also be associated with the tiny gap between the tube and soil. It has been learned that the outer layer of the Insentek sensor was made of chlorinated polyvinyl chloride with high chlorinated polyethylene (HCPE) coating, which helped improve the corrosion resistance to soil chemicals. To eliminate the impact of the gap on SWC measurement, it is suggested to adopt proper chemical material that has the same thermal expansion coefficient as soils.Table 4
Goodness-of-fit test for soil water content (% vol) between the Insentek and oven-drying methods.
TreatmentsLinear regression equationAverage SWC (% vol)Goodness-of-fitR2pInsentekOven-dryingRMSERPDIrrigation amount (mm)30y=0.9623x24.723.81.4361.560.7780.0145y=0.9671x26.225.31.9301.930.7840.0160y=0.9819x25.224.80.9774.040.9120.001Note:y is soil water content (SWC, % vol) of Insentek method; x is SWC of oven-drying method; RMSE is root mean square error (% vol); RPD is relative prediction deviation; R2 is coefficient of determination; and p is probability.


## 3.1. Comparison of Soil Water Content between Insentek and Microlysimetric Methods

SWC measured by Insentek sensor had a significant positive linear correlation with the data calculated by microlysimetric system (see Figure3). Coefficient of determination (R2) of the linear functions for sand, loam, and clay soils was between 0.995 and 0.997, and that for different soil bulk densities (1.1–1.5 g cm–3) was between 0.996 and 1.000. The correlation between microlysimetric and Insentek data was extremely significant (p<0.001). The soil physical property test showed that the wilting point of sandy and loam soil was significantly lower than that of clay, while the field holding capacity of loam and clay soil was higher than that of sandy soil (see Table 2). Consequently, the available SWC of loam soil is significantly greater than that of sandy and clay soils. According to the goodness-of-fit test, the RMSE between Insentek and microlysimetric SWC values is between 0.89 and 1.04% vol, and the RPD is between 5.6 and 6.8, larger than the threshold value (4.0), which indicates that the accuracy of Insentek sensor is high in laboratory condition. The best range of bulk density for the Insentek sensor in SWC measurement is 1.1–1.4 g cm–3. Under high bulk density of 1.5 g cm–3, the RMSE remains 1.48% vol, and RPD is 3.7, less than the threshold value.Figure 3
Linear correlation between soil water content (% vol) measured using the Insentek and microweighing lysimetric methods across (a) different soil textures (sand, loam, and clay) and (b) soil bulk densities (1.1 to 1.5 g cm–3).
(a)(b)Table 2
Goodness-of-fit test between soil water content (SWC, % vol) calculated based on Insentek sensor and oven-drying method across different soil textures and soil bulk densities.
TreatmentsMicroweighing lysimetric values (% vol)Goodness-of-fitR2pWPFCAverage SWCRMSERPDSoil textureSand3.46 c28.6 b18.09 c1.0446.780.9970.0001Loam7.26 b35.1 a22.71 b0.9277.990.9950.0001Clay17.7 a33.7 a25.35 a0.8975.610.9970.0001Soil bulk density (g cm–3)1.17.47 b25.1 c15.6 c0.9215.930.9970.00011.217.7 a35.1 a26.4 a1.1644.680.9980.00011.37.14 b33.7 a23.8 b0.8637.880.9990.00011.43.46 c28.8 b17.9 c0.32022.81.0000.00011.59.22 b28.2 b21.1 b1.4823.710.9960.0001Note: WP is wilt point; FC is field holding capacity;R2 is coefficient of determination; RMSE is root mean square error; RPD is relative prediction deviation, and values measured by Insentek sensor are reliable with RPD≥4.0; and p is probability. Different letters in each column indicate significant differences at p<0.05.So far, the standard oven-drying method is considered the most accurate method of SWC measurement [22]. Similarly, the weighing lysimetric method is regarded the most reliable approach to estimate crop evapotranspiration [6]. However, the oven-drying method is destructive, laborious, and time-consuming. That promotes the development of nondestructive, low-cost technique for SWC measurement. The Insentek sensor is an inexpensive and easy-to-use FDR probe, compared with TDR sensors. Before being applied to precision irrigation, Insentek sensors should be proven accurate and reliable in SWC measurement under a variety of field conditions. In laboratory condition, Insentek sensors have better performance in clay and loam soils compared to sandy soils, which is consistent with the results of Adla et al. [23]. The reason may be related to more reasonable range of available SWC for loam soils than sandy soils.


## 3.2. Comparison of Soil Water Content and Infiltration between Insentek and Oven-Drying Methods

SWC prior to irrigation was similar between Insentek and oven-drying methods (see Figure4). However, postirrigation SWC measured by Insentek sensors was larger than that of oven-drying method. Specifically, SWC in 0–10 cm soil layer decreased noticeably with time under 60 mm treatment. For 30 mm treatment, SWC in 10–20 cm layer remained stable within 1.0 d after irrigation and then gradually declined. As for SWC in 20–30 cm soil layer, it increased gradually and reached the maximum SWC values 2 d after irrigation under 30 mm treatment, while for SWC in 0–30 cm soil layer, it reached the maximum values within 1 d after irrigation under 45 and 60 mm treatments and remained stable for the remaining 6 d. Under 45 mm and 60 cm treatments, the mean SWC in 30–50 cm soil layers increased gradually with time and eventually reached the maximum values 3.0 d after irrigation. Averaged across the three irrigation treatments, postirrigation SWC of 10–20 cm layer measured by Insentek sensors is 7.0% higher than that of oven-drying method (33.3% vs. 31.1%). Similarly, SWC of 20–30 cm soil layer measured by Insentek sensors is 5.6% higher, and that of 30–50 cm layer is 4.8% higher (p<0.05), respectively, than that of the oven-drying method. Irrigation treatments showed nonsignificant effect (p<0.05) on SWC below 50 cm soil layer. Besides, SWC below 50 cm layer was similar between the Insentek and oven-drying methods.Figure 4
Dynamics of soil volumetric water content (% vol) for each 10 cm soil layer to a depth of 0–100 cm measured using the Insentek and oven-drying methods before and after irrigation events applied at the amounts of 30 mm (a, d), 45 mm (b, e), and 60 mm (c, f), respectively.Figure5 shows the differences in infiltration area between the Insentek and oven-drying methods. The maximum infiltration depth was detected 50 cm under the three irrigation treatments. Although the infiltration depth measured by the two methods was generally the same, the duration of infiltration (d) after irrigation measured by the oven-drying method was much shorter than that of the Insentek method, especially for the 30 mm treatment. Under the 30 mm treatment, the duration of infiltration for SWC ranging from 32 to 36% vol was 4.3 d by the Insentek method, 3.9 d longer than that of the oven-drying method (see Table 3).Figure 5
Contour plot of soil volumetric water content (% vol) in 0–100 cm soil layer measured using the Insentek and oven-drying methods before and after different irrigation events applied at the amounts of 30 mm (a, d), 45 mm (b, e), and 60 mm (c, f), respectively.Table 3
Infiltration area (cm d) within the 0–50 cm soil depth converted using pixel area.
TreatmentsIrrigation quota (mm)Range of SWC (% vol)Pixel area within 0–50 cm depthPercent of total pixel area (%)Infiltration depth (cm)Duration after irrigation (d)Infiltration area within 0–50 cm depth (cm d)Maximum value0–3610100507.0350Insentek method3032–360.686.8184.32428–321.6316.3215.85724–284.6146.1437.01614532–360.656.5311.72328–323.3133.1385.811624–283.6936.9487.01296032–360.717.1361.82528–324.4744.7427.015624–282.0620.6507.072Oven-drying method3032–360.0191.9120.46.728–321.0410.4192.73624–284.2842.8397.01504532–360.414.1251.11428–322.0520.5333.37224–284.1741.7456.91466032–360.232.3260.6828–323.9839.8415.613924–283.8138.1507.0133Note: maximum infiltration depth is 50 cm under the three irrigation treatments. Pixel area is calculated using the software ImageJ (NIH, Maryland, USA). Infiltration area is the product of infiltration depth (cm) and duration after irrigation (d).In this study, the difference in infiltration area between the Insentek and oven-drying methods was larger in the upper 0–10 cm depth, but it became smaller with soil depth. The result was consistent with the findings of Domínguez–Niño et al. [24]. In addition, other factors, such as soil temperature, were considered attributable to influence sensor accuracy [25]. Diurnal variations in soil temperature are heterogeneous across different soil depths [26] and are generally larger in the upper soil layer, resulting in disturbance in sensor measurement. To solve the problem, a temperature sensor and a moisture sensor have been installed in the Insentek sensors, and then, thousands of experiments were performed to develop proper temperature compensation models to improve the accuracy of the sensors [27]. In laboratory with a constant temperature, our results proved that soil texture and bulk density had little impact on the accuracy of Insentek sensors, except sandy soil with high bulk density. The result was consistent with Domsch et al. [28]. In general, sandy soil texture and high bulk density are associated with high thermal conductivity of soil [29], which negatively affects the soil dielectric constant [30], and thus SWC data.


## 3.3. Goodness-of-Fit Test between the Insentek and Oven-Drying Methods

Slopes for equations between Insentek and oven-drying methods were <1.0, indicating SWC measured by the oven-drying method was constantly lower, which was probably due to the destruction of soils causing increases in soil evaporation in upper soils and thus leading to lower SWC (see Table4). In soil sampling, great efforts were made to reduce destructive effects on soils, leaving soils within the measuring radius (15 cm) of the Insentek sensor unchanged. During the experiment, sampling locations were kept a distance of 20–30 cm away from the Insentek sensor. After soil sampling, holes were filled with original soils. However, some artificial factors such as compactness of soil and destruction of macro soil pores inevitably contributed to the variability of experimental results. We suggest that disturbance to the soils should be minimized when choosing the oven-drying method to calibrate SWC sensors. The overestimates of SWC and infiltration area may also be associated with the tiny gap between the tube and soil. It has been learned that the outer layer of the Insentek sensor was made of chlorinated polyvinyl chloride with high chlorinated polyethylene (HCPE) coating, which helped improve the corrosion resistance to soil chemicals. To eliminate the impact of the gap on SWC measurement, it is suggested to adopt proper chemical material that has the same thermal expansion coefficient as soils.Table 4
Goodness-of-fit test for soil water content (% vol) between the Insentek and oven-drying methods.
TreatmentsLinear regression equationAverage SWC (% vol)Goodness-of-fitR2pInsentekOven-dryingRMSERPDIrrigation amount (mm)30y=0.9623x24.723.81.4361.560.7780.0145y=0.9671x26.225.31.9301.930.7840.0160y=0.9819x25.224.80.9774.040.9120.001Note:y is soil water content (SWC, % vol) of Insentek method; x is SWC of oven-drying method; RMSE is root mean square error (% vol); RPD is relative prediction deviation; R2 is coefficient of determination; and p is probability.


## 4. Conclusions

We performed a comparative study using microweighing lysimeter and oven-drying methods to assess the accuracy of Insentek sensors in SWC measurement under laboratory and field conditions. The results showed that, in laboratory condition, the SWC of Insentek sensors had a significant (p<0.001) positive correlation to that of microweighing lysimeter across various soil textures and bulk densities, except sandy soil with high bulk density. In field condition, the Insentek method monitored the same infiltration depth as the oven-drying method did, but the duration of infiltration after irrigation was longer by the Insentek method. The accuracy of the Insentek sensor tested by the microlysimetric method in laboratory condition was higher than that of oven-drying method in field condition. The discrepancy in the SWC data between the Insentek and oven-drying methods was likely attributable to the tiny gap between the tube and soil. Also, the disturbance to soils in the oven-drying method may also contribute to variability in SWC data. Compared with the microlysimetric method, the oven-drying method with continuous disturbance in soil sampling may have drawbacks in assessing the accuracy of Insentek sensor in field conditions. We consider the microweighing lysimetric method an appropriate approach to assess the accuracy and reliability of the SWC measurement sensors.



---

*Source: 6665829-2021-04-23.xml*