The remote sensing of chlorophyll fluorescence is a rapidly advancing front in terrestrial vegetation science, with emerging capabilities in space-based methodologies and prospects for diverse applications (Mohammed et al. 2019). Chlorophyll fluorescence (ChlF) differs from the traditional reflectance-based vegetation indices as it refers to the red and far-red wavelengths (650–800 nm) emitted by chlorophyll
Chlorophyll fluorescence is widely used to monitor crop stress, while multi-frequent satellite observations allow for the retrieval of parameters such as pigment concentration and photosynthetic active radiation (PAR), which can be used for ChlF estimations. It is known that crop yield can be reduced by stress factors such as extreme temperatures, direct sunlight, and a shortage of water and nutrients (Baligar et al. 2001; Dąbrowska-Zielińska et al. 2011), which affect or block plants’ metabolism, growth, and development (Kranner et al. 2010). Compared to reflectance, induced fluorescence is considered to be a more accurate indicator of plant state and to be able to detect stress impacts at earlier growth stages (McFarlane et al. 1980; Kancheva 2008).
Vegetation indices derived from satellite images have, alternatively, been used to estimate plant stress over larger areas based on spectral reflectance (Bochenek et al. 2017; Gurdak & Grzybowski 2018). Unfortunately, most of the common vegetation indices measure only the “greenness” of the land cover and not photosynthesis itself. Consequently they have little or no sensitivity to short-term physiological changes in leaves (Zarco-Tejada et al. 2013). In the short term (hours or days), green leaves remain green but reduce photosynthesis when they are stressed. At longer time scales, when prolonged stress causes premature senescence, the stress can be monitored as a change in, for example, NDVI (Mathobo et al. 2017). ChlF is directly related to crop photosynthesis and is, therefore, recognized as a more reliable indicator of stress than VIs (Meroni et al. 2009; Li et al. 2018). Therefore, some efforts have been undertaken in research centers in the past years to prepare effective methods for deriving this variable from satellite data and applying it to crop estimations for wheat, barley, and maize in Europe and Asia (Lopez-Lozano et al. 2015; Wei et al. 2019).
On the other hand, among Essential Climatic Variables (ECV), remote-sensing based on the land surface temperature (LST) can be used to estimate plant water-vapor loss (gl, water vapor) (Miguel Costa et al. 2013). This is potentially an indicator of stomatal opening, which in turn is strongly correlated with photosynthesis (Chaves et al. 2003; Jones 1992). In monitoring individual Essential Climatic Variables, such as air temperature, precipitation derived from meteorological stations is also required for studying photosynthesis dynamics (Gobron et al. 2006).
Photosynthesis is one of the most important photochemical processes, however, it is difficult to measure at field scale. One promising approach is the use of ChlF, which is directly related to photosynthesis and can be measured at field scale with remote sensing techniques (Drusch et al. 2017). In 2019, the European Space Agency (ESA) announced that it has signed a contract with Thales Alenia Space to lead the Fluorescence Explorer (FLEX) satellite mission, which is scheduled for launch in 2023. Also, taking advantage of Sentinel-3's (launched in 2018) optical and thermal sensors will lead to an integrated package of measurements that can assess plant health. With the Sentinel-2 satellites, also in orbit, there is the unique opportunity for using the data of all three missions synergistically for vegetation studies. The FLEX mission, in tandem with Sentinel 3, aims to provide global maps of vegetation fluorescence that can reflect photosynthetic activity, and plant health and stress. In turn, this is important for better agricultural management and food security.
Maize (
The most important requirement for growing maize and sugar beets is the soil, which must contain a large supply of plant food, be rich in humus, and have the property of retaining a great deal of moisture. Climatic conditions, both air and surface temperature, sunshine, rainfall, and wind, all have an important bearing upon high yields of maize and sugar beets (Rolph 1873; Fernandez-Armesto 2011). Drought and moisture stress that occurs during the different development stages of maize and sugar beets may reduce the final grain yield to different degrees, with the extent of yield reduction depending not only on the severity of the stress, but also on the stage of the plants’ development (Wilson 1968; Claasen & Shaw 1970). Due to the plants’ water requirements, the most crucial phenological stages are at the stem elongation stage (maize) and at the rosette growth stage (sugar beets) (Doorenbos & Kassam 1979; Garrod 1974).
The objective of this paper is to examine the feasibility of estimating ChlF from VIs derived from the Sentinel-2 satellite. Monitoring crop stress and investigating the climatic variables that affect the phenological dynamics associated with temporal changes in VIs are very important. We evaluated of the usefulness of thirty three VIs that monitored the impact of various environmental conditions by finding how they related to the ground measurements of ChlF. To monitor the stress conditions of plants, a detailed temporal analysis of meteorological parameters was carried out in combination with the temporal variability of VIs derived from Sentinel-2 images. In order to investigate the thermal stress we also observed the land surface temperature derived from Sentinel-3 at 1km resolution, in conjunction with a meteorological database containing air temperature and precipitation. The outcome of the statistical analysis allowed the authors to find the relations between the chosen remote sensing-based indices and the ground measured ChlF. The authors were able to draw conclusions about the stress conditions that could impact the growth of maize and sugar beets, reflecting the changes in VI values, using temporal profiles for climate variables and land surface temperature.
The cropland test site is located in western Poland (Fig. 1), within the administrative area of Wielkopolska voivodeship (NUTS-2 PL41). Wielkopolska lies within the Odra river basin, however, 88% of the province's surface water drains into the Warta river basin. The flat terrain reaches elevation ranges of 60 to 90 meters. The western part of Wielkopolska is influenced by oceanic air masses that affect the mildness of the climate. On the other hand the eastern region of Wielkopolska is under the influence of a more distinctly continental climate. The average temperature for the year is approximately 8°C. The growing season is one of the longest in Poland: the southern plains are characterized by a growing season of around 228 days, while the northern region gradually declines to 216 days. Precipitation ranges from 500 to 550 mm per year (Gurdak & Grzybowski 2018).
The field site covered an area of 25 km × 25 km and was located in Wielkopolska. The coordinates are as follows: upper left 16°38′21″E, 52°10′31,500″ N; lower right 17°0′53: E, 51°57′27″ N. The site has been incorporated into JECAM (Joint Experiment of Crop Assessment and Monitoring), which is an international program for agricultural monitoring. JECAM was formed as part of the GEOGLAM (Group on Earth Observations Global Agricultural Monitoring) initiative and aims to carry out experiments to facilitate the fusion of satellite and field data for monitoring crop growth conditions and developing yield forecasting models. Its field experiments are carried out at a series of sites that represent the world's main cropping systems and rely on collecting in-situ vegetation biophysical parameters and meteorological measurements. In Poland, the JECAM site is characterized by a mixture of agricultural crops, with winter wheat, rape, maize, and sugar beets as the dominant species. The field pattern is composed of large fields that dominate the test site, intermixed with small fields. The crop calendar for winter cereals begins in September and lasts till July, while the spring/summer crop calendar starts in April and ends in October. The main crop types are winter/spring wheat, triticale, and barley; and winter rape and rye, maize, sugar beets, alfalfa, and potatoes. The analyzes were performed for maize and sugar beets.
In 2018, seven fields of maize and five fields of sugar beets were selected. In 2019, the number of plots increased to twelve for maize and nine for sugar beets. The area of the fields was between 5 ha and 50 ha, and formed a rough square shape. The maize was seeded at the beginning of May and harvested at the turn of September and October. The date for seeding sugar beets was the beginning of April, the harvest was in mid-October. Figure 1 presents the distribution of corn and sugar beet fields on which the ground measurements were conducted.
The ChlF ground measurements were conducted on the selected JECAM plots during the growing season. The actual locations of the plots of maize and sugar beet were investigated during field campaigns during July–August (2018–2019), taking into account the crop rotation system used in Poland. In-situ data were collected during two field campaigns: 7–9 August 2018 and 14–16 July 2019, at the stem elongation stage (maize) and at the rosette growth stage (sugar beets). These are the phases most susceptible to damage caused by drought and water stress (Doorenbos & Kassam 1979; Garrod 1974).
During 2018–2019, nineteen and fourteen crop fields were registered as maize and sugar beets, respectively, in the study area (Fig. 2). Ground measurement data and the related satellite-based indices were grouped into two sets: a training set for determining relations (23 points), comprising 70% of all collected data, and a test set for verifying the established relations (10 points), comprising 30% of all data. The regression equations obtained as a result of the correlation analysis were next applied to generate ChlF values at the test points. Finally, these values were compared with the ground ChlF measurements and the differences between these two datasets were computed and statistically assessed.
We used an OS5p+ Pulse Modulated Chlorophyll Fluorometer (OPTI-SCIENCES, USA), which registers maximum quantum yield (FV/FM) using the dark adapted test (Fig. 3). The dark adapted test is a measurement ratio that represents the maximum potential quantum efficiency of the Photosystem II if all capable reaction centers are open. This ratio is an estimate of the maximum portion of absorbed quanta used in the PSII reaction centers (Kitajima & Butler 1975). By dark adapting, one is allowing the re-oxidation of the PSII and the relaxation of the NPQ. The appropriate dark adaption time was between 20 and 35 minutes for best results. A series of ten dark adapted white clips were provided with the system to be used for dark adaption measurement. The Elementary Sampling Unit (ESU) for ten clips was a 30 × 30 m square, for the correct characterization of a 10 m Sentinel-2 pixel. The clips were placed on the leaves, with the black slider covering the cylindrical opening. After dark adaption, the end of the fiber optic bundle was placed in the cylindrical opening and the dark slide of the clip was opened allowing the sample to be exposed to the fiber optic bundle. FV/FM value ranges of 0.66 to 0.83 were obtained from ground measurement. The approximate optimal value range is 0.79 to 0.83 for most land plant species, with lowered values indicating plant stress (Maxwell & Johnson 2000). From the ten measurements carried out for each field, the average value was calculated as FV/FM.
Individual ranges for the electromagnetic spectrum are used for a detailed analysis of the state of plants. This is done by calculating the remote sensing vegetation indicators (VI); these are VIs that use various types of mathematical combinations of the respective spectral reflection indicators. VIs based on narrow ranges of the electromagnetic spectrum allow the content of individual substances in plants to be analyzed (Kycko 2017; Zagajewski et al. 2017).
The vegetation indices were selected and calculated from the Sentinel-2 Multispectral instrument. Cloudless satellite images were acquired on the nearest date to the ground measurement campaigns in order to ensure the comparability of the results (Sentinel-2B 2018/08/09 and Sentinel-2A 2019/07/25). The product type was Level-2A (after atmospheric correction). At the preliminary stage of the work thirty three vegetation indices, which characterized different aspects of crop condition and development, were derived from the Sentinel-2 data. These can be arranged into four groups (Table 1).
Vegetation indices calculated from Sentinel-2 satellite imagery
Assessment of the general condition of vegetation | CTVI | Corrected Transformed Vegetation Index |
|
Perry, 1984 |
DVI | Difference Vegetation Index | DVI= |
Richardson, 1977 | |
EVI | Enhanced Vegetation Index |
|
Huete, 1999 | |
GEMI | Global Environmental Monitoring Index |
|
Pinty, 1992 | |
GNDVI | Green Normalized Difference Vegetation Index |
|
Gitelson, 1998 | |
IRECI | Inverted Red Edge Chlorophyll Index |
|
Frampton et al., 2013 | |
MSAVI | Modified Soil Adjusted Vegetation Index |
|
Qi et al., 1994 | |
MSAVI2 | Modified Soil Adjusted Vegetation Index 2 |
|
Qi et al., 1994 | |
NDREI1 | Normalized Difference Red Edge Index 1 |
|
Gitelson And Merzlyak, 1994 | |
NDREI2 | Normalized Difference Red Edge Index 2 |
|
Barnes, 2000 | |
NDVI | Normalized Difference Vegetation Index |
|
Rouse, 1974 | |
NRVI | Normalized Ratio Vegetation Index |
|
Baret, 1991 | |
REIP | Red Edge Inflection Point |
|
Guyot And Barnet, 1988 | |
RVI | Ratio Vegetation Index |
|
Bannari et al., 1995 | |
SATVI | Soil Adjusted Total Vegetation Index |
|
Marsett, 2006 | |
SAVI | Soil Adjusted Vegetation Index |
|
Huete, 1988 | |
SLAVI | Specific Leaf Area Vegetation Index |
|
Lymburger et al., 2000 | |
SR | Simple Ratio Vegetation Index |
|
Birth, 1968 | |
TTVI | Thiam's Transformed Vegetation Index |
|
Thiam, 1997 | |
TVI | Transformed Vegetation Index |
|
Deering, 1975 | |
WDVI | Weighted Difference Vegetation Index | WDVI=RNIR-a*Rred | Richardson, 1977 | |
Assessment of photosynthetically active pigment | CLG | Chlorophyll Index Green |
|
Gitelson, 2003 |
CLRE | Red-edge-band Chlorophyll Index |
|
Gitelson, 2003 | |
MCARI | Modified Chlorophyll Absorption Ratio Index |
|
Daughtery, 2000 | |
MTCI | MERIS Terrestrial Chlorophyll Index |
|
Dash And Curran, 2004 | |
S2REP | Sentinel-2 Red-Edge Position Index |
|
Frampton et al., 2013 | |
Assessment of the amount of light used in photosynthesis | SIPI | Structure Insensitive Pigment Index |
|
Peñuelas et al., 1995 |
ZMI | Zarco-Tejada & Miller Index |
|
Zarco-Tejada et al., 2001 | |
Assessment of water content | DSWI | Disease Water Stress Index |
|
Galvão et al., 2005 |
MNDWI | Modified Normalized Difference Water Index |
|
Xu, 2006 | |
NDWI | Normalized Difference Water Index |
|
McFeeters, 1996 | |
NDWI2 | Normalized Difference Water 2 Index |
|
Gao, 1996 | |
NDII | Normalized Difference Infrared Index |
|
Hardisky et al., 1993 |
In order to analyze the thermal stress on plant growth we used the satellite temperature provided by the Copernicus Sentinel-3 mission. The temperature observations registered by a Sea and Land Surface Temperature Radiometer (SLSTR) instrument, mounted on-board both the Sentinel-3A and Sentinel-3B satellites, was downloaded. The Level-2 LST product from the satellite-based temperature data that provides land surface parameters, generated a 1 km wide measurement grid, which was used in the analyzes.
The temperature readings from the measurement points were made on the basis of the Sentinel-3 LST data. Next we carried out a series of descriptive statistical calculations, which were presented on box plots (mean (red cross), median (black line in green box), maximum and minimum (upper and lower whisker), outliers (black dots), 1st quartile and 3rd quartile (sides of green box)) on surface temperature data (Fig. 5) on fields with sugar beets and maize within the JECAM area. The 1km gridded Sentinel-3 images showing temperature were primarily resampled as ten meter resolution images in order to compare the results with the Sentinel-2 based vegetation indices and find the relation between temperature conditions and plant growth.
At the beginning of the Sentinel-3 mission's operational term, (i.e. July 2017), maize and sugar beet croplands were taken into account for validating surface temperatures registered by satellites. We analyzed twenty-five cloud-free Sentinel-3 satellite images of the study area during 2018–2019 (Table 2).
Sentinel-3 cloud-free satellite images of croplands in Wielkopolska. Actual state on 10. September 2019.
July | August | July | August | ||||||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
5 | 9 | 21 | 28 | 1 | 5 | 13 | 17 | 4 | 15 | 25 | 26 | 1 | 8 | 12 | 13 | 15 | 22 | 23 | 24 | 26 | 28 | 29 | 30 | 31 |
In the next stage, meteorological data were compiled from weather stations adjoining the study area. Meteorological data were collected from the closest station, which belonged to the national meteorological data network (the town of Kornik, 52°14′N, 17°6″E). Three meteorological parameters were analyzed in order to find anomalous weather periods that could affect the condition of the crops within study area: daily maximum and mean temperatures (measured at a height of 2 m from the ground) as well as the sum of the daily precipitation collected at the station from 2018 to 2019.
The statistical analyzes were conducted using Statistica 13 software. The acquired data had a normal distribution. The r-Pearson correlation analysis was used to calculate the correlation between the VIs and ChlF. Thereafter, linear polynomial regression was used to compute the ChlF from the VIs. In order to assess the accuracy, the following statistical analyses were used: r coefficient, mean absolute error (MAE), root mean squared error (RMSE), and p-value. Verification was carried out on the test set (10 points), comprising 30% of all data.
First, an analysis of the temporal variation of VIs was performed. Figure 4 illustrates the growth of maize and sugar beets during 2018–2019, and the differences in the plants’ conditions observed within the two-year period. A time series of vegetation indices began in mid-May with the germination phase and ended with the heading of the maize and the development of the beets’ root phases in July or August. Profiles were prepared for the commonly used indices that represented each of three chosen groups of VIs presented in Table 1. The average index values from the fields were calculated from the Sentinel-2 images and divided into the individual stages of crop development (Fig. 4).
The values of the NDVI and NDII indices increased during the germination and leaf development phases. During stem elongation and rosette growth the given indice values reached their peaks, after which the values decreased, achieving a level close to, or the same as, that at the beginning of growing season. The observed NDVI and NDII temporal curves are comparable to the runs noted by Masialeti et al. (2010) and Pan et al. (2015), thus, they appear to be characteristic for maize and sugar beet crops. For maize, the NDII value was 0.19 in 2018, and 0.29 in 2019. While the NDVI reached 0.66 and 0.73 in 2018 and 2019, respectively, at the stem elongation stage. For sugar beets at the rosette growth stage, we observed an NDII of 0.45 in 2018 and 0.38 in 2019, while the NDVI values reach 0.83 in 2018 and 0.75 in 2019. Estimates of the ratio of carotenoids to chlorophyll can be followed using the SIPI (structural independent pigment index) curves. The distinctly increased values observed for maize at heading and sugar beets at the beets’ root development stage in 2019 confirm that the vegetation was stressed (Kycko et al. 2019a). Assessing the land surface temperatures derived from the low-resolution satellite data, we noticed a consistency in the Sentinel-3 temperature registrations from all nineteen maize fields as well as from all ten sugar beet fields. The Sentinel-3 LST mean and median values extracted from included crop fields were highly consistent; for example, the average difference between the mean and median (Δmean, median) across all the Sentinel-3 satellite registrations of maize for 2018 was 0.42°C (Fig. 5), whereas the averaged Δmean, median from the satellite observations in 2019 was 0.17°C (Fig. 5). Next, the average difference between the mean and median Sentinel-3 LST for sugar beet in 2018 was 0.36°C (Fig. 5). In contrast, the land surface temperatures averaged Δmean, median in 2019 for sugar beet fields was 0.24°C (Fig. 5). The greater consistency of satellite temperature registrations during 2019 compared to 2018 could be the effect of long-term and stable drought conditions over the whole study area (Fig. 6, Fig. 7). Having determined the change in conditions over two years, we were able to draw conclusions about the influence of weather conditions on vegetation's state.
At the beginning of July 2018 we noticed higher maximum air temperatures, reaching 30°C, in contrast to the first days of July 2019, which had maximum air temperatures of around 20°C. We also observed corresponding changes in surface temperatures for the same period (Fig. 5). Spring and summer 2018 were characterized by short and intense periods of rain, with a total precipitation of 40–50 mm during one week; as well as periods of heat (Fig. 6) that may have significantly worsened drought conditions by hastening water evaporation and thermal stress on crops. During summer 2019 we noticed that precipitation was evenly distributed, with a total of 15–20 mm of rain every week (Fig. 7). On the other hand, we noted periodical cycles in surface temperatures. Thus the favorable meteorological conditions before the growing season in 2019 implied the stable growth of maize and sugar beets until August, when we noticed stressed vegetation, confirmed by clearly increased SIPI values (Fig. 4).
The low NDII and high SIPI values noted during the last phenological phases were associated with climatic factors such as an increased surface temperature and response to drought conditions.
At the next stage of the work, correlation analysis was performed comparing particular VIs derived from satellite images and those ChlF values measured on the ground in order to find if a significant relationship existed between these two parameters. Assessments of Fv/Fm from VIs derived from remote sensing data have been reported by several studies (Tan et al. 2012; Peng et al. 2017; Kycko et al. 2019a; Kycko et al. 2019b; Wei et al. 2019; Zagajewski et al, 2017; Zagajewski et al. 2018). For instance, some researchers compared the results of VIs to assess the Fv/Fm of maize, and reached the conclusion that out of the various kinds of VI that had a close relationship with Fv/Fm, SIPI performed best (Tan et al, 2012). If ground cover was significant, the impact of the background was significantly reduced, and Fv/Fm could be better estimated using NDVI (Wei et al. 2019). Kycko et al. (2019a) and Zagajewski et al. (2018) used in-situ hyperspectral remote sensing as an important technique to fulfill real-time monitoring of plant condition based on its superior performance in acquiring vegetation canopy information rapidly and non-destructively. In studies conducted by Peng et. al. (2017) an approach was demonstrated to estimate ChlF using spectral indices calculated from leaf reflectance measured using a hyperspectral radiometer, however, the reflectance spectra were resampled to simulate the spectral bands of the MSI aboard Sentinel-2.
Correlation analysis was performed between ground measured ChlF and VIs derived from Sentinel-2 images (Table 3).
Results of correlation analysis of ground measured ChlF and Vis (red color – results with the highest correlation coefficient)
Assessment of the general condition of vegetation | CTVI | 0.19 | 0.08 | 0.07 | 0.445 | 0.36 | 0.08 | 0.07 | 0.204 |
DVI | 0.24 | 0.07 | 0.07 | 0.315 | 0.29 | 0.08 | 0.08 | 0.320 | |
EVI | 0.61 | 0.03 | 0.03 | 0.005 | 0.45 | 0.05 | 0.04 | 0.106 | |
GEMI | −0.26 | 0.07 | 0.07 | 0.275 | −0.25 | 0.08 | 0.08 | 0.394 | |
GNDVI | 0.07 | 0.10 | 0.09 | 0.767 | 0.35 | 0.09 | 0.09 | 0.226 | |
IRECI | 0.64 | 0.02 | 0.02 | 0.003 | 0.22 | 0.08 | 0.08 | 0.447 | |
MSAVI | 0.17 | 0.08 | 0.07 | 0.478 | 0.37 | 0.07 | 0.06 | 0.188 | |
MSAVI2 | 0.17 | 0.08 | 0.07 | 0.478 | 0.37 | 0.07 | 0.06 | 0.188 | |
NDREI1 | 0.14 | 0.08 | 0.07 | 0.579 | 0.35 | 0.07 | 0.07 | 0.213 | |
NDREI2 | 0.06 | 0.10 | 0.09 | 0.815 | 0.35 | 0.07 | 0.07 | 0.223 | |
NDVI | 0.19 | 0.08 | 0.08 | 0.429 | 0.35 | 0.07 | 0.07 | 0.214 | |
NRVI | −0.19 | 0.08 | 0.08 | 0.429 | −0.35 | 0.07 | 0.07 | 0.214 | |
REIP | 0.15 | 0.08 | 0.07 | 0.538 | 0.31 | 0.08 | 0.07 | 0.284 | |
RVI | −0.17 | 0.08 | 0.07 | 0.478 | −0.37 | 0.07 | 0.07 | 0.188 | |
SATVI | 0.22 | 0.07 | 0.07 | 0.362 | 0.40 | 0.05 | 0.05 | 0.152 | |
SAVI | 0.56 | 0.04 | 0.03 | 0.012 | 0.32 | 0.08 | 0.07 | 0.269 | |
SLAVI | 0.37 | 0.06 | 0.06 | 0.124 | 0.26 | 0.08 | 0.07 | 0.374 | |
SR | 0.29 | 0.07 | 0.06 | 0.233 | 0.21 | 0.08 | 0.08 | 0.479 | |
TTVI | 0.19 | 0.08 | 0.07 | 0.445 | 0.36 | 0.07 | 0.07 | 0.478 | |
TVI | 0.55 | 0.04 | 0.04 | 0.015 | 0.21 | 0.08 | 0.08 | 0.478 | |
WDVI | 0.24 | 0.07 | 0.06 | 0.315 | 0.29 | 0.08 | 0.07 | 0.320 | |
Assessment of photosynthetically active pigment | CLG | 0.06 | 0.10 | 0.10 | 0.805 | 0.24 | 0.07 | 0.07 | 0.401 |
CLRE | 0.09 | 0.10 | 0.09 | 0.708 | 0.28 | 0.08 | 0.07 | 0.335 | |
MCARI | 0.35 | 0.06 | 0.06 | 0.141 | 0.26 | 0.07 | 0.06 | 0.377 | |
MTCI | −0.01 | 0.11 | 0.10 | 0.980 | 0.38 | 0.06 | 0.06 | 0.186 | |
S2REP | 0.46 | 0.06 | 0.05 | 0.046 | 0.43 | 0.05 | 0.04 | 0.125 | |
Assessment of the amount of light used in photosynthesis | SIPI | −0.68 | 0.02 | 0.02 | 0.001 | 0.24 | 0.07 | 0.07 | 0.401 |
ZMI | 0.54 | 0.04 | 0.04 | 0.017 | −0.13 | 0.09 | 0.08 | 0.665 | |
Assessment of water content | DSWI | 0.64 | 0.03 | 0.03 | 0.003 | 0.22 | 0.07 | 0.06 | 0.460 |
MNDWI | 0.54 | 0.04 | 0.04 | 0.017 | 0.09 | 0.09 | 0.09 | 0.751 | |
NDWI | −0.07 | 0.10 | 0.09 | 0.767 | −0.35 | 0.06 | 0.06 | 0.226 | |
NDWI2 | 0.37 | 0.06 | 0.06 | 0.120 | 0.33 | 0.07 | 0.06 | 0.251 | |
NDII | 0.65 | 0.03 | 0.03 | 0.002 | 0.31 | 0.07 | 0.07 | 0.279 |
The results obtained are supported by the phenological behavior of the maize and sugar beets, which is related to water content and the efficiency of photosynthesis (Rolph 1873; Plessis 2003). For maize, the highest correlation between ground measured ChlF and VIs appears for NDII (r=0.65), and at a negative correlation for SIPI (r=−0.68). Similar ranges of correctness were achieved by Kycko et. al. (2019b) when analyzing the influence of lead ions on the growth of pea plants. Tan et. al. (2012) reported that SIPI was the most sensitive remote sensing variable for monitoring Fv/Fm, with a correlation coefficient of r=−0.88. These reflectance indices should be effective for estimating ChlF in respective crops by using simple reflectance techniques (Peñuelas & Filella 1998; Kycko et al. 2014). High correlations were also obtained for VIs that contained various wavelengths. Some of the most significant ranges overlapped with the research by Zagajewski et al., (2018) for
The regression equations derived from the correlation analysis served to determine ChlF values on the basis of particular vegetation indices. The formulas for the linear regression equations used to estimate ChlF (FV/FM) are given in Table 4. Next, ChlF values, calculated based on the satellite VIs for control points, were compared with those measured on the ground in order to estimate the accuracy of the ChlF determination.
Linear regression equations for estimating ChlF (FV/FM)
MAIZE | NDII | FV/FM = 0.68053 + 0.25102 * NDII |
SIPI | FV/FM = 1.1527 − 0.3640 * SIPI | |
SUGAR BEETS | EVI | FV/FM = 0.60230 + 0.07402 * EVI |
S2REP | FV/FM = −8.753 + 0.01317 * S2REP |
The root mean square error (RMSE) for the predicted ChlF (FV/FM) from NDII was 0.033, and from SIPI was 0.032 for maize. While the RMSE for sugar beets was the same for both indices (EVI and S2REP) and equal to 0.038. Tan et al. (2012) achieved similar results measuring maize using a model established on the basis of SIPI, in which the RMSE was 0.082. However, EVI and NDVI noted by Wei et. al (2019), showed an RMSE below 0.148.
Reflectance indices such as structural independent pigment index (SIPI) and normalized difference infrared index (NDII), which are related to pigment and water content, may be useful in assessing ChlF. This adds new possibilities for indirectly assessing progressive leaf water stress. Average values of ChlF showed significant reductions under relatively serious drought, while simple VI means did not show any remarkable change until there was extreme drought. Our results confirmed the relations between ChlF measured at field sites and ChlF retrieved from satellite data. Therefore ChlF high-temporal satellite retrievals provide an opportunity for following changes in ChlF and finding the effects of drought stress on ChlF changes. The study of the relationships between reference chlorophyll fluorescence and vegetation indices derived from Sentinel-2 images led to the conclusion that there is quite a significant relation between these two types of data. Both indices – SIPI and NDII – derived from the Sentinel-2 during the maize vegetation season correlated well with ground measured fluorescence (r = −0.68 and r = 0.65 respectively). Moreover the temporal dynamics of SIPI and NDII reflect the surface temperature changes and meteorological conditions existing during the crucial phase of the growing season, especially at the end of July and August. This means that Sentinel-2 satellite data can be applied to monitoring stress conditions during the growing season. Nevertheless, it should be mentioned that the significant impact of drought conditions in 2018 and 2019 was observed for sugar beets, while analysis of the Sentinel-2 based vegetation indices and in-situ measured ChlF revealed lower correlation, at r=0.43 for S2REP and r=0.46 for EVI. The results confirmed that both VIs – EVI and S2REP – are reliable indicators for monitoring maize and sugar beet conditions related to chlorophyll content.
EVI and S2REP indices are desirable indicators for monitoring photosynthesis (or chlorophyll content), while the SIPI index (related to the amount of active light) and NDII (related to water content) better reflect stress conditions and ChlF. High SIPI values (increased carotenoids and decreased chlorophyll) are often an indicator of plant disease, which is associated with loss of chlorophyll in plants. Therefore, there is a negative correlation with the SIPI index, which provides information about stress conditions. The worst results in sugar beets may be caused by the fact that the signal acquired by the satellites was not homogeneous. Sugar beet leaves do not cover soil background signals, while the top canopy of maize cover the bare soil under them, especially in areas of dense canopy.
The results of this study indicate further research is needed, which should be based on a synergy of low and high resolution satellite data; these data will enable a more detailed analysis for estimating fluorescence and its relation to climatic conditions, environmental aspects, and indices derived from satellite images, taking into account the variability of other crop fields within the study area. In the near future, we expect that retrieval algorithms will become operational that allow the fluorescence for crop stress condition monitoring and crop estimation to be quantified.