Next Article in Journal
Variability in Water Capacity of Small-Leaved Linden Associated with Both the Presence of Honeydew and Various Sources of Pollution
Next Article in Special Issue
Further Molecular Diagnosis Determines Lack of Evidence for Real Seed Transmission of Tomato Leaf Curl New Delhi Virus in Cucurbits
Previous Article in Journal
Rhizoactinobacteria Enhance Growth and Antioxidant Activity in Thai Jasmine Rice (Oryza sativa) KDML105 Seedlings under Salt Stress
Previous Article in Special Issue
Yellow Leaf Disease Resistance and Melanaphis sacchari Preference in Commercial Sugarcane Cultivars
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Impact of Agricultural Mitigation Measures on the Spread of Sharka Disease in Sweet Cherry Orchards

by
Juan Pablo Gutiérrez-Jara
1,*,†,
Katia Vogt-Geisse
2,*,†,
Margarita C. G. Correa
1,†,
Karina Vilches-Ponce
3,
Laura M. Pérez
4 and
Gerardo Chowell
5
1
Centro de Investigación de Estudios Avanzados del Maule (CIEAM), Vicerrectoría de Investigación y Postgrado, Universidad Católica del Maule, Talca 3480112, Chile
2
Facultad de Ingeniería y Ciencias, Unidad Adolfo Ibáñez, Diagonal Las Torres 2640, Peñalolén, Santiago 7941169, Chile
3
Facultad de Ciencias Básicas, Universidad Católica del Maule, Avenida San Miguel 3605, Talca 3480112, Chile
4
Departamento de Física, Universidad de Tarapacá, Casilla 7D, Arica 1000000, Chile
5
School of Public Health, Georgia State University, Atlanta, GA 30303, USA
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Plants 2023, 12(19), 3442; https://doi.org/10.3390/plants12193442
Submission received: 6 September 2023 / Revised: 27 September 2023 / Accepted: 27 September 2023 / Published: 29 September 2023

Abstract

:
Sharka is a disease affecting stone fruit trees. It is caused by the Plum pox virus (PPV), with Myzus persicae being one of the most efficient aphid species in transmitting it within and among Prunus orchards. Other agricultural management strategies are also responsible for the spread of disease among trees, such as grafting and pruning. We present a mathematical model of impulsive differential equations to represent the dynamics of Sharka disease in the tree and vector population. We consider three transmission routes: grafting, pruning, and through aphid vectors. Grafting, pruning, and vector control occur as pulses at specific instants. Within the model, human risk perception towards disease influences these agricultural management strategies. Model results show that grafting with infected biological material has a significant impact on the spread of the disease. In addition, detecting infectious symptomatic and asymptomatic trees in the short term is critical to reduce disease spread. Furthermore, vector control to prevent aphid movement between trees is crucial for disease mitigation, as well as implementing awareness campaigns for Sharka disease in agricultural communities that provide a long-term impact on responsible pruning, grafting, and vector control.

1. Introduction

Sharka is one of the most severe diseases that affects stone fruit trees worldwide, caused by the Plum pox virus (PPV) of the genus Potyvirus from the family Potyviridae. Sharka symptoms range from slight leaf discoloration to tree death [1]. Although sensitivity to PPV can vary among Prunus cultivars, most of them show Sharka symptoms and a certain decrease in productivity [2]. The symptoms shown by infected plants can vary according to the virus strain, the spatial location, the Prunus species, the plant tissue, and the season, with spring being the best season to observe symptoms due to temperatures below 30 degrees Celsius [3,4]. Symptoms include other abnormal shapes, colors, or forms of leaves and fruits, premature drop of fruits, fruit lesions, etc. Asymptomatic infection may be more prevalent in warmer climates [5]. In general, infected plants can take several months to show symptoms and their appearance may be temporary and fade during the summer months [3,6].
Disease spread depends on biotic and abiotic vectors. Aphids are the leading biotic vectors of PPV from infested trees to healthy ones present in the orchard or surrounding [3]. More than 20 species of aphids can transmit the virus in a non-persistent manner, of which the species Myzus persicae is one of the most efficient in doing so [4]. The intensity of spread depends on several factors, including the efficacy of the aphid species in transmitting and its abundance of winged forms, the winds in the area, and the presence of susceptible cultivars around [5]. The virus can also be transmitted through vegetative propagation techniques in nursery trading such as grafting [3,4,5,7].
Unfortunately, infected trees cannot recover and even if the virus does not kill the plant, they remain PPV-positive. Previous research has reported that virus titers in PPV-positive trees can drop below the PPV detection levels, using ELISA, during warm weather [8]. Hence, infected trees serve as a reservoir for the virus if not eradicated, as well as nearby infected weeds and cultivars, due to the action of virus vectors. Consequently, PPV control relies on public and private prevention measures, such as laws and norms that regulate the trading of vegetative material susceptible to PPV and control illegal exchange trafficking; strict quarantine regulations for nursery stock; chemical vector control; proper eradication of infected material; and the control of weeds and susceptible vegetation [4,5].
Sharka disease was detected in Chile in 1992 and is now distributed in the main Prunus productive regions: Valparaíso, Metropolitana and O’Higgins; it was also recently detected in the Maule region [4]. Chile is the leading sweet cherry exporter worldwide and sweet cherries represent one of the fruits with more economic importance in the country; indeed, the 2022–2023 export season was expected to report nearly 3 billion USD to the country [9].
Mathematical modeling of disease transmission in plants has been developed to improve the understanding, forecasting, and warning, on different scales, of how infections spread in plants and plantations due to viruses, bacteria, and fungi [10]. In particular, epidemiological models of differential equations of the type SIR (susceptible-infectious-removed) and SEIR (susceptible-exposed-infectious-removed) were introduced to model the dynamics of vector-transmitted diseases for different types of plants; for example, citrus [11], banana [12], rice [13], and cereals [14]. Additionally, control techniques are analyzed in such SIR and SEIR models to evaluate the effects of roguing (removal–destruction of diseased plants) and replanting (healthy plants). The main results provide feasible control mechanisms for disease eradication, which could be achieved by implementing realistic roguing and replanting intensities [15,16,17].
In order to include variability in the transmission rates, SIR and SEIR models have included population dynamics and epidemiological aspects of vectors [18]. The interaction of two or more viruses and their coexistence in a host plant system has been modeled [14,18]. Seasonality has been represented by assuming periodic immigration and death rates, and biological control has been modeled by introducing an ecological predator into vector dynamics [19]. Stochasticity and delayed differential equations have been used in models to include uncertainty and time lags between vectors and disease dynamics [20].
The importance of using mathematical modeling to manage Sharka disease has already been highlighted [21]. In this work, we incorporate the impact of agricultural mitigation measures on the spread of Sharka disease in sweet cherry orchards and suggest strategies to prevent agricultural losses.
The proposed model explicitly incorporates risk perception of the farming community towards plant viruses that impacts the effective implementation of agricultural management strategies, which is more likely if individuals in the agricultural community are aware of the existence of a disease and their role in controlling it [22]. Campaigns that can help increase this awareness are a combination of agricultural extension programs that transfer consistent messages, communicate through farmers’ social referents and mass media, and implement group learning and individual communication with farm advisors, as well as new forms of knowledge exchange such as crowd-sourcing, since individuals may not even be aware of plant viruses [22,23]. In general, exposure to agricultural extension services that create awareness and the level of farmers’ education may be crucial for knowledgeable pest management [24].
The following section presents the mathematical model that expresses the dynamics of the Sharka virus, including variables of interest for the study of mitigation strategies. Subsequently, in Section 3, the results of the model are shown. Finally, in Section 4, the discussion and respective conclusions are opened.

2. Materials and Methods

Systems of differential equations have been a helpful tool for understanding the evolution of different phenomena. There are dynamics where the evolving process is subject to abrupt changes, represented in pulses. These perturbations have been called impulsive differential equations (IDE) and are used to model ecological, epidemiological, and other problems [25,26,27].
We present a mathematical model of impulsive differential equations to represent the dynamics of Sharka disease, including tree, vector, and human behavior. We denote the total tree population by N T and divide it according to its epidemiological states: susceptible ( S T ), exposed ( E T ), symptomatic infectious ( I T ), asymptomatic infectious ( A T ), and removed ( R T ). This last compartment refers to the trees that have been uprooted. Furthermore, we have divided the vector population ( N v ) into two states: susceptible ( S v ) and infectious ( I v ). Table 1 summarizes the definition of model variables.
The flow diagram in Figure 1 shows the dynamics of the disease in the tree population, where the segmented and dotted lines represent the instants (pulses) in which the continuous dynamics of the disease are altered. A susceptible tree ( S T ) can be infected in three ways: (i) through interaction with an infectious vector ( I v ) at a force of transmission
Ψ T = β T v , q r 2 S T I v N T ,
where β T v , q / r 2 represents the transmission rate, whose form can be interpreted as a dependence on the movement of the vector ( 1 / r 2 ) and a seasonal factor ( q = 0 , 1 , 2 , 3 ; for autumn, winter, spring, summer, respectively); (ii) through infected pruning equipment at certain instants (pulses, segmented line) at a force of infection
Φ T = β T p , q P P S T A T + I T N T ,
where the transmission rate β T p , q P / P , depends on a seasonal factor (q) and a dynamic variable P ( t ) , representing the agricultural worker’s risk perception of Sharka disease governed by the differential equation adopted from [28,29,30]
P = λ 1 ( P P ) + λ 2 ( I T + R T ) / N T
where λ 1 is the rate of resistance to change, P is the natural risk perception of the agricultural population and λ 2 is the speed of people’s reactions to the presence of symptomatic trees ( I T ) and removed trees ( R T ); (iii) through the introduction of infected grafts at certain instants (pulses, dotted line) at a rate λ e Δ T , where Δ T is the rate of introduction of grafts and λ e is the proportion of infected grafts. The average time that a tree remains exposed is 1 / α . The transition rate from an infected tree with symptoms to becoming asymptomatic is δ a , i and the reverse process occurs at a rate δ i , a . The average time it takes for a symptomatic and asymptomatic tree to be uprooted is given by 1 / γ I and 1 / γ A , respectively. 1 / Λ T is the average time it takes for a tree to be replaced by a susceptible (healthy) tree after being uprooted due to Sharka disease. Finally, the natural mortality rate of the tree is given by d T , and the recruitment (plantation) by Λ .
Figure 2 shows the dynamics of the model for the disease in the vector population. The force of infection of a vector when contacting an infectious tree is given by
Ψ v = β v q r 2 S v ( A T + I T ) / N T ,
where β v q / r 2 represents the rate of transmission, whose form can be interpreted as a dependence on the movement of the vector ( 1 / r 2 ) and a seasonal factor ( q = 0 , 1 , 2 , 3 ; for autumn, winter, spring, summer, respectively). An infected vector recovers at a rate γ v ; thus, 1 / γ v is the average time it takes for a vector to remain in its infectious state. Λ v and d v correspond to the recruitment and natural mortality rates, respectively. Finally, Θ q is an additional vector mortality rate, which occurs in certain instants (pulses) due to the implementation of preventive and vector control measures.
The model of impulsive differential equations given by system (5) follows the dynamics from Figure 1 and Figure 2. t = t i and t = t j correspond to the instants in which grafts are introduced and stone tree pruning is performed, respectively. t = t k are the moments of vector control, either by natural or chemical means. t = t n are the instants where awareness campaigns impact risk perception (P, see Equation (3)) of the agricultural community towards Sharka disease.
S T ˙ ( t ) = Λ ( β T v , q / r 2 ) S T I v / N T + Λ T R T d T S T E T ˙ ( t ) = ( β T v , q / r 2 ) S T I v / N T ( α + d T ) E T I T ˙ ( t ) = α E T + δ i , a A T ( δ a , i + γ I + d T ) I T A T ˙ ( t ) = δ a , i I T ( δ i , a + γ A + d T ) A T R T ˙ ( t ) = γ I I T + γ A A T Λ T R T S v ˙ ( t ) = Λ v q ( β v q / r 2 ) S v ( A T + I T ) / N T + γ v I v d v S v I v ˙ ( t ) = ( β v q / r 2 ) S v ( A T + I T ) / N T ( γ v + d v ) I v P ˙ ( t ) = λ 1 ( P P ) + λ 2 ( I T + R T ) / N T } t { t i , t j , t k , t n } S T ( t + ) = S T λ e Δ T ( P / P ) S T E T ( t + ) = E T + λ e Δ T ( P / P ) S T } t = t i , S T ( t + ) = S T β T p , q ( P / P ) S T ( A T + I T ) / N T E T ( t + ) = E T + β T p , q ( P / P ) S T ( A T + I T ) / N T } t = t j d v ( t + ) = [ 1 + Θ q ] ( P / P ) d v } t = t k . P ( t + ) = [ 1 + μ ] P } t = t n .

3. Results

In this section, we present simulations to show the impact of different risk factors and mitigation measures on the transmission of Sharka disease in cherry trees. Figure 3, Figure 4, Figure 5 and Figure 6 show the dynamics of different states of the trees with respect to time.
Figure 3, Figure 4 and Figure 5a and the black curves in Figure 6 each show a base case scenario, with parameters as in Table 2. In particular, we assume in this scenario that infection is produced only by initial grafting, when the orchard is established ( t i = 0 ), such that 1% of the total of 5% grafts made per unit of time were infected. Also, this base scenario considers pruning at the beginning of winter and summer ( t j = 90 and t j = 270 each year, respectively), with a higher transmission rate through pruning in winter than in summer; vector control measures applied at the beginning of spring and summer ( t k = 180 and t k = 270 each year, respectively), with a greater impact on vector mortality in spring than in summer; and awareness campaigns at the beginning of winter ( t n = 90 each year). Finally, vector–tree transmission rates with r = 2 are considered (see Equations (1) and (4)), as well as the fact that infected symptomatic trees were uprooted on average after two years of presenting symptoms ( γ I = 1 / 720 ) and asymptomatic trees were not detected ( γ A = 0 ). All other sub-figures are compared to this base case scenario.
Figure 3b,c study the effect of vector arrival in a disease-free orchard, where 1 % and 5 % of vectors arrive infected, respectively. In particular, Figure 3d compares the dynamics of the total number of infectious stone fruit trees ( I T + A T ) over time, for an initial arrival of 1 % , 5 % , 10 % , 20 % , and 50 % , showing, as expected, a clear increase in total tree infection with an increase in the prevalence of initial disease in vectors. In any case, the impact of the arrival of infected vectors is less than the infection produced by 1 % infected grafts (light blue curve in Figure 3a).
Figure 4b,c show the dynamics of the tree states for a scenario where the vector–tree transmission rates are not reduced (vector movement is not considered, r = 1 ), and a scenario where it is reduced by a factor of 1 / 9 (vector movement becomes more local and limits widespread transmission, r = 3 ), respectively. Figure 4d summarizes the dynamics of the total number of infectious trees ( I T + A T ) in each case: r = 1 (widespread transmission), r = 2 (base case transmission), r = 3 (local spread transmission). As expected, we observe a reduction in the prevalence of the disease in trees with less vector movements. This emphasizes the importance of preventing the development of winged-form aphids by early detection and control of overpopulated trees.
Figure 5 discusses the dynamics of tree states for different removal strategies for infectious trees. Figure 5a,b show the dynamics when only infected symptomatic trees are detected and removed after two years and one year on average, respectively. Figure 5c shows the dynamics when symptomatic and asymptomatic trees are detected and removed on average within one year after the onset of the disease. Figure 5d illustrates, in one plot, the curves of the total number of infectious trees ( I T + A T ) for the three scenarios described previously (from Figure 5a–c). We can observe a significant approximate reduction of 80% by the third year in the total number of infectious trees compared to the base case (red vs black curve) when detecting and uprooting not only symptomatic trees, but also asymptomatic trees early after the onset of the disease.
Figure 6 shows the impact that improving risk perception can have on the disease. We observe that after increasing the speed of reaction ( λ 2 ) and decreasing the resistance to change ( λ 1 )—parameters that appear in Equation (3) of the risk perception variable P ( t ) —the percentage of infected trees decreases by approximately 25% (red curve) by the third year (Figure 6a), compared to the base case (black curve). Regarding the effectiveness of prevention campaigns (increasing μ ), which are applied at instants t n , and keeping λ 1 and λ 2 fixed, we can observe from Figure 6b that the percentage of infected trees by the third year decreases by approximately 35% (black vs red curve). Finally, in the case that campaigns are effective (increasing μ at certain instants) and resistance to change improves along with the speed of reaction in people (decreasing λ 1 and increasing λ 2 , respectively, remaining in time), the decrease in infected trees is significant, reaching a reduction in infected trees of more than 50% by the third year (see Figure 6c, black vs red curve).

4. Discussion and Conclusions

We have proposed a mathematical model of impulsive differential equations to express the dynamics of the disease among trees and vectors. We obtained results that can serve as input for decision making on mitigation measures against the spread of the virus. Grafting, pruning, and vector control occur as pulses at certain instants and depend on a dynamic human behavior toward the presence of the disease. Human behavior is, in particular, affected by pulse disease awareness campaigns [22,23,24]. In addition, we consider the movement and seasonality of vectors that affect disease transmission. Although several aphid species have been identified as Sharka’s vectors [44], here we used Myzus persicae as the main aphid vector, considering that it is the main aphid present in Chilean Prunus orchards [45].
In the proposed model, three routes of Sharka infection are considered between trees: grafting, pruning, and through vectors (aphids). Our results show that grafting with infected material can be critical for the spread of the virus, even more so than the arrival of infectious vectors (see Figure 3a,d). Therefore, increasing people’s risk perception for the disease that reduces grafting with infectious vegetative material could significantly lower disease prevalence (see Figure 6 black curve vs blue and red curve). Hence, agricultural management strategies are suggested to be developed to monitor and test vegetative material for grafting, and awareness campaigns that motivate the agricultural community to purchase grafting material from certified sellers are proposed to be evaluated. Additionally, the arrival of infectious vectors is also a significant problem that increases disease prevalence (see Figure 3d). Although there are orchards with rigorous surveillance and management for PPV vectors, neighboring farms may not provide the same conditions. The latter become a constant reservoir of the disease, which can be easily transferred through vector mobility to Sharka-free orchards.
Climatic conditions, orchard maintenance, and vegetation around the property, among others, predetermine the movement of vectors. Our results describe the impact of vector movement on disease spread (see Figure 4). Thus, to avoid vector migration to other trees, preventing the development of winged-form vectors may be critical in keeping the disease controlled. Since vector overpopulation, as well as other environmental factors, favor the appearance of winged-form aphids and aphid movement [38,39], applying vector control measures is essential for mitigation and prevention measures.
Trees, in general, show symptoms of possible infection through their leaves. However, detecting symptoms by visual inspection with the naked eye alone is complex and certainly misses asymptomatic infectious trees. The results of our model show the need to detect infectious trees (symptomatic and asymptomatic) in the short term; otherwise, the spread of the disease within the orchard may increase significantly (Figure 5), which can produce considerable economic losses. Therefore, growers should have simple and feasible technological monitoring mechanisms to detect the virus in their orchards, especially to detect asymptomatic trees, which are a constant reservoir of the disease. From an agronomic point of view, the use of molecular tools to provide an early diagnostic of the disease in the orchard should be implemented routinely [46,47]. Another diagnostic tool that needs to be studied further is the use of drones to detect symptomatic trees, as these technologies are becoming widely used and have already been useful for other sweet cherry diseases [48].
The risk perception of individuals in agricultural communities towards this disease depends on multiple factors: social, psychological, and economic, among others. Our model considers a change in risk perception according to a dynamic law (see Equation (3)), in which we have tried to encompass these factors into two: the resistance to change while reducing risk perception ( λ 1 ) and the speed of reaction to disease prevalence ( λ 2 ), while increasing risk perception. Resistance to change can be associated with work, health, housing, education, food security, beliefs, etc., subject to the unique circumstances around farmers’ goals and agricultural contexts [22]. The reaction speed can be related to awareness campaigns, such as information via the media, educational interventions, and socialization, among others. In our model, risk perception impacts the efficiency and responsible implementation of vector control, pruning and grafting, and, hence, it impacts the dynamics of the disease (see system (5)). When a population experiences a lower resistance to change and a higher reaction speed, we can observe a decrease in disease prevalence (see Figure 6a). In addition, we have incorporated into our model the effect of awareness campaigns that help to increase risk perception (through λ 2 ). Our results show the impact of these campaigns on the number of infectious trees, reducing it considerably (see Figure 6b). In particular, when awareness campaigns have a long-lasting effect—reducing in the long-term the resistance to change and increasing the speed of reaction—we observe a reduction in the curve representing infectious trees by more than 50% of the prevalence expected without awareness campaigns in a three-year period (see Figure 6c). A higher risk perception for this disease may encourage growers to implement strategies to prevent its transmission. Therefore, we hope that after awareness campaigns that impact risk perception, aphid control [40], virus reservoirs such as weeds [49] and wild Prunus will be managed [50], and, last but not least, we hope that growers understand the importance of healthy plant material for the establishment of the Prunus orchard [51].
In the model presented through impulsive differential equations, environmental, economic, and dispersion factors, among others, were not involved, but may affect disease transmission dynamics. For example, it would be of interest to include aphid characteristics that can affect propagation, such as that winged aphids may facilitate plant infection [52] or that aphid abundance and feeding may be determined by the presence of natural predators [39]. Therefore, to evaluate the impact of these variables and others on the spread of the virus would enrich future work on the transmission and control of Sharka disease.

Author Contributions

Conceptualization, J.P.G.-J., K.V.-G., K.V.-P. and M.C.G.C.; methodology, J.P.G.-J. and K.V.-G.; software, J.P.G.-J.; validation, J.P.G.-J., K.V.-G., M.C.G.C., K.V.-P., L.M.P. and G.C.; formal analysis, J.P.G.-J., K.V.-G. and M.C.G.C.; investigation, J.P.G.-J., K.V.-G. and M.C.G.C.; resources, J.P.G.-J., K.V.-G., M.C.G.C. and K.V.-P.; data curation, J.P.G.-J., K.V.-G., M.C.G.C. and K.V.-P.; writing—original draft preparation, J.P.G.-J., K.V.-G., M.C.G.C. and K.V.-P.; writing—review and editing, J.P.G.-J., K.V.-G., M.C.G.C. and G.C.; visualization, J.P.G.-J., K.V.-G., M.C.G.C. and G.C.; supervision, J.P.G.-J., K.V.-G. and M.C.G.C.; project administration, J.P.G.-J., K.V.-G. and M.C.G.C.; funding acquisition, J.P.G.-J., K.V.-G., M.C.G.C. and K.V.-P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Universidad Católica del Maule (Fondos para Publicaciones en Revistas Open Access 2023) and Universidad Adolfo Ibáñez (Programa de Apoyo a la Investigación UAI 2023). Partially supported by ANID through Subvención a Instalación en la Academia, grant number SA77210040, and Fomento a la Vinculación Internacional, grant number FOVI220194.

Data Availability Statement

Not applicable.

Acknowledgments

We would like to thank the Agro-Math laboratory from the UCM-UAI, the EME lab and all its team and finally to thank the Universidad Católica del Maule, Universidad Adolfo Ibáñez, Tarapacá University, and Georgia State University that supported this research work.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Clemente-Moreno, M.J.; Hernández, J.A.; Diaz-Vivancos, P. Sharka: How do plants respond to Plum pox virus infection? J. Exp. Bot. 2014, 66, 25–35. [Google Scholar] [CrossRef] [PubMed]
  2. Rubio, M.; García-Ibarra, A.; Dicenta, F.; Martínez-Gómez, P. Plum pox virus (sharka) sensitivity in Prunus salicina and Prunus cerasifera cultivars against a Dideron-type isolate. Plant Breed. 2011, 130, 283–286. [Google Scholar] [CrossRef]
  3. Invasive Species Compendium. Plum Pox Virus (Sharka). Available online: https://www.cabi.org/isc/datasheet/42203 (accessed on 27 October 2022).
  4. Madariaga, V.M. La Enfermedad de Sharka: Descripción y Situación en Chile; Informativo Instituto de Investigación Agropecuarias (INIA) La Platina: Santiago, Chile, 2020; [en línea]; Available online: https://hdl.handle.net/20.500.14001/5011 (accessed on 13 August 2023).
  5. Herrera, G. Investigations of the Plum pox virus in Chile in the past 20 years. Chil. J. Agric. Res. 2013, 73, 60–65. [Google Scholar] [CrossRef]
  6. Llácer, G.; Cambra, M. Hosts and symptoms of Plum pox virus: Fruiting Prunus species. EPPO Bull. 2006, 36, 219–221. [Google Scholar] [CrossRef]
  7. Damsteegt, V.; Scorza, R.; Stone, A.; Schneider, W.; Webb, K.; Demuth, M.; Gildow, F. Prunus host range of Plum pox virus (PPV) in the United States by aphid and graft inoculation. Plant Dis. 2007, 91, 18–23. [Google Scholar] [CrossRef]
  8. Gougherty, A.; Nutter Jr, F. Impact of Eradication Programs on the Temporal and Spatial Dynamics of Plum pox virus on Prunus spp. in Pennsylvania and Ontario, Canada. Plant Dis. 2015, 99, 593–603. [Google Scholar] [CrossRef] [PubMed]
  9. D, A. Chile Launches 2022/23 Cherry Season, Expects $3 Billion in Exports to China, Produce Report. 2022. Available online: https://www.producereport.com/article/chile-launches-202223-cherry-season-expects-3-billion-exports-china (accessed on 14 August 2023).
  10. Jeger, M.; Madden, L.; Van Den Bosch, F. Plant virus epidemiology: Applications and prospects for mathematical modeling and analysis to improve understanding and disease control. Plant Dis. 2018, 102, 837–854. [Google Scholar] [CrossRef] [PubMed]
  11. Arias, J.H.; Gómez-Gardenes, J.; Meloni, S.; Estrada, E. Epidemics on plants: Modeling long-range dispersal on spatially embedded networks. J. Theor. Biol. 2018, 453, 1–13. [Google Scholar] [CrossRef]
  12. Tankam-Chedjou, I.; Touzeau, S.; Mailleret, L.; Tewa, J.J.; Grognard, F. Modelling and control of a banana soilborne pest in a multi-seasonal framework. Math. Biosci. 2020, 322, 108324. [Google Scholar] [CrossRef]
  13. Savary, S.; Nelson, A.; Willocquet, L.; Pangga, I.; Aunario, J. Modeling and mapping potential epidemics of rice diseases globally. Crop Prot. 2012, 34, 6–17. [Google Scholar] [CrossRef]
  14. Goleniewski, G. Modelling cultivar mixtures using SEIR compartmental models. Biom. J. 1996, 38, 281–297. [Google Scholar] [CrossRef]
  15. Fishman, S.; Marcus, R.; Talpaz, H.; Bar-Joseph, M.; Oren, Y.; Salomon, R.; Zohar, M. Epidemiological and economic models for spread and control of citrus tristeza virus disease. Phytoparasitica 1983, 11, 39–49. [Google Scholar] [CrossRef]
  16. Fishman, S.; Marcus, R. A model for spread of plant disease with periodic removals. J. Math. Biol. 1984, 21, 149–158. [Google Scholar] [CrossRef]
  17. Chan, M.S.; Jeger, M.J. An analytical model of plant virus disease dynamics with roguing and replanting. J. Appl. Ecol. 1994, 31, 413–427. [Google Scholar] [CrossRef]
  18. Jeger, M.; Van Den Bosch, F.; Madden, L.; Holt, J. A model for analysing plant-virus transmission characteristics and epidemic development. Math. Med. Biol. A J. IMA 1998, 15, 1–18. [Google Scholar] [CrossRef]
  19. Jackson, M.; Chen-Charpentier, B.M. Modeling plant virus propagation with seasonality. J. Comput. Appl. Math. 2019, 345, 310–319. [Google Scholar] [CrossRef]
  20. Chen-Charpentier, B. Stochastic modeling of plant virus propagation with biological control. Mathematics 2021, 9, 456. [Google Scholar] [CrossRef]
  21. Rimbaud, L.; Dallot, S.; Gottwald, T.; Decroocq, V.; Jacquot, E.; Soubeyrand, S.; Thébaud, G. Sharka Epidemiology and Worldwide Management Strategies: Learning Lessons to Optimize Disease Control in Perennial Plants. Annu. Rev. Phytopathol. 2015, 53, 357–378. [Google Scholar] [CrossRef]
  22. Ritter, C.; Jansen, J.; Roche, S.; Kelton, D.; Adams, C.L.; Orsel, K.; Erskine, R.J.; Benedictus, G.; Lam, T.J.G.M.; Barkema, H.W. Invited review: Determinants of farmers’ adoption of management-based strategies for infectious disease prevention and control. J. Dairy Sci. 2016, 100, 3329–3347. [Google Scholar] [CrossRef]
  23. Hilaire, J.; Tindale, S.; Jones, G.; Pingarron-Cardenas, G.; Bacnik, K.; Ojo, M.; Frewer, L. Risk perception associated with an emerging agri-food risk in Europe: Plant viruses in agriculture. Agric. Food Secur. 2022, 11, 21. [Google Scholar] [CrossRef]
  24. Abtew, A.; Niassy, S.; Affognon, H.; Subramanian, S.; Kreiter, S.; Tropea Garzia, G.; Martin, T. Farmers’ knowledge and perception of grain legume pests and their management in the Eastern province of Kenya. Crop Prot. 2016, 87, 90–97. [Google Scholar] [CrossRef]
  25. Benchohra, M.; Henderson, J.; Ntouyas, S. Impulsive Differential Equations and Inclusions; Contemporary Mathematics and Its Applications; Hindawi Publishing Corporation: London, UK, 2006. [Google Scholar]
  26. Córdova-Lepe, F.; Pinto, M.; González-Olivares, E. A new class of differential equations with impulses at instants dependent on preceding pulses. Applications to management of renewable resources. Nonlinear Anal. Real World Appl. 2012, 13, 2313–2322. [Google Scholar] [CrossRef]
  27. Gutiérrez-Jara, J.P.; Córdova-Lepe, F.; Muñoz-Quezada, M.T.; Chowell, G. Pesticide application, educational treatment and infectious respiratory diseases: A mechanistic model with two impulsive controls. PLoS ONE 2020, 15, e0243048. [Google Scholar] [CrossRef] [PubMed]
  28. Cabrera, M.; Córdova-Lepe, F.; Gutiérrez-Jara, J.P.; Vogt-Geisse, K. An SIR-type epidemiological model that integrates social distancing as a dynamic law based on point prevalence and socio-behavioral factors. Sci. Rep. 2021, 11, 10170. [Google Scholar] [CrossRef]
  29. Gutiérrez-Jara, J.P.; Vogt-Geisse, K.; Cabrera, M. Collateral effects of insecticide-treated nets on human and environmental safety in an epidemiological model for malaria with human risk perception. Int. J. Environ. Res. Public Health 2022, 19, 16327. [Google Scholar] [CrossRef] [PubMed]
  30. Gutiérrez-Jara, J.P.; Vogt-Geisse, K.; Cabrera, M.; Córdova-Lepe, F.; Muñoz-Quezada, M.T. Effects of human mobility and behavior on disease transmission in a COVID-19 mathematical model. Sci. Rep. 2022, 12, 10840. [Google Scholar] [CrossRef] [PubMed]
  31. Moreno, A.; Fereres, A.; Cambra, M. Quantitative estimation of plum pox virus targets acquired and transmitted by a single Myzus persicae. Arch. Virol. 2009, 154, 1391–1399. [Google Scholar] [CrossRef]
  32. Gao, S.; Xia, L.; Liu, Y.; Xie, D. A Plant Virus Disease Model with Periodic Environment and Pulse Roguing. Stud. Appl. Math. 2016, 136, 357–381. [Google Scholar] [CrossRef]
  33. Nicaise, V. Crop immunity against viruses: Outcomes and future challenges. Front. Plant Sci. 2014, 5, 660. [Google Scholar] [CrossRef]
  34. Farhan, S.S.; Darwesh, D.; Nsr, M.M. Grafting “Pioneer” plum cultivar on two different rootstocks 1-Effect of different rootstock and grafting technique on “Pioneer” plum cultivar seedlings. Egypt. J. Agric. Res. 2022, 100, 653–660. [Google Scholar]
  35. Tedesco, S.; Fevereiro, P.; Kragler, F.; Pina, A. Plant grafting and graft incompatibility: A review from the grapevine perspective. Sci. Hortic. 2022, 299, 111019. [Google Scholar] [CrossRef]
  36. Svanes, E.; Johnsen, F.M. Environmental life cycle assessment of production, processing, distribution and consumption of apples, sweet cherries and plums from conventional agriculture in Norway. J. Clean. Prod. 2019, 238, 117773. [Google Scholar] [CrossRef]
  37. Li, S.; Travadon, R.; Trouillas, F. Seasonal susceptibility of sweet cherry pruning wounds to Calosphaeria pulchella, Cytospora sorbicola and Eutypa lata. Plant Dis. 2023. [Google Scholar] [CrossRef] [PubMed]
  38. Bosch, M.J.; Ives, A.R. Identifying environmental factors affecting the production of pea aphid dispersal morphs in field populations. Ecol. Entomol. 2023, 48, 597–609. [Google Scholar] [CrossRef]
  39. Kansman, J.T.; Nersten, C.E.; Hermann, S.L. Smelling danger: Lady beetle odors affect aphid population abundance and feeding, but not movement between plants. Basic Appl. Ecol. 2023, 71, 1–8. [Google Scholar] [CrossRef]
  40. Hazir, A.; Yurtmen, M.; Fidan, H. Potential aphid (Hemiptera: Aphididae) vectors of Plum-pox virus (virus: Potyviridae) and status of Sharka disease in stone fruit orchards in the East Mediterranean region of Turkey. J. Agric. Sci. 2021, 27, 484–492. [Google Scholar] [CrossRef]
  41. Saleesha, F.M.A.; Kennedy, J.; Rajabaskar, D.; Geethalakshmi, V. Temperature dependent development of Myzus persicae Sulzer (Hemiptera: Aphididae) in Cauliflower (Brassica oleracea var. botrytis). Agric. Res. 2022, 11, 488–498. [Google Scholar] [CrossRef]
  42. Bernardi, D.; Garcia, M.S.; Botton, M.; Nava, D.E. Biology and fertility life table of the green aphid Chaetosiphon fragaefolli on strawberry cultivars. J. Insect Sci. 2012, 12, 28. [Google Scholar] [CrossRef]
  43. Muñoz-Quezada, M.T.; Lucero, B.A.; Gutiérrez-Jara, J.P.; Buralli, R.J.; Zúñiga-Venegas, L.; Muñoz, M.P.; Ponce, K.V.; Iglesias, V. Longitudinal exposure to pyrethroids (3-PBA and trans-DCCA) and 2, 4-D herbicide in rural schoolchildren of Maule region, Chile. Sci. Total Environ. 2020, 749, 141512. [Google Scholar] [CrossRef]
  44. Cinar, C.T.; Gazel, M.; Kaya, K.; Olmos, A.; Caglayan, K. Susceptibility of different prunus rootstocks to natural infection of plum pox virus-Turkey (PPV-T) in Central Anatolia. Physiol. Mol. Plant Pathol. 2022, 119, 101837. [Google Scholar] [CrossRef]
  45. Yang, Q.; Umina, P.A.; Wei, S.; Bass, C.; Yu, W.; Robinson, K.L.; Gill, A.; Zhan, D.; Ward, S.E.; van Rooyen, A.; et al. Diversity and Regional Variation of Endosymbionts in the Green Peach Aphid, Myzus persicae (Sulzer). Diversity 2023, 15, 206. [Google Scholar] [CrossRef]
  46. De Mori, G.; Savazzini, F.; Geuna, F. Chapter 10—Molecular tools to investigate Sharka disease in Prunus species. In Applied Plant Biotechnology for Improving Resistance to Biotic Stress; Poltronieri, P., Hong, Y., Eds.; Academic Press: Cambridge, MA, USA, 2020; pp. 203–223. [Google Scholar] [CrossRef]
  47. Espinoza, C.; Bascou, B.; Calvayrac, C.; Bertrand, C. Deciphering Prunus Responses to PPV Infection: A Way toward the Use of Metabolomics Approach for the Diagnostic of Sharka Disease. Metabolites 2021, 11, 465. [Google Scholar] [CrossRef] [PubMed]
  48. Chaschatzis, C.; Karaiskou, C.; Mouratidis, E.G.; Karagiannis, E.; Sarigiannidis, P.G. Detection and Characterization of Stressed Sweet Cherry Tissues Using Machine Learning. Drones 2022, 6, 3. [Google Scholar] [CrossRef]
  49. Kil, E.J.; Ho, P.T.; Fadhila, C.; Lal, A.; Vo, T.T.B.; Kim, M.; Lee, S. Plum pox virus: Diagnosis and spread inhibition by weed control. J. Plant Dis. Prot. 2021, 128, 1091–1099. [Google Scholar] [CrossRef]
  50. Collum, T.D.; Stone, A.L.; Sherman, D.J.; Damsteegt, V.D.; Schneider, W.L.; Rogers, E.E. Viral Reservoir Capacity of Wild Prunus Alternative Hosts of Plum Pox Virus Through Multiple Cycles of Transmission and Dormancy. Plant Dis. 2022, 106, 101–106. [Google Scholar] [CrossRef]
  51. Zhou, J.; Xing, F.; Wang, H.; Li, S. Occurrence, Distribution, and Genomic Characteristics of Plum Pox Virus Isolates from Common Apricot (Prunus armeniaca) and Japanese Apricot (Prunus mume) in China. Plant Dis. 2021, 105, 3474–3480. [Google Scholar] [CrossRef] [PubMed]
  52. Guo, H.; Zhang, Y.; Li, B.; Li, C.; Shi, Q.; Zhu-Salzman, K.; Ge, F.; Sun, Y. Salivary carbonic anhydrase II in winged aphid morph facilitates plant infection by viruses. Proc. Natl. Acad. Sci. USA 2023, 120, e2222040120. [Google Scholar] [CrossRef]
Figure 1. Schematics of the disease dynamics for trees. Table 1 and Table 2 describe, respectively, the variables and parameters used.
Figure 1. Schematics of the disease dynamics for trees. Table 1 and Table 2 describe, respectively, the variables and parameters used.
Plants 12 03442 g001
Figure 2. Schematics of the disease dynamics for vectors. Table 1 and Table 2 describe, respectively, the variables and parameters used.
Figure 2. Schematics of the disease dynamics for vectors. Table 1 and Table 2 describe, respectively, the variables and parameters used.
Plants 12 03442 g002
Figure 3. Effects of grafting with infected material and arrival of infectious vectors on disease progress. (ac) show the curves representing exposed trees ( E T , blue), symptomatic infectious trees ( I T , red), asymptomatic infectious trees ( A T , green), removed trees ( R T , yellow), and the total number of infectious trees ( I T + A T , light blue), for four scenarios: (a) infection is only produced by inserting 1% of infected grafts out of 5 % of grafts made initially ( t i = 0 ); (b,c) infection is only produced by the initial arrival of a percentage of infected vectors ( I v ( 0 ) ) and no grafting with infected material occurs ( Δ T = 0 ). Initial conditions were taken: in (a) S T ( 0 ) = 1 , E T ( 0 ) = I T ( 0 ) = A T ( 0 ) = R T ( 0 ) = 0 ; S v ( 0 ) = 1 , I v ( 0 ) = 0 , P ( 0 ) = P ; in (b) S v ( 0 ) = 0.99 , I v ( 0 ) = 0.01 and the remaining as in (a); in (c) S v ( 0 ) = 0.95 , I v ( 0 ) = 0.05 and the remaining as in (a). (d) shows the total number of infectious trees ( I T + A T ) for different initial numbers of infectious vectors ( I v ( 0 ) ), under the same conditions as in (b,c). All other parameters used in the simulations are as in Table 2.
Figure 3. Effects of grafting with infected material and arrival of infectious vectors on disease progress. (ac) show the curves representing exposed trees ( E T , blue), symptomatic infectious trees ( I T , red), asymptomatic infectious trees ( A T , green), removed trees ( R T , yellow), and the total number of infectious trees ( I T + A T , light blue), for four scenarios: (a) infection is only produced by inserting 1% of infected grafts out of 5 % of grafts made initially ( t i = 0 ); (b,c) infection is only produced by the initial arrival of a percentage of infected vectors ( I v ( 0 ) ) and no grafting with infected material occurs ( Δ T = 0 ). Initial conditions were taken: in (a) S T ( 0 ) = 1 , E T ( 0 ) = I T ( 0 ) = A T ( 0 ) = R T ( 0 ) = 0 ; S v ( 0 ) = 1 , I v ( 0 ) = 0 , P ( 0 ) = P ; in (b) S v ( 0 ) = 0.99 , I v ( 0 ) = 0.01 and the remaining as in (a); in (c) S v ( 0 ) = 0.95 , I v ( 0 ) = 0.05 and the remaining as in (a). (d) shows the total number of infectious trees ( I T + A T ) for different initial numbers of infectious vectors ( I v ( 0 ) ), under the same conditions as in (b,c). All other parameters used in the simulations are as in Table 2.
Plants 12 03442 g003
Figure 4. Effects of vector mobility on disease progress. (ac) show the curves representing exposed trees ( E T , blue), infectious symptomatic trees ( I T , red), infectious asymptomatic trees ( A T , green), removed trees ( R T , yellow), and the total number of infectious trees ( I T + A T , light blue) with respect to time, for r = 1 , r = 2 , and r = 3 , respectively. (d) shows the total number of infectious trees ( I T + A T ) with respect to time and for r = 1 (blue), r = 2 (black), and r = 3 (red). The initial conditions were taken as in Figure 3a and all other parameters not specified here and used in the simulations are as in Table 2.
Figure 4. Effects of vector mobility on disease progress. (ac) show the curves representing exposed trees ( E T , blue), infectious symptomatic trees ( I T , red), infectious asymptomatic trees ( A T , green), removed trees ( R T , yellow), and the total number of infectious trees ( I T + A T , light blue) with respect to time, for r = 1 , r = 2 , and r = 3 , respectively. (d) shows the total number of infectious trees ( I T + A T ) with respect to time and for r = 1 (blue), r = 2 (black), and r = 3 (red). The initial conditions were taken as in Figure 3a and all other parameters not specified here and used in the simulations are as in Table 2.
Plants 12 03442 g004
Figure 5. Effects of detection and removal of infected trees on disease progress. (ac) show the curves representing exposed trees ( E T , blue), infectious symptomatic trees ( I T , red), infectious asymptomatic trees ( A T , green), removed trees ( R T , yellow), and the total number of infectious trees ( I T + A T , light blue) with respect to time. (a,b) represent scenarios where only symptomatic trees ( I T ) were removed after two years and one year on average, respectively, after presenting symptoms ( 1 / γ T = 720 and 1 / γ T = 360 , respectively). (c) represents a scenario in which symptomatic and asymptomatic trees were removed after an average of one year after becoming infectious ( 1 / γ T = 360 and 1 / γ A = 360 , respectively). (d) shows scenarios (ac) for the total number of infectious trees ( I T + A T ). The initial conditions were taken as in Figure 3a and all other parameters not specified here and used in the simulations are as in Table 2.
Figure 5. Effects of detection and removal of infected trees on disease progress. (ac) show the curves representing exposed trees ( E T , blue), infectious symptomatic trees ( I T , red), infectious asymptomatic trees ( A T , green), removed trees ( R T , yellow), and the total number of infectious trees ( I T + A T , light blue) with respect to time. (a,b) represent scenarios where only symptomatic trees ( I T ) were removed after two years and one year on average, respectively, after presenting symptoms ( 1 / γ T = 720 and 1 / γ T = 360 , respectively). (c) represents a scenario in which symptomatic and asymptomatic trees were removed after an average of one year after becoming infectious ( 1 / γ T = 360 and 1 / γ A = 360 , respectively). (d) shows scenarios (ac) for the total number of infectious trees ( I T + A T ). The initial conditions were taken as in Figure 3a and all other parameters not specified here and used in the simulations are as in Table 2.
Plants 12 03442 g005
Figure 6. Effects of risk perception on disease progress. (ac) show the curve of the total number of infectious trees ( I T + A T ) with respect to time. (a) pictures three scenarios for different values of the resistance to change parameter, λ 1 , and the rate of reaction to disease, λ 2 , from Equation (3), where the black curve represents the smallest level of risk perception, and the red curve the largest. (b) shows three scenarios for a varying factor of increase in the risk perception, μ , at instants t n (see Equation (5)), such that the black curve represents the lowest and the red curve the highest effectiveness of awareness campaigns. (c) shows scenarios (a,b) combined. The initial conditions were taken as in Figure 3a, and all other parameters not specified here and used in the simulations are in Table 2.
Figure 6. Effects of risk perception on disease progress. (ac) show the curve of the total number of infectious trees ( I T + A T ) with respect to time. (a) pictures three scenarios for different values of the resistance to change parameter, λ 1 , and the rate of reaction to disease, λ 2 , from Equation (3), where the black curve represents the smallest level of risk perception, and the red curve the largest. (b) shows three scenarios for a varying factor of increase in the risk perception, μ , at instants t n (see Equation (5)), such that the black curve represents the lowest and the red curve the highest effectiveness of awareness campaigns. (c) shows scenarios (a,b) combined. The initial conditions were taken as in Figure 3a, and all other parameters not specified here and used in the simulations are in Table 2.
Plants 12 03442 g006
Table 1. Description of the variables of the model given in system (5).
Table 1. Description of the variables of the model given in system (5).
VariableDescriptionUnits *
S T ( t ) Susceptible trees at time tT
E T ( t ) Exposed trees at time tT
I T ( t ) Infectious symptomatic trees at time tT
A T ( t ) Infectious asymptomatic trees at time tT
R T ( t ) Removed (uprooted) trees at time tV
S v ( t ) Susceptible vectors at time tV
I v ( t ) Infectious vectors at time tV
P ( t ) Human risk perception to Sharka disease at time tR
* T:= tree, V:= vector, R:= Risk perception.
Table 2. Description of the parameters of the model given in system (5).
Table 2. Description of the parameters of the model given in system (5).
Param.DescriptionUnits *BaselineValueRef. +
Tree
β T v , q Transmission rate to tree through vector, on season q . d 1 [ 0 , 1 5 ] 1 5 [ 0 , 1 8 , 1 , 1 2 ] [31]
β T p , q Transmission rate to tree through pruning, on season q . d 1 [ 0 , 1 100 ] 1 100 [ 0 , 1 , 0 , 1 3 ]
1 / α Average virus-PPV incubation time.d 31 0.17 [32]
1 / γ I Average time until a symptomatic infectious tree is uprooted.d 1 720 C.A.
1 / γ A Average time until an asymptomatic infectious tree is uprooted.d 0C.A.
δ a , i Transition rate from symptomatic to asymptomatic infectious. d 1 [ 1 90 , 1 60 ] 1 70 [6]
δ i , a Transition rate from asymptomatic to symptomatic infectious. d 1 1 240 [6]
Λ T Replacement rate of removed trees. d 1 1 360 [33]
Λ Recruitment rate. T d 1 N T d T , N T = 1 1 7200
Δ T Grafting rate. d 1 5 100 [34,35]
λ e Proportion of infected grafts.Unitless 0 , 1 1 100 C.A.
d T Natural mortality rate d 1 [ 1 35 , 1 15 ] 1 365 1 7200 [36]
t i Time instants where grafting occurs.d 0[34,35]
t j Time instants where pruning occurs.d 90 ( 1 + 2 k ) , k N 0 [37]
Vector
β v q Transmission rate to vector, on season q . d 1 [ 0 , 13 100 ] 13 100 [ 1 , 1 , 1 , 1 ] [31]
1 / r 2 Vector–host transmission rate reduction due to vector movement.Unitless 1 4 [38,39]
1 / γ v Average time a vector remains infectious.d 0.1 [40]
Λ v Recruitment rate. V d 1 N v d v , N v = 1 1 30
d v Natural mortality rate. d 1 1 30 [41,42]
Θ q Mortality rate increase factor, on season q .Unitless [ 0 , 0 , 3 10 , 3 100 ] C.A.
t k Time instants where vector control occurs.d 180 ( 2 k + 1 ) ;
90 + 180 ( 2 k + 1 ) ,
k N 0 .
[33]
Human
λ 1 Rate of resistance to change risk perception. d 1 0 , 1 0.01[29,30,43]
λ 2 Per capita reaction rate to change risk perception. R d 1 0 , 1 0.01[29,30,43]
P Natural risk perception.R 0.5[29,30,43]
μ Risk perception increase factor.Unitless 0.1C.A.
t n Time instants where educational prevention campaigns occur.d 90 ( 1 + 4 k ) , k N 0 C.A.
* d:= day, T:= tree, V:= vector. ** q = 1(autumn), 2(winter), 3(spring), 4(summer), + C.A. = Chosen by author.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Gutiérrez-Jara, J.P.; Vogt-Geisse, K.; Correa, M.C.G.; Vilches-Ponce, K.; Pérez, L.M.; Chowell, G. Modeling the Impact of Agricultural Mitigation Measures on the Spread of Sharka Disease in Sweet Cherry Orchards. Plants 2023, 12, 3442. https://doi.org/10.3390/plants12193442

AMA Style

Gutiérrez-Jara JP, Vogt-Geisse K, Correa MCG, Vilches-Ponce K, Pérez LM, Chowell G. Modeling the Impact of Agricultural Mitigation Measures on the Spread of Sharka Disease in Sweet Cherry Orchards. Plants. 2023; 12(19):3442. https://doi.org/10.3390/plants12193442

Chicago/Turabian Style

Gutiérrez-Jara, Juan Pablo, Katia Vogt-Geisse, Margarita C. G. Correa, Karina Vilches-Ponce, Laura M. Pérez, and Gerardo Chowell. 2023. "Modeling the Impact of Agricultural Mitigation Measures on the Spread of Sharka Disease in Sweet Cherry Orchards" Plants 12, no. 19: 3442. https://doi.org/10.3390/plants12193442

APA Style

Gutiérrez-Jara, J. P., Vogt-Geisse, K., Correa, M. C. G., Vilches-Ponce, K., Pérez, L. M., & Chowell, G. (2023). Modeling the Impact of Agricultural Mitigation Measures on the Spread of Sharka Disease in Sweet Cherry Orchards. Plants, 12(19), 3442. https://doi.org/10.3390/plants12193442

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