Skip to main content

Comparison of visible imaging, thermography and spectrometry methods to evaluate the effect of Heterodera schachtii inoculation on sugar beets



Phenotyping technologies are expected to provide predictive power for a range of applications in plant and crop sciences. Here, we use the disease pressure of Beet Cyst Nematodes (BCN) on sugar beet as an illustrative example to test the specific capabilities of different methods. Strong links between the above and belowground parts of sugar beet plants have made BCN suitable targets for use of non-destructive phenotyping methods. We compared the ability of visible light imaging, thermography and spectrometry to evaluate the effect of BCN on the growth of sugar beet plants.


Two microplot experiments were sown with the nematode susceptible cultivar Aimanta and the nematode tolerant cultivar BlueFox under semi-field conditions. Visible imaging, thermal imaging and spectrometry were carried out on BCN infested and non-infested plants at different times during the plant development. Effects of a chemical nematicide were also evaluated using the three phenotyping methods. Leaf and beet biomass were measured at harvest. For both susceptible and tolerant cultivar, canopy area extracted from visible images was the earliest nematode stress indicator. Using such canopy area parameter, delay in leaf growth as well as benefit from a chemical nematicide could be detected already 15 days after sowing. Spectrometry was suitable to identify the stress even when the canopy reached full coverage. Thermography could only detect stress on the susceptible cultivar. Spectral Vegetation Indices related to canopy cover (NDVI and MCARI2) and chlorophyll content (CHLG) were correlated with the final yield (R = 0.69 on average for the susceptible cultivar) and the final nematode population in the soil (R = 0.78 on average for the susceptible cultivar).


In this paper we compare the use of visible imaging, thermography and spectrometry over two cultivars and 2 years under outdoor conditions. The three different techniques have their specific strengths in identifying BCN symptoms according to the type of cultivars and the growth stages of the sugar beet plants. Early detection of nematicide benefit and high yield predictability using visible imaging and spectrometry suggests promising applications for agricultural research and precision agriculture.


The rapid development of sensitive tools for plant phenotyping allows the assessment of very complex traits such as root morphology, biomass, leaf characteristic, yield related traits, biotic and abiotic response [1,2,3]. In most cases, phenotyping approaches are tested independently under a given scenario which does not facilitate the objective comparison of the methods tested. Often, the different methods are investigated at various scales (field or greenhouse) by following diverse protocols (cultivar, type and level of infestation, growth duration). Sugar beet is an interesting crop since the harvested organ develops vegetatively, thereby integrating environmental effects over time. It has recently been shown that beet development is reflected by aboveground development facilitating the use of shoot phenotyping procedures for yield estimation and disease effects [4]. On sugar beet, limited studies have been published which makes it difficult to evaluate the advantages and disadvantages of the different phenotyping approaches to characterize nematode symptoms.

Nematodes are soil borne parasites that occur naturally in soil. They cause annually up to 20% of yield losses in crops such as soybean, cotton, cereals, tuber crops, legumes, fruit and vegetables [5]. Sugar beet is a root crop which is widely cultivated in Europe and North America for sugar production. The sugar beet cyst nematode Heterodera schachtii (Schmidt) is a major threat and can cause severe beet damage and compromise the final yield. It has been demonstrated that there is a strong link between number of nematodes and crop performance such as shoot development and root biomass accumulation [6, 7].

In order to manage the damage caused by nematodes, dedicated strategies have been developed. A first approach consists of evaluating whether the level of infestation in the field is above a given economical threshold thereby justifying specific nematode control methods. However, soil sample analyses are expensive and technically difficult because of the cluster distribution of BCN in the field [8, 9]. Thus, many samples per hectare are required to achieve a reasonable estimation of the potential crop damage.

To reduce costs and increase the spatial resolution of BCN soil pressure evaluation, non-destructive methods have been developed [10]. It is worth noting that BCN occurs in patches in the field, has a low mobility, and causes diverse and rather generic visible aboveground symptoms, for example stunted growth, decreased chlorophyll content and canopy wilting [7, 11]. All this makes BCN an appropriate target for non-destructive phenotyping method development. Several remote sensing methods to detect stress caused by nematodes have already been successfully tested on a variety of crops such as potato, soybean or sugar beet. These methods are mainly based on imaging and non-imaging multi- and hyper spectral measurements, with the calculation of spectral vegetation indices (SVIs) [12,13,14].

Hillnhütter et al. [15] demonstrated the potential of normalized difference vegetation index (NDVI) to evaluate the symptoms caused by BCN on sugar beet plants under controlled conditions. Use of specific SVIs to predict the final beet yield and the nematode population in the soil has also been reported in field experiments [16]. Schmitz et al. [17] showed the ability of remote sensing thermography at field level to detect small changes in the canopy temperature of BCN-infested sugar beet. Thus, thermography and spectrometry appear to be suitable phenotyping methods for the detection of belowground symptoms caused by BCN. However, these systems require the use of expensive devices and complex data analysis methods.

Alternatively, visible imaging technology can be used for sugar beet phenotyping. Such a technology is cheaper than the aforementioned technologies, since it uses low cost sensors and the devices are easy to handle and calibrate [18]. The projected shoot area of the plants is usually calculated and used as a parameter to predict shoot biomass in different plant species [19,20,21,22]. Particularly in sugar beet, the use of visible images showed very promising results in discriminating, at an early plant developmental stage, BCN-infested and non-infested plants in the greenhouse [4]. In this study, the “digital canopy area” parameter calculated was a suitable proxy for shoot and root biomass estimation during the first 2 months of growth.

Beside the need to identify damage caused by nematodes and to evaluate the degree of infestation in the field, the use of phenotyping tools plays a role in agricultural research and development activities aiming at the discovery and development of new solutions for nematode control. In most of the cases, the evaluation studies aimed at evaluating the efficacy of the solutions by looking at the impact on the final yield. This implies that trials need to be kept up to harvest and last 3 months or longer. Using non-destructive measurements to get early insights regarding the efficacy of new solutions (compounds or cultivars) on the yield potential would allow to reduce the duration and costs of the trials, and to increase the testing cycles per year. Overall such new tools can have a substantial impact on the efficiency of compound screening or development of new cultivars.

The present study compares the ability of several traits (canopy area, canopy temperature and SVIs) obtained with three different phenotyping devices (visible imaging, thermography and spectrometry visible) to identify and characterize stress caused by BCN on sugar beet plants at the semi-field level. More specifically, the main objectives were to:

  1. (i)

    Compare the ability of the three phenotyping methods to detect stress generated by BCN on nematode susceptible and tolerant sugar beet plants,

  2. (ii)

    Evaluate the capability of the methods to predict sugar beet yield,

  3. (iii)

    Evaluate the potential of visible imaging to detect benefit of a contact nematicide.


Plant cultivation

Studies were conducted in 2014 and 2015 on a polytunnel area located in the Syngenta Research Centre in Stein (Switzerland). The area was equipped with a microplot system (Fig. 1a), which simulates real field conditions and allows to monitor the main environmental conditions. The experimental layout includes 70 microplots consisting of a pot in pot system. One 150 L plastic container (65 cm diameter and 60 cm depth) is nested inside of another, with both recessed in the ground up to the rim to reduce fluctuation of soil temperature.

Fig. 1
figure 1

a Overview of the semi-field platform with 70 microplots. b Top view of a microplot 416 °Cd. White dots represent the three “sowing locations”

The nematode susceptible cultivar Aimanta (Syngenta AG, Basel, Switzerland) and the nematode tolerant cultivar Bluefox (Syngenta AG, Basel, Switzerland) were used in 2014 and 2015 trials, respectively. The soil used was non-sterile sandy loam (56% sand, 31% silt, 11% clay, pH 7.7, 2% O.M.). A commercial seed treatment consisting of Thiram, Hymexazol, Thiamethoxame and Tefluthrine (6 + 14 + 60 + 8 gai/unit (one sugar beet unit is 100,000 seeds)) was applied to the seeds to avoid early insect attack and fungal disease. Six seeds were sown per microplot at three different locations (two seeds per location) (Fig. 1b). Two weeks after sowing, three seedlings were left in each microplot (one seedling kept per location) simulating a sowing density similar to the real sowing density adopted under field conditions (100,000 seeds per hectare).

Preparation of the soil and nematode inoculum

Cysts of H. schachtii were cultured at the Syngenta research centre in Stein. Cysts were coming from greenhouse oilseed rape plants cultured in loess soil. Infested soil was prepared by mixing the sandy loam soil with the amount of infested loess soil to reach a final level of 600 eggs and juveniles (J2) per 100 cm3 soil. Only the upper layer (corresponding to a volume of 40 L) of each microplot was infested.

In the 2015 experiment, an additional treatment including soil infested with cysts of H. schachtii + Fosthiazate nematicide (ISK Bioscience Corporation, Concord, OH, USA) was added [23]. Fosthiazate was applied as granules of Nemathorin 10G product (Syngenta AG, Basel, Switzerland) at the same time as the soil infestation with a final rate of 30 kg/ha. A randomized complete block design with ten replicates was used in the 2 years. Experimental settings and main crop management operations are reported in Table 1.

Table 1 Summary of the experimental settings and the crop management operations during the two microplot experiments

In both years, microplots were equipped with sensors to monitor air temperature ( and soil moisture ( Soil sensor technology is based on the microthermic measurements of soil moisture. Details of the environmental conditions are shown in Fig. 2. Thermal Time (TT) expressed in degree days (°Cd) was calculated using trait temperature as: TT = Σ if ≥ 0 \(\left({\frac{{T_{\max} + T_{\min}}}{2}} \right) - T_{base}\), with Tbase of 1.1 °C [24].

Fig. 2
figure 2

Evolution of daily average air temperature and soil moisture during 2014 and 2015 trials. Soil moisture is expressed as percentage of the field capacity (FC). The end of the lines corresponds to the respective harvest dates

In 2015, soil moisture during the first 40 days was higher than in 2014. In both years, soil moisture conditions were sufficient to allow a homogenous growth of the seedlings. It is worth noting that high air temperatures were observed during the sugar beet emergence (between 6 and 9 days after sowing) in 2015 (21 °C on average).

Evaluation of plants and nematodes

After harvest, fresh weight of beets was determined for each plant. Dry weight of the leaves was measured after a drying period of 72 h (70 °C). Final nematode population was assessed by sampling 1000 g of soil. Soil sampling was performed between five and 20 cm depth in the middle of each pot. All the soil samples were subsequently sent to an external lab (ClearDetection, Wageningen, NL) for analysis of the number of cysts and number of eggs and larvae per 100 cm3 of dry soil according to EPPO method 1/25 ( Plant growth stages (GS) were defined according to the BBCH scale [25].

Visible imaging

Canopy visible images were captured from seedling emergence up to 1300 °Cd every 2 or 3 days using a digital camera Canon S100 (Canon, Tokyo, Japan). The device was mounted on a mobile monopod and images were obtained from 1.8 m above plant canopy with a resolution of 0.0029 cm2 pixel−1. The monopod was held vertically in order to have the camera centred in the middle of the pot. To optimize image processing, photos were captured, when possible, under cloudy conditions early in the morning using the automatic settings of the camera. Fifteen minutes were necessary to acquire the 70 images which prevented any changes in the illumination and therefore also the necessity for white balancing. Raw pictures (Fig. 3a) were processed using ImageJ, the Java-based open-source image processing and analysis program (, following a workflow described by Joalland et al. [4] based on an image segmentation proposed by Woebbecke et al. [26]. This fast and non-invasive method allows to evaluate the “digital canopy area” (green area) at different times (Fig. 3b).

Fig. 3
figure 3

a Raw visible image taken from the top. b Image after processing 534 °Cd (susceptible cultivar Aimanta in 2014)


Thermal images were acquired using an infrared camera (Testo 885, Testo Ltd, UK). The thermal device was calibrated prior to taking pictures by setting up the emissivity to 96% and the reflected temperature compensation parameter to the current air temperature [27]. Pictures were then taken from the top of each pot in manual mode (autofocus off). Two images were automatically generated by the camera during the image acquisition; one thermal image, in which pixels correspond to temperature value, and one visible image (Fig. 4a, b). A macro was specifically built on ImageJ to extract the canopy temperature by combining both thermal and visible images.

Fig. 4
figure 4

a Visible and b thermal images obtained simultaneously with the thermal camera 856 °Cd in 2015

For each date and timing of measurement, the Vapour Pressure Deficit (VPD) was calculated according to the equation of Anderson [28] using air temperature and relative humidity. VPD reflects the ability of the air to hold water and it reflects the transpirational demand.


Data acquisition

Spectral measurements were performed several times during plant development (Table 2) with a non-imaging spectroradiometer (ASD FieldSpec® 4, Analytic Spectral Devices, Boulder, CO, USA) with a spectral range of 350–2500 nm. Spectra were acquired from the top of the microplots at a height of 1.20 m with a 25° field of view using a mobile dark box specifically designed for the microplot experiments. A 120 W halogen lamp (spotlight 120 W, Kent, Lyon, France) was used to provide constant optimum illumination of the canopy inside the box during the measurements. The dark box and halogen lamp combination made the conditions of measurement consistent between microplots and between days.

Table 2 Summary of phenotyping measurements during the two studies

Instrument optimization and reflectance calibration were performed using a Zenith Polymer® (SphereOptics, Germany) 99% reflectance target as white reference before the sample acquisition. Each sample scan represented an average of five reflectance spectra.

Spectral vegetation indices

For each date of measurement and each microplot, a selection of 123 published SVIs was computed to reduce the data dimension. SVIs were calculated using ratios of several bands at different ranges of the spectrum. For each measurement date, a correlation matrix was built for the 123 SVIs using control non-infested and control infested treatments. Indices highly inter-correlated to each other (Pearson’s correlation coefficient R > 0.8) were grouped resulting in 20 groups. One SVI was then selected by group which resulted in a final selection of 20 SVIs per date. For the final study, eight SVIs were selected out of 20 following a discriminant analysis between non-infested and nematode infested treatments to reflect the broad range of traits for which the SVIs were initially developed (Table 3).

Table 3 Selected SVIs, their respective equations, the aimed detection trait and references

Statistical data analysis

The program R [36] was used for analysis of the biological data. Beet fresh weight, leaf dry weight and canopy area of BCN infested and non-infested plants were tested for homogeneity of variance. They were then exposed to analysis of variance (ANOVA) at a probability level of 0.05 using the factor “nematode infestation”. Linear regression models were used to quantify the relationship between final beet fresh weight, nematode population and several phenotyping parameters. Regarding spectrometry data, a discriminant analysis was performed to identify, for each date of measurement, indices that allow to discriminate between control infested and non-infested treatments.


Plant fresh weight and nematode population

In both experiments, an artificial inoculation corresponding to 600 eggs and J2 per 100 cm3 of soil led to a moderate pressure similar to what can be expected in field situations. Such nematode pressure significantly affected the final beet fresh weight (Table 4). In 2014, beet biomass of nematode infested treatment was reduced by 32% compared to the non-infested treatment for the susceptible cultivar Aimanta, whereas in 2015, the final beet biomass reduction was 11% for the tolerant cultivar. Final average leaf dry biomass of infested plants (39.4 g plant−1) was significantly lower than that of non-infested plants (47 g plant−1) for the susceptible cultivar (−16%) (p < 0.05). This was not the case for the tolerant cultivar where no effect of BCN could be observed on the final leaf dry biomass.

Table 4 Effect of BCN on final beet fresh weight and leaf dry weight of sugar beet plants

As expected, almost no nematodes were found in the soil of non-infested treatments for both trials. Presence of negligible numbers of eggs and larvae can be explained by the non-sterile field soil that was used for these experiments. A larger number of nematodes was found in the infested pots in 2014 (on average 13,535 eggs and larvae per 100 cm3 soil) compared to 2015 (6950 eggs/larvae on average). The average Pf (Final nematode population)/Pi (Initial nematode population) ratio was 27 in 2014 and 11.5 in 2015.

Non-infested and infested treatments are displayed for 2014 and 2015 experiments

For both susceptible and tolerant cultivars, final aboveground biomass was strongly correlated with the belowground biomass (Fig. 5). Linear regression in 2014 and 2015 resulted in R2 of 0.82 and 0.74 respectively suggesting that leaf biomass is a good indicator of the beet biomass. The close relationship between above and belowground sugar beet biomass confirms the interesting use of non-destructive phenotyping tools to evaluate the status of the plant canopy over time. By measuring the canopy, the growth of the beet can be indirectly investigated.

Fig. 5
figure 5

Final leaf dry weight as a function of the final beet fresh weight (n = 381, R2 = 0.79, p < 0.01)

Early stress detection using visible imaging

Evolution of the canopy area of infested sugar beets presented similar patterns for susceptible and tolerant cultivars (Fig. 6). Canopy areas of both varieties were strongly affected by BCN during the first 600 °Cd. From 600 to 1000 °Cd, canopy area differences between infested and non-infested plants decreased due to a combination of leaf overlapping and plant recovering. After 1000 °Cd, differences between infested and non-infested treatments were not visible anymore using the canopy area parameter.

Fig. 6
figure 6

Evolution of the canopy area of infested susceptible and tolerant cultivars. Canopy area is expressed as a percentage of the non-infested treatment. Time periods where the difference in canopy area between infested and non-infested treatments are significant are represented on the figure (p < 0.05)

Canopy area allowed the detection of the nematode stress that was applied and the statistically significant discrimination between infested and non-infested plants from 230 °Cd (GS 13) to 880 °Cd (GS 33) for the nematode tolerant cultivar BlueFox and from 335 °Cd (GS 14) to 995 °Cd (GS 35) for the susceptible cultivar Aimanta (p < 0.05) (Fig. 6). This difference in the timing of the stress detection between the 2 years of trial was most likely caused by the low air temperatures during the first 2 weeks after sowing in 2014 (Fig. 2). This led to a slow and non-homogenous crop establishment and a low nematode pressure. The high variability in canopy area was reflected by the higher coefficient of variation observed (38% in 2014 and 18% in 2015). In summary, the visible imaging method was very sensitive in detecting the damaging effect of nematodes on aboveground plant growth already in very early stages (after the development of the third leaf 335 °Cd).

Canopy temperature evaluation using thermography

In 2014, canopy temperature differed significantly between treatments throughout the whole observed period (Table 5). Always, canopy temperature of nematode infested plants was significantly higher than the canopy temperature of the non-infested plants. At 1446 °Cd, average canopy temperature of the non-infested treatment was 18.6 versus 19.5 °C for the infested treatment. For the susceptible cultivar, differences between the two treatments increase with increasing VPD.

Table 5 Canopy temperature (°C) of non-infested and nematode infested treatments at different times during the season

In 2015, for the five dates of measurement, the canopy temperature of infested plants was on average 0.5 °C higher than the canopy temperature of the non-infested ones. However the difference was statistically significant only at 581 °Cd. There was no obvious correlation between the differences in canopy temperature of the two treatments and the VPD for the tolerant cultivar.

It can be stated that such canopy temperature differences observed between treatments are caused by nematode stress and not by the environmental condition variability on the platform. In fact, the randomized complete block design of the experiment was set up according to an air temperature gradient which prevented any effect of air temperature variability on the canopy temperature comparison between nematode infested and non-infested treatments. Soil moisture was similar for all the pots at each measurement date.

Nematode stress identification by a spectrometry approach

In Table 6, indices are grouped according to their relevance in assessing plant biomass, chlorophyll content, water status and general stress. Among the different indices, those related to the biomass, chlorophyll and general stress resulted in better detection of the nematode infestation and damage at the different stages of the crop development. In 2015, from 404 to 736 °Cd (GS 15 to GS 31), SVIs mainly related to plant biomass such as NDVI or leaf area such as MCARI2 were significantly affected by nematodes which confirmed the previous observation concerning the canopy area. At more advanced stages (GS 31 to GS 39), differences could be detected on both susceptible and tolerant cultivars using the CHLG and TGI respectively. The Health Index (HI), which was developed specifically for sugar beet, was particularly effective consistently across the 2 years of investigation and at different stages of crop development.

Table 6 Selection of SVIs that allowed to statistically discriminate non-infested and nematode infested treatments

In the last measurement (1618 °Cd), sugar beet plants displayed additional symptoms of general stresses with early leaf senescence which affected the identification of sole nematode effects.

Phenotyping parameters and final data

Most of the SVIs in Table 6 were significantly correlated with beet fresh weight and final nematode population in the soil. Correlation coefficients were always higher for the susceptible cultivar compared to the nematode tolerant cultivar (Table 7). On average NDVI, MCARI2 and CHLG were highly correlated with the beet fresh weight (R = 0.69) and the final nematode population (R = 0.78) for the susceptible cultivar. Cumulative canopy area which reflects the ability of the plants to absorb light over the season was significantly correlated with the final beet fresh weight in 2015 (R = 0.54, p < 0.1). In 2014, the weak correlation observed was not significant (R = 0.32). There was no significant correlation between cumulative canopy area and final BCN population in the soil. Canopy temperatures did not show any significant correlations with the final beet fresh weight and BCN population.

Table 7 Pearson’s correlation (R) between phenotyping variables at different dates and the fresh weight of the beet and final nematode population in the soil (n = 20, p < 0.1)

Practical application of visible imaging for nematicide research

In 2015, clear differences were observed in the canopy area between treatments during the first 35 days of plant development (Fig. 7). 244 °Cd (15 das), canopy area of fosthiazate treated plants was 29% higher than canopy area of the nematode infested treatment. At this date, average canopy areas of non-infested and fosthiazate treatments were statistically larger than the canopy area of nematode infested treatment (Fig. 8A). Evolution of canopy areas of non-infested and fosthiazate treated plants showed similar pattern. Both treatments showed significantly higher canopy area compared to the nematode-infested treatment from 244 to 560 °Cd (GS 16) (Fig. 7).

Fig. 7
figure 7

Canopy area of the non-infested and fosthiazate treatments as a percentage of the nematode infested treatment (n = 10). Only the first 35 days of growth are represented. From 244 to 560 °Cd both non-infested and fosthiazate treatments showed statistically significant higher canopy area than the nematode infested treatment (p < 0.05)

Fig. 8
figure 8

A Average canopy area and B final beet fresh weight of three treatments from 2015 trial (n = 10). Bars represent the standard error of the mean. Different letters indicate significant differences (p < 0.05, n = 10)

The same trend was observed in the final sugar beet yield. Fosthiazate treatment showed significant benefit in the final beet fresh weight compared to the nematode infested treatment (+14%) (Fig. 8B). Early canopy area differences reflected the final sugar beet yield.


The present study compares three phenotyping techniques in the same experimental settings across 2 years and two sugar beet cultivars. A multi-sensor approach was found to be of advantage for continuous crop phenotyping or monitoring in experimental and field settings [2, 37]. While more studies have been conducted under controlled greenhouse conditions [4, 15], this work was carried out outdoor by simulating conditions that are close to the real field situations comprising soil type, nematode infestation, plant density, duration of the crop cycle, plant canopy and root development. To the best of our knowledge, such a comparison of three phenotyping methods under outdoor conditions to characterize the sugar beet growth has not been published before.

During the 2 years of experimentation, the artificial nematode infestation successfully led to a yield reduction of 32% for the susceptible cultivar and 11% for the tolerant cultivar. These results are consistent with reports from other field and microplot studies [7, 38]. Thus the microplot settings used in this study were successful to simulate a realistic timing of nematode infestations and crop damage. BCN multiplication was 2.5 times higher for the susceptible cultivar compared to the tolerant one. This order of magnitude is consistent with respect to the definition of “nematode tolerance” given by Trudgill [39]. The fact that no differences could be observed in the final shoot biomass for the tolerant cultivar can be explained by the ability of the nematode tolerant cultivar to endure nematode damage and recover during the second part of the growing season [40, 41].

Visible imaging, thermography and spectrometry measurements enabled detection in a non-invasive, dynamic and objective manner of the effect of nematode infestation on sugar beet plants. Digital canopy area extracted from top-view visible images is a suitable tool to evaluate the effective plot-based canopy area. This parameter is taking account of different morphological components of the sugar beet such as the number of leaves, the area of the leaves and the plant architecture [42]. Canopy area appeared particularly suitable to dynamically characterize the early growth of the sugar beet plant from sowing to an advanced vegetative stage (GS 35). In a previous study carried out under greenhouse conditions, the “digital canopy area” parameter was identified as a proxy to estimate the shoot and root biomass of the sugar beets [4]. Such dynamic prediction of leaf biomass using visible images was also reported and used for high throughput phenotyping on other crops under greenhouse conditions [43, 44]. In the present study we demonstrated that the digital canopy area can be adapted to, and is effective in, outdoor conditions by looking at clusters of plants simulating the natural seed density expected in real field conditions. Overall, the top down visible imaging method showed its strength in the early evaluation of the degree of growth inhibition of the plant biomass. The nematode tolerant cultivar did not prevent BCN affecting the early plant development. Surprisingly, in the early growth stages, the canopy area reduction was higher for the tolerant cultivar compared to the susceptible cultivar which indicates that the tolerance mechanism does not prevent early nematode damage [45]. Benefits in very early plant growth (244 °Cd) observed with the use of fosthiazate showed the ability of a contact nematicide to protects the root by suppressing the first generation of J2 hatching from the cysts and to ensure yield benefit compared to the untreated plants [46].

Canopy temperature reflects plant water status, stomatal conductance and transpiration rate of the leaves [27, 47, 48]. It has been shown that nematodes strongly decrease water uptake of the roots which increases the stomatal resistance and consequently reduces the leaf evapotranspiration [49, 50]. In 2014, significantly higher canopy temperatures were observed for the nematode infested sugar beets compared to the non-infested plants. These results are consistent with previous observations made by Schmitz et al. [17] where a correlation between canopy temperature and nematode density was observed. Temperature difference between the two treatments increased with VPD. Infested susceptible plants had difficulties in cooling down their leaves when air conditions become constraining (VPD > 1.5). Most likely, infested plants were not able to keep up the high transpiration rate because of nematode damage at root level which compromised water uptake. The tolerant cultivar Bluefox behaved differently. Tolerance mechanisms allow sugar beet plants to maintain their transpiration rate even under high VPD.

Spectrometry measurements allowed the calculation of SVIs that reflected specific agronomical or physiological traits such as chlorophyll content, water content, biomass or photosynthesis rate [51,52,53,54]. The present study showed that specific SVIs allowed to differentiate between nematode infested and non-infested plants. Nematodes have an effect on different physiological parameters in both susceptible and tolerant cultivars. On the tolerant cultivar, most symptoms occur during the first 2 months of growth whereas on the susceptible cultivar, symptoms persist at more advanced growth stages since the plants are not able to recover from the infestation. The 2015 experiment helped to associate a type of BCN stress with the growth stages or time period where it occurs. The performance of indices related to the biomass and chlorophyll content was variable depending on the growth stages and among them MCARI2, 780/700 and TGI were the most promising. A close relationship has been demonstrated between the value of TGI index and the leaf chlorophyll content on a variety of crops [32, 55]. Such effect of BCN decreasing the leaf chlorophyll content was also reported by Schmitz et al. [11]. Nematode effect on the leaf water content was low in the tolerant cultivar which confirmed the limited effect of nematodes (also observed with thermography) in reducing transpiration rate on a tolerant cultivar. Two SVIs appeared suitable from early growth stages (GS 15) to advanced stages (GS 39) in detecting the stress caused by nematodes; HI and PRI. Health index (HI) uses two spectral regions centered on 700 nm and 534 nm. Reflectance near 700 nm is a feature of green vegetation and chlorophyll content whereas reflectance around 534 nm is an indicator of photosynthetic function [34, 52]. Thus, HI can be classified as a general stress index [1]. Photochemical reflectance index (PRI), based on reflectance at 531 and 570 nm, reflects the light use efficiency [56]. Although HI and PRI are not nematode specific, they appear suitable in detecting nematode stress over the whole season under semi-field conditions.

Correlations between SVIs and final sugar beet biomass demonstrated the ability of spectrometry in predicting final yield on both susceptible and tolerant cultivars. In particular, CHLG and MCARI2 were the best SVIs to predict final yield on susceptible cultivars. Close relationship between beet fresh weight and nematode incidence make the correlation between SVIs and BCN populations evident. Correlations were higher for the susceptible cultivar compared to the tolerant cultivar because of the larger range of beet fresh weight that was observed. Cumulative canopy area was also correlated with the final sugar beet yield suggesting a close relationship between the early plant growth and the final yield of the sugar beet. Such relationship between early phenotyping parameters and final yield is not so clear with other crops such as maize or wheat where the early plant growth does not always reflect the final plant yield as reported by Tekrony and Egli [57], Egli and Rucker [58] and Sankaran et al. [59]. In this respect, sugar beet appears a suitable crop for early yield prediction using phenotyping measurements.

Our results obtained with visible phenotyping showed that a slight delay in plant growth during the first 30 days had a significant effect on the final yield. Similar results were highlighted in a previous field study by sowing seeds of sugar beets at different timings to simulate a delay in the plant development [60]. The larger the canopy, the greater is the use of incident radiation. Olthof [61] observed higher damage on the plant when seeds were sown directly in nematode infested soil than when the infestation occurred 2 weeks after sowing. Early growth delay observed for the nematode infested plants could not be compensated during later growth of the crop. Thus, it appears crucial to avoid stress during the first growth stages of the sugar beet [62]. The fosthiazate effect in 2015 supports this point. In this experiment, fosthiazate nematicide was used as additional “positive control”. This nematicide acts by suppressing nematode hatching from the cysts and paralyzing juveniles and it is known to provide a strong root protection during the first month of the plant growth [46]. According to the rate that was applied (30 kg ha−1) and the concentration required for biological activity, it is likely that the fosthiazate effect in the soil stopped after 6–8 weeks [46, 63]. However, the early protection enabled a good development of the seedlings and insured yield benefit compared to the nematode infested plants. This result is of interest for crop protection research. Under moderate nematode pressure, protection of sugar beet plants against nematode damage should occur from sowing to 1200 °Cd (GS 25). Late nematode infestation did not significantly affect the plant growth of tolerant cultivars.

Given the complicated nature of the investigated nematode-plant interaction with belowground damage and unspecific symptoms displayed in the canopy, the phenotyping techniques evaluated here provided very encouraging results with potential applications in the area of sugar beet research. It can be stated that, visible imaging, thermography and spectrometry compared in the present investigation are complementary tools and are particularly suitable for automation. Field phenotyping platforms with multiple sensor systems will therefore be a valuable tool to improve crop performance via optimized management schedules [64, 65].


In conclusion, the study demonstrated that it was possible to use non-invasive and non-destructive technologies to characterize the dynamic of the plant growth and detect stress symptoms caused by BCN on nematode susceptible and tolerant sugar beets. While thermography only showed the ability to detect BCN stress on a susceptible cultivar, spectrometry and visible imaging technologies allowed the indirect observation of BCN damage on both susceptible and nematode tolerant cultivars and to give a prediction of the yield potential. In addition, the three different techniques have their specific strength at different points in time reflecting particular growth stages of the sugar beet. Visible imaging was the earliest stress indicator whereas spectrometry and thermography could identify the stress still when the canopy reached full coverage. Further applications of these tools could be developed for controlled environment and field situations. Under control conditions, canopy area has a great potential to be used as an early parameter to predict the degree of inhibition of the plant biomass caused by BCN and to quantify the degree of benefit from a new compound. Under field conditions visible image analysis, alone, may not be sufficiently specific to identify nematode damage because canopy area reduction can be caused by other types of stress. Therefore this technique would need to be combined with other approaches (e.g. spectrometry; thermography and/or soil sampling).



beet cyst nematode


days after sowing


growth stage


spectral vegetation indices


  1. Mahlein AK, Rumpf T, Welke P, Dehne HW, Plümer L, Steiner U, Oerke EC. Development of spectral indices for detecting and identifying plant diseases. Remote Sens Environ. 2013;128:21–30.

    Article  Google Scholar 

  2. Liebisch F, Kirchgessner N, Schneider D, Walter A, Hund A. Remote, aerial phenotyping of maize traits with a mobile multi-sensor approach. Plant Methods. 2015;11(1):9.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Walter A, Liebisch F, Hund A. Plant phenotyping: from bean weighing to image analysis. Plant Methods. 2015;11(1):1–11.

    Article  Google Scholar 

  4. Joalland S, Screpanti C, Gaume A, Walter A. Belowground biomass accumulation assessed by digital image based leaf area detection. Plant Soil. 2016;398(1):257–66.

    Article  CAS  Google Scholar 

  5. Luc M, Sikora RA, Bridge J. Plant parasitic nematodes in subtropical and tropical agriculture. Egham: CABI Bioscience; 2005.

    Book  Google Scholar 

  6. Seinhorst JW. The relation between nematode density and damage to plants. Nematologica. 1965;11:137–54. doi:10.1163/187529265X00582.

    Article  Google Scholar 

  7. Cooke DA. Beet cyst nematode (Heterodera schachtii Schmidt) and its control on sugar beet. Agric Zool Rev. 1987;2:135–83.

    Google Scholar 

  8. Evans K, Webster RM, Halford PD, Barker AD, Russell MD. Site-specific management of nematodes—pitfalls and practicalities. J Nematol. 2002;34:194–9.

    PubMed  PubMed Central  Google Scholar 

  9. Wyse-Pester DY, Wiles LJ, Westra P. The potential for mapping nematode distributions for site-specific management. J Nematol. 2002;34:80–7.

    PubMed  PubMed Central  Google Scholar 

  10. Hillnhütter C, Schweizer A, Kühnhold V, Sikora RA. Remote sensing for the detection of soil-borne plant parasitic nematodes and fungal pathogens. In: Oerke EC, Gerhards R, Menz G, Sikora RA, editors. Precision crop protection—the challenge and use of heterogeneity. Dordrecht: Springer; 2010. p. 151–65.

    Chapter  Google Scholar 

  11. Schmitz A, Tartachnyk II, Kiewnick S, Sikora RA, Kühbauch W. Detection of Heterodera schachtii infestation in sugar beet by means of laser-induced and pulse amplitude modulated chlorophyll fluorescence. Nematology. 2006;8:273–86.

    Article  Google Scholar 

  12. Heath WL, Haydock PPJ, Wilcox A, Evans K. The potential use of spectral reflectance from the potato crop for remote sensing of infection by potato cyst nematodes. Asp Appl Biol. 2000;60:185–8.

    Google Scholar 

  13. Laudien R. Entwicklung eines GIS-gestützten schlagbezogenen Führungsinformationssystems für die Zuckerwirtschaft (Development of a field- and GIS-based management information system for the sugar beet industry). Ph.D. thesis University of Hohenheim. Germany. 2005.

  14. Nutter FW, Tylka GL, Guan J, Moreira AJD, Marett CC, Rosburg TR, et al. Use of remote sensing to detect soybean cyst nematode-induced plant stress. J Nematol. 2002;34:222–31.

    CAS  PubMed  PubMed Central  Google Scholar 

  15. Hillnhütter C, Mahlein AK, Sikora RA, Oerke EC. Use of imaging spectroscopy to discriminate symptoms caused by Heterodera schachtii and Rhizoctonia solani on sugar beet. Precis Agric. 2012;13:17–32. doi:10.1007/s11119-011-9237-2.

    Article  Google Scholar 

  16. Hillnhütter C, Mahlein AK, Sikora RA, Oerke EC. Remote sensing to detect plant stress induced by Heterodera schachtii and Rhizoctonia solani in sugar beet fields. Field Crop Res. 2011;122(1):70–7.

    Article  Google Scholar 

  17. Schmitz A, Kiewnick S, Schlang J, Sikora RA. Use of high resolution digital thermography to detect Heterodera schachtii infestation in sugar beets. Commun Agric Appl Biol Sci. 2004;69(3):359–63.

    CAS  PubMed  Google Scholar 

  18. Li L, Zhang Q, Huang D. A review of imaging techniques for plant phenotyping. Sensors. 2014;14:20078–111.

    Article  PubMed  PubMed Central  Google Scholar 

  19. Sher-Kaul S, Oertli B, Castella E, Lachavanne JB. Relationship between biomass and surface area of six submerged aquatic plant species. Aquat Bot. 1995;51:147–54. doi:10.1016/0304-3770(95)00460-H.

    Article  Google Scholar 

  20. Smith SM, Garrett PB, Leeds JA, McCormick PV. Evaluation of digital photography for estimating live and dead aboveground biomass in monospecific macrophyte stands. Aquat Bot. 2000;67:69–77. doi:10.1016/S0304-3770(99)00085-6.

    Article  Google Scholar 

  21. Mizoue N, Masutani T. Image analysis measure of crown condition, foliage biomass and stem growth relationships of Chamaecyparis obtusa. For Ecol Manage. 2003;172:79–88. doi:10.1016/S0378-1127(02)00281-5.

    Article  Google Scholar 

  22. Tackenberg O. A new method for non-destructive measurement of biomass, growth rates, vertical biomass distribution and dry matter content based on digital image analysis. Ann Bot Lond. 2007;99:777–83.

    Article  Google Scholar 

  23. Koyanagi T, Imai O, Yoshida K. Development of a new nematicide, fosthiazate. J Pestic Sci. 1998;23:174–83.

    Article  CAS  Google Scholar 

  24. Holen CD, Dexter AGA. Growing degree day equation for early sugarbeet leaf stages. Res Ext Rep. 1996;27:152–7.

    Google Scholar 

  25. Meier U, Bachmann E, Buhtz H, Hack H, Klose R, Märländer B, Weber E. Phenological growth stages of beta beets (Beta vulgaris L. spp.). Nachrichtenbl Deut Pflanzenschutzd. 1993;45:37–41.

    Google Scholar 

  26. Woebbecke DM, Meyer GE, Von Bargen K, Mortensen DA. Color indices for weed identification under various soil, residue and lighting conditions. Trans ASAE. 1995;38:259–69.

    Article  Google Scholar 

  27. Oerke EC, Steiner U. Potential of digital thermography for disease control. In: Oerke EC, Gerhards R, Menz G, Sikora RA, editors. Precision crop protection—the challenge and use of heterogeneity. Dordrecht: Springer; 2010. p. 167–82.

    Chapter  Google Scholar 

  28. Anderson DB. Relative humidity or vapor pressure deficit. Ecology. 1936;17(2):277–82.

    Article  Google Scholar 

  29. Rouse JW, Haas RH, Schell JA, Deering DW. Monitoring vegetation systems in the Great Plains with ERTS. NASA Special Publ. 1974;351:309.

    Google Scholar 

  30. Haboudane D, Miller JR, Pattey E, Zarco-Tejada PJ, Strachan IB. Hyperspectral vegetation indices and novel algorithms for predicting green LAI of crop canopies: modeling and validation in the context of precision agriculture. Remote Sens Environ. 2004;90(3):337–52.

    Article  Google Scholar 

  31. Mistele B, Gutser R, Schmidhalter U, Mulla DJ. Validation of field-scaled spectral measurements of the nitrogen status in winter wheat. In: Proceedings of the 7th international conference on precision agriculture and other precision resources management, Hyatt Regency, Minneapolis, MN, USA, 25–28 July 2004. Precision Agriculture Center, University of Minnesota, Department of Soil, Water and Climate, p. 1187–95.

  32. Hunt ER, Daughtry CST, Eitel JU, Long DS. Remote sensing leaf chlorophyll content using a visible band index. Agron J. 2011;103(4):1090–9.

    Article  Google Scholar 

  33. Gitelson AA, Keydan GP, Merzlyak MN. Three-band model for noninvasive estimation of chlorophyll, carotenoids, and anthocyanin contents in higher plant leaves. Geophys Res Lett. 2006;33(11):L11402.

    Article  Google Scholar 

  34. Gamon JA, Penuelas J, Field CB. A narrow-waveband spectral index that tracks diurnal changes in photosynthetic efficiency. Remote Sens Environ. 1992;41:35–44.

    Article  Google Scholar 

  35. Clay DE, Kim KI, Chang J, Clay SA, Dalsted K. Characterizing water and nitrogen stress in corn using remote sensing. Agron J. 2006;98(3):579–87.

    Article  Google Scholar 

  36. R Development Core Team. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. 2008. ISBN 3-900051-07-0.

  37. Pfeifer J, Khanna R, Constantin D, Popovic M, Galceran E, Kirchgessner N, Walter A, Siegwart R, Liebisch F. Towards automatic UAV data interpretation for precision farming. In: CIGR-AgEng conference. Aarhus, Denmark. 2016.

  38. Herr L. Sugar beet diseases incited by Rhizoctonia spp. In: Sneh B, Jabaji-Hare S, Neate S, Dijst G, editors. Rhizoctonia species: taxonomy, molecular biology, ecology, pathology and disease control. Dordrecht: Springer; 1996. p. 341–9. doi:10.1007/978-94-017-2901-7_31.

    Chapter  Google Scholar 

  39. Trudgill DL. Resistance to and tolerance of plant parasitic nematodes in plants. Ann Rev Phytopathol. 1991;29(1):167–92.

    Article  Google Scholar 

  40. Wallace HR. A perception of tolerance. Nematologica. 1988;33:419–32.

    Article  Google Scholar 

  41. Davy de Virville J, Person-Dedryver F. Growth and respiratory activity of roots of various Triticeae tolerant or resistant to Heterodera avenae Woll. with or without infection by the nematode. Rev Nematol. 1989;12:379–86.

    Google Scholar 

  42. Milford GFJ, Pocock TO, Riley J. An analysis of leaf growth in sugar-beet. I. Leaf appearance and expansion in relation to temperature under controlled conditions. Ann Appl Biol. 1985;106:163–72. doi:10.1111/j.1744-7348.1985.tb03106.x.

    Article  Google Scholar 

  43. Hartmann A, Czauderna T, Hoffmann R, Stein N, Schreiber F. HTPheno: an image analysis pipeline for high throughput plant phenotyping. BMC Bioinform. 2011;12:148.

    Article  Google Scholar 

  44. Golzarian MR, Frick RA, Rajendran K, Berger B, Roy S, Tester M, Lun DS. Accurate inference of shoot biomass from high-throughput images of cereal plants. Plant Methods. 2011;7:1–11.

    Article  Google Scholar 

  45. Trudgill DL. Mechanisms of damage and of tolerance in nematode infested plants. Nematology from molecule to ecosystem. Dundee: European Society of Nematologists Invergourie; 1992. p. 133–45.

    Google Scholar 

  46. Woods SR, Haydock PJ, Edmunds C. Mode of action of fosthiazate used for the control of the potato cyst nematode Globodera pallida. Ann Appl Biol. 1999;135(1):409–15.

    Article  CAS  Google Scholar 

  47. Inoue Y, Kimball BA, Jackson RD, Pinter PJ, Reginato RJ. Remote estimation of leaf transpiration rate and stomatal resistance based on infrared thermometry. Agric For Meteorol. 1990;51(1):21–33.

    Article  Google Scholar 

  48. Cohen Y, Alchanatis V, Meron M, Saranga Y, Tsipris J. Estimation of leaf water potential by thermal imagery and spatial analysis. J Exp Bot. 2005;56:1843–52.

    Article  CAS  PubMed  Google Scholar 

  49. Haverkort AJ, Fasan T, Van de Waart M. The influence of cyst nematodes and drought on potato growth. 2. Effects on plant water relations under semi-controlled conditions. Neth J Plant Pathol. 1991;97(3):162–70.

    Article  Google Scholar 

  50. Jones HG. Application of thermal imaging and infrared sensing in plant physiology and ecophysiology. Adv Bot Res. 2004;41:107–63.

    Article  Google Scholar 

  51. Curran PJ, Dungan JL, Macler BA, Plummer SE. The effect of a red leaf pigment on the relationship between red edge and chlorophyll concentration. Remote Sens Environ. 1991;35(1):69–76.

    Article  Google Scholar 

  52. Gitelson AA, Merzlyak MN. Signature analysis of leaf reflectance spectra: algorithm development for remote sensing of chlorophyll. J Plant Physiol. 1996;148(3):494–500.

    Article  CAS  Google Scholar 

  53. Mahlein AK, Oerke EC, Steiner U, Dehne HW. Recent advances in sensing plant diseases for precision crop protection. Eur J Plant Pathol. 2012;133(1):197–209.

    Article  CAS  Google Scholar 

  54. Liebisch F, Küng G, Damm A, Walter A. Characterization of crop vitality and resource use efficiency by means of combining imaging spectroscopy based plant traits. In: Workshop on hyperspectral image and signal processing: evolution in remote sensing, 24–27 June 2014, Lausanne, Switzerland, IEEE International. p. 6.

  55. Constantin D, Rehak M, Akhtman Y, Liebisch F. Detection of crop properties by means of hyperspectral remote sensing from a micro UAV. 21. In: Workshop Computer-Bildanalyse in der Landwirtschaft and 3. Workshop Unbemannte autonom fliegende Systeme (UAS) in der Landwirtschaft. 6. und 7. Mai 2015, Braunschweig, Germany.

  56. Trotter GM, Whitehead D, Pinkney EJ. The photochemical reflectance index as a measure of photosynthetic light use efficiency for plants with varying foliar nitrogen contents. Int J Remote Sens. 2002;23(6):1207–12.

    Article  Google Scholar 

  57. Tekrony DM, Egli DB. Relationship of seed vigor to crop yield: a review. Crop Sci. 1991;31(3):816–22.

    Article  Google Scholar 

  58. Egli DB, Rucker M. Seed vigor and the uniformity of emergence of corn seedlings. Crop Sci. 2012;52(6):2774–82.

    Article  Google Scholar 

  59. Sankaran S, Khot LR, Espinoza CZ, Jarolmasjed S, Sathuvalli VR, Vandemark GJ, Miklas PN, Carter AH, Pumphrey MO, Knowles NR, Pavek MJ. Low-altitude, high-resolution aerial imaging systems for row and field crop phenotyping: a review. Eur J Agron. 2015;70:112–23.

    Article  Google Scholar 

  60. Scott RK, English SD, Wood DW, Unsworth MH. The yield of sugar beet in relation to weather and length of growing season. J Agric Sci. 1973;81:339–47. doi:10.1017/S0021859600059001.

    Article  Google Scholar 

  61. Olthof HA. Effect of plant age and transplanting damage on sugar beets infected by Heterodera schachtii. J Nematol. 1983;15:555–9.

    CAS  PubMed  PubMed Central  Google Scholar 

  62. Griffin GD. The relationship of plant age, soil temperature, and population density of Heterodera schachtii on the growth of sugarbeet. J Nematol. 1981;13(2):184–90.

    CAS  PubMed  PubMed Central  Google Scholar 

  63. Pantelelis I, Karpouzas DG, Menkissoglu-Spiroudi U, Tsiropoulos N. Influence of soil physicochemical and biological properties on the degradation and adsorption of the nematicide fosthiazate. J Agric Food Chem. 2006;54(18):6783–9.

    Article  CAS  PubMed  Google Scholar 

  64. Kirchgessner N, Liebisch F, Yu K, Pfeifer J, Friedli M, Hund A, Walter A. The ETH field phenotyping platform FIP-A cable-suspended multi-sensor system. Funct Plant Biol. 2016;44:154–68.

    Article  Google Scholar 

  65. Virlet N, Sabermanesh K, Sadeghi-Tehran P, Hawkesford M. Field scanalyser: an automated robotic field phenotyping platform for detailed crop monitoring. Funct Plant Biol. 2016;44:143–53.

    Article  Google Scholar 

Download references

Authors’ contributions

SJ, CS, HVV, FL, AG and AW designed the experiments. SJ performed the trials and the phenotyping measurements with the support of FL. SJ and HVV did the statistical analysis. SJ and CS drafted the manuscript supported by FL and AW. All authors read and approved the final manuscript.


We would like to thank Dr. Roger Hall for his careful review of the manuscript and for his valuable suggestions and corrections as well as Beat Reber and Tobias Plec for the technical support in setting up the trials. We also thank Brigitte Slaats and Tobias Straumann for growing and providing the nematode inoculum and Cliff Watrin and Rita Kuznia for their technical advice through the project phases.

Competing interests

The authors declare that they have no competing interests.

Availability of data and materials

The datasets generated and analysed during the current study are not publicly available because part of the study includes sensitive information related to industrial project, but are available from the corresponding author on reasonable request.

Consent for publication

All authors read and approved the final manuscript.

Ethics approval and consent to participate

Not applicable.


Not applicable.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Claudio Screpanti.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Joalland, S., Screpanti, C., Liebisch, F. et al. Comparison of visible imaging, thermography and spectrometry methods to evaluate the effect of Heterodera schachtii inoculation on sugar beets. Plant Methods 13, 73 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: