Acceso abierto

Remote sensing techniques to assess chlorophyll fluorescence in support of crop monitoring in Poland

Miscellanea Geographica's Cover Image
Miscellanea Geographica
Thematic Issue: “Innovation in geospatial and 3D data” focuses on the newest achievements in the field of Geodata, which are used in Geosciences and for various applications such as urban planning, territorial management, damage assessment, environmental monitoring, 3D city modelling, renewable energy assessment, land registry, heritage documentation.

Cite

Introduction

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 a pigments, which are measured a few nanoseconds after light absorption (Damm et al. 2015). Accompanied by a photochemical reaction and heat dissipation, chlorophyll fluorescence emission is one of three pathways used for light consumption in the photosystem (Genty et al. 1989). Chlorophyll fluorescence is light that is re-emitted at a longer wavelength after being absorbed by chlorophyll molecules. Variable chlorophyll fluorescence is only observed in chlorophyll a, in photosystem II. Plant physiology can be investigated by measuring the intensity and nature of the variable chlorophyll fluorescence, and by using the protocols that have been developed (Pilar 2013). In many ChlF parameters, Fv/Fm is the maximum quantum yield (Photochemical efficiency of PSII) and this is used to characterize the conversion efficiency of the light energy of the PS II reaction center; its numerical changes are of special significance. However, conventional methods of assessing Fv/Fm from field observations that involve site-specific, complicated parameterizations and calculations, make it difficult to apply over large agricultural areas (Damm et al. 2015). These shortcomings can be overcome through the complementary use of satellite imagery.

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 (Zea mays L.) is one of the most important grain crop in the world and is produced throughout the world under diverse environments. In developed countries, maize is consumed mainly as second-cycle produce, in the form of meat, eggs, and dairy products (Plessis 2003). In developing countries, maize is consumed directly and serves as a staple diet for some 900 million poor people. By 2025, maize will be the developing world's largest crop and between now and 2050 the demand for maize in the developing world is expected to double. Sugar beet (Beta vulgaris L. subsp. Vulgaris), on the other hand, is a root crop that is the world's second sugar source (after sugar cane): 24% of world production. During the years 2010–2019, Poland was in the top ten countries of sugar beet producers.

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.

Study area

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).

Figure 1

JECAM cropland site 25 km × 25 km, highlighted in yellow, with field sites where chlorophyll fluorescence were measured

Source: own elaboration

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.

Materials and methods

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.

Figure 2

Scheme of work

Source: own elaboration

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.

Figure 3

Chlorophyll fluorescence measurements on sugar beets with OS5p+ Pulse Modulated Chlorophyll Fluorometer

Source: photo by Maciej Bartold

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

1 2 3 4 5
Application Index Description Equation Reference
Assessment of the general condition of vegetation CTVI Corrected Transformed Vegetation Index CTVI=(NDVI+0.5)|NDVI+0.5|*NDVI+0.5 {\rm{CTVI}} = {{\left( {{\rm{NDVI}} + 0.5} \right)} \over {\left| {{\rm{NDVI}} + 0.5} \right|}}*\sqrt {{\rm{NDVI}} + 0.5} Perry, 1984
DVI Difference Vegetation Index DVI= aRNIR-Rred Richardson, 1977
EVI Enhanced Vegetation Index EVI=RNIRRredRNIR+C1*RredC2*Rblue+L {\rm{EVI}} = {{{R_{NIR}} - {R_{red}}} \over {{R_{NIR}} + {C_1}*{R_{red}} - {C_2}*{R_{blue}} + L}} Huete, 1999
GEMI Global Environmental Monitoring Index GEMI=n(10.25n)Rred0.1251Rred {\rm{GEMI}} = {\rm{n}}\left( {1 - 0.25{\rm{n}}} \right) - {{{R_{red}} - 0.125} \over {1 - {R_{red}}}} Pinty, 1992
GNDVI Green Normalized Difference Vegetation Index GNDVI=RNIRRgreenRNIR+Rgreen {\rm{GNDVI}} = {{{R_{NIR}} - {R_{green}}} \over {{R_{NIR}} + {R_{green}}}} Gitelson, 1998
IRECI Inverted Red Edge Chlorophyll Index IRECI=(RNIRRred)Rrededge1/Rrededge2 {\rm{IRECI}} = {{\left( {{R_{NIR}} - {R_{red}}} \right)} \over {{R_{rededge1}}/{R_{rededge2}}}} Frampton et al., 2013
MSAVI Modified Soil Adjusted Vegetation Index MSAVI=RNIRRredRNIR+Rred+L*(1+L) {\rm{MSAVI}} = {{{R_{NIR}} - {R_{red}}} \over {{R_{NIR}} + {R_{red}} + L}}*\left( {1 + L} \right) Qi et al., 1994
MSAVI2 Modified Soil Adjusted Vegetation Index 2 MSAVI2=12[2*R800+1(2*R800+1)8*(R800R670)] {\rm{MSAVI}}2 = {1 \over 2}\left[ {2*{R_{800}} + 1 - \sqrt {\left( {2*{R_{800}} + 1} \right) - 8*\left( {{R_{800}} - {R_{670}}} \right)} } \right] Qi et al., 1994
NDREI1 Normalized Difference Red Edge Index 1 NDREI1=R790R720R790+R720 {\rm{NDREI}}1 = {{{R_{790}} - {R_{720}}} \over {{R_{790}} + {R_{720}}}} Gitelson And Merzlyak, 1994
NDREI2 Normalized Difference Red Edge Index 2 NDREI2=R750R705R750+R705*R445 {\rm{NDREI}}2 = {{{R_{750}} - {R_{705}}} \over {{R_{750}} + {R_{705}}*{R_{445}}}} Barnes, 2000
NDVI Normalized Difference Vegetation Index NDVI=RNIRRredRNIR+Rred {\rm{NDVI}} = {{{R_{NIR}} - {R_{red}}} \over {{R_{NIR}} + {R_{red}}}} Rouse, 1974
NRVI Normalized Ratio Vegetation Index NRVI=RredRNIR1RredRNIR+1 {\rm{NRVI}} = {{{{{R_{red}}} \over {{R_{NIR}}}} - 1} \over {{{{R_{red}}} \over {{R_{NIR}}}} + 1}} Baret, 1991
REIP Red Edge Inflection Point REIP=700+40((R670+R7802)R700R740R700) {\rm{REIP}} = 700 + 40\left( {{{\left( {{{{R_{670}} + {R_{780}}} \over 2}} \right) - {R_{700}}} \over {{R_{740}} - {R_{700}}}}} \right) Guyot And Barnet, 1988
RVI Ratio Vegetation Index RVI=RredRNIR {\rm{RVI}} = {{{R_{red}}} \over {{R_{NIR}}}} Bannari et al., 1995
SATVI Soil Adjusted Total Vegetation Index SATVI=RNIRRredRNIR+Rred+L*(1+L)RSWIR2 {\rm{SATVI}} = {{{R_{NIR}} - {R_{red}}} \over {{R_{NIR}} + {R_{red}} + L}}*\left( {1 + L} \right) - {{{R_{SWIR}}} \over 2} Marsett, 2006
SAVI Soil Adjusted Vegetation Index SAVI=(1+L)(RNIRRred)RNIR+Rred+L {\rm{SAVI}} = {{\left( {1 + L} \right)\left( {{R_{NIR}} - {R_{red}}} \right)} \over {{R_{NIR}} + {R_{red}} + L}} Huete, 1988
SLAVI Specific Leaf Area Vegetation Index SLAVI=RNIRRred+RSWIR {\rm{SLAVI}} = {{{R_{NIR}}} \over {{R_{red}} + {R_{SWIR}}}} Lymburger et al., 2000
SR Simple Ratio Vegetation Index SR=RNIRRred {\rm{SR}} = {{{R_{NIR}}} \over {{R_{red}}}} Birth, 1968
TTVI Thiam's Transformed Vegetation Index TTVI=|NDVI+0.5| {\rm{TTVI}} = \sqrt {\left| {{\rm{NDVI}} + 0.5} \right|} Thiam, 1997
TVI Transformed Vegetation Index MSAVI2=12[120*(R750R550)200*(R670R550)] {\rm{MSAVI}}2 = {1 \over 2}\left[ {120*\left( {{R_{750}} - {R_{550}}} \right) - 200*\left( {{R_{670}} - {R_{550}}} \right)} \right] Deering, 1975
WDVI Weighted Difference Vegetation Index WDVI=RNIR-a*Rred Richardson, 1977
Assessment of photosynthetically active pigment CLG Chlorophyll Index Green CLG=RNIRRgreen1 {\rm{CLG}} = {{{R_{NIR}}} \over {{R_{green}}}} - 1 Gitelson, 2003
CLRE Red-edge-band Chlorophyll Index CLRE=R750R7101 {\rm{CLRE}} = {{{R_{750}}} \over {{R_{710}}}} - 1 Gitelson, 2003
MCARI Modified Chlorophyll Absorption Ratio Index MCARI=[(R700R670)0.2*(R700R550)]*(R700/R670) {\rm{MCARI}} = \left[ {\left( {{R_{700}} - {R_{670}}} \right) - 0.2*\left( {{R_{700}} - {R_{550}}} \right)} \right]*\left( {{R_{700}}/{R_{670}}} \right) Daughtery, 2000
MTCI MERIS Terrestrial Chlorophyll Index MTCI=R754R709R709R681 {\rm{MTCI}} = {{{R_{754}} - {R_{709}}} \over {{R_{709}} - {R_{681}}}} Dash And Curran, 2004
S2REP Sentinel-2 Red-Edge Position Index S2REP=705+35*((RNIR+Rred)/2)R705(R740R705) {\rm{S}}2{\rm{REP}} = 705 + 35*{{\left( {\left( {{R_{NIR}} + {R_{red}}} \right)/2} \right) - {R_{705}}} \over {\left( {{R_{740}} - {R_{705}}} \right)}} Frampton et al., 2013
Assessment of the amount of light used in photosynthesis SIPI Structure Insensitive Pigment Index SIPI=R800R450R800+R650 {\rm{SIPI}} = {{{R_{800}} - {R_{450}}} \over {{R_{800}} + {R_{650}}}} Peñuelas et al., 1995
ZMI Zarco-Tejada & Miller Index ZMI=R750R710 {\rm{ZMI}} = {{{R_{750}}} \over {{R_{710}}}} Zarco-Tejada et al., 2001
Assessment of water content DSWI Disease Water Stress Index DSWI=R802+R547R1657+R682 {\rm{DSWI}} = {{{R_{802}} + {R_{547}}} \over {{R_{1657}} + {R_{682}}}} Galvão et al., 2005
MNDWI Modified Normalized Difference Water Index MNDWI=RgreenRMIRRgreen+RMIR {\rm{MNDWI}} = {{{R_{green}} - {R_{MIR}}} \over {{R_{green}} + {R_{MIR}}}} Xu, 2006
NDWI Normalized Difference Water Index NDWI=RgreenRNIRRgreen+RNIR {\rm{NDWI}} = {{{R_{green}} - {R_{NIR}}} \over {{R_{green}} + {R_{NIR}}}} McFeeters, 1996
NDWI2 Normalized Difference Water 2 Index NDWI2=R857R1241R587+R1241 {\rm{NDWI2}} = {{{R_{857}} - {R_{1241}}} \over {{R_{587}} + {R_{1241}}}} Gao, 1996
NDII Normalized Difference Infrared Index NDII=R850R1650R580+R1650 {\rm{NDII}} = {{{R_{850}} - {R_{1650}}} \over {{R_{580}} + {R_{1650}}}} Hardisky et al., 1993

Source: own elaboration

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.

Year 2018 2019
Month July August July August
Day 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

Source: own elaboration

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.

Results and discussion

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).

Figure 4

Time series of vegetation indices for maize and sugar beets during 2018–2019

Source: own elaboration

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.

Figure 5

Sentinel-3 based land surface temperatures of maize and sugar beets during 2018–2019. The red crosses indicate mean, thin black lines median, green boxes 25%–75% ranges, black dots are outliers

Source: own elaboration

Figure 6

Maximum and mean temperatures, as well as total precipitation, noted in 2018 at the meteorological station in the town of Kornik

Source: own elaboration

Figure 7

Maximum and mean temperature, as well as total precipitation, noted in 2019 at the meteorological station in the town of Kornik

Source: own elaboration

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 Vegetation Index Maize Sugar beets
R coefficient MAE RMSE p-value R coefficient MAE RMSE p-value
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

Source: own elaboration

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 Bistorta vivipara: SIPI (r=−0.46) and ZMI (r=0.61). Alternatively Zagajewski et al. (2017) used RapidEye satellite imagery to assess how well satellite data can detect vegetation stress on the ground. However, a simple ratio, such as NDVI, is characterized by low correlations and seem to be insufficient for detecting drought in plant condition (Zarco-Tejada et al. 2013; Liu et al. 2018). Moreover, analyses made on spectra collected in the SWIR region (1400–1550 nm and 1850–2000 nm) using airborne or satellite platforms are problematic due to the effects on measured signals of atmospheric water vapor (Zagajewski et al. 2017; Hillnhütter et al. 2011). The least significant relation between ChlF and VIs exist for sugar beets. Despite this, it should be noted that the highest correlation appears for EVI (r=0.45) and S2REP (r=0.43). The impact of single stresses, such as water deficit or drought, could be detected. However, a similar decrease in fluorescence intensity hindered the differentiation between both stresses. This fact was an overall obstacle in our study and was also observed for combined stresses. One possible explanation for these observations is that all chosen stresses reduce plant photosynthetic efficiency and the differentiation becomes more complicated due to the interaction between the specific processes of the individual stresses (Leufen et al. 2014). The results of the 2013 study by Li et. al. show that sugar beets tolerate tissue dehydration better, and can be distinguished from sensitive sugar beets that rely on the efficiency of electron transport in PSII.

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)

CROP TYPE INDEX FORMULA
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

Source: own elaboration

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.

Conclusions

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.

eISSN:
2084-6118
Idioma:
Inglés
Calendario de la edición:
4 veces al año
Temas de la revista:
Geosciences, Geography, other