Complex Analysis of the Efﬁciency of Difference Reﬂectance Indices on the Basis of 400–700 nm Wavelengths for Revealing the Inﬂuences of Water Shortage and Heating on Plant Seedlings

the Abstract: A drought, which can be often accompanied by increased temperature, is a key adverse factor for agricultural plants. Remote sensing of early plant changes under water shortage is a prospective way to improve plant cultivation; in particular, the sensing can be based on measurement of difference reﬂectance indices (RIs). We complexly analyzed the efﬁciency of RIs based on 400–700 nm wavelengths for revealing the inﬂuences of water shortage and short-term heating on plant seedlings. We measured spectra of reﬂected light in leaves of pea, wheat, and pumpkin under control and stress conditions. All possible RIs in the 400–700 nm range were calculated, signiﬁcances of differences between experimental and control indices were estimated, and heatmaps of the significances were constructed. It was shown that the water shortage (pea seedlings) changed absolute values of large quantity of calculated RIs. Absolute values of some RIs were signiﬁcantly changed for 1–5 or 2–5 days of the water shortage; they were strongly correlated to the potential quantum yield of photosystem II and relative water content in leaves. In contrast, the short-term heating (pea, wheat, and pumpkin seedlings) mainly inﬂuenced light-induced changes in RIs. Our results show new RIs, which are potentially sensitive to the action of stressors.


Introduction
Plant cultivation is an important branch of the world's economy; the consumption of plant products increases from year to year. However, environmental conditions are rather variable [1][2][3]; the frequency and magnitude of fluctuations in these conditions can be stimulated under the development of global climate change. There are numerous stressors, including non-optimal temperatures, high intensities of light, excess precipitations, etc., which can decrease plant biomass and crop quality [4,5]. In particular, the factors strongly influence photosynthesis [6,7], which is a basis of plant productivity.
A soil drought is a key adverse factor for agricultural plants [5,8]. Water exchange and photosynthesis are important targets of the drought action [9,10]. Their changes are main reasons for the decrease in the plant productivity. The drought is often accompanied by the action of an increased temperature that can stimulate negative effects of the drought on the plant productivity [11]. It is important that even short-term water shortage (hours and days), which is induced by decreased irrigation [10,12] or the action of osmolitics [13,14], can induce numerous physiological and biochemical changes in plants. The short-term water shortage strongly influences photosynthetic processes and transpiration [10,12,14]; decreasing CO 2 assimilation and efficiency of photosynthetic machinery, inducing stomata in particular, drought-induced plant changes can be revealed by measurements of NDVI and WI at different time and spatial scales [54,76,[81][82][83]. Heating-induced changes in plants can also be revealed on the basis of the measurements of the indices [81,84].
However, further investigations into using RIs for estimating physiological processes in plants remain topical. Current RIs are calculated on the basis of analysis of only small amounts of potential combinations of spectral bands; at that, even small shifts of used wavelengths can modify sensitivities of the RIs to physiological processes (e.g., earlier we showed changes in relations between PRI and photosynthetic parameters under such shifts [59]). Solving this problem can be based on a complex analysis of the efficiency of all possible RIs, which are calculated based on all combinations of wavelengths of reflected light in a measured spectrum.
Only a few works have been devoted to this complex analysis [51,[85][86][87][88][89]. Correlograms, which are heatmaps of the dependence of the correlation coefficient between RI and analyzed physiological parameters on the values of two wavelengths (used for RI calculation), are the main results of these investigations. Their use can simultaneously show relations of all possible reflectance indices with investigated physiological parameters [51,85,89]. Potentially, similar heatmaps of the significance (p) and directions of differences between spectral indices under control and experimental conditions can be used to quickly reveal RIs that are sensitive to actions of stressors including drought and heating; however, this possibility requires an experimental analysis.
The aim of our work was to conduct the complex analysis of the efficiency of using difference reflectance indices for early revealing the influences of the short-term water shortage and heating on plant seedlings. Our investigation was focused on the analysis of RIs calculated on the basis of 400-700 nm wavelengths which form a photosynthetic activity spectrum [90,91]. Considering the influence of the water shortage and heating on photosynthesis [7,10,12], the spectral range was expected to be sensitive to action of these stressors. The investigation was performed in laboratory conditions because investigations in controlled conditions are widely used approach for primary revealing effective RIs for remote sensing (see, e.g., [36] for WI, [77] for NDWI, and [33,48] for PRI).

Plant Materials and Treatments
We used two-to three-week-old seedlings of pea (Pisum sativum L., cultivar Albumen) for an investigation of short-term water shortage influence on RIs (in accordance with our previous work [92]). Plants were cultivated in a sand substrate in a Binder KBW 240 (Binder GmbH, Tuttlingen, Germany) under a 16/8 h (light/dark) photoperiod at 24 • C. The water shortage was emulated by irrigation absence for five days of the experiment. Control plants were irrigated every two days. Measurement of reflectance and photosynthetic parameters was performed every day of the water shortage for experimental and control seedlings. Quantities of investigated seedlings were 6 for plants under water shortage and 6 for control plants.
The relative water content in the sand substrate was measured as the ratio of difference between fresh and dry weights to fresh weight; the dry weight of sand was measured after 2 h of heating in a TV-20-PZ-K thermostat (Kasimov Instrument Plant, Kasimov, Russia) (about 100 • C). This content was rapidly decreased after the termination of irrigation: about 10% of water under control conditions, less than 5% of water on the first day of water shortage, and less than 0.5% on the following days of the shortage. The method was also used for estimation of the relative water content and dry weight in the leaves which were measured every day of water shortage in separate groups of seedlings. First, second, and third leaves from each investigated plant were analyzed; the relative water content and dry weight were calculated per leaf. Quantities of investigated seedlings were 6 for all experimental and control variants.
We used two-to three-week-old seedlings of pea (Pisum sativum L., cultivar Albumen), wheat (Triticum aestivum L., cultivar Zlata), and pumpkin (Cucurbita pepo L., cultivar Mozoleevskaya) for the investigation of short-term heating influence on RIs (in accordance with [92]). The seedlings were hydroponically cultivated in the Binder KBW 240. The plants were heated in the thermostat for 30 min at 46.5 • C for pea and pumpkin and at 44 • C for wheat in accordance with our previous work [92]. Measurements of leaf plant reflectance and photosynthetic parameters for each seedling were performed on seedlings without heat (control plants) and in 1 h and one day after termination of heating (experimental plants). Quantities of control and heated seedlings were 5 for pea, 6 for wheat, and 6 for pumpkin.

Measurement of Reflected Light and Maximal Quantum Yield of Photosystem II
The measurement of the reflectance indices and photosynthetic parameters was performed on the second mature leaves of investigated plants (in accordance with our previous investigation [92]), which were kept at a fixed position ( Figure 1).

Calculation of Difference Reflectance Indices and Data Analysis
The calculation of difference reflectance indices and data analysis were based on several programs, which were developed using the Python 3.8 programming language. These programs solved the following tasks: Figure 1. Schema of measurements of maximal quantum yield of photosystem II (Fv/Fm) and intensity of reflected light in leaves of pea, wheat, and pumpkin seedlings. Leaves of investigated plants were kept at a fixed angle. Measuring light (ML) refers to standard pulses of a weak blue light (460 nm, 24 μmol m −2 s −1 , 2.5 μs pulse length). Saturation pulse (SP; red light, 630 nm, 10,000 μmol m −2 s −1 , 300 ms) was used for estimation of Fv/Fm before illumination by white actinic light (AL). MLs and SPs were generated by Dual-PAM-100, including blocks DUAL-DB and DUAL-E. A halogen lamp Osram Decostar was used as a source of white AL (about 630 μmol m −2 s −1 ); duration of illumination was 9 min. Reflected light (RL) was measured by an S-100 spectrometer. Angle of light incidence on leaf surface was about 30°, angle between leaf surface and input for RL was about 30°, and angle between leaf surface and direction of SP and ML was about 60°.
(i) Calculation of all possible difference RIs according to Equation (2): where I(R1) and I(R2) are intensities of reflected light at R1 and R2 wavelengths (in the spectral range from 400 to 700 nm). The intensities were averaged in a 3 nm range for increasing accuracy of the analysis. If R2 ≥ R1 then RI was not calculated. Additionally, light-induced changes in RIs (ΔRIs) were calculated according to Equation (3):  The photosynthetic parameters and intensity of reflected light were simultaneously measured in leaves, as described in our previous works [52,80,92,93]. A halogen lamp (3000 K, 20 W, 12 V; Osram Decostar, OSRAM GmbH, Munich, Germany) was the source of white actinic light. The intensity of leaf illumination was about 630 µmol m −2 s −1 . The distance from the halogen lamp to the leaf was about 15 cm; the angle between the leaf's surface and the direction of the actinic light was about 30 • . The white actinic light was turned on after the measurement of maximal quantum yield of photosystem II (Fv/Fm). The total duration of plant illumination by halogen lamp was about 9 min.
A Dual-PAM-100 PAM fluorometer (Heinz Walz GmbH, Effeltrich, Germany) was used to measure photosynthetic parameters. Plants were adapted under dark conditions for 15 min. The dark (Fo) and maximal (Fm) fluorescence yields in photosystem II (PSII) [24,94] were measured at the saturation pulse (630 nm, 10,000 µmol m −2 s −1 , 300 ms) before initiation of actinic light illumination and measurement of reflected light. In our investigation, we estimated the damage of the photosynthetic apparatus based on the value of Fv/Fm, which was calculated as: Fv The intensity of reflected light by plant leaves was measured by using an S100 spectrometer (SOLAR Laser Systems, Minsk, Belarus) connected with a fiber optic cable. The tip of the cable was equipped with the simplest collimator as a small black tube. The geometry of the experiment is shown in Figure 1. The angle between the leaf's surface and the input of the fiber-optic cable was about 30 • , and the distance from the leaf to the fiber-optic surface was about 1.5 cm. The spectral range of the S100 was 190-1050 nm, and the spectral resolution was about 1 nm. The white panel of a QPcard 101 calibration card ver. 3 (Argraph Corp., Carlstadt, NJ, USA) was used for the calibration procedure before measuring the spectra of reflected light in leaves.
The integration time was 5 s, and continuously repeating measurements were used. We used two time points for our analysis: after turning on the halogen lamp and after 7 min of illumination. The final intensity of reflected light at both time points was calculated as the difference between this intensity under illumination and the dark rate of the signal.

Calculation of Difference Reflectance Indices and Data Analysis
The calculation of difference reflectance indices and data analysis were based on several programs, which were developed using the Python 3.8 programming language. These programs solved the following tasks: (i) Calculation of all possible difference RIs according to Equation (2): where I(R 1 ) and I(R 2 ) are intensities of reflected light at R 1 and R 2 wavelengths (in the spectral range from 400 to 700 nm). The intensities were averaged in a 3 nm range for increasing accuracy of the analysis. If R 2 ≥ R 1 then RI was not calculated. Additionally, light-induced changes in RIs (∆RIs) were calculated according to Equation (3): where RI 1 was calculated based on I(R 1 ) and I(R 2 ), which were measured immediately after turning on the halogen lamp (first time point), and RI 2 was calculated based on I(R 1 ) and I(R 2 ), which were measured 7 min after turning on the halogen lamp (second time point). (ii) Calculation of significance (p) of differences between control and experimental plants for all calculated RIs (or ∆RIs); the nonparametric Mann-Whitney U test was used for estimation of p. Directions of changes were also estimated. Two-dimensional data arrays (significances and directions of changes for each RI as function of R 1 and R 2 ) were results of this stage of the analysis. The arrays were used for direct revealing the effective RIs (or ∆RIs) or for following construction of heatmaps. NDVI (Equation (4)) and WI (Equation (5)) were additionally calculated in the experiments because the indices could be used for estimating drought and heating action [54,76,[81][82][83][84]: where I(764-796) and I(664-676) are the intensities of reflected light averaged in 764-796 nm and 664-676 nm ranges, respectively [75].

Statistics
Considering the quantity of repetitions equaling 5-6, their distributions could be differed from the normal distribution. As a result, we preliminarily used the Shapiro-Wilk test for estimation of normality of distributions of calculated RIs, ∆RIs, NDVI, WI, Fv/Fm, and relative water content in control and experiments. The most of investigated parameters had the normal distribution (data not shown); however, there were parameters which had not the normal distribution and could not analyzed by methods of parametric statistical analysis (e.g., Student's t-test). For standardization of the analysis, we used the nonparametric Mann-Whitney U test calculated by using the Python 3.8 for analysis of all investigated parameters. Medians of investigated values were used as a nonparametric analog of averaged values.
In experiments with water shortage, linear correlation coefficients between calculated reflectance indices and investigated physiological parameters (Fv/Fm, relative water content) were calculated. Medians, which were separately calculated on the basis of control and experimental values for each day of the water shortage, were used for the calculation (n = 10).

Influence of Water Shortage on Relative Water Content and Maximal Quantum Yield of Photosystem II
On the first stage of the investigation, we analyzed influencing water shortage on the relative water content, Fv/Fm, and dry weight of leaves of pea seedlings. The water shortage was emulated by the absence of irrigation for five days; this water regime quickly decreased the water content in the sand substrate. Only pea seedlings were used in this experimental variant because our preliminary experiments showed that they were sensitive to the water shortage and suitable for investigation of the influence of water deficit on photosynthetic and spectral properties of leaves [54,92]. Figure 2a shows that the water shortage decreased the relative water content in leaves. The significant decrease was initiated in the second day of water shortage and was about 17% in the fifth day. The result showed that decrease of the soil relative water content induced decrease of the leaf relative water content for the investigated time interval. Figure 2b shows that significant decrease of Fv/Fm was observed for 2-5 days of the water shortage development. The strong decrease of Fv/Fm, which probably indicated damage of photosynthetic machinery, was observed on the fourth day, and especially the fifth day of the water shortage action. In contrast, Figure 2c shows that the dry weight was not significantly changed for five days under the irrigation absence. was about 17% in the fifth day. The result showed that decrease of the soil relative water content induced decrease of the leaf relative water content for the investigated time interval. Figure 2b shows that significant decrease of Fv/Fm was observed for 2-5 days of the water shortage development. The strong decrease of Fv/Fm, which probably indicated damage of photosynthetic machinery, was observed on the fourth day, and especially the fifth day of the water shortage action. In contrast, Figure 2c shows that the dry weight was not significantly changed for five days under the irrigation absence. Figure 2. Dynamics of the relative water content (a), photosynthetic maximal quantum yield (Fv/Fm) (b), and dry wight (c) in leaves of pea seedlings subjected to the water shortage and a scatter plot between the relative water content and Fv/Fm (d). The relative water content and dry wight were calculated per leaf. Control plants were irrigated every two days, irrigation was not used for five days in experiment. R 2 and R are determination and correlation coefficients. Medians are shown in the figure. * Value significantly differed from control (p < 0.05, Mann-Whitney U test).
After that, we investigated relations between the decrease of the leaf relative water content and the Fv/Fm decrease. Medians, which were separately calculated for values in Figure 2. Dynamics of the relative water content (a), photosynthetic maximal quantum yield (Fv/Fm) (b), and dry wight (c) in leaves of pea seedlings subjected to the water shortage and a scatter plot between the relative water content and Fv/Fm (d). The relative water content and dry wight were calculated per leaf. Control plants were irrigated every two days, irrigation was not used for five days in experiment. R 2 and R are determination and correlation coefficients. Medians are shown in the figure. * Value significantly differed from control (p < 0.05, Mann-Whitney U test).
After that, we investigated relations between the decrease of the leaf relative water content and the Fv/Fm decrease. Medians, which were separately calculated for values in control and experimental seedlings in each day of the water shortage, were used for calculation of correlation coefficient. Figure 2d shows that the relative water content in leaves and Fv/Fm were strongly related at the water shortage (R was about 0.82). The result showed that the water loss was probable to induce the photosynthetic damage at the emulated water shortage. Thus, this experimental variant could be used for further investigation of the influence of the water shortage on the RIs of pea leaves.

Influence of Water Shortage on Difference Reflectance Indices
Further, we investigated influencing the water shortage on difference reflectance indices in leaves of pea seedlings (Figures 3 and 4). Heatmaps for significances and directions of water shortage-induced changes in RIs and ∆RIs, which were calculated on the basis of wavelengths in the 400-700 nm range, were constructed for each day of the water shortage action. 50% of all RIs was changed for 4-5 days of the water shortage. Percentages of changed ΔRIs were lower, e.g., about 15% was changed for 4-5 days of the soil water shortage.    ∆RI (the light-induced change in RI) was calculated as RI 2 − RI 1 . RI 1 was calculated based on I(R 1 ) and I(R 2 ) measured immediately after turning on the halogen lamp (first time point), and RI 2 based on I(R 1 ) and I(R 2 ) measured 7 min after turning on lamp (second time point). The Mann-Whitney U test was used for p calculations. Figure 3 shows heatmaps of significances and directions of water shortage-induced changes in absolute values of RIs in pea leaves in different days of the water shortage. Two large spectral regions with significant differences between control and experimental seedlings (with increased and with decreased RIs) and several small regions were observed in the first day of the water shortage. In the second day, regions with decreased RIs were absent (excluding separate pixels or small groups of pixels) and areas of spectral regions with increased RIs were strongly reduced. Since the third day of the water shortage, areas of spectral regions with significantly changes in RIs were increased. On the fifth day of the water shortage, areas of the spectral regions with significant changes in absolute values of RIs formed the most part of total area of the heatmap (about 71%). Figure 4 shows heatmaps of significances and directions of water shortage-induced changes in ∆RIs (the light-induced change in RIs) in pea leaves in different days of the water shortage. There were only small spectral regions and separate pixels in the heatmaps with significant differences between control ∆RIs in experimental and control plants in the first, second, and third days of the water shortage. Positions of the regions and pixels were varied in different days of the water shortage. Since the fourth day of the water shortage, large spectral regions with positive and negative changes in ∆RIs were observed. In the fifth day of the water shortage, areas of the spectral regions with significant changes in ∆RIs were about 30% from total area of the heatmap.
The results showed that the water shortage action could influence the number of RIs; changes in their absolute values were more expressive than changes in ∆RIs. Several spectral regions with significant differences between RIs (or ∆RIs) were observed in specific days of the water shortage. For example, a large spectral region with decreased RIs (R 1 was about 640-690 nm, R 2 was about 500-650 nm) was revealed in the first day (Figure 3a) but it was absent in the second day (Figure 3b). In contrast, there were spectral regions with significant differences which were observed for several days of the water shortage. For example, a large spectral region with increased RIs (R 1 was about 560-640 nm, R 2 was about 520-600 nm) was observed since the third day of the water shortage (Figure 3c-e).

Analysis of Efficiencies of Difference Reflectance Indices for Revealing Water Shortage-Induced Changes in Seedlings
In the next stage of our work, we directly analyzed the two-dimensional data arrays, which included significances and directions of changes for each RI as function of R 1 and R 2 . Revealing RIs, which were changed on initial stages of the water shortage and remained changed after that, was main task of this analysis because the indices could be potentially considered as markers of action of drought on plants. Table 1 shows that there were only a small portion of RIs (less than 2%), which were significantly changed for 1-5 or 2-5 days of the water shortage. In contrast, about 25% of the total quantity of RIs was changed for 3-5 days of the water shortage and more than 50% of all RIs was changed for 4-5 days of the water shortage. Percentages of changed ∆RIs were lower, e.g., about 15% was changed for 4-5 days of the soil water shortage. RIs, which were significant changed for 1-5 and 2-5 days of the water shortage, were the most interesting for analysis because they were sensitive to the earliest action of the water shortage. Table 2 shows these RIs and their correlations with the relative water content and Fv/Fm in leaves of pea seedlings. Only absolute values of RIs were analyzed. Additionally, Table 2 shows NDVI and WI, which are widely used indices for estimating the drought action on plants [54,76,[81][82][83].  Table 2 shows that there were five reflectance indices, which were significantly increased since the first day of the water shortage, and another 41 reflectance indices, which were significantly increased since the second day of the water shortage. The most of these RIs were strongly correlated with Fv/Fm and the relative water content in leaves of pea seedlings. NDVI was also strongly related to the physiological parameters; however, it was significantly changed since the third day of the water shortage. Correlation coefficients between WI and investigated physiological parameters were lower than ones for other RIs (e.g., the coefficient between WI and the relative water content was 0.63). The result was in a good accordance with literature data [76] which showed that average correlation coefficient between WI and the relative water content was about 0.66. Significant change in WI was observed in the fifth day of the water shortage. Thus, revealed RIs seemed to be more effective markers of the water shortage action in investigated plants than NDVI and WI.
Dynamics of NDVI, WI, and several revealed RIs in control and experimental pea seedlings were separately investigated. We analyzed the RI(621,442) calculated on the basis of reflected light intensities at 621 and 442 nm, RI(500,471) calculated on the basis of reflected light intensities at 500 and 471 nm, and RI(692,421) calculated on the basis of reflected light intensities at 692 and 421 nm. RI(621,442) was selected because it was significantly changed for 1-5 days of the water shortage. RI(500,471) and RI(692,421) were selected because they were calculated on the basis of different spectral regions (green-blue and red-blue regions, respectively) and were strongly correlated to physiological parameters. Figure 5 shows that all selected RIs had relatively stable dynamics in control plants; in contrast, the NDVI and, especially, WI dynamics were more varied under control conditions. A magnitude of RI(692,421) changes was more than the magnitude of NDVI changes (about 0.2 and about 0.3, respectively). The magnitudes of RI(621,442) and RI(500,471) changes were lower (about 0.015 and about 0.05, respectively); however, the magnitudes were comparable to the magnitudes of other widely-used reflectance indices (e.g., PRI [59]).

Influence of Short-Term Heating on Maximal Quantum Yield of Photosystem II and Difference Reflectance Indices
It is known [11] that development of soil drought is often accompanied by increased temperature actions; at that, durations of these actions (daily temperature changes and/or direct heating of leaves induced by sunlight) can be strongly shorter than the duration of soil drought. Considering this point, we additionally investigated influencing the short-term heating on Fv/Fm and reflectance indices. Earlier, we showed that the shortterm heating could strongly influence photosynthetic parameters and PRI in pea, wheat, and pumpkin seedlings [92]. As a result, wheat and pumpkin seedlings were additionally investigated in the experiment with the short-term heating. Experimental Fv/Fm and reflectance indices were measured in 1 h and one day after the action of high temperature on seedlings of pea, wheat, and pumpkin.
It was shown that Fv/Fm was significantly decreased in 1 h after the short-term heating in pea, wheat, and pumpkin seedlings ( Figure 6). In pea, Fv/Fm did not recover in one day after the heating. In wheat and pumpkin, Fv/Fm partially recovered in this time. Figure 7 shows that heatmaps of significances and directions of heating-induced changes in absolute values of RIs. Two large spectral regions with significant differences between control and experimental seedlings (with increased and decreased RIs) and several small regions were observed in pea and pumpkin leaves in 1 h after the heating. The spectral regions with increased RIs were strongly decreased in one day after the heating (especially, in pumpkin leaves); the spectral regions with decreased RIs were practically absent after this time interval. In wheat seedlings, there were only several small spectral regions (and separate pixels) with significantly increased absolute values of RIs in 1 h and one day after the short-term heating. Spectral regions and separate pixels with decreased absolute values of RIs were absent in wheat seedling. magnitudes were comparable to the magnitudes of other widely-used reflectance indices (e.g. PRI [59]).  several small regions were observed in pea and pumpkin leaves in 1 h after the heating The spectral regions with increased RIs were strongly decreased in one day after the heating (especially, in pumpkin leaves); the spectral regions with decreased RIs were practically absent after this time interval. In wheat seedlings, there were only severa small spectral regions (and separate pixels) with significantly increased absolute values of RIs in 1 h and one day after the short-term heating. Spectral regions and separate pixels with decreased absolute values of RIs were absent in wheat seedling.   Figure 8 shows that heatmaps of significances and directions of heating-induced changes in ∆RIs (the light-induced change in RIs). Several large spectral regions with significant differences between control and experimental seedlings (with increased and with decreased RIs) and several small regions were observed in pea, wheat, and pumpkin leaves in 1 h after heating. Some spectral ranges were specific for different investigated species. For example, a large spectral region with decreased ∆RIs (R 1 was about 525-575 nm, R 2 was about 460-540 nm) was observed in wheat seedlings in 1 h after the short-term heating but it was absent in pea and pumpkin seedlings. In contrast, some spectral ranges with significant changes in ∆RIs were observed in all investigated plants. For example, a large spectral region with decreased ∆RIs (R 1 was about 570-610 nm, R 2 was about 510-570 nm) was observed in all pea, wheat, and pumpkin seedlings in 1 h after the short-term heating.
The large areas of the spectral regions with changed ∆RIs were observed in one day after the short-term heating in pea seedlings. Some spectral regions were like the regions in 1 h after the heating; for example, the spectral region with decreased ∆RIs (R 1 was about 570-610 nm, R 2 was about 510-570 nm) or spectral region with increased ∆RIs (R 1 was about 675-700 nm, R 2 was about 400-650 nm). Only small spectral areas with significant changes in ∆RIs were observed in one day after the short-term heating in pumpkin seedlings. Finally, only separate pixels in the heatmap showed significant changes in ∆RIs in one day after the short-term heating in wheat seedlings. The results seemed to be similar to different recoveries of photosynthetic damages in different investigated species (the Fv/Fm recovery absence in pea seedlings and the partial recovery in wheat and pumpkin seedlings, Figure 6).    I(R 1 )+I(R2) , where I(R 1 ) and I(R 2 ) are intensities of reflected light at different narrow spectral bands (averaged for 3 nm). ∆RI (the light-induced change in RI) was calculated as RI 2 -RI 1 . RI 1 was calculated based on I(R 1 ) and I(R 2 ) measured immediately after turning on the halogen lamp (first time point), and RI 2 based on I(R 1 ) and I(R 2 ) measured 7 min after turning on lamp (second time point). The Mann-Whitney U test was used for p calculations.

Analysis of Efficiencies of the Water Shortage-Sensitive Reflectance Indices for Revealing Heating-Induced Changes
Finally, we analyzed the question: Could the revealed RIs, which were sensitive to water shortage-induced changes (see Table 2), be significantly changed after the short-term heating? Only pea seedlings were investigated to increase accuracy of comparison between water shortage-and heating-induced changes in absolute values of RIs.
The two-dimensional data arrays, which included significances and directions of heating-induced changes for each RI as function of R 1 and R 2 , were directly analyzed. It was shown that 31 RIs from Table 2 were significantly increased in 1 h after the shortterm heating and 2 RI was significantly increased in one day after that (Table 3). NDVI was significantly decreased in 1 h after the short-term heating; however, the significant change was absent in one day after the heating. Significant changes in WI were absent. Table 3. Changes in absolute values of NDVI, WI, and the water shortage sensitive RIs (from Table 2) after the short-term heating. The RIs from Table 2, which were not significantly changed after the short-term heating, were not shown.
The results showed that there were RIs which were sensitive to both the water shortage and heating; however, the indices were mainly sensitive to plant changes, which were observed immediately after action of the increased temperature (in 1 h after the heating). Long-term changes in plants (changes, which were observed in one day after the heating) did not influence absolute values of 44 from 46 RIs, which were sensitive to the soil drought (see Tables 2 and 3).
There are numerous reflectance indices (e.g., NDVI, PRI, WI, NDWI, and NIR V [33][34][35][36]39,42,48,[74][75][76][77][78][79]) which are calculated on the basis of narrow spectral bands of reflected light and used for plant remote sensing. Using RIs for analysis of reflected light spectra has the number of preferences including the possibility to measure based on both multi-and hyperspectral cameras, simple algorithms for calculating RIs, the use of relative units for RIs and minimization of errors related to light intensity, and the possibility to reveal specific physiological processes or specific physiological parameters related to specific wavelengths. The reflectance indices can be used for remote sensing of water shortageand heating-induced changes in plants (e.g., NDVI, WI, NDWI, or PRI [54,76,77,[81][82][83]92]).
However, current RIs are calculated on the basis of analysis of only small amounts of potential combinations of spectral bands; at that, even small shifts of the used wavelengths can modify the sensitivities of RIs to physiological processes (e.g., see [59]). As a result, the next step in development of methods of using RIs is their complex analysis, which includes estimation of efficiencies of all possible RIs, which are calculated based on all combinations of wavelengths of reflected light in the measured spectrum [51,[85][86][87][88][89]. Construction of two-dimensional heatmaps, which show the efficiency of each possible RI (e.g., correlation coefficients between the RI and physiological parameters in correlograms [51,[85][86][87][88][89]) as function of two wavelengths used for the calculation of this RI, is actively used for the complex analysis. In the current work, we analyze the similar heatmaps; however, significances and directions of water shortage-and heating-induced changes in RIs are estimated in our complex analysis. There are several important points which were shown in the current work.
First, we show that the heatmaps showing significances and directions of stressorinduced changes in RIs and ∆RIs (400-700 nm wavelength range) can be used for the complex revealing RIs and ∆RIs, which are sensitive to the short-term water shortage and heating (see Figures 3,4,7 and 8). In the heatmaps there are both spectral regions and separate pixels with significantly increased or decreased RIs and ∆RIs. It shows that both correlograms [51,[85][86][87][88][89] and heatmaps of significances and directions of changes (our results) can be used for the analysis of efficiencies of reflectance indices as markers of stressor actions.
Second, the analysis of influencing the emulated water shortage on pea seedlings shows that water shortage-induced changes in absolute values of RIs are more expressive than the changes in ∆RIs (the light-induced change in RIs): the changes in RIs are initiated earlier (Figures 3 and 4) and spectral regions with these changes are larger (Table 1). There are both "specific" (significant changes in RIs are only observed at specific days after the water shortage initiation) and "universal" (significant changes in RIs are observed for several days of the water shortage action) spectral regions with drought-sensitive RIs in the heatmaps (Figure 3).
Further analysis shows that there are five RIs increased for 1-5 days of water shortage and other 41 RIs increased for 2-5 days of the shortage (Table 2). These RIs are more sensitive to water shortage then NDVI (increased for 3-5 days of water shortage) and WI (increased in 5th day), which are widely used for the remote sensing of drought-induced changes in plants [54,76,[81][82][83]; the most of the RIs are strongly correlated to Fv/Fm and the relative water content in leaves. All 46 RIs are based on wavelengths related to absorption spectra for carotenoids and chlorophyll a and b [95]. It is additionally shown that at least some revealed RIs (e.g. RI(621,442), RI(500,471), and RI(692,421)) are more stable under control conditions than NDVI and WI ( Figure 5). It is an important property because agedependent dynamics of absolute value of reflectance index can strongly limit the RI using (e.g., changes in absolute value of PRI under control conditions can limit the index using for estimation of water shortage action on plants [92,96]). Considering absorption spectra for carotenoids and chlorophyll a and b [95], it can be additionally speculated that RI(621,442) is mainly related to the contents of carotenoids and chlorophyll a, RI(500,471) is mainly related to the contents of carotenoids and chlorophyll b, and RI(692,421) is probably related to the contents of chlorophyll a and b. Finally, 32 from these 46 RIs are also sensitive to short-term heating (Table 3), which is important considering the relations between drought and increased temperature [11].
Third, the analysis of influencing short-term heating on pea, wheat, and pumpkin seedlings shows that the heating influences absolute values of RIs in pea and pumpkin seedlings in a similar manner; in contrast, there are only minimal changes in RIs in wheat seedlings ( Figure 7). However, the analysis of ∆RIs (the light-induced change in RIs) shows that there are many spectral regions with significant changes in ∆RIs, which are similar in pea, wheat, and pumpkin seedlings at 1 h after the short-term heating (Figure 8). In one day after heating, the regions are minimal in wheat and pumpkin seedlings, but they remain in pea seedlings. Both results are in a good accordance with photosynthetic damage because Fv/Fm is significantly decreased 1 h after the short-term heating in pea, wheat, and pumpkin seedlings; in contrast, Fv/Fm partially recovers in wheat and pumpkin seedlings and is decreased in pea seedlings in one day after the heating.
It can be supposed that the differences in efficiency of RIs and ∆RIs can be related to the sensitivity of absolute values of the indices and their light-induced changes to different groups of processes. Changes in absolute values in RIs are probable to be mainly related to changes in content and composition of photosynthetic pigments [24,58], which can be affected by water shortage [97,98]. In contrast, ∆RIs are light-induced changes in RIs for 7 min of illumination. This means that ∆RIs should be related to light-induced physiological changes in this time interval. It is known that there are numerous photosynthetic changes, which are formed in this time interval (activation of electron transport chain [99], the formation of different components of non-photochemical quenching [100,101], changes in the electrochemical gradient of H + [100,102], etc.), and influence plant reflectance (e.g., [35,42,80]). Damaging the processes probably results in plant heating [92]; this damage can explain changes in ∆RIs in our experiments.
Thus, our results show that the complex analysis of efficiencies of RIs (and ∆RIs), which is based on the calculation of all reflectance indices (and their light-induced changes) in the measured spectral range of reflected light, on following the estimation of significances and directions of differences between RIs (and ∆RIs) in control and experimental plants, and on the construction of heatmaps of these significances and directions, which can be used for revealing new reflectance indices, and that are sensitive to stressors (by example of the 400-700 nm spectral range and the water shortage and heating).
It should be noted that our results are shown in laboratory experiments. It means that efficiencies of the specific RIs, which were sensitive to short-term water shortage in our work, require additional verification in field conditions in future. Potentially, efficiencies of revealed reflectance indices can be modified by complicated geometry of canopy in comparison to fixed leaf, fluctuations in environmental conditions (e.g., changes in light intensity or wind gusts), spatial heterogeneity in the investigated field (e.g., different plant growth density), etc. Additionally, it cannot be excluded that RIs, which are sensitive to early plant physiological changes under water shortage, will not be sensitive to physiological changes under long-term drought. As a result, further verification of the revealed specific indices in all noted cases is important problem of future investigations.
In contrast, it can be expected that the general proposed approach based on complex analysis of significance of changes in RIs and ∆RIs can be used for both leaf and canopy levels, different species, various stressors, different spatial and temporal scaled, etc. This means that the approach seems to be perspective for development of methods of remote sensing of action of stressors on agricultural plants in fields by multi-and hyperspectral imaging that will increase food security and support a sustainable agriculture.

Conclusions
In the current work, we performed complex analysis of the efficiency of using difference reflectance indices for early revealing the influences of the short-term water shortage and heating on plant seedlings. In laboratory experiments, it was shown that calculation of all possible RIs and ∆RIs based on 400-700 nm spectral range and construction of heatmaps of significances of differences between control and experimental seedlings could be used for revealing indices, which were sensitive to plant changes induced by water shortage and heating. Using this approach, we showed that there were a number of RIs that were more sensitive to action of water shortage than NDVI and WI; some of the indices were also sensitive to heating action. It was additionally shown that absolute values of reflectance indices were more sensitive to the action of water shortage. In contrast, light-induced changes in the indices were more sensitive to the action of heating.
New specific reflectance indices, which are sensitive to stressors in our work, can be potentially used for early revealing action of water shortage on plants; i.e., they can contribute to development of methods of plant remote sensing. However, a check of their efficiency in changeable environmental conditions, in different spatial and temporal scales, in different plant species, and under action of different stressors is an important future task because it cannot be excluded that efficiencies of the RIs may be strongly modified in these cases.
In contrast, it is very probable that the general proposed approach based on complex analysis of significance of changes in RIs and ∆RIs can be used for revealing reflectance indices, which will be effective tool for plant remote sensing at all noted cases. Therefore, another important task of future investigations is further analysis of efficiency of this approach for different experimental and field conditions.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.