Next Article in Journal
SELWAK: A Secure and Efficient Lightweight and Anonymous Authentication and Key Establishment Scheme for IoT Based Vehicular Ad hoc Networks
Previous Article in Journal
Distributed Airflow Sensing Based on High-Spatial-Resolution BOTDA and a Self-Heated High-Attenuation Fiber
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Real-Time Assessment of Mandarin Crop Water Stress Index

1
College of Water Conservancy and Civil Engineering, South China Agricultural University, No. 483, Wushan Road, Tianhe District, Guangzhou 510642, China
2
College of Engineering, National Center for International Collaboration Research on Precision Agricultural Aviation Pesticides Spraying Technology (NPAAC), South China Agricultural University, No. 483, Wushan Road, Tianhe District, Guangzhou 510642, China
3
Department of Agricultural Engineering, University of Cape Coast, Cape Coast PMB, Ghana
4
Department of Agricultural Engineering, University of Ghana, Accra P.O. Box LG 77, Ghana
5
Soil, Water and Environmental Engineering Division, Soil Research Institute of Ghana, Kumasi PMB, Ghana
6
Soil Science Division, Cocoa Research Institute of Ghana (Ghana COCOBOD), New Tafo-Akim P.O. Box 8, Ghana
*
Author to whom correspondence should be addressed.
Sensors 2022, 22(11), 4018; https://doi.org/10.3390/s22114018
Submission received: 28 January 2022 / Revised: 14 February 2022 / Accepted: 21 February 2022 / Published: 26 May 2022
(This article belongs to the Section Optical Sensors)

Abstract

:
The use of plant-based indicators and other conventional means to detect the level of water stress in crops may be challenging, due to their difficulties in automation, their arduousness, and their time-consuming nature. Non-contact and non-destructive sensing methods can be used to detect the level of water stress in plants continuously and to provide automatic sensing and controls. This research aimed at determining the viability, efficiency, and swiftness in employing the commercial Workswell WIRIS Agro R infrared camera (WWARIC) in monitoring water stress and scheduling appropriate irrigation regimes in mandarin plants. The experiment used a four-by-three randomized complete block design with 80–100% FC water treatment as full field capacity and three deficit irrigation treatments at 70–75% FC, 60–65% FC, and 50–55% FC. Air temperature, canopy temperature, and vapor pressure deficits were measured and employed to deduce the empirical crop water stress index, using the Idso approach (CWSI(Idso)) as well as baseline equations to calculate non-water stress and water stressed conditions. The relative leaf water content (RLWC) of mandarin plants was also determined for the growing season. From the experiment, CWSI(Idso) and CWSI were estimated using the Workswell Wiris Agro R infrared camera (CWSIW) and showed a high correlation (R2 = 0.75 at p < 0.05) in assessing the extent of water stress in mandarin plants. The results also showed that at an altitude of 12 m above the mandarin canopy, the WWARIC was able to identify water stress using three modes (empirical, differential, and theoretical). The WWARIC’s color map feature, presented in real time, makes the camera a suitable device, as there is no need for complex computations or expert advice before determining the extent of the stress the crops are subjected to. The results prove that this novel use of the WWARIC demonstrated sufficient precision, swiftness, and intelligibility in the real-time detection of the mandarin water stress index and, accordingly, assisted in scheduling irrigation.

1. Introduction

Detecting crop water stress on time and precisely is very vital for both small- and large-scale precision agricultural systems, in order to avoid drought-related problems such as poor plant growth and development, reduction of yield, and yield quality [1,2,3]. In recent times, water scarcity issues are becoming widespread on a global scale, due to global warming resulting from the frequent emission of greenhouse gases, deforestation, and other environmental degradation factors, leading to erratic changes in climatic conditions. Rapid urbanization, massive industrial development, high population growth, and the increasing demand for fresh water for both domestic and commercial purposes are all significant threats to irrigated agriculture [4,5]. If unchecked, this scenario will bring about the production of crops under limiting water (water-stressed) conditions, with the resulting output being detrimental [2,6,7].
The most frequently used indicators for measuring water stress in tree crops are water potentials (Ψ) and stomatal conductance (gs), but these indicators are time consuming, very arduous, cannot be easily automated, and can only offer a point measurement in a leaf, a single soil location, a branch, or a tree [8,9,10]. Their measurements may not be true representations of the actual water status of an entire orchard field without a massive financial commitment, because several monitoring sites will be required in very large orchard setups. Furthermore, the information obtained from water potential (Ψ) measurements may not be consistent, because water potential (Ψ) decreases in plant species with isohydric characteristics (for example, in mandarin plants), due to their ability to preserve an effective stomatal regulation. This phenomenon impedes the noticeable decrease in water potential (Ψ) during limiting water (drought) conditions or during periods of high evaporative demand [10,11]. Change into right sign
Infrared thermometry has been employed from space to study the variation in plant canopy temperature as the condition of plant water varies [9,12,13]. However, the unavailability of geographic and time-based resolution of satellite images has limited the wider acceptance of the technology [14,15]. The research in [15] and [16] proved that equipping unmanned aerial vehicles (UAVs) with multispectral (IR and near-IR) cameras and correlating their data with data from the crop water stress index (CWSI) to study water status in olive and mandarin orchards was very effective. In commercial orchards, less difficult and automatic remote sensing technologies are required to rapidly and accurately estimate the water status of the plants, due to the large farm sizes and the estimated yield reduction resulting from water deficit [9,17,18,19]. It is worth noting that remote sensing technologies provide better observations in large established farms and easier means for managing the exact levels of plant water stress, thereby assisting subsequent operational decisions [20,21]. Several research projects have used less expensive and non-complex means of measuring canopy temperature (from infrared thermometers) and other related indicators such as crop water stress index (CWSI), the water potential index (Ψ), and the stomatal conductance index (Ig) to monitor and study crop water level for several years [9,22,23].
The crop water stress index (CWSI) is a standardized index used to measure crop water stress and to minimize the possible impacts of environmental conditions affecting the interactions between plant canopy temperature and air temperature [12,24,25,26]. The measurement of the CWSI is established by the relationship between the actual canopy temperature (Tc) data and the theoretical canopy temperature data, obtained from two different thresholds, namely a lower limit (when plant transpiration is maximum under full irrigation when stomata are completely opened) and an upper limit (when plants are not transpiring because they have been subjected to conditions of extreme water stress leading to complete closure of stomata) [9,26,27]. Furthermore, mandarin plant’s water stress can be determined by estimating the canopy temperature [28,29,30,31] at any given time. Nonetheless, the sensitivity of mandarins to air vapor pressure is very high, thereby minimizing the rate of transpiration under high VPD conditions and making the use of only canopy temperature very challenging [9].
Citrus water requirement is difficult to predict, due to differences in water consumption at different stages of growth, seasonal change, soil types, ground covers, ages and sizes of tree, and locations of cultivated plants. Different authors have described the different volumes of water needed by citrus. Mature citrus use about 100 to 511 L of water daily in summer [32,33,34,35]. In recent times, the CWSI model has been used to identify the water status in most crops. However, the procedures employed during the calculations of CWSI values in assessing plant water stress are very arduous and tedious, and involve a good deal of technical know-how, resulting in low farmer-researcher compliance [9,27,36,37].
To minimize the problems of farmer-researcher compliance and related issues, the Workswell Wirirs Agro R infrared camera has been designed to effectively evaluate crop water stress promptly, independent of secondary weather data, complex mathematical and scientific formulae, and workable operating system software. The WWARIC is designed to monitor crop water stress status in diverse farm situations, irrespective of the prevailing environmental conditions or the type of plants cultivated. This tool is very useful in on-farm water deficit estimation, irrespective of plant type, plants’ growth stages, and farm conditions. The aim of this experiment is to study the effectiveness and feasibility of using the WWARIC to monitor water stress and to assist in the appropriate scheduling of irrigation regimes for mandarin plants.

2. Materials and Methods

2.1. Experimental Site

The field experiment was conducted on a 576 m2 (0.142 acres) mandarin orchard that was already established at the South China Agricultural University Citrus Experimental Orchard near Yangcun town, Boluo county, Huizhou city of the Guangdong Province in the Peoples’ Republic of China, which has geographical coordinates of 36°16′2″ North and 115°28′32″ East, a mean annual rainfall of 1955 mm, and a mean annual temperature of 31 °C. The area is typical of a subtropical monsoon area, with the Tropic of Cancer crossing the city. The climate is mild and humid, with frequent rains from April to August [38].
The orchard trees were planted in 2014 at a spacing of 2.5 m within a block and 3.6 m between blocks. A four-by-three randomized complete block design (RCBD) was used for this field experiment. Each plot consisted of three sub-blocks with three mandarin plants each were studied, constituting a plant population of 36 distributed over 12 experimental blocks (plots). Each treatment was replicated three times. To minimize the possible effects on the treatments (especially the purposefully stressed treatments) of rainfall and evaporation of water from the soil, the demarcated field was covered with black woven polyethene mulch to ascertain, during the study period, whether the loss of irrigation water was only through transpiration and deep percolation. Figure 1 shows the layout of the experimental plot.

2.2. Soil Data

The sampled soil collected and analyzed revealed that the type of soil at the study field was loamy sand made up of 72.1% sand, 19.4% silt, and 8.5% clay. Soil samples were collected and tested before treatments. The soil’s physical and chemical properties are shown in Table 1.

2.3. Water Stress Treatments

Irrigation scheduling was based on recharging the volume of water in the soil back to the previous soil moisture content at field capacity before usage by the plant for photosynthesis and transpiration. Four different irrigation treatments, i.e., full irrigation (80–100% FC) and regulated deficit irrigation schedules (70–75% FC, 60–65% FC and 50–55% FC) were applied throughout the study period.
Water was supplied by means of a gravity-driven drip system with an adjusted flow rate of 6 L/h. Prior to each irrigation scheduling, the soil gravimetric water content was measured and expressed in terms of volumetric water content (θv). These data were compared with the soil water content (θv) data obtained by measuring with an IMKO TRIME-TDR probe that was inserted horizontally along the walls of ditches dug closed to the plant root zone (about 0.5 m from the base of the plant) at different depth intervals. The amount of irrigation water applied was estimated from the equation used in [39,40]:
d n = FC SMC 100 × BD × D
where FC = the field capacity of the soil, BD = the bulk density of the soil, D = the rootzone depth, and SMC = the prevailing soil moisture content on a volume basis, dn = amount of irrigation water applied.

2.4. Fertilizer Application, Pest and Disease Management

In May 2021, the soil nutrient was improved by adding 4.5 kg h−1 of N-P2O5-K2O fertilizer (20-10-10 of nitrogen, phosphorus, and potassium, respectively) as a top-dressing 30 cm from the plant stem into the soil, prior to treatment allocation. The farm was often sprayed with a mixture of Fenproximate, Azadirachtin, Bacillus thuringiensis, and other organic pesticides and insecticides to minimize the destruction caused by pests and insects, through a rotary center-pivot overhead mist sprayers mounted on metal pipes about 5 m high at 5 m intervals across the length and breadth of the field.

2.5. CWSI Workswell Wiris Agro R Infrared Camera

2.5.1. Camera Specification, Data Collection, and Processing

The WWARIC is designed to be mounted on UAVs for assessing the crop water stress index (CWSI) as well as the temperature. The camera, endowed with a powerful operating system (WIRIS OS), ensures full accessibility and utilization of all the camera’s operational functions during instantaneous data streaming and assessment, for onboard operations or offboard use. The WWARIC can be operated simply through the MavLink protocol, the RJ-45 connector, the S. Bus system, the CAN bus system, or a triggers support system. With a “plant biomass coverage” index, this camera displays the percentage of a real-time calculated mass of vegetation in the RGB setup. The CWSI, estimated with the WWARIC, is normalized within the limit of 0 to 1 (or 0–100%, with 100% being very stressed conditions while movement towards 0 represents less stressed conditions), representing the extent of stress on the crop in relation to pixel value. This information is vital in developing yield maps and making precise water management decisions in large-scale agriculture environments.
The Workswell Wiris Agro R infrared camera has a long-wave infrared band sensor with a specific sensor resolution of 640 by 512 pixels and a focal-plane array (FPA) active sensor with a size of 1088 × 10−5 by 8705 × 10−6 m. The field of view of the WWARIC lens is 45° and the standard temperature responsiveness is 0.05 °C (50 mK). The camera comes with four color maps: crop pallets, cropstep pallets, water pallets, and waterstep pallets that are used to manage and assess the CWSI. The WWARIC can perform digital zooming of images up to 14 times continuously (i.e., 1–14×). The camera uses the advanced Workswell CorePlayer application software package for offline CWSI data analysis, which comes with two licenses. The 3D mapping software works perfectly with Agisoft and Pix4D softwares. Operationally, the 10× optical zoom RGB has a resolution of 1920 × 1080 pixels (full high-definition), a 1 3 inch sensor, automatic white offset, a powerful varied range, backlight counterbalance, exposure, and gamma control systems. The WWARIC has 3.0 to 350 mm focal length and 10× optical zoom that supports vibration and a view angle as low as 6.9° (for ultra zoom) with a maximum of 58.2° (for extra-wide viewing). The WWARIC has an auto-focus function that synchronizes a direct focus zoom, an in-built high-speed solid-state drive with 128 GB capacity for storing captured images and videos, an exterior secure digital card slot, and a USB 2.0 port.
The WWARIC stores images and videos in the following layouts: JPEG images, TIFF images, digital camera full high-definition JPEG images, and digital camera h.264 encode video high-definition recordings and full-frame video recordings (raw data recordings). With an average mass of 4.30 kg and dimensions of 8.30 cm long, 8.50 cm wide, and 6.80 cm high, the camera can be fixed on a 2 × 1 4 − 20 UNC support system connected to a tripod stand, an unmanned aerial vehicle, or a satellite system. The camera is encased in a robust aluminum plate that can endure external temperature fluctuations, and which serves as a protector for the internal parts against being damaged, to ensure long-term measurement stability. The operating and storage temperatures are, respectively, a minimum of −10 °C and a maximum of 50 °C and a minimum of 30 °C and a maximum of 60 °C. The camera has an input power source of 9–36 V DC, with a coaxial cable of 2 × 6.4 mm and an outer shell GND that dissipates 12 W of power.

2.5.2. Camera Measurement Functions

The Workswell Wiris Agro R infrared camera has an inbuilt four-color map with full radiometric (temperature) information. The camera can be operated rapidly to capture videos and images simultaneously as thermal imagery, digital imagery, or digital video imagery. During operation, the camera images can be viewed as picture-in-picture (PIP), full screen RGB with segmentation, or dual-screen. There are three micro-HDMI video outputs in the formats: 1280 × 720 pixels (720p), an aspect ratio of 16:9, and a micro-HDMI video output.

2.6. Field Measurements

2.6.1. Image Acquisition

The WWARIC, purchased from Albrechtova Vrchu, Prague, Czech Republic, was mounted on a metallic pole, 12 m high, to capture images of mandarins cultivated under different soil water contents at various field capacities to estimate the water levels in the mandarin plants. The CWI-640-545-128 camera specifications are as follows: it was released on 20 January 2021 with revision number 210120, serial number 1010-CWI-201002, version 1.3.8, and resolution 640 × 512. The RTC temperature is 37.6 °C, with CPU and IR core temperatures of 46.3 °C and 45.8 °C, respectively. Data gained by the WWARIC were processed with CorePlayer software version 1.4.1. The CWSIW was compared with the CWSI(Idso) and the relative leaf water content (RLWC) to determine precision through correlation.

2.6.2. Canopy Temperature Measurement

The mandarin canopy temperature was calculated with the help of a Raytek® portable, noncontact, infrared thermometer having an adjusted emissivity of 0.95 Wm−2 with a 0.2 °C (0.5 °F) display resolution. The measurement was done by holding the IR thermometer 0.30 m directly above the mandarin leaf, with the optical laser spot directed vertically to the targeted leaf [20,41]. The temperature readings were recorded once every 7 days after July 2021, when the first irrigation dose was administered until harvesting. The mandarin leaf temperatures of three sampled plants per treatment were recorded from four directions, with north as the initial point of reference, and the average temperatures of each treatment were estimated. The temperature recordings were usually measured on very clear sky sunny days between 11:00 A.M. and 14:00 P.M., on fully developed sunlit mandarin leaves at the tip of the branches, as proposed by [42,43,44].

2.6.3. Estimating the Relative Leaf Water Content (RLWC)

The RLWC was estimated every week by weighing freshly plucked sampled mandarin leaves to deduce their fresh weight (FW), subsequently soaking and bottling them in distilled water at 25 °C for 144 h to ensure full turgidity, and then weighing them again to determine the weight when fully turgid (TW). The turgid leaves were cleaned, oven-dried at 80 °C for 48 h and reweighed to deduce the dried weight (DW). The relative leaf water content (RLWC) was then estimated according to the equation used by [20,45]:
RLWC = FW DW TW DW × 100
where RLWC (%) = the relative leaf water content, FW = the weight of fresh leaves, DW = the weight of dried leaves, and TW = the weight of turgid leaves.

2.6.4. Estimating SPAD and the Nitrogen Sufficiency Index (NSI)

The chlorophyll content of each leaf was measured with a portable Konica Minolta SPAD 502 chlorophyll meter. The SPAD measurements were made on new sampled and completely opened leaves every week, from the beginning of irrigation. The measured SPAD values were converted into the nitrogen sufficiency index (NSI) by finding the ratios of the measured leaf chlorophyll concentration in each plant for each treatment to the measured chlorophyll concentration of a non-stressed reference plant, on the same day of measurement as described by [46] in the following equation:
NSI = SPAD M SPAD R
where NSI is the nitrogen sufficiency index, SPADM is the measured chlorophyll content in the observed plant, and SPADR is the measured chlorophyll content in the fully irrigated reference mandarin crop.

2.6.5. Estimating the CWSI Using an Empirical Method (the Idso Method)

The CWSI can be determined experimentally by studying the differences in plant canopy temperature and air temperature (Tc  Ta), as well as considering the vapor pressure deficit (VPD) of the air on a well-irrigated plant on a clear, sunlit, unclouded day with maximum plant environmental settings [27,36,37]. In calculating the CWSI, two focal points were designated and used as a lower-referenced point (the non-water stressed baseline dependent on VPD) and an upper-referenced point (the non-transpiring canopy/water stressed baseline not dependent on VPD) [27,36,37,47,48,49,50].
The CWSI empirical value is estimated using the equation applied in [36]:
CWSI Idso = ( T C T a ) D 2   D 1 D 2
where Tc is the mandarin leaf cover temperature (°C); Ta is the orchard air temperature (°C); D2 = the lower reference point ( ( T C T a ) × VPD for well-watered mandarin plants, i.e., non-stressed conditions obtained from 80–100 FC treatment); and D1 = the upper reference point (the difference in canopy and air temperature difference for highly watered mandarin plants (extremely water stressed). The upper reference point was determined from canopy temperature data obtained from fully stressed mandarin plants (50–55% FC treatment) during the peak of the summer period when the maximum transpiration was envisaged [47,48]. The weather data used to calculate the CWSI by Idso procedures are presented in Table 2.
This CWSI empirical method, or CWSI(Idso), expresses the plant water level within the range of 0 to 1, indicating no or less water stress levels and high-water stress status, respectively. However, negative values may also be obtained, indicating the presence of excess water.

2.6.6. Estimating the CWSI with the Workswell WIRIS Agro R Infrared Camera

The WWARIC is designed to estimate the crop water stress index via three different methods. The CWSI model differential (CWSId) mode represents the fundamental method that employs the measurement of air temperature (Tair) as the only parameter during estimation. The CWSId value is estimated from the differences in the crop canopy temperature and the prevailing air temperature. The second method, the CSWI model empirical (CWSIe) mode, is based on measuring the temperatures of two reference points, Twet (representing no stress or a 0% stress level) and Tdry (representing fully stressed or a 100% stress level). The CWSI model theoretical (CWSIt) mode is the third model that estimates the CWSI using the following parameters: the category of plant, the canopy temperature (°C), the air temperature (°C), the relative humidity (%), the baseline slopes, and the baseline intercept.

2.7. Data Analysis

Statistical analyses were conducted using the IBM Statistical Package for Social Sciences version 21 software. One-way ANOVA was employed to analyze the differences in the water stress parameters studied (the leaf temperature (LT°C), the relative leaf water content (RLWC) (%), the leaf chlorophyll content (LCC) (mg g−1), and the water stress index, CWSI(Idso)) within treatments and blocks. The treatment means were separated by the Tukey HSD method to indicate significant differences among the treatment means, at 0.001, 0.01, and 0.05 levels of significance. Linear regression analysis was used to show the relation between soil volumetric water content and soil water content, measured with the time domain reflectometer, as well as the water stress indices obtained by both the Idso approach and the WWARIC.

3. Results

3.1. Soil Water Content Measurement

This study showed that the estimated volumetric soil water content and the TDR readings were strongly correlated. The correlation coefficient (R2) was found to be R2 = 0.93 at p < 0.05, as shown in Figure 2. The total irrigation water supplied to the mandarins from July 2021 to 30 November 2021, was 184.64 mm, 134.02 mm, 116.40 mm, and 100.14 mm at 80–100% FC, 70–75% FC, 60–65% FC, and 50–55% FC, respectively.

Variations in Soil Moisture Content

Fluctuations in weather conditions resulted in a corresponding variation in soil moisture content throughout the entire research period; hence, there were variations in the depths of the irrigation water supplied. The study showed that the irrigation water supplied was highest for all treatments during mid-summer (from 3 August 2021 to 31 August 2021) when transpiration was very high due to high daily temperatures. In contrast, the depth of the irrigation water supplied was low in July (when irrigation was started) for all the treatments but increased after July. In autumn, there was a decrease in the depth of the irrigation water supplied, resulting from the decrease in temperatures. The irrigation amount supplied and the weekly variation in irrigation water are presented in Table 3 and Figure 3, respectively.

3.2. Plant-Based Indicators

3.2.1. Relative Leaf Water Content (RLWC)

The relative leaf water content increased with increasing volume of water available to the plant. The research showed that there was a gradual increase in RLWC across all the treatments from the beginning of irrigation treatment, notwithstanding that there were few weeks where there was little decline in the values observed. The estimated average relative leaf water content (%) during the study period was 93.93, 86.79, 84.59, and 82.88 for treatments 80–100% FC, 70–75% FC, 60–65% FC, and 50–55% FC, respectively. The observed relative leaf water content for the study period is presented in Figure 4.

3.2.2. Leaf Temperature (°C)

The plant canopy temperature (leaf temperature) varied within treatments and between treatments. The average seasonal values recorded were 34.49, 35.64, 36.36, and 37.01 for the fully irrigated plants (80–100% FC) and the deficit irrigated treatments (70–75% FC, 60–65% FC, and 50–55% FC), respectively.

3.2.3. Leaf Chlorophyll (mg g−1) Content and the Nitrogen Sufficiency Index

The weekly SPAD dataset measured during the period were highly significant. The observed SPAD values increased with the increasing soil water content at field capacity. The average SPAD values estimated were 72.3, 69.4, 68.5 and 65.4 mg g−1 for treatments 80–100% FC, 70–75% FC, 60–65% FC, and 50–55% FC, respectively. The average NSI values observed were 1, 0.948, 0.935, and 0.894 for treatments 80–100% FC, 70–75% FC, 60–65% FC, and 50–55% FC, respectively.

3.3. Mandarin Crop Water Stress Index (CWSI) Based on Idso Procedures and Baseline Equations

The CWSI values varied under the different irrigation treatments imposed on the plants from July 2021 to November 2021. The vapor deficit pressure was estimated from the relationship developed by [51] and utilized by [26,52], as follows;
VPD = 0.6108 × exp   ( 17.27   T a T a + 237.3 ) × 100 RH 100
where VPD   = vapor deficit pressure Ta = the air temperature and RH = the relative humidity.
The CWSI empirical value was subsequently calculated from the relationship between VPD and the corresponding canopy and the air temperature difference (Tc − Ta ) was calculated from the relationship described previously in Equation (4).
Weekly disparities were observed in both summer and autumn with regard to VPD, as were differences in the canopy temperature and the air temperatures measured during the period. The average seasonal Tc − Ta (°C) values recorded in summer 2021 were 2.23, 2.41, 2.57, and 2.70 for treatments 80–100% FC, 70–75% FC, 60–65% FC, and 50–55% FC, respectively, with their corresponding average VPD (kPa), monitored for the same period, at 1.31, 1.47, 1.58, and 1.69, respectively. In autumn, the observed mean seasonal Tc − Ta (°C) values measured for the treatments 80–100% FC, 70–75% FC, 60–65% FC, and 50–55% FC were 2.21, 2.36, 2.50, and 2.63, respectively, with their corresponding VPD at 1.28, 1.44, 1.53, and 1.65, respectively. Figure 5A–D shows the weekly disparities in Tc − Ta, as well as the VPD during the study period.
The baseline equations for the mandarins were established as Tc Ta = 0.56 × VPD + 4.05 for the non-water stress baseline and Tc Ta = 4.05 for non-transpiring (i.e., extremely water stressed conditions) baseline equations. During the establishment of these baseline equations, all data collected were put together to get a fair representation of the weather data available. Correlating VPD and Tc  Ta datasets for the non-transpiring treatments gave rise to the lower baseline equation and VPD against Tc  Ta data values for the extremely stressed treatments resulted in the development of the upper baseline equation. Tc  Ta had a significant correlation with VPD at R2 = 0.58 for the non-stressed treatments while it was weakly correlated to VPD at R2 = 0.005 for the highly stressed treatments. These are shown in Figure 6A,B.
The computed seasonal CWSI(Idso) values ranged from −0.21–0.79, with average seasonal values of 0.14, 0.33, 0.48, and 0.61 recorded for treatments 80–100% FC, 70–75% FC, 60–65% FC, and 50–55% FC, respectively, in summer, while the average seasonal values for autumn were 0.16, 0.32, 0.46, and 0.58 for the treatments 80–100% FC, 70–75% FC, 60–65% FC, and 50–55% FC, respectively. These values are shown in Figure 7A,B.
Both the VPD and the measured canopy-air temperature differences were significantly different (at p < 0.001) among the treatments and within treatments; however, both parameters were strongly correlated to the CWSI, as determined from the Idso procedures at R2 = 0.67 and R2 = 0.79, respectively, as shown in Figure 8A,B.

3.4. Mandarin CWSI Estimated from the Workswell WIRIS Agro R Infrared Camera

The results obtained from the WWARIC in estimating the mandarin CWSI followed a similar pattern and showed consistency with the results when the CWSI by Idso procedures were adopted to calculate the mandarin CWSI. Figure 9, Figure 10 and Figure 11 show a digital image of the mandarin orchard captured using the WWARIC, the isothermal graph, and the CWSI image, as well as a 3D graph gained from the WWARIC after analysis with the WIRIS core player version 1.4.1 during a data collection session in the autumn of 2021.
The outcome of the estimated CWSI from using the WWARIC on the four water application levels was established in autumn 2021 by categorizing the CWSI values into a 0.25 interval range, comparable to the procedures proposed by [24,53,54]. The CWSI recorded during the measurement session under the 80–100% FC treatment was within the CWSI range of 0.00 < CWSIW  0.25, which is shown by a deep green color, while treatments under the 70–75% FC water stress level were within the range 0.25 < CWSIW  0.5, represented with a light green color. In addition, the CWSI for mandarins under the 60–65% FC water stress level was within the 0.5 < CWSIW  0.75 range, which is represented by a yellow color, while the CWSI for mandarins under the 50–55% FC water stress level was within the range of 0.75 < CWSIW  1.00, which is represented with a light brown color. Nevertheless, the CWSI value > 1, indicated with the dark brown color on the TIFF image, shows the mulch cover, but not necessarily water stress within the canopy, as seen in Figure 10.

3.5. CWSI Based on the Idso Method and the WWARIC

The outcome of this research showed a high relationship between the CWSI estimate based on the Idso procedure, by relying on the field measurements of the mandarin canopy temperatures measured via the portable Raytek® handheld, non-contact, infrared thermometer, and the CWSI estimated from using the WWARIC, during the study period with a correlation coefficient of R2 = 0.75, as seen in Figure 12.

3.6. The CWSI and Plant-Based Water Stress Indicators (RLWC, Leaf Temperature, and Chlorophyll Content)

3.6.1. The CWSI Calculated by the Idso Method and Plant-Based Water Stress Indicators (RLWC, Leaf Temperature, and Chlorophyll Content)

The research revealed a positive connection between relative leaf water content, the chlorophyll content, and the soil water content, while an inverse relationship existed between the soil water level and the measured leaf temperatures. Generally, increasing the soil water content resulted in an increase in both the RWLC and the chlorophyll content, and vice versa. Mandarins subjected to no stress (80–100% FC) recorded the highest values for both the RWLC and the chlorophyll content, and vice versa. However, mandarin leaves subjected to this same treatment (i.e., 80–100% FC) recorded the lowest temperature measurements. Nonetheless, all three of the studied water stress indicators (RLWC, LT, and LCC) were significantly correlated with the CWSI estimated by the Idso approach at R2 = 0.68, 0.51, and 0.80, respectively, as seen in Figure 13A–C. Nevertheless, the calculated CWSI(Idso) was inversely related to both the RWLC and the leaf chlorophyll content, but directly related to the leaf temperature. An increase in the CWSI(Idso) corresponded to lower values of the RWLC and the measured leaf chlorophyll, and vice versa. On the other hand, the CWSI calculated by the Idso procedure increased with the increasing leaf temperatures recorded.

3.6.2. The CWSI Estimated from the WWARIC and Plant-Based Water Stress Indicators (RLWC and Leaf Temperature)

The research showed that the CWSI estimated with the WWARIC was significantly correlated with the RLWC and the leaf temperature, with R2 values of 0.77 and 0.79, respectively. This correlation is shown in Figure 14A,B. The mandarins’ RLWC, and the estimated CWSI from the WWARIC on mandarins from the field, depicted a negative relationship. It was observed that at high values of mandarins’ RLWC, there were corresponding lower values of the CWSI estimated from the WWARIC, and vice versa. Furthermore, the recorded leaf temperatures varied directly with the values of the estimated CWSI gained with the WWARIC, leading to a positive significant relationship. It was observed that at low measured leaf temperatures, the corresponding CWSI values gained from the camera were also low, and vice versa.

4. Discussion

The crop water stress index has been useful in numerous agronomic researches to determine plant stress caused by water and/or soil-water deficit conditions [12,13,25,26,52,55]. This concept was first established and proposed by Idso and his research team, who adopted and worked with environmental factors, such as prevailing vapor deficit pressures and air temperatures [36]. The Idso method has been subsequently adopted, modified, and applied in other related studies to monitor plant-water stress related conditions. The results from these previous studies indicated a very strong relationship between the CWSI calculated using the Idso approach and the CWSI calculated using other procedures [25,41,56,57,58]. Nonetheless, the arduous nature of obtaining the required metrological data, the difficulty in computations, and the high complexity, advanced technicality, and time involved in estimating the crop water stress index by the Idso approach made it challenging at the farm level to assess crop water and soil-water conditions to efficiently control drought and yield loss; hence, the need arose for most agronomists to resort to using remotely sensed data obtained from very high-resolution cameras to resolve these problems [49,52,56]. These studies showed that the baseline equations established for mandarins ( Tc Ta = 0.56 × VPD + 4.02 ) and ( Tc Ta = 4.02 ) for the non-water stress baseline and for the non-transpiring (extremely water stressed conditions) baseline equations, respectively, were similar to the equations developed by [59]. Nevertheless, although the linear regressions developed for the non-water stress baseline from the Tc  Ta and VPD datasets were significant at R2 = 0.58 (p < 0.05), the observed data were scattered. This outcome was consistent with the findings from [59,60,61], which observed that differences in the canopy-air temperatures and the VPD values measured for mandarin to develop the non-water stress baseline have a more scattering characteristic when compared to values of Tc  Ta and VPD recorded for developing non-water stress baseline equations for different orchard crops, such as sweet lime and pistachio. These authors attributed this phenomenon to variations in the mandarin canopy temperature.
The results further revealed that the CWSI computed using the Idso approach increased with an increasing soil-water deficit. The average estimated CWSI from the Idso approach for the period of study was in the order of 80–100% FC < 70–75% FC < 60–65% FC < 50–55% FC, which corresponded to CWSI values of 0.15 < 0.32 < 0.47 < 0.60. These figures are not different from the results of [27,52,54,61,62], where it was observed that the CWSI in plants increases with a decreasing availability of water in the soil (or increasing water stress). Furthermore, the CWSI calculated from the Idso approach was correlated with both the VPD and the difference in canopy-air temperature. This is consistent with the results from [62,63], where it was argued that the CWSI calculated by the Idso approach was a more sensitive index to the VPD and the canopy-air temperature difference, compared to other indices such as degree of non-stress and degrees above canopy threshold. Furthermore, the inverse relationship between estimated CWSI by the Idso approach and the RWLC from our studies is similar to the results from the works of [64,65,66], in their studies of soybean and corn. The leaf temperature was also directly related to the calculated CWSI, which is also consistent with studies from [59,62,67]. The effects of water stress on the leaf chlorophyll content were observed during the study period. This result showed a strong correlation between the calculated CWSI and the measured leaf chlorophyll content. This correlation is not different from the reports in [66,68,69,70,71], where it was observed that increasing crop water stress decreases the leaf chlorophyll content in a plant.
The strong association between the estimated CWSI from the field data and the CWSI gained by the WWARIC is an indication that the camera is capable of estimating the level of stress a particular plant is subjected to in real time and also in any plant setting. This finding is comparable to the studies in [24], where the WWARIC was used to successfully study water stress in tomato grown under greenhouse conditions. In estimating conventional CWSI by any method, a meteorological dataset is required before computations can be done. This requirement limits its use in modern-day precision agriculture, where the Internet of Things and big data systems, as well as artificial intelligence, make life easy. It is with this effect that the WWARIC has been endowed with all the features making it possible and easy to measure the correct prevailing air and canopy temperatures, humidity, and other parameters simultaneously to estimate water stress precisely in plants under any plant condition, from small fields to vast open fields and/or inside greenhouses in real time Other commercial infrared cameras and sensors do not have these features, providing the WWARIC with an advantage over them as far as in situ, instantaneous, and precise CWSI measurements are concerned. Furthermore, the works of [72,73,74,75,76] indicated that most thermal infrared cameras and sensors require surrogate calibrations before temperature data can be retrieved and utilized for subsequent determination of the CWSI. The WWARIC needs no calibration before it can measure the temperature and determine the CWSI of the plant under consideration. In addition, the in-built color map feature allows for easy identification of the temperature and the CWSI values of plants within specific regions of the field. It is worth noting that all three modes of the camera, differential, empirical and theoretical, were able to detect water stress in the mandarin orchard at heights above 12 m. The WWARIC is affected by the clouds, as the presence of a cloudy weather reduces the reflections on the plant canopy, thereby obscuring clear image acquisition; hence, it is advisable to use the camera on a clear sunny day, taking the measurements for calculating the CWSI with the Idso procedure.

5. Conclusions

The use of thermal indices based on canopy temperature for crop water stress surveillance and irrigation management has been established by numerous scientists in the last few decades, following the introduction of portable infrared thermometers. Since then, scientists have designed various scientific equipment and tools for standardizing the effects of some variable meteorological parameters, including VPD, wind speed, solar radiation, air temperature, etc. This novel research utilized a commercial WWARIC to remotely monitor mandarin CWSI precisely and instantaneously in real time, to swiftly model and schedule irrigation appropriately. Using the WWARIC, the predictable CWSI in mandarin development at each water deficit treatment level was shown by the following color codes: deep green representing the treatment 80–100% FC (little or no stress) within a range of 0 < CWSI 0.25 ; light green showing 75–70% FC (mild to medium water stress level) for the range 0 .25 < CWSI 0.5 ; yellow representing 65–60% FC (medium to high water stress level) within the range 0 .5 < CWSI 0.75 ; and light brown depicting 55–50% FC (high to extremely high water stress level) for the range 0 .75 < CWSI 1 . In contrast, dark brown colors, seen on the TIFF image, represent complete drought with no water present ( CWSI > 1 ) , which in this case was the image of the ground plastic mulch cover, since there was not a total drought situation in the study field. In this research, all three pixel modes (empirical, theoretical, and differential) of the WWARIC were used to establish the CWSI at a height of 12 m. The CWSIW values found in this study were consistent with the results obtained from the CWSII and with the relative leaf water content (RLWC), as well as with the leaf chlorophyll content of the mandarin plant. However, the WWARIC was designed to estimate CWSI (a standardized index established to compute water stress) while subduing the effects of environmental factors that affect the relationship between stresses, canopy temperature, and air temperature, unlike other infrared sensors and cameras. This makes this WWARIC a better option in determining the CWSI in plants, while minimizing the effects of several environmental factors. The CWSI obtained using the WWARIC is suitable for monitoring crop water conditions in real time to overcome drought and making immediate management decisions to reschedule irrigation earlier to avoid plant death and yield loss, while maximizing water use efficiency for mandarin production. This novel use of the WWARIC has demonstrated sufficient precision, swiftness, and intelligibility in the real-time measurement of mandarin water stress, depicting estimated CWSI in specific color ranges and values. When adopted and used in precision agriculture, this camera will save significant time, money, and other resources.

Author Contributions

Conceptualization, S.A.A.; data curation, K.E.A., A.A.A., H.W. and S.Q.; formal analysis, S.A.A., E.A.A., E.A., F.I. and R.O.D.; funding acquisition, J.L. and Y.L.; investigation, S.A.A., K.E.A., A.A.A., H.W. and S.Q.; methodology, S.A.A., K.E.A., A.A.A., H.W. and S.Q.; project administration, S.A.A. and J.L.; software, S.A.A., K.E.A., E.A.A., E.A., F.I. and R.O.D.; supervision, J.L. and Y.L.; validation, S.A.A.; writing—original draft, S.A.A.; writing—review and editing, S.A.A., K.E.A., A.A.A., H.W., S.Q., E.A.A., E.A., F.I. and R.O.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported with funding from the Key Area Research and Development Program of the Guangdong Province (2019B020214003).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Acknowledgments

We would like to acknowledge the efforts of everyone who contributed to the success of this article, especially the funding institution, the field workers, and the laboratory crew.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Li, J.; Chen, J. Citrus Fruit-Cracking: Causes and Occurrence. Hortic. Plant J. 2017, 3, 255–260. [Google Scholar] [CrossRef]
  2. Blanco, V.; Blaya-Ros, P.J.; Torres-Sánchez, R.; Domingo, R. Influence of Regulated Deficit Irrigation and Environmental Conditions on Reproductive Response of Sweet Cherry Trees. Plants 2020, 9, 94. [Google Scholar] [CrossRef] [Green Version]
  3. Saitta, D.; Consoli, S.; Ferlito, F.; Torrisi, B.; Allegra, M.; Longo-Minnolo, G.; Ramírez-Cuesta, J.M.; Vanella, D. Adaptation of citrus orchards to deficit irrigation strategies. Agric. Water Manag. 2021, 247, 106734. [Google Scholar] [CrossRef]
  4. Kuscu, H.; Karasu, A.; Oz, M.; Demir, O.A.; Turgut, I. Effect of irrigation amounts applied with drip irrigation on maize evaporation, yield, water use efficiency and net return in a sub-humid climate. Turk. J. Field Crops 2013, 18, 13–19. [Google Scholar]
  5. Al-Ghobari, H.M.; Dewidar, A.Z. Integrating deficit irrigation into surface and subsurface drip irrigation as a strategy to save water in arid regions. Agric. Water Manag. 2018, 209, 55–61. [Google Scholar] [CrossRef]
  6. Fernandes-Silva, A.; Oliveira, M.; Paço, T.A.; Ferreira, I. Deficit Irrigation in Mediterranean Fruit Trees and Grapevines: Water Stress Indicators and Crop Responses. In Irrigation in Agroecosystems; Ondrašek, G., Ed.; IntechOpen: London, UK, 2018; pp. 52–85. [Google Scholar]
  7. Faghih, S.; Zamani, Z.; Fatahi, R.; Liaghat, A. Effects of deficit irrigation and kaolin application on vegetative growth and fruit traits of two early ripening apple cultivars. Biol. Res. 2019, 52, 43. [Google Scholar] [CrossRef] [Green Version]
  8. Fernández, J.E. Plant-based methods for irrigation scheduling of woody crops. Horticulturae 2017, 3, 35. [Google Scholar] [CrossRef] [Green Version]
  9. Romero-Trigueros, C.; Gambín, J.M.B.; Tortosa, P.A.N.; Cabañero, J.J.A.; Nicolás, E.N. Determination of Crop Water Stress Index by Infrared Thermometry in Grapefruit Trees Irrigated with Saline Reclaimed Water Combined with Deficit Irrigation. Remote Sens. 2019, 11, 757. [Google Scholar] [CrossRef] [Green Version]
  10. Ihuoma, S.O. The Use of Spectral Reflectance Data to Assess Plant Stress and Improve Irrigation Water Management. Ph.D. Thesis, McGill University, Montreal, Canada, 2020. [Google Scholar]
  11. Jones, H.G. Irrigation scheduling: Advantages and pitfalls of plant based methods. J. Exp. Bot. 2004, 55, 2427–2436. [Google Scholar] [CrossRef] [Green Version]
  12. Drechsler, K.; Kisekka, I.; Upadhyaya, S. A comprehensive stress indicator for evaluating plant water status in almond trees. Agric. Water Manag. 2019, 216, 214–223. [Google Scholar] [CrossRef]
  13. Berni, J.A.J.; Zarco-Tejada, P.J.; Sepulcre-Cantó, G.; Fereres, E.; Villalobos, F. Mapping canopy conductance and CWSI in olive orchards using high resolution thermal remote sensing imagery. Remote Sens. Environ. 2009, 113, 2380–2388. [Google Scholar] [CrossRef]
  14. Bastiaanssen, W.G.M.; Pelgrum, H.; Soppe, R.W.O.; Allen, R.G.; Thoreson, B.P.; de, C. Teixeira, A.H. Thermal-infrared technology for local and regional scale irrigation analyses in horticultural systems. Acta Hortic. 2006, 792, 33–46. [Google Scholar]
  15. Berni, J.A.; Zarco-Tejada, P.J.; Suárez, L.; Fereres, E. Thermal and narrowband multispectral remote sensing for vegetation monitoring from an unmanned aerial vehicle. IEEE Trans. Geosci. Remote Sens. 2009, 47, 722–738. [Google Scholar] [CrossRef] [Green Version]
  16. Stagakis, S.; González-Dugo, V.; Cid, P.; Guillén-Climent, M.L.; Zarco-Tejada, P.J. Monitoring water stress and fruit quality in an orange orchard under regulated deficit irrigation using narrow-band structural and physiological remote sensing indices. ISPRS J. Photogramm. Remote Sens. 2012, 71, 47–61. [Google Scholar] [CrossRef] [Green Version]
  17. Martínez-Gimeno, M.A.; Bonet, L.; Provenzano, G.; Badal, E.; Intrigliolo, D.S.; Ballester, C. Assessment of yield and water productivity of clementine trees under surface and subsurface drip irrigation. Agric. Water Manag. 2018, 206, 209–216. [Google Scholar] [CrossRef]
  18. Romero-Trigueros, C.; Nortes, P.A.; Alarcón, J.J.; Hunink, J.E.; Parra, M.; Contreras, S.; Droogers, P.; Nicolás, E. Effects of saline reclaimed waters and deficit irrigation on Citrus physiology assessed by UAV remote sensing. Agric. Water Manag. 2017, 183, 60–69. [Google Scholar] [CrossRef] [Green Version]
  19. Zhang, F.; Zhou, G. Estimation of vegetation water content using hyperspectral vegetation indices: A comparison of crop water indicators in response to water stress treatments for summer maize. BMC Ecol. 2019, 19, 18. [Google Scholar] [CrossRef] [Green Version]
  20. Ihuoma, S.O.; Madramootoo, C.A. Sensitivity of spectral vegetation indices for monitoring water stress in tomato plants. Comput. Electron. Agric. 2019, 163, 104860. [Google Scholar] [CrossRef]
  21. Stimson, H.C.; Breshears, D.D.; Ustin, S.L.; Kefauver, S.C. Spectral sensing of foliar water conditions in two co-occurring conifer species: Pinus edulis and Juniperus monosperma. Remote Sens. Environ. 2005, 96, 108–118. [Google Scholar] [CrossRef]
  22. Jones, H.G.; Stoll, M.; Santos, T.; de Sousa, C.; Chaves, M.M.; Grant, O.M. Use of infrared thermography for monitoring stomatal closure in the field: Application to grapevine. J. Exp. Bot. 2002, 53, 2249–2260. [Google Scholar] [CrossRef]
  23. Merlot, S.; Mustilli, A.C.; Genty, B.; North, H.; Lefebvre, V.; Sotta, B.; Vavasseur, A.; Giraudat, J. Use of infrared thermal imaging to isolate Arabidopsis mutants defective in stomatal regulation. Plant J. 2002, 30, 601–609. [Google Scholar] [CrossRef] [Green Version]
  24. Alordzinu, K.E.; Li, J.; Lan, Y.; Appiah, S.A.; Al Aasmi, A.; Wang, H. Rapid Estimation of Crop Water Stress Index on Tomato Growth. Sensors 2021, 21, 5142. [Google Scholar] [CrossRef] [PubMed]
  25. Poblete-Echeverría, C.; Espinace, D.; Sepúlveda-Reyes, D.; Zuñiga, M.; Sanchez, M. Analysis of crop water stress index (CWSI) for estimating stem water potential in grapevines: Comparison between natural reference and baseline approaches. Acta Hortic. 2017, 1150, 189–194. [Google Scholar] [CrossRef]
  26. Ru, C.; Hu, X.; Wang, W.; Ran, H.; Song, T.; Guo, Y. Evaluation of the Crop Water Stress Index as an Indicator for the Diagnosis of Grapevine Water Deficiency in Greenhouses. Horticulturae 2020, 6, 86. [Google Scholar] [CrossRef]
  27. Jackson, R.D.; Idso, S.B.; Reginato, R.J.; Pinter, P.J. Canopy temperature as a crop water stress indicator. Water Resour. Res. 1981, 17, 1133–1138. [Google Scholar] [CrossRef]
  28. García-Tejero, I.; Durán-Zuazo, V.H.; Muriel-Fernández, J.L.; Martínez-García, G.; Jiménez-Bocanegra, J.A. Benefits of low-frequency irrigation in citrus orchards. Agron. Sustain. Dev. 2011, 31, 779–791. [Google Scholar] [CrossRef]
  29. Stagno, F.; Roccuzzo, G.; Allegra, M.; Intrigliolo, F.; Parisi, R.; Cirelli, G.; Consoli, S.; Barbagallo, S. Deficit irrigation strategies: Preliminary assessment on a sicilian young orange orchard. Acta Hortic. 2015, 1065, 1713–1718. [Google Scholar] [CrossRef]
  30. Consoli, S.; Stagno, F.; Vanella, D.; Boaga, J.; Cassiani, G.; Roccuzzo, G. Partial root-zone drying irrigation in orange orchards: Effects on water use and crop production characteristics. Eur. J. Agron. 2017, 82, 190–202. [Google Scholar] [CrossRef]
  31. Consoli, S.; Stagno, F.; Roccuzzo, G.; Cirelli, G.L.; Intrigliolo, F. Sustainable management of limited water resources in a young orange orchard. Agric. Water Manag. 2014, 132, 60–68. [Google Scholar] [CrossRef]
  32. Wright, G.C. Irrigating Citrus Trees; University of Arizona: Tucson, AZ, USA, 2000. [Google Scholar]
  33. Swain, S. Growing Citrus in the North Bay; University of California: Berkeley, CA, USA, 2012. [Google Scholar]
  34. Carr, M.K.V. The Water Relations and Irrigation Requirements of Citrus (Citrus Spp.): A Review. Exp. Agric. 2012, 48, 347–377. [Google Scholar] [CrossRef]
  35. Khedkar, D.D. Minimal Volume Irrigation Scheduling for Citrus Cultivation in Warud Tehsil, Maharashtra, India. Agric. Agric. Sci. Procedia 2015, 4, 365–371. [Google Scholar] [CrossRef]
  36. Idso, S.B.; Jackson, R.D.; Pinter, P.J.; Reginato, R.J.; Hatfield, J.H. Normalizing the stress degree-day parameter for environmental variability. Agric. Meteorol. 1981, 24, 45–55. [Google Scholar] [CrossRef]
  37. Idso, S.B. Non-water stressed baselines: A key to measuring and interpreting plant water stress. Agric. Meteorol. 1982, 27, 59–70. [Google Scholar] [CrossRef]
  38. Cai, L.; Xu, Z.; Ren, M.; Guo, Q.; Hu, X.; Hu, G.; Wan, H.; Peng, P. Source identification of eight hazardous heavy metals in agricultural soils of Huizhou, Guangdong Province, China. Ecotoxicol. Environ. Saf. 2012, 78, 2–8. [Google Scholar] [CrossRef]
  39. Fahad, M.; Wajid, S.A.; Ahmad, A.; Cheema, M.J.M. Response of wheat cultivars to deficit irrigation under semiarid conditions of Faisalabad. Int. J. Agric. Biol. 2019, 21, 1004–1012. [Google Scholar]
  40. Khan, S.; Rasool, A.; Irshad, S.; Hafeez, M.B.; Ali, M.; Saddique, M.; Asif, M.; Hasnain, Z.; Naseem, S.; Alwahibi, M.S.; et al. Potential soil moisture deficit: A useful approach to save water with enhanced growth and productivity of wheat crop. J. Water Clim. Change 2021, 12, 2515–2525. [Google Scholar] [CrossRef]
  41. Cao, Z.; Wang, Q.; Zheng, C. Best hyperspectral indices for tracing leaf water status as determined from leaf dehydration experiments. Ecol. Indic. 2015, 54, 96–100. [Google Scholar] [CrossRef]
  42. Fuchs, M. Infrared measurement of canopy temperature and detection of plant water stress. Theor. Appl. Clim. 1990, 42, 253–261. [Google Scholar] [CrossRef]
  43. Ihuoma, S.O.; Madramootoo, C.A. Narrow-band reflectance indices for mapping the combined effects of water and nitrogen stress in tomato crops. Biosyst. Eng. 2020, 192, 133–143. [Google Scholar] [CrossRef]
  44. Rashid, A.; Stark, J.C.; Tanveer, A.; Mustafa, T. Use of Canopy Temperature Measurements as a Screening Tool for Drought Tolerance in Spring Wheat. J. Agron. Crop Sci. 1999, 182, 231–237. [Google Scholar] [CrossRef]
  45. Pieczynski, M.; Marczewski, W.; Hennig, J.; Dolata, J.; Bielewicz, D.; Piontek, P.; Wyrzykowska, A.; Krusiewicz, D.; Strzelczyk-Zyta, D.; Konokpa-Postupolska, D.; et al. Down-regulation of CBP80 gene expression as a strategy to engineer a drought-tolerant potato. Plant Biotechnol. J. 2013, 11, 459–469. [Google Scholar] [CrossRef] [PubMed]
  46. Zhu, J.; Tremblay, N.; Liang, Y. A corn nitrogen status indicator less affected by soil water content. Agron. J. 2011, 103, 890–898. [Google Scholar] [CrossRef]
  47. Jamshidi, S.; Zand-Parsa, S.; Kamgar-Haghighi, A.A.; Shahsavar, A.R.; Niyogi, D. Evapotranspiration, crop coefficients, and physiological responses of citrus trees in semi-arid climatic conditions. Agric. Water Manag. 2020, 227, 105838. [Google Scholar] [CrossRef]
  48. Jamshidi, S.; Zand-Parsa, S.; Niyogi, D. Assessing Crop Water Stress Index of Citrus Using In-Situ Measurements, Landsat, and Sentinel-2 Data. Int. J. Remote Sens. 2021, 42, 1893–1916. [Google Scholar] [CrossRef]
  49. Katsoulas, N.; Elvanidi, A.; Ferentinos, K.P.; Kacira, M.; Bartzanas, T.; Kittas, C. Crop reflectance monitoring as a tool for water stress detection in greenhouses: A review. Biosyst. Eng. 2016, 151, 374–398. [Google Scholar] [CrossRef]
  50. Paltineanu, C.; Septar, L.; Moale, C. Crop Water Stress in Peach Orchards and Relationships with Soil Moisture Content in a Chernozem of Dobrogea. J. Irrig. Drain. Eng. 2013, 139, 20–25. [Google Scholar] [CrossRef]
  51. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration: Guidelines for Computing Crop Water Requirements; FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998. [Google Scholar]
  52. Tanriverdi, C.; Atilgan, A.; Degirmenci, H.; Akyuz, A. Comparasion of Crop Water Stress Index (CWSI) and Water Deficit Index (WDI) by Using Remote Sensing (RS). Infrastruct. Ecol. Rural Areas 2017, 879–894. [Google Scholar] [CrossRef]
  53. Magney, T.S.; Vierling, L.A.; Eitel, J.U.; Huggins, D.R.; Garrity, S.R. Response of high frequency Photochemical Reflectance Index (PRI) measurements to environmental conditions in wheat. Remote Sens. Environ. 2016, 173, 84–97. [Google Scholar] [CrossRef]
  54. Veysi, S.; Naseri, A.A.; Hamzeh, S.; Bartholomeus, H. A satellite based crop water stress index for irrigation scheduling in sugarcane fields. Agric. Water Manag. 2017, 189, 70–86. [Google Scholar] [CrossRef]
  55. Yang, Y.; Ling, P.; Fleisher, D.H.; Timlin, D.J.; Reddy, V.R. Non-contacting techniques for plant drought stress detection. Trans. ASABE 2008, 51, 1483–1492. [Google Scholar] [CrossRef]
  56. Espinoza, C.Z.; Khot, L.R.; Sankaran, S.; Jacoby, P.W. High resolution multispectral and thermal remote sensing-based water stress assessment in subsurface irrigated grapevines. Remote Sens. 2017, 9, 961. [Google Scholar] [CrossRef] [Green Version]
  57. Gerhards, M.; Schlerf, M.; Mallick, K.; Udelhoven, T. Challenges and Future Perspectives of Multi- Hyperspectral Thermal Infrared Remote Sensing for Crop Water-Stress Detection: A Review. Remote Sens. 2019, 11, 1240. [Google Scholar] [CrossRef] [Green Version]
  58. Zhang, F.; Zhou, G. Deriving a light use efficiency estimation algorithm using in situ hyperspectral and eddy covariance measurements for a maize canopy in Northeast China. Ecol. Evol. 2017, 7, 4735–4744. [Google Scholar] [CrossRef] [PubMed]
  59. Gonzalez-Dugo, V.; Zarco-Tejada, P.J.; Fereres, E. Applicability and limitations of using the crop water stress index as an indicator of water deficits in citrus orchards. Agric. For. Meteorol. 2014, 198–199, 94–104. [Google Scholar] [CrossRef]
  60. Testi, L.; Goldhamer, D.; Iniesta, F.; Salinas, M. Crop water stress index is a sensitive water stress indicator in pistachio trees. Irrig. Sci. 2008, 26, 395–405. [Google Scholar] [CrossRef] [Green Version]
  61. Sepaskhah, A.R.; Kashefipour, S.M. Evapotranspiration and crop coefficient of sweet lime under drip irrigation. Agric. Water Manag. 1995, 27, 331–340. [Google Scholar] [CrossRef]
  62. DeJonge, K.C.; Taghvaeian, S.; Trout, T.J.; Comas, L.H. Comparison of canopy temperature-based water stress indices for maize. Agric. Water Manag. 2015, 156, 51–62. [Google Scholar] [CrossRef]
  63. Taghvaeian, S.; Comas, L.; DeJonge, K.C.; Trout, T.L. Conventional and simplified canopy temperature indices predict water stress in sunflower. Agric. Water Manag. 2014, 144, 69–80. [Google Scholar] [CrossRef]
  64. Dahlgren, R.P.; Vanderbilt, V.C.; Daughtry, C.S.T. Estimates of Leaf Relative Water Content from Optical Polarization Measurements. In Proceedings of the 2017 American Geophysical Union Fall Meeting, New Orleans, LA, USA, 11–15 December 2017; Volume 2017. Abstract A21B-2160. [Google Scholar]
  65. Swain, S.; Rundquist, D.; Arkebauer, T.; Narumalani, S.; Wardlow, B. Non-invasive estimation of relative water content in soybean leaves using infrared thermography. Isr. J. Plant Sci. 2012, 60, 25–36. [Google Scholar] [CrossRef]
  66. Pirzad, A.; Shakiba, M.R.; Zehtab-Salmasi, S.; Mohammadi, S.A.; Darvishzadeh, R.; Samadi, A. Effect of water stress on leaf relative water content, chlorophyll, proline and soluble carbohydrates in Matricaria chamomilla L. J. Med. Plants Res. 2011, 5, 2483–2488. [Google Scholar]
  67. Carroll, D.A.; Hansen, N.C.; Hopkins, B.G.; DeJonge, K.C. Leaf temperature of maize and Crop Water Stress Index with variable irrigation and nitrogen supply. Irrig. Sci. 2017, 35, 549–560. [Google Scholar] [CrossRef]
  68. Sanchez, R.A.; Hall, A.J.; Trapani, N.; de Hunau, R.C. Effects of water stress on the chlorophyll content, nitrogen level and photosynthesis of leaves of two maize genotypes. Photosynth. Res. 1983, 4, 35–47. [Google Scholar] [CrossRef] [PubMed]
  69. Chutia, J.; Borah, S.P. Water Stress Effects on Leaf Growth and Chlorophyll Content but Not the Grain Yield in Traditional Rice (Oryza sativa Linn.) Genotypes of Assam, India II. Protein and Proline Status in Seedlings under PEG Induced Water Stress. Am. J. Plant Sci. 2012, 3, 971–980. [Google Scholar] [CrossRef] [Green Version]
  70. Khayatnezhad, M.; Gholamin, R. The effect of drought stress on leaf chlorophyll content and stress resistance in maize cultivars (Zea mays). Afr. J. Microbiol. Res. 2012, 6, 2844–2848. [Google Scholar] [CrossRef]
  71. Shahenshah; Isoda, A. Effects of Water Stress on Leaf Temperature and Chlorophyll Fluorescence Parameters in Cotton and Peanut. Plant Prod. Sci. 2010, 13, 269–278. [Google Scholar] [CrossRef]
  72. Jensen, A.M.; McKee, M.; Chen, Y. Procedures for processing thermal images using low-cost microbolometer cameras for small unmanned aerial systems. In Proceedings of the 2014 IEEE Geoscience and Remote Sensing Symposium, Quebec City, QC, Canada, 13–18 July 2014; IEEE: Piscataway, NJ, USA, 2014; pp. 2629–2632. [Google Scholar]
  73. Khanal, S.; Fulton, J.; Shearer, S. An overview of current and potential applications of thermal remote sensing in precision agriculture. Comput. Electron. Agric. 2017, 139, 22–32. [Google Scholar] [CrossRef]
  74. Maes, W.; Huete, A.; Steppe, K. Optimizing the processing of UAV-based thermal imagery. Remote Sens. 2017, 9, 476. [Google Scholar] [CrossRef] [Green Version]
  75. Aragon, B.; Johansen, K.; Parkes, S.; Malbetéau, Y.; Al-Mashharawi, S.; Al-Amoudi, T.; Andrade, C.F.; Turner, D.; Lucieer, A.; McCabe, M.F. A calibration procedure for field and UAV-based uncooled thermal infrared instruments. Sensors 2020, 20, 3316. [Google Scholar] [CrossRef]
  76. Stutsel, B.; Johansen, K.; Malbéteau, Y.M.; McCabe, M.F. Detecting Plant Stress Using Thermal and Optical Imagery from an Unoccupied Aerial Vehicle. Front. Plant Sci. 2021, 12, 734944. [Google Scholar] [CrossRef]
Figure 1. Field layout and experimental design.
Figure 1. Field layout and experimental design.
Sensors 22 04018 g001
Figure 2. Correlation between soil volumetric water content and TDR readings.
Figure 2. Correlation between soil volumetric water content and TDR readings.
Sensors 22 04018 g002
Figure 3. Variations in soil water content.
Figure 3. Variations in soil water content.
Sensors 22 04018 g003
Figure 4. Seasonal variations in relative leaf water content.
Figure 4. Seasonal variations in relative leaf water content.
Sensors 22 04018 g004
Figure 5. (A) Weekly Tc − Ta (°C) for summer, 2021; (B) Weekly Tc  Ta (°C) for summer, 2021; (C) Weekly VPD (kPa) for autumn, 2021; (D) Weekly VPD (kPa) for autumn, 2021.
Figure 5. (A) Weekly Tc − Ta (°C) for summer, 2021; (B) Weekly Tc  Ta (°C) for summer, 2021; (C) Weekly VPD (kPa) for autumn, 2021; (D) Weekly VPD (kPa) for autumn, 2021.
Sensors 22 04018 g005
Figure 6. Relationship between Tc  Ta (°C) and VPD (kPa) for developing (A) non-water stressed and (B) non-transpiring baseline reference points in relation to fully irrigated and extremely stressed mandarin trees at noon on a very bright, clear day.
Figure 6. Relationship between Tc  Ta (°C) and VPD (kPa) for developing (A) non-water stressed and (B) non-transpiring baseline reference points in relation to fully irrigated and extremely stressed mandarin trees at noon on a very bright, clear day.
Sensors 22 04018 g006
Figure 7. (A) CWSI for summer 2021; (B) CWSI for autumn 2021.
Figure 7. (A) CWSI for summer 2021; (B) CWSI for autumn 2021.
Sensors 22 04018 g007
Figure 8. Correlation between (A) CWSI(Idso) and VPD (kPa); (B) CWSI(Idso) and Tc  Ta (°C).
Figure 8. Correlation between (A) CWSI(Idso) and VPD (kPa); (B) CWSI(Idso) and Tc  Ta (°C).
Sensors 22 04018 g008
Figure 9. (A) digital imagery of the mandarin orchard. (B) Isothermal panel showing histogram of pixels and their corresponding CWSI values from the WWARIC.
Figure 9. (A) digital imagery of the mandarin orchard. (B) Isothermal panel showing histogram of pixels and their corresponding CWSI values from the WWARIC.
Sensors 22 04018 g009aSensors 22 04018 g009b
Figure 10. Isothermal imagery of the mandarin orchard depicting various water stress levels and their corresponding range of CWSI values recorded by the WWARIC. (Subset (A) is the isothermal image of mandarin under low water stress level (B). isothermal image of mandarin subjected to mild water stress (C). isothermal image of mandarin under moderate water stress (D). isothermal image of mandarin under high water stress.
Figure 10. Isothermal imagery of the mandarin orchard depicting various water stress levels and their corresponding range of CWSI values recorded by the WWARIC. (Subset (A) is the isothermal image of mandarin under low water stress level (B). isothermal image of mandarin subjected to mild water stress (C). isothermal image of mandarin under moderate water stress (D). isothermal image of mandarin under high water stress.
Sensors 22 04018 g010
Figure 11. 3D graph of the thermal data obtained by the WWARIC.
Figure 11. 3D graph of the thermal data obtained by the WWARIC.
Sensors 22 04018 g011
Figure 12. Correlation between CWSI calculated from Idso procedures and CWSI from WWARIC.
Figure 12. Correlation between CWSI calculated from Idso procedures and CWSI from WWARIC.
Sensors 22 04018 g012
Figure 13. Correlation between (A) CWSI(Idso) and RLWC (%) and (B) CWSI(Idso) and leaf temperature (°C); (C) CWSI(Idso) and LCC (mg g−1).
Figure 13. Correlation between (A) CWSI(Idso) and RLWC (%) and (B) CWSI(Idso) and leaf temperature (°C); (C) CWSI(Idso) and LCC (mg g−1).
Sensors 22 04018 g013
Figure 14. Correlation between (A) CWSI(WWARIC) and RLWC (%) and (B) CWSI(WWARIC) and leaf temperature (°C).
Figure 14. Correlation between (A) CWSI(WWARIC) and RLWC (%) and (B) CWSI(WWARIC) and leaf temperature (°C).
Sensors 22 04018 g014
Table 1. General soil physical and chemical properties.
Table 1. General soil physical and chemical properties.
Soil Physical Properties
Depth (cm)Textural AnalysisField Capacity (%)Permanent Wilting Point (%)Bulk Density (gcm−3)
Sand (%)Clay (%)Silt (%)Textural Class
0–3072.8020.007.20Loamy sand22.60 ± 0.110.35 ± 0.11.47 ± 0.1
30–6073.8019.207.02Loamy sand24.20 ± 0.110.30 ± 0.11.48 ± 0.1
60–9069.70 19.1011.20Loamy sand22.10 ± 0.110.27 ± 0.11.53 ± 0.1
Average72.1019.438.47Loamy sand22.97 ± 0.110.31 ± 0.11.49 ± 0.1
Soil Chemical Properties
ParameterpHO.M (g/kg)Total N (g/kg)Total K (g/kg)Total P (g/kg)Alkalized N (mg/kg)Available P (mg/kg)Available K (mg/kg)
Values 4.908.830.422.270.3739.8185.9172.89
Table 2. Meteorological data during mandarin seasonal development.
Table 2. Meteorological data during mandarin seasonal development.
SeasonSummerAutumn
Air temperature (℃)32.528.0
Relative humidity (%)79.375.0
Vapor pressure deficit (kPa)1.51.5
Note: (1) Values shown are average seasonal values (22 June 2021 to 30 November 2021). (2) Summer (22 June 2021 to 22 September 2021); autumn (23 September 2021 to 20 December 2021).
Table 3. Irrigation water supplied for mandarin in different seasons in 2021.
Table 3. Irrigation water supplied for mandarin in different seasons in 2021.
SeasonIrrigation Period (Days)Irrigation Amount (mm)
80–100% FC70–75% FC60–65% FC50–55% FC
SummerJune (9 days)----
July (31 days)24.0517.4415.0312.63
August (31 days)57.2041.5035.7030.00
September (22 days)24.3918.6015.6013.70
Total for summer105.6477.5466.3356.33
AutumnSeptember (7 days)9.406.005.805.94
October (31 days)32.8024.0421.2017.96
November (30 days)36.8026.4423.0719.91
Total for autumn79.0056.4850.0743.81
Total irrigation amount184.64134.02116.40100.14
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Appiah, S.A.; Li, J.; Lan, Y.; Darko, R.O.; Alordzinu, K.E.; Al Aasmi, A.; Asenso, E.; Issaka, F.; Afful, E.A.; Wang, H.; et al. Real-Time Assessment of Mandarin Crop Water Stress Index. Sensors 2022, 22, 4018. https://doi.org/10.3390/s22114018

AMA Style

Appiah SA, Li J, Lan Y, Darko RO, Alordzinu KE, Al Aasmi A, Asenso E, Issaka F, Afful EA, Wang H, et al. Real-Time Assessment of Mandarin Crop Water Stress Index. Sensors. 2022; 22(11):4018. https://doi.org/10.3390/s22114018

Chicago/Turabian Style

Appiah, Sadick Amoakohene, Jiuhao Li, Yubin Lan, Ransford Opoku Darko, Kelvin Edom Alordzinu, Alaa Al Aasmi, Evans Asenso, Fuseini Issaka, Ebenezer Acheampong Afful, Hao Wang, and et al. 2022. "Real-Time Assessment of Mandarin Crop Water Stress Index" Sensors 22, no. 11: 4018. https://doi.org/10.3390/s22114018

APA Style

Appiah, S. A., Li, J., Lan, Y., Darko, R. O., Alordzinu, K. E., Al Aasmi, A., Asenso, E., Issaka, F., Afful, E. A., Wang, H., & Qiao, S. (2022). Real-Time Assessment of Mandarin Crop Water Stress Index. Sensors, 22(11), 4018. https://doi.org/10.3390/s22114018

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