Analysis of Airborne Optical and Thermal Imagery for Detection of Water Stress Symptoms

: High-resolution airborne thermal infrared (TIR) together with sun-induced ﬂuorescence (SIF) and hyperspectral optical images (visible, near- and shortwave infrared; VNIR/SWIR) were jointly acquired over an experimental site. The objective of this study was to evaluate the potential of these state-of-the-art remote sensing techniques for detecting symptoms similar to those occurring during water stress (hereinafter referred to as ‘water stress symptoms’) at airborne level. Flights with two camera systems (Telops Hyper - Cam LW , Specim HyPlant ) took place during 11th and 12th June 2014 in Latisana, Italy over a commercial grass ( Festuca arundinacea and Poa pratense ) farm with plots that were treated with an anti-transpirant agent (Vapor Gard ® ; VG) and a highly reﬂective powder (kaolin; KA). Both agents affect energy balance of the vegetation by reducing transpiration and thus reducing latent heat dissipation (VG) and by increasing albedo, i.e., decreasing energy absorption (KA). Concurrent in situ meteorological data from an on-site weather station, surface temperature and chamber ﬂux measurements were obtained. Image data were processed to orthorectiﬁed maps of TIR indices (surface temperature ( T s ), Crop Water Stress Index (CWSI)), SIF indices (F 687 , F 780 ) and VNIR/SWIR indices (photochemical reﬂectance index (PRI), normalised difference vegetation index (NDVI), moisture stress index (MSI), etc.). A linear mixed effects model that respects the nested structure of the experimental setup was employed to analyse treatment effects on the remote sensing parameters. Airborne T s were in good agreement ( ∆ T < 0.35 K) compared to in situ T s measurements. Maps and boxplots of TIR-based indices show diurnal changes: T s was lowest in the early morning, increased by 6 K up to late morning as a consequence of increasing net radiation and air temperature ( T air ) and remained stable towards noon due to the compensatory cooling effect of increased plant transpiration; this was also conﬁrmed by the chamber measurements. In the early morning, VG treated plots revealed signiﬁcantly higher T s compared effect on the plant while photosynthetic efﬁciency of light reactions remained unaffected. VNIR/SWIR indices as well as SIF indices were highly sensitive to KA, because of an overall increase in spectral reﬂectance and thus a reduced absorbed energy. On the contrary, the TIR domain was highly sensitive to subtle changes in the temperature regime as induced by VG and KA, whereas VNIR/SWIR and SIF domain were less affected by VG treatment. The beneﬁt of a multi-sensor approach is not only to provide useful information about actual plant status but also on the causes of biophysical, physiological and photochemical changes.


Introduction
In the context of climate change and an increasing global water scarcity, water deficit stress (normally shortened to water stress) is one of the most critical abiotic stressors to plant growth. In order to increase the quantity and quality of food production using a reduced amount of water, the detection and quantification of plant stresses is of major interest for agriculture in general and precision farming in particular [1]. Remote sensing provides powerful tools in different spectral domains for spatiotemporal monitoring of water stress [2].
Plant response to water stress is expressed by a variety of physiological changes (e.g., stomatal behaviour, and leaf water content), biophysical changes (energy balance, leaf and canopy structure, and biomass and yield) as well as photochemical processes [3][4][5]. Accordingly, various attempts have been made to detect these changes using remotely sensed signatures. To date, the main remote sensing imaging techniques for plant water stress detection are thermal imaging (TIR; 8-14 µm), visible, near-and shortwave infrared reflectance (VNIR/SWIR; 0.4-2.5 µm), and sun-induced fluorescence (SIF; 0.69 and 0.76 µm).
TIR imaging has been well studied for water stress detection, e.g., [6][7][8]. The underlying principle is that plant temperature rises with increasing water stress in comparison to a well-watered plant due to decreasing evaporative cooling through stomatal closing [9]. Since stomatal closure is one of the first responses to water stress, plant temperature as measured by TIR sensors can be used to detect water stress pre-visually [7,10,11]. Plant temperature, however, is not solely governed by the plant water supply but also by the actual micro-meteorological conditions (i.e., solar radiation, wind speed, leaf boundary layer resistance and vapour pressure deficit (VPD)). Therefore, as an alternative to using absolute temperature, several temperature-based indices were developed over the last few decades with the aim of compensating for varying meteorological conditions. For example, the prominent Crop Water Stress Index (CWSI) [12] by means of artificial references does not require any additional meteorological measurements to be calculated. The usefulness of temperature-based indices has recently been demonstrated in several airborne studies [13][14][15][16]. Currently, the main limitation of all temperature-based (i.e., absolute temperature and temperature-based indices) approaches arises from the use of broadband infrared cameras with erroneous temperature retrieval by assuming an emissivity value that is pre-determined (e.g., ε = 0.97 for vegetation) and constant over the spectral range from 8 to 14 µm. For example, an error of 1% emissivity results in an absolute temperature error of 1 K [7]. However, very recently hyperspectral TIR airborne imagers such as Telops Hyper-Cam LW, Itres TASI-600, or Specim AisaOWL have become available. These devices allow for stable temperature and emissivity separation (TES) and very accurate temperature retrieval (i.e., <0.5 K) by measuring the emitted radiation in many narrow bands [17]. The importance of airborne (hyperspectral) TIR remote sensing lies in the possibility of bridging the gap between ground-based thermography e.g., [11,18] and proposed satellite missions, such as ECOSTRESS (Ecosystem Spaceborne Thermal Radiometer Experiment on Space Station, [19]), HyspIRI (Hyperspectral Infrared Imager, [20]) or the hyperspectral mission concept HiTeSEM (High-resolution Temperature and Spectral Emissivity Mapping, [21]).
In the VNIR/SWIR domain, various spectral vegetation indices (VI) have been developed, each for a specific purpose such as canopy water content (MSI; [22]), greenness or fractional vegetation cover (NDVI; [23]), or photosynthetic activity (PRI; [24]). As these plant parameters are somehow related to plant water status (e.g., photosynthetic activity, chlorophyll and water content may decrease under water stress conditions [3,25,26]), these VIs may serve as indicators of water stress [13][14][15]. Among the large variety of available indices, the PRI (photochemical reflectance index; [24]) is related to non-photochemical heat dissipation, which may be linked to water availability and photosynthetic efficiency. In fact, PRI is heavily influenced by canopy effects and pigment content, but Suarez et al. [27][28][29] suggest that the use of PRI in combination with radiative transfer modelling, accounting for these effects, can effectively provide very good means for monitoring water stress in crops and natural vegetation at airborne level. However, the ability of the PRI to be used for water stress detection is not conclusive at a small-scale experimental ground and airborne level [13,14,18].
At the same time, remote sensing of sun-induced fluorescence (SIF) has become increasingly popular over the last decade with a variety of studies ranging from ground-based experiments to airborne campaigns and even towards satellite missions (see Meroni et al. [30] for a comprehensive review). In principle, radiative energy absorbed by leaf chlorophyll is processed along three competing pathways: (i) conversion of photosynthetic active radiation (PAR) to sugars through photosynthesis, (ii) the re-emission of non-used energy through chlorophyll fluorescence or (iii) dissipation of heat [31]. Hence, these three processes are in competition with each other; variation in one of the processes affects the others [32]. For example, plant photosynthetic efficiency is reduced under environmental stress conditions (i.e., water or nutrient shortage, heat stress, and other types of stress) due to plants protective mechanisms (e.g., leaf rolling reduces the plants surface and therefore the PAR (Photosynthetically Active Radiation) absorbance; stomatal closure reduces water loss and CO 2 uptake). In this context, SIF is expected to be a direct indicator of photosynthetic efficiency and plant stress, although further studies are needed to establish a consistent basis for robust assessment [31,33]. The HyPlant sensor facilitates sub-nanometre airborne acquisitions from the red (0.68 µm) towards the far red (0.78 µm) spectral range. Current experiments using the HyPlant sensor demonstrated the capability of quantitative plant stress detection at airborne level using both red and far red chlorophyll fluorescence peaks [32].
These preceding studies suggest that plant water stress symptoms may be detectable by means of the three approaches (TIR temperature indices, VNIR/SWIR VIs and SIF indices). However, since they are based on different models and physiological processes, the sensitivity and suitability of these approaches may vary depending on the target and environmental conditions. Therefore, to examine the detectability of water stress symptoms by the three approaches, we conducted a field-scale experiment specifically designed to tentatively evaluate the potential of state-of-the-art remote sensing at airborne level. Hyper-Cam LW as well as HyPlant sensors were flown over a commercial grass farm. Water stress symptoms were simulated by treating grass surfaces with two different chemical agents and comparing them to untreated Control (CR) plots.
The first agent, the anti-transpirant Vapor Gard ® (VG; Miller Chemical and Fertilizer, Hanover, PA 17331, USA), is composed of di-1-p-menthene, a natural terpene polymer. The emulsion surrounds the leaves with a thin film and is supposed to reduce the plants water transpiration by limiting stomatal conductance. Thus, VG is recommended for farming to reduce water loss and prevent water stress in times of limited water availability [34][35][36]. Because of decreasing stomatal conductance, transpiration is also reduced and consequently plant surface temperature is expected to be significantly increased in comparison to a CR plot. Furthermore, VG does not only have an effect on the permeability of water but also reduces CO 2 uptake rates [37]. The second agent, kaolin (KA) is a highly reflective white powder and can be dissolved in water and sprayed on plants. The white color of KA increases the albedo of the plant surface and therefore reduces the absorbed light energy (APAR).
Thus, we assume that VG treated plants (in comparison to CR plants) have: (i) decreased transpiration, (ii) increased leaf temperature, and (iii) decreased photosynthetic activity, while keeping relatively high leaf water and chlorophyll content (in case of enough soil moisture). We assume further that KA treated plants (in comparison to CR plants) have: (iv) decreased absorbed radiation resulting in decreased transpiration under energy-limited conditions (in case of sufficient soil water), but also (v) little change in stomatal conductance, leaf temperature, photosynthetic efficiency, leaf water and chlorophyll content.
However, it is still unknown how well they allow the detection of plant water stress symptoms at airborne level. Thus, the overall aim of this study was to evaluate the capability of selected TIR and VNIR/SWIR hyperspectral remote sensing indices, as well as SIF remote sensing for detection of plant water stress symptoms. The specific objectives were: (i) examining if different temperature-based indices (e.g., CWSI, T s , T s -T air ), traditional VNIR/SWIR indices (e.g., PRI, NDVI, MSI), and SIF indices reveal differences over grass plots treated with chemical agents (VG and KA); (ii) an assessment of why the different indices may or may not have changed from treatments with respect to the underlying physiological processes; (iii) specifically examining diurnal changes in temperature-based index values regarding different treatments.

Experimental Design
In June 2014, a field-scale experiment was conducted on a 2 ha commercial grass (Festuca arundinacea and Poa pratense at the ratio of 90% to 10%, respectively) farm near to Latisana, northeast Italy (45 • 46 46.1"N, 13 • 00 50.5"E) ( Figure 1). The grass was cut at a height of 30-50 mm. All plots were well watered with an irrigation of 40 mm per week. The soil texture was silty clay. Plant characteristics were modified by two chemical treatments (VG and KA) to simulate symptoms similar to those occurring during water stress. A replicated block design with three replicates was adopted ( Figure 1). The average plot size was 12 × 9 m for CR and VG and 9 × 9 m for KA plots. In the evening of 10th June 2014, tagged plots were treated by 15 L of an emulsion of 3% VG and 15 L of 60 g L −1 KA. assume further that KA treated plants (in comparison to CR plants) have: (iv) decreased absorbed radiation resulting in decreased transpiration under energy-limited conditions (in case of sufficient soil water), but also (v) little change in stomatal conductance, leaf temperature, photosynthetic efficiency, leaf water and chlorophyll content. However, it is still unknown how well they allow the detection of plant water stress symptoms at airborne level. Thus, the overall aim of this study was to evaluate the capability of selected TIR and VNIR/SWIR hyperspectral remote sensing indices, as well as SIF remote sensing for detection of plant water stress symptoms. The specific objectives were: (i) examining if different temperature-based indices (e.g., CWSI, Ts, Ts-Tair), traditional VNIR/SWIR indices (e.g., PRI, NDVI, MSI), and SIF indices reveal differences over grass plots treated with chemical agents (VG and KA); (ii) an assessment of why the different indices may or may not have changed from treatments with respect to the underlying physiological processes; (iii) specifically examining diurnal changes in temperaturebased index values regarding different treatments.

Experimental Design
In June 2014, a field-scale experiment was conducted on a 2 ha commercial grass (Festuca arundinacea and Poa pratense at the ratio of 90% to 10%, respectively) farm near to Latisana, northeast Italy (45°46′46.1″N, 13°00′50.5″E) ( Figure 1). The grass was cut at a height of 30-50 mm. All plots were well watered with an irrigation of 40 mm per week. The soil texture was silty clay. Plant characteristics were modified by two chemical treatments (VG and KA) to simulate symptoms similar to those occurring during water stress. A replicated block design with three replicates was adopted ( Figure 1). The average plot size was 12 × 9 m for CR and VG and 9 × 9 m for KA plots. In the evening of 10th June 2014, tagged plots were treated by 15 L of an emulsion of 3% VG and 15 L of 60 g L −1 KA.
During the experiment, airborne and in situ data were collected. Airborne data were acquired with two camera systems (Telops Hyper-Cam LW, Specim HyPlant) during 11th and 12th June 2014. Concurrent in situ meteorological data from an on-site weather station, surface temperature and chamber flux measurements were obtained.  During the experiment, airborne and in situ data were collected. Airborne data were acquired with two camera systems (Telops Hyper-Cam LW, Specim HyPlant) during 11th and 12th June 2014. Concurrent in situ meteorological data from an on-site weather station, surface temperature and chamber flux measurements were obtained.

Overview
The airborne flight campaign was organised within the framework of an annual FLEX-EU campaign. The experimental site was overpassed three times during 11th and 12th June from both Hyper-Cam LW and HyPlant airborne sensor, respectively. TIR images were acquired using Telops Hyper-Cam LW on 11th June 2014 in the early morning at 09:18, late morning at 10:48 and at midday near solar noon at 12:51 (all Central Europe Solar Time, CEST). Unfortunately, synchronised overpasses of both airborne sensor were not possible. Thus, the nearest available three HyPlant imagery acquisitions were selected. The first flight occurred on 11th June at 14:52, whereas, second and third overpasses were performed one day later (12th June) at 12:02 and 13:40 (all CEST), respectively.
Hyper-Cam LW sensor was flown at a survey height of 1430 m above ground resulting in a ground sampling distance (GSD) of 0.5 m per pixel. HyPlant sensor was flown at 600 m height resulting in 1 m GSD. To obtain a consistent dataset, HyPlant images were co-registered to Telops Hyper-Cam LW spatial resolution (0.5 m). Nearest neighbour interpolation was used for resampling. Over each plot a circular region of interest (ROI, size: 32 pixels) was defined and the corresponding image values were extracted and stored in a database. The ROI were defined to be circular around the central pixel of each plot, in order to prevent any boundary effects. Each ROI covers a grass surface area of approximately 8 m 2 .

Thermal Images
TIR airborne data were recorded using a Hyper-Cam LW (Telops Inc., Québec City, Québec, Canada, www.telops.com) mounted on a customized airborne platform. The base instrument is a FTIR (Fourier transformed infrared) spectrometer consisting of a 320 × 256 pixel MCT (mercury cadmium telluride) focal plane array (FPA) detector with up to 0.25 cm −1 spectral resolution in the spectral range of 7.7 to 11.5 µm (1298.7 to 869.6 cm −1 ). The camera system is provided with two internal blackbodies with a known emissivity (ε ≥ 0.99) for radiometric calibration. The airborne platform combines the Image Motion Compensator (IMC) mirror, an Inertial Measurement Unit (IMU) consisting of GPS and Inertial Navigation System (INS), as well as a visible bore-sighted camera mounted on a base plate [17]. The IMC mirror compensates for airborne forward motion and guarantees gapless image acquisition. Furthermore, the IMC balances airborne pitch and roll, while the stabilization platform reduces airborne vibrations and airborne yaw. The IMU recordings enable ortho-rectification and geo-referencing of the collected data. TIR images were acquired with a spectral sampling distance of 6.65 cm −1 with 8 cm −1 FWHM (Full Width at Half Maximum) in the spectral domain from 878 to 1277 cm −1 resulting in 60 spectral bands.
In order to retrieve LST (Land Surface Temperature) from Hyper-Cam LW data cubes several radiometric and geometric processing steps were performed. First, raw interferograms were transformed to spectral radiances (W m −2 sr −1 cm −1 ) using a Fourier transformation, followed by a two-point calibration based on the internal blackbodies with known temperature and emissivity (ε ≥ 0.99). This process was conducted using Telops Reveal calibrate software resulting in single image cubes in units of spectral radiance. Second, the single image cubes were orthorectified and geo-referenced using the information of the IMU, before they were mosaicked to single flight lines. This geometric correction was performed using Telops Airborne software. Third, the so-called "Blackbody Fit" approach [38] was performed to retrieve LST from spectral radiances. This simple approach needs a minimum of a priori knowledge to fit the measured spectral radiances of each pixel to a Planck-curve with a known temperature. Following Schlerf et al. [17] the fit was performed in the spectral range between 880 cm −1 (11.4 µm) and 912 cm −1 (11.0 µm) and with an assumed emissivity of 97% (ε = 0.97), which is in good agreement with common findings about plant emissivity e.g., [18,[39][40][41]. The best fit was then found by iteratively changing the temperature of the Planck-curve and using a non-linear least square curve fitting approach to assess the goodness of fit [42].
Any atmospheric distortions within the optical path (i.e., upwelling and downwelling radiance as well as transmittance) was neglected for the following reasons: (i) because of the low survey altitude, the atmosphere imposes limited effects on the recorded signal; (ii) the small-scale study site limited horizontal variations in atmospheric gas concentrations and (iii) only relative differences in temperature (i.e., temperature-based indices) were studied.
To validate the accuracy of LST retrieval, airborne temperature (T airborne ) were compared to ground measured temperatures (T ground ) of a water body. T ground was measured at the surface of a water filled pool (366 cm in diameter) along a profile (i.e., along the diameter of the pool) of five thermocouples. Thermocouples have a certificated accuracy (ISO 17025) of 0.5 K. T airborne was extracted from corresponding TIR images as an average of four pixels within the pool.
Since absolute plant temperatures are very sensitive to rapidly changing meteorological conditions and therefore are not very suitable to detect water stress symptoms and diurnal changes, two additional temperature-based indices were calculated: (i) the difference between plant surface and air temperature (T s -T air ) based on the approach of Stress Degree Day (SDD) [43] and (ii) the in-scene Crop Water Stress Index (CWSI) approach following the idea of using artificial wet and dry reference surfaces [44]. CWSI was retrieved using the following equation: where T s is the plant surface temperature, T wet for the lower boundary of plant surface temperature, assuming a leaf with stomata fully open and a potential transpiration rate of 100% and T dry is the upper boundary compared to a non-transpiring leaf with stomata completely closed. Since the fleece material artificial reference surfaces located next to the treatment plots turned out unsuitable for this study (i.e., wet reference surface was warmer than coolest CR pixel and dry reference surface was up to 15 K warmer compared to the warmest pixel of the VG treatment), CWSI was determined by an empirical approach. The lower limit (T wet ) was calculated using the coolest 5% of the CR treatment pixel [45], assuming that the CR plots were well watered. The upper boundary was calculated by adding 5 K to the current T air (i.e., T dry = T air + 5 K). This approach has been previously suggested by Meron et al. [46] and Cohen et al. [6].

Hyperspectral Optical Images
HyPlant airborne sensor was developed at Forschungszentrum Jülich (Jülich, Germany) in cooperation with Specim Spectral Imaging Ltd. (Oulu, Finland, www.specim.fi). HyPlant was designed as an airborne demonstrator of ESA's Fluorescence Explorer (FLEX) satellite mission. The push-broom imager consists of two modules. First, the dual-channel module is measuring reflected surface radiance in the wavelength range of 0.38-0.97 µm with a spectral resolution of 4 nm (FWHM) at a sampling interval of 1.7 nm and from 0.97-2.5 µm with 13.3 nm FWHM at 5.5 nm sampling interval, respectively. Second, the fluorescence module was especially designed to retrieve the emitted fluorescence signal within the atmospheric oxygen bands, O 2 -A (760 nm, F 760 ) and O 2 -B (687 nm, F 687 ). Therefore, the fluorescence module provides a very high spectral resolution of 0.25 nm (FWHM) covering the spectral range from 670-780 nm resulting in 1024 spectral bands. Therefore, HyPlant is the only available airborne imaging spectrometer offering the detection of the two broadband chlorophyll fluorescence emission peaks at 687 and 760 nm. Further technical details, sensor calibration and image processing are documented in Rascher et al. [33].
The emitted sun-induced chlorophyll fluorescence F 760 and F 687 were retrieved from HyPlant's fluorescence module images using a Frauenhofer Line Depth (FLD) approach. In detail, the fluorescence retrieval implemented for HyPlant is based on the 3FLD approach as initially introduced by Maier et al. [47] and the iFLD method as proposed by Alonso et al. [48]. The method was fundamentally refined and adapted to high performance spectrometers and for airborne fluorescence retrieval. In essence it is complemented with the simulation of atmospheric components using MODTRAN5 [49] in combination with the MODTRAN5 interrogation technique as introduced by Verhoef and Bach [50], and using an empirical constraint based on non-fluorescing reference surfaces (see Wieneke et al. [51] for a recent detailed description of the retrieval method).
Additionally, several reflectance-based narrowband VIs (Table 1) related to plant physiology, structure and water content were calculated from the HyPlant dual-channel module to determine their ability to detect plant water stress symptoms induced by chemical treatments. PRI [24] was calculated using the average of three spectral bands in HyPlant closest to 531 nm and 570 nm, respectively (centre wavelength ± 1 band). As simple ratio indices related to water content we applied the Water Index (WI) [52], the Leaf Water Index (LWI) [53] and the Moisture Stress Index (MSI) [22], which are sensitive in the domain of the water absorption bands at 970 nm, 1450 nm and 1600 nm respectively. Additionally, the Simple Ratio (SR) [54] index and the Normalised Difference Vegetation Index (NDVI) were evaluated [23]. All water content and greenness related indices were calculated with a spectral window of centre wavelength ± 4 HyPlant spectral bands.

Meteorological Data
A weather station was set up close to the experimental plots within a non-treated grass area. The station was equipped with a four-band net radiometer (NR01-L Campbell Sci., Logan, UT, USA), an air temperature and humidity sensor (HMP45AC, Vaisala, Helsinki, Finland), measurement of soil water content using Time-domain Reflectrometry (TDR CS616, Campbell Sci., USA), and a canopy temperature sensor (Calex Electronics, Leighton Buzzard, UK). All of the sensors were connected to a data logger (CR1000, Campbell Sci., USA) and acquired at 0.1 Hz, then averaged half hourly. As vapour pressure deficit (VPD) is a more sensitive indicator of water vapour condition than relative humidity (RH) and therefore describes the interaction of plants with the intervening atmosphere more precisely, VPD (hPa) was calculated using Equation (2) [55]: where e s is the saturation vapour pressure in mbar, L is the latent heat of vaporization (2.5 × 10 6 J kg −1 ), Rv is the specific gas constant for water vapour (461 J K −1 kg −1 ) and T is the current air temperature (K). RH is the relative humidity (%). Non-steady-state flow-through a chamber system consisting of transparent and non-transparent chambers was used to measure net ecosystem exchange (NEE) and ecosystem respiration (Reco) fluxes, respectively. Both chambers had a volume of 0.3 m 3 and their dimensions were 0.78 × 0.78 × 0.50 m. Chambers were equipped with fans, temperature sensor and a vent to equilibrate pressure changes during measurements [56]. Water vapour fluxes were calculated based on measurements taken with transparent chamber made from 3 mm thick PLEXIGLAS (Clear 0A000 GT, Evonik Industries, Weiterstadt, Germany) [57,58]. Gas concentration changes in the chambers were measured with an infrared gas analyser (LI-840, LICOR, Lincoln, NE, USA) installed in the portable control box equipped with a pump and CR1000 data logger (Campbell Sci., Logan, UT, USA). Air was circulated between the chamber and the analyser in a closed loop with the flow rate of 0.7 L min −1 . During measurements chambers were fixed to the preinstalled soil frames (one per experimental plot) inserted to soil to 5 cm depth on 5th of June 2014. The chamber closure time used for NEE and H 2 O fluxes estimation did not exceed 1 min to avoid overheating of the chamber headspace and 2 min for the non-transparent chamber. The chambers were un-cooled in order not to reduce H 2 O fluxes due to condensation on cooling items in accordance with Chojnicki et al. [57] and Acosta et al. [58]. The chamber measurements started at 9:30 in the morning and were continued until 17:00 (CEST), but only measurements taken until the 3rd overpass at 13:40 were considered in analyses corresponding to the Hyper-Cam LW overpasses. There were 10 flux measurements taken on VG and CR plots within this time window.

Chamber Flux Measurements
H 2 O and CO 2 fluxes were calculated based on a gas concertation changes over the closure time using the linear regression type as described in Juszczak et al. [56]. In order to avoid underestimation of the fluxes due to a gas saturation, fluxes were calculated from the first 30-40 s of measurements for data with the highest regression slopes in accordance with Hoffmann et al. [59]. Coefficient of determination (R 2 ) for the fluxes considered in analyses exceeds 0.9.
Gross Ecosystem Productivity (GEP), indicating the amount of CO 2 assimilated by grass in photosynthesis, was calculated as a sum of absolute NEE and Reco values from the two consecutive measurements of NEE and Reco (Reco was measured just after NEE).
Incoming PAR was measured continuously by a PAR quantum sensor (SKP215, Skye Instruments, Llandrindod Wells, UK) installed on a weather station located approximately 100 m from the experimental plots. Considering that chamber flux measurements were taken at CR and VG plots at different times and under different PAR and temperature conditions, the H 2 O fluxes were normalised over PAR, dividing H 2 O fluxes by PAR (both in µmols m −2 s −1 ).

Statistical Analysis
For statistical analysis, we used R [60] and the R package nlme [61] to perform a linear mixed effects model to describe the nested structure of our experimental setup. An interaction term between the variables treatment and flight time was introduced as fixed effects, because we were rather interested in the comparison of the treatments at each flight time and less in the differences between the treatments accumulated over all three overpasses. As random effects, we set the variable ID (flight: treatment: plot), which describes the nesting of the variables plot within treatment and treatment within flight as well as the repeated measure problem of the three overpasses. To take care of the autocorrelation of neighbouring pixels within each plot, we included a Gaussian spatial correlation structure term in the model for the residuals, which considered a distance (d) dependent correlation following a Gaussian distribution for a certain range (r) (exp{−[(r 1−2 )/d] 2 }). The overall number of observation is n = 864, 3 treatments with 3 repetitions each by 3 flights by 32 pixels per plot. Visual inspection of residual plots did not reveal any obvious deviations from homoscedasticity or normality. For pairwise comparison we used a student t-statistic with a 0.05 Tukey-adjusted significance level, which is implemented in the R package lsmeans [62]. Figure 2 shows the diurnal variations of environmental conditions (net radiation (Rn), air temperature (T air ), surface temperature (T s ), the difference of T s and T air , and vapour pressure deficit (VPD) as well as soil water content (SWC)) for 11th and 12th of June from sunrise (05:30 CEST) until sunset (21:00 CEST). The meteorological data were retrieved from the weather station measurements (see Section 2.2). Measurements were taken over an untreated site and can therefore be considered as a Control (CR) plot. In fact, no major differences were observed between the 11th and 12th of June. For both dates, Rn steeply increased from 0 Wm −2 to a maximum of 670 Wm −2 around solar noon, before it decreased again. T air and T s also rapidly rose on the 11th June until solar noon from 19.7-32.   Figure 3 shows the diurnal course of H2O flux (solid lines) and PAR (dashed lines) for CR (green) and VG (orange) treatments observed by chamber flux measurements. During the day, H2O flux and PAR increased until solar noon. The observed variations in the measured fluxes might reflect the differences between plots (measurements were taken on three CR and three VG plots). The average H2O flux rates (n = 10, both for CR and VG) for the time between 09:30 and 13:40 (CEST) was 20% higher on average for CR plots (10.6 ± 2.61 mmols H2O m −2 s −1 ) than for VG plots (8.44 ± 3.21 mmols H2O m −2 s −1 ), although this difference was not significant (p = 0.11, Figure 3 lower left). Considering that measurements were taken at different time and under different PAR and temperature conditions, the H2O fluxes were normalised over PAR (Figure 3 lower right). After normalisation with PAR the difference between VG and CR transpiration was significant (p < 0.05 *), which proves that VG substantially reduced transpiration.

Chamber Flux Measurements
The average GEP (n = 10, both for CR and VG) for the same time between 09:30 and 13:40 (CEST) was 21% higher (p < 0.01 **) for CR plots (35.91 ± 5.15 µmols CO2 m −2 s −1 ) than for VG plots (28.51 ± 9.59 µmols CO2 m 2 s −1 ).  Figure 3 shows the diurnal course of H 2 O flux (solid lines) and PAR (dashed lines) for CR (green) and VG (orange) treatments observed by chamber flux measurements. During the day, H 2 O flux and PAR increased until solar noon. The observed variations in the measured fluxes might reflect the differences between plots (measurements were taken on three CR and three VG plots). The average H 2 O flux rates (n = 10, both for CR and VG) for the time between 09:30 and 13:40 (CEST) was 20% higher on average for CR plots (10.6 ± 2.61 mmols H 2 O m −2 s −1 ) than for VG plots (8.44 ± 3.21 mmols H 2 O m −2 s −1 ), although this difference was not significant (p = 0.11, Figure 3 lower left). Considering that measurements were taken at different time and under different PAR and temperature conditions, the H 2 O fluxes were normalised over PAR (Figure 3 lower right). After normalisation with PAR the difference between VG and CR transpiration was significant (p < 0.05 *), which proves that VG substantially reduced transpiration.

Accuracy of Temperature Images
Image temperatures retrieved from the TIR airborne sensor (Telops Hyper-Cam LW) were compared with ground temperatures of the water filled pool at each overpass (Table 2). Overall, the airborne image temperatures agree very well with the ground temperatures with less than 0.35 K averaged absolute error, which is within the accuracy of the thermocouples (0.5 K, see above). The observed error increased with increasing temperature (i.e., ∆T = 0.19 at 09:18 and ∆T = 0.49 at 12:51). In comparison, using a single band of the airborne hyperspectral scanner (AHS), Sobrino et al. [63] were only able to retrieve surface temperatures with an accuracy of 2 K. The average GEP (n = 10, both for CR and VG) for the same time between 09:30 and 13:40 (CEST) was 21% higher (p < 0.01 **) for CR plots (35.91 ± 5.15 µmols CO 2 m −2 s −1 ) than for VG plots (28.51 ± 9.59 µmols CO 2 m 2 s −1 ).

Accuracy of Temperature Images
Image temperatures retrieved from the TIR airborne sensor (Telops Hyper-Cam LW) were compared with ground temperatures of the water filled pool at each overpass (Table 2). Overall, the airborne image temperatures agree very well with the ground temperatures with less than 0.35 K averaged absolute error, which is within the accuracy of the thermocouples (0.5 K, see above). The observed error increased with increasing temperature (i.e., ∆T = 0.19 at 09:18 and ∆T = 0.49 at 12:51). In comparison, using a single band of the airborne hyperspectral scanner (AHS), Sobrino et al. [63] were only able to retrieve surface temperatures with an accuracy of 2 K.

Temperature-Based Indices
Maps of temperature-based indices ( Figure 4) as well as boxplots ( Figure 5) show diurnal changes within the experimental site and treatments. It can be noted that the relative relationships (order) among the three treatments were consistent for all indicators and all flights ( Figure 5, Appendix A, Table A1 for mean (x) and standard deviation (sd) of each index per flight and treatment).
For the CR plots, T s was lowest in the early morning (09:18) with 301.43 ± 0.21 K and increased by more than 6 K to 307.53 ± 0.44 K in the late morning (10:48). From late morning to midday (12:51), T s remained almost constant (307.99 ± 0.14 K). T s -T air was also lowest in the early morning (0.97 ± 0.21 K), increased to late morning (3.5 ± 0.44 K) and then slightly decreased by midday (2.3 ± 0.14 K). Finally, CWSI values also increased from early (0.08 ± 0.04) to late morning (0.44 ± 0.2), but decreased almost towards initial values during midday (0.11 ± 0.04). Diurnal courses of VG and KA treatment plots were similar to those of CR plots.
A distinct pattern between all temperature-based indices was recognized ( Figure 5), with KA was showing the lowest values in comparison to CR, whereas VG had the highest for all three flights (e.g., T s (KA) = 301.03 K, T s (CR) = 301.43 K and T s (VG) = 302.42 K at early morning). A pairwise comparison of the treatments (see Tukey's HSD (Honestly Significant Difference) in Appendix A, Table A2) revealed that in the early morning VG and CR were significantly different for T s and T s -T air , but not for CWSI (p = 0.59). However, VG plots could be easily distinguished from CR and surroundings by visual interpretation of the CWSI map at early morning ( Figure 4). Additionally, the boxplots of CWSI also showed a distinct difference between CR and VG treatment at early morning. Thus, all temperature-based indices showed considerably increased values for the VG treatment at early morning. By contrast, for the late morning and midday flights VG had no significant effect (p > 0.05) to any temperature-based index. Similar results were observed for KA, where for all temperature-based indices only at late morning KA was significantly different to CR.  Figure 6 shows the mean reflectance spectra from HyPlant's dual-channel module (0.38-2.5 µm) per treatment and flight. The reflectance spectra for VG plots were very similar to CR, which is also confirmed by the boxplots (Figure 7) as well as by the p-values (see Appendix A, Table A3). Furthermore, no visual differences between VG and CR were recognized in the HyPlant images from the maps (Figure 4). In contrast, KA mean reflectance spectra had an overall higher reflectance for all wavelengths. In particular, in the visible spectral range (0.38-0.7 µm) KA displayed up to 10% higher (about 5 times larger) reflectance values, underpinning the properties of KA to increase plant albedo. To test for any structural changes that may have occurred because of plant growth or stress (e.g., leaf rolling, leaf angle) in the course of the campaign, we evaluated the structural indices NDVI and SR that were derived from VNIR/SWIR spectra over the three acquisitions. Both structural indices were not significantly different at the three airborne overpasses (Figure 7), as shown by a pairwise Tukey's HSD post-hoc test (see Appendix A, Table A3). The similarity of all three NDVI and SR values suggests stable data collection (good repeatability). Again, NDVI and SR behaved in a similar manner over both chemical agents (Figure 7). They reacted highly significant to the fact that the KA increased overall reflectance (see above) in comparison to CR plots (p < 0.001 ***, see Appendix A Table A3). In contrast, VG had no significant effect on the grass surface reflectance in comparison to CR (VG and CR are almost identical, Figure 6). Further, VG had no significant effect on NDVI and SR as reported by very small mean differences between VG and CR (e.g., 0.01 at early morning and even 0 for late morning and midday for NDVI, see Appendix A, Table A1).  Figure 6 shows the mean reflectance spectra from HyPlant's dual-channel module (0.38-2.5 µm) per treatment and flight. The reflectance spectra for VG plots were very similar to CR, which is also confirmed by the boxplots (Figure 7) as well as by the p-values (see Appendix A, Table A3). Furthermore, no visual differences between VG and CR were recognized in the HyPlant images from the maps (Figure 4). In contrast, KA mean reflectance spectra had an overall higher reflectance for all wavelengths. In particular, in the visible spectral range (0.38-0.7 µm) KA displayed up to 10% higher (about 5 times larger) reflectance values, underpinning the properties of KA to increase plant albedo.  Figure 6 shows the mean reflectance spectra from HyPlant's dual-channel module (0.38-2.5 µm) per treatment and flight. The reflectance spectra for VG plots were very similar to CR, which is also confirmed by the boxplots (Figure 7) as well as by the p-values (see Appendix A, Table A3). Furthermore, no visual differences between VG and CR were recognized in the HyPlant images from the maps (Figure 4). In contrast, KA mean reflectance spectra had an overall higher reflectance for all wavelengths. In particular, in the visible spectral range (0.38-0.7 µm) KA displayed up to 10% higher (about 5 times larger) reflectance values, underpinning the properties of KA to increase plant albedo. To test for any structural changes that may have occurred because of plant growth or stress (e.g., leaf rolling, leaf angle) in the course of the campaign, we evaluated the structural indices NDVI and SR that were derived from VNIR/SWIR spectra over the three acquisitions. Both structural indices were not significantly different at the three airborne overpasses (Figure 7), as shown by a pairwise Tukey's HSD post-hoc test (see Appendix A, Table A3). The similarity of all three NDVI and SR values suggests stable data collection (good repeatability). Again, NDVI and SR behaved in a similar manner over both chemical agents (Figure 7). They reacted highly significant to the fact that the KA increased overall reflectance (see above) in comparison to CR plots (p < 0.001 ***, see Appendix A Table A3). In contrast, VG had no significant effect on the grass surface reflectance in comparison to CR (VG and CR are almost identical, Figure 6). Further, VG had no significant effect on NDVI and SR as reported by very small mean differences between VG and CR (e.g., 0.01 at early morning and even 0 for late morning and midday for NDVI, see Appendix A, Table A1). To test for any structural changes that may have occurred because of plant growth or stress (e.g., leaf rolling, leaf angle) in the course of the campaign, we evaluated the structural indices NDVI and SR that were derived from VNIR/SWIR spectra over the three acquisitions. Both structural indices were not significantly different at the three airborne overpasses (Figure 7), as shown by a pairwise Tukey's HSD post-hoc test (see Appendix A, Table A3). The similarity of all three NDVI and SR values suggests stable data collection (good repeatability). Again, NDVI and SR behaved in a similar manner over both chemical agents (Figure 7). They reacted highly significant to the fact that the KA increased overall reflectance (see above) in comparison to CR plots (p < 0.001 ***, see Appendix A Table A3). In contrast, VG had no significant effect on the grass surface reflectance in comparison to CR (VG and CR are almost identical, Figure 6). Further, VG had no significant effect on NDVI and SR as reported by very small mean differences between VG and CR (e.g., 0.01 at early morning and even 0 for late morning and midday for NDVI, see Appendix A, Table A1). The performance of reflectance-based VIs to detect plant (water) stress was tested using PRI, as well as water content related indices (i.e., WI, LWI, MSI). As KA greatly changed the reflectance in the green spectral range (Figure 6), the PRI was also greatly affected in comparison to CR. No significant differences between VG and CR plots could be recognized for PRI. However, PRI values are slightly increased in the VG treatment ( Figure 7). Furthermore, the water absorption bands at 0.97 and 1.45 µm, which were used in the WI and LWI, respectively, were also highly sensitive to the KA treatment compared to CR. In fact, the KA plots could be easily distinguished from the surrounding areas looking at PRI, NDVI and LWI maps (Figure 4). Just MSI was less sensitive and only indicated significantly differences for late morning and midday.

VNIR/SWIR Indices and Sun-Induced Fluorescence (SIF)
As a fundamental difference to reflectance-based VI, SIF is an emitted radiation. Therefore, SIF is not directly affected by changing the reflectance properties of grass surface by chemical agent, but rather indirectly by induced plant photochemical changes. F687 shows slightly negative values for KA. However, the overall pattern of F687 and F760 displays consistently lower values in KA plots ( Figure  7), except for F760 at the first HyPlant overpass. The differences between KA and CR were up to 1.6 mW nm −1 sr −1 m −2 for F687 and 0.7 mW nm −1 sr −1 m −2 for F760 during the second HyPlant overpass. Thus, The performance of reflectance-based VIs to detect plant (water) stress was tested using PRI, as well as water content related indices (i.e., WI, LWI, MSI). As KA greatly changed the reflectance in the green spectral range (Figure 6), the PRI was also greatly affected in comparison to CR. No significant differences between VG and CR plots could be recognized for PRI. However, PRI values are slightly increased in the VG treatment ( Figure 7). Furthermore, the water absorption bands at 0.97 and 1.45 µm, which were used in the WI and LWI, respectively, were also highly sensitive to the KA treatment compared to CR. In fact, the KA plots could be easily distinguished from the surrounding areas looking at PRI, NDVI and LWI maps (Figure 4). Just MSI was less sensitive and only indicated significantly differences for late morning and midday.
As a fundamental difference to reflectance-based VI, SIF is an emitted radiation. Therefore, SIF is not directly affected by changing the reflectance properties of grass surface by chemical agent, but rather indirectly by induced plant photochemical changes. F 687 shows slightly negative values for KA. However, the overall pattern of F 687 and F 760 displays consistently lower values in KA plots (Figure 7), except for F 760 at the first HyPlant overpass. The differences between KA and CR were up to 1.6 mW nm −1 sr −1 m −2 for F 687 and 0.7 mW nm −1 sr −1 m −2 for F 760 during the second HyPlant overpass. Thus, F 687 and F 760 of KA plots are significantly different in comparison to CR (see Appendix A Table A3). In contrast, VG has no effect on F 687 and F 760 compared to CR.

Discussion
As shown by the results of the meteorological measurements (i.e., SWC, VPD; Figure 2) no actual plant water stress occurred during the experiment due to sufficient water supply from the soil. Nevertheless, the VG treatment had a modest but clearly measurable effect on H 2 O flux relative to PAR (reduction of relative transpiration of grass by 20% on average) as well as on CO 2 flux (reduction by nearly 21%). This proved that real, but rather mild effects on transpiration and CO 2 uptake arose from the VG treatment.
The diurnal change of T s can be explained by the effects of plant transpiration. At late morning, as a consequence of increasing net radiation and VPD (Figure 2), plant transpiration and thus, evaporative cooling increased and prevented a rise in T s . Due to the same process, T s -T air also increased from early to late morning. However, T s -T air slightly decreased from late morning to midday ( Figure 5) because T s remained constant due to higher transpiration rate while T air continued to rise due to increased net radiation ( Figure 2). Furthermore, the diurnal changes in CWSI demonstrated the same effect as described for T s -T air . Since the upper boundary (T dry ) was determined by adding 5 K to the current T air , T dry increased in the same manner as T air resulting in overall lower CWSI values at midday. Negative CWSI values in KA resulted from the fact that T s was lowest for KA and the lower boundary (T wet ) was determined by the coolest 5% of CR plots, which were partly warmer than KA.
Indeed, the mild physiological manipulations through the chemical agents (i.e., 20% reduced transpiration in VG) only induced very small effects to temperature-based indices. VG reduced plant transpiration and induced symptoms of water stress, i.e., an overall increase in T s compared to CR plots. Contrary, the KA treatment highly increased the plant albedo, and thus reduced the overall energy uptake by the plant causing a decrease in T s . In this study, only absolute temperature differences of less than 1.0 K were observed between the treatments. Nevertheless, these minor differences could be distinguished by TIR remote sensing, as observed by a distinct pattern for all temperature-based indices ( Figure 5). These results demonstrate that TIR remote sensing indices were sensitive to small T s differences induced by the chemical treatments. In comparison, a recent study has demonstrated that temperature-based indices can detect water stress from airborne multispectral TIR data under distinct water deficit conditions with large temperature difference of up to 4.7 K [13,14]. However, our results indicate that the TIR-based indices can be used to detect minor or moderate level of water stress. It is well known that the leaf temperature can be increased by various stressors including water deficit and diseases [64,65]. Therefore, TIR imaging of high spatial and spectral resolution greatly contributes to the detection and/or quantification of the physiological anomaly caused by abiotic and biotic stresses.
Basically, canopy reflectance spectra in VNIR/SWIR are affected by water content, pigment content and composition (chlorophyll, carotenoid, xanthophyll), leaf internal structure, and canopy geometry e.g., Reference [66]. Thus, in turn, the stress symptoms are detected via the spectral change associated to the response of these physiological traits to environmental stressors. In this study, we assumed little change in water and chlorophyll content as well as in structure among CR, KA and VG treatments. The response of all analysed VIs (PRI, NDVI, WI, LWI, MSI) that are closely related to either structure, water and pigment, or physiological functioning showed little difference between CR and VG. This fact suggests the consistency of these VIs in various conditions. On the other hand, the significant differences found between CR and KA imply that care should be taken in the application of these VIs to some specific leaf surface conditions (e.g., after yellow-sand dust and volcanic ash fall). In such cases, the surface reflectance spectra could be altered (as in KA treatment) without any physiological changes. Generally, the applicability of normalised VIs is based on the consistent relationship of multiple bands to the internal physiological changes. Accordingly, alteration of the relationship between spectral bands (i.e., shape of spectra) caused by external factors would strongly hamper the consistency of VIs. In fact, KA treatment increased plant albedo, but did not increase reflectance at the same ratio in all wavebands. Thus, the significant difference in KA was not caused by physiological changes, but rather indicated abnormal alteration of reflectance spectra. Since reflectance spectra is altered unusually by plant diseases such as powdery mildew and rust diseases e.g., Reference [67], similar unusual changes of VIs as in KA would be detected by hyperspectral measurement. Therefore, our results from the unique experimental design provide useful insights on the consistency and usefulness of hyperspectral reflectance measurements, especially for discrimination of causes of detected changes.
According to previous research [30], SIF decreases with a reduction of photosynthesis under high light and stressed conditions (e.g., shortage in water availability). However, in our experiment, SIF indices showed no significant difference between VG and CR treatments while the photosynthetic efficiency should have been lowered in VG plots because of reduced H 2 O and CO 2 fluxes (as proved by chamber flux measurements). This fact suggests that either the stress was too subtle to reduce photosynthesis, and thus had no influence on the efficiency of photosynthetic light reaction or that the used SIF indices were not sufficiently sensitive to detect the slight photosynthesis reduction. On the other hand, the significant reduction in SIF indices (F 687 , F 760 ) found under KA treatment is proportionally related to the reduced amount of APAR (Absorbed Photosynthetically Active Radiation) by increasing surface reflectance, and thus lowered the photosynthetic rate. Nevertheless, it is not clear whether the reduction of SIF was directly related to the change of photosynthetic photon use efficiency because the fraction of APAR (fAPAR) as well as water and chlorophyll content remained unchanged. In addition, one has to be careful about appropriate conclusions regarding the attenuation of SIF and incident PAR at leaf surfaces. It is still not obvious how SIF can be related to photochemistry and identifying a unique relationship between SIF and photochemical efficiency is challenging [31]. The slightly negative values in F 687 originate from uncertainties in the F 687 retrieval due to the largely modified surface reflectance by the KA treatment. However, our results clearly show that the subtle change in SIF indices for KA treatments was detected by the HyPlant sensor. Therefore, optical hyperspectral techniques bear the ability to detect both rapid photochemical change (e.g., photon use efficiency) and integrated physiological change (e.g., photosynthetic capacity), and to identify the causes for the changes by using multiple wavelengths. Since remotely sensed SIF is affected by the non-linear interactions of photochemical, physiological and biophysical factors, further experimental studies in combination with process-based modelling e.g., Reference [68], are needed for both scientific and industrial applications.
In fact, the interpretation of the results presented here with regard to real water stress occurrence in natural conditions is limited due to: (i) logistical challenges related to synchronous flight plan of both airborne sensor systems HyPlant and HyperCam-LW which complicated the comparison of TIR versus VNIR/SWIR and SIF water stress indices, and (ii) instead of a real water stress occurrence only water stress symptoms were simulated by the chemical agents (i.e., VG and KA). Nevertheless, the study allowed us to deduce new findings about the sensitivity and the interplay of different indices in relation to water stress symptoms as induced by chemical agents at the airborne level.
In general, biotic and abiotic plant stressors (e.g., water stress, heat stress, and diseases) often occur simultaneously and cause similar plant physiological responses (e.g., increase in plant temperature, reduced photosynthetic efficiency, and change in canopy structure). Therefore, multi-sensor and multi-temporal approaches have a great potential to obtain useful information not only about the current plant status but also on the causes of biophysical, physiological and photochemical changes. In particular, airborne remote sensing with its high spatial and temporal resolutions can bridge the gap between in situ and satellite observations.

Conclusions
In this paper, we examined the performance of different airborne hyperspectral remote sensing approaches (i.e., TIR, VNIR/SWIR and SIF) for detecting the subtle effects caused by chemical agents over a commercial grass farm. The two chemical agents (i.e., VG and KA) induced symptoms similar to those that can be expected from the occurrence of water stress. From the VG treatment, real but rather mild effects on plant transpiration and CO 2 uptake arose. From the KA treatment, plant albedo was increased.
Although a direct comparison of the spectral domains regarding their performance of real water stress detection is not reliable, the most important finding of this study is that the TIR domain was highly sensible to subtle changes in the temperature regime as induced by VG and KA, whereas VNIR/SWIR and SIF domains were less affected by the artificial water stress symptoms. Furthermore, the assessment of the observed effects of the treatments on the different indices is summarised in Table 3. First, T s was increased in VG plots due to decreased plant transpiration. Similarly, also the diurnal changes in temperature-based indices were caused by plant transpiration. Second, SIF indices were not affected in VG plots because the changes in photosynthetic efficiency were too subtle. Third, VNIR/SWIR indices were also not affected in VG plots due to unchanged leaf biochemical components. Fourth, KA treatment increased plant albedo, reduced APAR resulting in a reduction of T s , in a reduction of SIF, and in a change of VNIR/SWIR indices (see Table 3 for details). Finally, a multi-sensor approach not only provides useful information about the current plant status but also on the causes of biophysical, physiological and photochemical changes. T s was increased due to reduced transpiration.
SIF indices remained unchanged probably due to too subtle changes in photosynthetic efficiency.

KA
Indices were highly sensitive to an overall increase in reflectance and corresponding reduction of APAR (Absorbed Photosynthetically Active Radiation).
T s was reduced due to a decrease in absorbed radiation.
SIF indices were reduced due to decreased overall available absorbed energy (APAR).