Next Article in Journal
Multi-Year Mapping of Major Crop Yields in an Irrigation District from High Spatial and Temporal Resolution Vegetation Index
Next Article in Special Issue
Design and Calibration of a Low-Cost SDI-12 Soil Moisture Sensor
Previous Article in Journal
BeeKeeper 2.0: Confidential Blockchain-Enabled IoT System with Fully Homomorphic Computation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Performance Assessment of Five Different Soil Moisture Sensors under Irrigated Field Conditions in Oklahoma

1
Department of Biosystems and Agricultural Engineering, Oklahoma State University, Stillwater, OK 74078, USA
2
Department of Plant and Soil Sciences, Oklahoma State University, Stillwater, OK 74078, USA
3
USDA-ARS Grazinglands Research Laboratory, El Reno, OK 73036, USA
*
Author to whom correspondence should be addressed.
Submission received: 9 October 2018 / Revised: 29 October 2018 / Accepted: 2 November 2018 / Published: 5 November 2018
(This article belongs to the Special Issue Sensors and Systems for Smart Agriculture)

Abstract

:
Meeting the ever-increasing global food, feed, and fiber demands while conserving the quantity and quality of limited agricultural water resources and maintaining the sustainability of irrigated agriculture requires optimizing irrigation management using advanced technologies such as soil moisture sensors. In this study, the performance of five different soil moisture sensors was evaluated for their accuracy in two irrigated cropping systems, one each in central and southwest Oklahoma, with variable levels of soil salinity and clay content. With factory calibrations, three of the sensors had sufficient accuracies at the site with lower levels of salinity and clay, while none of them performed satisfactorily at the site with higher levels of salinity and clay. The study also investigated the performance of different approaches (laboratory, sensor-based, and the Rosetta model) to determine soil moisture thresholds required for irrigation scheduling, i.e., field capacity (FC) and wilting point (WP). The estimated FC and WP by the Rosetta model were closest to the laboratory-measured data using undisturbed soil cores, regardless of the type and number of input parameters used in the Rosetta model. The sensor-based method of ranking the readings resulted in overestimation of FC and WP. Finally, soil moisture depletion, a critical parameter in effective irrigation scheduling, was calculated by combining sensor readings and FC estimates. Ranking-based FC resulted in overestimation of soil moisture depletion, even for accurate sensors at the site with lower levels of salinity and clay.

1. Introduction

Irrigated agriculture, a major contributor to the United States (U.S.) economy, plays a vital role in supplying the demand for food, feed, and fiber. Although only 27% of all croplands in the U.S. are irrigated, this sector is responsible for nearly 50% of crop revenues [1]. Sustaining high levels of food production through irrigated agriculture requires large amounts of water. In 2010, irrigation was the second largest consumer of freshwater withdrawals in the U.S., accounting for approximately 33% (approximately 159 million m3 year−1) of the total water withdrawals [2]. Irrigation water sources, however, are usually limited in amount and are subject to increasing competition. In addition, more variability in precipitation patterns is expected due to climate change, which may threaten the availability of irrigation water supplies [3,4]. These challenges create the need to optimize irrigation management and avoid over- or under-irrigation. Over-irrigation, in addition to wasting water and valuable nutrients, can create favorable conditions for pests and diseases, increase energy costs, and reduce the lifespan of irrigation infrastructure. It can also result in erosion of topsoil and contamination of downstream water resources due to movement of water-soluble chemicals [5]. In contrast, under-irrigation reduces crop yield and negatively impacts economic viability of agricultural production.
Several advanced technologies are available to assist with achieving and implementing optimized irrigation management, including weather stations, air- and spaceborne remote sensing platforms, computer models, plant feedback sensors, and soil moisture sensors [6,7]. Soil moisture sensors, in particular, can be used effectively to improve irrigation management [8]. As a tool for irrigation scheduling, these sensors have been shown to increase crop yields while conserving water [8,9,10,11]. For example, Zotarelli et al. [12] showed that users who manage irrigation with soil moisture sensors applied 15 to 51% less irrigation water compared to fixed-time irrigation plan and observed a crop yield increase of 11 to 26% in Florida, U.S. In addition, sensors can provide continuous estimate of soil moisture conditions in a nondestructive way at a reasonable cost and usually require little maintenance over their lifetime [13]. Soil moisture sensors include tensiometers, neutron gauges, electromagnetic sensors, electrical resistance sensors, and heat dissipation sensors, to name a few [14]. Among these different types, electromagnetic sensors have been widely used by producers for irrigation scheduling.
Despite their numerous advantages, electromagnetic sensors are sensitive to soil salinity and clay content. The impact of soil salinity on sensor readings of soil volumetric water content (θv) (m3 m−3) has been highlighted in several studies [15,16,17]. For example, Wyseure et al. [16] reported that θv error was acceptable at soil bulk electricity conductivity (EC) (dS m−1) levels below 2.0 dS m−1, and Schwartz et al. [18] found that θv estimates were not affected at bulk EC levels below 2.8 dS m−1. These thresholds are exceeded in many irrigated areas in arid/semi-arid regions, where there is a great need for improving irrigation management using sensor technologies. The results from prior studies on the impact of clay content have been somewhat variable. Rüdiger et al. [19] observed overestimation error in θv that increased with clay content. In contrast, Fares et al. [20] observed underestimation of θv for electromagnetic sensors due to high clay content, which was more prevalent at lower soil moisture content. Mittelbach, et al. [21] reported both under- and overestimation errors at different depths of a clay loam soil in Switzerland. In light of these variable results, and since high salinity and clay content conditions are encountered in many agricultural fields, there is a need to undertake further field studies to investigate the accuracy of electromagnetic soil moisture sensors under varying levels of salinity and clay content.
The goal of this study was to evaluate the performance of soil moisture sensors for irrigation scheduling purposes under low and high salinity/clay content conditions. Specific objectives were to (1) assess the performance of five different commercially available electromagnetic sensors in estimating θv in situ under soils with variable salt and clay content, (2) compare the accuracy of several approaches of determining soil moisture thresholds used in irrigation scheduling, and (3) investigate the accuracy of estimated soil moisture depletion based on sensor readings and different threshold approaches.

2. Materials and Methods

2.1. Sensor Description

Five commercially available electromagnetic sensors were evaluated in this study: TDR315, CS655, GS1, SM100, and CropX.

2.1.1. TDR315

The TDR315 (Acclima Inc., Meridian, ID, USA) is a recently commercialized sensor for agricultural applications [18]. This sensor operates on principles of Time Domain Reflectometry (TDR) that estimates the soil apparent permittivity (Ka) (unitless) at relatively higher frequencies (3.5 GHz), which are less sensitive to bulk EC compared to lower frequency electromagnetic techniques [22]. Conventional TDR sensors have a problem sustaining high frequency signals because of signal attenuation in the sensor’s coaxial cables. The TDR315 addresses this issue by embedding all the electronics required for pulse generation and waveform acquisition in a compact circuit within the probe head. The data are transmitted digitally via SDI-12 (Serial Data Interface at 1200 baud) protocol, which is an asynchronous, ASCII, serial communications protocol and can support a cable length of up to 60 m. The sensor shares the same advantages of the conventional TDRs, but, it is more portable, affordable, and convenient to use [18]. These sensors have a planar three-conductor transmission line, each 15 cm long, and transmit the incident pulse in the center rod and two exterior grounds. The TDR315 reports volumetric water content (θv) (%) based on a proprietary dielectric mixed model which estimates Ka using Topp equation (Equation (1)) [17]. The sensor also reports soil temperature (°C), bulk relative permittivity (unitless), bulk EC (µS cm−1), and soil pore water EC (µS cm−1). DataSnap SDI-12 data-loggers from the same manufacturer were used with TDR315 sensors to collect data on hourly basis.
θ v = 4.3 × 10 6 × ( K a 3 ) 5.5 × 10 4 × ( K a 2 ) + 2.92 × 10 2 × ( K a ) 5.3 × 10 2

2.1.2. CS655

The CS655 sensor (Campbell Scientific, Inc., Logan, UT, USA) is a water content reflectometer. An electronic pulse is sent from the probe head and reflected at the end of the rods (12 cm in length). Upon detecting the returned pulse, another pulse is sent. Then, the probe records the frequency of these pulses and inverses the frequency as period in microseconds (µs). This period is impacted by the velocity of electromagnetic pulse, which is influenced by Ka [17,23]. The probe estimates θv from Ka using the Topp equation (Equation (1)) [17]. Apart from θv, other measured parameters include period average (µs), soil’s relative dielectric permittivity (unitless), bulk EC (dS m−1), and soil temperature (°C). Like the TDR315, the CS655 communicates with a data-logger using an SDI-12 interface. To collect hourly data, CR1000 data-loggers (Campbell Scientific, Inc., Logan, UT, USA) were used in this study.

2.1.3. GS1

The GS1 sensor (METER Group, Inc., Pullman, WA, USA) estimates θv by generating an electromagnetic field to measure the dielectric constant of the surrounding medium. This sensor uses capacitance and frequency domain technology and operates at 70 MHz [24]. It provides oscillating waves to the sensor rods that charge in response to the dielectric of the material. The sensor quantifies the charge and provides a raw value (RV) that is strongly correlated with θv (Equation (2)). The GS1 has a rugged design and is capable of remaining in the soil for a long time. It has a two-rod design, with each rod measuring 5.5 cm in length. Hourly data were collected throughout the cropping season using EM5B analog data-loggers (METER Group, Inc., Pullman, WA, USA).
θ v = 3.62 × 10 4 ( RV ) 0.554

2.1.4. SM100

The WaterScout SM100 sensor (Spectrum Technologies, Aurora, IL, USA) has two electrodes functioning as a capacitor, with surrounding soil acting as the dielectric. The capacitor is driven by an 80 MHz oscillator and converts the soil’s dielectric permittivity to an output signal, which correlates with θv. Watchdog 1400 data-loggers from the same manufacturer were used to collect hourly data.

2.1.5. CropX

The CropX sensor (CropX Ltd., Tel Aviv, Israel) integrates soil moisture sensing and a cellular communication package. The sensor electrodes are built into a helical wing attached to a central shaft for installation with reduced soil disturbance. The sensor measures soil moisture based on the amplitude domain reflectometry. When the amount of water changes in the soil, the sensor measures the change in amplitude differential due to changes in dielectric permittivity, which directly correlates to changes in water content. CropX is a multiprobe sensor that measures θv at 20 and 46 cm depths in the soil.

2.2. Study Sites

The study took place during the 2017 crop growing season. Two sites were selected for sensor installation, one with lower salinity and lower clay content (LSLC) located in central Oklahoma and the other in southwest Oklahoma with higher salinity and higher clay content (HSHC). Figure 1 shows the location of the study sites overlaid on the map of long-term mean annual precipitation across Oklahoma, obtained from Daly et al. [25]. The LSLC site had a Pond Creek fine sandy loam soil (fine-silty, mixed, superactive, thermic Pachic Argiustolls) while the HSHC site had a Hollister silty clay loam soil (fine, smectitic, thermic Typic Haplusterts). The EC of the soil solution (1:1 soil–water ratio) was 1.2 dS m−1 at LSLC compared to 7.0 dS m−1 at HSHC. Table 1 provides additional information on soil characteristics at each site. In addition to variations in soils, the two sites were different in crop types, irrigation systems, and climatic conditions. Corn (Zea mays L.) was planted at the LSLC site under a center-pivot irrigation system, while the HSHC site was under furrow-irrigated cotton (Gossypium hirsutum L.). Key meteorological parameters for each site are given in Table 2.

2.3. Experimental Setup

Four replications of TDR315, CS655, GS1, and SM100 and two replications of CropX were installed on 7/20/2017 and 7/27/2017 at LSLC and HSHC sites, respectively. The sensors were used with manufacturer-provided data-loggers and calibrations because the results obtained in this manner would best represent the conditions that irrigators and farm managers would face in the field [26]. Therefore, the raw θv readings reported by the sensors were used in analysis without any alteration [13]. It should be noted that developing and utilizing site-specific calibrations can significantly improve accuracies if the required technical and financial resources are available to users. All sensors were installed at a depth of 20 cm from the soil surface. The top 20 cm is important for plant water uptake as root distribution of plants is denser in this layer than deeper in the soil profile [27].
At each replication, a pit was dug between two rows of crops to install the soil moisture sensors. Physical properties of soil in each pit were determined in the Soil Physics Laboratory at Oklahoma State University (OSU) by taking undisturbed soil cores (diameter = 2.5 cm, length = 5.1 cm) using the Sample Ring Kit (Model C, Eijkelkamp Soil & Water, Inc., Giesbeek, The Netherlands) on the day of sensor installation. Soil textural information (particle size distribution) were determined by hydrometer following the protocol proposed by Ashworth, et al. [28]. Additionally, four replications of soil samples were taken at each site on the installation day to measure soil salinity. The salinity test was done by Soil, Water and Forage Analytical Laboratory at OSU using the 1:1 soil water extraction method [29].
Sensors were inserted horizontally into the side wall of the pit (undisturbed soil) so that the rods of the sensors were on top of each other (vertical orientation) and the middle point of the sensor rods was directly under the crop row. The θv readings are often impacted by the sensor installation procedure [30], so extra care was taken to maintain minimal disturbance to the surrounding soil while inserting the rods. The spacing between the sensors was determined based on the volume of influence of individual sensors plus an additional distance to eliminate any possible interference. This spacing between the sensors were varied from 10 to 18 cm depending on the volume of influence of sensors’ electromagnetic field. Then, wires were run below and away from the sensors for some distance to avoid creating any preferential flow channels. After that, the wires were run through PVC pipes to the data-logger encasement. The CropX sensors were installed using the spiral auger provided by the manufacturer to minimize soil disturbance.
The excavated soil was collected in different buckets for different soil layers and carefully used to backfill the pits, attempting to recreate the original bulk density. Precipitation amounts were recorded by a tipping bucket rain gage (model TE525-L, Campbell Scientific, Inc., Logan, UT, USA) at the LSLC site, whereas, these measurements were collected from an Oklahoma Mesonet weather station located 678 m to the southwest of the sensor installation location at the HSHC site [31]. Gravimetric soil samples (diameter = 3 cm, height = 5.1 cm) were collected using a Giddings soil sampling probe (Giddings Machine Company, Windsor, CO, USA) to estimate reference θv (θref) (m3 m−3) throughout the crop growing season. On each sampling date, four gravimetric samples were taken at each site and the probe was centered at the sensor installation depth (20 cm). If there was an irrigation and/or precipitation event around the sampling dates, extra care was taken not to compact the areas above the sensors. Soil samples were put in plastic bags immediately after sample collection and kept out of sunlight to minimize evaporation. All soil samples of known volumes were oven-dried at 105 °C for 24 h and used to determine bulk density.

2.4. Soil Moisture Thresholds

Efficient irrigation management requires knowledge of two important soil moisture thresholds that indicate water availability for plant consumption [5]. These thresholds are field capacity (FC) and wilting point (WP). The FC is often estimated as the water retained at a soil matric potential of −33 kPa, although research has shown that this can result in underestimation of FC and −10 kPa may provide a more suitable approximation [32]. The WP is often estimated as the water retained at −1500 kPa [33]. These values can be different depending on soil texture, crop type, and other factors.
In this study, FC and WP were determined using three different approaches: laboratory, sensor-based, and the Rosetta model [34]. Undisturbed soil cores extracted from each site were used in laboratory tests where FC was determined at −33 kPa using the Tempe cell method and WP at −1500 kPa using the pressure plate method [35]. The sensor-based approach was based on ranking of the collected data following the procedure proposed in Hunt et al. [36]. This method uses sensor readings to estimate FC and WP as the 95th and the 5th percentiles of all θv values collected during the study period. This method assumes that the hydrologic conditions during the measurement period result in θv values, which span from values lower than WP to values higher than FC. The Rosetta model uses hierarchical pedotransfer functions to estimate van Genuchten water retention parameters [34]. In this study, three different FC-WP outputs were generated from the Rosetta model by providing different types and combination of input data. The three types of input data included (i) only the textural class of soils at study sites, (ii) textural information (percentages of sand, silt, and clay), and, (iii) textural information and bulk density. Estimated FC and WP from all methods described above were compared with those reported in the U.S. Department of Agriculture’s Web Soil Survey at each study site [37]. In addition to FC and WP, the available water content (AWC), which is the difference between FC and WP, was calculated and compared with values obtained from different methods described above [38].
To optimize irrigation management based on soil moisture sensing, sensor readings must be converted to soil moisture depletion (SMD) (m3 m−3). In this study, SMD was calculated as the difference between FC and θv:
S M D ( i ) = θ F C θ v ( i )
where, SMD(i) is the soil moisture depletion at the ith time-step, θFC is the θv at FC (constant) (m3 m−3), and θv(i) is the θv at the ith time-step. In estimating SMD, θv(i) values were obtained from sensor readings and θFC values were based on two different approaches, resulting in two SMD estimates for each sensor at each site. The two θFC approaches were the laboratory and the ranking methods explained above. The results were compared against SMD estimates based on θref (gravimetric measurements) and laboratory θFC. After SMD is estimated, it can be multiplied by the root zone depth to obtain an estimate of irrigation requirement in units of water depth.

2.5. Statistical Analysis

To evaluate the performance of the selected sensors, θv readings of sensors were compared with θref values. Four statistical parameters, namely root mean square error (RMSE), RMSE-observations standard deviation ratio (RSR), mean bias error (MBE), and index of agreement (k) were estimated according to the following equations.
RMSE = 1 n i = 1 n ( P i O i ) 2
RSR = RMSE STDEV O ( i ) = i = 1 n ( O i P i ) 2 i = 1 n ( O i O ¯ ) 2
MBE = 1 n i = 1 n ( P i O i )
k = 1 [ i = 1 n ( P i O i ) 2 i = 1 n ( | P i O ¯ | + | O i O ¯ | ) 2 ]
where, n is the sample size, i is the index of sample pairs, P is the sensor reading (predicted), O is the θref (observed), and O ¯ is the mean of all θref values.
The accuracy categories outlined in Fares et al. [20] were adopted in this study for interpreting RMSE values. These categories include good (RMSE ≤ 0.01 m3 m−3), fair (0.01 ≤ RMSE ≤ 0.05 m3 m−3), poor (0.05 ≤ RMSE ≤ 0.10 m3 m−3), and very poor (RMSE ≥ 0.10 m3 m−3). The RSR provides benefits of incorporating error index statistics and it includes a normalization factor applicable to various constituents [39]. The RSR varies from a value of zero indicating zero RMSE and a perfect model simulation to a large positive value. The performance of a model is determined by different categories of RSR: very good model fit (0.00 ≤ RSR ≤ 0.50), good model fit (0.50 ≤ RSR ≤ 0.60), satisfactory model fit (0.60 ≤ RSR ≤ 0.70), and unsatisfactory model fit (RSR > 0.70). However, these categories are based on simulations running on a monthly time-step. Moriasi et al. [39] noted that the acceptable range of RSR would increase in magnitude when using smaller time-steps, which was the case in this study. The MBE measures the average difference between sensor-estimated θv and θref. A MBE of zero indicates the predicted and observed values are unbiased. A positive value of MBE means sensor is overestimating θv, and negative MBE indicates underestimation [40]. The index of agreement (k) was used to determine how well the sensor-estimated θv agreed with θref [41]. The value of k can range from zero to one, with one representing the highest level of agreement and zero representing complete disagreement [42].
In addition to the above statistical parameters, Pearson correlation coefficients (r) were calculated for pairwise sensor comparisons to evaluate the similarity in their temporal variations throughout the study period. Closely correlated temporal patterns have a r value near one, while this parameter is near zero in case of uncorrelated patterns [43]. Finally, linear regression models were fitted to sensor-estimated θv and θref using the Minitab statistical software (version 17.3) (Minitab, Inc., State College, PA, USA) [44]. These linear models and the reported intercepts and slopes for each sensor can be used as field calibration equations in future applications at the study sites.

3. Results and Discussion

3.1. Sensor Performance

The fluctuations in θv were similar across all sensors at both study sites (Figure 2). All sensors responded to most irrigation and precipitation events. In some cases, there was little or no change in θv following a watering event, mainly because the amount of water received was not large enough to reach sensor installation depth. The results of performance evaluation (statistical indicators) are summarized in Table 3. In general, all sensors performed better at the LSLC. At this site, the RMSE was the lowest for CS655 (0.019 m3 m−3), followed by TDR315 (0.028 m3 m−3) and GS1 (0.048 m3 m−3). These values belong to the fair accuracy category defined in Fares et al. [20], suggesting that CS655, TDR315, and GS1 can be implemented for effective irrigation scheduling under conditions similar to those of LSLC. The RMSE values obtained in this study were smaller than the RMSE values of 0.105 and 0.049 m3 m−3 reported by Singh et al. [45] for the CS655 and TDR315 in a loam soil, respectively. Adeyemi et al. [46] found a similar RMSE of 0.020 m3 m−3 for TDR315 and 0.050 m3 m−3 for GS1 in a sandy loam soil under laboratory conditions. The RMSE of CropX was 0.051 m3 m−3, which is in the poor category. The SM100’s RMSE was very poor (0.110 m3 m−3).
The MBE and RSR revealed similar patterns in sensor performance at the LSLC, with the CS655 performing the best, followed by TDR315, GS1, CropX, and SM100. The MBE indicated that all sensors overestimated θv at LSLC. This overestimation can also be observed in Figure 3 as most of the points were above the 1:1 line. Overestimation of θv by CS655 was observed by Kisekka et al. [47] and Michel et al. [48] too. Adeyemi et al. [46] found that TDR315 and GS1 underestimated θv in sandy loam soil, but, with increasing clay content, the underestimation became overestimation. The RSR ranged from 0.53 for CS655 to 3.00 for SM100 at LSLC site. According to categories defined by Moriasi et al. [39], the CS655 had a good model fit whereas all other sensors were classified as having unsatisfactory model fit. But as mentioned previously, running a model on temporal resolution higher than monthly would warrant less strict performance rating. Therefore, higher RSR values are expected in this study because of hourly time-step analysis. This trend was also observed in a study by Wyatt et al. [49], which produced high RSR values at daily time-step.
All sensors had larger RMSE at the HSHC site compared to LSLC (Table 3). However, the magnitude of the increase in RMSE was not uniform and changed from a slight increase for CropX to over an eight-fold increase for CS655. The CropX sensor had the smallest RMSE, followed by TDR315, GS1, CS655, and SM100. The values of RMSE belonged to the poor accuracy category in case of CropX and TDR315 and very poor category for other sensors according to classifications in Fares et al. [20], suggesting that none of the sensors can be implemented for effective irrigation scheduling under conditions similar to those of HSHC. In addition, the variability of readings among the replications of the same sensors increased at HSHC; the average standard deviation (SD) ranged from 0.021 m3 m−3 for TDR315 to 0.050 m3 m−3 for CS655. At LSLC, the average SD varied from 0.011 m3 m−3 for TDR315 to 0.023 m3 m−3 for GS1. The average SD of θref was 0.015 m3 m−3 at LSLC and 0.010 m3 m−3 at HSHC.
High clay content and elevated levels of salinity seem to be the main reasons behind lower sensor accuracies at the HSHC site. Adeyemi et al. [46] concluded that the errors in TDR315 and GS1 would increase with an increase in soil salinity level. In addition, Wyseure et al. [16] reported that the error in TDR sensors would remain within reasonable limits if the bulk EC is kept less than 2 dS m−1. The bulk EC at HSHC, however, was well over this threshold. The MBE estimates were larger at HSHC than LSLC and showed that all sensors except CropX overestimated θv. This is also evident in Figure 3. Most of previous studies have reported overestimation error for TDR sensors under saline conditions. This is mainly due to the fact that in saline soils, the dielectric permittivity measured by TDR increases and therefore θv is overestimated as mentioned in Dalton [15]. However, Schwartz et al. [18] found that TDR315 underestimated θv in a saline Pullman clay loam soil. The RSR values followed a pattern similar to other error indicators at HSHC, having the smallest value of 1.34 for CropX and the largest value of 5.66 for SM100.
Some noise in θv readings of the TDR315 at HSHC can be seen in Figure 2b. Schwartz et al. [18] reported that TDR315 sensors were insensitive to bulk EC up to 2.8 dS m−1 and corresponding pore water EC up to 7.3 dS m−1. The bulk EC and pore water EC exceeded these thresholds at HSHC on many days at the beginning of the study period. This might have caused signal attenuation that induced noise in θv readings. This noise was quantified using standard deviation (SD) in θv among the replications. At the beginning of the growing season, the SD had a range of zero to 0.099 m3 m−3 and an average of 0.021 m3 m−3 at HSHC for TDR315, which is much larger when compared to the range of 0.002 to 0.043 m3 m−3 and average of 0.012 m3 m−3 at LSLC during the same period. The observed noise was reduced later in the growing season, probably due to decrease in soil EC because of leaching of salts by irrigation water.
The hourly bulk EC estimates from TDR315 and CS655 were in agreement with soil EC determined in the laboratory and showed the significant difference between the two study sites (Figure 4). Both sensors reported small bulk EC at LSLC with similar ranges of 0.1 to 0.4 dS m−1. At HSHC, however, the bulk EC was significantly larger with ranges of 1.1 to 3.4 and 0.9 to 3.0 dS m−1 based on TDR315 and CS655 sensors, respectively.
In utilizing soil moisture sensors for irrigation management, obtaining a complete time series is as important as taking accurate readings. In this study, CropX and CS655 had significant data gaps for different reasons. On average, 41% of the CropX data were missing at LSLC compared to less than one percent at HSHC. Several correspondences with the manufacturer revealed that the potential reason behind this issue could be the tall corn canopy at LSLC, which can block the transmitted signals. Upon recommendation from the manufacturer, extension antennas were installed on CropX sensors at LSLC. The observed crop height was 2.16 m and the extension antennas were installed in such a way that the tops of the antennae were 1.91 m from the ground. However, this modification did not help with the apparent transmission problem.
The CS655 had 21% missing data at HSHC. Sugita et al. [50] conducted a reliability test on CS655 and found that the sensor was missing 64% of the measurements when exposed to high salinity levels (bulk EC = 1.2–2.1 dS m−1). The bulk EC at HSHC was larger than the values reported in Sugita et al. [50]. In addition to high salinity, the HSHC site had relatively high clay content (38.7%). The clay particles have highly charged surface areas which increase dielectric losses and cause the apparent permittivity (Ka) values to go outside the acceptable range of Topp equation [17]. The combined effect of higher soil salinity and clay content results in the attenuation of the electromagnetic signal from the sensor [18].Therefore, the sensor fails to report θv in case of Ka 42 and θv ≥ 0.52 m3 m−3 as the internal logical test rejects these data.
Linear regression equations were developed to estimate θref based on sensor-estimated θv (Table 4). These equations can be used to get more accurate θv readings in areas matching this study’s local conditions. At LSLC, the regression models were all statistically significant at α = 0.05, with r2 values ranging from 0.57 for CropX to 0.85 for CS655. Although SM100 had low accuracy, the high r2 value (0.84) indicates that this sensor had high degree of correlation with the reference values. At the HSHC site, the linear regression model for CS655 was not statistically significant. Models of other sensors were significant and had r2 values varying from 0.73 to 0.85.

3.2. Correlations between Sensors

In general, the Pearson’s correlation coefficients (r) of θv readings were larger at LSLC than HSHC (Table 5). At this site, the strongest correlation (r = 0.99) was between TDR315 and CS655 and the weakest was between CropX and SM100 (r = 0.79). The correlation coefficients for CropX were smallest among all sensors at the LSLC site, ranging from 0.79 to 0.81. Despite being the least accurate sensor, SM100 had strong correlation with the top two accurate sensors, i.e., TDR315 and CS655. This indicates that SM100 closely followed the temporal changes in θv of more accurate sensors. At HSHC, the correlation between TDR315 and GS1 was the strongest (r = 0.97). The SM100 also had strong correlations with TDR315, GS1, and CropX. On the other hand, CS655 had weak correlations with other sensors.
The strong correlation between sensors with different accuracies suggests that the response of less accurate sensors to soil moisture fluctuations was similar to those of more accurate sensors. The differences in θv readings were relatively constant over the study period (offset error). This provides an opportunity for potential utilization of less accurate sensors in some limited applications where the user is only interested in determining the movement of the water front in the soil profile. One example of this application is leaching salts below the root zone. In this case, the user needs to ensure water front has moved below the bottom of the root zone. Another example is preventing deep percolation to ensure applied water remains within the root zone and that soluble chemicals are not transported to shallow groundwater resources.

3.3. Soil Moisture Thresholds

At LSLC, the FC and WP estimated in the laboratory were similar to the output of the Rosetta model based on textural class, textural information, and textural information plus bulk density (Table 6). Thresholds obtained from USDA’s Web Soil Survey (USDA-WSS) were slightly larger than the results of the laboratory and Rosetta methods. However, the estimates based on the ranking of sensor readings were significantly larger than those of the other methods. The FC and WP values were larger at HSHC compared to LSLC irrespective of the method used because of larger clay content in the soils. The FC values from the Rosetta model and the USDA-WSS were either similar or slightly smaller than those obtained with the laboratory approach. All ranking estimates of FC were significantly larger than those with laboratory approach except for CropX, which was slightly larger. In the case of WP, estimates from the Rosetta model were significantly smaller than those with the laboratory approach, while USDA-WSS reported a similar value. Ranking method estimates were significantly larger except for CropX. The differences between AWC estimates of the ranking and laboratory methods were smaller than the differences in the FC and WP estimates of the same methods at both sites, mainly because overestimations in FC and WP estimates of the ranking method were of similar magnitudes and thus cancelled out to a large extent.
Results of this study reveal that the Rosetta model is capable of accurately estimating soil moisture thresholds even with minimal input data (textural classes). The USDA-WSS also performed satisfactorily, despite the fact that it is based on coarse soil surveys. However, the ranking method resulted in significant overestimation of FC when compared to laboratory estimates, ranging from 59 to 117% at the LSLC and from 6 to 94% at HSHC site. The difference between WP estimates of the ranking and laboratory methods varied from 100 to 283% at LSLC and from −14 to 129% at HSHC. A potential reason behind this poor performance could be that the full range of soil moisture conditions was not experienced at both sites during the period of study. However, this situation could be the case in many irrigated areas, since producers attempt to replenish soil moisture well before it reaches WP to avoid water stress and yield loss. Another reason behind the poor performance of the ranking method is the error in sensor readings, especially at HSHC, where most sensors overestimated soil moisture due to high clay content and elevated salinity levels.
Variations in hourly SMD are presented in Figure 5. In this figure, dots represent observed SMD based on θref and laboratory-determined FC, while lines represent sensor SMD based on sensor θv and FC from two methods: laboratory and ranking. The Rosetta model was not considered here because the FC values obtained from the model were similar to those of the laboratory. At LSLC, observed SMD values were zero except on two sampling dates in early September. This is because this site was under full to slightly over-irrigation at most times during the study period. The only exception for the same period was in September when crop water demand outpaced irrigation application. Possible underestimation of θFC in the laboratory method may have contributed to zero SMD on most measurement dates too. In this study, a soil matric potential of −33 kPa was used to measure θFC. But as mentioned before, this value can be as high as −10 kPa in sandy loam soil, resulting in a larger θFC and consequently a larger SMD estimate. Sensor SMDs based on laboratory-FC had similar patterns, indicating no depletion during the study period except in the month of September (Figure 5a). On the other hand, sensor SMDs based on ranking-FC showed significant depletions at most times, reaching values as large as 0.15 m3 m−3 (Figure 5b). This increase in SMD is mainly due to overestimation of FC in the ranking method, since the same sensors readings were used in both SMD approaches.
At the HSHC site, the observed SMD indicated a larger depletion, especially during early September to early October. This pattern was expected since this site was under a low-frequency (7–10 days) flood irrigation regime that was not able to meet cotton water demand during the hot and dry month of September. At this site, sensor SMDs based on laboratory-FC showed no depletion except for CropX and TDR315. The SMD estimates of CropX were larger and the SMD estimates of TDR315 were smaller than observed SMD. This is because CropX underestimated θv, while TDR315 overestimated this parameter. The overestimation errors of the other sensors were so large that their θv readings were above laboratory-FC at all times, resulting in no depletion. The sensor SMDs based on ranking-FC were significantly larger than those based on laboratory-FC, except for CS655. This was because of the overestimation of FC by the ranking method. Hence, depletion was calculated at most times. The SMDs of CS655 were similar to the observed SMD, since the overestimation errors in θv readings and ranking-FC were similar in magnitude.

4. Conclusions

The performance of five types of commercially available soil moisture sensors was evaluated at two fields with significantly different salinity levels and clay contents. The sensors included TDR315, CS655, GS1, SM100, and CropX. The accuracy of each sensor was determined by comparing its readings with gravimetric measurements of soil water content obtained at several times during the study period. In general, all sensors responded to wetting and drying events. The TDR315, CS655, and GS1 sensors had acceptable accuracies for managing irrigations at the site with low salinity and low clay content (LSLC) based on root mean square error (RMSE). However, none of the sensors performed satisfactorily at the site with high salinity and high clay content (HSHC), with RMSE estimates that were up to eight times larger compared to the values at LSLC. In addition, high levels of noise were observed in TDR315 due to high salinity and out-of-range responses and consequently missing readings in case of CS655 sensor. A potential solution for using soil moisture sensors in irrigation scheduling under such conditions is the use of site-specific calibrations.
For practical irrigation scheduling, sensor readings must be used in conjunction with soil moisture thresholds of field capacity (FC) and wilting point (WP) in order to estimate soil moisture depletion (SMD) and consequently irrigation requirement. In this study, FC and WP values determined in the laboratory using undisturbed soil cores were compared against those obtained from three independent approaches: the Rosetta model, the ranking of sensor readings, and the values reported in the U.S. Department of Agriculture’s Web Soil Survey (USDA-WSS). The Rosetta model was capable of providing estimates similar to those of the laboratory approach, regardless of the type and number of input data used in the model. The USDA-WSS approach resulted in acceptable estimates of FC and WP. The ranking method, however, significantly overestimated FC and WP at both sites, even for accurate sensors. The ranking method did not perform well in estimating SMD either, except for one sensor at the HSHC site where the overestimation error in FC was similar to overestimation error in soil water content and canceled each other out. The results of this study show that two major conditions are required before the ranking method can be used effectively in estimating soil moisture thresholds: sensor readings that are the basis of calculations must be accurate; and, the full range of moisture conditions from below WP to above FC must be experienced during the data collection period.
This study contributes to the existing knowledge on sensor-based irrigation scheduling through quantifying the accuracies of five widely-used soil moisture sensors as impacted by soil clay content and salinity, as well as investigating the effectiveness of different soil moisture threshold estimation approaches for agricultural irrigation applications. The results highlight the wide range of accuracies that exist among soil moisture sensors and methods for determining soil moisture thresholds. Such a wide range creates major challenges in utilizing soil moisture sensors for irrigation scheduling applications. As new sensors are being developed frequently, studies like this need to be conducted under variable field conditions to evaluate the performance of the new sensors and to provide guidelines on how they can be used for irrigation scheduling purposes.

Author Contributions

S.D. conducted the study, collected and analyzed field data, and wrote the paper. S.T. conceived of this study, collected field data, and helped in analysis and writing the paper. T.E.O. helped in designing and planning the experimental setup of the study and provided valuable comments in writing and revising the paper. D.M., P.G., and J.L.S. were instrumental in acquiring instrumentation and helped in writing and revising the paper.

Funding

This manuscript is based upon work supported by the Natural Resources Conservation Service, U.S. Department of Agriculture, under No. 69-3A75-16-013. Additional funding was provided by the Agricultural Research Service, U.S. Department of Agriculture, under Agreement No. 58-3070-5-007 and by the U.S. Geological Survey 104 (b) program administered through the Oklahoma Water Resources Center.

Commercial Endorsement Disclaimer

The use of trade, corporation, and firm names in this article is for informational purpose only and is at convenience of the reader. Oklahoma State University and US Department of Agriculture neither approve nor endorse the use of products from these companies.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. United States Department of Agriculture (USDA). Farm and Ranch Irrigation Survey (2013); National Agricultural Statistics Service, Ed.; USDA: Washington, DC, USA, 2014; Volume 3, p. 266.
  2. Maupin, M.A.; Kenny, J.F.; Hutson, S.S.; Lovelace, J.K.; Barber, N.L.; Linsey, K.S. Estimated Use of Water in the United States in 2010; US Geological Survey: Reston, VA, USA, 2014.
  3. Fischer, G.; Tubiello, F.N.; Van Velthuizen, H.; Wiberg, D.A. Climate change impacts on irrigation water requirements: Effects of mitigation, 1990–2080. Technol. Forecast. Soc. Chang. 2007, 74, 1083–1107. [Google Scholar] [CrossRef] [Green Version]
  4. Fishman, R.M. Climate Change, Rainfall Variability, and Adaptation through Irrigation: Evidence from Indian Agriculture. Unpublished work. 2012.
  5. Datta, S.; Taghvaeian, S.; Stivers, J. Understanding Soil Water Content and Thresholds for Irrigation Management; Oklahoma Cooperative Extension Service: Stillwater, OK, USA, 2017. [Google Scholar]
  6. Broner, I. Irrigation scheduling. In Service in Action; No. 4.708; Colorado State University Extension: Fort Collins, CO, USA, 2005. [Google Scholar]
  7. Martin, E. Methods of Measuring for Irrigation Scheduling—WHEN; Arizona Cooperative Extension: Tucson, AZ, USA, 2009. [Google Scholar]
  8. Martin, E.C.; Pegelow, E.J.; Stedman, S. Comparison of Irrigation Scheduling Methods in Cotton Production; College of Agriculture, University of Arizona: Tucson, AZ, USA, 1995. [Google Scholar]
  9. Fisher, D.K.; Hanks, J.E.; Pringle, H.L., III. Comparison of Irrigation Scheduling Methods in the Humid Mid-South. Irrigation Association. 2009. Available online: https://www.irrigation.org/IA/FileUploads/IA/Resources/TechnicalPapers/2009/ComparisonOfIrrigationSchedulingMethodsInTheHumidMid-South.pdf (accessed on 21 August 2018).
  10. Sui, R. Irrigation Scheduling Using Soil Moisture Sensors. J. Agric. Sci. 2017, 10, 1. [Google Scholar] [CrossRef]
  11. Kebede, H.; Fisher, D.K.; Sui, R.; Reddy, K.N. Irrigation Methods and Scheduling in the Delta Region of Mississippi: Current Status and Strategies to Improve Irrigation Efficiency. Am. J. Plant Sci. 2014, 5, 50005. [Google Scholar] [CrossRef]
  12. Zotarelli, L.; Scholberg, J.M.; Dukes, M.D.; Muñoz-Carpena, R.; Icerman, J. Tomato yield, biomass accumulation, root distribution and irrigation water use efficiency on a sandy soil, as affected by nitrogen rate and irrigation scheduling. Agric. Water Manag. 2009, 96, 23–34. [Google Scholar] [CrossRef]
  13. Cardenas-Lailhacar, B.; Dukes, M. Precision of soil moisture sensor irrigation controllers under field conditions. Agric. Water Manag. 2010, 97, 666–672. [Google Scholar] [CrossRef]
  14. Yoder, R.; Johnson, D.; Wilkerson, J.; Yoder, D. Soilwater sensor performance. Appl. Eng. Agric. 1998, 14, 121–133. [Google Scholar] [CrossRef]
  15. Dalton, F. Development of time-domain reflectometry for measuring soil water content and bulk soil electrical conductivity. In Advances in Measurement of Soil Physical Properties: Bringing Theory into Practice; Topp, G.C., Daniel, W.D., Green, R.E., Eds.; Soil Science Society of America: Madison, WI, USA, 1992; pp. 143–167. [Google Scholar]
  16. Wyseure, G.; Mojid, M.; Malik, M. Measurement of volumetric water content by TDR in saline soils. Eur. J. Soil Sci. 1997, 48, 347–354. [Google Scholar] [CrossRef]
  17. Topp, G.C.; Davis, J.; Annan, A.P. Electromagnetic determination of soil water content: Measurements in coaxial transmission lines. Water Resour. Res. 1980, 16, 574–582. [Google Scholar] [CrossRef]
  18. Schwartz, R.C.; Evett, S.R.; Anderson, S.K.; Anderson, D.J. Evaluation of a Direct-Coupled Time-Domain Reflectometry for Determination of Soil Water Content and Bulk Electrical Conductivity. Vadose Zone J. 2016, 15, 8. [Google Scholar] [CrossRef]
  19. Rüdiger, C.; Western, A.W.; Walker, J.P.; Smith, A.B.; Kalma, J.D.; Willgoose, G.R. Towards a general equation for frequency domain reflectometers. J. Hydrol. 2010, 383, 319–329. [Google Scholar] [CrossRef]
  20. Fares, A.; Abbas, F.; Maria, D.; Mair, A. Improved calibration functions of three capacitance probes for the measurement of soil moisture in tropical soils. Sensors 2011, 11, 4858–4874. [Google Scholar] [CrossRef] [PubMed]
  21. Mittelbach, H.; Lehner, I.; Seneviratne, S.I. Comparison of four soil moisture sensor types under field conditions in Switzerland. J. Hydrol. 2012, 430, 39–49. [Google Scholar] [CrossRef]
  22. Robinson, D.; Jones, S.B.; Wraith, J.; Or, D.; Friedman, S. A review of advances in dielectric and electrical conductivity measurement in soils using time domain reflectometry. Vadose Zone J. 2003, 2, 444–475. [Google Scholar] [CrossRef]
  23. Chávez, J.L.; Varble, J.L.; Andales, A.A. Performance evaluation of selected soil moisture sensors. In Proceedings of the 23rd Central Plains Irrigation Conference, Burlington, CO, USA, 22–23 February 2011; pp. 29–38. [Google Scholar]
  24. Environment, M. GS1: Soil Moisture Sensor Operations Manual; METER Group, Inc.: Pullman, WA, USA, 2015; p. 22. [Google Scholar]
  25. Daly, C.; Halbleib, M.; Smith, J.I.; Gibson, W.P.; Doggett, M.K.; Taylor, G.H.; Curtis, J.; Pasteris, P.P. Physiographically sensitive mapping of climatological temperature and precipitation across the conterminous United States. Int. J. Climatol. 2008, 28, 2031–2064. [Google Scholar] [CrossRef] [Green Version]
  26. Leib, B.G.; Jabro, J.D.; Matthews, G.R. Field evaluation and performance comparison of soil moisture sensors. Soil Sci. 2003, 168, 396–408. [Google Scholar] [CrossRef]
  27. Brutsaert, W. Daily evaporation from drying soil: Universal parameterization with similarity. Water Resour. Res. 2014, 50, 3206–3215. [Google Scholar] [CrossRef]
  28. Ashworth, J.; Keyes, D.; Kirk, R.; Lessard, R. Standard procedure in the hydrometer method for particle size analysis. Commun. Soil Sci. Plant Anal. 2001, 32, 633–642. [Google Scholar] [CrossRef]
  29. Zhang, H.; Kress, M.F.; Johnson, G.V. Procedures Used by OSU Soil, Water, and Forage Analytical Laboratory; Division of Agricultural Sciences and Natural Resources, Oklahoma State University: Stillwater, OK, USA, 2002. [Google Scholar]
  30. Aguilar, J.; Rogers, D.; Kisekka, I. Irrigation Scheduling Based on Soil Moisture Sensors and Evapotranspiration. Kansas Agric. Exp. Stn. Res. Rep. 2015, 1, 20. [Google Scholar] [CrossRef] [Green Version]
  31. McPherson, R.A.; Fiebrich, C.A.; Crawford, K.C.; Kilby, J.R.; Grimsley, D.L.; Martinez, J.E.; Basara, J.B.; Illston, B.G.; Morris, D.A.; Kloesel, K.A. Statewide monitoring of the mesoscale environment: A technical update on the Oklahoma Mesonet. J. Atmos. Ocean. Technol. 2007, 24, 301–321. [Google Scholar] [CrossRef]
  32. Van Lier, Q.d.J. Field capacity, a valid upper limit of crop available water? Agric. Water Manag. 2017, 193, 214–220. [Google Scholar] [CrossRef]
  33. Tolk, J.A. Soils, Permanent Wilting Points. In Encyclopedia of Water Science; Taylor & Francis: New York, NY, USA, 2003. [Google Scholar]
  34. Schaap, M.G.; Leij, F.J.; Van Genuchten, M.T. Rosetta: A computer program for estimating soil hydraulic parameters with hierarchical pedotransfer functions. J. Hydrol. 2001, 251, 163–176. [Google Scholar] [CrossRef]
  35. Dane, J.H.; Hopmans, J.W. Water retention and storage—Introduction. In Methods of Soil Analysis Part 4; Dane, J.H., Topp, G.C., Eds.; Soil Science Society of America: Madison, WI, USA, 2002; pp. 671–717. [Google Scholar]
  36. Hunt, E.D.; Hubbard, K.G.; Wilhite, D.A.; Arkebauer, T.J.; Dutcher, A.L. The development and evaluation of a soil moisture index. Int. J. Climatol. 2009, 29, 747–759. [Google Scholar] [CrossRef] [Green Version]
  37. NRCS. Web Soil Survey. 2009. Available online: http://www.websoilsurvey.ncsc.usda.gov/app/ (accessed on 29 October 2017).
  38. Cassel, D.; Nielsen, D. Field Capacity and Available Water Capacity. In Methods of Soil Analysis Part 1; Klute, A., Ed.; Soil Science Society of America: Madison, WI, USA, 1986; pp. 901–926. [Google Scholar]
  39. Moriasi, D.N.; Arnold, J.G.; Van Liew, M.W.; Bingner, R.L.; Harmel, R.D.; Veith, T.L. Model evaluation guidelines for systematic quantification of accuracy in watershed simulations. Trans. ASABE 2007, 50, 885–900. [Google Scholar] [CrossRef]
  40. Addiscott, T.; Whitmore, A. Computer simulation of changes in soil mineral nitrogen and crop nitrogen during autumn, winter and spring. J. Agric. Sci. 1987, 109, 141–157. [Google Scholar] [CrossRef]
  41. Willmott, C.J. On the validation of models. Phys. Geogr. 1981, 2, 184–194. [Google Scholar] [CrossRef]
  42. Mishra, A.; Vu, T.; Veettil, A.V.; Entekhabi, D. Drought monitoring with soil moisture active passive (SMAP) measurements. J. Hydrol. 2017, 552, 620–632. [Google Scholar] [CrossRef]
  43. Cosh, M.H.; Jackson, T.J.; Bindlish, R.; Prueger, J.H. Watershed scale temporal and spatial stability of soil moisture and its role in validating satellite estimates. Remote Sens. Environ. 2004, 92, 427–435. [Google Scholar] [CrossRef]
  44. Montgomery, D.C.; Peck, E.A.; Vining, G.G. Introduction to Linear Regression Analysis; John Wiley & Sons: Hoboken, NJ, USA, 2012; Volume 821. [Google Scholar]
  45. Singh, J.; Lo, T.; Rudnick, D.R.; Dorr, T.J.; Burr, C.A.; Werle, R.; Shaver, T.M.; Muñoz-Arriola, F. Performance assessment of factory and field calibrations for electromagnetic sensors in a loam soil. Agric. Water Manag. 2018, 196, 87–98. [Google Scholar] [CrossRef]
  46. Adeyemi, O.; Norton, T.; Grove, I.; Peets, S. Performance evaluation of three newly developed soil moisture sensors. In Proceedings of the CIGR-AgEng Conference, Aarhus, Denmark, 26–29 June 2016; pp. 26–29. [Google Scholar]
  47. Kisekka, I.; Aguilar, J.; Lamm, F.; Rogers, D. Using Soil Water and Canopy Temperature to Improve Irrigation Scheduling for Corn. In Proceedings of the 2014 Irrigation Association Conference, Phoenix, AZ, USA, 3–6 November 2014. [Google Scholar]
  48. Michel, A.; Brown, H.; Gillespie, R.; George, M.; Meenken, E. Automated measurement of crop water balances under a mobile rain-exclusion facility. Agron. N. Z. 2015, 45, 39–46. [Google Scholar]
  49. Wyatt, B.M.; Ochsner, T.E.; Fiebrich, C.A.; Neel, C.R.; Wallace, D.S. Useful Drainage Estimates Obtained from a Large-Scale Soil Moisture Monitoring Network by Applying the Unit-Gradient Assumption. Vadose Zone J. 2017, 16, 16. [Google Scholar] [CrossRef]
  50. Sugita, M.; Kubota, A.; Higuchi, M.; Matsuno, A.; Tanaka, H. Continuous soil moisture monitoring under high salinity conditions by dielectric sensors: A reliability test. Tsukuba Geoenviron. Sci. 2016, 12, 17–22. [Google Scholar]
Figure 1. Experimental study site locations.
Figure 1. Experimental study site locations.
Sensors 18 03786 g001
Figure 2. Time series of sensor-estimated θv along with point measurements of θref at (a) lower salinity and lower clay content (LSLC) and (b) higher salinity and higher clay content (HSHC) sites. Error bars for θref represent standard error of mean. The FC and WP limits were determined in the laboratory.
Figure 2. Time series of sensor-estimated θv along with point measurements of θref at (a) lower salinity and lower clay content (LSLC) and (b) higher salinity and higher clay content (HSHC) sites. Error bars for θref represent standard error of mean. The FC and WP limits were determined in the laboratory.
Sensors 18 03786 g002
Figure 3. Sensor-estimated θv vs. θref at LSLC and HSHC sites.
Figure 3. Sensor-estimated θv vs. θref at LSLC and HSHC sites.
Sensors 18 03786 g003
Figure 4. Time series of sensor-estimated Bulk electricity conductivity (EC) at (a) LSLC and (b) HSHC sites.
Figure 4. Time series of sensor-estimated Bulk electricity conductivity (EC) at (a) LSLC and (b) HSHC sites.
Sensors 18 03786 g004
Figure 5. Time series of hourly soil moisture depletion (SMD) estimated based on sensor readings of θv and FC estimates from laboratory (a) and ranking (b) methods at LSLC site and laboratory (c) and ranking (d) methods at HSHC site. Dots represent SMD estimated based on θref and FC estimates from laboratory method.
Figure 5. Time series of hourly soil moisture depletion (SMD) estimated based on sensor readings of θv and FC estimates from laboratory (a) and ranking (b) methods at LSLC site and laboratory (c) and ranking (d) methods at HSHC site. Dots represent SMD estimated based on θref and FC estimates from laboratory method.
Sensors 18 03786 g005
Table 1. Soil properties at study sites.
Table 1. Soil properties at study sites.
SiteSoil TextureParticle Size DistributionEC ¥θv (m3 m−3)Ksat
% Sand% Silt% Clay(dS m−1)Sat. FC §WP *(mm day−1)
LSLCFine sandy loam72.214.413.41.20.340.170.05390.0
HSHCSilty clay loam23.537.838.77.00.390.320.2132.4
¥ Electrical conductivity. Saturation level; § Field capacity at −33 kPa; * Wilting point at −1500 kPa; Saturated hydraulic conductivity.
Table 2. Twenty-year (1997–2016) average annual and study period (July to October 2017) meteorological parameters obtained from Oklahoma Mesonet weather network.
Table 2. Twenty-year (1997–2016) average annual and study period (July to October 2017) meteorological parameters obtained from Oklahoma Mesonet weather network.
ParameterAnnualStudy Period
LSLCHSHCLSLCHSHC
Total Prec. 1 (mm)752616451340
Mean Rs 2 (MJ m−2)17.117.719.921.8
Minimum Tair 3 (°C)9.410.018.118.9
Maximum Tair (°C)22.124.130.631.2
Mean Tair (°C)15.416.823.924.8
Minimum RH 4 (%)41.937.645.544.7
Mean VPD 5 (kPa)0.91.01.01.1
Mean U2 6 (m s−1)2.52.53.02.5
1 Precipitation; 2 Daily accumulation of solar radiation; 3 Daily air temperature; 4 Daily relative humidity. 5 Daily vapor pressure deficit; 6 Daily wind speed at 2.0 m above the ground.
Table 3. Performance indicators of soil moisture sensors.
Table 3. Performance indicators of soil moisture sensors.
IndicatorsTDR315CS655GS1SM100CropX
LSLCHSHCLSLCHSHCLSLCHSHCLSLCHSHCLSLCHSHC
RMSE (m3 m−3)0.0280.0640.0190.1650.0480.1220.1100.2330.0510.055
RSR0.761.550.533.991.312.973.005.662.531.34
MBE (m3 m−3)0.0200.0530.0080.1600.0420.1210.1080.2330.045−0.049
k0.850.690.940.300.690.410.440.260.580.75
Table 4. Parameters and the p-values of the linear regression equation: θref = Slope × (sensor θv) + Intercept.
Table 4. Parameters and the p-values of the linear regression equation: θref = Slope × (sensor θv) + Intercept.
SiteSensorInterceptSloper2p-Value
LSLCTDR315−0.0170.9750.800.001
CS6550.0360.7710.85<0.001
GS10.0170.7370.700.005
SM100−0.0330.7470.840.001
CropX−0.0521.0300.570.018
HSHCTDR3150.0560.6830.850.001
CS655−0.0560.7740.200.267 Ŧ
GS1−0.1080.9710.730.007
SM100−0.1650.8730.790.003
CropX0.1370.6560.850.001
Ŧ The linear regression model was not statistically significant.
Table 5. Pearson correlation coefficients among installed sensors at study sites.
Table 5. Pearson correlation coefficients among installed sensors at study sites.
LSLC
TDR315CS655GS1SM100CropX
TDR3151.00
CS6550.991.00
GS10.970.991.00
SM1000.950.950.921.00
CropX0.790.810.810.791.00
HSHC
TDR315CS655GS1SM100CropX
TDR3151.00
CS6550.501.00
GS10.970.571.00
SM1000.900.480.901.00
CropX0.860.420.850.781.00
Note all correlation coefficients were significant at p = 0.05.
Table 6. Estimates of field capacity (FC), wilting point (WP), and available water content (AWC) (all in m3 m−3) obtained from various methods.
Table 6. Estimates of field capacity (FC), wilting point (WP), and available water content (AWC) (all in m3 m−3) obtained from various methods.
MethodLSLCHSHC
FCWPAWCFCWPAWC
Laboratory 10.170.060.110.320.210.09
Rank-TDR315 20.270.160.110.490.290.20
Rank-CS655 20.270.120.150.510.430.08
Rank-GS1 20.310.160.150.500.370.13
Rank-SM100 20.370.230.140.620.480.14
Rank-CropX 20.280.170.110.340.180.16
Rosetta-TC 30.170.060.110.310.120.19
Rosetta-TI 40.170.070.100.290.140.15
Rosetta-TBD 50.150.070.080.260.140.12
USDA-WSS 60.210.120.090.290.210.08
1 Laboratory measurement; 2 Ranking method performed for each sensor; 3 Rosetta model using soil textural class only; 4 Rosetta model using soil textural information (% sand, silt, and clay); 5 Rosetta model using textural information and bulk density; 6 USDA’s Web Soil Survey.

Share and Cite

MDPI and ACS Style

Datta, S.; Taghvaeian, S.; Ochsner, T.E.; Moriasi, D.; Gowda, P.; Steiner, J.L. Performance Assessment of Five Different Soil Moisture Sensors under Irrigated Field Conditions in Oklahoma. Sensors 2018, 18, 3786. https://0-doi-org.brum.beds.ac.uk/10.3390/s18113786

AMA Style

Datta S, Taghvaeian S, Ochsner TE, Moriasi D, Gowda P, Steiner JL. Performance Assessment of Five Different Soil Moisture Sensors under Irrigated Field Conditions in Oklahoma. Sensors. 2018; 18(11):3786. https://0-doi-org.brum.beds.ac.uk/10.3390/s18113786

Chicago/Turabian Style

Datta, Sumon, Saleh Taghvaeian, Tyson E. Ochsner, Daniel Moriasi, Prasanna Gowda, and Jean L. Steiner. 2018. "Performance Assessment of Five Different Soil Moisture Sensors under Irrigated Field Conditions in Oklahoma" Sensors 18, no. 11: 3786. https://0-doi-org.brum.beds.ac.uk/10.3390/s18113786

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop