Loading [MathJax]/jax/output/SVG/jax.js
Theory article

Dynamic analysis of a delayed differential equation for Tropidothorax elegans pests

  • Received: 23 August 2023 Revised: 11 October 2023 Accepted: 15 October 2023 Published: 30 October 2023
  • In this paper, we establish an infectious disease model of Tropidothorax elegans to study the impact of them on plants. Our model involves the time delay for Tropidothorax elegans to hatch eggs, which is influenced by temperature. Second, we theoretically analyze the existence and the stability of the equilibrium and the normal form near the Hopf bifurcating critical point. Next, we choose three groups of parameters for numerical simulations to verify theoretical analysis of our model. Then, based on numerical simulations, we give bioanalysis which are consistent with the patterns of Tropidothorax elegans pests, such as dying off in large numbers of adults during the winter and one or two generations a year.

    Citation: Tingru Yang, Yuting Ding. Dynamic analysis of a delayed differential equation for Tropidothorax elegans pests[J]. Electronic Research Archive, 2023, 31(11): 6947-6963. doi: 10.3934/era.2023352

    Related Papers:

    [1] Rina Su . Dynamic analysis for a class of hydrological model with time delay under fire disturbance. Electronic Research Archive, 2022, 30(9): 3290-3319. doi: 10.3934/era.2022167
    [2] Miao Peng, Rui Lin, Zhengdi Zhang, Lei Huang . The dynamics of a delayed predator-prey model with square root functional response and stage structure. Electronic Research Archive, 2024, 32(5): 3275-3298. doi: 10.3934/era.2024150
    [3] Mengting Sui, Yanfei Du . Bifurcations, stability switches and chaos in a diffusive predator-prey model with fear response delay. Electronic Research Archive, 2023, 31(9): 5124-5150. doi: 10.3934/era.2023262
    [4] Xinrui Yan, Yuan Tian, Kaibiao Sun . Effects of additional food availability and pulse control on the dynamics of a Holling-($ p $+1) type pest-natural enemy model. Electronic Research Archive, 2023, 31(10): 6454-6480. doi: 10.3934/era.2023327
    [5] Qinghua Zhu, Meng Li, Fang Han . Hopf bifurcation control of the ML neuron model with Hc bifurcation type. Electronic Research Archive, 2022, 30(2): 615-632. doi: 10.3934/era.2022032
    [6] Yujia Xiang, Yuqi Jiao, Xin Wang, Ruizhi Yang . Dynamics of a delayed diffusive predator-prey model with Allee effect and nonlocal competition in prey and hunting cooperation in predator. Electronic Research Archive, 2023, 31(4): 2120-2138. doi: 10.3934/era.2023109
    [7] Fengrong Zhang, Ruining Chen . Spatiotemporal patterns of a delayed diffusive prey-predator model with prey-taxis. Electronic Research Archive, 2024, 32(7): 4723-4740. doi: 10.3934/era.2024215
    [8] San-Xing Wu, Xin-You Meng . Hopf bifurcation analysis of a multiple delays stage-structure predator-prey model with refuge and cooperation. Electronic Research Archive, 2025, 33(2): 995-1036. doi: 10.3934/era.2025045
    [9] Yuan Xue, Jinli Xu, Yuting Ding . Dynamics analysis of a diffusional immunosuppressive infection model with Beddington-DeAngelis functional response. Electronic Research Archive, 2023, 31(10): 6071-6088. doi: 10.3934/era.2023309
    [10] Xueyong Zhou, Xiangyun Shi . Stability analysis and backward bifurcation on an SEIQR epidemic model with nonlinear innate immunity. Electronic Research Archive, 2022, 30(9): 3481-3508. doi: 10.3934/era.2022178
  • In this paper, we establish an infectious disease model of Tropidothorax elegans to study the impact of them on plants. Our model involves the time delay for Tropidothorax elegans to hatch eggs, which is influenced by temperature. Second, we theoretically analyze the existence and the stability of the equilibrium and the normal form near the Hopf bifurcating critical point. Next, we choose three groups of parameters for numerical simulations to verify theoretical analysis of our model. Then, based on numerical simulations, we give bioanalysis which are consistent with the patterns of Tropidothorax elegans pests, such as dying off in large numbers of adults during the winter and one or two generations a year.



    In the 21st century, the trend of global warming is becoming more and more pronounced, and industrial gases such as carbon dioxide are a major culprit in global warming. Forests play a key role in carbon sequestration to mitigate the process of global warming, and the capacity of forests to sequester carbon is assessed in [1,2]. To achieve the double carbon goal, one aspect is to protect forests from forest diseases and pests. Therefore, many scholars are studying ways to protect the forest to avoid the damage of forest diseases and pests. The key is to master the patterns of forest pest and disease outbreaks.

    To find the patterns of forest pest and disease, studying the effects of factors such as meteorological factors and human intervention on the outbreaks of forest diseases and pests is significant. Various models of forest diseases and pests were studied [3,4,5,6]. Many scholars provided theoretical knowledge for modeling analysis [7,8,9]. In the work of Yao et al. [4], the relation of environment meteorological conditions with the area of forest diseases and pest was clearly studied by using the method of abnormal index analysis. Kiewra et al. [5] established the weather research and forecasting model to show the influence of different meteorological variables. Likewise, Wang et al. [6] used a MaxEnt model and incorporated annual meteorological and human activity factors to study the outbreaks of pine wilt disease. Obviously, meteorological condition is an important factor in outbreaks of forest pest and disease.

    In addition, human intervention is also a key factor affecting forest pests and diseases. In References [10,11,12,13], scholars studied many managements to control forest pest and disease, such as using the GIS-based system which is a geographic information system to project the outbreaks of pest, using the skill of genetically engineered baculoviruses for pest, using chemical control and biological control and so on. These are efforts made by people to control pests.

    In this paper, we focus on pests caused by the Tropidothorax elegans. This is a class of Lygaeidae with an extremely wide distribution. The appearance and biological habits of Lygaeidae are studied [14,15]. All pests of Lygaeidae are a class of stinging pests and harm plants through their mouthparts. They reproduce one or two generations per year, get through the winter as adults and lay eggs in the following year. During the period of damage, both Tropidothorax elegans nymphs and adults suck the leaves, stem and roots of plants. The sucked parts of the leaves will cause the plant to be infected with sooty blotch under the action of some spores, forming a black spot on the leaves, which will greatly affect the photosynthesis of the plant and prevent the plant from absorbing nutrients. This mode of transmission is similar to infectious diseases, and many scholars studied the transmission patterns of many types of infectious diseases and predicted the development of the infectious disease by establishing the model of infectious diseases [16,17,18]. Some other scholars gave different measures based on the model analysis to block the spread of the virus [19,20,21,22].

    After gathering preliminary knowledge of Tropidothorax elegans and the current research status, we establish a model about Tropidothorax elegans and arrange the rest of the context as follows. In Section 2, according to the effect of Tropidothorax elegans on plants, we establish an infectious disease model on Tropidothorax elegans. In Section 3, we discuss the stability of the equilibrium point of the model and the existence of the Hopf bifurcation. In Section 4, we calculate the normal form of the Hopf bifurcation. In Section 5, we use numerical simulations to verify the correctness of bioanalysis. Finally, we draw the corresponding conclusions in Section 6.

    Based on the stinging characteristics of Tropidothorax elegans, we assume that a susceptible plant becomes an infected plant when it is sucked by Tropidothorax elegans. Therefore, we divide all plants into two groups. One group is susceptible plants (S(t), 104 plants), which have not been harmed by Tropidothorax elegans. The other group is infected plants (I(t), 104 plants), i.e., from susceptible plants that are sucked by Tropidothorax elegans (S(t)I(t)). The infected plants have two types of changes: somes become healthy after treatment in time, and the others will die without timely treatment. Besides, the recovered plants can be reinfected again. We denote the number of Tropidothorax elegans as X(t) (106 pests). Based on above information, the relationship between the variables can be given in Figure 1.

    Figure 1.  SIX Model diagram.

    In addition, the hatching time of Tropidothorax elegans is also an important factor. We make the hatching time of Tropidothorax elegans as the time delay in our study. Subject to resource constraints, we view the reproductive capacity of Tropidothorax elegans as a decreasing function related to the number of prehatch Tropidothorax elegans. Thus, we get α(t)=α1(1X(tτ)N), where N is the environmental capacity of Tropidothorax elegans and τ is the time delay of the pest hatching period (year). We establish the model as follows:

    {dS(t)dt=α2γX(t)S(t)+θI(t),dI(t)dt=γX(t)S(t)θI(t)β2I(t),dX(t)dt=α1X(t)(1X(tτ)N)β1X(t), (2.1)

    where α1,α2,β1,β3,γ,θ,N are nonnegative parameters. The specific descriptions are given in Table 1. In this table, all variables and parameters are nonnegative.

    Table 1.  The descriptions of variables and parameters.
    Symbol Descriptions
    S Number of susceptible plants
    I Number of infected plants
    X Number of Tropidothorax elegans
    α1 Reproductive ability of Tropidothorax elegans
    α2 Replenishment rate of plants
    β1 Mortality rate of Tropidothorax elegans
    β2 Mortality rate of infective plants
    γ Transition rate from S to I
    θ Transition rate from I to S
    N Environmental capacity of Tropidothorax elegans in the area

     | Show Table
    DownLoad: CSV

    In this section, the system (2.1) is considered. When α2=0, we calculate a nonisolated boundary solution E0=(S0,0,0), where S0 is an arbitrary value. It means that there is no local infestation of Tropidothorax elegans. Therefore, there are no infected plants. However, we focus on the presence of Tropidothorax elegans. Thus, We consider the following assumption:

    (A1)α1β1>0.

    When (A1) holds, we can obtain the unique and positive equilibrium of system (2.1):

     E=(S,I,X), (3.1)

    where S=α1α2(β2+θ)γN(α1β1)β2, I=α2β2, X=(α1β1)Nα1.

    Since α1 represents the reproductive capacity of Tropidothorax elegans and β1 represents the mortality rate of Tropidothorax elegans, (A1) implies that the reproductive capacity of insects is much greater than the mortality rate of insects in fact.

    In this section, we consider the nonnegativity of system (2.1).

    Theorem 3.1. If S(θ)0, I(θ)0, X(θ)0(θ[τ,0]), the solution S(t),I(t),X(t) of system (2.1) with τ0 is nonnegative when t>0.

    Proof. First, we prove X(t)0 when t0 under the initial condition of system (2.1). We assume that X(t)0 is not always nonnegative for t0 and make t1 be the time that X(t1)=0 and X(t1)<0 hold.

    According to the third equation of system (2.1), we can obtain X(t1)=0. The two conclusions are contradictory. Therefore, X(t)0 when t>0.

    Then, we make t2 be the time that I(t2)=0 and I(t2)<0 hold. According to the second equation of system (2.1), we can obtain I(t2)=γX(t2)S(t2). If I(t2)<0 holds, then S(t2)<0 and t3<t2. Let S(t3)=0, S(t3)=0.

    According to the first equation of system (2.1), we can obtain S(t3)=α2+θI(t3). Because I(t3)>0, we obtain S(t3)>0. The two conclusions are contradictory. So S(t3)0 when t>0. Thus, the solution to system (2.1) is nonnegative when t>0.

    In this section, we consider the stability and the existence of Hopf bifurcation of E. We can obtain the characteristic equation of system (2.1) at E as follows:

     [λ+(α1β1)eλτ][λ2+(θ+β2+γN(α1β1)α1)λ+γ(α1β1)β2α1]=0. (3.2)

    When τ=0, it becomes:

     [λ+(α1β1)][λ2+(θ+β2+γN(α1β1)α1)λ+γ(α1β1)β2α1]=0. (3.3)

    Equation (3.3) has three roots: let λ1, λ2 be the roots of the part of quadratic equation in Eq (3.3), λ1+λ2=(θ+β2+γN(α1β1)α1), λ1λ2=γ(α1β1)β2α1 and λ3=β1α1. (A1) holds and other parameters are positive. Thus, Eq (3.3) has three negative roots and E is locally asymptotically stable for τ=0.

    Next, we consider the existence of bifurcation periodic solutions near the equilibrium E for τ>0. Let λ=iω(ω>0) be a root of characteristic Eq (3.2). We only consider the equation as follows:

     λ+(α1β1)eλτ=0. (3.4)

    Then, substituting λ=iω(ω>0) into Eq (3.4) and separating the real and imaginary parts, we obtain:

    cos(ωτ)=0,sin(ωτ)=ωα1β1. (3.5)

    Adding the square of the two equations in Eq (3.5), we obtain ω=α1β1, and the expression of τ:

     τ(j)=1(α1β1)(π2+2jπ),j=0,1,2,. (3.6)

    When τ=τ(j)(j=0,1,2,), characteristic Eq (3.4) have a pair of pure imaginary roots λ=±iω. Calculating the transversality conditions, we obtain:

    Re(dλdτ)1τ=τ(j)=1(α1β1)2>0,j=0,1,2,.

    Theorem 3.2. When (A1) holds, the unique equilibrium E exists and is positive, and system (2.1) undergoes Hopf bifurcation when τ=τ(j)(j=0,1,2,) near E, where τ(j) is given by Eq (3.6). Then, when τ[0,τ(0)), the equilibrium E is locally asymptotically stable and unstable when τ>τ(0).

    In this section, we derive the normal form of Hopf bifurcation for the system (2.1) by using the multiple time scales method. To reflect the actual situation, we focus on the delay in hatching and study the patterns of Tropidothorax elegans pests. We consider the time-delay τ as a bifurcation parameter. Let τ=τc+εμ, where τc=τ(j)(j=0,1,2,) is the Hopf bifurcation critical value, μ is perturbation parameter and ε is dimensionless scale parameter. When τ=τc, the Eq (3.4) has eigenvalue λ=iω(ω>0), at which system (2.1) undergoes a Hopf bifurcation at equilibrium E.

    To translate the equilibrium point E=(S,I,X) to the origin and simplify the expression, we write ˜S as S, ˜I and ˜X are the same. At the same time, to normalize the delay, we let ttτ rescale the time. Then, we obtain following equations:

    {˙S=τ[γ(XS+XS+SX)+θI],˙I=τ[γ(XS+XS+SX)(θ+β2)I],˙X=τ[α1Xα1N(XX(t1)+XXX(t1)X)]. (4.1)

    We write Z(t)=(S(t),I(t),X(t))T,Z(t1)=(S(t1),I(t1),X(t1))T. The Eq (4.1) can be written as:

     ˙Z=τAZ+τBZ(t1)+τF(Z,Z(t1)), (4.2)

    where F(Z,Z(t1))=(γSX,γSX,α1NXX(t1))T,

    A=(γN(α1β1)α1θα1α2(θ+β2)N(α1β1)β2(γNα1β1)α1(θ+β2)α1α2(θ+β2)N(α1β1)β2000),B=(00000000β1α1).

    Then, we make h the eigenvector corresponding to eigenvalue λ=iωτ of the linear section of Eq (4.2) and h the eigenvector corresponding to eigenvalue λ=iωτ of adjoint matrix of the linear section of Eq (4.2). Then, h and h should satisfy h,h=¯hTh=1. According to the above conditions, we obtain:

    h=(h1h2h3)=(VmωmUV+Cm),h=(001). (4.3)

    where m=α12α2(β2+θ)N(α1β1)β2, U=γ(α1β1)N, V=β2i+ω, C=(iω+θ+β2)α1.

    The solution of Eq (4.2) is assumed as follows:

     Z(t)=Z(T0,T1,T2,)=+k=1εkZk(T0,T1,T2,), (4.4)

    where

    Z(t)=Z(T0,T1,T2,)=(S(T0,T1,T2,),I(T0,T1,T2,),X(T0,T1,T2,))T,
    Zk(t)=Zk(T0,T1,T2,)=(Sk(T0,T1,T2,),Ik(T0,T1,T2,),Xk(T0,T1,T2,))T.

    The derivative, with regard to t, is transformed into

    ddt=T0+εT1+ε2T2+=D0+εD1+ε2D2+,

    where Di is differential operator, and

    Di=Ti,i=0,1,2,.

    From Eq (4.4), we can obtain

     ˙Z=εD0Z1+ε2D1Z1+ε3D2Z1+ε2D1Z2+ε3D0Z3+. (4.5)

    Then, using the Taylor expansion to expand X(T01,ε(T01),ε2(T01),) at X(T01,T1,T2,), we obtain

     X(t1)=εX1,τc+ε2X2,τc+ε3X3,τcε2D1X1,τcε3D2X1,τcε3D1X2,τc+, (4.6)

    where Xj,τc=Xj(T01,T1,T2,),j=1,2,. Combing in Eqs (4.4)–(4.6) and separating the coefficients of ε, ε2 and ε3 on both sides of the equation, the expressions of the coefficient before ε, ε2 and ε3 are shown from (4.7)–(4.9):

    {D0S1+τc(γSX1+γXS1θI1)=0,D0I1τc[γSX1+γXS1(θ+β2)I1]=0,D0X1+τc[α1X1α1XN(X1+X1,τc)β1X1]=0. (4.7)
    {D0S2+τc(γSX2+γXS2θI2)=D1S1γτcX1S1μ(γSX1+γXS1θI1),D0I2τc[γSX2+γXS2(θ+β2)I2]=D1I1+γτcX1S1+μ[γSX1+γXS1(θ+β2)I1],D0X2τc[α1X2α1XN(X2+X2,τc)β1X2]=D1X1+α1τcN(XD1X1,τcX1X1,τc)+μ[α1X1α1XN(X1+X1,τc)β1X1]. (4.8)
    {D0S3+τc(γSX3+γXS3θI3)=D2S1D1S2γτc(X1S2+X2S1)μγ(SX2+XS2+X1S1)+μθI2,D0I3τc[γSX3+γXS3(θ+β2)I3]=D2I1D1I2+γτc(X1S2+X2S1)+μγ(SX2+XS2+X1S1)+μ(θ+β2)I2,D0X3τc[α1X3α1XN(X3+X3,τc)β1X3]=D2X1D1X2α1τcN[X1X2,τcX1D1X1,τc+X2X1,τcX(D2X1,τc+D1X2,τc)]+μ[α1X2α1N(X1X1,τc+X(X2+X2,τcD1X1,τc))β1X2]. (4.9)

    The solution of Eq (4.7) can be assumed as:

     Z1=GheiωτcT0+ˉGˉheiωτcT0, (4.10)

    where h is given by Eq (4.3). We substitute Eq (4.10) into the right side of the Eq (4.7) and write the coefficient of eiωτcT0 as m1. According to the solvability condition h,m1=0, the expression of GT1 is obtained as follows:

     GT1=MμG, (4.11)

    where

    M=a[α1XN(1+eiωτc)(α1β1)],

    with

    a=(τcα1XNeiωτc1)1.

    Since μ is a perturbation term, its effect on higher order is smaller. Thus, we only consider its effect on the linear part and ignore its effect on higher order. We assume that the solution of Eq (4.8) is as follows:

     S2=g1e2iωτT0G2+¯g1e2iωτT0¯G2+l1G¯G, I2=g2e2iωτT0G2+¯g2e2iωτT0¯G2+l2G¯G, X2=g3e2iωτT0G2+¯g3e2iωτT0¯G2+l3G¯G, (4.12)

    where

    g1=k1k2,g2=2iω2iω+β2k1k2,g3=α1h23eiωτcN(2iωα1+β1)α1X(1+e2iωτc),l1=2θk3+(γSk4k3)β2γXβ2,l2=2k3β2,l3=k4,

    with

    k1=γSg3+γh1h3,k2=(2iω+γX+2θωi2iω+β2)1,
    k3=γ(h3¯h1+¯h3h1),k4=α1h3¯h3(eiωτc+eiωτc)(α1β1)N2α1X,

    and h1,h2,h3 is given by Eq (4.3).

    Then we substitute Eq (4.12) into the right side of Eq (4.9), and the coefficient vector of eiωτcT0 is denoted by m2. According to the solvability condition h,m2=0, the expression of GT2 can be obtained as follows:

     GT2=α1τcHG2¯G, (4.13)

    where

    H=d[h3l3(1+eiωτc)+¯h3g3(e2iωτc+eiωτc)],

    with

    d=[h3(τcα1XeiωτcN)]1,

    and h3,g3,l3 are given by Eqs (4.3) and (4.12).

    Let GG/ε. Then, the deduced third-order normal form of Hopf bifurcation of system (2.1) is:

     ˙G=MμG+α1τcHG2¯G, (4.14)

    where M is given in Eq (4.11) and H is given in Eq (4.13).

    Substituting G=reiθ into Eq (4.14), the following normal form of Hopf bifurcation in polar coordinates is obtained:

    {˙r=Re(M)μr+Re(H)α1τcr3,˙θ=Im(M)μ+Im(H)α1τcr2. (4.15)

    Based on the normal form of Hopf bifurcation in polar coordinates, we only need to consider the first equation in system (4.15). Thus, the following theorem holds:

    Theorem 4.1. When Re(M)μRe(H)α1τc<0, the system (4.15) exits a semitrivial fixed point r=Re(M)μRe(H)α1τc, and the system (2.1) has a periodic solution.

    1) If Re(M)μ<0, the periodic solution reduced on the center manifold is unstable.

    2) If Re(M)μ>0, the periodic solution reduced on the center manifold is stable.

    In this section, we choose three groups of parameters to simulate the effects of temperature changes and artificial interventions on the model. Then, we give the corresponding biological explanations based on the simulations. Finally, we discuss what humans should do to reduce the harm caused by Tropidothorax elegans.

    First, we choose the first group of parameters to describe the situation that Tropidothorax elegans grow freely. Second, considering the context of global warming, we choose the second group of parameters to simulate the effects of global warming on Tropidothorax elegans pests. At higher latitudes, the increase in temperature will lead to a decrease in the number of adults that freeze to death in winter, and the best temperature for the survival of Tropidothorax elegans can be reached more easily, so the incubation time can be shortened [23]. At lower latitudes, increased temperatures will lead to a decrease in the mortality rate of adults, so there will be more Tropidothorax elegans laying eggs in next year. Thus, globally higher temperature leads to lower mortality and higher reproduction of Tropidothorax elegans around the world. Finally, we take the third group of parameters to simulate the effect of timely artificial control on Tropidothorax elegans pests. Timely artificial control can increase the mortality rate of Tropidothorax elegans [24]. Thus, increasing the recovery rate of infected plants and reducing the mortality rate of infected plants.

    In addition, α2 is the replenishment rate of plants, β2 is the mortality rate of infective plants and θ is the transition rate from I(t) to S(t). In fact, because people find plants dead and promptly replenish them with the same number of plants, we consider the replenishment rate to be the same as the mortality rate of plants. The infected plants have two types of changes: some become healthy after treatment in time and the others will die without timely treatment. Thus, we consider α2+θ=1.

    Based on the above analyses, three reasonable groups of parameters are taken to describe the above three situations as follows:

    1)N=4.5,α1=15,α2=0.05,θ=0.95,β1=0.8,β2=0.05,γ=0.29;2)N=4.5,α1=20,α2=0.10,θ=0.90,β1=0.6,β2=0.02,γ=0.29;3)N=4.5,α1=15,α2=0.02,θ=0.98,β1=0.9,β2=0.02,γ=0.29;

    First, we show the simulated results under the first group of parameters:

    N=4.5,α1=15,α2=0.05,θ=0.95,β1=0.8,β2=0.05,γ=0.29.

    We calculate the equilibrium E=(S,I,X)=(0.815195239,1,0.94) and critical τc=0.11135 of system (2.1). Thus, for the initial functions [S(θ),I(θ),X(θ)]=[0.81,0.99,0.93](θ[τc,0]), we choose τ=0, τ=0.08<τc and τ=0.18>τc for simulations. Based on Theorem 3.2, the equilibrium E is locally asymptotically stable when τ=0 and τ=0.08<τc. The result of simulation with τ=0.08 is shown in Figure 2, which is consistent with theoretical analyses.

    Figure 2.  Equilibrium E of the system (2.1) is locally asymptotically stable when τ=0.08.

    In Figure 2, S,I,X become stable after a short period of variations. In the early stage, the number of Tropidothorax elegans increases rapidly, and peak after a period of time. Since then, the number of Tropidothorax elegans begins to decline. We suspect this is because limited resources act as a constraint on their growth. After a short period of intraspecific competition, the number of Tropidothorax elegans becomes relatively stable. This also means that when the incubation time of the Tropidothorax elegans is small enough, the Tropidothorax elegans will remain locally forever. This may happen at low latitude regions where temperature is high enough.

    For τ=0.18>τc, we calculate the normal form of Hopf bifurcation and obtain Re(M)>0, Re(H)<0 from Eq (4.15). According to the Theorem 4.1, the periodic solution is stable when μ>0. The solution is shown in Figure 3.

    Figure 3.  When τ=0.18, there exists a stable periodic solution near equilibrium E.

    Clearly, S,I,X all form the stable periodic solutions in Figure 3. This means that this pest has become an endemic disease and is difficult to eradicate with τ=0.18. In Figure 3(c), the infestation period of Tropidothorax elegans is about one month and outbreaks once a year which is consistent with the habits of the Tropidothorax elegans (http://museum.ioz.ac.cn/topic_detail.aspx?id=7358). In addition, the number of Tropidothorax elegans tends to zero for a period of time. Then, it gradually increases again, repeating the trend of the first cycle. The biological meaning is that Tropidothorax elegans spend the winter as adults, the number of Tropidothorax elegans decreases considerably due to the low temperature. However, since Tropidothorax elegans go through the winter everywhere, such as near plants, on the backs of plant leaves, under rocks or in soil burrows, it is difficult to kill. Therefore, the number of the Tropidothorax elegans will increase again in next year because of their strong reproductive ability, repeating the experience of the previous year. Besides, the external environment is stable, so the environmental capacity rarely changes and the peak value remains constant.

    For the group of parameters (2):

    N=4.5,α1=20,α2=0.10,θ=0.90,β1=0.6,β2=0.02,γ=0.29.

    We calculate the equilibrium E=(S,I,X)=(0.789983015,1,0.97) and τc=0.08093. We know that E is locally asymptotically stable for τ(0,τc) and unstable for τ>τc. Similar to the solution of the first group of parameters when τ=0.18, the periodic solution of the second group of parameters with τ=0.12>τc is also stable when μ>0. We choose the initial functions [S(θ),I(θ),X(θ)]=[0.78,0.99,0.96](θ[τc,0]) for simulation and the result of simulation is shown in Figure 4.

    Figure 4.  When τ=0.12, there exists a stable periodic solution near equilibrium E.

    The overall trend is the same as Figure 3, the actual situation of these parameters is: under the background of global warming, temperature has increased. This leads to a lower natural mortality rate and a higher reproduction capacity of Tropidothorax elegans, which are more suitable for their survival. The time of the number of Tropidothorax elegans tends to zero is significantly shorter in Figure 4(c). That means that the frequency of infestation is significantly higher and the damage to the forest is greater. Thus, the increase of temperature will lead to the increase of damage frequency and damage degree of Tropidothorax elegans.

    Finally, we simulate the third group of parameters:

    N=4.5,α1=15,α2=0.02,θ=0.98,β1=0.9,β2=0.02,γ=0.29.

    We calculate the equilibrium E=(S,I,X)=(0.81520,1,0.94) and τc=0.11135 of system (2.1). Compared with the first group of parameters, we choose τ=0.18>τc and the initial functions [S(θ),I(θ),X(θ)]=[0.81,0.99,0.93](θ[τc,0]) for simulation. The result of simulation is shown in Figure 5.

    Figure 5.  When τ=0.18, there exists a stable periodic solution near equilibrium E.

    The Figure 5 corresponds to the actual situation: after people found this pest, people will kill pests in time. Obviously, the more timely the manual control, the longer it takes the number of Tropidothorax elegans to tend to zero, which shows that the manual control has suppressed the harm of the Tropidothorax elegans to the plants. When we choose the prevention methods, we should choose the optimal control method according to the different characteristics of the pest outbreak respectively.

    Based on the research [11,12,27], the advantage of biological control is that it is highly selective and safer to the environment. However, its killing effect is slow and the economic cost is high. Compared with it, the killing effect of chemical control is faster and the economic cost is low, but it will bring more environmental pollution and other hidden problems. For example, when chemical control is carried out, pests and their natural enemies are killed in large quantities, at which time the effectiveness of biological control is greatly reduced, but it does achieve the purpose of eliminating pests. With time going by, the remaining chemical reagents are reduced, and the ability to kill pests is weakened. Then, biological control is relied on to kill pests [28]. Therefore, in the case of sudden outbreak of Tropidothorax elegans, we can use chemical control to prevent and control the bug urgently, and in the later stage of pest management, we can use biological control technology to reduce the pollution.

    In addition, we compare the curves above three groups together as shown in Figure 6:

    Figure 6.  Comparison of three sets of parameters.

    We can see the similarities of these curves:

    1) The trend of the change of the number of infected plants and the number of susceptible plants is opposite, this is because that the total number of plants is a constant. which is consistent with the actual situation.

    2) In the stage of increasing number of Tropidothorax elegans, after the number of Tropidothorax elegans peaking, the number of susceptible and the number of infected plants peak. There is a delay between the time which the number of Tropidothorax elegans peaks with the time which the number of susceptible plants and infected plants peak. We analyze that this is because the replenishment rate of the plants is smaller than the reproductive capacity of Tropidothorax elegans and the mortality rate of the plants is also smaller than the mortality rate of Tropidothorax elegans, so the change of the number of plants is smaller and slower, while the change of Tropidothorax elegans will be more flexible.

    3) In the first and third groups of parameters, we choose different parameters and the same τ=0.18. In the Figure 6, we find that the two curves almost coincide. However, the time-delay of the second group of parameters is different from other groups, and the curves are different. This means that the effect of parameters variations on the patterns of Tropidothorax elegans pests is less than the effect of temperature on the patterns of Tropidothorax elegans pests. Therefore, we believe that the speed of reproduction determines the outbreak of Tropidothorax elegans.

    Some differences in the Figure 6:

    1) When the time-delay is the same, we compare the curves for the first group of parameters with the curves for the third group of parameters. The peak of the curves of the third group is lower than the peak of the curves of the first group in Figure 6(a), (b). In fact, this represents a reduction in the number of mortality plants through manual control during the outbreak. In addition, we compare the curves of the first and third group of parameters in Figure 6(c), we find that the breakout period of the curve of the artificial intervention is later than the breakout period of the first set of parameters. This means that artificial intervention can effectively postpone the outbreak of the pest. Therefore, we should do a good job of patrolling in the forest in time, so that the insect pests can be detected and treated in time to ensure the safety of the plants.

    2) When the time-delay is different, we compare the first group of parameters with τ=0.18 and the second group of parameters with τ=0.12 in Figure 6, and find that the curves of the second group of parameters change more fastly. As the temperature rises, the incubation time of Tropidothorax elegans will be shortened and the number of infected plants and susceptible plants will also change rapidly. Due to resource constraints and more intense intraspecific competition, the maximum number of Tropidothorax elegans in the region has decreased.

    In this paper, considering the characteristics of Tropidothorax elegans to sting plants, we have constructed a SIX model with a time delay for Tropidothorax elegans to hatch eggs, which is influenced by temperature and manual control. We have studied the stability of the equilibrium and the existence of Hopf bifurcation. Then, we have analyzed the stability and bifurcating direction of the Hopf bifurcating periodic solution by calculating the normal form with the multiple time scales method.

    Based on three groups of parameters, we have conducted numerical simulations. First, we have simulated the change of temperature through the first and the second groups of parameters to verify that Tropidothorax elegans pests become more frequent with higher temperature and that the laws of outbreaks is once or twice a year, which is consistent with actual situation. Second, we have simulated the effects of human intervention through the first and third groups of parameters and have compared them to conclude that human intervention can delay the outbreak of Tropidothorax elegans pests and reduce the number of Tropidothorax elegans at the peak of the outbreak. Further more, based on the analyzed characteristics of Tropidothorax elegans, we could choose different control methods and apply the right medicine to the right situation. In the future, we will further break down the effects of different control methods on this pest, focusing on two methods: pesticide control and natural enemy control, so that we could remove pests more efficiently and protect the environment at the same time.

    The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.

    This study was funded by Fundamental Research Funds for the Central Universities of China (No. 2572022DJ06) and College Students Innovations Special Project funded by Northeast Forestry University of China (No. DC-2023179).

    The authors declare that they have no competing interests.



    [1] W. Yang, Z. Min, M. X. Yang, J. Yan, Exploration of the implementation of carbon neutralization in the field of natural resources under the background of sustainable development-an overview, Int. J. Environ. Res. Public Health, 19 (2022), 14109. https://doi.org/10.3390/ijerph192114109 doi: 10.3390/ijerph192114109
    [2] H. Chen, C. C. Tie, Y. Teng, Z. Geng, J. M. Zhu, Research on forest carbon sequestration assessment and forest management plan based on TOPSIS method and grey correlation analysis, Int. J. Environ. Pollut., 70 (2021), 130–147. https://doi.org/10.1504/IJEP.2021.132005 doi: 10.1504/IJEP.2021.132005
    [3] B. Yang, W. Q. Cao, C. M. Tian, Visual analysis of impact factors of forest pests and diseases, J. Visualization, 22 (2019), 1257–1280. https://doi.org/10.1007/s12650-019-00593-2 doi: 10.1007/s12650-019-00593-2
    [4] X. H. Yao, Y. P. Xu, B. S. Yuan, L. Y. Han, X. L. Yao, H. H. Han, Study of meteorological assessment techniques of forest diseases and pests in Xiaolong mountain (in Chinese), J. Northwest A, F Univ. (Nat. Sci. Ed.), 10 (2012), 64–68. https://doi.org/10.13207/j.cnki.jnwafu.2012.10.018 doi: 10.13207/j.cnki.jnwafu.2012.10.018
    [5] D. Kiewra, M. Kryza, M. Szymanowski, Influence of selected meteorological variables on the questing activity of Ixodes ricinus ticks in Lower Silesia, SW Poland, J. Vector Ecol., 39 (2014), 138–145. https://doi.org/10.1111/j.1948-7134.2014.12080.x doi: 10.1111/j.1948-7134.2014.12080.x
    [6] W. Wang, W. T. Peng, X. Y. Liu, G. He, Y. L. Cai, Spatiotemporal dynamics and factors driving the distributions of pine wilt disease-damaged forests in China, Forests, 13 (2022), 261. https://doi.org/10.3390/f13020261 doi: 10.3390/f13020261
    [7] R. Castro, S. Lanucara, V. Piccione, G. Pioggia, G. Modica, M. A. Ragusa, MEDALUS model evolutions and prospects case study Sicily, in International Conference on Computational Science and Its Applications, Springer Nature Switzerland, (2023), 310–326. https://doi.org/10.1007/978-3-031-37114-1_21
    [8] X. M. Chen, D. M. Luo, Dynamical analysis of an almost periodic multispecies mutualism system with impulsive effects and time delays, Filomat, 37 (2023), 551–565. https://doi.org/10.2298/FIL2302551C doi: 10.2298/FIL2302551C
    [9] Z. Khan, J. Khan, S. Saifullah, A. Ali, Dynamics of hidden attractors in fourdimensional dynamical systems with power law, J. Funct. Spaces, 2022 (2022), 3675076. https://doi.org/10.1155/2022/3675076 doi: 10.1155/2022/3675076
    [10] H. A. Wood, R. R. Granados, Genetically engineered baculoviruses as agents for pest control, Annu. Rev. Microbiol., 45 (1991), 69–87. https://doi.org/10.1146/annurev.mi.45.100191.000441 doi: 10.1146/annurev.mi.45.100191.000441
    [11] V. Mondello, A. Songy, E. Battison, C. Pinto, C. Coppin, P. Trozel-Aziz, et al., Grapevine trunk diseases: A review of fifteen years of trials for their control with chemicals and biocontrol agents, Plant Dis., 102 (2018), 1189–1217. https://doi.org/10.1094/PDIS-08-17-1181-FE doi: 10.1094/PDIS-08-17-1181-FE
    [12] T. Shimales, E. Mendesil, B. Zewdie, B. Ayalew, K. Hylander, A. J. M. Tack, Management intensity affects insect pests and natural pest control on Arabica coffee in its native range, J. Appl. Ecol., 60 (2023), 911–922. https://doi.org/10.1111/1365-2664.14410 doi: 10.1111/1365-2664.14410
    [13] W. X. Yao, S. Guo, J. Wang, C. L. Chen, F. H. Yu, X. Li, et al., Droplet deposition and pest control efficacy on pine trees from aerial application, Pest. Manage. Sci., 78 (2022), 3324–3336. https://doi.org/10.1002/ps.6959 doi: 10.1002/ps.6959
    [14] Z. Arslangundogdu, E. Hizal, S. Acer, First record of oxycarenus lavaterae (Fabricius, 1787) (Heteroptera, Lygaeidae) in Turkey, Appl. Ecol. Environ. Res., 16 (2018), 1305–1311. https://doi.org/10.15666/aeer/1602_13051311 doi: 10.15666/aeer/1602_13051311
    [15] C. Q. Yang, Y. Fan, C. G. Du, A preliminary study on the Tropidothorax elegans (in Chinese), Insect Knowl., 25 (1988), 36–38.
    [16] W. S. Avusuglo, R. Mosleh, T. Ramaj, A. Li, S. S. Sharbayta, A. A. Fall, et al., Workplace absenteeism due to COVID-19 and influenza across Canada: a mathematical model, J. Theor. Biol., 572 (2023), 111559. https://doi.org/10.1016/j.jtbi.2023.111559 doi: 10.1016/j.jtbi.2023.111559
    [17] L. L. Chang, W. Gong, Z. Jin, G. Q. Sun, Sparse optimal control of pattern formations for an SIR reaction-diffusion epidemic model, SIAM J. Appl. Math., 82 (2022), 1764–1790. https://doi.org/10.1137/22M1472127 doi: 10.1137/22M1472127
    [18] Y. Tatsukawa, M. R. Arefin, S. Utsumi, K. Kuga, J. Tanimoto, Stochasticity of disease spreading derived from the microscopic simulation approach for various physical contact networks, Appl. Math. Comput., 431 (2022), 127328. https://doi.org/10.1016/j.amc.2022.127328 doi: 10.1016/j.amc.2022.127328
    [19] G. Q. Sun, X. Ma, Z. Z. Zhang, Q. H. Liu, B. L. Li, What is the role of aerosol transmission in SARS-Cov-2 Omicron spread in Shanghai, BMC Infect. Dis., 22 (2022), 880. https://doi.org/10.1186/s12879-022-07876-4 doi: 10.1186/s12879-022-07876-4
    [20] X. Ma, G. Q. Sun, Z. H. Wang, Y. M. Chu, Z. Jin, B. L. Li, Transmission dynamics of brucellosis in Jilin province, China: Effects of different control measures, Commun. Nonlinear Sci., 114 (2022), 106702. https://doi.org/10.1016/j.cnsns.2022.106702 doi: 10.1016/j.cnsns.2022.106702
    [21] G. Q. Sun, H. T. Zhang, L. L. Chang, Z. Jin, H. Wang, S. Q. Ruan, On the dynamics of a diffusive foot-and-mouth disease model with nonlocal infections, SIAM J. Appl. Math., 82 (2022), 1587–1610. https://doi.org/10.1137/21M1412992 doi: 10.1137/21M1412992
    [22] T. A. Biala, Y. O. Afolabi, A. Q. M. Khaliq, How efficient is contact tracing in mitigating the spread of COVID-19? A mathematical modeling approach, Appl. Math. Model., 103 (2022), 714–730. https://doi.org/10.1016/j.apm.2021.11.011 doi: 10.1016/j.apm.2021.11.011
    [23] C. A. Deutsch, J. J. Tewksbury, M. Tigchelaar, D. S. Battisti, S. C. Merrill, R. B. Huey, et al., Increase in crop losses to insect pests in a warming climate, Science, 361 (2018), 916–919. https://doi.org/10.1126/science.aat3466 doi: 10.1126/science.aat3466
    [24] H. Puebla, P. K. Roy, A. Velasco-Perez, M. M. Gonzalez-Brambila, Biological pest control using a model-based robust feedback, IET Syst. Biol, 12 (2018), 233–240. https://doi.org/10.1049/iet-syb.2018.5010 doi: 10.1049/iet-syb.2018.5010
    [25] L. James, F. D. Hugh, Interspecific hybridization and genetic divergence in Milkweed Bugs (Oncopeltus: Hemiptera: Lygaeidae), Evolution, 37 (1983), 583–591. https://doi.org/10.1111/j.1558-5646.1983.tb05575.x doi: 10.1111/j.1558-5646.1983.tb05575.x
    [26] A. Duarte, N. Borralho, P. Cabral, M. Caetano, Recent advances in forest insect pests and diseases monitoring using UAV-Based data: A systematic review, Forests, 13 (2022), 911. https://doi.org/10.3390/f13060911 doi: 10.3390/f13060911
    [27] L. Z. Ji, Z. Wang, X. W. Wang, Forest insect pest management and forest management in China: An overview, Environ. Manage., 48 (2011), 1107–1121. https://doi.org/10.1007/s00267-011-9697-1 doi: 10.1007/s00267-011-9697-1
    [28] J. N. Liu, Q. Qi, B. Liu, S. J. Gao, Pest control switching models with instantaneous and non-instantaneous impulsive effects, Math. Comput. Simul., 205 (2023), 926–938. https://doi.org/10.1016/j.matcom.2022.10.027 doi: 10.1016/j.matcom.2022.10.027
  • Reader Comments
  • © 2023 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(1237) PDF downloads(62) Cited by(0)

Figures and Tables

Figures(6)  /  Tables(1)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog