Remote Sensing for Plant Water Content Monitoring: A Review

: This paper reviews the different remote sensing techniques found in the literature to monitor plant water status, allowing farmers to control the irrigation management and to avoid unnecessary periods of water shortage and a needless waste of valuable water. The scope of this paper covers a broad range of 77 references published between the years 1981 and 2021 and collected from different search web sites, especially Scopus. Among them, 74 references are research papers and the remaining three are review papers. The different collected approaches have been categorized according to the part of the plant subjected to measurement, that is, soil (12.2%), canopy (33.8%), leaves (35.1%) or trunk (18.9%). In addition to a brief summary of each study, the main monitoring technologies have been analyzed in this review. Concerning the presentation of the data, different results have been obtained. According to the year of publication, the number of published papers has increased exponentially over time, mainly due to the technological development over the last decades. The most common sensor is the radiometer, which is employed in 15 papers (20.3%), followed by continuous-wave (CW) spectroscopy (12.2%), camera (10.8%) and THz time-domain spectroscopy (TDS) (10.8%). Excluding two studies, the minimum coefﬁcient of determination (R 2 ) obtained in the references of this review is 0.64. This indicates the high degree of correlation between the estimated and measured data for the different technologies and monitoring methods. The ﬁve most frequent water indicators of this study are: normalized difference vegetation index (NDVI) backscattering coefﬁcients (10.8%), spectral reﬂectance (8.1%), reﬂection coefﬁcient (8.1%) and dielectric constant (8.1%).


Introduction
Water is one of the most important factors that have a significant influence on crop production. There is a direct relationship between biomass production and water consumed through transpiration [1]. Analyzing specific parts of a plant, leaf water content (LWC) is a crucial indicator of essential physiological processes [2], such as stomatal conductance (g s ), transpiration, photosynthesis and respiration [3], which is closely related to forest fire ignition and fire propagation [4]. Furthermore, irrigation currently represents the largest consumption of water on the planet, reaching about 20% of the total amount of consumed freshwater and two-thirds of the total amount destined for human beings [5]. In addition to these factors, climate change, which has been increasingly consolidated over the last decades, has induced water scarcity in many parts of the world, leading leaf optical depth (τ) and leaf surface area (LA) [58,59] and NDVI [67]) that also agree with LWC.
CWC can be also estimated by monitoring the water content of trunks or stems. Several studies use invasive methods to monitor CWC, consisting in the insertion of needles or probes inside the trunk or stems [71][72][73][74]76,81]. It is noteworthy to highlight that these approaches cannot be strictly considered as remote sensing techniques in spite of their widespread use in this research field. Other studies describe methods that need contact with the trunk or stems to develop the monitoring function [75,[78][79][80]. Among the main sensors that estimate the tree water status include coaxial probes called portable dielectric probes (PDP) [71][72][73][74][75], volumetric soil water content (VSWC) sensors inserted into tree trunks [76], sap-flow measurements based on sensors that measure sap-flow rate (SFR) (g·h -1 ) [78][79][80] or sap-flux density (SFD) (m 3 ·m -2 ·h -1 ) [77,81], CW vector network analyzer (VNA) spectroscopy [82,83], portable low-power cost-effective FMCW radars [84] and portable nuclear magnetic resonance (NMR) sensors [85]. In relation to the parameters measured by these sensors, PDPs, by using the complex reflection coefficient, estimate the complex xylem dielectric constant (XDC) to be correlated with the xylem SFD, the xylem sap chemical composition and xylem water potential (XWP). The VSWC sensors estimate XDC to be correlated with the stem water potential (SWP). Sap-flow measurements constitute an estimation of the tree water status. The most widely used techniques in order to measure sap-flow rate and SFD are the heat balance (HB) [78][79][80] and the heat pulse (HP) [77,81] methods. CW-VNA spectrometers measure the trunk reflection coefficient correlated with the VSWC (humidity probe) and with the trunk diameter variation measured by a dendrometer. Portable FMCW radars acquire the signal reflected by the xylem tissue to be correlated with the VSWC. Finally, NMR sensors measure the amplitude of the NMR signal, which scales linearly and quantitatively with the total amount of water contained within the sensor.
In this paper, we reviewed the different remote sensing techniques and approaches found in the literature, as well as their applications, for assessing plant water status. The scope of this paper covers a broad range of 77 scientific works published between the years 1981 and 2021, including mostly journal papers and some books and conference proceedings, and collected from different search web sites, especially Scopus. Different technologies and sensors used to conduct remote sensing are described in Section 2. Section 3 introduces the most relevant results of the collected references. Section 4 presents the discussion of results. Finally, Section 5 concludes this paper.

Sensors and Techniques
Throughout the review process, a wide variety of passive and active sensors used to monitor plant water status have been identified. Passive sensors measure radiation reflected by the tested object coming from a source different from the instrument or emitted by the object itself. Sunlight is the most frequent source of radiation detected by this type of devices. The most typical passive sensors include radiometers and spectrometers mounted on satellites and aircraft, cameras and portable spectrometers. Conversely, active sensors illuminate the object with their own electromagnetic radiation and then measure the radiation that is reflected, backscattered or transmitted through this object. The most common active sensors are radars, lasers, scatterometers, VSWC probes, PDPs, SFR and SFD gauges and some kinds of spectroscopy techniques, such as THz-TDS, THz-QTDS, CW spectroscopies, FTIR spectroscopy and air-coupled broadband ultrasonic spectroscopy. In the case of ultrasonic spectroscopy and sap-flow measurements, the electromagnetic radiation is replaced by acoustic waves and thermal energy respectively. In the following subsections, the most relevant sensors collected in this paper are briefly outlined.

Radiometer
This is an instrument that measures the radiant flux emitted or reflected by a given object and usually works in a portion of the spectrum selected among the visible, infrared or microwave bands [86]. Sometimes, it covers multiple bands simultaneously and is called multichannel radiometer. Satellites and manned aircraft are typically used as platforms.
In the case of microwave radiometers, it is very common to express the power in terms of an equivalent temperature, called brightness temperature ( Figure 1) and defined as the physical temperature of a blackbody that would radiate the same power at the considered wavelength. The output of a microwave radiometer is actually a voltage value proportional to its antenna temperature (T A ) and the final goal of the measurement is to relate this antenna temperature to the brightness temperature of the object [87]. In the literature, different methods to retrieve SMC and VWC on the basis of T B have been found, such as the direct correlation between water content and T B [9] or T B normalized by the land surface temperature [9], the calculation of specific indices sensitive to water content [9,19,22], the application of retrieval models [18,23] and the implementation of decision-making algorithms [21]. Regarding radiometers operating in the infrared and visible ranges, all the devices found in the literature are imaging radiometers, which include the scanning capability to provide a two-dimensional array of pixels for producing an image. They are often called imagers or scanners, and the scanning can be carried out mechanically or electronically by using an array of detectors [86]. From the obtained images, imaging radiometers allow the measurement of the radiance in a specific direction towards the object or under a specific angle of view and, finally, the calculation of the reflectance in order to obtain the suitable spectral indices.
A large number of terms are usually used to describe the same type of radiometer leading to confusion. In this work, the following criterion has been followed: when the device operates in a single frequency band, it is simply called radiometer, for instance, the radiometer on board the soil moisture active passive (SMAP) satellite at 1.41 GHz [23]; when it works in several discontinuous frequency bands of bandwidth larger than approximately 10 nm, multichannel or multispectral radiometer is the term used to define it. Most multispectral radiometers have four basic spectral bands, such as blue, green, red, and NIR bands. Some of them have additional spectral bands in the SWIR and thermal infrared (TIR) regions of the spectrum. Two examples of this kind of scanner are the operational land imager (OLI) and the thermal infrared sensor (TIRS) on board the Landsat-8 satellite [14,26]; finally, when the instrument is designed to obtain imagery over hundreds of narrow and continuous spectral bands with typical bandwidths of 10 nm or less, it is called hyperspectral radiometer or spectrometer, which is described in the following subsection.

Spectrometer
It is a device that detects, measures and analyzes the spectral content of the incident electromagnetic radiation. Conventional imaging spectrometers use gratings or prisms to discriminate the spectral content on the basis of its wavelength [86]. The spectrometers collected in this paper are of two types: imaging spectrometers mounted on manned aircraft, such as airborne visible/infrared imaging spectrometer (AVIRIS) [40] and airborne visible/infrared imaging spectrometer-next generation (AVIRIS-NG) [41], and portable spectrometers.
Imaging spectrometry consists in the acquisition of images in hundreds of contiguous spectral bands such that for each pixel a radiance spectrum or spectral radiance measured in W·sr −1 ·m −2 ·nm −1 can be derived [88]. Imaging spectrometers are designed to measure the spectral radiance coming from an object as a function of two spatial and one spectral dimension, that is, a 3-D dataset usually called data cube. Afterwards, from these radiance values, the reflectance in each pixel is obtained in order to calculate the required spectral indices. Figure 2 shows the data cube provided by a spectrometer under the hyperspectral imaging technique. The third dimension in the cube represents the spectral content of each pixel. Hyperspectral sensors divide the spectrum into several narrow spectral bands, instead of a single wide band. On account of this, these instruments can provide high resolution spectral information from the object. However, a reduction in the width of the spectral bands also means a reduction of the signal received by the sensor. With the objective of maintaining a specific signal-to-noise ratio (SNR), the pixel size has to be enlarged and the spatial resolution decreases, leading to a trade-off between spatial and spectral resolution. In order to improve the spatial resolution of the images, airborne imagery is usually adopted [90]. Remote sensing imaging spectrometers may be categorized by the method they use to achieve the spatial and spectral discrimination. Whiskbroom, pushbroom, framing and windowing are the main techniques to discriminate the spatial information, while filtering, dispersive and interferometric techniques are the most commonly used methods to obtain the spectral discrimination [88,91,92].
Additionally, portable spectrometers that estimate the water content of leaves with destructive and non-destructive techniques are widely used in the reviewed literature because they allow the user to transport the instrument from the laboratory to the sample, instead of taking the sample to the spectrometer. There are different types of portable spectrometers, such as optical or molecular spectrometers, elemental spectroscopy, mass spectrometry (MS), portable gas chromatography-mass spectrometry (GC-MS), ion-mobility spectrometry, NMR spectrometers, small electron spin resonance (ESR) instruments and millimeter wave spectrometers among others [93]. The spectrometers collected in this review are optical covering the visible, NIR and SWIR spectral bands.

Camera
Over the last decade, UAVs, also known as drones, unmanned aircraft systems (UAS) or remotely piloted aircraft systems (RPAS), have been increasingly employed. Recent studies have indicated that the application of these devices is revolutionizing a large number of research areas, such as forestry studies, spatial ecology, ecohydrology and other fields related to environmental monitoring. This revolution has been possible thanks to the technological progress and miniaturization of sensors, airframes, and software. UAV sensing technology has enabled scientists to obtain custom cost-effective imagery at high spatial resolution ranging from 1 cm to 1 m. Before UAVs were consolidated, external companies or institutions were the responsible entities for providing datasets by means of the utilization of sensors mounted on satellites and manned aircraft. Nowadays, research groups purchase or build their own remote sensors without the support of these entities. Although UAV remote sensors are not very different from those on board satellites and manned aircraft, small differences may exist in performance, data processing, calibration, spatial and spectral resolution, measurement geometries and measurement timings [94].
In the heart of this miniaturization, cameras have burst onto the scene as devices combined with UAVs to monitor plant water content. The four types of cameras found in the literature covered by this review are conventional RGB, multispectral, hyperspectral and TIR cameras. A conventional digital camera possesses an image sensor that converts electromagnetic energy into electrical charges. Although the image sensor in most digital cameras is a charge-coupled device (CCD), some of them use complementary metal oxide semiconductor (CMOS) technology. Its spectral range is limited to the three visible bands, that is, red, green and blue bands, and it automatically combines the spectral bands into an image [95]. Multispectral and hyperspectral cameras are the miniaturized devices analogous to the multispectral imaging radiometers and hyperspectral spectrometers respectively, explained in the two previous subsections. Similarly to these devices, the multispectral and hyperspectral cameras of this review operate in the VIS/IR frequency range. Finally, TIR cameras are described in the following subsection.

Thermal Remote Sensors
This type of sensors measure the radiation emitted by an object due to its thermal state, instead of measuring the radiation reflected by the object after being illuminated by solar radiation. All objects with a temperature higher than the absolute zero (−273.15 • C) emit electromagnetic radiation. Due to the infrared atmospheric windows, thermal remote sensors usually operate in the 3-5 µm and 8-14 µm frequency ranges [86]. In the particular case of this review, all the used thermal sensors operate in the second range of 8-14 µm.
The thermal infrared sensors collected in this review are TIR radiometers [9,27], which have already been described in Section 2.1, conventional portable infrared thermometers [27] and TIR cameras. A TIR camera is a device similar to a conventional RGB camera but with a spectral range of approximately 7.5-13.5 µm. Figure 3 shows the simplified crop water stress index (CWSIsi) of a cotton field, which is estimated on the basis of the cotton canopy temperature measured by a TIR camera on board a drone. Its calculation procedure is detailed in Section 3.2.

Radar
It is an electromagnetic sensor, usually in the range of microwave frequencies, used to detect the position, distance and velocity of an object from the radar placement. It consists in radiating electromagnetic energy towards the object and monitoring the reflected signal from it [86]. Depending on the type of signal used to monitor the object, radars can be categorized into two main groups: continuous-wave (CW) and pulsed radars. Specifically, the sensors found in this review are synthetic aperture radars (SAR) mounted on satellites ( Figure 4a) and manned or unmanned aircraft, which use pulsed waveform, and frequencymodulated continuous-wave (FMCW) radars (Figure 4b). (a) Soil moisture content retrieval of the Ugan-Kuqa River Delta Oasis expressed in relative water content (%) by using the Sentinel-1 C-band SAR data [96]; (b) Portable FMCW radar from Anteral S.L. applied to a grapevine trunk [84]. Copyright 2021 IEEE.

Synthetic Aperture Radar
SAR systems, such as the C-band SAR on board the Sentinel-1 satellite, were created in an attempt to improve the spatial resolution at long distances. For a real aperture radar system, the minimum dimension of an object that can be detected (R s ) is calculated as [97]: where λ is the wavelength, h is the height of observation and φ is the diameter of the antenna. Therefore, when the height of observation reaches a considerable value, as in space platforms, achieving a suitable resolution would require using a very large antenna. This limitation can be overcome by using an artificially synthesized virtual antenna that is the fundamental concept of SAR systems. Doppler effect, namely the change in the frequency of the received signal as a result of the relative movement between the sensor and the observed object, is used to synthesize the adequate antenna and obtain high-resolution image data of the surface under test [97]. SAR technology features relevant characteristics compared to other imaging techniques used in remote sensing. Firstly, as happens with other technologies, such as microwave radiometers, radars, scatterometers and thermal sensors among others, it enables the monitoring of the Earth's surface independently of the natural light, that is, at night or on cloudy days. Another important advantage of SAR technology is the penetration capacity that allows the penetration into vegetation to estimate the SMC beneath the canopy at different depths. The specific parameters used to monitor the soil and canopy moisture content are the backscattering coefficients, which depend on specific physical properties of the observed surfaces, such as the dielectric constant (ε), closely related to their amount of water. Each pixel of an image obtained by means of radar imaging represents the backscattering coefficient of that area on the ground that reflects the electromagnetic energy, and the greater the stored value, the more intense the received signals will be [97].

Frequency-Modulated Continuous-Wave Radar
FMCW radars emit continuous radio frequency (RF) power during the measurement, changing their operating frequency by means of a frequency synthesizer that offers several modulation ramps. This type of sensor divides the transmitted power into two signals, the reference signal used as local oscillator (LO) and the RF signal, which is radiated toward the target, reflected, and sent back to the receiver to be mixed with the reference signal. At the mixer input, the time delay between both signals can be estimated from their frequency difference and used to determine the distance D to the reflecting object by applying the following equation [84]: where c is the speed of light in vacuum, ∆ f is the frequency difference measured between both signals, T r is the modulation wave period or ramp time, and BW is the bandwidth of the transmitted signal. For a medium different from vacuum, c must be replaced by the group velocity V g at which the electromagnetic wave propagates in the medium. The depth resolution ∆D, defined as the minimum separation in depth between two targets of equal cross section that makes it possible to identify them as separate targets, can be expressed as follows [84]: The application of this technology on the monitoring of plant water content consists in measuring the power reflected on the plant for different degrees of hydration. The part of the plant more commonly subjected to measurement is its trunk [84].

Scatterometer
It is a side-looking radar sensor designed to measure σ with high radiometric resolution at the cost of sacrificing spatial resolution [98]. Radiometric resolution is the capacity of a detector to detect different input signal levels. The more sensitive the sensor, the higher the radiometric resolution is. Scatterometers are sensors complementary to SAR devices. These sensors, unlike scatterometes, optimize the spatial resolution at the expense of the radiometric resolution. Therefore, while satellite SAR systems are characterized by a spatial resolution of about 1 m, with a low radiometric resolution related to the speckle noise, satellite scatterometers have a spatial resolution of about 10 km, with a radiometric resolution to measure σ within a tenth of a decibel. Although these remote sensors were initially designed to estimate the speed and direction of wind across ocean surfaces, its use has also been extended to monitor continental surfaces. One of the applications of this continental monitoring, which has been included among the different collected references [42], is the vegetation water stress detection. The application of this technology on the monitoring of plant water content is very similar to SAR technology, consisting in the measure of the backscattering coefficients of the observed surfaces.

Terrestrial Laser
The monitorization of CWC by means of laser devices is mainly based on the terrestrial laser scanning (TLS) technology, also referred to as terrestrial light detection and ranging (LiDAR) or topographic LiDAR. This technology measures range (XYZ coordinates) and reflectance of the surface points of the object under test after emitting laser pulses toward them. Modern TLS devices, such as Leica P20, P40 and P50 (Leica Geosystems, Heerbrugg, Switzerland) can measure about 10 6 points per second with an accuracy lower than 1 mm within a range of 10 m (Figure 5a). Lidar instruments can be mounted on a tripod over the ground in order to capture the shape of objects in the surrounding (TLS technology), or placed on moving platforms such as aircraft, helicopters, cars or vessels. If the scanning instrument is placed on board a flying platform, the technology is called airborne Lidar or airborne laser scanning. On the contrary, if the laser is mounted on a car, van or boat, the technology is referred to as mobile laser scanning, terrestrial mobile mapping or, more commonly, mobile mapping [99]. The technology more frequently used in order to estimate the water status of a plant by means of laser devices is TLS. The found references related to this technology are based on the NDII calculation of individual leaves and on the subsequent correlation of these NDII values with EWT, which is the amount of water in leaf per unit leaf surface area. Once this model is created and validated, different trees in their entirety are monitored in order to obtain NDII point clouds and, therefore, using the model, EWT point clouds ( Figure 5b) [37,38].

Soil Moisture Sensors
There are basically two types of soil moisture sensors: those that measure the soil water potential, such as tensiometers and granular matrix sensors, and those based on dielectric techniques that measure VSWC. Tensiometers measure the soil matric potential by means of a porous tip, that is, the suction pressure required for a plant to obtain water from the soil. Granular matrix sensors measure the electrical resistance between two electrodes in a porous matrix, which is a function of the soil matric potential [100,101]. Concerning the dielectric techniques, they can be divided into time-domain and frequency-domain approaches. Time domain reflectometry (TDR) estimates the apparent permittivity of the soil by measuring the propagation time of a pulsed electromagnetic signal along a waveguide consisting of two or three parallel rods. The time the wave takes to reach the end of the waveguide and back is directly related to the dielectric properties of the soil and, therefore, to its soil moisture content [102]. Time domain transmissiometry (TDT) is a variant of TDR that measures the transmission time, instead of the reflection time, of a pulse along a transmission line consisting of a looped rod. As with TDR, there is no complex waveform analysis [103]. In the case of the techniques based on the frequency domain, two types of approaches can be adopted: frequency domain reflectometry (FDR) and capacitance sensors. Although FDR sensors, as in the case of TDR, are constituted by a waveguide composed of two or three parallel rods, a frequency-swept signal is transmitted along the waveguide. Depending on the type of FDR sensor, some establish a relationship between permittivity and resonance frequency and others relate permittivity with reflection and transmission coefficients. Unlike TDR, FDR enables the estimation of the complex dielectric constant. Finally, capacitance techniques estimate the soil dielectric constant by measuring the charge time of a capacitor that uses the soil as the dielectric material between its conductors [102].
The use of soil moisture probes in this review is clearly intended to provide groundtruth data in order to correlate this information with the one provided by the specific remote sensing technique. Only a few VSWC sensors have been used as the main sensing technique to monitor plant water content. This is the case of commercial VSWC sensors inserted into tree trunks ( Figure 6) [76] or a type of FDR sensor specially designed to monitor the xylem water content inside the trunk [71][72][73][74][75]. This last approach, called PDP, is based on open-ended coaxial probes at microwave frequencies that measure the complex reflection coefficient at the end of the semirigid transmission line. The xylem water content is estimated by means of the relationship between the magnitude and phase of Γ and the complex dielectric constant of the medium at the probe tip. Concerning the use of soil moisture sensors as tools to provide ground-truth data, it is worth highlighting the soil moisture sensing controller and optimal estimator (SoilSCAPE) network [104] developed to measure soil moisture at a high temporal resolution by using a combination of wireless sensor networks and soil sensors. SoilSCAPE has been widely used in the calibration/validation of many satellite missions.

Sap-Flow Rate and Sap-Flux Density Gauges
A review of the literature on this kind of sensors has allowed the identification of different methods to carry out sap-flow measurements. These sensors are based on the application of heat in the sapwood, which is transported by the moving sap. Thanks to this movement, sap flow can be estimated by monitoring temperature changes around the heater. The first classification of these measurement devices can be made between those that measure the total sap flow in a plant stem or stem section, which is called sap-flow rate (SFR) (g·h -1 ), and those that measure the amount of sap passing through a specific surface per unit of time, which is called sap-flux density (SFD) (m 3 ·m -2 ·h -1 ) [77].
Among the SFR approaches, it is possible to distinguish the stem heat balance (SHB) and trunk sector heat balance (THB) methods. The SHB method can be used in both woody and herbaceous stems as small as about 4 mm in diameter. Commercial heat balance gauges can fit stems of diameters ranging from 2-125 mm. The SHB technique uses a flexible heater wrapped around the stem and, in turn, enclosed in a layer of cork. These two layers are covered by two others, a layer of foam insulation and an outer weather-resistance shield. By measuring the radial (between the inner and outer cork surfaces) and axial (above and below the heater) temperature gradients by means of thermocouples and applying the heat balance to the stem by means of these temperature gradients, the mass flow rate of sap is calculated in g·h -1 [105]. The THB method is used in tree trunks with diameters larger than 120 mm. As in the previous method, the heat balance is applied to a trunk sector. However, rather than superficially spherical, the heat is applied internally using five stainless steel electrode plates inserted into the wood. Using four thermocouples inserted into the sapwood and positioned at the same height as the plates, and four other thermocouples 100 mm below, it is possible to apply the heat balance in order to estimate the mass flow rate of sap [105].
SFD techniques can be divided into continuous heat SFD methods and heat-pulse SFD methods. Continuous heat approaches relate SFD to a temperature difference measured between a constantly heated needle and different unheated needles intentionally inserted in specific positions around the first one. Thermal dissipation (TD) and heat field deformation (HFD) methods are the most important continuous techniques [77]. In the heat-pulse method, short pulses of heat are applied by means of a needle inserted in the sapwood and the sap-flux density is calculated by measuring the propagation velocity of the heat pulses along the stem. The most relevant heat-pulse techniques are compensation heat-pulse velocity (CHP), T max , heat ratio (HR), calibrated average gradient (CAG) and sapflow+ among others [77]. A limitation of heat-pulse approaches is that they need to be applied to stems with diameters large enough to be able to accommodate the measurement needles. The smallest stem diameters found in the literature range from 5-6 mm [81].

Electromagnetic Spectroscopy
Spectroscopy is the study of the interaction between matter and electromagnetic radiation as a function of the wavelength. Currently, there are different spectroscopy techniques depending on their particular application. The most relevant spectroscopy techniques used to estimate plant water status are: terahertz time-domain spectroscopy (THz-TDS), terahertz continuous-wave (THz-CW) spectroscopy, terahertz quasi time-domain spectroscopy (THz-QTDS), continuous-wave vector network analyzer (CW-VNA) spectroscopy and Fourier transform infrared (FTIR) spectroscopy. Table 1 presents a summary of the different electromagnetic spectroscopy techniques used in this review, including the setup components, if the generated signal is pulsed or continuous and the operating frequency.

Terahertz Time-Domain Spectroscopy
The operating principle of THz-TDS is to generate periodic femtosecond pulses of electromagnetic radiation by means of a photoconductive antenna that acts as the system transmitter and to receive these pulses after passing through the sample by using another photoconductive antenna that acts as the system receiver [106]. The properties of the sample can be estimated by comparing the results between the sample measurement and the reference measurement. The reference measurement, which is carried out for each sample, consists in performing a measurement between the transmitter and receiver without the sample to record the characteristics of the measurement setup. Once the measurement has been completed, it is possible to obtain its frequency dependence by applying the Fourier transform (FT) to the time-domain data [61]. The central component of the THz-TDS is a laser, which emits femtosecond pulses of light. These short pulses are used to excite the photoconductive antennas of the transmitter and receiver. At the transmitter antenna, by means of the application of a static bias voltage to the antenna semiconductor, the photoexcited carriers are accelerated giving rise to short radiated THz pulses. The detection of these pulses is performed by an unbiased photoconductive antenna simultaneously excited by the laser pulses that generate free carriers in the semiconductor material. These carriers are accelerated by the THz incoming electric field generating a photocurrent that feeds a transimpedance amplifier followed by a lock-in amplifier [61,106,107].

Terahertz Continuous-Wave Spectroscopy
THz-CW spectroscopy uses a similar setup to the one of THz-TDS, but the photoconductive antennas are excited by a CW laser beam, instead of light pulses. This beam is the result of photomixing two laser beams slightly detuned against each other. When the photoconductive antenna of the transmitter is excited by this optical signal, a THz radiation at the difference frequency of the two lasers is generated. Similar to THz-TDS, the amplitude of the THz signal after passing through the sample is an indicator of the plant water status. The development of this spectroscopy technique by means of photomixing has the advantage of generating a THz beam with high spectral resolution. This mixing process provides a quasi-monochromatic THz beam, whose bandwidth could be less than 1 kHz. This technique is particularly appropriate when the required frequency band must fit the specific absorption line of the samples under test [61,106].

Terahertz Quasi Time-Domain Spectroscopy
THz-QTDS setup is almost identical to the one of THz-TDS, but replacing the femtosecond laser by a cheap multimode laser diode. Since a laser diode can generate light at different wavelengths simultaneously, it is possible to make use of a large number of difference frequencies by using only one multimode laser diode, unlike THz-CW spectroscopy, which is performed with two CW lasers. The final result is a signal that looks like a train of THz pulses. This similarity with the THz-TDS time-domain waveform is where the name quasi time-domain spectroscopy comes from. The main objective is to achieve a spectroscopy technique similar to THz-TDS, but more compact and cost-effective [61].

Continuous-Wave Vector Network Analyzer Spectroscopy
CW-VNA spectroscopy makes use of a vector network analyzer, apparatus that provides the scattering parameters of the sample under test for a large bandwidth. This spectroscopy technique works at microwave and terahertz frequencies. Generally, a CW signal at microwave frequency is generated by a microwave source included in the VNA. This signal is multiplied by nonlinear diodes (e.g., Schottky diodes) in order to convert it up to higher microwave or terahertz frequencies. These nonlinear diodes used to perform the frequency upconversion are included in modules called VNA extenders. Finally, by means of an antenna, attached at the end of the extender waveguide, the beam is directed toward the antenna of the receiver extender after passing through the sample. At the receiver, a downconversion process reverse to the previous one is carried out. High dynamic range (higher than 100 dB), high spectral resolution (sub-kHz for the best systems) and large bandwidth are the main advantages of this spectroscopy technique. On the contrary, the need for a precise calibration, its upper frequency limit and the possible appearance of the standing wave effect between the antennas are the most relevant drawbacks [106].

Fourier Transform Infrared Spectroscopy
FTIR spectroscopy is based on the Michelson interferometer, consisting in a polychromatic infrared source, a beam splitter, a fixed mirror and a moving mirror ( Figure 7) [106]. The beam splitter is designed to transmit half of the incident light toward the moving mirror and to reflect the other half toward the fixed mirror. The fixed and moving mirrors reflect the radiation back to the beam splitter. Again, half of this reflected radiation is transmitted and half is reflected at the beam splitter. Due to these second reflections, two light beams will appear at the detector path, which are recombined into a single light beam that leaves the interferometer, interacts with the sample and strikes the detector. The optical path difference (OPD) between the two beams due to the displacement of the moving mirror results in interference patterns. For instance, for a specific wavelength, depending on the mirror position, the interference between both beams will be constructive or destructive. When this recombined light beam passes through the sample, a part of the radiation at certain frequencies is absorbed by it and does not reach the detector. Therefore, the final interferogram is composed of all frequencies except for those absorbed by the sample. Finally, in order to convert the OPD-dependent raw data into a frequency-dependent function, the Fourier transform is applied.

Air-Coupled Broadband Ultrasonic Spectroscopy
The development of new broadband ultrasonic transducers, capable of working in gaseous media (air-coupled transducers), has enabled the application of this spectroscopy technique to the measurement of materials without the need to use any coupling fluid to attach the transducer to the surface of the material [55]. The experimental setup ( Figure 8) consists of two specially designed air-coupled piezoelectric transducers operating at MHz frequencies and positioned facing each other at a distance of a few centimeters (approximately 0.3-1.2 MHz and 2 cm for the measurement of LWC). A high voltage (100-400 V) microsecond pulse is applied to the transmitter transducer, which converts it into an ultrasonic pulse that propagates through the air. Passing through the sample, the ultrasonic signal reaches the receiver transducer, which converts it into an electrical one again. Then, it is amplified and filtered (typically a low-pass filter at 10 MHz). Finally, an oscilloscope digitizes this electrical signal, averages several waveforms to reduce the high frequency noise and applies the fast Fourier transform (FFT) transferring the data to a computer. Before performing the measuring process, a measurement without the sample is carried out in order to calibrate the setup. Afterwards, the final measurement is performed by positioning a leaf between the transducers for a few seconds at normal incidence. From this measurement, the magnitude and phase of the transmission coefficient is obtained in the frequency domain. In the particular case of LWC estimation, other parameters, such as the standardized thickness resonance frequency (f/f 0 ), the quality factor of this thickness resonance (Q), the macroscopic effective elastic constant in the leaf thickness direction (c 33 ), the attenuation coefficient (α) and the variation of the attenuation coefficient with the frequency (n), are calculated on the basis of this transmission coefficient and correlated with the RWC. In relation to the f/f 0 parameter, f represents the frequency value associated with the maximum transmittance for any RWC value of a leaf and f 0 is the frequency value at which the maximum transmittance is reached when RWC = 1. Although this spectroscopy technique is a non-destructive and non-invasive approach, all references found in the literature that monitor RWC of leaves by means of this technique make use of leaves detached from the plant instead of in-vivo measurement. For this reason, these references have been classified in this work as destructive approaches.

Global Navigation Satellite System Reflectometry
The most common configuration of a GNSS-R system is the same as that used by a bistatic radar. It consists in measuring the direct and reflected GPS signals by means of right-hand (RHCP) and left-hand (LHCP) circularly polarized antennas respectively. The direct signal is used to calibrate the system and the reflected signal contains information about the Earth's surface. One important application of this technology is the assessment of SMC and VWC from ground-based [29,110], airborne [16] and spaceborne [111] platforms.

Different Approaches Based on Targets
This section presents a brief description of the most relevant approaches on monitoring of plant water content found in the literature. Results are presented according to the part of the plant subjected to measurement, namely, soil, canopy, leaves or trunks. Moreover, approaches applied to leaves are divided into destructive techniques, in which the leaves under test are detached from the plant before measurement, and non-destructive techniques, in which the remote sensing measurements are carried out in-vivo in the plant itself. Concerning the monitoring techniques applied to trunks or stems, they are separated considering invasive approaches that make use of probes, needles and electrodes inserted in the sapwood, and non-invasive approaches, which monitor the plant water status remotely by means of a sensor.

Soil Techniques
The approaches found in the literature that try to estimate SMC [9][10][11][12][13][14][15][16][17] use radiometers [9,11,[13][14][15], radars [10][11][12][13]15], GNSS-R instruments [16] and/or cameras [17] on board satellites, manned and unmanned aircraft, and/or UAVs, commonly known as drones, operating in the range of microwave, infrared and/or visible. Table 2 shows a summary of the references that use soil techniques in this review including the indicator used to estimate the water content, the indicator operating frequency, the fitting curve used between the estimated and measured data, the coefficient of determination (R 2 ), the estimation error, the employed technology, the target or part of the plant under test and the platform used by the sensor to develop the monitoring function.  When radiometers are the remote sensors chosen to monitor SMC, satellites [13][14][15] and manned aircraft [9,11] are the most common platforms. Radiometers on board manned aircraft usually work at microwave frequencies and use the brightness temperature in order to obtain a SMC map of the area under study. Fang et al. [11], in addition to an L-band (1.413 GHz) radiometer mounted on a manned aircraft, use a high-resolution L-band (1.26 GHz) SAR mounted on an unmanned aircraft. By means of a soil moisture disaggregation algorithm that integrates soil moisture retrievals provided by the radiometer and high-resolution observations provided by the radar, it is possible to achieve soil moisture estimates at a finer scale. For a spatial resolution of 800 m, the estimated SMC is well correlated to the in-situ soil moisture measurements with a coefficient of determination (R 2 ) ranging from 0.628 (1.5 Kg/m 2 ≤ VWC ≤ 2.5 Kg/m 2 ) to 0.794 (VWC ≥ 2.5 Kg/m 2 ) accompanied by an unbiased root mean square error (RMSE) ranging from 0.025 to 0.091 m 3 /m 3 . Macelloni et al. [9] make use of two radiometers on board an ultralight manned aircraft, a multichannel microwave radiometer at L (1.4 GHz), C (6.8 GHz) and X (10 GHz) bands, and a thermal infrared radiometer (8-14 µm) to normalize T B by the thermal infrared temperature of the surface. For the case of bare soils, the correlation coefficient (R) of the linear regression model between the estimated and ground measured SMC at L band reaches a value of 0.84 with a standard error of estimate (SEE) of 3%. In addition, an experimental study for estimating the SMC at different depths was carried out by those authors revealing a strong sensitivity of L band to SMC within a layer of 5 cm regardless of the surface roughness, and of X band to SMC within a layer of 2.5 cm on smooth surfaces. The correlation coefficients of the linear regression models between the normalized brightness temperature (T N ) and ground measured SMC at L and X bands are 0.94 and 0.8 respectively. Huang et al. [13] conduct a soil moisture retrieval model combining two satellite sensors, the Sentinel-1 C-band SAR to estimate the SMC and the Landsat-8 VIS/IR OLI to eliminate the influence of the vegetation canopy on the radar backscattering coefficients by means of three vegetation indices: NDVI, VWC and EVI. The highest coefficient of determination between the estimated and in-situ soil moisture is 0.76 (RMSE = 2.04%), which is reached under the VWC vegetation index. Wang et al. [14] provide two combination approaches of the PDI and TVDI indices, the joint model and the combined model. These indices are calculated on the basis of the remote sensing data acquired from the Landsat-8 multichannel radiometers: OLI and TIRS. The RMSE of the joint/combined models between the estimated and ground measured SMC at different soil depths are: 1.49%/1.48% at 0-20 cm, 1.46%/1.48% at 0-30 cm and 1.86%/1.90% at 0-40 cm. Practically for all measurements, RMSE values of the joint and combined models are lower than the ones generated by the individual models. Tao et al. [15] propose an approach similar to the one of [14]. They developed a vegetation backscattering model to retrieve vegetation-covered SMC based on the backscattering coefficient of the RADARSAT-2 C-band SAR and field measurements. In this case, the vegetation scattering is estimated by using multichannel VIS/IR imaging radiometers, particularly four 16-m wide field imagers (WFI) on board the GaoFen-1 satellite that enable the calculation of the NDVI index. The best fit between the measured and estimated SMC features an R 2 of 0.806 with an RMSE of 0.043 m 3 /m 3 .
Concerning the use of radars to monitor SMC, all the references except [10,12] have already been described previously. Eweys et al. [10] use the same SAR previously described in [15] to develop a cumulative distribution function (CDF) matching approach for quantifying the soil and vegetation contributions to the satellite signal. This is achieved by means of the comparison of two signals acquired at the same location when the soil is bare and covered with vegetation on the condition that both signals must have the same cumulative probability. The linear regression model between the measured and retrieved SMC in the presence of vegetation shows an R 2 of 0.64 and an RMSE of 0.03 m 3 /m 3 . Chatterjee at al. [12] make use of the same SAR described in [13] to develop an empirical SMC retrieval model by combining the sensor data with multiple ancillary datasets (e.g., terrain, land cover, soil properties) and correlating the estimated SMC with field measurements performed at the United States Climate Reference Network (USCRN) stations within a layer of 5 cm. In terms of calibration, R 2 = 0.941 (RMSE = 0.029 m 3 /m 3 ) is the highest coefficient of determination achieved by one of the three developed models. In terms of validation, the coefficient of determination decreases to 0.676 with an RMSE of 0.065 m 3 /m 3 .
Regarding the use of GNSS-R instruments to estimate SMC, Munoz-Martin et al. [16] evaluate the reflectivity of the land surface by using an airborne-based microwave interferometer reflectometer (MIR) at L1 (1575.42 MHz) and L5 (1176.45 MHz) bands, and at different integration times. The soil moisture retrieval is obtained by applying an artificial neural network (ANN) to the reflectivity values. The ANN soil moisture data are correlated with a 20-m resolution downscaled estimation approach developed by combining the data of the soil moisture and ocean salinity (SMOS) mission together with the Sentinel-2 NDVI data and the European centre for medium-range forecast (ECMWF) land surface temperature. The maximum correlation coefficient (R = 0.99) is reached for an integration time of 5 s with a SEE value of 0.016 m 3 /m 3 .
In relation to the SMC estimation by using cameras, the collected approaches lack in accuracy. Putra and Nita [17] use aerial photography by means of a conventional digital camera on board a drone to calculate different spectral indices based on the red, green and blue bands. The results show that the correlation coefficient between the ground measured SMC and the RGB bands, the digital elevation model (DEM) and the different spectral indices is below 0.5, where the triangular greenness index (TGI) features the highest value (R = 0.32).

Canopy Techniques
Concerning the approaches aimed at estimating the CWC, in addition to the same sensors used to monitor the SMC (i.e., radiometers [18][19][20][21][22][23][24][25][26][27], radars [28], GNSS-R instruments [29] and/or cameras [26,[30][31][32][33][34][35]), other remote inspection devices have been adopted, such as NIR and SWIR lasers under the TLS technique [36][37][38], portable spectrometers [39] or mounted on manned aircraft [40,41] operating in the visible and infrared spectrum, and a microwave rotating pencil-beam scatterometer installed on the ISS [42]. Table 3 shows a summary of the references that use canopy techniques in this review.  Hunt et al. [18] compare independent estimates of VWC between the Windsat polarimetric microwave multichannel radiometer on board the Coriolis satellite and the moderate resolution imaging spectroradiometer (MODIS), which is a multispectral VIS/IR imaging radiometer aboard the Terra and Aqua NASA's satellites. Results show that there is a linear relationship between both estimates with an R 2 = 0.82 and an RMSE = 0.96 Kg/m 2 . The spectral index used to estimate CWC by using MODIS and, in turn, VWC is NDII. Calvet et al. [19] use the Nimbus-7 scanning multichannel microwave radiometer (SMMR) temperature-corrected data at 6.6 and 37 GHz to derive the water status evolution of the Amazon forest canopy. For this purpose, and after applying atmospheric corrections to the data for water vapor, clouds and rain, the moisture ratio (MR) vegetation humidity index is defined and correlated with variations of stomatal resistance (R 2 = 0.85 with RMSE = 10.91 s·m -1 for the wet season). Jackson et al. [20] map and monitor VWC for corn and soybean canopies by using two spaceborne multispectral VIS/IR imaging radiometers, the thematic mapper (TM) and the enhanced thematic mapper plus (ETM+) on board the Landsat-5 and Landsat-7 satellites respectively. From the obtained images, the NDVI and NDWI indices are calculated and mathematically related to VWC, calculated with gravimetric techniques at the different sampling fields. In relation to the validation of the mathematical models, the minimum RMSEs are obtained with the NDWI index with values for corn and soybean canopies lower than 0.576 and 0.171 Kg/m 2 respectively. In addition, CWC estimation is a crucial factor for early detection of fire risk [21]. Thus, these authors develop a decision-making system in which the moisture of forest canopy is estimated by using the brightness temperature provided by the manned flying platform IL-18 equipped with microwave radiometers operating at 1.43, 13.3 and 37.5 GHz. Firstly, by means of a regional hydrological model, the low moisture areas of the forest canopy are identified.
Secondly, the remote sensing platform is used only in the local fire-dangerous areas identified in the first stage. As a final result, the spatial distribution of the forest fire occurrence probability (FFOP) in Siberia is presented. Chakraborty et al. [22] propose a new vegetation index (VI), the normalized diurnal difference vegetation water content (nddVWC), which is calculated as the difference between VWC at predawn (VWCp) and VWC at afternoon (VWCa) divided by LAI. VWC estimates are obtained by using a radiative transfer model that enables the retrieval of these variables by means of the brightness temperature at 6.9, 10.7 and 18.7 GHz of the advanced microwave scanning radiometer (AMSR-E) aboard the Aqua satellite. The nddVWC index is significantly (p = 0.05) correlated with rainfall, particularly 3 days after the meteorological phenomenon happens. The drought severity index (DSI) is also correlated with nddVWC with a strong negative correlation coefficient of −0.82. Baur et al. [23] combine the brightness temperature of two satellite radiometers, an L-band (1.4 GHz) radiometer on board the soil moisture active passive (SMAP) satellite and the multichannel advanced microwave scanning radiometer 2 (AMSR2) at 6.9 GHz (C band) and 10.7 GHz (X band) onboard the GCOM-W1 satellite. By using these brightness temperatures and by means of the inversion of the tau-omega model, they estimate the vegetation optical depth (VOD) and the single scattering albedo (ω), both indicators of VWC. VOD quantifies the total attenuation within the canopy, which is normalized by the vegetation height derived from LiDAR measurements performed by the geoscience laser altimeter system (GLAS) sensor on board ICESat. This normalized attenuation (k e ), which enables the comparison of areas around the globe with different vegetation height, is the sum of the absorption (k a ) and scattering (k s ) losses, where k s = ωk e . Once these coefficients are calculated, different global maps for time averaged k e , k a and k s are obtained for the L, C and X frequency bands. In addition, seasonal time series of VOD for four selected areas of interest and using a 15-day moving average are drawn. Pan et al. [24] [26], Bhatti et al. conduct a variable rate irrigation (VRI) study for maize and soybean fields in eastern Nebraska by using UAS and satellite imagery, specifically a multispectral VIS/IR camera mounted on a drone, the Landsat-7 ETM+ radiometer and the Landsat-8 OLI and TIRS. From the remote sensing data and incorporating soil water measurements, a spatial evapotranspiration and water balance model is developed to make prescriptions for the VRI treatments. Four types of irrigation treatments are compared: VRI using satellite imagery (VRI-L), VRI using UAS imagery (VRI-U), uniform (U) irrigation, and rainfed (R) irrigation. VRI treatments are compared with U and R treatments according to the crop yield and the mean prescribed seasonal irrigation depth (I p ). For instance, in 2017 the maize yield under VRI-L treatment (12.2 Mg/ha) was significantly larger than the one under R treatment (11.6 Mg/ha), mainly due to the irrigation. In the same year, the I p under VRI-L treatment (76 mm) was significantly greater than the one under U treatment (51 mm). Pereira et al. [27] estimate the water deficit (WD) of a commercial sugarcane area by using the T S -T a and T L -T a temperature differences, where T S is the land surface temperature, T L is the leaf temperature and T a is the air temperature. T L is measured by a portable infrared thermometer and T S is obtained by using a spectral image from Landsat-8 TIRS. The T L -T a difference, which is found to be positive for WD periods and negative for water surplus (WS) periods, correlates better with WD than the T S -T a difference.
Only one collected study report the use of radars to monitor CWC or VWC [28]. In this reference, Ma et al. assess the ability of the radar included in the NASA's passive and active L-and S-band sensor (PALS) to estimate VWC of corn and soybean fields by means of the backscattering coefficients (σ hh , σ vv , σ vh and σ hv ), polarization ratios, and radar vegetation index (RVI). After correlating these parameters with field VWC measurements, L-band σ hh /σ vv is found to be the most correlated for corn with R = 0.81 (RMSE = 0.53 Kg/m 2 ), while for soybean the L-band σ hh /σ hv parameter reaches the highest correlation with R = 0.9 (RMSE = 0.12 Kg/m 2 ).
Regarding the use of GNSS-R instruments to estimate CWC, Camps et al. [29] analyze the influence of the vegetation on the GPS L1 coarse/acquisition (C/A) code signal by installing a dual-input GPS receiver with a dual-polarization antenna under a beech forest for more than a year. The carrier-to-noise (C/N 0 ) data are compared with different groundtruth datasets, observing the highest correlation for the NDVI indicator with a maximum R 2 value of 0.92 (RMSE = 0.49 dB) for the right-hand circular polarization.
When the sensor chosen to monitor CWC or VWC is a camera, in all cases except in [33] the sensor is mounted on an unmanned aerial vehicle, commonly referred to as drone. Caruso et al. [30] make use of two cameras, a conventional RGB camera and a multispectral NIR-RG camera, to estimate LAI, tree height, canopy diameter and canopy volume of olive trees that are either irrigated or rainfed. LAI is estimated by means of the NDVI index, while the rest of the parameters are obtained from the digital surface model (DSM) by using automatic aerial triangulation. The coefficient of determination of the linear correlation between LAI and NDVI is 0.78. In addition, the water stress integral (WSI) calculated in the field from the SWP is linearly correlated with the monthly canopy volume increment yielding an R 2 value of 0.99. Chen et al. [31] estimate the water content of different parts of the cotton plant by using a VIS/IR multispectral camera, by means of which they calculate several spectral indices normalized by the canopy temperature (T c ). This type of vegetation index, which is called vegetation water supply index (VSWI = VI/T c ), enables an increase in the correlation coefficient in comparison with the single vegetation index. But the highest correlation and accuracy in predicting cotton plant water content are achieved by means of multivariate linear regression models, reaching an R 2 value of 0.998 for the model calibration between LWC, measured by gravimetric techniques, and several VSWIs. The RMSE for the validation of this model reaches a value of 0.023. Bian et al. [32] make use of the same multispectral camera as the one of [31] together with a TIR camera in order to simplify the calculation of the crop water stress index (CWSI), giving rise to the CWSIsi indicator. CWSIsi is calculated by means of the TIR camera, while the multispectral sensor is intended to calculate the NDVI, transformed chlorophyll absorption in reflectance index (TCARI) and optimized soil-adjusted vegetation index (OSAVI) spectral indices. CWSIsi is calculated as the ratio between T c -T wet and T dry -T wet , where T c is the average canopy temperature, T wet is the lower boundary temperature of the canopy and T dry is the upper boundary temperature of the canopy. T c , T wet and T dry are obtained from the canopy temperature histograms, which, in turn, are obtained from the UAV thermal infrared images (Figure 3). Particularly for this work, T wet is the mean of the lowest 0.5% of canopy temperatures and T dry is the mean of the highest 0.5% of canopy temperatures. In comparison with other calculation methods of CWSI and the three previous spectral indices, CWSIsi features higher correlation with the stomatal conductance (R 2 = 0.66), transpiration rate (R 2 = 0.592) and soil volumetric water content at the depth of 45 cm (R = 0.812). Huo et al. [33] retrieve the surface temperature and emissivity (ε λ ) of twelve potted wheat plants with different water contents by means of a hyperspectral TIR camera mounted on a tripod in order to detect water deficits and their health condition. Comparing the temperature of the different plants, changes in the canopy structure are observed and the health of each plant may be inferred. An increase in emissivity is noticed for the first three measurements as the plants' water content decreases. But amazingly, between the third and fourth measurements, the period of highest drought in the experiment, the emissivity decreases. This trend is related to changes in the plants' internal structure, physical and biochemical properties, and normal photosynthesis processes. Zhang et al. [34] combine two different cameras, a high-resolution RGB camera and a TIR camera, mounted on two different drones to monitor water stress in maize crops with different irrigation treatments. By using both cameras, a T c map is extracted and, from this temperature map, three water stress indicators are obtained, such as the degrees above non-stress (DANS), the standard deviation of canopy temperature (CTSD) and the canopy temperature coefficient of variation (CTCV). These three indices together with T c are correlated with field measurements of ground-truth T c , stomatal conductance, LAI and SWC. The linear regression model between the estimated and ground measured T c features a R 2 value of 0.94 (RMSE = 0.7 • C). In relation to the linear regression models between the previous water stress indicators and the field measurements, the highest R 2 is obtained between the DANS indicator and LAI with a value of 0.77 (p < 0.01), which is the same as the one between T c and LAI. Finally, by means of a TIR camera, Quebrajo et al. [35] calculate the temperature difference between the vegetation cover of sugar beet plants in a plot with large spatial variability in terms of soil properties and the prevailing air temperature at the time of flight. This temperature difference is used to obtain the WSI and CWSI water stress indices. The highest linear correlation (R 2 = 0.94) is reached between WSI and plant sugar content evaluated experimentally by means of the manual harvesting of six samples per plot.
Concerning the following sensor to monitor CWC, three references [36][37][38] use terrestrial LiDAR to estimate LWC vertical profiles of different tree species. Zhu et al. [36] estimate the LWC of the tree canopy by calculating the backscatter coefficient of the target, defined as the backscatter cross section of the target normalized by the illuminated target area, which is correlated with measured LWC reaching a maximum R 2 value of 0.66 (RMSE = 0.004 g/cm 2 and p < 0.0001) for five plant species. Once the model is created, its validation provides a maximum R 2 value of 0.85 (RMSE = 0.0016 g/cm 2 ) correlating the estimated and measured LWC of individual leaves. Then, this model is applied to the whole tree in order to obtain a LWC vertical profile. They employ an SWIR (1550 nm) laser. Elsherif et al. [37,38] make use of a similar measurement procedure, but they estimate the EWT by using the NDII index by means of two lasers at NIR (808 nm) and SWIR (1550 nm) wavelengths. In [37], the highest correlation between EWT and NDII is obtained in the October datasheet for four different tree species reaching a maximum R 2 value of 0.93 (p < 0.05). The RMSE and RRMSE for the validation of this model reach values of 0.00065 g/cm 2 and 4.5% respectively. Similarly in [38], the highest R 2 for four different tree species corresponds to a value of 0.94 with a relative error for the model validation of 6.3%.
Hyperspectral techniques, commonly known as spectrometers, are also used to monitor CWC. Among the different types of spectrometers employed for this purpose, portable spectrometers and imaging spectrometers mounted on manned aircraft are the most relevant. Kycko et al. [39], by means of the calculation of hyperspectral vegetation indices by applying a field-portable spectrometer (350-2500 nm) to the canopy, assess the health condition of alpine sward dominant species (Agrostis rupestris, Festuca picta, and Luzula alpino-pilosa) of the Tatra National Park (TPN), which are strongly influenced by tourists. All the polygons under test are divided into: areas subjected to trampling, reference areas not subjected to trampling and recultivated areas with limited access and protected by biodegradable jute mats. Each species is compared in each of the three areas correlating the different spectral indices with field measurement. In the specific case of CWC estimation, the water band index (WBI) and NDWI indicators correlated with T L -T a are used to assess the plant water stress. The highest correlation coefficients are obtained for Luzula alpinopilosa with values of 0.88 and 0.92 for WBI and NDWI respectively. Serrano et al. [40] assess the ability of AVIRIS to estimate canopy RWC of three different chaparral communities. For this purpose, several vegetation indices are calculated from the surface reflectance values obtained from the AVIRIS images at different wavelengths. Water index (WI) and NDWI are the most correlated indicators with canopy RWC, which is measured by gravimetric techniques applied to the leaves of the different species. The highest R 2 for WI and NDWI is obtained when the ceanothus and chamise chaparral species represent more than 70% of the total vegetation cover, reaching a value of 0.88 and 0.86 respectively (p < 0.001). Miller et al. [41], in addition to AVIRIS, use the AVIRIS-NG sensor mounted on a manned aircraft to quantify the response of common urban tree species and turfgrass to the 2012-2016 California drought by means of multiple endmember spectral mixture analysis (MESMA), vegetation indices and continuum removed absorption features. Although they do not include correlations between vegetation indices and field measurements, variation analyses and correlations among the own indices are performed. All tree species and turfgrass feature four or more spectral indices with lower mean values during the drought. Most spectral indicators are correlated with one another with R values larger than 0.4 in the majority of cases.
Concerning the last sensor to monitor CWC, van Emmerik et al. [42] use a rotating pencil-beam scatterometer installed on the ISS to observe changes in the radar backscatter over the Amazon and relate these changes with tree water deficit measured with dendrometers. Predawn backscatter is the most sensitive to vegetation water stress. Both VV and HH backscatter decrease between 0.5 and 1 dB when the tree water deficit increases a maximum value of 600 mm (for the Dipteryx odorata species).
In the three THz-TDS approaches found in the literature, the Landau-Lifshitz-Looyenga model [112] is employed to estimate the LWC. This model relates the dielectric constant of a leaf to the dielectric constants and volumetric concentrations of its components, that is, water, dry solid tissue and air. In addition, Gente et al. [43] implement a numerical algorithm in order to iteratively minimize the difference between the transmission coefficient measured by THz-TDS and the one obtained by using the previous model. When the simulation of the algorithm converges, the LWC is obtained. This process is applied to different barley leaves with different degrees of hydration and the results are compared with gravimetric techniques obtaining a RMSE value of approximately 2%. Jördens et al. [44], after calculating the volumetric concentrations of the leaf components by a destructive technique, apply the model in order to estimate the dielectric constant and, in turn, the refractive index (RI) and the absorption coefficient of coffee leaves with different degrees of hydration. For this purpose, the influence of the scattering due to the roughness of the leaf surface is considered by means of a Rayleigh roughness factor. Although no correlation analyses are presented, a good agreement between measurement (0.9 THz) and simulation is obtained. Furthermore, a long-term study was performed by subjecting a Coffea arabica plant to water stress for 21 days. By applying the developed model, at the beginning of the study the estimated water content of the leaf under test is about 70%. At the end of the experiment, the water content is reduced to 15%. Finally, Singh et al. [45], by following an iterative algorithm similar to the one of [43], develop the first three-dimensional (3D) water mapping of an agave leaf. In order to achieve this 3D distribution, each leaf is cut transversely to produce thin slices (~600 µm) that are imaged using THz-TDS in transmission mode pixel by pixel. In addition to this, by averaging the water content of several adjacent slices, they obtain 12 water content values along the leaf that are subsequently compared with gravimetric methods. However, no correlation analyses between them are presented.  The following destructive technique applied to leaves is the CW spectroscopy. Chuah et al. [46], by means of a Wiltron scalar network analyzer and a slotted waveguide, measure the dielectric constant of rubber and oil palm leaves at X band by following the waveguide thin sheet technique. Measurements of the leaves' dielectric constant are compared with their gravimetric moisture content and with dielectric constant values obtained from two theoretical models. Although no correlation analyses are performed, the theoretical model from Ulaby and El-Rayes [113] is found to be able to give good estimates of the dielectric constant for the two types of leaf samples. Zahid et al. [47] measure the transmission scattering parameters of coffee arabica, pea shoot and baby spinach leaves in the frequency range of 0.75 to 1.1 THz by using a VNA and two frequency extenders in order to monitor their water content for four consecutive days. On each day, 10 rounds of weight measurements are performed in order to calculate the gravimetric water content of the leaves. From the measurements obtained from the VNA in different domains, 25 significant features are considered to develop three machine learning algorithms, such as support vector machine (SVM), K-nearest neighbor (KNN) and decision-tree (D-Tree). The most accurate results are achieved for the SVM algorithm with an overall accuracy of 98.8%, 97.15% and 96.82% for coffee arabica, pea shoot and baby spinach leaves respectively. Finally, Sancho-Knapik et al. [48] compare two RWC measurement techniques of poplar leaves: reflectivity at 1730 MHz by means of a VNA and a microwave digital cordless telephony (DCT) patch antenna, and ratio of reflectances between 1300 and 1450 nm (R 1300 /R 1450 index) measured by a portable NIR spectrometer. The correlation between these measurements and gravimetric methods provides R 2 values of 0.94 (SEE = 0.0219 and p < 0.0001) and 0.98 (SEE = 0.00394 and p < 0.0001) for the NIR and microwave techniques respectively.
In relation to the references that use FTIR spectroscopy, only the one performed by Fabre et al. [49] has been collected. This study is the first one to publish the relationship between LWC and reflectance in the 3-15 µm range. For this purpose, they use sorghum, sunflower and cherry tree leaves with different degrees of hydration. In the 3-5.5 µm range, the reflectance for all samples decreases exponentially when the LWC increases. In the 8-15 µm range and with LWC lower than 30%, all samples except for the adaxial side cherry tree leaf (almost unaffected by drying) follow a similar exponential trend but with much lower variation levels. At higher LWC values, reflectances remain unchanged.
Concerning portable spectrometers, one of the two collected references has already been described previously [48]. In the second reference, Arevalo-Ramirez et al. [50] assess the feasibility of estimating the water content of Pinus radiata and Eucalyptus globulus leaves at different dehydration stages without prior knowledge of their reflectance in the SWIR region (1000-2500 nm), which is not available in most vision-based sensors. To this end, 70% of the reflectance values measured in 90 leaves per species in the 350-2500 nm range is used to train two machine-learning regressors, Kernel Ridge Regressor (KRR) and Random Forest Regressor (RFR). The implemented algorithms based on these regressors enable the reconstruction of the leaf reflectance in the SWIR region with a minimum RMSE lower than 3%, value obtained from the KRR-based algorithm validation with the remaining 30% of the data. Finally, by using the reflectance reconstructed by the KRR-based algorithm, 16 spectral indices are calculated with a relative percentage error lower than 10% for seven of them.
Finally, the last technique found in the literature to monitor LWC is the air-coupled broadband ultrasonic spectroscopy. Six references that use this approach have been collected. The performed measurement process is the same for all references and is thoroughly described in Section 2.11. For this reason, only the correlation results of each work between the gravimetric RWC and the corresponding parameter will be presented. Sancho-Knapik et al. [51], by using piezoelectric transducers in the 0.3-1.2 MHz range, correlate the RWC of Populus x euramericana and Prunus laurocerasus leaves with f/f 0 , adjusting this relationship to a cubic function with an R 2 value of 0.99 (p < 0.0001) for both species. Sancho-Knapik et al. [52] correlate the RWC of Prunus laurocerasus, Ligustrum lucidum, Populus x euroamericana and Platanus hispanica leaves with ∆f/f 0 and c 33 [54] correlate the RWC of Quercus muehlenbergii Engelm leaves with the percentage of loss of thickness (PLT), the leaf density (ρ), f/f 0 and 1/Q with R 2 values of 0.93 (linear fitting curve with p < 0.0001), 0.88 (second-order fitting curve with p < 0.002), 0.99 (third-order fitting curve with p < 0.0001) and 0.97 (linear fitting curve with p < 0.0001) respectively. Sancho-Knapik et al. [55] compare the air-coupled ultrasonic spectroscopy in the 0.3-1.2 MHz range with a CW microwave (2.25-2.55 GHz) technique, similar to the one presented in [48], to estimate the LWC of grapevine leaves. For the microwave technique and from the reflectance values, they define the existing area (S) indicator. A high linear correlation is obtained between RWC and S with a R 2 value of 0.995, (p < 0.0001). Concerning the ultrasonic technique, the highest R 2 is reached between RWC and f/f 0 with a value of 0.99 (third-order fitting curve with p < 0.0001). Finally, Fariñas et al. [56], after performing ultrasonic measurements in the 0.15-1.6 MHz range to obtain the transmission coefficient of Viburnum tinus leaves, apply two different machine-learning algorithms to estimate RWC, convolutional neural networks (CNN) and random forest (RF). In the evaluation process, both methods are tested successfully with R values of 0.92 (RMSE = 0.04) and 0.84 (RMSE = 0.06) respectively.
Concerning the THz-TDS approaches, Li et al. [57] use the partial least squares (PLS) water-content prediction model to estimate the LWC of soybean leaves from THz-TDS measurements. The correlation coefficient between the water content predicted by the PLS model and the measured water content is 0.9153 (RMSE = 0.0526%). In order to validate this estimation method, three long-term experiments are performed to estimate: changes in the water content of soybean plants subjected to a gradient of soil water availability, changes in the leaf moisture of soybean plants growing in different substrates and changes in the stomatal conductance of soybean leaves before and after abscisic acid treatment. Pagano et al. [58] use a THz generation technology similar to THz-TDS, but replacing the femtosecond laser and the photoconductive antennas with a THz quantum cascade laser (QCL). By applying the THz-QCL technique in transmission at 2.55 THz, they estimate the leaf water mass of six species by means of the product of the leaf optical depth (τ) and the leaf surface area (LA). Leaf optical depth is obtained from the THz transmittance measurements and LA by using an RGB digital camera and an image processing software. The highest R 2 between τ·LA and the gravimetric leaf water mass (M w ) is obtained for the Corylus avellana tree leaves with a value of 0.99. Baldacci et al. [59] apply the same measurement technique as the one of [58] to Vitis vinifera leaves at the same THz frequency. The linear regression model between τ·LA and M w features a R 2 value of 0.95. Castro-Camus et al. [60] use THz-TDS and the Landau-Lifshitz-Looyenga model employed by [43][44][45] to monitor the water loss of Arabidopsis thaliana cauline leaves under water restriction in two substrates with different water retention capacity. No correlation analyses are presented. Gente et al. [61] present a review on monitoring of leaf water content with THz and sub-THz waves, in which most of the collected references have already been considered in this work. Finally, Gente and Koch [62] use THz-TDS in the 150-300 GHz range with an innovative setup that enables the automatic measurement of silver fir (Abies alba) needles during several weeks. Several silver fir seedlings under different drought stress treatments are arranged in a circular arrangement and one needle of these small trees is placed periodically in the same position between the transmitter and the receiver by using a small and weak magnet. After THz measurements, the gravimetric RWC of each seedling is obtained. A significant R value of −0.98 (p < 0.001) is obtained between THz transmission expressed in percentage and gravimetric RWC for all measured seedlings.
Gente et al. [63] propose the first outdoor measurement system for continuous invivo monitoring of LWC by using THz-QTDS in transmission at 100-200 GHz. For this purpose, a corn leaf is inserted into the portable spectrometer and measured every 10 min during two days and nights. The effective water thickness of the leaf is calculated from the THz transmission values and represented as a function of time. Although no correlation analyses are presented, higher transmission values correspond to lower water content in the leaf and vice versa.
In relation to CW spectroscopy, three references have been collected. Van Emmerik et al. [64] use a microstrip line resonator to measure in-vivo corn leaves through a sampling window (9 mm × 9 mm) cut in the ground plane. Teflon blocks are placed under the sensor and on top of the leaf to have the sample inside a material of known dielectric constant. This resonator is connected to a portable cable and antenna analyzer that can works as VNA. For each leaf, the S 11 parameter, which depends on the water content of the sample, is measured at 2.1-4.1 GHz. Particularly, the resonant frequency, at which the magnitude of S 11 is minimum, depends on the dielectric constant of the leaf and, therefore, on its water content. Although no correlation analyses are presented, they show the relationship between the resonance frequency variation (∆f r ) and the gravimetric RWC for individual leaves with different degrees of hydration. Observing this relationship, it is possible to state that as the corn plant dries down and the RWC of the leaf under test decreases, ∆f r increases following a quadratic trend. The resonance measurements also show a dynamic vertical profile in the dielectric constant of the leaves, depending on their relative height in the plant. Afsharinejad et al. [65] apply THz-CW spectroscopy at 0.1-1.1 THz to measure the transmission coefficient in live leaves of seven species (cabbage, mint, radish, strawberry, coffee, rhubarb, and lettuce) to study the range of variation in THz transmission among leaves of the same plant, depending on their position, height on a stem and age, among others, and among leaves of the same species, but belonging to different individual plants.
Given the large variability in the transmission coefficient of leaves on the same plant and among plants of the same species, it is necessary to collect a sufficiently large number of leaves in each plant to make accurate estimates of THz transmission possible. Considering all the plants of this study, the standard deviation (SD) of the leaf transmission at 0.4 THz is >0.3. The following reference, written by Gente et al. [66], could also be classified as a canopy technique. However, due to the employed technology, frequency range and plant size, it has been categorized as a leaf technique. They implement a setup, which consists of a 35 GHz Gunn oscillator, a detector and two HDPE lenses, to measure in transmission the amount of radiation absorbed by 22 potted plants with different irrigation treatments. For this purpose, the transmitter is kept in a fixed position and the detector can be moved an angular range of 270 • behind the plant by using a motorized arm. From the measurements with and without the plant, it is possible to calculate the amount of radiation absorbed by the plant. A linear regression model between the plant absorbance and the plant water content measured by gravimetric techniques is determined with a R 2 value of 0.83 (RMSE = 0.17 g).
In relation to the use of portable spectrometers to monitor LWC non-destructively, only one reference, the one written by Ribera-Fonseca [67], has been collected. From the reflectance value measured by a portable VIS/NIR spectrometer on 24 leaves of 12 blueberry shrubs subjected to two different irrigation treatments (full irrigation and mild water deficit), the NDVI and WI spectral indices are calculated and correlated with field measurements, such as the midday leaf water potential (Ψ MD ) and the stomatal conductance. The relationship between Ψ MD and both indicators is described by a linear model with R 2 values (p < 0.01) of 0.61 and 0.60 for NDVI and WI respectively. In the case of the stomatal conductance, the relationship is exponential with R 2 values (p < 0.01) of 0.58 and 0.69 for NDVI and WI respectively. This study reveals that the NDVI and WI indices can be used as non-destructive predictors of water stress in blueberry shrubs, but further studies should be performed to confirm these results.
Although, it is conducted on fruits, it is worth mentioning the study of Litvischenko et al. [68] who designed a device similar to a portable FMCW radar to measure the reflection coefficient of sunflowers in the 25.86-37.51 GHz frequency range in order to determine the ripeness state of sunflower seeds in the field, which depends on their water content. For this purpose, the ratio between the reflection coefficient of the sunflower under test and the reflection coefficient of ripe sunflower seeds with a moisture content of 7% is calculated. Simultaneously, the moisture content of the seeds is measured with a capacitive hydrometer. The relationship between both magnitudes can be adjusted by a quadratic polynomial with an R 2 value of 0.99 and an RMSE for the model calibration of 0.27% in the volumetric moisture content estimation of the seeds.
The last type of sensor, usually called porometer, is based on the measurement of g s (mmol·m -2 ·s -1 ), an indicator of the plant's water status and, in the end, of its drought stress level. Vilagrosa et al. [69] analyze the response of Pistacia lentiscus and Quercus coccifera seedlings to an intense drought period in terms of g s , vulnerability to cavitation, leaf dieback and survival. They correlate the maximal g s registered throughout the day (g s-max ) with the predawn water potential (Ψ PD ) measured by means of a pressure chamber.
For the Q. coccifera species, the found fitting curve is exponential with an R 2 value of 0.91, while for the P. lentiscus species the estimated Ψ PD follows a 1/g s-max trend with an R 2 value of 0.98. Medrano et al. [70] present a review on how certain species, grapevines among others, respond to progressive drought. For this purpose, the different studies relate photosynthetic parameters of the plant to the maximum light-saturated g s . Particularly, the correlation with Ψ PD is presented, obtaining a quadratic regression model with an R 2 value of 0.97 (p < 0.001) for potted grapevines.

Techniques Applied to the Trunk or Stems
Although all the collected techniques that are applied to the trunk or stems of the plant are non-destructive, some of them use invasive methods [71][72][73][74]76,81] to monitor its water status, which consist in the insertion of needles or probes inside the trunk. Concerning the non-invasive methods, half of them need contact with the trunk or stems to be able to develop the monitoring function [75,[78][79][80].
Among the main sensors that estimate tree water status by invasive techniques we can find coaxial PDPs [71][72][73][74], VSWC sensors inserted into tree trunks [76] and sap-flow measurements based on sensors that measure SFD [81]. Table 5 shows a summary of the references that use trunk techniques specifying if the used technique is invasive or non-invasive.  The measurement procedure of the PDP invasive techniques [71][72][73][74] is similar in all of them. It consists in inserting an open-ended coaxial probe in the xylem tissue of the trunk in order to estimate the real and imaginary parts of its dielectric constant from the direct measurement of the xylem reflection coefficient. The value of this dielectric constant is an indicator of the plant water status. Although no correlation analyses are presented, different magnitudes simultaneously measured in each study are represented in order to find a relationship among their trends. Xylem dielectric constant of several tree species at P (0.5 GHz), L (1.25 GHz) and C (5 GHz) bands, XWP, vapor pressure deficit (VPD) and xylem SFD are the main concurrent magnitudes analyzed in most PDP studies. In short, xylem SFD is well correlated with VPD, VPD is negatively correlated with XDC with apparent decorrelation at low VPDs (e.g., after rain or during cloudy days), XDC responds directly to changes in XWP and temporal variations in XDC correspond to changes in xylem SFD (water content), sap chemistry and VPD. Stott et al. [76] insert four types of commercial VSWC sensors into the sapwood of peach trees in order to estimate their dielectric constant at 70 MHz, 78 MHz, 3.5 GHz and 4 GHz, depending on the used sensor, and correlate it with the SWP. Results improve when the sensor waveguides are in contact with a greater percentage of sapwood. An increment of 0.5 MPa in SWP due to irrigation is associated with a maximum increment of 0.39 in dielectric constant, obtained for the 78 MHz sensor. A maximum R 2 value of 0.29 between both magnitudes is reached for the GHz sensors. Finally, Miner et al. [81] design and implement a device based on a new HP technique, the T m Ratio method, in order to measure the SFD of potted corn and sunflower plants grown in a greenhouse. Instead of presenting SFD data, the product between the apparent heat pulse velocity (V) and the stem area (A) is obtained. Results of this product are compared to the HP T max method and correlated with the gravimetric transpiration measured by digital balances. While for sunflower plants R 2  Finally, among the main non-invasive techniques applied to the plant trunk or stems we can include coaxial PDPs [75], sap-flow measurements based on sensors that measure SFR [78][79][80], CW-VNA spectroscopy [82,83], portable FMCW radars [84] and portable NMR sensors [85].
Franchois et al. [75], by using a non-invasive contact PDP, estimate the dielectric constant of branches, parts of trunk and needles of two conifers (Abies Normanniana and Picea Omorika) at different tree heights in the 1-10 GHz frequency range. Results confirm that wood possesses an anisotropic nature, the longitudinal complex dielectric constant is approximately 1.5-3 times higher than the transverse component and the phloem layer and the new needles have the highest complex dielectric constant.
Sakuratani [78] use SHB to estimate SFR of potted soybean and sunflower plants. The estimated SFR is compared to the transpiration rate determined by two methods: digital balances and chamber method by using ventilated platinum resistance psychrometers. Although no correlation analyses are presented, there is a good agreement between both magnitudes for the two methods. Dugas [79] measure the xylem mass flow rate of cotton plant by using a SHB sensor in field experiments. In order to evaluate the accuracy of the used method, it is previously applied to potted cotton plants of a glasshouse, relating SFR with the plant transpiration measured by successive mass measurements. The RMSE between both magnitudes reaches a value of 24.6 g·plant -1 ·d -1 , which is 8.6% of the average value of the measured transpiration values. Repeating the process in the field and evaluating the transpiration by means of a weighing lysimeter, a RMSE value of 69 g·plant -1 ·d −1 is obtained. Finally, Nakano and Iwasaki [80] measure the SFR in xylem and phloem separately by applying a combination of the SHB method and girdling to Meiwa kumquat trees. After studying the influence of the gauge temperature on the measurement accuracy, they keep this temperature constant to estimate that for a nonbearing branch the phloem sap flows basipetally at 0.27 g·h -1 . For a bearing branch, it is estimated that the phloem sap flows acropetally at 0.30 g·h -1 .
Concerning the CW-VNA spectroscopy, Torres et al. [82] and Santesteban et al. [83] implement the first monitoring system of plant water status by using a non-destructive non-invasive contactless technique applied to the plant trunk. In both references, a VNA, two plane-convex lenses and a frequency extender are used to measure the reflection (S 11 parameter) of a grapevine trunk in the time domain at 0.14-0.22 THz. For this purpose, the potted plant is settled inside a grown chamber in order to control the external conditions, mainly temperature and lighting. In both studies, two similar experiments are performed: irrigation cycles and day/night cycles. In these measurements, the S 11 module of a chosen trunk point is recorded every 5-10 min, identifying the value corresponding to the trunk reflection in each recording. Although no correlation analyses are presented, concurrent magnitudes, such as the trunk thickness displacement (dendrometer) and the volumetric water content of the soil (VSWC sensor), are measured and compared with excellent agreement. In addition, Santesteban et al. [83] perform a third experiment based on the girdling of the phloem, by means of which the S 11 module changes totally during the day and night periods (inverted response).
Quemada et al. [84] use a technique similar to the one of [82,83], but replacing the VNA with a portable low-power cost-effective FMCW radar from Anteral S.L. Instead of recording the S 11 module for each measurement, the fast Fourier transform (FFT) of the signal I-j Q is calculated, where I and Q are the outputs of the radar's IQ mixer. As in the two previous studies, the potted grapevine is settled inside a grown chamber to control the external conditions and a VSWC sensor is used to provide the truth data. With the objective of assessing the ability of this technique to monitor plant water content, three different experiments are performed: irrigation cycles, day/night cycles and comparison between irrigation cycles at different temperatures. For the first two experiments, regression models between the reflected signal and the soil volumetric water content are developed. For the first drying period after the first irrigation, both magnitudes are adjusted by a quadratic polynomial with an R 2 value of 0.995. For the day/night cycles, choosing two daytime and two nighttime periods between the 40th and 44th six-hour cycles, the regression model is linear with a R 2 value of 0.983.
Finally, Windt and Blümler [85] use portable NMR sensors to estimate xylem sap flow in a poplar tree (Populus nigra) and changes in the trunk absolute water content of an oak tree (Quercus robur) on the basis of the NMR signal (A NMR ). These authors compare the results obtained from the NMR sensor with variations in the trunk diameter measured with two dendrometers positioned above and below the NMR sensor, assessing the same pattern for both parameters during a 6-day period.

Comparison between Different Remote Sensing Approaches According Their Characteristics
The scope of this paper covers a broad range of 77 references published between 1981 and 2021 and collected by using different search web sites, especially Scopus. Among all of them, 74 references are research papers and the remaining three references are review papers [61,70,77]. Because each of these review papers includes several studies and some of them have already been considered in the results of this work, they have been removed from this section.
From the collected results, it is possible to assert that the use of remote sensing to monitor plant water content is well established around different countries over the last decades. This expansion has been possible thanks to the development of a wide variety of cutting-edge sensors and platforms. According to the year of publication, the number of published papers has increased exponentially over time, as shown in Figure 9, with two papers published in the decade 1981-1990 (2.7%), seven in the decade 1991-2000 (9.5%), seven in the decade 2001-2010 (9.5%) and 57 in the decade 2011-2020 (77%). Technological development over the last decades has enabled the improvement and irruption of new techniques to monitor plant water content. Together with the year of publication, Figure 9 includes the technologies employed to develop the monitoring function. The most used sensor is the radiometer, which is employed in 15 papers (20.3%), followed by CW spectroscopy (12.2%), camera (10.8%) and THz-TDS (10.8%).
In relation to the references that use techniques applied to leaves, THz-TDS, CW and ultrasonic spectroscopies are the most common technologies. Reflection, transmission and absorption coefficients are the most typical estimators used to develop the monitoring function, reaching high R 2 values (0.995) [55]. Finally, among the references that make use of techniques applied to the plant trunk or stems, PDPs are the most frequent sensors. It is worth emphasizing that very few R 2 values are reported in the collected literature. The highest R 2 value (0.995) is achieved by using a portable FMCW radar at 21.4-24.8 GHz [84]. Figure 11 shows the number of publications in which each water estimator is used, considering only the estimators that have been used two or more times. The NDVI vegetation index is the most common estimator, being used in nine studies (12.2%). Backscattering coefficients are the second most used indicator with eight uses (10.8%) followed by the spectral reflectance, reflection coefficient and dielectric constant, employed in 6 publications each (8.1%). Figure 11. Number of publications in which each water estimator is used, considering only the estimators that have been used two or more times. Note that each publication can use more than one water indicator. NDVI = Normalized difference vegetation index, σ = Backscattering coefficients, R λ = Spectral reflectance, Г = Reflection coefficient, ε = Dielectric constant, NDII = Normalized difference infrared index, NDWI = Normalized difference water index, MSI = Moisture stress index, δ = Absorption coefficient, SFR = Sap-flow rate, T B = Brightness temperature, EVI = Enhanced vegetation index, τ = Leaf optical depth, LA = Leaf surface area, EWT = Equivalent water thickness, NMDI = Normalized multi-band drought index, T S = Land surface temperature, GI = Greenness index, TCARI = Transformed chlorophyll absorption in reflectance index, WBI = Water band index, WI = Water index, RI = Refractive index. Figure 12 shows the multiple correspondence analysis (MCA) conducted on 72 collected studies (two studies [29,68] are not included in the MCA) and described by 6 variables (estimator, technology, target, platform, adjustment and R 2 coefficient, categorized as high, R 2 > 0.5; low R 2 < 0.5 or NA for unknown R 2 ) in order to bring up the underlying structures in the dataset. MCA shows that different studies cluster separately according to the target of measurement. Thus, methodologies with the soil moisture as target are characterized by using the NDVI, NDWI, MR, TGI and R λ estimators, with a linear adjustment and the use of cameras and radars with satellite as platform. When the crop canopy is the target, the adjustment is quadratic, several estimators are used, R 2 coefficient can be low and, in addition of radars, radiometers and spectrometers, scatterometers and a combination of sensors are used.
Separately from the studies directed to canopy and soil measurements, the studies targeted to leaves and trunks are grouped. Thus, when the target is leaf water content, the adjustment can be cubic or exponential, based on different electromagnetic spectroscopies and using several estimators. This cluster also includes methods based on g s measured with porometers. Finally, trunk based methods are characterized by technologies based on electromagnetic spectroscopies in combination with radars and specific methods such as SFR gauges, PDPs or NMR.
In addition, recent research supports that the efficiency of the evaluation of the plant water status by direct measurements is inefficient and costly, and could be improved for implementing it in commercial agriculture given the amount of data points required for a successful interpolation and, thereby, the use of remotely acquired indices provides reliable information for informed management [116]. These technologies also highlight the connection between the soil to plant water status with the primary and secondary metabolism in wine grapes evidencing their potential use in precision viticulture [115]. Finally, these authors conclude that the promptness and efficiency of proximal sensing can be transformed into realistic utilization, which can be significantly beneficial in large fields.

Conclusions
From the 77 studies analyzed in this review, it is possible to assert that the use of remote sensing to monitor plant water content is well established around different countries over the last decades, allowing farmers to control the irrigation management and to avoid unnecessary periods of water shortage and a needless waste of valuable water. According to the year of publication, the number of published papers has increased exponentially over time, mainly due to the technological development of a wide variety of cutting-edge sensors and platforms in recent years. The revision of the literature shows that the four most common sensors or technologies used for the estimation of crop water status are radiometer, CW spectroscopy, camera and THz-TDS. These technologies are mainly applied to canopy, followed by soil, leaves and trunks. This work also highlights the difficulty of comparing different approaches for crop water status assessment given the diversity of employed technologies, operating frequencies, monitoring targets, monitoring methods and estimating error. NDVI, backscattering coefficients, spectral reflectance, reflection coefficient and dielectric constant were the preferred water indicators in the reviewed studies. Finally, data shows the high degree of correlation between the estimated and measured data for the different technologies and monitoring methods. Altogether, this study evidences the potential use that remote sensing currently has in water management in crop production systems.