Assessing the Potential of Downscaled Far Red Solar-Induced Chlorophyll Fluorescence from the Canopy to Leaf Level for Drought Monitoring in Winter Wheat

: Solar-induced chlorophyll ﬂuorescence (SIF) from ground, airborne, and satellite-based observations has been increasingly used in drought monitoring recently due to its close relationship with photosynthesis. SIF emissions respond rapidly to droughts, relative to the widely used vegetation indices (VIs), thus indicating their potential for early drought monitoring. The response of SIF to droughts can be attributed to the confounding effects of both the physiology and canopy structure. In order to reduce the reabsorption and scattering effects, the total emitted SIF (SIFtot) was proposed and served as a better tool to estimate GPP compared with the top-of-canopy SIF (SIFtoc). However, the response time and response magnitude of SIFtot to droughts and its relationships with the environmental parameters and soil moisture (SM) (i.e., the knowledge of drought monitoring using SIFtot) remains unclear. Here, the continuous ground data of F760toc (SIFtoc at 760 nm) from a nadir view that was downscaled to F760tot (SIFtot at 760 nm), NIR v , and the NDVI, SM, meteorological, and crop growth parameters were measured from four winter wheat plots with different intensities of drought (well-watered, moderate drought, severe drought, and extreme drought) over 2 months. The results indicated that F760tot was more closely correlated with the SM than the VIs at short time lags but weaker at longer time lags. The daily mean values of F760tot and NIR v were able to distinguish the differences between different drought levels, and F760tot responded quickly to the onset of drought, especially for the moderate drought intensity. These ﬁndings demonstrated that F760tot has potential for early drought monitoring and may contribute to mitigating the risk of agricultural drought.


Introduction
Droughts are one of the most extreme and least understood climatic events in many regions of the world that has been responsible for significant socioeconomic and ecological consequences [1]. Climate change increases the global frequency and severity of droughts, thus exacerbating yield uncertainty and posing several risks to agriculture [2]. Consequently, there is an urgent need for obtaining a better understanding of drought detection. Historically, drought monitoring approaches have evolved from in situ station-based measurements (e.g., the Palmer Drought Severity Index (PDSI)) toward the continuous observations and monitoring of key drought-related variables at regional and global scales using remote sensing [3]. Among them, vegetation indices (VIs) have been widely used for estimating drought severity across space and serve as good indicators to monitor drought [4]. However, it is noted that greenness-based VIs are incapable of capturing dynamic changes of drought responses in shorter timescales, since these indices relate to canopy greenness directly and potentially decouple from photosynthetic functioning [5].
Solar-induced fluorescence (SIF) is an optical signal (650-850 nm) emitted from plant chlorophyll and has been used as an indicator of the actual functional state of plant photosynthesis. Recent research has demonstrated that SIF is promising for estimating gross primary production (GPP) [6], monitoring water stress [7], tracking phenological changes [8], predicting crop yield [9], and estimating leaf nitrogen content [10], among other uses. With respect to stress detection, various studies have explored the impacts of droughts on the red (e.g., F685) or far red (e.g., F760) SIF observed by the ground, airborne, and satellite platforms. For example, Pérez-Priego et al. [11] conducted the field experiments to acquire ratios between the 757 nm (out) and 760 nm (in) bands and revealed the good relationship between fluorescence and droughts at the canopy scale. The far-red SIF was more sensitive to VPD [7] and the short-term standard precipitation index (SPI) than the Normalized Difference Vegetation Index (NDVI) [12]. More importantly, the site-observed F760 was significantly correlated with the root zone soil moisture (SM) [13], and it outperformed the VIs for early drought detection in a short time lag, especially over a closed canopy [14].
However, the top-of-canopy SIF (SIFtoc) is only a fraction of the emitted SIF photons due to the absorption and scattering effects of SIF emission in the canopy scale, thus influencing its relationship with GPP and showing the confounding effects of the leaf optical properties and canopy structures on SIF [15]. Therefore, many efforts have been made to reduce the influence of normalizing the solar zenith angle (SZA) of SIF observations or downscaling SIF from the canopy to leaf level. Some studies [16][17][18] have used the bidirectional reflectance distribution function model to correct the angular effect of SIF and demonstrated a stronger relationship between SIF and GPP. Alternatively, the importance of downscaling SIF from the canopy to leaf level has attracted increasing attention recently. Generally, at the canopy scale, the amount of SIFtoc is frequently defined as SIF toc = PAR·fPAR·Φ F ·f esc (1) where PAR represents the photosynthetically active radiation received, fPAR is the fraction of PAR absorbed by vegetation, Φ F is the fluorescence efficiency, and f esc is the fluorescence escape ratio (i.e., the probability of a photon escaping from the canopy in the direction of the sensor). Based on SIFtoc and f esc , the term of SIFtot (total canopy SIF emission) can be defined as SIF toc = SIF tot ·f esc (2) The core of downscaling SIF from the canopy to leaf level is to estimate f esc . For example, a model for retrieving the spectral shape of the leaf fluorescence spectrum was proposed by Romero et al. [19]. Liu et al. [20] utilized the Random Forest (RF) model to estimate SIF at the photosystem level in both the red and far-red bands and improved the correlations between SIF and absorbed photosynthetically active radiation (APAR). A simple approach to estimate the f esc of SIF at the far-red band was introduced based on the near-infrared reflectance of vegetation (NIR v ) and the fraction of absorbed photosynthetically active radiation (fPAR) [21]. Generally, the O 2 -A band (~760 nm) is most commonly used for retrieving SIF results, or in other words, F760 (SIF at 760 nm). However, the fact is that only a fraction of the emitted SIF photons is observed from the top of the canopy, which is therefore called F760toc. Using the above method, F760toc can be used to calculate the total emitted SIF (SIFtot) at 760 nm (i.e., F760tot), thus reflecting the total emitted SIF from the entire canopy. In this regard, some studies [22] suggested that SIFtot exhibited better performance in predicting GPP compared with SIFtoc. In contrast to these findings, Dechant et al. [23] actually found a decreased performance of F760tot compared with F760toc for GPP estimation.
With respect to stress detection, SM is used to indicate drought-induced water stress, and the relationship between SIF and SM is complex. Generally, the physiological response of plants to a drought includes the stomatal and non-stomatal limitations, thus leading to the variation of photosynthetic electron transport, non-photochemical quenching (NPQ), and SIF. As the SM continuously declines until a certain threshold, the stomatal and nonstomatal limitations lead to a decrease both in photosynthesis and fluorescence. Based on many canopy-scale studies [5,13,14], a decrease in F760 can be observed under water stress, and F760 has a higher correlation relationship with SM than VI. Specifically, the behavior of the canopy-scale SIF under a drought is the ensemble of the changes in the physiological reaction and canopy structure [23], and among them, Xu et al. [5] suggested that structural factors dominated the spatial response of SIF to water stress over the physiological reaction, while F760tot represents the total emitted SIF at the leaf level and should be less affected by canopy structures. Therefore, it is unreasonable to speculate on the performance of F760tot for drought monitoring based on the performance of F760toc. That aside, only a few studies, such as Liu et al. [24] demonstrated that the value of SIFtot apparently declined compared with the multi-year average under a drought, and there appeared to be similar monitoring effects for both SIFtoc and SIFtot regarding droughts. This study used satellite SIF data to obtain these findings without ground experiments. Moreover, the response time and magnitude of SIFtot (F760tot) for droughts and its relationships with the environmental parameters and SM (i.e., the potential of F760tot in drought detection) remain unclear. If this is figured out, the results might provide a deep insight into the utility of F760tot in drought monitoring, thus contributing to mitigating the risk of agricultural droughts.
To explore the monitoring effects of F760tot on droughts, here, the continuous ground measurements to four plots under different levels of drought for a wheat crop were carried out. An automated field spectroscopy system was used to collect the time series of canopy F760toc, and then F760toc was downscaled to F760tot for analyzing its relationship with the environmental and crop growth parameters and addressing the following three objectives: (1) explore the detailed responses between F760tot with F760toc and VIs in responding to different intensities of drought; (2) reveal the relationships between PAR and the growth parameters with seasonal F760tot; (3) determine the relationships of F760tot, F760toc, and VIs with SM across the growth season of wheat.

Ground Measurements
A dataset comprising the ground spectral, canopy structure, and vegetation physiological measurements of winter wheat, as well as meteorological and SM measurements, was acquired at the Fangshan Comprehensive Experimental Station (39 • 35 N, 115 • 42.5 E) of Beijing Normal University, following the experimental settings of different drought levels for four plots. It was used to evaluate the performance of F760tot for drought monitoring.
Winter wheat was sown in October 2016 in rows in four plots (4 m × 4 m) and harvested in June 2017. To protect the wheat from freezing, the same irrigation was conducted on four plots before the green-up stage. After the green-up, the four different irrigation treatments in Table 1 were applied in order to generate four drought levels: wellwatered treatment, moderate drought, severe drought, and extreme drought, corresponding to the order of P1, P2, P3, and P4. Among the four plots, only the irrigations of P1 met the needs of the winter wheat during growth. The winter wheat in other plots suffered different levels of water stress. During the study period, there was no rainfall, and therefore the experiment was not affected by precipitation. Eventually, the prolonged drought levels of four plots were generated, and the period between 28 March 2017 and 19 May 2017 was selected as the time window for field measurements. Details on the measurements including the parameters, methods or devices, and intervals are summarized in Table 2 as follows. The spectral measurements were conducted using an automated field spectroscopy system [14], which was improved based on the FluoSpec system [25]. As displayed in Figure 1, the core components of this system were an Ocean Optics HR2000+ spectrometer (wavelength: 680-775 nm; FWHM:~0.13 nm) and a fiber switch (OceanOptics, Inc., Dunedin, FL, USA) to measure the up-and downwelling radiation fluxes in sequence (nadir view,~4.5 m height, hemispherical-conical configurations, footprint: 0.45 m). One of the fibers which pointed toward the wheat canopy to collect upwelling radiance was a bare fiber (FOV: 25 • ), and the other fiber attached to a cosine corrector (CC-3, OceanOptics, Inc., Dunedin, FL, USA) was pointed toward the sky to collect incident light (FOV: 180 • ). In order to reduce the influence of dark current, the HR2000+ spectrometer was placed in an incubator with a fixed temperature of 23 • C.
Another key component was the stepping motor with a round trip that could be rotated over the 4 plots in 6-min intervals. Moreover, the surface soil moisture at the 20-cm soil layer, incident photosynthetically active radiation (PAR, 400-700 nm), and meteorological parameters were collected automatically using different devices. Additionally, some critical growth parameters of wheat such as the leaf area index (LAI), chlorophyll content (Chl), and relative water content of leaf (RWC) were manually measured in a coarser temporal resolution. Specifically, the growth parameters were measured six times during the field campaign.

F760toc Retrieval
Before SIF retrieval, pre-processing, including dark current correction, a Savitzky-Golay (SG) filter, radiometric and wavelength calibration, and quality control were conducted for solar irradiance and canopy radiance, which were then converted to irradiance (mW/m 2 /nm) and radiance (mW/m 2 /sr/nm). Because of the limited height between the sensor and canopy surface (<2 m), the effect of atmospheric radiation was negligible and therefore not corrected in this study.
Various methods have been developed to decouple SIF from the total at-sensor radiance from spectral measurements. Among them, major SIF retrieval methods were derived from the Fraunhofer line discrimination (FLD) principle [26]. For example, several different FLD-based algorithms including the three-band FLD (3FLD), the improved FLD (iFLD), and the pFLD method have been developed, focusing on the telluric O 2 absorption bands or solar Fraunhofer lines. These methods have proven to be robust and simple to apply but may overestimate the values of SIF emissions [27]. Hence, by assuming that the SIF and reflectance were polynomials or other appropriate mathematical functions of wavelengths, the Spectral Fitting Methods (SFM), proposed by Meroni et al. [28], utilized a set of contiguous channels to estimate SIF. Apart from these FLD-based algorithms, statistical methods based on singular vector decomposition (SVD) or principal component analysis (PCA) were proposed to enhance FLD SIF retrieval from satellites [29,30] but were applied to extract ground-based SIF [31].
Here, only F760toc was retrieved due to the limitation of our spectrometer. According to the analysis by Chang et al. [32], the FLD-based retrievals of SIF were noisier and more sensitive when assuming the SIF spectral shape than the O 2A -based method, and the authors recommended that the SVD or SFM method using a reduced fitting window (759.5-761.5 nm) was robust for far-red SIF retrievals across different sky conditions. Therefore, based on this reduced fitting window, the SFM algorithm was selected to estimate F760toc at the canopy level. The basic equation was where L(λ) is the radiance upwelling from the vegetation canopy at ground level, E(λ) is the total solar irradiance incident on the target, r MOD (λ) and F MOD (λ) are the mathematical functions used to describe and model (MOD) the two key variables, and they are all distributed at each wavelength.

Estimation of fPARchl
According to Equations (1) and (2) in Section 1, fPAR was essential for the escape probability to calculate SIFtot, which was not measured directly in the field experiment. Therefore, the statistical model based on regression analysis between the fPAR and VIs was selected as the method to estimate fPAR. Pioneering works [33,34] have shown that fPAR (i.e., fPAR canopy or fPAR chl ) could be approximated with the NDVI, enhanced vegetation index (EVI), wide dynamic range vegetation index (WDRVI), and so on. Here, the SCOPE 2.0 model [35] was utilized to illustrate the relationships between different fPAR values and VIs.
Details about the SCOPE input parameters are listed in Table 3 (i.e., different leaf chlorophyll contents (Cab) and LAI levels to cover the most common vegetation conditions of winter wheat, as well as different soil moisture contents (SMCs) and SZA to represent different drought levels and solar-view geometries, respectively). Consequently, 10,530 different samples of spectral reflectance (1-nm intervals), aPAR, and aPAR byCab were generated. Different VIs containing the Normalized Difference Vegetation Index (NDVI), Red Edge NDVI (NDVI 705 ), WDRVI, and Optimized Soil-Adjusted Vegetation Index (OS-AVI) are calculated in Table 4, where the average values of 684-686 nm, 704-706 nm, and 749-751 nm were used for the Red, Red edge , and NIR bands, respectively. As depicted in Figure 2, the fPAR chl estimation of NDVI 705 outperformed the original counterparts with the highest R 2 (0.921). The linear fPAR chl -NDVI 705 function was expressed as: The PAR above the canopy (PAR above ) was directly measured using three quantum sensors (S-LIA-M003, Onset Computer Corporation, Bourne, MA, USA). Consequently, the APAR chl could be calculated as follows:

F760tot Calculation
SIFtoc may be downscaled to SIFtot of all leaves or all photosystems by considering the absorption and scattering effects of SIF emissions within the canopy or the leaf [15], respectively. Based on the concept of recollision probability (or 'p-theory') [40], several studies aimed at describing the radiative transfer of SIF emissions within the canopy and the SIF escape probability (f esc ) from the leaf to canopy level. For example, Yang and van der Tol [41] and Liu et al. [20] showed that f esc could be modeled as follows: where δ FC (λ) is the canopy scattering of SIF, BRF(λ, Ω) is the directional canopy reflectance for the direction Ω, i 0 is the proportion of incident photons intercepted by the canopy, and w is the leaf-scattering coefficient or leaf single-scattering albedo. However, these parameters cannot, in general, be accurately estimated, so different available variables were chosen to represent the components of δ FC (λ), i 0 , and w. For BRF(λ, Ω), it was usually assumed to be consistent with NIR v = NDVI * NIR [21], R Veg = R Obs −R Soil [22], or Red V = Red tot ×NDVI 2 [42]. For i 0 , it appeared to play a similar role to fPAR [42] or be calculated as 1 − EXP(−0 .5 × L e /cos θ S ) [22]. For w, it was expected to be relatively stable (e.g., 0.853~0.888) at the far red band and 0.287~0.444 at the red band, as suggested by Liu et al. [20]. In this study, NIR v for BRF, fPAR for i 0 , and 0.87 for w were combined to calculate F760tot: Based on f esc , F760 tot can be calculated as:

Data Analysis
The pre-processing of different measurements was introduced in the previous sections, and data analyses were then performed to address the scientific questions raised previously. For all different measurement data, the SIF data were dominant, and only the time steps of the meteorological data and PAR, SM, and vegetation parameters that corresponded to the SIF data were selected in order to ensure direct comparability. The ratio of F760/PAR (referred to as AF760toc and AF760tot hereafter) was also calculated to remove the influence of solar irradiance.
The temporal resolutions of the measurements are displayed in Table 2. The hourly and daily mean values (from 8:00 a.m. to 5:00 p.m.) of the different parameters were calculated. For the hourly scale, the hourly mean F760 and F760/PAR as well as the meteorological data were calculated along with the PAR. The fPAR derived from the linear regression and APAR was also averaged hourly to calculate F760tot. For the daily scale, the paired daily data mentioned above were also calculated depending on the average hourly mean values.
First, the utility of F760tot for drought monitoring in comparison with F760toc, F760/PAR, and VIs (NDVI, NIR v ) was evaluated during most of the growing stages. Secondly, to reveal the effect of different parameters on F760tot under a long-running drought, correlation analyses of F760tot with PAR and other vegetation parameters such as LAI and Chl were conducted for the four plots. Third, the goal of this part is to discuss the detailed differences in the relationship of the root zone soil moisture (20 cm SM) with F760tot, AF760, and the VIs. Furthermore, the lag responses of 20 cm SM to F760tot were also explored with the development of wheat. All correlation analyses were performed based on the Pearson correlation coefficient, and the statistical significance of the Pearson correlation coefficient was evaluated with a two-sided t-test at a confidence level of 95%.

Differences of F760tot, F760toc, and VIs in Responding to Different Drought Levels
In Section 2.1, a dataset containing different types of measurements was obtained for subsequent analysis. As illustrated in Figure 3, there was a steady increase in these three parameters, including the VPD, Ta, and SR, during the growing season, among which Ta appeared to have the greatest increment ( Figure 3b). Specifically, these three parameters showed the daily variability, and the peaks of Ta were well correlated with those of the VPD   In order to discuss the responses of F760tot to different drought levels over four plots, here, the seasonal patterns of F760tot in comparison with F760toc, the NDVI, and NIR v were analyzed to illustrate the detailed relationships of these parameters with drought development. As shown in Figure 4, the drought levels of the four plots (P1 > P2 > P3 > P4 from wet to drought) was successfully established, thus supporting exploring the utility of F760tot for drought monitoring according to Liu et al. [14]. It was clearly observed that the F760tot, F760toc, and VIs showed different seasonal patterns; that is, the mean values of F760tot as well as F760toc were generally consistent with a trend of going up slightly first (27 March to 24 April) and then going down slightly (24 April to 20 May), and some valleys for F760tot and F760toc were perfectly matched with the low values of SR as mentioned above (Figure 4a-c), indicating the potential links between them. The VI values were characterized as a smoothly declining trend during the growing season ( Figure 4d,e). The different seasonal patterns suggested that the SIF emissions of the wheat were susceptible to the influences of various factors compared with the VIs.
Given the different drought levels, there were similarities and differences between the performance of F760tot as well as F760toc and the VIs in the four plots. As for the similarities, the values of all variables of the four plots exhibited a difference in gradient consistent with the drought level (P1 > P2 > P3 > P4) in the growing season after irrigation, such as in the period of the second rectangle (11)(12)(13)(14)(15)(16)(17)(18)(19)(20). These results indicate that all variables had the potential to respond to the irrigation-induced variations in SM during sustained and prolonged droughts. However, on the contrary, the F760tot values of P4 were higher than its values for P1 in the period of the second rectangle (Figure 4b), which may have been related to premature yellowing of the wheat leaves during long-term droughts. As for the differences, these mainly involved that the response time of F760tot, F760toc, and NIR v to droughts were obviously ahead of the NDVI in the period of the first rectangle (18)(19)(20)(21)(22)(23)(24)(25)(26)(27)(28)(29), thus indicating the rapid response of F760tot, F760toc, and NIR v to the onset of a drought, and among them, the response patterns of F760toc and F760tot to a drought were quite similar (i.e., responding simultaneously to the second irrigation (22 April) on a daily scale).
The newly proposed NIR v index appeared to respond more rapidly to droughts than the NDVI, but the magnitudes of its values of P1 and P2 were opposed to the different drought levels between them, and the order of NIR v values matched with the drought levels until the third irrigation. That aside, a rapid response difference between F760tot, F760toc, and the NDVI for the irrigation-induced variation of SM was also observed. For instance, a sudden rise happened with F760tot and F760toc of the four plots after irrigation, especially in the second irrigation (22 April to 20 May), while the values of the NDVI showed a smooth trend, indicating an accumulative effect of the drought. Conclusively, it could be inferred that both F760tot and F760toc had an edge over the VIs in the early detection of droughts.
Additionally, the seasonal patterns of AF76toc and AF76tot (SIF/PAR, normalized by PAR) under different intensities of droughts were also analyzed. Overall, SIF/PAR showed similar trends of change with those of SIF ( Figure 5); that is, its values for the four plots were generally directly proportional to the drought levels. Specifically, as shown in Figure 5a,b, the fluctuation of SIF/PAR was reduced after eliminating the effect of PAR on SIF. In the first rectangle of two irrigations, the AF760 values of the four plots were displayed in the order of P2 > P1 > P3 > P4 for AF760toc and P1 > P2 > P3 > P4 for AF760tot. In the second rectangle of the four irrigations, the AF760 values of the four plots were shown in the order of P1 > P2 > P3 > P4 for AF760toc and P1 > P4 > P2 > P3 for AF760tot. In order to further investigate the impact of F760tot on droughts, the drought-induced decreases, especially those caused by SM, in the NDVI, NIR v , F760toc, and F760tot between P1 and the other plots are explored in Figure 6. The theoretical basis for this analysis was that the field experiment took the SM as the control condition, and other properties such as the illumination, temperature, and wheat of the four plots were basically similar. Therefore, the differences for P1 and the other plots directly represented the impact of SM droughts on different variables. As shown in Figure 6, the drought-induced decreases of the four variables gradually increased along with the aggravation of the drought, ranging from 40% to 90%. In the same intensities of droughts, the drought-induced decreases of the four variables also gradually increased as the drought progressed, except for F760tot in an extreme drought. As for the NDVI, it exhibited a later response to the onset of a drought than other variables and decreased more dramatically for the long-term drought than SIF, particularly in extreme droughts. As for NIR v , it seemed to have a compromise in performance between the SIF and NDVI in early drought monitoring, and among these four variables, the drought-induced decreases of NIR v and the NDVI were gradually higher than those of SIF, whereas in a moderate drought, the NIR v values of P2 were conversely higher than those of P1. Given the SIF variables, the F760tot values decreased more than other variables when in response to the onset of a moderate drought, and the decreases of F760tot were slightly less than those of F760toc.

Relationships between PAR and Growth Parameters with Seasonal F760tot
SIF was increasingly used for crop stress monitoring, but it remained challenging to reveal the effects of the environmental and growth parameters on SIF emissions during long-running droughts. For example, more and more researchers have proven that the PAR is the primary environmental parameter driving SIF variations at different time scales (Yang et al., 2015). As shown in Figure 7, a relatively strong correlation between F760toc and the PAR of four plots (r value = 0.736 for P1) was found, and the correlation coefficients declined with the drought levels (P1-P4, r value: 0.7~0.4). After being obtained by normalizing F760toc for f esc , F760tot appeared to have a stronger correlation relationship with the PAR (r value:~0.8 for P1), and compared with F760toc, there was only a minimal drop of the r value from 0.799 (P1) to 0.75 (P4), which seemed to suffer fewer effects from the drought for the relationship between F760tot and the PAR. In terms of the APAR, the better correlations between them of P1-P3 were found, whereas in P4, F760tot and the APAR exhibited a weak positive relationship. Figure 7. Relationship between hourly PAR, APAR, and F760toc (blue) or F760tot (golden) in the field campaign collected from four plots (P1-P4). The blue or golden lines with band denote the regression line and 95% confidence interval, respectively. Significance was performed according to two-tailed test (ns: p > 0.05. * p < 0.05. ** p < 0.01. *** p < 0.001).
Apart from the PAR, LAI and Chl are two growth parameters that will dramatically reduce when crops experience water stress, and these two parameters also influence the emission of SIF. Therefore, partial correlation analysis was then applied to explore the relationship between the SIF and VIs with LAI and Chl in Figure 8. Here, the LAI and Chl measurements seemed to be insufficient for rigorous analysis, but it is noted that our analyses were considered in the global perspective and focused on the results that had been tested for significance.
First, the average correlation of LAI with SIF (r values > 0.6) was greater than that of Chl with a significant positive correlation, indicating that LAI could have had much more of an impact on SIF emissions than Chl over the growing season under different drought levels. Then, as depicted in Figure 8a for LAI, the impact of LAI on F760toc usually exceeded that of F760tot from P1 to P3. This suggested that SIF was less prone to vegetation canopy structures after downscaling to the leaf level, whereas for P4 in an extreme drought, a similar impact from LAI was observed on F760toc and F760tot. As with AF760toc and AF760tot, the former was more related to LAI from P1 to P2, and the latter was more linked to LAI from P1 to P2. Among them, AF760toc showed the most significant correlation with LAI (r value: 0.92, p < 0.001). As depicted in Figure 8b for Chl, F760tot was more related to Chl than F760toc from P1 to P3, indicating that the former was closely integrated with the plant physiology. F760tot showed the most significant correlation with Chl, whereas the impact of Chl on F760toc was less than that of F760tot for P4 under an extreme drought. As for AF760toc and AF760tot, the impact of Chl on AF760toc was usually superior to that of AF760tot for P1, P2, and P4 instead of P3. Figure 8. Partial correlation coefficients (r) of (a,b) with F760toc, F760tot, AF760toc, AF760tot, NDVI, and NIR v based on daily measurements. Significance was performed according to two-tailed test (ns: p > 0.05. * p < 0.05. ** p < 0.01. *** p < 0.001).

Relationships of Root Zone Soil Moisture with F760tot
As mentioned in the previous section, the PAR is the common driving force for the dynamics of SIF, whereas in Figure 7, such a driving effect declined for four plots along with the drought levels, thus indicating the involvement of SM in SIF variations. Pearson analysis was then applied to further investigate the relationships between F760tot with SM in different time scales (Figure 9). On the daily and hourly scales, it was obvious that the impacts of SM on F760tot and the VIs were becoming close, along with an increase in drought stress. For P1 with no drought, SM seemed to have little influence on F760tot and the VIs. As the drought intensified from P2 to P4, the response of F760tot and the VIs tended to be different. The correlations between the NDVI and NIR v (r values: 0.68 and 0.75) with SM were positively much greater than those with SIF. As for SIF, only AF760toc exhibited the most positive relationship with SM (r values: 0.54 for P4), and the reason for this may lie in that there were strong correlation relationships between AF760toc and the VIs. The previous analysis was conducted using the whole data from the growing season. Furthermore, taking P4 as example, the relationships between F760tot and the VIs with SM in the 20-cm layer under different time windows from 3 days to 28 days ( Figure 10) were explored, aiming to reveal the distribution of correlation coefficients under a drought. There were similarities and differences between the performances of F760tot and the VIs with SM of P4 (extreme drought). It was confirmed that the correlation coefficients of different parameters showed the larger distribution of numerical data in the short time window and became concentrated along with the increase of time steps. Meanwhile, different variables containing SIF, SIF/PAR, and the VIs had their unique patterns. For example, the relationship of F760tot and F760toc with SM was weak, and AF760 and the VIs showed an increasing and positive trend correlation coefficient with SM.
The correlations of SM in the 20-cm layer (SM) with lagged SIF (F760toc and F760tot) as well with lagged VIs (NDVI and NIR v ) were also explored in Figure 11. In order to figure out their relationships under different time windows (from 1 week to 4 weeks), the r values using data from P4 representing an extreme drought were calculated. First, the correlations of SIF-SM and VI-SM presented a steady increasing trend and had little amplitude fluctuation with the increase in the time windows, and the correlations of the latter exceeded that of the former using the data over the growing season, suggesting that the VIs were more applicable for the drought lasting for a long time scale. In the time window of 1 week (Figure 11a), the r values of SIF-SM had a similar pattern to the trigonometric curve, and the r values of VI-SM increased in the range of 8 days and then kept consistent (r values > 0.75). Specifically, the r values of SIF-SM were higher than that of VI-SM when the time lag was~5 days (r values > 0.7), and the latter exceeded that of the former when the time lag was~9 days. This further indicated that SIF outperformed the VIs in early drought monitoring. In other panels, especially in Figure 11c,d, it is obvious that the VIs had a strong relationship with SM in different time lags, while F760tot turned out to decouple with SM in time lags of about 3 or 4 weeks.

Discussions
Here, the responses of F760tot to different intensities of droughts are investigated, and then the relationships between F760tot and other variables with PAR and growth parameters, especially with root zone soil moisture over a seasonal timescale, are explored. The discussions are presented below.

Pros and Cons of F760tot for Drought Monitoring
The first main result lied in that F760tot and F760toc were better than the NDVI at early drought monitoring, which is consistent with previous research such as that by Lee et al. [7] and Liu et al. [14]. On this basis, it was further found that both F760tot and F760toc had the advantage of early crop drought monitoring, and they exhibited similar responses to the onset of a drought (Figures 4b,c, 10 and 11). SIFtot was downscaled from the canopy to leaf level and outperformed SIFtoc at estimating the GPP [15]. In general, it should be less influenced by canopy structures and more closely related to the plant physiology, which was consistent with the results of Figure 8. The reason for the similar response of these two is that it was incapable of capturing the onset of a drought earlier in time compared with SIFtoc as observed in nadir viewing, because the latter was representative among different observation directions due to its effectiveness in removing the directional effects and reducing the uncertainties induced by varying sun-target-sensor geometries [18].
However, although the response times of F760tot and F760toc to droughts were similar, further studies ( Figure 6) demonstrated that F760tot had the potential for detecting moderate droughts because it had more of a decrease than other variables in moderate droughts. Specifically, SIFtot could represent the overall chlorophyll fluorescence emissions, which could reflect more of a decrease of fluorescence emissions under the same level of drought, therefore possibly improving early monitoring for moderate droughts.
Another main result was that the NDVI and NIR v outperformed F760tot in long-term drought monitoring over the growing season. As depicted in Figures 9-11, the r values between the VIs and SM up to~0.75 were higher than that of F760tot and F760toc, regardless of the different time windows or whether the water stress hysteresis effects were considered. Specifically, it was obvious that due to the weakness of F760tot and F760toc signals and the influence of many factors, the values of F760tot and F760toc fluctuated dramatically during the growing season compared with the VIs in Figure 4 and even showed a slow, regular up and down movement in the time window of 1 week in Figure 11a. The values of the NDVI, in contrast, maintained a smooth, continuous curve in the seasonal scale in Figure 4 and tended to rise during the first 8 days and stabilize in the rest (Figure 11a). Additionally, according to Figures 4 and 5 and previous studies [25], although the relationship of F760toc with PAR gradually decreased with the increasing severity of the drought, a pronounced day-to-day variation in the SIF pattern was highly correlated with that of PAR, so the daily accumulated values of PAR varied in accordance with the fluctuation of SIF during the growing season, thus leading to the weak correlations between SIF with SM.
Here, the drought is defined by the moisture deficits of the root zone. More and more studies have proven that the SIF emissions were sensitive to drought stress, which could be attributed to the decrease in carboxylation capacity and the increase in non-photochemical quenching. For example, compared with VPD, SM is the dominant driver of dryness stress on ecosystem production that is characterized by SIF [43]. Similarly, using SIF as a proxy for GPP, Dang et al. [44] demonstrated that SM plays a more important role, and SIF is very sensitive to aridity gradients in arid and semi-arid ecosystems. The results coincided with our reports. Based on this sensitivity, Yu et al. [45] integrated the relationship of the photosynthesis rate and photosystem II operating efficiency into SUCROS97 and improved the prediction of biomass under a drought. Aside from that, Cannière et al. [46] proposed a framework for parametrizing the water stress function using SIF in the AgroC crop growth model and improved the estimates of actual evapotranspiration during the stressed periods.
In terms of the two results, it was noted that F760tot could respond rapidly to droughts, and the NDVI was more sensitive to long-term water deficits. The main reason for this is that SIF is strongly correlated with photosynthesis and serves as an indicator of photosynthetic activity [47], thus prompting SIF to respond near-instantaneously to rapid adjustments in the photosynthetic machinery [48] and highlighting its potential to improve the present drought detection approach [49]. From the physiological perspective, there are three possible fates for solar radiation absorbed by plant chlorophyll molecules, including photosynthesis (PQ), heat (NPQ), and emitted fluorescence, and the relationship between them is generally competitive but actually neither unique nor simple [50]. For example, the relationships between steady state fluorescence and photochemistry under low and high light (or light-statured conditions with increasing irradiance and moisture) stress were a negative and positive correlation, respectively [51]. Under prolonged drought conditions, a decrease in SM (i.e., a decrease in soil water potential) would induce a decrease in plant water potential in order to absorb as much water as possible from the soil. Moreover, plants need to decrease transpiration to limit their water loss by regulating the opening and closing of the stomata on the leaf surface. However, the closing of the stomata prevents plants from taking up carbon dioxide (CO 2 ), thus limiting the dark reaction of photosynthesis and reducing the photosynthetic rate. It was expected that drought stress significantly decreased the leaf photosynthesis, and here, the SIF emissions also reduced along with the drought levels ( Figure 4). Therefore, the relationship between photosynthesis and fluorescence has a positive correlation. Plants can decrease the rate of photosynthesis rapidly in response to a drought, which similarly reduces SIF values quickly based on the positive and strong correlation between them. In contrast, the NDVI was mainly influenced by many factors, such as the chlorophyll content, LAI, and other indirect drivers such as the total plant cover, biomass, plant, and SM. These factors tend to decline gradually under prolonged drought conditions and lead to the late response of NDVI to SM changes reflecting the accumulative effect [52].
There is one important variable of SIF normalized by PAR (AF760) that provides a better understanding of their contributions to F760. It can be found that in extreme droughts (P4), the correlation between AF760 and SM was enhanced at different time scales and time windows compared with F760 (Figures 9-11). This finding proves the necessity of normalizing SIF by PAR in drought monitoring. Moreover, the correlation between AF760toc and the growth parameters (Chl and LAI) was enhanced (R 2 from 0.4 to 0.85 and R 2 from 0.45 to 0.65, respectively), as well as the degree of enhancement improving along with the severity of the drought, which may explain the better correlation between AF760toc and SM than that of F760. This coincides with the findings by Sun et al. [53] that SIF normalized by PAR can improve the net impacts of droughts on SIF.

Evaluation of NIR v for Drought Monitoring
NIR v , a newly proposed index which aims to address the mixed pixel problem and aspects of photosynthetic capacity [54], is highly correlated with SIF and can be used as the proxy of GPP (Wang et al., 2021) [55]. It has been demonstrated that NIR v and SIF exhibit similar performances for estimating GPP [56], while their differences have not been systematically explored under drought stress [57]. Generally, SIF is strongly driven by PAR or APAR (Figure 7), and NIR v is theoretically independent of PAR or APAR. SIF is less impacted by the soil background, and similarly, NIR v also has the potential for disentangling the vegetation signal from the soil's brightness and capturing the depth distributions of the canopy's photosynthetic capacity. Therefore, many studies have been conducted. For example, Wang et al. [58] demonstrated that NIR v and SIF were complementary for GPP estimation; that is, NIR v outperformed the other proxies in nonevergreen sparsely vegetated areas, whereas SIF performed best in high-productivity areas. The main reason for this lies in that the inherent noise associated with SIF retrieval [59] hindered the capturing of GPP dynamics. Aside from that, other NIR v -based variables such as NIR v,Rad , defined as the product of observed NIR radiance and the NDVI, has great potential for estimating the GPP from satellite remote sensing data [60].
As for drought monitoring, there are many indices for monitoring droughts in different scales, and each index has its advantages and limitations. For example, Wei et al. [61] provided a systematic study on evaluating the single drought indices, combined drought indices, and in situ drought indices for monitoring droughts across China. The results indicated that different drought indices have specific performance for drought monitoring in different land use types in China, and this research provided an effective way to explore the potential of different drought indices for monitoring droughts. Orimoloye et al. [62] utilized the Moderate Resolution Imaging Spectroradiometer (MODIS) EVI products to calculate the Vegetation Condition Monitoring Index (VCI), which was proven to be suited for monitoring drought disasters.
Meanwhile, for NIR v , for instance, during the 2012 North American drought, NIR v rapidly declined in parallel with the GPP, while the NDVI showed little deviation compared with non-drought years [63]. During the 2018 drought in northwestern Europe, strong reductions in SIF and NIR v were observed, corresponding to the occurrence of reduced precipitation [64]. Liu et al. [24] selected a severe drought event in 2016 over the Hulun Buir grasslands as an example and explored the potential of SIF and NIR v for drought detection and monitoring. Specifically, many parameters, including the VIs, self-calibrating PDSI, precipitation, SM, surface water storage, and land surface temperature were acquired to detect and characterize the drought event. Meanwhile, according to the composition of SIFtoc, the satellite SIF was used to calculate SIF yield (normalized by APAR), the physiological SIF emission yield, and SIFtot. The results indicated that NIR v outperformed other VIs in drought monitoring, with the largest decrease among the VIs, which was consistent with our study. For lag time analysis, similar results showing that SIF had a shorter lag time than the VIs in responding to meteorological droughts were found. Considering the relationship between NIR v and SM, Buitink et al. [65] used two satellite-derived vegetation indices, NIR v and Vegetation Optical Depth (VOD), and two SM networks for the 2018 agricultural drought in a lowland area. The results indicated that NIR v had the potential to infer the critical SM content. In this study, it was noted that in Figure 4, the NIR v values of four plots appeared to separate in the order of P2 > P1 > P3 > P4 in the first rectangle, thus indicating its greater potential for early drought monitoring than a VI. Furthermore, NIR v had a stronger correlation relationship with SM and the growth parameters in extreme droughts regardless of the different time windows or time lags (Figures 9-11). This also highlights its ability for long-term drought monitoring.

Conclusions
Here, the continuous ground measurements of SIF and the environmental and crop growth parameters of four plots under different levels of droughts for wheat crops were carried out, thus systematically exploring the potential of F760tot for drought monitoring. The results indicated that: (1) F760tot was capable of distinguishing the differences in different drought levels and responded quickly to the onset of moderate droughts compared with other variables, which appeared to have the greatest decrease; (2) compared with F760toc, F760tot appeared to be more related to the physiology and was subjected to the canopy structure less, but these relationships varied in extreme droughts; and (3) F760tot contained effective information on SM in terms of the correlations in short time lags, and the VIs were more strongly correlated with SM in the longer time lags. The results of this study demonstrate that F760tot was more sensitive to moderate droughts that usually appeared in the early stages of drought stress in plants, which may be attributed to the representation of the total emitted SIF and less influence from the canopy structure.
Conclusively, F760tot has great potential for early drought detection and monitoring. Our findings may contribute to mitigating the risk of agricultural droughts. In addition to these findings, because of the reabsorption effect, further works can be conducted to explore the performance of the total emitted SIF at the red band for drought monitoring. Aside from that, structural equation modeling (SEM) is a statistical technique that can be used to reveal the relationships of the observed and latent variables, and further studies are needed to use this method to illustrate the potential relationships between SIF and other variables, thus contributing to understand its driving factors.