Next Article in Journal
Influence of Cactus Pear Mucilage-Based Edible Coating on Marketability and Edibility Parameters of Minimally Processed Loquat Fruits
Next Article in Special Issue
Evidence of Resistance to QoI Fungicides in Contemporary Populations of Mycosphaerella fijiensis, M. musicola and M. thailandica from Banana Plantations in Southeastern Brazil
Previous Article in Journal
Use of Active Sensors in Coffee Cultivation for Monitoring Crop Yield
Previous Article in Special Issue
Efflux Pumps and Multidrug-Resistance in Pyricularia oryzae Triticum Lineage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Fungicide Application Timing Based on Soybean Rust Prediction Model on Application Technology and Disease Control

by
Matheus Mereb Negrisoli
1,*,
Flávio Nunes da Silva
1,
Raphael Mereb Negrisoli
1,
Lucas da Silva Lopes
1,
Francisco de Sales Souza Júnior
1,
Bianca Rezende de Freitas
2,
Edivaldo Domingues Velini
1 and
Carlos Gilberto Raetano
1
1
Department of Plant Protection, School of Agriculture, Sao Paulo State University, 3780 Avenida Universitária, Botucatu 01049-010, SP, Brazil
2
Department of Crop Science, School of Agriculture, Sao Paulo State University, 3780 Avenida Universitária, Botucatu 01049-010, SP, Brazil
*
Author to whom correspondence should be addressed.
Agronomy 2022, 12(9), 2119; https://doi.org/10.3390/agronomy12092119
Submission received: 12 July 2022 / Revised: 3 August 2022 / Accepted: 8 August 2022 / Published: 7 September 2022

Abstract

:
The application of remote sensing techniques and prediction models for soybean rust (SBR) monitoring may result in different fungicide application timings, control efficacy, and spraying performance. This study aimed to evaluate the applicability of a prediction model as a threshold for disease control decision-making and to identify the effect of different application timings on SBR control as well as on the spraying technology. There were two experimental trials that were conducted in a 2 × 4 factorial scheme: 2 cultivars (susceptible and partially resistant to SBR); and four application timings (conventional chemical control at a calendarized system basis; based on the prediction model; at the appearance of the first visible symptoms; and control without fungicide application). Spray deposit and coverage at each application timing were evaluated in the lower and upper region of the soybean canopy through quantitative analysis of a tracer and water-sensitive papers. The prediction model was calculated based on leaf reflectance data that were collected by remote sensing. Application timings impacted the application technology as well as control efficacy. Calendarized system applications were conducted earlier, promoting different spray performances. Spraying at moments when the leaf area index was higher obtained poorer distribution. None of the treatments were capable of achieving high spray penetration into the canopy. The partially resistant cultivar was effective in holding disease progress during the crop season, whereas all treatments with chemical control resulted in less disease impact. The use of the prediction model was effective and promising to be integrated into disease management programs.

1. Introduction

Soybean rust (SBR) causes significant crop yield losses throughout the world [1,2]. The disease is caused by Phakopsora pachyrhizi Sydow and is controlled mostly by fungicide application at pre-determined scheduled timings of the soybean growth stages, usually regardless of disease incidence and pressure level [1,3]. Due to constant fungicide spraying over the seasons and several times at the same season, a large number of fungi populations that are resistant to different chemical groups of fungicides have also been reported [4,5], significantly reducing fungicide efficacy over time [6].
Most soybean cultivation is conducted in extensive field areas, which hardens disease monitoring conductance. Nowadays, specialists monitor SBR in the field through extensive scouting that is based on disease symptom recognition. However, monitoring is usually absent, and farmers choose to spray on a calendarized system basis (i.e., at a pre-determined period) as a guarantee of crop yield. Besides input losses, these practices can lead to fungicide resistance selection pressure, poor spraying quality at periods that are not appropriate to achieve the best spray distribution, as well as possible wrong timing, which also leads to lower efficacy [1,6,7,8,9]. Other techniques are being applied to the integrated disease management (IDM) of the soybean–P. pachyrhizi pathosystem, such as resistant or partially resistant cultivars [10,11], planting date restrictions [1,2], and biological control [12]. Nonetheless, disease monitoring is considered the basis of any control method that is applied.
One of the alternatives to aid in disease control decision-making is the use of remote sensing technologies to facilitate data acquisition in wider fields as well as the obtainment of reliable information. It means that these technologies can precisely detect the disease remotely and at a faster pace since it allows the monitoring of wider fields depending on the sensor and where it is based [13,14,15]. Studies have reported the possibility of identifying SBR by remote sensing techniques [16,17,18], as well as for other diseases and crops [15,19]. Moreover, researchers have been applying all of this information in data modeling so that remote sensing could serve as the background information in the construction of disease prediction models as decision support systems [14,20]. These models are expected to improve the application timing, control effectiveness, as well as spraying performance regarding its uniform distribution in the crop.
Fungicide application timing strongly influences disease control, and it may vary depending on the fungicide mode of action and translocation capacity [8,21]. Besides reaching the target at the moment of most fungicide susceptibility, the application timing may influence the application technology in terms of product distribution, coverage, and penetration into the canopy [9,22]. Different timings also represent variations in soybean canopy density, especially due to leaf area index (LAI) and vegetation density [23] which play a major role as a barrier to the spray reaching the interior of the canopy. Most fungicide applications target the lower region of the soybean canopy since this is where SBR starts its development [1]. Therefore, applications at moments of higher LAI tend to have reduced penetration capacity and worse spray distribution [24].
According to Müller et al. [8], effective monitoring along with fungicide application as soon as the disease is identified in the field is a key factor to mitigate excessive and unnecessary application. However, since most systemic fungicides that are used for SBR control have both curative and preventive action modes, it is unknown how the application timing will affect control efficacy. The hypothesis is that SBR monitoring through prediction models that are based on remote sensing can help identify the first appearance of the disease in the field, and, thereby, improve application timing accuracy. Moreover, different application timings will potentially influence product distribution across soybean canopy regions.
The goal of this study was to evaluate the applicability of prediction models as the threshold for disease control decision-making and to identify the effect of different application timings on soybean rust control as well as on the spraying technology. This study is a continuation of the research that was conducted by Negrisoli et al. [18], using a prediction model to determine the fungicide application timing and how it differs from conventional methods.

2. Materials and Methods

There were two experiment replications that were carried out in the field over the 2020/2021 crop season, in different experimental areas of Botucatu, SP, Brazil (Field 1: 22°48′48″ S and 48°25′37″ W; Field 2: 22°49′38″ S and 48°25′40″ W) (Figure 1). In both experimental fields, the no-tillage system was adopted, and all sowing operations, phytosanitary management, and evaluations were carried out homogeneously.
The trials were sown on 8 December 2020, spaced at 0.45 m between rows, and with a population of 196,528 plants ha−1 in Field 1 and 211,806 plants ha−1 in Field 2. The crop was fertilized with 250 kg ha−1 of the commercial formula 02-20-20 (NPK) homogeneously throughout the experimental area.
The experiments were carried out in a randomized complete block design and the treatments were distributed in a 2 × 4 factorial scheme: 2 cultivars (Brevant DS6217 IPRO, susceptible to SBR; TMG 7063 IPRO, partially resistant to SBR); and 4 application timings that were based on different parameters as decision-making of control: conventional chemical control, spraying at a scheduled basis at the soybean reproductive stage R1 and R1 + 15 days; application timing (A1) defined based on the disease prediction model and A1 + 15 days; application timing (A2) at the appearance of the first visible symptoms and A2 + 15 days; and control treatments without fungicide application (Table 1), and all in four repetitions. The descriptions of each cultivar are shown in Table 2.
The definition of the application timings “A1” and “A2” were based on weekly disease severity monitoring starting at the vegetative growth stage, V6, in each plot. Firstly, the definition of A1 was conducted using an SBR prediction model that was proposed by Negrisoli et al. [18], based on the leaf reflectance, as the decision support system of control, as described subsequently herein.

2.1. Leaf Reflectance Assessment

This assessment was conducted every five days after the V6 vegetative growth stage, by randomly collecting 5 leaflets per plot from the lower region of the canopy for leaf reflectance analysis. The UV 2700 non-imaging spectrophotometer (Shimadzu, Kyoto, Japan) was used coupled to an integrated base ISR-603: Integrating Sphere Attachment, analyzing the reflectance of each leaflet in the range of 270 to 1000 nm with an interval of 3.0 nm, as described by Negrisoli et al. [18]. Each sample (leaflet) was evaluated separately, generating the reflectance values of each sample in the previously mentioned spectral range interval.
The spectral curves were reduced from 270–1000 nm to 270–900 nm for high noise levels reduction [27] and the Savitzky–Golay filter was applied using 11 central points and a third-degree polynomial [28,29]. The data were used to calculate a list of 19 vegetation indices (Vis) that was representative of the disease effect on the crop and to allow disease severity classification and, therefore, to predict or detect diseased plants among the samples [16,18,30]. Vegetation indices were chosen to be used instead of the full spectra length so that other spectral sensors besides the hyperspectral ones could be used to acquire the data that were required to run this model [18], besides being able to supply valuable information for disease detection and plant stress identification [14].
The Vis were calculated for each sample at every reflectance evaluation date and this database was used to supply the prediction model that was based on the Support Vector Machine (SVM) algorithm, which is programmed to classify into four classes: “healthy”, “low severity”, “moderate severity”, and “high severity”. At the moment when the samples from susceptible (T5) and partially resistant (T6) cultivars that averaged 1 plant per plot were classified at “low severity” (diseased plant), the application was thereafter immediately conducted. All the Vis calculations, evaluation methodologies, and prediction model construction are fully described by [18]. The formulas and list of Vis that were used are also fully described by Negrisoli et al. [18].
For the definition of “A2” application timing, 10 leaflets were collected from the lower region of the canopy in each plot for visual assessment of the symptoms. The samples were taken to a laboratory for better visualization of the fungal structures and symptoms under a stereoscopic microscope. Susceptible and resistant reactions may result in different SBR symptoms. A susceptible reaction results in tan to light-brown lesions (TAN reaction type) and the presence of uredinia throughout the leaf, while in partially resistant cultivars the lesions are characterized by reddish-brown lesions (RB-reaction type) and a lower quantity of uredinia [3,11]. The spraying was conducted at the time when the first SBR symptoms were detected in the plots of susceptible (T7) and partially resistant cultivars (T8).
The standard application timing was conducted at the end of the vegetative growth and the beginning of the reproductive stage (R1) of susceptible (T3) and partially resistant cultivars (T4). This application timing is common throughout the country and was used as the standard timing parameter.

2.2. Fungicide Sprayings

The fungicide Elatus®, Syngenta (azoxystrobin 60 g L−1 + benzovindiflupyr 30 g L−1) was used for SBR control at a dose of 0.250 kg ha−1, following the manufacturer’s recommendations. The spraying was carried out with a CO2 pressurized backpack sprayer with a 2.0-m boom that was equipped with four flat fan spray nozzles (Teejet XR110-02), displacement speed of 5 km h−1, and a working pressure of 200 kPa, providing a spraying volume of 150 L ha−1.

2.3. Quali-Quantitative Analysis of Spraying

The spray deposit capacity as well as the percentage of coverage of each treatment were evaluated in both first and second applications (Table 1). For the quantitative analysis of the deposit, the food dye Brilliant Blue marker (Duas Rodas Industrial LTDA, São Bernardo do Campo, Brazil) was applied immediately before fungicide spraying to avoid the influence of the fungicide formulation on the optical reading of the marker by spectrophotometry. The Brilliant Blue marker dye was solubilized in distilled water at a concentration of 3000 mg L−1. After spraying, three samples that were composed of 10 leaflets each were randomly collected from the upper and lower region of the plant canopy in the central lines of each plot, also evaluating the ability of the treatments to penetrate the crop canopy.
The samples were taken to the laboratory and processed according to Palladini et al. [31]. Each sample received 100 mL of distilled water, stirred for 15 s, and the resulting solution was transferred to plastic containers. The quantification of deposits was performed in a spectrophotometer (Shimadzu VIS 1601 PC) with an absorbance reading at a wavelength of 630 nm [31]. After the tracer was washed, the leaf area of each sample was evaluated with the aid of a benchtop leaf area meter (LICOR, model LI-3100). The readings of known concentrations of the dye were correlated to the absorbance values that were obtained in the spectrophotometer and the calibration curve was constructed, obtaining the dye concentration in the target in mg L−1. Finally, the volume that was found in the target was established by correlating the dye concentration in the samples washing solution with the dye concentration in the final spray solution, presenting the data in μL cm−2.
The qualitative analysis of the spraying was carried out through the evaluation of the spray coverage using water-sensitive papers (WSP) (26 × 76 mm) that were distributed in each plot. There were two WSPs that were used attached to a metal rod in the central lines of the plots, one located horizontally at the top and the other at the bottom part of the canopy. Spray coverage and deposit evaluations were performed at the same time.
After spraying, the WSPs were placed in Petri dishes to prevent moisture absorption and taken to the laboratory for analysis. These samples were scanned at a resolution of 600 dpi and the digitalized images were analyzed by the “GOTAS” software (Embrapa®) to obtain the percentage of surface coverage.
The leaf area index (LAI) of the plots was also measured in the period between 50 to 70 days after emergence (DAE) (Field 1) and 50–88 DAE (Field 2), corresponding to the range of all application timings in each field. Therefore, at 50, 60, and 70 DAE of Field 1, and 50, 65, and 88 DAE of Field 2, 10 plants of each cultivar were randomly collected across each experimental field trial. The whole plants were taken to the laboratory, completely defoliated, and the leaves were measured by a benchtop leaf area meter (LICOR, model LI-3100). The total leaf area that was measured of each sample was converted to m2 and divided by the mean number of plants m−1 (Field 1 = 9.93 plants m−1 of the susceptible cultivar; 7.75 plants m−1 of the partially resistant cultivar; Field 2 = 10.34 plants m−1 of the susceptible cultivar; 8.72 plants m−1 of the partially resistant cultivar) [32]. The LAI is an important factor to understand spraying quality behavior as well as the capacity of each operational parameter according to the leaf density as a barrier to spraying.

2.4. Assessment of Disease Severity and Control Efficacy

Starting at the V6 growth stage, 10 leaflets were collected weekly per plot from the lower third of the plant and taken to the laboratory for better visualization of fungal structures and symptoms under a stereoscopic microscope. The disease severity (%) was estimated based on the SBR diagrammatic scale that was proposed by Godoy et al. [33], based on visual observation of the symptoms. The Area Under the Disease Progress Curve (AUDPC) was calculated according to Campbell and Madden [34] using the mean values of disease severity that were obtained in the plots and on the respective evaluation dates.
In addition, the disease severity was indirectly assessed by the defoliation level. The evaluations started at 85 DAE, the moment when the highest level of severity began, and no defoliation was still detected. The evaluation was carried out using the ASD FieldSpec Dual Spectroradiometer (Analytical Spectral Devices, Boulder, CO, USA), with a spectral range from 350 to 1070 nm and 7.5° of field of view. A white panel with approximately 100% reflectance (Spectralon) was used as a reference for calibration. Calibration and optimization of the equipment were conducted every 10 min.
The evaluations were conducted weekly under intense sunlight at around 11:00 a.m., performing three readings in the central region of each plot in each experimental trial, maintaining uniform height and equipment position. The reflectance data that were obtained were smoothed by Savitsky–Golay filtering as described above. The reflectance database was then used to calculate the VIs that was the most correlated to the effect of SBR on crop defoliation. For that, the calculation of the LAI through the normalized difference vegetation index (NDVI) and Beer–Lambert law [35] stood out and was used as a reference for evaluating crop defoliation by the disease. Tan et al. [35] proposed the calculation of the LAI through the integration of the Beer–Lambert Law and the NDVI of the samples, also considering its leaf orientation values. Here, the soybean crop was considered as a middle-type plant with leaf orientation values ranging from 30 to 60°, and the whole model is fully described by the authors [35].

2.5. Evaluation of the Effect of Control on Crop Yield

At the end of the crop season, each plot was individually harvested to determine the influence of disease control effectiveness on soybean crop yield. A total of one meter of the three center rows of each plot was manually harvested. The production was weighed on a precision scale to determine the weight of a thousand seeds (TSW) (g) and the crop yield (kg ha−1) of each treatment, considering humidity correction to 13% [36].

2.6. Data Analysis

Statistical analysis for the construction of the disease prediction model was performed according to the procedure that was already described, using the same script for data processing [18]. The prediction model was executed during the experiment to aid in the decision-making of the control and it was cross-validated (n = 10) with the original database that was used for the construction of this prediction model to confirm the disease severity classes identification [18].
The results were analyzed in the factorial scheme that was described and submitted to the analysis of homogeneity and normality. The data were submitted to analysis of variance (ANOVA) using the F test and, when significant, compared by Tukey’s test, both at 5% of significance.
Since a significant difference was found between the experimental field replications (p < 0.01), both areas were compared separately. In both fields, the spray deposit and coverage mean values were compared separately for each region of the canopy (upper and lower). The percentage of control was calculated based on the AUDPC of the control without application of each cultivar (T1 and T2). All statistical analyses and models were conducted using the R 3.6.3 software [37].

3. Results

3.1. Soybean Rust Detection and Application Timings

The periods of disease detection in the field trials according to each treatment are shown in Figure 2. In Field 1, SBR’s first symptoms were detected in both cultivars (Treatments 7 and 8) concomitant to the scheduled application (R1) at 50 DAE, and, therefore, T7 and T8 were sprayed on the same date as T3 and T4. The model was able to detect plants that were classified as “low severity” 6 days after that and in both cultivars when symptoms were about 0.2% severity. In the second experimental replication (Field 2), disease symptoms were detected in susceptible and partially resistant cultivars (T7 and T8) at 70 DAE, as well as the model was also able to detect diseased plants in the susceptible cultivar (T5), which had more characteristic symptoms of the disease at this time. For the partially resistant cultivars, which demonstrated significantly slower disease progression, the model detected diseased plants 7 days after that. The calendarized application (R1) was sprayed at the same time as in Field 1 (50 DAE).
At both Fields 1 and 2, a peak of LAI was observed around 60–65 DAE, occurring concurrently with the second application of the treatments (Figure 3). This period corresponds to the moment when higher foliar density is found, offering greater deposit and spray penetration challenges.

3.2. Spray Deposit

Consistent results were observed between the field trials and regarding the effect of the application timings on the qualitative and quantitate spraying variables. No significant effect was found in the interaction of the factors for any deposit values in both field trial replications.
In the first application of Field 1, a greater spray deposit was observed in the upper region of the crop canopy when spraying based on the prediction model (0.85 µL cm−2), compared to when based on symptoms (0.65 µL cm−2) and calendarized (0.55 µL cm−2) that were sprayed previously (Figure 4). In the lower region of the canopy, a significant difference was found only between the cultivars, where higher spray deposits (p < 0.01) were found in partially resistant cultivars (0.30 µL cm−2) than in the susceptible cultivars (0.09 µL cm−2) (Figure 5).
In the second application, the application timings affected the deposit values at both the upper and lower canopy regions, disregarding the cultivar that was used. Greater spray deposit was found when spraying based on the prediction model at the upper (0.98 µL cm−2) and lower (0.21 µL cm−2) canopy regions (Figure 6). The prediction model application timings (T5 and T6) occurred at a different moment than the other two treatments.
Furthermore, the cultivar also affected the quantity of the deposits that were found in the lower region of the canopy, in which the partially resistant cultivar (0.17 µL cm−2) significantly overcame the susceptible cultivar (0.05 µL cm−2) (Figure 7). Irrespective of the application timing, an expressive reduction in the spray deposit was found again at the lower region of the canopy, demonstrated by an uneven spray distribution across the canopy regions and reduced the capacity of droplets penetration into the lower regions.
The second experimental field trial presented a similar trend in deposit values. In the first application, a significantly greater deposit was found in the upper region in treatments that were sprayed at the first symptoms appearance (0.82 µL cm−2) and prediction model decision (0.84 µL cm−2) (Figure 8). These applications happened 20 days after the calendarized application (0.66 µL cm−2), once more demonstrating better spray deposit in the region with the later application. Moreover, a higher mean deposit value was found in the lower region of partially resistant cultivars (0.13 µL cm−2) than those found in the susceptible cultivars (0.05 µL cm−2) (Figure 9).
In the second application of Field 2, significant differences were found only in the lower region of the canopy. The application timings that were based on the prediction model and symptoms that were applied on the same date, obtained significantly higher deposition than the scheduled application (Figure 10). The higher mean deposit values were also found in the lower region of the canopy of partially resistant cultivars (Figure 11).

3.3. Spray Coverage

Finally, the qualitative evaluations that were based on the percentage of spray coverage were also consistent between the experimental trials. For Field 1, no significant difference was found in the first application, irrespective of application timing, cultivar, or canopy region. In the second application, a higher percentage of coverage was observed in the lower region of the canopy according to the cultivar, in which there was greater coverage in the partially resistant cultivars (42.3%) compared to the susceptible ones (31.8%) (Figure 12).
On the other hand, a significant difference was found in the interaction between the factors in the first application of Field 2 only in the lower region of the canopy (Table 3). A higher percentage was found at the lower-positioned targets of partially resistant cultivar treatments (21.9%) compared to the susceptible cultivars (5.14%) when spraying at a calendarized timing. Moreover, the coverage of the partially resistant cultivars were also significantly higher compared to the other application timings. In the second application, another difference was found in the lower region of the canopy as a function of application timings, in which the scheduled spraying timing presented significantly lower values compared to the others (Figure 13).

3.4. Effect of Application Timings on SBR Control and Crop Defoliation

The disease severity (AUDPC) and the control efficacy in both experimental fields are shown in Table 4. In general, a significant reduction in the disease severity was observed in all treatments with the fungicide application compared to the control. In addition, there was also a significant difference between the cultivars, in which the partially resistant cultivar presented lower AUDPC, regardless of the application timing. For the percentage of control, a difference was found only in Field 2, where a higher disease severity was observed in all the treatments and, therefore, with higher disease pressure. In this field trial, the percentage of control with the scheduled application (R1) in the partially resistant cultivar was significantly lower than in the others. Although the percentage of control seems much lower in the partially resistant cultivars, it is possibly due to the lower severity that was found even for the control treatment without applications.
The disease progress curves were considerably similar between both experimental fields (Figure 14). It is possible to observe greater development starting at 42 days after the first application (DAA) and rapid growth after this moment. On the other hand, greater severity progress was found at Field 2, where other treatments were also affected by the disease and promoted greater disease development. Overall, susceptible cultivars showed greater development, especially the control without fungicide application.
For the level of defoliation as an indirect severity evaluation, the spectroradiometer proved to be effective in representing the leaf stand level of the treatments (Figure 15). As the crop moved towards the end of the season, a clear reduction in the LAI was observed in both field trials regardless of the treatment.

3.5. Effect of SBR on Crop Yield

Regarding the effect of the disease on the crop yield (kg ha−1), no significant differences were found in the interaction of factors. However, the application timings affected the crop yield in Field 2, in which the control treatment without application presented a lower crop yield (p < 0.05) compared to the other application timings (Table 5). The numerical difference was kept similar in Field 1.
Likewise, a lower TWS was observed for the control treatment without control, significantly (p < 0.05) for Field 1 (Table 5). Furthermore, the effect of the cultivars was also observed for TWS in Field 2, in which the susceptible cultivar presented lower TWS (160.9 g) than the partially resistant cultivar (178.8 g) (Figure 16).

4. Discussion

This study evaluated the applicability of an SBR prediction model and the effect of different fungicide application timings on the spraying quality as well as on disease control and crop yield. Based on the disease epidemiology, it is important to provide proper SBR control as soon as possible, aiming to avoid the quick dispersal of spores and the emergence of epidemics in the field [2,38]. In this scenario, adequate disease monitoring is key to improving the application timing accuracy, spraying quality, and control efficacy. Besides, remote sensing and the usage of prediction models as support decision systems showed to be a valuable tool for IDM of SBR.
The prediction model that was applied here was able to identify plants with severity levels as low as 0.2% severity. The disease detection matched with the first appearance of symptoms, which can be considered an advantage since remote sensing may be applied in a simpler and faster way than extensive scouting in the field [15], especially with innovations in technology such as hyperspectral cameras and drone imaging of wide fields in lesser times [13]. The possibility of identifying the disease as quickly as possible can provide greater chances of a successful disease management program.
The application timings impacted the spraying quality, resulting in significant differences in the spray deposit and coverage at different parts of the crop canopy. A difference in the quantity of product that was deposited in the leaves, as well as the area that was covered by the spraying, may play a significant effect on the control of the biological agent, interfering with the control effectiveness and epidemic management. For instance, Berger-Neto et al. [39] reported higher white mold (Sclerotinia sclerotiorum) incidence in soybean with treatments that produced lower spray deposits, especially in lower canopy regions.
In our study, a marked difference was found mainly with later application timings. It promoted greater retention of spray deposit on the upper region of the canopy, therefore demonstrating a greater barrier for spray penetration. This outcome was also observed as the percentage of coverage. In Field 1, spraying that was based on the prediction model was conducted 6 days later than the others, which already promoted slightly greater retention in the upper section. None of the treatments were able to promote proper deposition and coverage in the lower region of the canopy. Meanwhile, the second application was conducted when a reduction of LAI was already started due to crop defoliation [23], influencing the penetration of the spray into the canopy and, therefore, promoting better deposition.
Nonetheless, the recommended timing of application at the end of the vegetative growth stage and reproductive growth stage is due to the possibility of still reaching the entire vertical profile of the plants and, therefore, achieving a better spray distribution, besides targeting the period of greater disease development [38]. However, it was seen that even spraying on R1 did not provide a uniform distribution. Furthermore, most fungicides that are recommended for SBR control act preventively and curatively, while most applications aim to act preventively to avoid any disease incidence and proliferation throughout the field [1]. Oppositely, when sprayed too far away from the first disease incidence in the field it is possible that the residual of the fungicide may not still be as active as needed, leading to lower efficacy [40].
Müller et al. [8] reported similar results of soybean rust severity at different application timings. Higher SBR severity (AUDPC) was observed when the application was the furthest from the first disease incidence. So, an application after R3 produced the best control when the disease incidence happened only after this period, whereas an application at R1 resulted in the best control efficacy when the disease was first seen before this growth stage and, therefore, closer to the disease incidence. Although an earlier application during the vegetative growth stage may provide better fungicide distribution in the canopy [24], the residual effect of the fungicide may not last until the period of higher disease incidence [40].
Moreover, Müller et al. [8] found that applications that were conducted closer to the first disease observations resulted in higher crop yields compared to other treatments that were sprayed after a longer time. This information also reassures the importance of proper disease monitoring to reduce both unnecessary fungicide applications in the field and to improve application timing accuracy and disease epidemic control. This is especially true when considering a country with proper environmental conditions for SBR development, along with fungicide application as the main source of control that is majorly characterized by several applications of only a few chemical choices (triazole, strobilurin, and carboxamides) that contributes to the selection of resistant populations [1,5,6].
In this study, especially for Field 2 where the disease started only 20 days after the calendarized application, at least one application would have been saved considering the whole crop season and the 15-day interval between applications. A smaller number of applications could reduce the chemical usage and, therefore, decrease fungal exposure to these fungicides, while a more accurate application having the potential to cease fungal development at determinant points that can increase control efficacy and decrease the chance of fungal survival in the field [8]. Finally, a better distribution throughout the crop canopy could also ensure a lower risk of over- or under-application rate, which could also lead to fungicide resistance [41]. Therefore, the benefits of sensor-based prediction models might come from the reduction of the number of applications or the fungicide efficacy that it might achieve.
The structural differences of the cultivars also played a major role in the spraying technology. The partially resistant cultivar (TMG7063) visually had less inter-row canopy closure compared to the susceptible cultivar (DS6217), in addition to a lower LAI. Therefore, the architecture format of the cultivar allowed a better penetration into the canopy, represented by a better spray distribution. Thereby, in most spraying that was conducted, better spray deposition and coverage were observed in the lower section of the plants.
A study unveiled the dynamic spray deposition according to different soybean cultivars’ architecture, reporting a significant increase in the spray coverage and overall penetration capacity when spraying in cultivars with lower height, LAI, and numbers of branches [22]. An increase in the coverage of almost 96% at the lower canopy regions was found for these cultivars.
The coverage and spray deposit results were very similar and consistent. The significant reduction in spraying reaching capacity into the innermost regions of the canopy stood out in all cases. Therefore, it highlights the need to adjust the technology according to the application timing, especially based on the LAI and canopy closure [22]. Alternatives such as using specialized spray nozzles or higher application volumes can be very responsive to provide significant improvements [9,39,42,43].
The results that were obtained from the application technology corroborate the disease severity that was found. The treatments with higher spray deposit and coverage also promoted slightly lesser disease severity, such as those that were observed when spraying based on the prediction model and first symptoms. Therefore, in general, the spraying decision-making that was based on remote sensing data improved or at least maintained the control effectiveness level of conventional methods. Furthermore, greater disease development was observed 42 days after the first application (DAA). Therefore, the earlier applications may have promoted a lower control rate due to greater distance from the peak of disease severity and lower fungicide residual.
With regards to the disease severity, the partially resistant cultivar was also considered to be highly effective in terms of suppressing the level of disease growth, with a marked progress curve reduction and lower severity indexes, even when without application. Along with the aforementioned cultivar structure that allowed a better spray distribution, all these features contribute to disease management and are, therefore, considered an excellent tool for the IDM.
To date, seven soybean resistance loci have been identified [44,45]. Vittal et al. [11] identified varied infection capacities of P. pachyrhizi among different soybean genotypes that were tested, reporting reduced hyphae development in resistant soybean cultivars. Resistant (immune reaction) and partially resistant cultivars are known to produce fewer pustules, lesions, and longer latent periods [10].
Another important factor is that, in this study, the disease incidence happened later than the usual first appearance, which can be visually identified in the curve progress. In this study, the later application timings at the disease onset provided better control. It is possible that in regions where there is higher disease pressure, monitoring by remote sensing can even anticipate the application compared to the calendarized applications. Improvements in disease monitoring to supply information to the decision support system have the potential to improve control efficacy as well as to even reduce the number of applications in regions of lower disease pressure [13,46,47,48].
Moreover, the use of remote sensing to assess the disease severity and defoliation was also considered successful, as reported by other authors [35,49]. Although no significant difference was found between the treatments, the LAI that was calculated through the integration of NDVI and the Beer–Lambert law [35] properly represented the defoliation levels, allowing better visualization of the effect of the disease on the crop leaf mass throughout the crop season.
The fungicide application promoted a reduced disease effect on the soybean crop yield, allowing greater yields irrespective of the cultivar or application timings. However, the differences that were found for each treatment control and spraying quality were not fully correlated to the crop yield, despite the proximity of the values that were found. For example, in both field repetitions, the calendarized application which obtained poorer spraying quality and higher disease severity, also achieved slightly higher crop yield than the others, even though it was only numerical. The crop yield is mostly affected when disease pressure is high, such as in the values that were observed in the control treatment without application. Besides, the effect of SBR on crop yield is most pronounced when its progression occurs at pod formation and filling [44,50]. Other injuries than visual symptoms may also play an important role in affecting the crop yield, such as on the carbon exchange rates [51].
Overall, the model was considered an effective tool and showed promising results to be used as a tool in the integrated management of SBR, since the detection periods were similar to those that were based on visual diagnosis and with potential for maintenance or improvement of disease control. In addition, improvements in assessment techniques are reported as a result of time savings compared to visual assessments. For cultivars with lower disease severity, or situations where there is a lower risk of epidemic or lower disease severity and incidence, it could be considered even more effective. It is important to state that the applicability of remote sensing techniques and data obtainment to supply these types of prediction models is difficult since many variables can interfere with the optical analyses, such as climate conditions and the presence of other injuries [13,19]. More studies are encouraged to be conducted to improve the prediction model database as well as to implement it in IDM programs.
Our results denote the practicality of using remote sensing and prediction models in the integrated management of the disease. Positive results were found, indicating the possibility of integrating it among the currently present management techniques. In addition, the timing of application impacted the application technology, resulting in significant differences in spray deposition and coverage, which have a high potential to also interfere with the effective control and management of epidemics. The application at lower LAI promoted better spray distribution and better SBR control efficacy. Finally, partially resistant cultivars also played a major role in the SBR control and as a powerful tool of the IDM.

Author Contributions

Conceptualization, M.M.N., R.M.N. and C.G.R.; Data curation, E.D.V. and C.G.R.; Formal analysis, M.M.N. and R.M.N.; Funding acquisition, E.D.V. and C.G.R.; Investigation, M.M.N., F.N.d.S., L.d.S.L., F.d.S.S.J. and B.R.d.F.; Methodology, F.N.d.S., L.d.S.L., F.d.S.S.J. and B.R.d.F.; Project administration, M.M.N. and C.G.R.; Resources, E.D.V. and C.G.R.; Supervision, C.G.R.; Validation, M.M.N., R.M.N. and E.D.V.; Visualization, F.N.d.S. and R.M.N.; Writing—original draft, M.M.N.; Writing—review & editing, M.M.N. and C.G.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the São Paulo Research Foundation (FAPESP), grant #2018/26486-0 and #2018/24869-0, by the National Council for Scientific and Technological Development, grant #142443/2018, and in part by Coordination for the Improvement of Higher Education Personnel (CAPES) (Grant number 001).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank Edivaldo Velini and Caio Carbonari for their assistance with equipment and technical knowledge on the methodology, Sergio A. Rodrigues for reviewing in the statistical analysis, and the School of Agriculture, Sao Paulo State University for assistance in trials seeding and conductance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Godoy, C.V.; Seixas, C.D.S.; Soares, R.M.; Marcelino-Guimarães, F.C.; Meyer, M.C.; Costamilan, L.M. Asian soybean rust in Brazil: Past, present, and future. Pesqui. Agropecuária Bras. 2016, 51, 407–421. [Google Scholar] [CrossRef]
  2. Langenbach, C.; Campe, R.; Beyer, S.F.; Mueller, A.N.; Conrath, U. Fighting Asian soybean rust. Front Plant Sci. 2016, 7, 797. [Google Scholar] [CrossRef] [PubMed]
  3. Juliatti, F.C.; Azevedo, L.A.S.; Cristina, J. Strategies of chemical protection for controlling soybean rust. In Soybean; Kasai, M., Ed.; IntechOpen: London, UK, 2017; Volume 1, pp. 35–62. [Google Scholar]
  4. Hossain, M.M.; Yamanaka, N. Pathogenic variation of Asian soybean rust pathogen in Bangladesh. J. Gen. Plant Pathol. 2019, 85, 90–100. [Google Scholar] [CrossRef]
  5. Müller, M.A.; Stammler, G.; De Mio, L.L.M. Multiple resistance to DMI, QoI and SDHI fungicides in field isolates of Phakopsora pachyrhizi. Crop Prot. 2021, 145, 105618. [Google Scholar] [CrossRef]
  6. Dalla Lana, F.; Paul, P.A.; Godoy, C.V.; Utiamada, C.M.; da Silva, L.H.C.; Siqueri, F.V.; Forcelini, C.A.; Jaccoud-Filho, D.D.S.; Miguel-Wruck, D.S.; Borges, E.P.; et al. Meta-analytic modeling of the decline in performance of fungicides for managing soybean rust after a decade of use in Brazil. Plant Dis. 2018, 102, 807–817. [Google Scholar] [CrossRef]
  7. Ozkan, H.E.; Zhu, H.; Derksen, R.C.; Guler, H.; Krause, C. Evaluation of various spraying equipment for effective application of fungicides to control Asian soybean rust. Asp. Appl. Biol. 2006, 77, 423–431. [Google Scholar]
  8. Mueller, T.A.; Miles, M.R.; Morel, W.; Marois, J.J.; Wright, D.L.; Kemerait, R.C.; Hartman, G.L. Effect of fungicide and timing of application on soybean rust severity and yield. Plant Dis. 2009, 93, 243–248. [Google Scholar] [CrossRef]
  9. Cunha, J.P.A.R.; Juliatti, F.C.; Reis, E.F. Tecnologia de aplicação de fungicida no controle da ferrugem asiática da soja: Resultados de oito anos de estudos em minas gerais e goiás. Biosci. J. 2014, 30, 950–957. [Google Scholar]
  10. Hartman, G.L.; Miles, M.R.; Frederick, R.D. Breeding for resistance to soybean rust. Plant Dis. 2005, 89, 664–666. [Google Scholar] [CrossRef]
  11. Vittal, R.; Paul, C.; Hill, C.B.; Hartman, G.L. Characterization and quantification of fungal colonization of Phakopsora pachyrhizi in soybean genotypes. Phytopathology 2014, 104, 86–94. [Google Scholar] [CrossRef]
  12. Twizeyimana, M.; Hartman, G.L. Effect of selected biopesticides in reducing soybean rust (Phakopsora pachyrhizi) development. Plant Dis. 2019, 103, 2460–2466. [Google Scholar] [CrossRef] [PubMed]
  13. Corley, C.D.; Pullum, L.L.; Hartley, D.M.; Benedum, C.; Noonan, C.; Rabinowitz, P.M.; Lancaster, M.J. Disease prediction models and operational readiness. PLoS ONE 2014, 9, e91989. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Martinelli, F.; Scalenghe, R.; Davino, S.; Panno, S.; Scuderi, G.; Ruisi, P.; Dandekar, A.M. Advanced methods of plant disease detection. A review. Agron. Sustain. Dev. 2015, 35, 1–25. [Google Scholar] [CrossRef]
  15. Zhang, J.; Huang, Y.; Pu, R.; Gonzalez-Moreno, P.; Yuan, L.; Wu, K.; Huang, W. Monitoring plant diseases and pests through remote sensing technology: A review. Comput. Electron. Agric. 2019, 165, 104943. [Google Scholar] [CrossRef]
  16. Cui, D.; Zhang, Q.; Li, M.; Zhao, Y.; Hartman, G.L. Detection of soybean rust using a multispectral image sensor. Sens. Instrum. Food Qual. Saf. 2009, 3, 49–56. [Google Scholar] [CrossRef]
  17. Furlanetto, R.H.; Nanni, M.R.; Mizuno, M.S.; Crusciol, L.G.T.; da Silva, C.R. Identification and classification of Asian soybean rust using leaf-based hyperspectral reflectance. Int. J. Remote Sens. 2021, 42, 4177–4198. [Google Scholar] [CrossRef]
  18. Negrisoli, M.M.; Negrisoli, R.M.; da Silva, F.N.; Lopes, L.L.; Souza Júnior, F.S.; Velini, E.D.; Carbonari, C.A.; Rodrigues, S.A.; Raetano, C.G. Soybean rust detection and disease severity classification by remote sensing. Agron. J. 2022. [Google Scholar] [CrossRef]
  19. Bajwa, S.G.; Rupe, J.C.; Mason, J. Soybean disease monitoring with leaf reflectance. Remote Sens. 2017, 9, 127. [Google Scholar] [CrossRef]
  20. Donatelli, M.; Magarey, R.D.; Bregaglio, S.; Willocquet, L.; Whish, J.P.; Savary, S. Modelling the impacts of pests and diseases on agricultural systems. Agric. Syst. 2017, 155, 213–224. [Google Scholar] [CrossRef]
  21. Twizeyimana, M.; Hartman, G.L. Sensitivity of Phakopsora pachyrhizi isolates to fungicides and reduction of fungal infection based on fungicide and timing of application. Plant Dis. 2017, 101, 121–128. [Google Scholar] [CrossRef]
  22. Müller, M.; Rakocevic, M.; Caverzan, A.; Boller, W.; Chavarria, G. Architectural characteristics and heliotropism may improve spray droplet deposition in the middle and low canopy layers in soybean. Crop Sci. 2018, 58, 2029–2041. [Google Scholar] [CrossRef]
  23. Tagliapietra, E.L.; Streck, N.A.; Da Rocha, T.S.M.; Richter, G.L.; Da Silva, M.R.; Cera, J.C.; Guedes, J.V.C.; Zanon, A.J. Optimum leaf area index to reach soybean yield potential in subtropical environment. Agron. J. 2018, 10, 932–938. [Google Scholar] [CrossRef] [Green Version]
  24. Negrisoli, M.M.; Raetano, C.G.; de Souza, D.M.; Souza, F.M.S.; Bernardes, L.M.; Junior, L.D.B.; Rodrigues, D.M.; Sartori, M.M.P. Performance of new flat fan nozzle design in spray deposition, penetration and control of soybean rust. Eur. J. Plant Pathol. 2019, 155, 755–767. [Google Scholar] [CrossRef]
  25. Brevant Sementes—Cultivares de Soja: DS6217IPRO. Available online: https://www.brevant.com.br/produtos/soja/ds6217ipro.html (accessed on 26 January 2022).
  26. Tropical Melhoramento Genético (TMG)—Cultivares de Soja: TMG 7063 IPRO. Available online: http://www.tmg.agr.br/ptbr/cultivar/tmg-7063-ipro (accessed on 26 January 2022).
  27. Abdel-Rahman, E.M.; Mutanga, O.; Odindi, J.; Adam, E.; Odindo, A.; Ismail, R. A comparison of partial least squares (PLS) and sparse PLS regressions for predicting yield of Swiss chard grown under different irrigation water sources using hyperspectral data. Comput. Electron. Agric. 2014, 106, 11–19. [Google Scholar] [CrossRef]
  28. Savitzky, A.; Golay, M.J. Smoothing and differentiation of data by simplified least squares procedures. Analytical. Chem. 1964, 36, 1627–1639. [Google Scholar] [CrossRef]
  29. Bohnenkamp, D.; Behmann, J.; Mahlein, A.K. In-field detection of yellow rust in wheat on the ground canopy and UAV scale. Remote Sens. 2019, 11, 2495. [Google Scholar] [CrossRef]
  30. Abdi, H.; Williams, L.J.; Valentin, D. Multiple factor analysis: Principal component analysis for multitable and multiblock data sets. Wiley Interdiscip. Rev. Comput. Stat. 2013, 5, 149–179. [Google Scholar] [CrossRef]
  31. Palladini, L.A.; Raetano, C.G.; Velini, E.D. Choice of tracers for the evaluation of spray deposits. Sci. Agric. 2005, 62, 440–445. [Google Scholar] [CrossRef]
  32. Jonckheere, I.; Fleck, S.; Nackaerts, K.; Muys, B.; Coppin, P.; Weiss, M.; Baret, F. Review of methods for in situ leaf area index determination: Part I. Theories, sensors and hemispherical photography. Agric. For. Meteorol. 2004, 121, 19–35. [Google Scholar] [CrossRef]
  33. Godoy, C.V.; Koga, L.J.; Canteri, M.G. Diagrammatic scale for assessment of soybean rust severity. Fitopatol. Bras. 2006, 31, 63–68. [Google Scholar] [CrossRef]
  34. Campbell, C.L.; Madden, L.V. Introduction to Plant Disease Epidemiology; John Wiley & Sons: New York, NY, USA, 1994; 532p. [Google Scholar]
  35. Tan, C.W.; Zhang, P.P.; Zhou, X.X.; Wang, Z.X.; Xu, Z.Q.; Mao, W.; Yun, F. Quantitative monitoring of leaf area index in wheat of different plant types by integrating NDVI and Beer-Lambert law. Sci. Rep. 2020, 10, 1–10. [Google Scholar] [CrossRef] [PubMed]
  36. Brasil. Ministério da Agricultura, Pecuária e Abastecimento. Regras Para Análise de Sementes; Ministério da Agricultura, Pecuária e Abastecimento—Secretaria de Defesa Agropecuária: Brasília, Brazil, 2009; 395p.
  37. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2013. [Google Scholar]
  38. Kelly, H.Y.; Dufault, N.S.; Walker, D.R.; Isard, S.A.; Schneider, R.W.; Giesler, L.J.; Hartman, G.L. From select agent to an established pathogen: The response to Phakopsora pachyrhizi (soybean rust) in North America. Phytopathology 2015, 105, 905–916. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Berger-Neto, A.; Jaccoud-Filho, D.S.; Wutzki, C.R.; Tullio, H.E.; Pierre, M.L.C.; Manfron, F.; Justino, A. Effect of spray droplet size, spray volume and fungicide on the control of white mold in soybeans. Crop Prot. 2017, 92, 190–197. [Google Scholar] [CrossRef]
  40. Balardin, R.; Madalosso, M.; Debortoli, M.; Lenz, G. Factors affecting fungicide efficacy in the tropics. In Fungicides; Carisse, O., Ed.; InTech Open: London, UK, 2010; pp. 23–38. [Google Scholar]
  41. Van Den Bosch, F.; Paveley, N.; Shaw, M.; Hobbelen, P.; Oliver, R. The dose rate debate: Does the risk of fungicide resistance increase or decrease with dose? Plant Pathol. 2011, 60, 597–606. [Google Scholar] [CrossRef]
  42. Ferguson, J.C.; Hewitt, A.J.; O’Donnell, C.C. Pressure, droplet size classification, and nozzle arrangement effects on coverage and droplet number density using air-inclusion dual fan nozzles for pesticide applications. Crop Prot. 2016, 89, 231–238. [Google Scholar] [CrossRef]
  43. Sharpe, S.M.; Boyd, N.S.; Dittmar, P.J.; Macdonald, G.E.; Darnell, R.L.; Ferrel, J.A. Spray penetration into a strawberry canopy as affected by canopy structure, nozzle type, and application volume. Weed Technol. 2017, 32, 80–84. [Google Scholar] [CrossRef]
  44. Childs, S.P.; Buck, J.W.; Li, Z. Breeding soybeans with resistance to soybean rust (Phakopsora pachyrhizi). Plant Breed. 2018, 137, 250–261. [Google Scholar] [CrossRef]
  45. Hartman, G.L.; Hill, C.B.; Twizeyimana, M.; Miles, M.R.; Bandyopadhyay, R. Interaction of soybean and Phakopsora pachyrhizi, the cause of soybean rust. CAB Rev. Perspect. Agric. Vet. Sci. Nutr. Nat. Resour. 2011, 6, 59. [Google Scholar] [CrossRef]
  46. Verreet, J.A.; Klink, H.; Hoffmann, G.M. Regional monitoring for disease prediction and optimization of plant protection measuares: The IPM wheat model. Plant Dis. 2000, 84, 816–826. [Google Scholar] [CrossRef]
  47. Sapak, Z.; Salam, M.U.; Minchinton, E.J.; MacManus, G.P.V.; Joyce, D.C.; Galea, V.J. POMICS: A simulation disease model for timing fungicide applications in management of powdery mildew of cucurbits. Phytopathology 2017, 107, 1022–1031. [Google Scholar] [CrossRef]
  48. Gama, A.B.; Peres, N.A.; Singerman, A.; Dewdney, M.M. Evaluation of disease alert systems for postbloom fruit drop of citrus in Florida and economic impact of adopting the Citrus Advisory System. Crop Prot. 2022, 155, 105906. [Google Scholar] [CrossRef]
  49. Yuan, H.; Yang, G.; Li, C.; Wang, Y.; Liu, J.; Yu, H.; Yang, X. Retrieving soybean leaf area index from unmanned aerial vehicle hyperspectral remote sensing: Analysis of RF, ANN, and SVM regression models. Remote Sens. 2017, 9, 309. [Google Scholar] [CrossRef] [Green Version]
  50. Kawuki, R.S.; Tukamuhabwa, P.; Adipala, E. Soybean rust severity, rate of rust development, and tolerance as influenced by maturity period and season. Crop Prot. 2004, 23, 447–455. [Google Scholar] [CrossRef]
  51. Kumudini, S.; Prior, E.; Omielan, J.; Tollenaar, M. Impact of Phakopsora pachyrhizi infection on soybean leaf photosynthesis and radiation absorption. Crop Sci. 2008, 48, 2343–2350. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Study area map of both field trial replications.
Figure 1. Study area map of both field trial replications.
Agronomy 12 02119 g001
Figure 2. Graphical representation of fungicide application timings along with qualitative and quantitative spraying evaluations as a function of disease detection that was proposed for each treatment.
Figure 2. Graphical representation of fungicide application timings along with qualitative and quantitative spraying evaluations as a function of disease detection that was proposed for each treatment.
Agronomy 12 02119 g002
Figure 3. Leaf area index (LAI) of Field Trials 1 (A) and 2 (B) that were analyzed during the period interval of the fungicide applications.
Figure 3. Leaf area index (LAI) of Field Trials 1 (A) and 2 (B) that were analyzed during the period interval of the fungicide applications.
Agronomy 12 02119 g003
Figure 4. Mean deposit values (µL cm−2) that were collected in the upper region of the soybean canopy at different application timings in the first application of Field 1. The means followed by the same letter did not differ by the Tukey test at 5% probability (p < 0.05).
Figure 4. Mean deposit values (µL cm−2) that were collected in the upper region of the soybean canopy at different application timings in the first application of Field 1. The means followed by the same letter did not differ by the Tukey test at 5% probability (p < 0.05).
Agronomy 12 02119 g004
Figure 5. Mean deposit values (µL cm−2) that were collected in the lower region of the soybean canopy of different soybean cultivars in the first application of Field 1. The means followed by the same letter did not differ by the Tukey test at 5% probability (p < 0.05).
Figure 5. Mean deposit values (µL cm−2) that were collected in the lower region of the soybean canopy of different soybean cultivars in the first application of Field 1. The means followed by the same letter did not differ by the Tukey test at 5% probability (p < 0.05).
Agronomy 12 02119 g005
Figure 6. Mean deposit values (µL cm−2) that were collected in the upper (A) and lower (B) regions of the soybean canopy at different application timings in the second application of Field 1. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05) in each comparison.
Figure 6. Mean deposit values (µL cm−2) that were collected in the upper (A) and lower (B) regions of the soybean canopy at different application timings in the second application of Field 1. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05) in each comparison.
Agronomy 12 02119 g006
Figure 7. Mean deposit values (µL cm−2) that were collected in the lower region of the soybean canopy of different soybean cultivars in the second application of Field 1. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Figure 7. Mean deposit values (µL cm−2) that were collected in the lower region of the soybean canopy of different soybean cultivars in the second application of Field 1. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Agronomy 12 02119 g007
Figure 8. Mean deposit values (µL cm−2) that were collected in the upper region of the soybean canopy at different application timings in the first application of Field 2. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Figure 8. Mean deposit values (µL cm−2) that were collected in the upper region of the soybean canopy at different application timings in the first application of Field 2. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Agronomy 12 02119 g008
Figure 9. Mean deposit values (µL cm−2) that were collected in the lower region of the soybean canopy of different soybean cultivars in the first application of Field 2. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Figure 9. Mean deposit values (µL cm−2) that were collected in the lower region of the soybean canopy of different soybean cultivars in the first application of Field 2. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Agronomy 12 02119 g009
Figure 10. Mean deposit values (µL cm−2) that were collected in the upper region of the soybean canopy at different application timings in the second application of Field 2. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Figure 10. Mean deposit values (µL cm−2) that were collected in the upper region of the soybean canopy at different application timings in the second application of Field 2. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Agronomy 12 02119 g010
Figure 11. Mean deposit values (µL cm−2) that were collected in the lower region of the soybean canopy of different soybean cultivars in the second application of Field 2. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Figure 11. Mean deposit values (µL cm−2) that were collected in the lower region of the soybean canopy of different soybean cultivars in the second application of Field 2. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Agronomy 12 02119 g011
Figure 12. Mean percentage of coverage (%) in water-sensitive papers that were located at the lower region of the soybean canopy, according to different soybean cultivars in the second application of Field 1. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Figure 12. Mean percentage of coverage (%) in water-sensitive papers that were located at the lower region of the soybean canopy, according to different soybean cultivars in the second application of Field 1. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Agronomy 12 02119 g012
Figure 13. Mean percentage of coverage (%) in water-sensitive papers that were located at the lower region of the soybean canopy, according to different application timings in the second application of Field 2. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Figure 13. Mean percentage of coverage (%) in water-sensitive papers that were located at the lower region of the soybean canopy, according to different application timings in the second application of Field 2. The means followed by the same letter did not differ by the Tukey’s test at 5% probability (p < 0.05).
Agronomy 12 02119 g013
Figure 14. Soybean rust disease progress curve at Field 1 and Field 2 trial repetitions, according to the different soybean cultivars and application timings (treatments). Note: Susc_: susceptible soybean cultivar (DS6217); Resist_: partially resistant soybean cultivar (TMG7063). Control: without fungicide application; Calend: calendarized application (reproductive growth stage R1); Model: application timing based on the prediction model; Sympt: application timing based on the appearance of the first symptoms.
Figure 14. Soybean rust disease progress curve at Field 1 and Field 2 trial repetitions, according to the different soybean cultivars and application timings (treatments). Note: Susc_: susceptible soybean cultivar (DS6217); Resist_: partially resistant soybean cultivar (TMG7063). Control: without fungicide application; Calend: calendarized application (reproductive growth stage R1); Model: application timing based on the prediction model; Sympt: application timing based on the appearance of the first symptoms.
Agronomy 12 02119 g014
Figure 15. Defoliation assessment based on the leaf area index (LAI) through the integration of NDVI and Beer–Lambert law, according to the spectral curves of the susceptible and partially resistant soybean cultivars under soybean rust effect, across experimental evaluation periods of Field 1 (A) and Field 2 (B). Notes: Susc_: susceptible soybean cultivar (DS6217); Resist_: partially resistant soybean cultivar (TMG7063). Control: without fungicide application; Calend: calendarized application (reproductive growth stage R1); Model: application timing based on the prediction model; Sympt: application timing based on the appearance of the first symptoms. LAI calculated according to Tan et al. [35].
Figure 15. Defoliation assessment based on the leaf area index (LAI) through the integration of NDVI and Beer–Lambert law, according to the spectral curves of the susceptible and partially resistant soybean cultivars under soybean rust effect, across experimental evaluation periods of Field 1 (A) and Field 2 (B). Notes: Susc_: susceptible soybean cultivar (DS6217); Resist_: partially resistant soybean cultivar (TMG7063). Control: without fungicide application; Calend: calendarized application (reproductive growth stage R1); Model: application timing based on the prediction model; Sympt: application timing based on the appearance of the first symptoms. LAI calculated according to Tan et al. [35].
Agronomy 12 02119 g015
Figure 16. Mean values of thousand seeds weight (TSW) (g) as a function of soybean rust effect over different soybean cultivars in the experimental Field 2. The means followed by the same letter did not differ by Tukey’s test at 5% probability (p < 0.05).
Figure 16. Mean values of thousand seeds weight (TSW) (g) as a function of soybean rust effect over different soybean cultivars in the experimental Field 2. The means followed by the same letter did not differ by Tukey’s test at 5% probability (p < 0.05).
Agronomy 12 02119 g016
Table 1. Description of the treatments according to the soybean cultivar and application timing.
Table 1. Description of the treatments according to the soybean cultivar and application timing.
TreatmentsCultivarFirst Application TimingSecond Application Timing
1DS6217 IPRO (Susceptible)--
2TMG 7063 (Partially resistant)--
3DS6217 IPRO (Susceptible)Calend *Calend + 15
4TMG 7063 (Partially resistant)CalendCalend + 15
5DS6217 IPRO (Susceptible)ModelModel + 15
6TMG 7063 (Partially resistant)ModelModel + 15
7DS6217 IPRO (Susceptible)SymptSympt + 15
8TMG 7063 (Partially resistant)SymptSympt + 15
* Calend: calendarized application at reproductive stage R1; Model: spraying according to SBR prediction model; Sympt: spraying at the first symptoms observations; “+15”: second application conducted at 15 days after the first application; “-”: no application (control treatments).
Table 2. Description of the cultivars that were adopted in the experimental trials.
Table 2. Description of the cultivars that were adopted in the experimental trials.
DescriptionBrevant DS6217 IPRO 1TMG 7063 2
P. pachyrhizi susceptibilitySusceptiblePartially resistant (Inox)
Maturity groups6.26.3
Growth habitIndeterminateIndeterminate
TraitsIntacta RR2 PRO®Intacta RR2 PRO®
Data from Brevant 1 [25] and Tropical Melhoramento Genético 2 [26].
Table 3. Mean percentage coverage (%) on water-sensitive papers in the first application of Field 2, according to different application timings and soybean cultivars in each canopy region.
Table 3. Mean percentage coverage (%) on water-sensitive papers in the first application of Field 2, according to different application timings and soybean cultivars in each canopy region.
Application TimingCultivar
SusceptiblePartially Resistant
Calendarized5.14 aB21.97 aA
Prediction model6.57 aA2.06 bA
Symptoms6.74 aA4.27 bA
Cause of VariationFP
Application timing (AT)2.8850.087 NS
Cultivar ©0.9220.352 NS
AT x C3.9580.042 *
NS: Not significant; * significant a p ≤ 0.05 by F test. Means followed by the same letter did not differ according to Tukey’s test at 5% probability (p < 0.05). Lowercase letters compare between the means of application timings at each cultivar level (lines). Uppercase letters compare between the means of cultivars at each application timing (columns). Each statistical comparison was conducted separately for each canopy region (upper and lower).
Table 4. Mean values of disease severity (AUDPC) and control (%) according to the application timings and susceptible and partially resistant soybean cultivars, for each field trial replication.
Table 4. Mean values of disease severity (AUDPC) and control (%) according to the application timings and susceptible and partially resistant soybean cultivars, for each field trial replication.
FieldApplication TimingAUDPCControl (%)
CultivarCultivar
SusceptibleP. ResistantSusceptibleP. Resistant
Field 1Control200.80 aA40.50 aB--
Calendarized68.90 bA18.10 aB64.4051.50
Prediction model54.60 bA15.90 aB71.6059.20
Symptoms63.00 bA15.00 aB67.9058.90
Causes of VariationFPFP
App. Timing (AT)53.450<0.001 ***0.5760.574 NS
Cultivar ©177.450<0.001 ***3.7720.071 NS
AT x C26.450<0.001 ***0.0420.959 NS
Field 2Control519.20 aA65.70 aB--
Calendarized138.90 bA32.50 aB72.50 aA47.80 bB
Prediction model142.50 bA22.40 aB71.70 aA65.70 aA
Symptoms191.00 bA27.10 aB62.10 aA56.74 abA
Causes of VariationFPFP
App. timing (AT)43.076<0.001 ***143.823<0.001 ***
Cultiv©(C)188.949<0.001 ***8.5770.008 *
AT x C28.402<0.001 ***5.0190.009 *
NS: Not significant; * significant a p ≤ 0.05; *** significant a p ≤ 0.01 by F test. Means followed by the same letter did not differ according to Tukey’s test at 5% probability (p < 0.05). Lowercase letters compare between the means of application timings at each cultivar level (lines). Uppercase letters compare between the means of cultivars at each application timing (columns). Each statistical comparison was conducted separately for each canopy region (upper and lower).
Table 5. Mean values of soybean crop yield (kg ha−1) and thousand seeds weight (TSW) (g), according to soybean rust effect of different application timings and soybeans cultivars, for each field trial repetition.
Table 5. Mean values of soybean crop yield (kg ha−1) and thousand seeds weight (TSW) (g), according to soybean rust effect of different application timings and soybeans cultivars, for each field trial repetition.
Application TimingCrop YieldTSW
kg ha−1g
Field 1Field 2Field 1Field2
Control2393.7142799.833 b2393.712799.83 b
Calendarized3143.0543462.89 a3143.053462.89 a
Prediction model2738.0613055.73 ab2738.063055.73 ab
Symptoms2817.0343169.037 ab2817.033169.04 ab
F value2.713 NS4.092 *5.510 ***1.804 NS
CV (%)19.0512.316.478.54
NS: Not significant; *** significant at p ≤ 0.01; * significant at p ≤ 0.05 by F test. Means followed by the same letter did not differ according to Tukey’s test at 5% probability (p < 0.05).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Negrisoli, M.M.; Silva, F.N.d.; Negrisoli, R.M.; Lopes, L.d.S.; Souza Júnior, F.d.S.; Freitas, B.R.d.; Velini, E.D.; Raetano, C.G. Impact of Fungicide Application Timing Based on Soybean Rust Prediction Model on Application Technology and Disease Control. Agronomy 2022, 12, 2119. https://doi.org/10.3390/agronomy12092119

AMA Style

Negrisoli MM, Silva FNd, Negrisoli RM, Lopes LdS, Souza Júnior FdS, Freitas BRd, Velini ED, Raetano CG. Impact of Fungicide Application Timing Based on Soybean Rust Prediction Model on Application Technology and Disease Control. Agronomy. 2022; 12(9):2119. https://doi.org/10.3390/agronomy12092119

Chicago/Turabian Style

Negrisoli, Matheus Mereb, Flávio Nunes da Silva, Raphael Mereb Negrisoli, Lucas da Silva Lopes, Francisco de Sales Souza Júnior, Bianca Rezende de Freitas, Edivaldo Domingues Velini, and Carlos Gilberto Raetano. 2022. "Impact of Fungicide Application Timing Based on Soybean Rust Prediction Model on Application Technology and Disease Control" Agronomy 12, no. 9: 2119. https://doi.org/10.3390/agronomy12092119

APA Style

Negrisoli, M. M., Silva, F. N. d., Negrisoli, R. M., Lopes, L. d. S., Souza Júnior, F. d. S., Freitas, B. R. d., Velini, E. D., & Raetano, C. G. (2022). Impact of Fungicide Application Timing Based on Soybean Rust Prediction Model on Application Technology and Disease Control. Agronomy, 12(9), 2119. https://doi.org/10.3390/agronomy12092119

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop