Research article

Modeling Citrus Huanglongbing transmission within an orchard and its optimal control

  • Received: 30 October 2019 Accepted: 23 December 2019 Published: 31 December 2019
  • Citrus Huanglongbing (HLB) is the most devastating citrus disease worldwide. In this paper, a deterministic dynamical model is proposed to explore the transmission dynamics of HLB between citrus tree and Asian citrus psyllid (ACP). Using the theory of dynamical system, the dynamics of the model are rigorously analyzed. The results show that the disease-free equilibrium is globally asymptotically stable when the basic reproduction number R0<1, and when R0>1 the system is uniformly persistent. Applying the global sensitivity analysis of R0, some parameters that have the greatest impact on HLB transmission dynamics are obtained. Furthermore, the optimal control theory is applied to the model to study the corresponding optimal control problem. Both analytical and numerical results show that: (1) the infected ACP plays a decisive role in the transmission of HLB in citrus trees, and eliminating the ACP will be helpful to curtail the spread of HLB; (2) optimal control strategy is superior to the constant control strategy in decreasing the prevalence of the diseased citrus trees, and the cost of implementing optimal control is much lower than that of the constant control strategy; and (3) spraying insecticides is more effective than other control strategies in reducing the number of ACP in the early phase of the transmission of HLB. These theoretical and numerical results may be helpful in making public policies to control HLB in orchards more effectively.

    Citation: Fumin Zhang, Zhipeng Qiu, Balian Zhong, Tao Feng, Aijun Huang. Modeling Citrus Huanglongbing transmission within an orchard and its optimal control[J]. Mathematical Biosciences and Engineering, 2020, 17(3): 2048-2069. doi: 10.3934/mbe.2020109

    Related Papers:

    [1] Xiao-Dong Yang, Ze-Fan Liao, Bin Shu, Ai-Jia Chen . Blockchain-based multi-authority revocable data sharing scheme in smart grid. Mathematical Biosciences and Engineering, 2023, 20(7): 11957-11977. doi: 10.3934/mbe.2023531
    [2] Jiayuan Zhang, Rongxin Guo, Yifan Shi, Wanting Tang . An anti-impersonation attack electronic health record sharing scheme based on proxy re-encryption and blockchain. Mathematical Biosciences and Engineering, 2024, 21(6): 6167-6189. doi: 10.3934/mbe.2024271
    [3] Zilong Liu, Jingbing Li, Jing Liu . Encrypted face recognition algorithm based on Ridgelet-DCT transform and THM chaos. Mathematical Biosciences and Engineering, 2022, 19(2): 1373-1387. doi: 10.3934/mbe.2022063
    [4] Ching-Chun Chang, Chang-Tsun Li . Algebraic secret sharing using privacy homomorphisms for IoT-based healthcare systems. Mathematical Biosciences and Engineering, 2019, 16(5): 3367-3381. doi: 10.3934/mbe.2019168
    [5] Kaimeng Chen, Chin-Chen Chang . High-capacity reversible data hiding in encrypted images based on two-phase histogram shifting. Mathematical Biosciences and Engineering, 2019, 16(5): 3947-3964. doi: 10.3934/mbe.2019195
    [6] Nyamsuren Vaanchig, Zhiguang Qin . Public key encryption with temporary and fuzzy keyword search. Mathematical Biosciences and Engineering, 2019, 16(5): 3914-3935. doi: 10.3934/mbe.2019193
    [7] Qi Wang, John Blesswin A, T Manoranjitham, P Akilandeswari, Selva Mary G, Shubhangi Suryawanshi, Catherine Esther Karunya A . Securing image-based document transmission in logistics and supply chain management through cheating-resistant visual cryptographic protocols. Mathematical Biosciences and Engineering, 2023, 20(11): 19983-20001. doi: 10.3934/mbe.2023885
    [8] Han-Yu Lin, Tung-Tso Tsai, Hong-Ru Wu, Miao-Si Ku . Secure access control using updateable attribute keys. Mathematical Biosciences and Engineering, 2022, 19(11): 11367-11379. doi: 10.3934/mbe.2022529
    [9] Shufen Niu, Wei Liu, Sen Yan, Qi Liu . Message sharing scheme based on edge computing in IoV. Mathematical Biosciences and Engineering, 2023, 20(12): 20809-20827. doi: 10.3934/mbe.2023921
    [10] Yang Zhao, Xin Xie, Xing Zhang, Yi Ding . A revocable storage CP-ABE scheme with constant ciphertext length in cloud storage. Mathematical Biosciences and Engineering, 2019, 16(5): 4229-4249. doi: 10.3934/mbe.2019211
  • Citrus Huanglongbing (HLB) is the most devastating citrus disease worldwide. In this paper, a deterministic dynamical model is proposed to explore the transmission dynamics of HLB between citrus tree and Asian citrus psyllid (ACP). Using the theory of dynamical system, the dynamics of the model are rigorously analyzed. The results show that the disease-free equilibrium is globally asymptotically stable when the basic reproduction number R0<1, and when R0>1 the system is uniformly persistent. Applying the global sensitivity analysis of R0, some parameters that have the greatest impact on HLB transmission dynamics are obtained. Furthermore, the optimal control theory is applied to the model to study the corresponding optimal control problem. Both analytical and numerical results show that: (1) the infected ACP plays a decisive role in the transmission of HLB in citrus trees, and eliminating the ACP will be helpful to curtail the spread of HLB; (2) optimal control strategy is superior to the constant control strategy in decreasing the prevalence of the diseased citrus trees, and the cost of implementing optimal control is much lower than that of the constant control strategy; and (3) spraying insecticides is more effective than other control strategies in reducing the number of ACP in the early phase of the transmission of HLB. These theoretical and numerical results may be helpful in making public policies to control HLB in orchards more effectively.


    Graphene nanoribbon (GNR) is a narrow strip of the patterned graphene sheet that has gained considerable research interest for its distinct electrical, thermal, structural, magnetic, and mechanical properties [1,2,3,4,5,6,7,8,9,10,11]. The properties like, linear E-k relationship, long mean free path (MFP), superior mobility, ideally negative temperature coefficient of resistance, and a width-dependent transport gap have made GNR a stronger alternative to the conventional copper for nano-interconnect applications [12,13]. Although GNR has very high carrier mobility [14], it is limited by different scattering mechanisms. Depending on the geometry, GNRs exhibit metallic or semiconducting property [15,16]. The transverse width of GNRs should be substantially smaller than that of the longitudinal length [17,18]. Depending upon the position of Carbon atoms at the edges, GNRs are classified as armchair and zigzag nanoribbons as depicted in Figure 1a, b. The difference in orientation at these edges of GNRs is 30° [18,19,20]. Subsequently, no contribution from GNR has been observed to energy states near the Fermi level of Carbon atoms as hydrogen atoms passivize the edge of the GNRs [19,20]. N is illustrated in Figure 1 as the number of dimers for armchair GNR as well as zigzag lines in ZZ-GNRs [12]. The width, w of the GNR has a direct relationship with the integer N, which is represented in Figure 1a, where it has been presented that, width and length of the unit cell of armchair GNRs are w=(N+1)a/2 w and aT=3a. However, in the case of ZZ-GNR, these are, w and a, w and a (w=a3+3Na/2) respectively, which has been shown in Figure 1b. The value of N is (3p+2) when an armchair GNR possesses the metallic property, and it is 3p or (3p+1), the nature is semiconducting [21]. It should be noted that p has been assumed to be an integer at this point (p = 1, 2, 3…).

    Figure 1.  GNR structure with (a) armchair and (b) zigzag edges (Reprinted (adapted) from Ref. [21] with permission of American Chemical Society).

    The bandgap shrinks with the increase of the width of the semiconducting GNR, and the metallic GNR has a large value of N [12]. The ZZ-GNR's Brillouin zone has partially flat bands due to their edge states and consequently shows zero bandgap with metallic behavior [19,20,22]. Furthermore, it is build up with three kinds of arrangements of carbon atoms; consequently, a significant difference in heat transfer capability is observed [23]. Therefore, the ZZ-metallic GNRs have attracted wide attention among researchers in the field of on-chip interconnect design [12,19,24].

    The intercalation is a process of implanting molecules or ions between the large van der Waals gap of GNR layers. It has been reported in the literature that intercalation doping in bulk material shows a significant distinction from the layered materials like GNR [25,26]. Subsequently, intercalation doping in layered material upshots significant improvement in the device and interconnect-based applications [27,28,29,30,31]. Intercalation doping has been proved as a very powerful way to tune the attributes of MLGNRs for few reasons like, (a) the associated changes of the material properties are irreversible, (b) through this process, highest possible doping can be achieved, and (c) the amount of intercalation is controllable, and periodically run-time control has been achieved through electrochemical voltage.

    With the intercalation doping, both the in-plane and interlayer (c-axis) resistances of MLGNR have been improved. Major improvement in Fermi energy, mobility, and mean free path has been achieved through Arsenic Pentafluoride (AsF5-), Lithium (Li-), Ferric Chloride (FeCl3), and Molybdenum Pentachloride (MoCl5) intercalation doping. Lithium (Li-) intercalated ultrathin graphite structure is fabricated, and improvement in Fermi level and in-plane conductivity has been experimentally observed [30]. The reported values of inter-layer conductivity (σc) of pristine and intercalation doped graphene are presented in the literature [25,26,28,29,32] as displayed in Table 1.

    Table 1.  Properties of pristine and intercalation doped MLGNRs [14,15,32].
    Properties Pristine (undoped) FeCl3-intercalated AsF5-intercalated Li-intercalated MoCl5-intercalated
    Ef (eV) 0.2 0.6 0.68 1.5 3.4
    δm (nm) 0.34 0.575 0.394 0.37 0.96
    C-axis conductivity (Ω·cm)−1 0.033 0.24 1.0 1.8 × 104 0.33

     | Show Table
    DownLoad: CSV

    The rest of the article is organized as follows. Section works on modelling and analysis of GNR based interconnects presents a detailed review on modeling and analysis of GNR based interconnects. The practical works on GNR are presented in section practical works on GNR, followed by section conclusions.

    A number of research articles have been published which demonstrate the performance of GNR interconnect in nano-scale high-speed CMOS-based design. Naeemi and Meindl [33] proposed a physics-based conductance model for graphene nanoribbon (GNR) interconnect for the first time. They also proposed a charge transport model of zigzag and armchair graphene nanoribbons [34]. Murali et al. have proposed a temperature-aware resistivity model of graphene nanoribbon interconnect [35]. Further, Naeemi and Meindl [36] have proposed a complete compact physics-based circuit model for GNR interconnect with an accurate calculation of the number of conduction channels (Nch) in a single layer GNR (SLGNR) using the close form approximation method. The equivalent single-conductor (ESC) model of MLGNR interconnect is illustrated in Figure 2 [13].

    Figure 2.  Basic equivalent single-conductor (ESC) model of MLGNR interconnect (Reprinted from Ref. [13] with permission of IEEE).

    Nasiri et al. proposed a compact formula to calculate temperature aware number of conduction channels of graphene nanoribbons with different Fermi energy [37]. However, this approach limits the computation for a given temperature and Fermi energy only. Figure 3 shows the design route of the MLGNR as an on-chip interconnect.

    Figure 3.  Schematic representation of graphene to MLGNR modeling approach.

    In [28], Banerjee et al. have proposed the conductance model of GNR with tight-binding approximation and obtained a linear E-k relationship using Landauer formalism. Here, the delay analysis in GNR interconnects has been presented as well as the performance comparison with other interconnect materials, such as copper (Cu), tungsten (W), and CNT, has been performed. Nasiri et al. have presented stability analysis in GNR interconnects [38]. The relative stability of multilayer GNR (MLGNR) interconnects is analyzed using the Nyquist stability criterion. Their analysis shows that by increasing the length and width, MLGNR interconnects become more stable. The crosstalk stability analysis using multilayer GNR interconnects is reported in [39]. Ramesh Kumar et al. have reported the delay, and stability analysis of multilayer GNR and multi-walled CNT interconnects [40].

    Kumar et al. first proposed different contact-based modeling (top-contact and side-contact) of graphene nanoribbon-based interconnect [41]. Following the route, Nisad et al. in [42] proposed a 2D resistive model for top-contact and side-contact MLGNR interconnects. It is observed that the side-contact MLGNR (SC-MLGNR) performs better than top-contact MLGNR (TC-MLGNR) in terms of delay, power-delay product (PDP), energy-delay product (EDP), etc. reported in [40,41,42,43]. But from the fabrication point of view, the TC-MLGNR is much better than SC-MLGNR reported in [41,42]. Subsequently, to improve the performance of TC-MLGNR, intercalation becomes important. The conductance modeling of intercalation doped TC-MLGNR has been presented in [22]. Intercalation doping with Arsenic pentafluoride (AsF5) in MLGNRs has been carried out by Xu et al., and substantial improvement in Fermi level consequently MFP and been observed [44].

    A self-heating model has been proposed by Jiang et al. [45] on Lithium (Li-) and Ferric chloride (FeCl3-) intercalation doped ultra-thin graphite arrangement has been fabricated as TC- interconnect and the improvement in Fermi level, and the in-plane conductivity has been experimentally observed. Das et al. proposed an electro-thermal RF model [46,47] for intercalation doped TC- as well as SC-MLGNR considering the integrated effect of distributed interconnect system based on the thermophysical property of multilayer GNR studied in the literature [48,49,50,51]. They have demonstrated superior RF performance of intercalation doped TC-MLGNR interconnect as compared to its pristine counterpart for different chip operating temperatures. Rakheja et al. [52] have compared the energy dissipation, and overall performance of graphene spin interconnects in a non-local spin-torque (NLST) circuit against the conventional CMOS-based circuit. They also reported the prospects of graphene-based interconnect that could be the replacement of copper in future microchips. The spin injection probability of electron in graphene wires has been further explored by them in [53]. The improvement of carrier transport with the help of the selective hydrogenation mechanism of edges for sub-30 nm graphene nanoribbons (GNRs) has been reported [54]. After improvement of the edges roughness, the GNRs exhibit improved transport properties. The carrier mobility has been improved up to 50% when the carrier density is cm−2 at room temperature.

    Pan et al. [55] reported the advantage of GNR interconnect over the traditional Cu-based nano-interconnect for next-generation VLSI interconnect. They have designed a 32-bit adder using MLGNR interconnect and observed a 70% improvement in EDP. They have also demonstrated the design of the ARM processor using MLGNR interconnect and achieved 22% and 15% improvement in EDP and clock frequency, respectively. It should be noted that the benefits observed in ARM processors are highly dependent on the quality, such as the MFP and the edge roughness of the GNR.

    The simultaneous switching noise (SSN) and IR-drop analysis for graphene nanoribbon-based power distribution network reported by Das et al. [56]. This paper reports the peak SSN in Cu-based interconnects 41–23% larger than that in GNR. Consequently, the peak IR-drop in GNR is 38–34% lesser than that of traditional copper wires. The propagation delay is highly affected due to SSN and IR-drop. GNR shows up to 64.49% lesser impact in the SSN generated delay in comparison with the Cu-based power networks. The reduced thickness model of MLGNR interconnect has been proposed by Bhattacharya et al. [57]. It has been demonstrated that keeping the width constant, reduced thickness of MLGNR has a lesser impact on crosstalk delay and overshoot/undershoot. The IR-drop induced temperature aware delay fault model of CNT and GNR based power distribution networks is proposed by Das et al. and Bhattacharya et al. [58,59]. It has been observed in these works that the delay faults can be reduced using CNT as well as GNR based power interconnects at longer lengths as compared to the conventional Cu-based power interconnects. The crosstalk and gate oxide reliability analysis in GNR interconnects have been investigated by Sahoo et al. [60] and Das et al. [61]. It has been observed that electric field acceleration rises significantly at a higher temperature. This set a substantial impact on gate oxide reliability on MLGNR driven interconnect system. This thermal impact of gate oxide reliability has been reported by Das et al. [62].

    The contact resistance-related issues in GNR based interconnect is very crucial in sub-nanometer technology. It may be top-contact or side-contact. The top-contact resistance for multi-layered GNR interconnect is ~4.3 kΩ reported in [41]. Leong et al. [63] reported the graphene-based devices with nickel-etched-graphene contacts, which shows considerably low contact resistance.

    Signal transmission and integrity analysis on GNR based interconnect reported by Zhao et al. [64]. Signal integrity analysis of single and multi-layered GNR interconnects has been performed based on their circuit equivalent models, with crosstalk effects characterized theoretically. It has been concluded in their work that longer SLGNR interconnect with larger width results in higher crosstalk noise than that of MLGNR or Cu-based interconnects. Cui et al. [65] have analyzed signal transmission characteristics in multilayer GNR (MLGNR) interconnects for 22 nm and 14 nm technology nodes. An ESC model has been derived, and transient response is studied in their work. Another signal transmission analysis along GNR interconnects has been reported by Xu et al. [66]. Sarkar et al. have reported the UHF behavior of GNR based interconnects. They have proposed a numerical model of GNR as an on-chip interconnect to calculate the surface impedance with anomalous skin effect in high-frequency using linear dispersion relation of graphene [67,68]. Chen and Fratini et al. [69,70] have reported the distress of low-bias resistivity due to phonon-based scattering (e.g., acoustic phonon scattering (AP) and remote interfacial phonon (RIP)) with temperature. These works imprinted a significant impact on temperature-dependent models and circuit analysis of GNR as an on-chip interconnect [37,45,47,59,62].

    The world's first graphene-based integrated circuit (IC) was developed by IBM in the year of 2011. An article on graphene-based IC design by the IBM research group is reported by Lin et al. [3]. T. J. Watson Research Center (Microelectronics Research Laboratory) and AMO Nanoelectronics group had demonstrated the possibility of creating the graphene-based transistor [71]; this evolutionary research took the technology a step further by binding the transistor and other electronics on a single chip to build a full-fledged integrated circuit. In this work, the wafer-scaled graphene-based circuit has been reported, where graphene channel-based FET, interconnect, and inductor, were monolithically integrated on a single SiC wafer. The IC was operated as a broadband RF-mixer operative up to 10 GHz frequency range. The graphene-based system exhibited outstanding thermal stability between 300 and 400 K chip operating temperature. This integration gives an inference to a platform created by Meyer et al. in 2007 with suspended graphene [72]. The pioneering work done by IBM [3] created a perpetual impact on future research on graphene-based technology for integrated circuit applications.

    Heinz and co-workers modeled as well as fabricated the monolayer graphene resonator on SiO2 substrate to demonstrate that graphene can transduce the motion of profound resonant circuit with minimal damping [73], the comprehensive experimental arrangement of the resonator along with schematics, circuit diagram, and SEM image are portrayed in Figure 4. This study demonstrated the superiority of graphene over CNT, evaluating the reproducible electrical properties and a larger surface area for capturing the incoming mass flux. In this work, monolayer graphene flakes have been fabricated on locating on Si/SiO2. Subsequently, atomically precise bottom-up fabrication of GNRs on a gold substrate has been demonstrated by Fasel and his team [74], which formed GNRs. This method produces GNRs with well-controlled widths and various shapes. The controllable preparation of graphene has been pioneered by Müllen and co-workers [75,76] with bottom-up organic synthesis so that it can be produced on a large scale as an on-chip VLSI interconnect. On top of that, recently, Fischer and his team demonstrated an ingenious method in Science [77] to produce metallic GNRs based on the atomically precise bottom-up synthesis.

    Figure 4.  Experimental setup of graphene resonator on SiO2 substrate: (a) SEM image of monolayer graphene-based resonators, (b) SEM image of suspended graphene nanoribbon, (c) Circuit diagram of graphene-based resonator where lock-in amplifier detects the current through graphene, and (d) Schematic of suspended graphene (Reprinted (adapted) from Ref. [73] with permission of Springer Nature).

    Chen et al. [78] fabricated graphene-based interconnects with CMOS-based mixed-signal application, as shown in Figure 5. The graphene layers have been moved to the substrate with the help of the poly-methyl-methacrylate (PMMA) assisted method, which has been demonstrated through the progression flow chart of Figure 5a. Subsequently, graphene-based interconnects and Ti-based via have been developed with CMOS-compatible fabrication methodology. Cr/Au has been used as the contact metal of graphene-based interconnects, as shown in Figure 5b. An optical image has been captured demonstrating the graphene-based interconnects on the top layer of the ring oscillator array in CMOS based process (Figure 5c, d).

    Figure 5.  Process flow of on-chip graphene interconnect fabrication. (a) Ring oscillator circuit with graphene-based interconnect, (b) CVD-grown graphene transferred to substrate, (c) integration of graphene with CMOS circuitry of ring oscillator, (d) optical image of (a), (e) optical image of single graphene interconnect as a subset of (a). (Inset) AFM image of Si/SiO2 substrate-based graphene stripes (Reprinted from Ref. [78] with permission of IEEE).

    Han et al. [79] of the IBM research group reported the graphene-based on-chip RF receiver circuit, which has been integrated with a high-performance 3-stage RF amplifier, filter, and down-converter mixer module. This work revealed the use of the graphene-based system for high-end RF communication of 4.3 GHz bandwidth. Among a number of methods which have been reported to obtain intercalation doping in MLGNR, specifically chemical and electro-chemical intercalation method have been subjugated in practical implementation. Each of them is further subdivided into the categories depending upon the base of the intercalation species. As an example, when the chemical vapor transport method has been adopted for the chemical intercalation method, it is categorized as gas-phase intercalation; on the other hand, wet chemical methods are classified as liquid phase intercalation. Solid-state electrolyte intercalation and liquid electrolyte intercalation are the categories of electro-chemical intercalation doping. Intermittently molten metal (salt) intercalation doping (Li in few-layer graphene) has been reported in experimental reports [25,80], whereas Crommie et al. demonstrated intercalation of C60 into a single fold in graphene in [81]. Figure 6 provides a summary of the field of intercalation doping of 2D materials.

    Figure 6.  Summary of the intercalation of 2D materials, intercalation types, properties, and applications (Republished from Ref. [27] with permission of Royal Society of Chemistry).

    Recently, Jiang et al. demonstrated the CVD-grown and FeCl3- intercalation-doped MLGNR for superior interconnect application [30]. MLGNR layer was first transferred to the Si- substrate. Subsequently, the Ni/Au-based alloy contacts were built for the measurement of resistivity, as shown in the flow of Figure 7. FeCl3 intercalation doping was carried out at 633 K and ~1.4 atmospheric pressure for 10 h in an Argon atmosphere, as shown in Figure 7i, j [30].

    Figure 7.  (a) Summery of ML-graphene fabrication, (b) ML-graphene transferred on SiO2 substrate followed by oxygen plasma etching, (c) FeCl3 intercalation at 633 K, (d) Four-probe test structures, SEM images of (e) the test structure with patterned GNR, (f–h) GNRs. AFM image of (i) pristine ML-graphene and (j) FeCl3 intercalated ML-graphene shows 20% increase in thickness post-doping of 10 h; The Raman spectrum of (k, n) pristine one, (l, o) FeCl3-intercalated one and (m) FeCl3 (Reprinted (adapted) from Ref. [30] with permission of American Chemical Society).

    A detailed literature review on graphene-based nano-interconnects is presented in this article. At first, the modeling works on GNR interconnects with different contact models are discussed. Several analysis works are reviewed thoroughly for the applicability of GNR as a nano-interconnect for future generation VLSI circuits. Next, the GNR growth technologies are presented along with different process technology. As prospective mend of nano-interconnect application, intercalation doped GNR has been emphasized and discussed with its experimental and model development outline in this review. In a nutshell, this article provides a scientific survey in favor of pristine and intercalation doped layered graphene nanoribbons as a potential candidate among nano-interconnects of future integrated circuit-based applications.

    This work has been partially supported by SMDP-C2SD, DeitY, MCIT, India.

    The authors declare no conflict of interest.



    [1] J. M. Bové, Huanglongbing: a destructive, newly-emerging, century-old disease of citrus, J. Plant. Pathol., 88 (2006), 7-37.
    [2] R. A. Taylor, E. A. Mordecai, C. A. Gilligan, J. R. Rohr, L. R. Johnson, Mathematical models are a powerful method to understand and control the spread of Huanglongbing, Peer J., 19 (2016), 2642.
    [3] J. M. Bové, M. E. Rogers, Huanglongbing-control workshop: summary, Acta. Hort., 1065 (2015), 869-889.
    [4] A. J. Ayres, J. B. Jr, J. M. Bové, The experience with Huanglongbing management in Brazil, Acta. Hort., 1065 (2015), 55-61.
    [5] K. Chen, C. Li, Survey of citrus Huanglongbing research, J. Zhejiang Agr. Sci., 56 (2015), 1048-1050. (in Chinese)
    [6] C. Chiyaka, B. H. Singer, S. E. Halbert, J. G. M. Jr, A. H. C. van Bruggen, Modeling Huanglongbing transmission within a citrus tree, Proc. Natl. Acad. Sci. USA, 109 (2012), 12213-12218.
    [7] M. S. Chan, M. J. Jeger, An analytical model of plant virus disease dynamics with roguing and replanting, J. Appl. Ecol., 31 (1994), 413-427.
    [8] M. J. Jeger, L. V. Madden, F. van den Bosch, The effect of transmission route on plant virus epidemic development and disease control, J. Theor. Biol., 258 (2009), 198-207.
    [9] M. S. Sisterson, Effects of insect-vector preference for healthy or infected plants on pathogen spread: Insights from a model, J. Econ. Entomol., 101 (2008), 1-8.
    [10] X. Z. Meng, Z. Q. Li, The dynamics of plant disease models with continuous and impulsive cultural control strategies, J. Theor. Biol., 266 (2010), 29-40.
    [11] X. Z. Meng, Z. T. Song, L. S. Chen, A new mathematical model for optimal control strategies of integrated pest management, J. Biol. Syst., 15, (2007), 219-234.
    [12] W. C. Zhao, J. Li, T. Q. Zhang, X. Z. Meng, T. H. Zhang, Persistence and ergodicity of plant disease model with Markov conversion and impulsive toxicant input, Commun. Nonlinear Sci., 48 (2017), 70-84.
    [13] L. J. Xia, S. J. Gao, Q. Zou, J. P. Wang, Analysis of a nonautonomous plant disease model with latent period, Appl. Math. Comput., 223 (2013), 147-159.
    [14] S. J. Gao, L. J. Xia, J. L. Wang, Z. J. Zhang, Modelling the effects of cross protection control in plant disease with seasonality, Int. J. Biomath., 10 (2017), 1-24.
    [15] S. J. Gao, D. Yu, X. Z. Meng, F. M. Zhang, Global dynamics of a stage-structured Huanglongbing model with time delay, Chaos. Soliton. Fract., 117 (2018), 60-67.
    [16] X. S. Zhang, J. Holt, J. Colvin, Mathematical models of host plant infection by helper-dependent virus complexes: why are helper viruses always avirulent, Phytopathology, 90 (2000), 85-93.
    [17] A. U. Awan, M. Ozair, Stability analysis of pine wilt disease model by periodic use of insecticides, J. Biol. Dynam., 10 (2016), 506-524.
    [18] M. A. Khan, R. Khan, Y. Khan, S. Islam, A mathematical analysis of Pine Wilt disease with variable population size and optimal control strategies, Chaos Soliton Fract., 108 (2018), 205-217.
    [19] K. Jacobsen, J. Stupiansky, S. S. Pilyugin, Mathematical modeling of citrus groves infected by Huanglongbing, Math. Biosci. Eng., 10 (2013), 705-728.
    [20] X. S. Zhang, J. Holt, J. Colvin, Synergism between plant viruses: a mathematical analysis of the epidemiological implications, Plant Pathol., 50 (2001), 732-746.
    [21] K. S. Lee, D. Kim, Global dynamics of a pine wilt disease transmission model with nonlinear incidence rates, Appl. Math. Model, 37 (2013), 4561-4569.
    [22] M. Ozair, X. Y. Shi, T. Hussain, Control measures of pine wilt disease, Comp. Appl. Math., 35 (2016), 519-531.
    [23] M. A. Khan, K. Ali, E. Bonyah, K. O. Okosun, S. Islam, A. Khan, Mathematical modeling and stability analysis of Pine Wilt Disease with optimal control, Sci. Rep., 7 (2017), 3115.
    [24] M. A. Khan, R. Khan, Y. Khan, S. Islam, A mathematical analysis of Pine Wilt disease with variable population size and optimal control strategies, Chaos Soliton Fract., 108 (2018), 205-217.
    [25] J. A. Lee, S. E. Halbert, W. O. Dawson, C. J. Robertson, J. E. Keesling, B. H. Singer, Asymptomatic spread of Huanglongbing and implications for disease control, Proc. Natl. Acad. Sci. USA, 112 (2015), 7605-7610.
    [26] R. G. dA. Vilamiu, S. Ternes, G. A. Braga, F. F. Laranjeira, A model for Huanglongbing spread between citrus plants including delay times and human intervention, AIP Conf. Proc., 1479 (2012), 2315-2319.
    [27] Wikipedia, Available from: https://en.wikipedia.org/wiki/Citrus_greening-disease.
    [28] M. X. Deng, Forming process and basis and technological points of the theory emphasis on control citrus psylla for integrated control Huanglongbing, Chin. Agric. Sci. Bull., 25 (2009), 358-363.
    [29] O. Diekmann, J. A. P. Heesterbeek, J. A. J. Metz, On the definition and the computation of the basic reproduction ratio R0 in models for infectious diseases in heterogeneous populations, J. Math. Biol., 28 (1990), 365-382.
    [30] P. Dreessche, J. Watmough, Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission, Math. Biosci., 180 (2002), 29-48.
    [31] X. Q. Zhao, Dynamical Systems in Population Biology (Second Edition), Springer-Verlag, New York, 2017.
    [32] H. L. Smith, Monotone Dynamical System: An Introduction to the Theory of Competitive and Cooperative Systems, American Mathematical Society, Providence, 1995.
    [33] C. Castillo-Chavez, H. Thieme, Asymptotically autonomous epidemic models, Math. Popul. Dyn.: Anal. Heterog., 1 (1995), 33-50.
    [34] H. R. Thieme, Persistence under relaxed point-dissipativity with an application to an epidemic model, SIAM J. Math. Anal., 24 (1993), 407-435.
    [35] R. Gamkrelidze, L. S. Pontrjagin, V. G. Boltjanskij, The Mathematical Theory of Optimal Processes, Macmillan Company, 1964.
    [36] S. M. Blower, H. Dowlatabadi, Sensitivity and uncertainty analysis of complex models of disease transmission: An HIV model, as an example, Int. Stat. Rev., 62 (1994), 229-243.
    [37] J. C. Helton, J. D. Johnson, C. J. Sallaberry, C. B. Storlie, Survey of sampling-based methods for uncertainty and sensitivity analysis, Reliab. Eng. Syst. Saf., 91 (2006), 1175-1209.
    [38] S. Marino, I. B. Hogue, C. J. Ray, D. E. Kirschner, A methodology for performing global uncertainty and sensitivity analysis in systems biology, J. Theor. Biol., 254 (2008), 178-196.
    [39] M. A. Sanchez, S. M. Blower, Uncertainty and sensitivity analysis of the basic reproductive rate: tuberculosis as an example, Am. J. Epidemiol., 145 (1997), 1127-1137.
    [40] S. Lenhart, J. T. Workman, Optimal Control Applied to Biological Models, Chapman and Hall/CRC, Boca Raton, FL, 2007.
    [41] X. F. Yan, Y. Zou, Optimal and sub-optimal quarantine and isolation control in SARS epidemics, Math. Comput. Model., 47 (2008), 235-245.
    [42] F. B. Agusto, M. A. Khan, Optimal control strategies for dengue transmission in pakistan, Math. Biosci., 305 (2018), 102-121.
  • This article has been cited by:

    1. Shahnawaz Ayoub, Rajesh Kumar Rameshbhai Savaliya, Purushothaman G, Nitin Girdharwal, Mrinal Verma, Mohit Tiwari, 2023, Intelligent Optimization Algorithm based Data Storage and Management: An Empirical Study, 979-8-3503-4664-0, 1430, 10.1109/ICEARS56392.2023.10085038
    2. Shu Yuan, Yunfeng Hua, Jun Ye, 2023, A Highway Engineering Evaluation System Model Based on Blockchain and Trusted Sharing, 979-8-3503-1341-3, 1, 10.1109/EASCT59475.2023.10392648
    3. Albandari Alsumayt, Nahla El-Haggar, Lobna Amouri, Zeyad M. Alfawaer, Sumayh S. Aljameel, Smart Flood Detection with AI and Blockchain Integration in Saudi Arabia Using Drones, 2023, 23, 1424-8220, 5148, 10.3390/s23115148
    4. Princy Diwan, Ramgopal Kashyap, Brijesh Khandelwal, Blockchain assisted encryption scheme for intellectual share estimation using medical research data, 2024, 36, 1532-0626, 10.1002/cpe.7896
    5. Shuangquan Chen, Jinguo Li, Kai Zhang, Aoran Di, Mengli Lu, Privacy-Preserving Breast Cancer Prediction Based on Logistic Regression, 2024, 67, 0010-4620, 2667, 10.1093/comjnl/bxae035
    6. Jiayuan Zhang, Rongxin Guo, Yifan Shi, Wanting Tang, An anti-impersonation attack electronic health record sharing scheme based on proxy re-encryption and blockchain, 2024, 21, 1551-0018, 6167, 10.3934/mbe.2024271
    7. Soumyaranjan Panda, Rajni Dubey, Biswajit Jena, Vikas Pareek, Lung-Wen Tsai, Sanjay Saxena, Federated Learning Frameworks in Smart e-Healthcare: A Systematic Literature Review with Bias Evaluation, 2025, 15684946, 112747, 10.1016/j.asoc.2025.112747
  • Reader Comments
  • © 2020 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(5210) PDF downloads(494) Cited by(17)

Figures and Tables

Figures(11)  /  Tables(2)

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog