Modeling seasonal variation for mosquito-borne disease in the tropical monsoon environment

Mosquitoes play an important role in the spread of mosquito-borne diseases. Considering the sensitivity of mosquitoes’ aquatic stage to the seasonal shift, in this paper, we present a seasonally forced mosquito-borne epidemic model by incorporating mosquitoes’ aquatic stage (eggs, larvae, and pupae) and seasonal shift factor, which is a periodic discontinuous differential system. Firstly, some sufficient conditions for the existence and uniqueness of a disease-free solution are obtained. Further, we define the basic reproduction number R0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\mathcal{R}_{0}$\end{document}, and obtain the stability of the disease-free solution when R0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\mathcal{R}_{0}$\end{document} is less than one. And, if R0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$\mathcal{R}_{0}$\end{document} is greater than one, the mosquito-borne disease is uniformly persistent and the model admits a positive periodic solution. Finally, some numerical simulations are given to illustrate the main theoretical results. In addition, simulation results also imply that ignoring the effects of seasonal succession can overestimate or underestimate mosquito-borne disease trends.


Introduction
In recent years, the spread of mosquito-borne diseases (such as dengue fever, malaria, Zika, chikungunya, yellow fever, and so on) has been characterized by high morbidity, high mortality, rapid growth in the number of cases, etc., which has become the major public health concern in the tropical and sub-tropical regions of the world. For example, in 2017, there were estimated 219 million cases of malaria in 87 countries, of which about 435,000 deaths [16]. According to the Philippine ministry of health, in 2019, more than 387,000 dengue cases (including 1452 deaths) have been registered in the Philippines as of November 2, up 99.3% from the same period last year [14].
Mathematical modeling is an important tool to understand the epidemiology of a mosquito-borne disease and assist in controlling the mosquito-borne disease. The first mosquito-borne disease model was developed to study the strategies to control malaria by Ross [21] in 1911. Subsequently, many scholars extended this model by considering acquired immunity in human population [5,24], age structure in human population [27,29], spatial heterogeneity [1,34], the incubation period of pathogen in the mosquito [26,32,33], the life-cycle of mosquito [11,12,35], and so on. In addition, many kinds of dynamical models of mosquito-borne disease have been developed to put forward relevant strategies for disease control. We refer to some of them [3,7,8] and the references therein.
It is well known that mosquito-borne diseases are influenced by complex demographic, environmental, and social factors. Temperature and rainfall are key environmental factors that determine the population size and the range of mosquitoes because they are critical to the length of mosquitoes' life cycle, the rate of reproduction and development of mosquitoes [20]. Therefore, many scholars have established periodic systems to study the effects of seasonality on mosquito-borne disease transmission (see [9,17,31]). Specifically, Mathieu et al. [4] established a vector-host model accounting for seasonal fluctuation in vector density and found that the periodic fluctuations of the vector population in phase with temperature variations suggest a strong climate effect on the vector density and, in turn, on the transmission of vector-borne. Lou et al. [15] presented and analyzed a mathematical malaria model with periodic environment. They obtained that dynamic behaviors of the model are determined by its basic reproduction number, and found that even small changes in temperature could change the prevalence of malaria.
To the best of our knowledge, in the existing models of seasonal mosquito-borne diseases, the effects of seasonality on the transmission of mosquito-borne disease are mainly characterized by a periodic continuous differential system. However, the periodic continuous differential system may not be appropriate for mosquito-borne disease modeling in regions with tropical monsoon climate, such as Mumbai (India), Colombo (Sri Lanka), Bangkok (Thailand), Ho Chi Minh City (Viet Nam), Phnom Penh (Cambodia), and so on mainly because there are only a dry season (Nov-Apr) and a rainy season (May-Oct) throughout the year in these areas, and the rainfall in the rainy season is significantly higher than that in the dry season. In fact, seasonal rainfall can increase the abundance of mosquitoes with aquatic larval stages, where reproduction depends on the availability of breeding sites [2]. Therefore, it is of great significance to consider a discontinuous periodic differential model to characterize the seasonal outbreak of mosquito-borne disease caused by seasonal change between rainy and dry seasons.
Based on the above discussion, in this paper, we propose a mosquito-borne disease transmission model with seasonal variation, which is a periodic discontinuous differential system. Applying the calculation method of the basic reproduction number for the continuous periodic differential system given in Wang and Zhao [30] and the theoretical analysis method of the discontinuous periodic differential system given in Tang et al. [23], we calculate the basic reproduction number for this model and discuss the uniform persistence of the disease and threshold dynamics. The organization of this paper is as follows. In the next section, model formulation, some hypotheses, and the well-posedness of this model are given. In Sect. 3, the basic reproduction number is calculated, and we show that it is the threshold value for uniform persistence and global extinction of the disease. In Sect. 4, we give some numerical simulations and discussion to illustrate our theoretical results, and we give some concluding remarks in the last section.

Model formulation and preliminaries
According to the mosquito-borne epidemic compartmental models in previous studies [22,28], the total human population at time t, denoted by N h (t), is divided into three classical classes: the susceptible S h (t), the infected I h (t), and the recovered R h (t). For the mosquito population, we consider the aquatic mosquitoes (comprising eggs, larvae, and pupae) and adult mosquitoes, where the adult population at time t, denoted by N m (t), is split into two classes: susceptible mosquitoes S m (t) and infected mosquitoes I m (t), the aquatic population at time t is denoted by A(t). Further, the following basic assumptions are introduced for our model.
(H 1 ) Considering the influence of seasonal factors (such as humidity, temperature, etc.) on mosquitoes' reproduction rate and breeding sites, we assume that the recruitment rate of the aquatic mosquitoes is governed by a logistic equation, in which these coefficients are periodic functions on account of seasonal effects. (H 2 ) For human population, let Λ h , μ h , d h , γ , and η be the human recruitment rate, natural death rate, disease induced death rate, recovery rate, and lose immunity rate, respectively. Denote the seasonal forced transmission coefficient from infected mosquitoes to a susceptible human by β(t), which is affected by seasonal succession. (H 3 ) For mosquito population, let k(t), b(t), μ a (t), and μ m (t) represent the carrying capacity of aquatic mosquitoes, mosquito birth rate, aquatic mosquito death rate, and adult mosquito death rate, respectively. The per capita rate of maturation from the aquatic form to the adult one is denoted by φ(t). Let α(t) denote the transmission coefficient from the infected human to susceptible mosquitoes. (H 4 ) Let ω be given as the period of disease transmission. Let J 1 and J 2 denote the dry season and the rainy season, respectively, where J 1 = [nω, nω + θω) and J 2 = [nω + θω, (n + 1)ω). Here, n ∈ Z + and θ ∈ (0, 1) is the fraction of the rainy season to the whole infection cycle. According to the assumptions above and taking into account the interaction between human and mosquito population, we give the following ordinary differential equations: with the following piecewise constant functions: and where J 1 = [nω, nω + θω) and J 2 = [nω + θω, (n + 1)ω). Assume In view of the biological background of model (1), we only consider the solution of model (1) with initial conditions (S h (0), I h (0), R h (0), A(0), S m (0), I m (0)) ∈ R 6 + = {(x 1 , . . . , x 6 ) : x i ≥ 0, i = 1, . . . , 6}. Now, we can give the following results about the existence, uniqueness, and boundedness of solutions of model (1).

Theorem 1 For any initial value condition
for all t ≥ 0, and this solution is ultimately bounded.
The global existence and uniqueness of the solution can be followed by Theorem 2.1 in [23]. The positivity and boundedness can be easily obtained, here omitted.
Consider the following two disease-free subsystems of model (1): and On the global dynamics of models (2) and (3), we have the following lemmas.

Lemma 1 Model
(2) has a unique positive globally asymptotically stable equilibrium The proof of Lemma 1 is obvious, so it is omitted.
Proof By the similarly method of Theorem 3 and Theorem 4 given in Ref. [10], it can be obtained that when (3) is persistent and has one positive periodic solution (A * (t), S * m (t)). Now, we will prove global uniform attractiveness of this periodic solution.
Let (A(t), S m (t)) be any positive solution of model (3). Defining a function V (t) = |A(t) -A * (t)| + |S m (t)-S * m (t)| and calculating the upper right derivative of V (t) along the solutions of model (3), it follows that Obviously, there exist positive constants 1 , 2 and T 1 ≥ 0 such that Integrating both sides of the above inequality on interval [T 1 , t] yields that . By Barbalat's lemma [13], we conclude that lim t→∞ |A(t) -A * (t)| = 0, lim t→∞ |S m (t) -S * m (t)| = 0. This completes the proof.

Basic reproduction number and threshold dynamics
In this section, we introduce the threshold quantity for the piecewise periodic model (1) and analyze its properties. First, we need to define the next infection operator for model (1), which arises from the combination with the idea in Ref. [30] for periodic ordinary differential models. Let and Then model (1) can be rewritten as where x), and , respectively. Then, by simple computations, it follows that Now, we extend t ∈ R + to t ∈ R and introduce some new notations. When t ∈ +∞ n=-∞ (J 1 ∪ J 2 ) = (-∞, +∞) and Obviously, F and V are 2 × 2 piecewise continuous periodic matrices with period ω in R, and -V (t) is cooperative. Let Y(t, s), t ≥ s, be the evolution operator of the linear system Then where E is the 2 × 2 identity matrix. Therefore, the monodromy matrix Φ -V (t) of model (6) equals Y(t, 0), t ≥ 0. According to Remark 3.5 in Sect. III.7 of Ref. [6], there exist K 1 > 0 and K 2 > 0 such that Obviously, there exists a constant K 3 > 0 such that F(t) < K 3 followed by the boundedness of F(t). Therefore, Let C ω be the ordered Banach space of all ω-periodic continuous functions from R to R 2 with the maximum norm · c , and C + ω = {I ∈ C ω : I(s) ≥ 0, s ∈ R}. Assume that I(s) ∈ C + ω is the initial distribution of infectious individuals in this periodic environment, then F(s)I(s) is the rate of new infectious produced by the infected individuals who were introduced at time s, and Y(t, s)F(s)I(s) represents the distributions of those infected individuals who were newly infected at time s and remain in the infected compartment at time t for t ≥ s.
which denotes the distribution of accumulative new infections at time t produced by all those infected individuals I(s) introduced at previous time to t. Define a linear operator L : Similar to Lemma 3.1 given in Ref. [23], the linear operator L is well defined.

Lemma 3
The operator L is positive, continuous, and compact on C ω . Therefore, following Wang and Zhao in [30], we call L the next infection operator and define the basic reproduction number R 0 for model (1) by where ρ(L) is the spectral radius of L. Similar to the analysis of Ref. [23], it is easy to follow the argument in [29] to characterize R 0 . Set U λ (t, s)(t ≥ s, s ∈ R) be the evolution operator of the following auxiliary ω-periodic switching model: Then U(ω, 0, λ) = Φ F-V (ω). Following the ideas in Ref. [30], we have the following result.
Next, we explore the threshold condition, which determines the extinction and persistence of the disease. 0) is globally asymptotically stable, where Θ(t) is given by (4).
It follows from an epidemiological concept that if R 0 < 1 the number of infected human and infectious mosquitoes tends to zero with the increase of time and the disease is extinct. However, the disease-free periodic solution is unstable, that is, no matter how small the number of infected hosts or infectious vectors in the initial state, the disease will occur repeatedly in the region and become endemic due to its long-term prevalence. In the following, we show that when R 0 > 1, the disease is uniformly persistent.
Proof Define From model (1), it is easy to see that X and X 0 are positively invariant and ∂X 0 is also a relatively closed set in X.
Let P : X → X be the Poincaré map associated with model (1), that is, P(P 0 ) = u(ω, P 0 ) = u 2 ω, θω, u 1 (θω, 0, P 0 ) , where u(t, P 0 ) is a solution of model (1) under the initial value P 0 ∈ R 5 + , and u i (t, t * , P * ) (i = 1, 2) is the solution semi-flow of the following model: It follows from the continuity of solution of model (1) with respect to the initial value P 0 that P is compact. Moreover, by Theorem 1, we obtain that P is point dissipative on X. Define It is a positive invariant set of P in ∂X 0 . Now, we claim Clearly, if φ 2 b 2 > μ m1 (φ 1 + μ a1 ) and lim inf t→∞ Θ(t) > 0, model (12) has a globally stable ω-periodic solution (Λ h /μ h , A * (t), S * m (t)). Therefore, the map P has a unique globally attractive fixed point restricted in M ∂ , which is M 1 = (Λ h /μ h , 0, 0, A * (t), S * m (0), 0). Define W s (M 1 ) = {x 0 : P n (x 0 ) → M 1 , n → ∞}, which is said to be a stable set of M 1 . Since R 0 > 1, the stable set W s (M 1 ) of M 1 satisfies that W s (M 1 ) ∩ X 0 = ∅. It follows from Lemmas 1 and 2 that {M 1 } is globally attractive in M ∂ when φ 2 b 2 > μ m1 (φ 1 + μ a1 ) and lim inf t→∞ Θ(t) > 0. Therefore, {M 1 } is isolated in M ∂ , and hence in X. Furthermore, {M 1 } is also invariant and {M 1 } does not form a cycle in M ∂ , and hence in X. By Theorem 1.3.1 in [34], we finally obtain that P is uniformly persistent with respect to (X 0 , ∂X 0 ), which implies that model (1) is uniformly persistent. The proof is completed.
From the above expressions, we can see that α, β, k, and φ are proportional to R 0 , while μ m and μ a are inversely proportional to R 0 . Using the sensitiveness analytical method [18], it is easily obtained that μ m is the most sensitive parameter about the basic reproduction number R 0 , followed by α, β, k and φ, b, and μ a are the least. Therefore, in order to effectively control the spread of the disease, during the rainy season, we can spray adult and larval insecticides, reduce contact rates, and reduce mosquito breeding sites at the same time. On the other hand, in view of pesticides causing environmental pollution, we can avoid the use of pesticides as much as possible in the dry season and control the spread of disease by reducing contact rates and reducing mosquito breeding sites.
In order to investigate the effect of the seasonal succession on the disease-free periodic solution of model (1), we choose parameters β = 2.15 × 10 -7 , α = 4.65 × 10 -7 , i.e., there is no seasonal variation. By direct calculation, we obtain the basic reproduction number R 0 ≈ 0.9143 of the autonomous model, which is close to the basic reproduction number R 0 ≈ 0.9177 < 1 of model (1) with rainy season and dry season. The solution curves of two  Fig. 1(d), respectively. By contrast of the trend of the black dotted line and the red-blue solid line, we find that the black dotted line undergoes a little bit of growth and then goes down very quickly and goes to zero; while the red-blue solid line has a relatively long growth and then gradually approaches zero after several disease cycles. These results indicate that even when the basic reproduction numbers are close to each other, it is easy to cause misjudgment for disease control departments without considering the seasonal succession. At the same time, we choose the above rainy season parameters for the autonomous model: β = 5.15 × 10 -7 , α = 5.65× 10 -7 , i.e., there is no seasonal variation. The simulation results are denoted by green chain dotted line in Fig. 1(d). Comparing the seasonal succession of rainy and dry, we find the basic reproduction number R 0 ≈ 1.5597 and the disease will not be extinct. That is to say, the disease will be overestimated without considering the dry season. The numerical simulation implies that considering the mosquito-borne epidemic model with seasonal change can more accurately describe the epidemic law of the mosquito-borne disease.
Furthermore, we choose rainy season parameters β 1 = 1.35 × 10 -6 , α 1 = 2.15 × 10 -6 and dry season parameters β 2 = 5.5 × 10 -7 , α 2 = 8.75 × 10 -7 in model (1). Clearly, φ 2 b 2 > μ m1 (φ 1 + μ a1 ), and by numerical calculations, one has lim inf t→∞ Θ(t) > 0 and obtains the basic reproduction number R 0 ≈ 2.2678 > 1. According to Theorem 4, model (1) is uniformly persistent, as shown in Fig. 2. In addition, the plots in Fig. 2(a) and Fig. 2(b) show that model (1) has a positive periodic solution, which is in accordance with the obtained conclusion in Corollary 1. It seems that this positive periodic solution is globally attractive from the trend of solution curves with different initial values in Fig. 2(c) and Fig. 2(d).  Finally, in order to investigate the effect of seasonal succession on the persistence of mosquito-borne disease, we take the autonomous model parameters β = 8.25 × 10 -7 , α = 1.23 × 10 -6 , i.e., there is no seasonal succession. By calculating, we obtain the basic reproduction number R 0 ≈ 2.2784, which is close to the basic reproduction number R 0 ≈ 2.2678 > 1 in Fig. 2, where the seasonal succession is considered. This is shown in the black dotted line and the red-blue solid line in Fig. 3. It is not hard to find that the black dotted line undergoes a little bit of growth and then gradually goes to one value, while the red-blue solid line has a relatively big growth and then gradually forms a stable state of periodic fluctuation. Numerical simulations indicate that although the basic regeneration numbers of the two cases are very close, the prevalence of the disease is completely different. Further, we take only the rainy season parameters β = 1.35 × 10 -6 , α = 2.15 × 10 -6 in the autonomous model and get R 0 ≈ 3.8224 in this case. A similar result can be obtained, as shown by the green chain dotted line in Fig. 3. Theoretical results and numerical simulations show that neglecting the effect of seasonality on mosquitoes can lead to overestimating or underestimating the epidemic trend of mosquito-borne infectious diseases and cause waste of medical resources to a certain extent.

Conclusion
In this paper, we establish a mosquito-borne disease transmission model with seasonal succession, which is different from the existing work due to the discussion of the dynamics of transmission of mosquito-borne diseases in certain areas. These areas have a tropical monsoon climate, year-round high temperatures (suitable for mosquito growth), and only two seasons a year: the dry season and the rainy season. The rainfall difference between the two seasons is very obvious. Therefore, we use piecewise constant functions to describe the seasonal succession rather than the continuous periodic function. All in all, we can conclude our work and findings in the following two aspects: (a) The main results obtained in this work: combined with the analysis methods in Refs. [30] and [23], the existence and uniqueness of the disease-free periodic solutions of this model are discussed, and the basic regeneration number of our model is obtained. On the premise of the existence of disease-free periodic solution, the basic regeneration number is the threshold that determines the extinction and persistence of the disease. Specifically, if R 0 < 1, then the disease-free periodic solution of this model is globally asymptotically stable; if R 0 > 1, the disease in the model is persistent and the model has a positive periodic solution. (b) The effect of seasonal succession on the dynamics of the disease: in order to investigate the effect of seasonal succession on the dynamics of mosquito-borne disease, we give two other simulations to illustrate the effect of seasonal succession on the disease-free periodic solution and the persistence of the disease. These are shown in Fig. 1(d) and Fig. 3, respectively. Numerical simulations show that regardless of whether the disease is extinct or persistent, ignoring the effects of seasonal transitions, especially the cycle alternation of the dry and rainy seasons, on the total population of mosquito populations and their behavior inevitably overestimates or underestimates the disease trend. Therefore, the establishment of a reasonable seasonal mosquito-borne disease model to accurately predict the epidemic trend of the disease has very important guiding significance for the reasonable allocation and use of limited medical resources.