Chlorophyll Fluorescence Kinetics May Be Useful to Identify Early Drought and Irrigation E ﬀ ects on Photosynthetic Apparatus in Field-Grown Wheat

: To assess the reliability and sensitivity of non-invasive optical methods to detect the early e ﬀ ects of water deﬁcit in the ﬁeld, we analyzed the time-series of non-invasive measurements obtained in a dry season in a representative collection of wheat genotypes grown in small-plot ﬁeld trials, in non-irrigated and irrigated variants. Despite a progressive water deﬁcit and signiﬁcant yield loss, the measurements indicated very minor changes in chlorophyll content or canopy cover. This corresponded well to the insigniﬁcant di ﬀ erences in spectral reﬂectance normalized di ﬀ erence vegetation index (NDVI) values. On the other hand, we identiﬁed the signiﬁcant and rapid response of fast ﬂuorescence kinetics data following the onset of irrigation. Analysis of parameters showed the main e ﬀ ects of drought were associated with changes in the amplitude of the I–P phase of the OJIP transient, indicating changes at the level of photosystem I and beyond. Statistical analyses identiﬁed the integrative parameter performance index PI tot as the most sensitive parameter, which well-reﬂects the di ﬀ erences in responses of the genotypes to water deﬁcit. Our results suggest that focusing on photosynthetic functions detected by the rapid chlorophyll ﬂuorescence records can provide more accurate information on the drought stress level, compared to the structural data obtained by absorbance or reﬂectance measurements.


Introduction
Climate change and related environmental stress factors play an increasingly important role in the performance, vulnerability, and productivity of crop plants on a global scale, with important economic impacts [1,2] and associated effects on food security, along with the expected increase of the world population [3]. This has led to increased pressure on researchers to contribute towards increased food production and quality, through development and innovation in relevant fields [4,5].
Considering the sustainability of food production, crop breeding plays a central role in ensuring food security by developing new varieties which are higher yielding, disease resistant, drought tolerant, or regionally adapted to different environments and growing conditions [6]. As conventional plant-breeding progress will not be enough to cope with the increase in demand for food production as a result of population growth [7], there exists an urgent need to accelerate breeding through improved genotyping and phenotyping methods and by efficiently using the available genetic diversity in crop germplasm [4]. New technologies are expected to simultaneously monitor genotypic and phenotypic data in the near future [8,9]. In addition to simple morphological and growth parameters, a larger contribution of physiological traits and responses measured in field conditions are required [10,11].
Field conditions are diverse with little control of environmental factors, and field experiments provide important information which is often different to that obtained in controlled environments [12]. Plants experience a range of stresses throughout their life cycle, caused by a variety of factors. Phenotyping of drought stress effects is particularly challenging, as the soil environment plays a crucial role and is difficult to simulate under controlled conditions [13]. Phenotyping under field environmental conditions remains a bottleneck for future breeding advances. Therefore, it seems necessary to develop a high-throughput phenotyping tools with improved spatial and temporal resolution, which may unlock new prospects for non-destructive field-based phenotyping in plants for a large number of traits, including those related to physiological, biotic, and abiotic stresses [14].
Assessment of the photosynthetic phenomena and apparatus can play a very important role in plant phenotyping, as it is considered to be a good indicator of overall performance or stress sensitivity of crops and genotypes. Chlorophyll fluorescence techniques have many practical applications, providing non-invasive, fast, and economic methods for the measurement of photosynthetic performance [15,16]. They can provide a large set of data on the function of PSII photochemistry obtained in only a few seconds, providing information about plant vitality, which may serve as a screening criterion. Chlorophyll fluorescence emitted from plant leaves provides information on the effects of the environmental factors on the plant's photosynthetic apparatus and the overall status of the photosynthetic systems within the leaf. Chlorophyll fluorescence can be used to measure the impact of different stress factors on photosynthesis, enabling us to detect the effects of stress before any visible external symptoms appear and before they might be detected by other methods [17,18].
Wheat (Triticum sp.) is the one of most dominant crop species in temperate countries, being used for human food and livestock feed. Wheat is the third most important crop to humankind, in terms of global production, and is a staple of many diets around the world [19]. Therefore, it is important to understand the different factors (e.g., soil fertility, availability of water, climate, diseases, or pests) affecting its crop photosynthetic activity and yield. In wheat, enormous genetic and phenotypic diversity can be observed, including a diversity in photosynthetic responses [20]; however, its efficient utilization in breeding is limited by a lack of efficient tools for testing the photosynthetic response in field collections.
In our previous studies, we identified specific parameters derived from the fast fluorescence kinetics which are useful to assess the level of nitrogen deficiency [21], effects of shading [22], or different PSII thermostability in different wheat genotypes [23]. In addition, in pot experiments, we have demonstrated that some fluorescence parameters, such as performance index, are sensitive to progressive dehydration in wheat, and thus may serve to recognize more sensitive or resistant genotypes [24]. The sensitivity of the fast fluorescence parameters to drought or osmotic stress has also been identified in barley [18], sorghum [25], rice [26], and maize [27]. There are several studies which have confirmed these findings, including a few studies using field conditions [28,29], but there is still insufficient information on the responses of PSII photochemistry on moderate, slowly progressing drought stress in field conditions and the specific effects of irrigation.
In this study, we focused on the effects of water deficit and additional water supply on the photosynthetic apparatus of wheat, indicated by the time-series of selected parameters derived from the chlorophyll fluorescence transient and additional leaf and canopy parameters, as recorded in a dry season in a collection of wheat genotypes. We compared the responses of the fluorescence parameters to water deficit and irrigation with the effects of stress on additional parameters measured at the leaf and canopy level. We also tested the ability of the techniques and parameters to recognize the different responses of the genotypes to stress and irrigation. The main objective of the present study was to identify the most suitable parameter(s) for the detection of drought stress effects and the responses of wheat genotypes under a drought stress scenario typical of a moderate climate zone.

Experimental Setup and Climate Conditions
Research Experiments were carried out with a field-grown wheat collection (field trials, Genebank of National Agricultural and The seeds were sown manually directly into the soil in autumn, into experimental plots of area 1.5 m 2 . The previous crop was field peas (Pisum sativum subsp. arvense). In the autumn, plants were fertilized with 330 kg ha −1 NPK (15-15-15), 150 kg ha −1 AMOFOS fertilizer (12% N, 52% P 2 O 5 ), and 120 kg ha −1 potassium chloride (60% K 2 O). Spring fertilization was realized as 110 kg ha −1 of ammonium nitrate with limestone (27% N). Furthermore, in autumn and spring, spraying against weeds was performed using optimal doses.
The experiments were significantly affected by the weather (Table 1), as the growing season (2016-17) was dry. In particular, the spring season was very dry and warm; except for April, which was normal. The crop moisture index (CMI), a drought index considering both temperature and precipitation to estimate the soil water balance [30], was calculated for the spring vegetation period, in order to recognize significant water deficits in field conditions.
The experimental plots were organized in two blocks: one block with a drip irrigation system (hereinafter indicated as the irrigated variant) and another block without irrigation (non-irrigated variant). The blocks were separated by a buffer (strip) of conventional wheat.
Drip irrigation was initiated in May, when the water reserves in soil decreased and the water deficit was recognized. Irrigation was then performed on regularly, once a week, with a dose of~20 mm in the irrigated variant.

Chlorophyll Fluorescence Measurements
Chlorophyll a fluorescence measurements were performed both on control and stressed variants from April to June, in regular intervals. Intact flag leaves of wheat plants were adapted to darkness for 15 min using leaf clips. Chlorophyll a fluorescence was measured using a Handy PEA (Plant Efficiency Analyzer, Hansatech Instruments, Kings Lynn, UK). After the adaptation of leaves to darkness, a single strong 1 s light pulse (3500 µmol m −2 s −1 ) was applied, provided by three light-emitting diodes (650 nm). The fast fluorescence kinetics (F 0 to F M ) was recorded from 10 µs to 1 s. For each genotype and variant, at least 10 measurements were applied. The measured data were analyzed by the JIP test [31,32]. The state of PSII photochemistry was analyzed using the parameters listed in Table 2. Table 2. Chlorophyll fluorescence parameters used in the study.

Parameter Definition Formula
Fo, Fm Minimum and maximum fluorescence intensity Amplitude of variable fluorescence in I-P phase; probability of electron transport from PSII acceptor side to the PSI electron acceptors

Measurements of Chlorophyll Content
Chlorophyll content was measured on the surface of the youngest fully developed leaf using a SPAD-502 chlorophyll meter (Konica-Minolta, Tokyo, Japan). The SPAD index values corresponding to the leaf chlorophyll content per leaf area were calculated based on the amount of light transmitted by the leaf in two wavelength regions (650 and 940 nm). The mean SPAD values were calculated as the average of ten measurements on the middle part of the leaf.

Assessment of Reflectance-Based Vegetation Index (NDVI)
The ASD FieldSpec 4 (Analytical Spectral Devices Inc., USA) broadband spectroradiometer was used to measure the spectral reflectance curves at wavelengths of 350-2500 nm. The measurements were performed directly in the field, where the sensor was held at a height of 60 cm above the plant. Measurements were performed in multiple replicates on one genotype.
In this study, we analyzed a main vegetation index-the normalized difference vegetation index (NDVI) (Equation (1)-obtained from the differential absorption of red and near-infrared reflectance bands [33], calculated as follows:

Assessment of the Crop Water Status
As the technical possibilities and organization of the measurements did not enable realization of conventional invasive assessment of leaf or plant water content, we monitored the changes of the water-sensitive bands of the spectral reflectance records, as measured by the ASD FieldSpec 4 (Analytical Spectral Devices Inc., Boulder, CO, USA) broadband spectroradiometer. In this experiment, as a measure of the canopy water status, we selected the water band index (WBI) (Equation (2) [34], based on near infra-red spectra: Although the WBI data could not be used for precise estimation of the water content, they provided rough information about the presence or absence of the differences in leaf water status between two variants.

Analysis of the Canopy Absorbance
The canopy absorbance was analyzed using the values of the light intensity measured above the canopy and near the soil level inside the canopy using a portable light sensor system Plant Canopy Analyzer (LICOR Inc, Lincoln, NE, USA).

Analysis of the Yield Components
In the spikes of individual plants, we assessed the number of grains per spike, thousand grain mass (g), and total grain mass per plant (g), representing grain yield.

Statistical Analyses
Statistical analysis was performed using two-way analysis of variance (ANOVA) followed by the post hoc Tukey HSD test (p < 0.05) using the Statistica version 9.0 software (Statsoft Inc., Tulsa, OK, USA). The factors analyzed were Irrigation (irrigated vs. non-irrigated variant) and Genotype (25 genotypes). The data presented in graphs represent mean value ± standard error. Measurements of all parameters in two variants were realized in 25 genotypes, performed nine times during the spring wheat season (March to June). Ten individuals of each genotype were analyzed using non-invasive methods. Analysis of statistical significance was performed also for a period of three analyses, before and after.

Results and Discussion
The relevance of drought is often underestimated for temperate climate regions of Central Europe, which is not generally characterized as a drought-prone region. However, water deficit frequently occurs, and drought is one of the most important climatic extremes, in terms of economic damage, along with its increasing trends associated with climate change [35,36]. A typical example was the season (Table 1) in which the presented study was performed, in which the autumn and winter period, with moderate water deficit, was followed by a dry and warm spring (partly alleviated by the weather in April), leading to the onset of soil water deficit, as shown by the values of CMI drought index ( Figure 1A). The onset of irrigation was immediately reflected in the values of the water band index (WBI; Figure 1B), which has been previously shown to be an indicator of dynamic changes in the water supply in plants [37][38][39]. The significant increase of WBI due to irrigation applied approximately four days before the measurements indicated that the non-irrigated plants had been subjected to a water deficit and, in turn, the plant water status of irrigated plants was significantly improved. On the other hand, the measurements of the light absorbance realized on the leaf and canopy level indicated that the water deficit was not associated with a decrease of the chlorophyll content in upper leaf levels, as shown by SPAD values (Figure 1C), and that the canopy structure was not significantly disturbed ( Figure 1D), providing a full leaf cover in water-deficient plots comparable to that in the irrigated variant, except during very late growth stages. This indicates that the plants were exposed to a slowly progressing moderate drought stress, which is typically associated with the absence of chlorophyll depletion [40] and, mostly, stomatal limitation of photosynthesis [41].
Agronomy 2020, 10, x FOR PEER REVIEW 6 of 14 associated with the absence of chlorophyll depletion [40] and, mostly, stomatal limitation of photosynthesis [41]. On the other hand, the analyses of yield components after the harvest ( Figure 2) indicated a decrease of grain yield by almost 50% on average. On the other hand, the analyses of yield components after the harvest ( Figure 2) indicated a decrease of grain yield by almost 50% on average. This indicates that, despite the moderate symptoms, the drought stress experienced was significant, which can be associated with a high economic impact. As the irrigation was initiated in a period when the number of plants and spikes were similar in both variants, the yield decrease was caused mostly by the decrease of grain number and grain weight (which was, however, less affected). The observed effects of the progressive moderate drought on yield components were fully in accordance with the results of other studies [42][43][44]. Based on the results presented above, we can conclude that we observed a drought stress scenario very typical for the temperate climate of Central Europe [36], with significant yield loss but without severe physiological symptoms observed at the leaf or canopy level. In such conditions, the assessment of simple traits such as leaf chlorophyll content or visual scoring may be insufficient in the breeding process. Therefore, there is an urgent need for fast and reliable methods to recognize the physiological effects of drought stress in a collection of genotypes in the field. parameters were calculated using the records measured in wheat leaves in two variants of irrigation. The points represent the mean values ± standard error of 25 genotypes, 10 individuals per genotype (n = 250). An asterisk indicates significant differences between non-irrigated and irrigated variants (Tukey HSD post hoc tests, p < 0.05).
On the other hand, the analyses of yield components after the harvest ( Figure 2) indicated a decrease of grain yield by almost 50% on average. As a reference, we tested the trends of the NDVI parameter ( Figure 3A). This indicates that, despite the moderate symptoms, the drought stress experienced was significant, which can be associated with a high economic impact. As the irrigation was initiated in a period when the number of plants and spikes were similar in both variants, the yield decrease was caused mostly by the decrease of grain number and grain weight (which was, however, less affected). The observed effects of the progressive moderate drought on yield components were fully in accordance with the results of other studies [42][43][44]. Based on the results presented above, we can conclude that we observed a drought stress scenario very typical for the temperate climate of Central Europe [36], with significant yield loss but without severe physiological symptoms observed at the leaf or canopy level. In such conditions, the assessment of simple traits such as leaf chlorophyll content or visual scoring may be insufficient in the breeding process. Therefore, there is an urgent need for fast and reliable methods to recognize the physiological effects of drought stress in a collection of genotypes in the field.
As a reference, we tested the trends of the NDVI parameter ( Figure 3A). This parameter represents the main and most frequently used spectral reflectance indicator, with broad application in both remote and proximal sensing [45]. However, the NDVI values followed the trends observed in leaf chlorophyll or canopy absorbance and significant effects under irrigation were observed only in the very late growth stage of the wheat, which may reflect the early onset of leaf senescence in water-limited plants.
The delayed response of NDVI to the onset of drought has been demonstrated in different species, as reported by several authors [46][47][48]. For example, Liu et al. [48] reported the delay of NDVI response by 18-20 days in grasses. This was in accordance with our study, in which the NDVI parameter was not suitable to detect the early effects of irrigation and differences between the irrigated and non-irrigated variants under the current drought stress scenario.
The second and the main method tested in our experiments was the analysis of fast fluorescence kinetics with analysis of the fluorescence transient using the model of Strasser [32,49]. Unlike the spectral reflectance, the fluorescence transient reflects the dynamic changes in the structure and functions of PSII photochemistry [50]; thus, providing completely different information, compared to methods based on the leaf absorbance or reflectance.
The fluorescence emitted by the chlorophyll molecules of the PSII antenna can serve as an intrinsic probe of the process of performing the successive steps of excitation energy use. The different steps and phases of the chlorophyll a fluorescence rise (O-J, J-I, and I-P) can be linked with the efficiencies of electron transfer in PSII, PSI, and between the two photosystems [49].
Continuous measurement of fast fluorescence kinetics produces two types of parameters: (1) The universal parameters provided by other fluorescence techniques (e.g., the PAM technique), such as F o , F m , and F v /F m ; and (2) model-specific parameters such as specific quantum yields, efficiencies, probability parameters, and performance indices [18]. Parameters such as F v /F m have been broadly used and proved to be sensitive to very severe drought [51] or whether the drought stress is accompanied with heat stress [52]. In our experiments, we observed some fluctuations in the F v /F m values, which were attributed to differences in the actual weather during measurements. On the other hand, a similar F v /F m was found in irrigated and non-irrigated variants. This is in accordance with the previous results, in which a very low sensitivity of the F v /F m parameter was found under moderate drought stress conditions [24,53]. The trend of the F 0 parameter was slightly different. We observed a small, but consistent, increase of F 0 in the drought-stressed variant. A significant but different trend was observed by Roostaei et al. [28], who observed a lower value in drought-stressed plants. As the fluorescence signal F 0 may depend on the chlorophyll content, the decrease may be associated with a loss of chlorophyll as a result of the drought stress, which was not the case in our experiment. In turn, the higher values of F 0 found in our study may be caused by a slower relaxation of basal fluorescence in drought-stressed plants, which is more typical of heat stress conditions [54].
The group of specific "JIP-test" parameters can be also divided into two groups: (1) parameters with a specific physiological meaning; and (2) integrative parameters, represented mostly by performance indices. The parameter area ( Figure 3D), related to the total pool of electron carriers at the PSII acceptor side [49], was found to be quite sensitive to drought stress; however, we also identified differences between two variants before irrigation was applied, which indicates the high sensitivity of this parameter to various factors, which limits its use. A similar situation was also observed in the case of the parameter related to the number of reaction centers, RC/ABS ( Figure 3E). Moreover, the differences between the two variants were relatively small.
The parameter V K /V J reflects the functional state of the oxygen-evolving complex, where an increase of this parameter indicates the damage of a fraction of the oxygen evolving complex (OEC) [23,55]. However, as expected, we observed only a small fluctuation of this parameter, mostly non-significant differences between the irrigated and non-irrigated variants.
In case of biophysically well-defined parameters, the most significant effects of irrigation were found in ψ REo , which represents the efficiency/probability with which a PSII-trapped electron is transferred to PSI acceptors [50]. This parameter is estimated as an amplitude of the I-P step and, in some works, has been denoted by ∆V I-P . The parameter has been found to be sensitive to ozone [56,57], nitrogen deficit [58], and low light [21,56], but also to water deficit [22,59] and combined heat and drought stress [60,61]. The response at the level of the I-P step under irrigation was clearly evident ( Figure 3G) and seems to be a specific stress-related response. Based on our previous results, the I-P step can be interpreted as the relative size of the pools of final PSI electron acceptors [62] and the change in I-P amplitude can be associated with the partitioning of part of PSI RCs into the cyclic electron flow, as an acclimation response to drought stress [21].
Performance indices are proposed to combine information on the status and activity of PSII and the efficiencies of specific electron transport reactions in the thylakoid membrane during the O-J-I-P rise, altogether providing a sensitive tool containing useful information beyond that provided by F v /F m , which has been broadly used to assess the effects of stress [63]. Most frequently used is the performance index calculated on absorption basis (P Iabs ), combining information on the number of reaction centers, efficiency of energy trapping, and electron transport from PSII to the plastoquinone pool [32]. More recently, the total performance index has been introduced [64], which contains also the efficiency with which an electron from the plastoquinone pool is transferred to the final PSI acceptors.
This study confirmed that PI tot and PI abs are very sensitive indicators of the physiological status of wheat under field conditions; however, the effect of drought stress was much more significant in case of the total performance index (PI tot ). The changes that occurred in both performance indices may be interpreted as evidence for the considerable modulation of PSII function during the season; however, in our specific drought stress scenario, limitation of the electron transport chain at the PSI level represented the main effects determining the statistical significance of the differences in PI tot between non-irrigated and irrigated variants. The idea to preferably use PI tot was also proposed by Stirbet et al. [63], as this parameter is related to the function of the "whole" linear electron transport, whereas PI abs is related only to electron transport to the plastoquinone (PQ) pool. Moreover, it has been shown that PI tot is less sensitive to diurnal changes and mid-day depression, compared to PI abs [21], which is important for practical applications, as the field measurements of a high number of genotypes can hardly be done within a narrow time window.
Considering the possible practical applications of the method, in addition to the sensitivity of response of the parameters, the differences in responses of individual genotypes need to be recognized. As the presentation of the values of multiple parameters in numerous genotypes and individual dates would be difficult and inefficient here, we address these issues through a summary statistical analysis of the significance of the factor genotype, in addition to analysis of the treatment effects (Table 3).
We also tested the interactions between the genotypes and treatment (irrigation). In the first step, the analysis was focused on the statistical significance of the differences between the two blocks in the period before the irrigation was initiated, in which no difference was expected. However, as blocks in the field are never completely identical, we found some significant effects in the parameters AREA and RC/ABS. This indicates that these parameters are very sensitive to minor side effects, which are difficult to explain. The sensitivity of these parameters to various environmental factors has been previously identified, including sensitivity to minor diurnal or developmental changes [18,22,26]. It has also been found that the number of active reaction centers (RC/ABS) changes during leaf development [65]. Similarly, a large seasonal variation in the fluorescence transient, with an apparent effect on the area above the fluorescence curve (indicated by the parameter AREA) has been found in cereal crop species [66]. In some other parameters (e.g., F v /F m and V k /V J ), developmental variation was found only under conditions of severe leaf senescence [67]. We could only guess that there were coincidences in the minor factors (e.g., microclimate, soil temperature, and so on) which led to the slightly faster development of new leaves in one of the blocks in early spring, which was then detected by the sensitive fluorescence parameters (i.e., RC/ABS and AREA) during the early spring measurements. Such sensitivity is not desirable and hence, the suitability of these overly sensitive parameters for the assessment of drought effects in the next period was considered limited. The second two-way analysis of variance assessed the significance of the factor "irrigation" and "genotype" in the period of three weeks after irrigation was initiated; that is, we assessed the statistical significance of the early effects. The statistical significance of the differences between two variants in three measurement cycles after irrigation were initiated by identifying the significant results in five parameters (PI tot , ψ REo , F 0 , AREA, and RC/ABS), whereas we did not find any statistical significance in three parameters (PI abs , F v /F m , and NDVI); in these parameters; therefore, we did not assess the effect of genotype and interaction.
As all the parameters were assessed in the same design and number of repetitions, we could also compare the power of the factors and interactions by comparing the F values. Although all presented effects in the five parameters were statistically significant, the power of the factors was not the same. We expect that a suitable parameter should be sensitive to drought/irrigation and able to recognize differences between the genotypes. A high interaction level for the parameter is also very convenient, as it means that the genotype differences in response to drought are shown by the values of this parameter. In this respect, the best values were found in the PI tot and AREA parameters. However, the significant difference found in the period before irrigation excluded the parameter AREA. A high F-value for interaction between the factor genotype and irrigation indicated that the patterns of AREA and, especially, PI tot were not fully identical but, in some genotypes, the differences were higher than in others. In the parameter F 0 , we identified a high F-value for the factors irrigation and genotype; however, the interaction was only moderate. This may be associated with the relatively low variation of F 0 values between the variants observed in our study. Similarly, a low relative variation between the variants, together with the high variation between the dates observed in ψ REo , might provide a reason for the lower F-values found in this parameter.
Thus, our statistical examination confirmed that the analysis of fast fluorescence kinetics may be used to recognize the early effects of moderate drought stress and irrigation on the photosynthetic apparatus of wheat. The integrative parameter performance index, PI tot , was identified as the most suitable parameter for the comparison of genotype-related differences in drought response.

Conclusions
Analyses of the time-series of non-invasive measurements obtained in a dry season within a representative collection of wheat genotypes grown in small-plot field trials in non-irrigated and irrigated variants showed that, despite a progressive water deficit and significant yield loss, the measurements of leaf and canopy absorbance showed very minor changes in chlorophyll content or canopy structure in the non-irrigated variant, compared to the irrigated variant. This corresponded well to the insignificant differences in spectral reflectance NDVI values. On the other hand, we identified the significant and rapid response of fast fluorescence kinetics indicators following the onset of irrigation. Analysis of multiple fluorescence parameters indicated that the main effects of drought/irrigation were associated with changes in amplitude of the I-P phase in the OJIP transient, which documents the main changes of the linear electron transport chain at the level of photosystem I and beyond. The statistical analyses identified the integrative parameter performance index PI tot , containing the information on the I-P amplitude as the most sensitive parameter, as well-reflecting the differences in responses of the considered genotypes to water deficit, which was not the case when using conventional parameters (e.g., F v /F m ) or many kinetic parameters. Our results demonstrated that focusing on photosynthetic functions, as detected by rapid chlorophyll fluorescence records, can provide more accurate information on the drought stress level, compared to the structural data obtained by absorbance or reflectance measurements.