# Modeling the Parasitic Filariasis Spread by Mosquito in Periodic Environment

Authors: Yan Cheng; Xiaoyun Wang; Qiuhui Pan; Mingfeng He
Journal: Computational and Mathematical Methods in Medicine
 (2017)

---

## Abstract

In this paper a mosquito-borne parasitic infection model in periodic environment is considered. Threshold parameterR0 is given by linear next infection operator, which determined the dynamic behaviors of system. We obtain that when R0<1, the disease-free periodic solution is globally asymptotically stable and when R0>1 by Poincaré map we obtain that disease is uniformly persistent. Numerical simulations support the results and sensitivity analysis shows effects of parameters on R0, which provided references to seek optimal measures to control the transmission of lymphatic filariasis.


---

## Body

## 1. Introduction

Lymphatic filariasis is a parasitic disease caused by filarial nematode worms and is a mosquito-borne disease that is a leading cause of morbidity worldwide. Lymphatic filariasis affects 120 million humans in tropical and subtropical areas of Asia, Africa, the Western Pacific, and some parts of the Americas [1]. It is estimated that 40 million people are chronically disabled by lymphatic filariasis, making lymphatic filariasis the leading cause of physical disability in the world [2]. There are some clinical manifestations for infective individuals, such as acute fevers, chronic lymphedema, elephantiasis, and hydrocele [3].W. bancrofti parasites, which account for 90% of the global disease burden, dwell in the lymphatic system, where the adult female worms release microfilariae (mf) into the blood. Mf are ingested by biting mosquitoes as a blood meal of a mosquito, through several developmental stages, that is, first into immature larvae and then L3 larvae. Infective stage larvae L3 actively escape from the mosquito mouthparts entering another human host at the next blood meal through skin [4]. These L3 larvae subsequently develop into worms in humans and the process continues. So in order to remove lymphatic filariasis from the society, not only are the infected persons to be recovered but also the infected vectors are to be killed or removed.Mathematical models are powerful tools in disease control and may provide a powerful strategic tool for designing and planning control programs against infectious diseases [5]. Since 1960s, simple mathematical models of infection have been in existence for filariasis and provided useful insights into the dynamics of infection and disease in human populations [6–8]. Michael et al. describe the first application of the moment closure equation approach to model the sources and the impact of this heterogeneity for microfilarial population dynamics [9]. Simulation model for lymphatic filariasis transmission and control [10, 11] suggests that the impact of mass treatment depends strongly on the mosquito biting rate and on the assumed coverage, compliance, and efficacy; sensitivity analysis showed that some biological parameters strongly influence the predicted equilibrium pretreatment mf prevalence. References [12–14] take into account the complex interrelationships between the parasite and its human and vector hosts and provide the management decision support framework required for defining optimal intervention strategies and for monitoring and evaluating community-based interventions for controlling or eliminating parasitic diseases. Gambhir and Michael have shown a joint stability analysis of the deterministic filariasis transmission model [15]. All such models have proved to be of great value in guiding and assessing control efforts [16, 17].Environmental and climatic factors play an important role for the transmission of vector-borne diseases and are researched in many articles [18, 19]. For lymphatic filariasis, proper temperature and humidity are more beneficial for mosquito population to give birth and propagate. For example, in temperate climates and in tropical highlands, temperature restricts vector multiplication and the development of the parasite in the mosquito, while in arid climates precipitation restricts mosquito breeding. Therefore, the transmission of lymphatic filariasis exhibits seasonal behaviors especially in the northern areas [20, 21]. Nonautonomous phenomenon in infectious disease often occurs, and basic reproductive number of periodic systems is described as the spectral radius of the next infection operator [22].But the dynamics system considers the periodic environment between human and mosquito is little. How to make a comprehensive understanding of the role of periodic environment in the transmission of lymphatic filariasis and how to control the transmission of lymphatic filariasis efficiently are problems that provide motivation for our study. For the limitation of ecology environmental resources such as food and habitat, it is reasonable to adopt logistic growth for mosquito population. Nonautonomous logistic equations have been studied [23–28]. Based on above works and [29–34], we investigate a simple lymphatic filariasis model in periodic environment:(1)Sh′t=Λt-β1tShtImt1+α1tSht-μ1tSht+υtIht,Ih′t=β1tShtImt1+α1tSht-μ1tIht-υtIht,Sm′t=rtSmt1-SmtKt-β2tSmtIht,Im′t=β2tSmtIht-μ2tImt.In view of the biological background, system (1) has initial values(2)Sh00>0,Ih00>0,Sm00>0,Im00>0,where Sh(t) and Ih(t) separately denote the densities of the susceptible and the infective individuals for human population at time t; Sm(t) and Im(t) represent the densities of the susceptible and the infected individuals for mosquito population at time t, respectively. It is easy to see that Nh(t)=Sh(t)+Ih(t) and Nm(t)=Sm(t)+Im(t) are size of human population and mosquito population, respectively. Λ(t) is the recruitment rates of human host at time t; μ1(t) and μ2(t) are the death rate of human host and infected mosquito, including the natural death rate and disease-induced death rate; β1(t) and β2(t) denote the contact rate of infected mosquito to humans or infected humans to mosquito; α1(t) is the force of infection saturation at time t; υ(t) is the recovery rate of infectious human host at time t; r(t) and K(t) are the intrinsic growth rate and the carrying capacity of environment for mosquito population at time t, respectively.In view of the biological background of system (1), we introduce the following assumptions:(H1)
All coefficients are continuous, positiveω-periodic functions;(H2)
∫0ωr(t)dt>0.The organization of this paper is as follows. In Section2, some preliminaries are given and compute the basic production number. In Section 3, we will study the globally asymptotical stability of the disease-free periodic solution and the uniform persistence of the model. In Section 4, simulations and sensitive analysis are given to illustrate theoretical results and exhibit different dynamic behaviors.


## 2. Basic Reproduction Number

Denote(3)fL=supt∈0,ωft,fM=inft∈0,ωft,where f(t) is a continuous ω-periodic function.Let(Rk,R+k) be the standard ordered k-dimensional Euclidean space with a norm ·. For u,v∈Rk, we denote u≥v if u-v∈R+k,  u>v if u-v∈R+k∖{0}, and u≫v if u-v∈Int(R+k), respectively.LetA(t) be a continuous, cooperative, irreducible, and ω-periodic k×k matrix function; we consider the following linear system:(4)dxtdt=Atxt.Denote ΦA(t) be the fundamental solution matrix of (4) and let ρ(ΦA(ω)) be the spectral radius of ΦA(ω). Then by the Perron-Frobenius theorem, ρ(ΦA(ω)) is the principle eigenvalue of ΦA(ω) in the sense that it is simple and admits an eigenvector V∗≫0.Lemma 1 (see [35]).
Letp=1/ωln⁡ρ(ΦA(ω)), where A(t) is a continuous, cooperative, irreducible, and ω-periodic k×k matrix function. Then system (4) gives a solution x(t)=eptv(t), where v(t) is a positive ω-periodic function.When system (1) gives disease-free solution, obviously Ih(t)≡0 and Im(t)≡0. So we get the following subsystem:(5)Sh′t=Λt-μ1tSht,(6)Sm′t=rtSmt1-SmtKt.From Lemma 2.1 of [33] and Lemma 2 of [23] we obtain the following lemma.Lemma 2.
(i) System (5) has a unique positive ω-periodic solution Sh∗(t) which is globally asymptotically stable. (ii) System (6) has a globally uniformly attractive ω-periodic solution Sm∗(t).So, according to Lemma2, system (1) has a unique disease-free periodic solution (Sh∗(t),0,0,Sm∗(t)).In the following, we use the generation operator approach to define the basic reproduction number of (1). We check the assumptions (A1)–(A7) in [22] and denote x=(Ih(t),Im(t),Sh(t),Sm(t))T and (7)
F
t
,
x
=
β
1
t
S
h
t
I
m
t
1
+
α
1
t
S
h
t
β
2
t
S
m
t
I
h
t
0
0
,
V
-
t
,
x
=
μ
1
t
I
h
t
+
υ
t
I
h
t
μ
2
t
I
m
t
β
1
t
S
h
t
I
m
t
1
+
α
1
t
S
h
t
+
μ
1
t
S
h
t
r
t
S
m
2
t
K
t
+
β
2
t
S
m
t
I
h
t
,
V
+
t
,
x
=
0
0
Λ
t
+
υ
t
I
h
t
r
t
S
m
t
.So system (1) can be written as the following form:(8)x′t=Ft,xt-Vt,xt≡ft,xt,where V(t,x)=V-(t,x)-V+(t,x). From the expressions of F(t,x) and V(t,x), it is easy to see that conditions (A1)–(A5) are satisfied. We will check (A6) and (A7).Obviously,x∗(t)=(0,0,Sh∗(t),Sm∗(t)) is disease-free periodic solution of system (8). We define (9)Mt=∂fit,x∗t∂xj3≤i,j≤4,where fi(t,x∗(t)) and xi are the ith component of f(t,x(t)) and x, respectively. So we can get (10)Mt=-μ1t00rt-2rtKtSm∗t3≤i,j≤4.For Sm∗(t) is the globally uniformly attractively ω-periodic solution of (6), (11)∫0ωrt1-Sm∗tKtdt=0.Hence,(12)exp⁡∫0ωrt-2rtKtSm∗tdt=exp⁡-∫0ωrtKtSm∗tdt<1.It is easy to see that ρ(ΦM(ω))<1, and condition (A6) holds.Further, we define(13)Ft=∂Fit,x∗t∂xj1≤i,j≤2,Vt=∂Vit,x∗t∂xj1≤i,j≤2.Fi(t,x∗(t)) and Vi(t,x∗(t)) are the ith component of F(t,x∗(t)) and V(t,x∗(t)). So we obtain that(14)Ft=0β1tSh∗t1+α1tSh∗tβ2tSm∗t0,Vt=μ1t+υt00μ2t.Obviously ρ(Φ-V(ω))<1; thus condition (A7) holds.LetY(t,s) be 2×2 matrix solution of the following initial value problem: (15)dYt,sdt=-VtYt,s∀t≥s;Ys,s=I.I is identity matrix. Let Cω be the ordered Banach space of all ω-periodic functions from R→R2, which is equipped with maximum norm ·∞ and the positive cone Cω+={ϕ∈Cω:ϕ(t)≥0,∀t∈R}. By the approach in [22], we consider the following linear operator L:Cω→Cω. Suppose that ϕ(s)∈Cω is the initial distribution of infectious individuals in this periodic environment. F(s)ϕ(s) is the distribution of new infections produced by the infected individuals who were introduced at time s, and Y(t,s)F(s)ϕ(s) represents the distributions of those infected individuals who were newly infected at time s and remain in the infected compartment at time t. Then (16)ψt=∫-∞0Yt,sFsϕsds=∫0+∞Yt,t-aFt-aϕt-adadenotes the distribution of accumulative new infections at time t produced by all those infected individuals ϕ(s) introduced at previous time to t. (17)Lϕt=∫0+∞Yt,t-aFt-aϕt-ada,∀t∈R,ϕ∈Cω.As in [22], L is the next infection operator, and the basic reproduction number of system (1) is given by (18)R0=ρL,where ρ(L) is the radius of L. Next we show that R0 serves as a threshold parameter for the local stability of the disease-free periodic solution.Theorem 3 (see Wang and Zhao [22], Theorem 2.2).
Assume that (A1)–(A7) hold; then the following statements are valid:(i)
R0=1 if and only if ρ(ΦF-V(ω))=1;(ii)
R0>1 if and only if ρ(ΦF-V(ω))>1;(iii)
R0<1 if and only if ρ(ΦF-V(ω))<1.So the disease-free periodic solution(Sh∗(t),0,0,Sm∗(t)) is asymptotically stable if R0<1 and unstable if R0>1.


## 3. Global Stability of Disease-Free Periodic Solution

Denote(19)Ω=Sh,Ih,Sm,Im:Sh>0,Ih≥0,Sm≥0,Im≥0,0<Sh+Ih≤ΛLμ1M<+∞,0≤Sm+Im≤Mm∗Δ<+∞.Ω is a positively invariant set with respect to system (1) and a global attractor of all positive solutions of system (1).(20)Nh′t=Λt-μ1tNht≤ΛL-μ1MNht,where ΛL=supt>0Λ(t) and μ1M=inft>0μ1(t). So it is easy to obtain Nh(t)≤ΛL/μ1M. (21)Nm′t=rtSmt1-SmtKt-μ2tImt≤rt+μ2tSmt-μ2tNmt≤Mm-μ2tNmt,where Mm=supt∈[0,ω)(r(t)+μ2(t))Sm(t).From the third equation of (1), for all t≥0 we have(22)Sm′t≤rtSmt1-SmtKt;by the comparison principle and Lemma2, we obtain(23)limt→∞sup⁡Smt≤limt→∞sup⁡Sm∗t≤Mm∗,where Sm∗(t) is the globally uniformly attractively positive ω-periodic solution and Mm∗=maxt∈[0,ω]Sm∗(t). So, for any small ϵ existing a t0, for all t≥t0 we have(24)Smt≤Sm∗t+ε≤Mm∗+ϵ.So we obtain(25)Nm′t≤supt≥0rt+μ2tMm∗+ϵ-μ2tNmt,and limt→∞sup⁡Nm(t)≤(Mm∗+ϵ)Δ, where Δ=supt>0(r(t)+μ2(t))/inft>0μ2(t). For ϵ small enough, Nm(t)≤Mm∗Δ.Theorem 4.
IfR0<1, the disease-free periodic solution (Sh∗(t),0,Sm∗(t),0) is globally asymptotically stable. And if R0>1, it is unstable.Proof.
By Theorem3 we obtain that if R0<1, (Sh∗(t),0,Sm∗(t),0) is locally stable. Next we prove that when R0<1 the disease-free solution (Sh∗(t),0,Sm∗(t),0) has global attractivity.
WhenR0<1 and by (iii) of Theorem 3, we have ρ(ΦF-V(ω))<1. So there exists a small enough constant ε1>0 such that ρ(ΦF-V+ε1N(ω))<1, where (26)Nt=0β1t1+α1tSh∗t+ε1β2t0.From Lemma 2 and nonnegativity of the solutions, for any ε1>0 there exists t1>0 such that Sh(t)≤Sh∗(t)+ε1 and Sm(t)≤Sm∗(t)+ε1, so for all t>t1 we have(27)Ih′t≤β1tSh∗t+ε1Imt1+α1tSh∗t+ε1-μ1tIht-υtIht,Im′t≤β2tSm∗t+ε1Iht-μ2tImt.Considering the auxiliary system(28)Ih′t~=β1tSh∗t+ε1Im′t~1+α1tSh∗t+ε1-μ1tIh′t~-υtIh′t~,Im′t~=β2tSm∗t+ε1Ih′t~-μ2tIm′t~.From Lemma 1, it follows that there exists a positive ω-periodic solution v1(t) such that J(t)≤eptv1(t), where J(t)=(Ih(t),Im(t))T and p=1/ωln⁡ρ(ΦF-V+ε1N(ω))<0. Then limt→∞J(t)=0; that is, limt→∞Ih(t)=0 and limt→∞Im(t)=0.
Moreover, from the equations ofSh(t),Sm(t), we get(29)limt→∞Sht=Sh∗t,limt→∞Smt=Sm∗t.Hence, disease-free periodic solution of system (1) is globally attractive. This completes the proof.Define(30)X=Sh,Ih,Sm,Im:Sh>0,Ih≥0,Sm≥0,Im≥0,X0=Sh,Ih,Sm,Im∈X:Ih>0,Im>0.We have (31)∂X0=X∖X0=Sh,Ih,Sm,Im∈X:IhIm=0.From system (1), it is easy to see that X and X0 are positively invariant, and ∂X0 is also a relatively closed set in X.LetP:X→X be the Poincaré map associated with system (1), satisfying (32)Px0=uω,x0,∀x0∈X;u(t,x0) is the unique solution of system (1) satisfying initial condition u(0,x0)=x0. P is compact for the continuity of solutions of system (1) with respect to initial value, and P is point dissipative on X.We further define(33)M∂=Sh0,Ih0,Sm0,Im0∈∂X0:PmSh0,Ih0,Sm0,Im0∈∂X0∀m>0,where Pm=P(Pm-1) for all m>1 and P1=P. Now, prove (34)M∂=Sh0,0,Sm0,0:Sh0>0,Sm0≥0.Obviously {(Sh,0,Sm,0):Sh>0,Sm≥0}⊆M∂.IfM∂∖{(Sh,0,Sm,0):Sh>0,Sm≥0}≠∅, then there exists at least a point (Sh0,Ih0,Sm0,Im0)∈M∂ satisfying Ih0>0 or Im0>0. We consider two possible cases.IfIh0=0 and Im0>0, then it is clear that from system (1) Im(t)≥0 for any t>0. From the second equation of system (1) and Sh>0, we obtain (35)Iht=Ih0e-∫0tμ1s+υsds+∫0tβ1sShsIms1+α1sShse∫stμ1τ+υτdτds>0,for all t>0.IfIm0=0 and Ih0>0, then Ih(t)=Ih0e-∫0t[μ1(τ)+υ(τ)]dτ>0. From the third equation of system (1) and Sm>0, we obtain (36)Imt=Im0e-∫0tμ2sds+∫0tβ2sSmsIhse∫stμ2τdτds>0,for all t>0. Hence, for any case, it follows that (Sh(t),Ih(t),Sm(t),Im(t))∉∂X0, so (Sh0,Ih0,Sm0,Im0)∉M∂. This leads to a contradiction; there exists one fixed point E0=(Sh∗(t),0,Sm∗(t),0) of P in M∂.In the following, we will discuss the uniform persistence of the disease, andR0 serves as a threshold parameter for the extinction and the uniform persistence of the disease.Theorem 5.
IfR0>1, then system (1) is uniformly persistent. There exists a positive constant ε, such that for all initial conditions (1) satisfies(37)limt→∞inf⁡Iht≥ε,limt→∞inf⁡Imt≥ε.When R0>1, system (1) admits at least one positive periodic solution.Proof.
From Theorem3, if R0>1 then we obtain ρ(ΦF-V(ω))>1. For an arbitrary small constant η>0, that ρ(ΦF-V-ηN(ω))>1, N(t) is the same as in Theorem 3. From assumption (H2), we obtain any small enough ε>0, ∫0ω[r(t)-α(t)ε]dt>0. Consider perturbed equations(38)Sεh′t=Λt-εβ1tSεht1+α1tSεht-μ1tSεht,(39)Sεm′t=rtSεmt1-SεmtKt-εβ2tSεmt.Using Lemma 2 in [25] and Lemma 1 of [27], we obtain (38) and (39) that admit globally uniformly attractive positive ω-periodic solutions Sεh∗(t) and Sεm∗(t). For the continuity of solutions with respect to ε, and for η>0 there exists ε1>0 for all t∈[0,ω]; thus we have(40)Sε1m∗t>Sm∗t-η,Sε1h∗t>Sh∗t-η.
Denotex0=(Sh0,Ih0,Sm0,Im0)∈X0, according to the continuity of the solution with respect to the initial condition; there exists δ for given ε1, for all x0∈X0 with x0-E0<δ; it follows u(t,x0)-u(t,E0)<ε1 for all t∈[0,ω].
Following, we prove(41)limm→∞sup⁡dPmx0,E0≥δ.We suppose the conclusion is not true; then following inequality holds:(42)limm→∞sup⁡dPmx0,E0<δ,for some x0∈X0. Without loss of generality, we can assume that(43)dPmx0,E0<δ∀m≥0.So we obtain(44)ut,Pmx0-ut,E0<ε1∀m≥0,t∈0,ω.For any t≥0,  t=mω+t′, where t′∈[0,ω] and m=[t/ω] is the greatest integer less than or equal to t/ω, so we have(45)ut,x0-ut,E0=ut′,Pmx0-ut′,E0<ε,∀t≥0.Hence, it follows that 0≤Ih(t)≤ε1 and 0≤Im(t)≤ε1 for all t≥0. Then from the first and third equations of (1),(46)Sh′t≥Λt-ε1β1tSht1+α1tSht-μ1tSht,Sm′t=rtSmt1-SmtKt-ε1β2tSmt.By the comparison principle, we obtain for any t≥0(47)Sht≥Sε1ht,Smt≥Sε1mt.Consider (38); there exists t1>0; for all t>t1 we have(48)Sε1mt>Sε1m∗t-η,Sε1ht>Sε1h∗t-η.By (38) and (48) we obtain(49)Smt>Sm∗t-η,Sht>Sh∗t-η.Then for all t>t1 we have (50)Iht≥β1tSh∗t-ηImt1+α1tSh∗t-η-μ1tIht-υtIht,Imt≥β2tSm∗t-ηIht-μ2tImt.Consider the following auxiliary system:(51)Iht¯=β1tSh∗t-ηImt¯1+α1tSh∗t-η-μ1tIht¯-υtIht¯,Imt¯=β2tSm∗t-ηIht¯-μ2tImt¯.From Lemma 1, it follows that there exists a positive ω-periodic function v2(t) such that (51) has a solution J(t)=v2(t)ep1t, where p1=1/ωln⁡ρΦF-V-ηNω. For ρΦF-V-ηNω>1, (52)limt→∞Iht=+∞,limt→∞Imt=+∞.This leads to a contradiction.That is to say,M∂∖{(Sh,0,Sm,0):Sh>0,Sm≥0}=∅ and {M1} is globally attractive in M∂, and all orbit in M∂ converges to {M1}. By [22], we obtain that P is weakly uniformly persistent with respect to (X0,∂X0). All solutions are uniformly persistent with respect to (X0,∂X0); thus we have limt→∞Ih(t)≥ε,  limt→∞Im(t)≥ε.


## 4. Sensitivity Analysis and Prevention Strategy

We conducted numerical simulation to this model and computed the reproductive numbersR0. It was confirmed that using the basic reproduction number of the time-averaged autonomous systems of a periodic epidemic model overestimates or underestimates infection risks in many other cases. Bacaer and Guernaoui give methods to compute R0, such as method of discretization of the integral eigenvalue [36] and Fourier series method for general periodic case and sinusoidal case and application of Floquet Theory method [37]. In [22] Wang and Zhao propose that in order to compute R0 we only need to compute the spectrum of evolution operator of the following system (53):(53)dwdt=-Vt+Ftλw,w∈Rn,λ∈0,∞;here system (53) is ω-periodic equation, and W(t,s,λ) is the evolution operator of system (53) with t≥s, s∈R. By Perron-Frobenius theorem ρ(W(ω,0,λ)) is an eigenvalue of W(t,0,λ), t≥0. Next, using Theorem 2.1 in [22] to compute R0 numerically, R0 serves as threshold parameter in periodic circumstances.Firstly, by the means of the software Matlab we computeR0. We choose parameters Λ(t)=0.6+0.4sin⁡(2πt/12), μ1(t)=0.5+0.1sin⁡(2πt/12), μ2(t)=0.8+0.1sin⁡(2πt/12), β1(t)=0.6+0.1sin⁡(2πt/12), β2(t)=0.7+0.1sin⁡(2πt/12), α1(t)=0.2+0.1sin⁡(2πt/12), υ(t)=0.02+0.03sin⁡(2πt/12), r(t)=0.5+0.4sin⁡(2πt/12), K(t)=0.9+0.3sin⁡(2πt/12). By numerical calculations, we obtain R0=0.9243<1; then the disease will be extinct; see Figure 1(a). If we choose β1(t)=0.9+0.1sin⁡(2πt/12), β2(t)=1.2+0.1sin⁡(2πt/12), then R0=1.4662>1; the disease is permanent; see Figure 1(b). The evolution trajectory in spaces (Sh,Ih) and (Sm,Im) are in Figures 2(a) and 2(b), respectively.Figure 1
Plot the evolution tendency of four populations. (a) Fixed parametersβ1(t)=0.6+0.1sin⁡(2πt/12),  β2(t)=0.7+0.1sin⁡(2πt/12); then R0=0.9243<1; (b) Parameters β1(t)=0.9+0.1sin⁡(2πt/12),  β2(t)=1.2+0.1sin⁡(2πt/12); then R0=1.4662>1.
(a) (b)Figure 2
WhenR0=1.4662, we graph the trajectory of two populations in spaces (Sh,Ih) and (Sm,Im), respectively.
(a) (b)In order to perform sensitivity analysis of parametersβ1(t), β2(t), K(t), and α1(t), we fix all parameters as in Figure 1, except that we choose the composite functions as follows: (54)β1t=β01+0.1sin⁡2πt12,β2t=β02+0.1sin⁡2πt12,Kt=k0+0.3sin⁡2πt12,α1t=α0+0.1sin⁡2πt12,where β01=1/12∫012β1(t)dt,  β02=1/12∫012β2(t)dt,  k0=1/12∫012K(t)dt, and α0=1/12∫012α1(t)dt.We first fix other parameters and detect the effect of parameters ofk0 and α0 on R0. From Figure 3(a), we see that with the increase of α0, R0 decreases, and the gradient also decreases, so this strengthens the psychological hint of susceptible human individuals to be benefit for the extinction of the disease. In Figure 3(b), with the increasing of k0 the sensitivity of R0 increases. That is to say, the carrying capacity of environment for mosquito is bigger and the disease is widespread more easily, so decreasing the circumstance fit survival for mosquitoes, such as contaminated pool or puddle and household garbage, is a necessary method for the extinction of disease.Figure 3
Sensitivity analysis of the basic reproductionR0 with parameter k0 or α0.
(a) (b)Next, we consider the combined influence of parametersβ10 and β20 on R0; in Figure 4 we can see that the basic reproduction number R0 may be less than 1 when β10 and β20 are small; the smaller β20 the more sensitive the effect on R0.Figure 4
Sensitivity analysis of the basic reproductionR0 with parameters β01 and β02.In Figure5, the basic reproduction number R0 is affected by β10 and k0; with the increasing of k0 the sensitivity of R0 increases; if we fix β10 as a constant the case will be similar to Figure 3(b). And the similar trend of β10 on the sensitivity of R0, so in the season in which temperature and humidity are more beneficial for mosquito population to give birth and propagate taking measures to avoid more bites is necessary.Figure 5
Sensitivity analysis of the basic reproductionR0 with parameters β01 and k0.


## 5. Conclusion

In this paper, we have studied the transmission of lymphatic filariasis; lymphatic filariasis is a mosquito-borne parasitic infection that occurs in many parts of the developing world. In order to systematically investigate the impact that vector genus-specific dependent processes may have on overall lymphatic filariasis transmission, we, according to the nature characteristic of lymphatic filariasis and considering the logistic growth in periodic environments of mosquito, model the transmission of lymphatic filariasis. The dynamic behavior of system (1) is determined by the threshold parameter R0; when R0<1 disease-free periodic solution is globally asymptotically stable and when R0>1 disease is uniformly persistent. We also give some numerical simulations which support the results we prove, confirming that R0 serves as a threshold parameter. Sensitivity analysis show effects of parameters on R0, which contribute to providing a decision support framework for determining the optimal coverage for the successful prevention programme.



---

*Source: 4567452-2017-02-08.xml*