Original Research ARTICLE
A New Optical Remote Sensing Technique for High-Resolution Mapping of Soil Moisture
- 1Department of Environmental Science, The University of Arizona, Tucson, AZ, United States
- 2Department of Electrical and Computer Engineering, Purdue University Northwest, Hammond, IN, United States
- 3School of Plant Sciences, The University of Arizona, Tucson, AZ, United States
- 4Department of Earth and Atmospheric Science, Saint Louis University, St. Louis, MO, United States
- 5Department of Civil, Environmental and Geo-Engineering, University of Minnesota, Minneapolis, MN, United States
- 6Arizona Experiment Station, The University of Arizona, Tucson, AZ, United States
- 7Department of Plants, Soils and Climate, Utah State University, Logan, UT, United States
The recently developed OPtical TRApezoid Model (OPTRAM) has been successfully applied for watershed scale soil moisture (SM) estimation based on remotely sensed shortwave infrared (SWIR) transformed reflectance (TRSWIR) and the normalized difference vegetation index (NDVI). This study is aimed at the evaluation of OPTRAM for field scale precision agriculture applications using ultrahigh spatial resolution optical observations obtained with one of the world's largest field robotic phenotyping scanners located in Maricopa, Arizona. We replaced NDVI with the soil adjusted vegetation index (SAVI), which has been shown to be more accurate for cropped agricultural fields that transition from bare soil to dense vegetation cover. The OPTRAM was parameterized based on the trapezoidal geometry of the pixel distribution within the TRSWIR-SAVI space, from which wet- and dry-edge parameters were determined. The accuracy of the resultant SM estimates is evaluated based on a comparison with ground reference measurements obtained with Time Domain Reflectometry (TDR) sensors deployed to monitor surface, near-surface and root zone SM. The obtained results indicate an SM estimation error between 0.045 and 0.057 cm3 cm−3 for the near-surface and root zone, respectively. The high resolution SM maps clearly capture the spatial SM variability at the sensor locations. These findings and the presented framework can be applied in conjunction with Unmanned Aerial System (UAS) observations to assist with farm scale precision irrigation management to improve water use efficiency of cropping systems and conserve water in water-limited regions of the world.
Introduction and Motivation
Soil moisture is a key hydrologic state variable that links land surface and atmospheric processes (Babaeian et al., 2019). Detailed knowledge about the state of SM and its spatial and temporal dynamics is of crucial importance for crop production to avoid water stress, but also for the mitigation of adverse environmental impacts due to over-irrigation as well as for the conservation of water resources, which is a pressing issue in water-limited regions of the world (Tuller et al., 2019). For example, the water scarcity projections are dire for much of the rapidly growing southwestern United States, for which the Colorado River is the major water source for over 36 million people and for close to 6 million acres of irrigated farmland (Owen, 2017). This escalating water crises demands development and adoption of transformative technologies for precision irrigation management.
Today, airborne remote sensing (RS) techniques with Unmanned Aerial Systems (UAS) provide an exceedingly powerful means for high temporal and spatial resolution SM observations (Stark et al., 2015). However, most of the currently employed RS techniques for SM estimation have been developed and evaluated for coarse-resolution (several tens of kilometers) satellite observations in the optical (Sadeghi et al., 2017), thermal (Shafian and Maas, 2015), and microwave (Kerr et al., 2001; Entekhabi et al., 2010) electromagnetic domains. Despite significant advances in large-scale SM estimation, for example, methods based on microwave observations (e.g., SMAP, SMOS) have limited applicability at smaller scales such as cropped fields, which renders them unsuitable for farm-level precision irrigation management. This advocates the need for the development and implementation of high-resolution RS techniques amenable for field scale SM monitoring and mapping.
Sadeghi et al. (2017) proposed the physically-based OPtical TRApezoid Model (OPTRAM) for estimation of spatiotemporal surface soil moisture dynamics based on the pixel distribution within the Normalized Difference Vegetation Index (NDVI) and shortwave infrared transformed reflectance (TRSWIR) space. More specifically, the OPTRAM estimates surface soil moisture based on the physical relationship between TRSWIR and the vegetation cover (Effati et al., 2019). The OPTRAM has been successfully validated with satellite data (i.e., Sentinel-2, Landsat-8, and MODIS) for several watersheds in the U.S. with vastly different climatic conditions, surface topologies, and vegetation covers via comparison of the model's surface soil moisture estimates with ground reference measurements (Sadeghi et al., 2017; Babaeian et al., 2018). Mananze et al. (2019) applied the OPTRAM for agricultural drought monitoring in Mozambique. The advantages of the OPTRAM are twofold—it does not require thermal data such as traditional triangle or trapezoid models, and it can be universally parameterized for a given location because the TRSWIR-soil moisture relationship is not affected by ambient environmental factors (e.g., air temperature and wind speed).
Motivated by the successful application of the OPTRAM at the watershed scale, we hypothesize that the model is also applicable to high spatial resolution visible near infrared (Vis-NIR) and shortwave infrared (SWIR) observations as the boundaries of the optical trapezoid are expected to be more distinct for the rather homogeneous surface conditions of agricultural fields when compared to the heterogeneous land surfaces of watersheds. Here, we test this hypothesis with ultrahigh spatial resolution Vis-NIR, and SWIR observations obtained with the TERRAPHENOTYPING–REFERENCE (TERRA-REF) platform, one of the world's largest field phenotyping robots located in Maricopa (AZ), and postulate that if the OPTRAM is applicable for SM mapping with ultrahigh spatial resolution optical reflectance data, it is also amenable for parametrization with high-resolution data captured with UAS's, and thus can be applied for farm scale precision irrigation management.
Optical Trapezoid Model
The volumetric SM content (cm3 cm−3) is obtained from the trapezoidal geometry of the transformed SWIR reflectance (TRSWIR)—vegetation index feature space (Sadeghi et al., 2017) as conceptually shown in Figure 1. For this study, we replaced the previously applied NDVI with the SAVI because it is generally more sensitive to sparse vegetation cover (i.e., the initial and development stages of agricultural crops) as it corrects for the influence of soil brightness when vegetation density is low. This leads to independence of the OPTRAM parameters from soil type. The volumetric SM is calculated as:
Where i and s are the intercept and slope of the dry- and wet-edge, respectively (Figure 1), n is the soil porosity, L is the soil brightness correction factor, and RRed, RNIR, and RSWIR denote the reflectances within the Red (636–673 nm), NIR (760–823 nm), and SWIR (2,110–2,290 nm) spectral bands captured with the TERRA-REF Vis-NIR and SWIR imaging systems. The average spectral frequencies of the applicable bands were used to derive RRed, RNIR, and RSWIR. The value of L varies with the extent of the green vegetation cover. For dense vegetation L = 0, and for bare soil L = 1. Generally, L = 0.5 is applied as the default value. When L = 0, then the SAVI equals the NDVI.
Ultrahigh Resolution Optical Observations and Reference Moisture Measurements
The optical observations (i.e., Vis-NIR, and SWIR) for calculation of the SAVI, TRSWIR, and SM (Equations 1–3) were acquired at ultrahigh spatial resolution with the TERRA-REF platform (Figure 2a). This state-of-the-art instrument located in Maricopa, Arizona consists of a 30-ton steel gantry with suspended instrument box that contains various imaging systems and sensors, including hyperspectral Vis-NIR and SWIR cameras (Figure 2b) with their specifications listed in Table 1. The robot autonomously moves along two 200-meter steel rails and continuously images the soil and crop below, providing reflectance measurements while the land surface transitions from bare soil to full plant canopy (Burnette et al., 2018).
Figure 2. Aerial view of the TERRA-REF field with marked locations (black squares) of the TDR sensor nests and the East-West swats (yellow rectangles) captured by the Vis-NIR and SWIR cameras (a), the TERRA-REF steel gantry with suspended instrument box (b), close up of the instrument box with various imaging systems and sensors (c), and the True TDR-315 sensors installed in duplicate in 2, 10, and 50 cm depths at the three sensor nest locations (d).
The hyperspectral imagers (Headwall Photonics, Bolton, MA, USA) are push-broom sensors that acquire 939 (Vis-NIR) and 275 (SWIR) spectral bands within the 380–1,000 nm and 900–2,500 nm electromagnetic frequency ranges, respectively (Figure 2c). The distance between the soil surface and the hyperspectral cameras was 2-m. The Vis-NIR hyperspectral images for calculation of the SAVI were radiometrically calibrated using real-time down-welling irradiance captured with a spectral irradiance meter (Ocean Optics, Largo, FL, USA) mounted on top of the gantry, which covers the range from 337 to 823 nm. The SWIR images were calibrated with white and gray reference panels with known spectral reflectance values placed at the 2-m distance below the imaging system (Figure 3). Prior to irradiance-based calibration, the dark current was subtracted from the raw image to correct for detector artifacts. Hyperspectral bands over each domain were averaged and used for analysis.
Figure 3. Reference panels exhibiting various spectral reflectances used for radiometric calibration of Vis-NIR and SWIR observations.
The OPTRAM-estimated SM values were then evaluated based on soil moisture reference measurements with state-of-the-art time domain reflectometry (TDR) sensors (True TDR-315, Acclima, Inc., Meridian, ID, USA) installed at three locations within the scanner field at 2, 10, and 50 cm depths (Figure 2d). The TDR-315 sensor houses the entire measurement circuitry, including a microprocessor, within the sensor head and communicates with a datalogger for transfer of processed data via the SDI-12 protocol. The TDR-measured soil moisture data were recorded at 15-min intervals with CR1000 dataloggers (Campbell Scientific, Inc., Logan, UT, USA). The data used for this study are from an energy sorghum phenotyping experiment conducted from September 1 to October 28, 2018.
Results and Discussion
TRSWIR-SAVI Feature Space for Soil Moisture Estimation With OPTRAM
The TRSWIR-SAVI trapezoidal space was generated via integration of five scenes collected at the initial stage (i.e., bare soil), during sorghum development, and at the final growth stage (i.e., before harvest) from the beginning of September to the end of October, 2018. The advantages of the integrated trapezoid are increased computational efficiency and independence of dry- and wet-edge parameters from the growth stage. The dry- and wet-edges and their associated parameters (Figures 1 and 4) were determined via manual fitting. This is justifiable because of the very distinct geometry exhibited by the measured TRSWIR-SAVI relationship during the sorghum experiment (Figure 4). Furthermore, Babaeian et al. (2018) demonstrated with a thorough sensitivity analysis that the selection of the dry- and wet-edge locations does not need to be overly precise for the OPTRAM to accurately estimate surface soil moisture. It should be noted that the more homogeneous soil and vegetation cover conditions in agricultural fields in conjunction with ultrahigh spatial resolution observations alleviate problems with OPTRAM parameterization occasionally experienced for large scale, coarse-resolution satellite observations, such as oversaturation of the wet-edge or uncertainty associated with image pixels that neither belong to the soil nor to vegetation (e.g., surface water bodies, roads, or buildings) (Sadeghi et al., 2017).
Figure 4. Pixel distribution within the integrated TRSWIR-SAVI trapezoidal space for the 2018 sorghum experiment. The blue and red solid lines represent the manually fitted wet- and dry-edges, respectively. The black, green, red, yellow, and blue point clouds correspond to the observations from Sept. 15, Sept. 28, Oct. 9, Oct. 18, and Oct. 28, respectively.
Soil Moisture Estimation With OPTRAM
To evaluate the accuracy of the OPTRAM SM estimates, they were compared with the TDR reference measurements at the three sensor nest locations. The SM contents measured with the TDR sensors throughout the growth period ranged from 0.029 to 0.465 cm3 cm−3. The OPTRAM SM was obtained by averaging the pixels directly covering the TDR sensor locations (i.e., an area of 30 times 50 cm). It is interesting to note that because of the ultrahigh resolution, the data volume generated for such a small area with the TERRA-REF Vis-NIR and SWIR imaging systems exceeds several gigabytes for a single scene. Figure 5 depicts the correlation between OPTRAM-determined SM and the TDR reference measurements at 2, 10, and 50 cm depths, with correlation coefficients (R) ranging from 0.66 to 0.83 and root mean square errors (RMSE) from 0.045 to 0.057 cm3 cm−3, which is considered reasonably accurate for remote sensing of SM (Entekhabi et al., 2014). Because of the limited penetration depth of Vis-NIR and SWIR electromagnetic radiation it was expected that the OPTRAM SM estimates better match the near-surface (2-cm) TDR reference measurements. However, there is also an obvious correlation between the OPTRAM estimates and root-zone SM (Figure 5), which contradicts the common presumption that surface soil moisture is decoupled from the root-zone (Qiu et al., 2014, 2016; Tayfur et al., 2014). Our results support findings in Sadeghi et al. (2019) and Koster et al. (2018), who established links between near-surface SM and the SM in deeper depths, and are in line with results shown in Wu and Dickinson (2004) and Ford et al. (2014), who demonstrate that root-zone SM can be accurately inferred from SWIR remote sensing observations. It is interesting to note that the OPTRAM SM estimates can be applied to directly calculate crop water consumption with a new analytical approach based on inversion of the linearized Richard's equation (Sadeghi et al., 2019). Furthermore, OPTRAM SM time series could be coupled with numerical water flow models to simulate root-zone moisture dynamics considering plant water uptake.
Figure 5. OPTRAM SM estimates compared with the reference TDR measurements at the three sensor nest locations.
Figure 6 depicts the OPTRAM soil moisture distributions across sensor nest location 1 at ultrahigh spatial resolution (i.e., 2-mm pixel size) for the early stage of sorghum development (September 28), close to maturity (October 9), and short before harvest (October 18). The maps clearly indicate the OPTRAM's capability to capture soil moisture variations when parameterized with ultrahigh resolution Vis-NIR and SWIR observations, not only for single scenes, but also for various soil moisture states related to irrigation and plant development stage. This instills confidence that the proposed method is amenable for parameterization with high-resolution UAS observations, which opens new avenues for farm scale precision irrigation management to increase crop water use efficiency while alleviating the risk of environmental contamination and contributing to conservation of strained water resources in water-limited regions of the world.
Figure 6. An example for the applicability of the OPTRAM to estimate moisture variations when parameterized with ultrahigh resolution Vis-NIR and SWIR observations. Scenarios for wet and dry conditions are shown.
This contribution should be viewed as a proof of concept for the applicability of the OPtical TRApezoid Model (OPTRAM) for estimation of SM when parameterized with remotely-sensed ultrahigh resolution Vis-NIR and SWIR reflectance data. A comparison of the OPTRAM SM estimates with reference moisture data obtained with state-of-the-art TDR sensors reveals the promising potential of the model to be applied in conjunction with rapidly evolving UAS imaging capabilities. This opens new avenues for UAS-based farm-level precision irrigation management. Such technological advances are much- needed to combat the water crisis in the arid and semiarid regions of the world that is expected to further escalate in view of the rapidly growing human population and a changing global climate. The next step is to test OPTRAM with UAS observations and apply the obtained SM information for farm-level irrigation management.
Data Availability Statement
The datasets generated for this study are available from the corresponding author on request.
EB, MT, MS, and SJ conceptualized the study and contributed to data analyses and writing of the manuscript. SW, MN, JD, RW, and DL contributed to sensor installation, field data collection, and analyses. PS, MM, and VS contributed to the calibration and analyses of the TERRA-REF Vis-NIR and SWIR data.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The authors gratefully acknowledge support from the University of Arizona (UA) College of Agriculture and Life Sciences (CALS) Innovation Venture Investment Project (iVIP). The Vis-NIR and SWIR data utilized for this study are from the TERRA-REF experiment, funded by the Advanced Research Projects Agency-Energy (ARPA-E), U.S. Department of Energy, under Award Number DE-AR0000594. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof.
Babaeian, E., Sadeghi, M., Franz, T. E., Jones, S., and Tuller, M. (2018). Mapping soil moisture with the OPtical TRApezoid Model (OPTRAM) based on long-term MODIS observations. Remote Sens. Environ. 211, 425–440. doi: 10.1016/j.rse.2018.04.029
Babaeian, E., Sadeghi, M., Jones, S. B., Montzka, C., Vereecken, H., and Tuller, M. (2019). Ground, proximal and satellite remote sensing of soil moisture. Rev. Geophys. 57, 530–616. doi: 10.1029/2018RG000618
Burnette, M., Willis, C., Kooper, R., Maloney, J. D., Ward, R., Shakoor, N., et al. (2018). “TERRA-REF data processing infrastructure,” in Proceedings of the Practice and Experience on Advanced Research Computing (Pittsburgh, PA). doi: 10.1145/3219104.3219152
Effati, M., Bahrami, H.-A., Gohardoust, M. R., Babaeian, E., and Tuller, M. (2019). Application of satellite remote sensing for estimation of dust emission probability in the Urmia Lake Basin in Iran. Soil Sci. Soc. Am. J. 83, 993–1002. doi: 10.2136/sssaj2019.01.0018
Entekhabi, D., Njoku, E. G., O'Neill, P. E., Kellogg, K. H., Crow, W. T., Edelstein, W. N., et al. (2010). The soil moisture active passive (SMAP) mission. Proc. IEEE. 98, 704–716. doi: 10.1109/JPROC.2010.2043918
Entekhabi, D., Yueh, S., O'Neill, P. E., Kellogg, K. H., Allen, A., Bindlish, R., et al. (2014). SMAP Handbook-Soil Moisture Active Passive: Mapping Soil Moisture and Freeze/Thaw from Space. Pasadena, CA: Jet Propulsion Laboratory.
Kerr, Y. H., Waldteufel, P., Wigneron, J. P., Martinuzzi, J., Font, J., and Berger, M. (2001). Soil moisture retrieval from space: the Soil Moisture and Ocean Salinity (SMOS) mission. IEEE Trans. Geosci. Remote Sens. 39, 1729–1735. doi: 10.1109/36.942551
Mananze, S., Pôças, I., and Cunha, M. (2019). Agricultural drought monitoring based on soil moisture derived from the optical trapezoid model in Mozambique. J. Appl. Remote Sens. 13:024519. doi: 10.1117/1.JRS.13.024519
Qiu, J., Crow, W. T., and Nearing, G. S. (2016). The impact of vertical measurement depth on the information content of soil moisture for latent heat flux estimation. J. Hydrometeorol. 17, 2419–2430. doi: 10.1175/JHM-D-16-0044.1
Qiu, J., Crow, W. T., Nearing, G. S., Mo, X., and Liu, S. (2014). The impact of vertical measurement depth on the information content of soil moisture times series data. Geophy. Res. Lett. 41, 4997–5004. doi: 10.1002/2014GL060017
Sadeghi, M., Babaeian, E., Tuller, M., and Jones, S. B. (2017). The optical trapezoid model: a novel approach to remote sensing of soil moisture applied to Sentinel-2 and Landsat-8 observations. Remote Sens. Environ. 198, 52–68. doi: 10.1016/j.rse.2017.05.041
Sadeghi, M., Tuller, M., Warrick, W. A., Babaeian, E., Parajuli, K., Gohardoust, M. R., et al. (2019). An analytical model for estimation of land surface net water flux from near-surface soil moisture observations. J. Hydrol. 570, 26–37. doi: 10.1016/j.jhydrol.2018.12.038
Stark, B., McGee, M., and Chen, Y. Q. (2015). “Short wave infrared (SWIR) imaging systems using small Unmanned Aerial Systems (sUAS),” in Proceedings of the International Conference on Unmanned Aircraft Systems (ICUAS) (Denver, CO). doi: 10.1109/ICUAS.2015.7152328
Tayfur, G., Zucco, G., Brocca, L., and Moramarco, T. (2014). Coupling soil moisture and precipitation observations for predicting hourly runoff at small catchment scale. J. Hydrol. 510, 363–371. doi: 10.1016/j.jhydrol.2013.12.045
Tuller, M., Babaeian, E., Jones, S. B., Montzka, C., Vereecken, H., and Sadeghi, M. (2019). The paramount societal impact of soil moisture. EOS 100. Available online at: https://eos.org/editors-vox/the-paramount-societal-impact-of-soil-moisture
Keywords: remote sensing, high-resolution, soil moisture, OPTRAM, precision irrigation, TERRA-REF
Citation: Babaeian E, Sidike P, Newcomb MS, Maimaitijiang M, White SA, Demieville J, Ward RW, Sadeghi M, LeBauer DS, Jones SB, Sagan V and Tuller M (2019) A New Optical Remote Sensing Technique for High-Resolution Mapping of Soil Moisture. Front. Big Data 2:37. doi: 10.3389/fdata.2019.00037
Received: 21 June 2019; Accepted: 15 October 2019;
Published: 05 November 2019.
Edited by:Rasmus Houborg, Planet Labs Inc, United States
Reviewed by:Michael Cosh, University of San Diego, United States
YangQuan Chen, University of California, Merced, United States
Copyright © 2019 Babaeian, Sidike, Newcomb, Maimaitijiang, White, Demieville, Ward, Sadeghi, LeBauer, Jones, Sagan and Tuller. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Markus Tuller, email@example.com