# American Institute of Mathematical Sciences

## Persistence of mosquito vector and dengue: Impact of seasonal and diurnal temperature variations

 1 Department of Mathematics and Statistics, San Diego State University, San Diego, CA 92182, USA, Computational Science Research Center, San Diego State University, San Diego, CA 92182, USA, Viral Information Institute, San Diego State University, San Diego, CA 92182, USA 2 Department of Natural Science in the Center for General Education, Chang Gung University, Guishan, Taoyuan 333, Taiwan, Community Medicine Research Center, Chang Gung Memorial Hospital, Keelung Branch, Keelung 204, Taiwan

Received  August 2020 Revised  January 2021 Published  February 2021

Dengue, a mosquito-borne disease, poses a tremendous burden to human health with about 390 million annual dengue infections worldwide. The environmental temperature plays a major role in the mosquito life-cycle as well as the mosquito-human-mosquito dengue transmission cycle. While previous studies have provided useful insights into the understanding of dengue diseases, there is little emphasis put on the role of environmental temperature variation, especially diurnal variation, in the mosquito vector and dengue dynamics. In this study, we develop a mathematical model to investigate the impact of seasonal and diurnal temperature variations on the persistence of mosquito vector and dengue. Importantly, using a threshold dynamical system approach to our model, we formulate the mosquito reproduction number and the infection invasion threshold, which completely determine the global threshold dynamics of mosquito population and dengue transmission, respectively. Our model predicts that both seasonal and diurnal variations of the environmental temperature can be determinant factors for the persistence of mosquito vector and dengue. In general, our numerical estimates of the mosquito reproduction number and the infection invasion threshold show that places with higher diurnal or seasonal temperature variations have a tendency to suffer less from the burden of mosquito population and dengue epidemics. Our results provide novel insights into the theoretical understanding of the role of diurnal temperature, which can be beneficial for the control of mosquito vector and dengue spread.

Citation: Naveen K. Vaidya, Feng-Bin Wang. Persistence of mosquito vector and dengue: Impact of seasonal and diurnal temperature variations. Discrete & Continuous Dynamical Systems - B, doi: 10.3934/dcdsb.2021048
##### References:

show all references

##### References:
A schematic diagram of the dengue transmission model. $A, M_s, M_e, M_i$: aquatic, susceptible, exposed, and infected mosquitos. $H_s, H_e, H_i, H_r$: susceptible, exposed, infected, and recovered humans. Solid arrows represent birth, maturation, infection, transfer, death, while dashed arrows indicate the effects of time dependent periodic environmental temperature, $T(t)$
Best-fit curves provided by the experimental data [44] for $\delta (T)$ (oviposition rate), $\mu_a(T)$ (aquatic phase mortality rate), $\theta (T)$ (mosquito emergence rate from acuatic phase), and $\mu_m(T)$ (mosquito mortality rate)
Best-fit curve provided by the data generated from the previous estimates [17] for the transmission probability from human to mosquito, $\beta_m (T)$
Temperature profile of $T_m(t)$ [Left], $\Psi_d(t)$ [Middle], and $T(t)$ [Right]. Parameters used are $T_0 = 25 ^oC, \epsilon_m = 5 ^oC, \tau_m$ = 365 day, $\phi_m = 0$, $\epsilon_d = 5 ^oC$, $\tau_d = 1$ day, and $\phi_d = 0$
Mosquito reproduction number ($\mathcal{R}^M$) [Left] and infection invasion threshold ($\mathcal{R}^0$) [right] for different values of the mean temperature ($T_0$) with amplitudes of seasonal temperature and diurnal temperature fixed at $\epsilon_m = 5$ $^oC$ and $\epsilon_d = 5$ $^oC$, respectively. For comparison purposes, $\mathcal{R}^M$ and $\mathcal{R}^0$ for the constant temperature (i.e., $\epsilon_m = \epsilon_d = 0$ $^oC$) are also plotted
Mosquito reproduction number ($\mathcal{R}^M$) [left column] and infection invasion threshold ($\mathcal{R}^0$) [right column] for different values of the amplitudes of seasonal temperature ($\epsilon_m$) with the amplitude of diurnal temperature fixed at $\epsilon_d = 5$ $^oC$ and the mean temperature fixed at $T_0 = 16$ $^oC$ [top row], $T_0 = 28$ $^oC$ [middle row], and $T_0 = 38$ $^oC$ [bottom row]
Mosquito reproduction number ($\mathcal{R}^M$) [left column] and infection invasion threshold ($\mathcal{R}^0$) [right column] for different values of the amplitudes of diurnal temperature ($\epsilon_d$) with the amplitude of seasonal temperature fixed at $\epsilon_m = 5$ $^oC$ and the mean temperature fixed at $T_0 = 16$ $^oC$ [top row], $T_0 = 28$ $^oC$ [middle row], and $T_0 = 38$ $^oC$ [bottom row]
Model parameters
 Parameter Description Value Reference $k$ Fraction of female larvae from eggs 0.5 (0-1) [18,26] $b$ Per capita biting rate 0.1 [6,26] $\mu_h$ Natural death rate of humans 4.22$\times 10^{-5}$ d$^{-1}$ Calculated, [16] $1/\gamma_h$ Intrinsic period 10 days [6,16,18,26] $\alpha_h$ Human recovery rate 0.1 d$^{-1}$ [18,26] $\delta_m$ In $\delta(t)$ 9.531 Data fitting $\delta_h$ In $\delta(t)$ 22.55 Data fitting $N_{\delta}$ In $\delta(t)$ 7.084 Data fitting $a_{\delta}$ In $\delta(t)$ 0 Data fitting $\epsilon_{\delta}$ In $\delta(t)$ $10^{-6}$ Data fitting $a_{0\mu_a}$ In $\mu_a(t)$ 2.914 Data fitting $a_{1\mu_a}$ In $\mu_a(t)$ -0.4986 Data fitting $a_{2\mu_a}$ In $\mu_a(t)$ 0.03099 Data fitting $a_{3\mu_a}$ In $\mu_a(t)$ -0.0008236 Data fitting $a_{4\mu_a}$ In $\mu_a(t)$ 7.975$\times 10^{-6}$ Data fitting $a_{0\theta}$ In $\theta(t)$ 8.044$\times 10^{-5}$ Data fitting $a_{1\theta}$ In $\theta(t)$ 11.386 Data fitting $a_{2\theta}$ In $\theta(t)$ 40.1461 Data fitting $\epsilon_{\theta}$ In $\theta(t)$ $10^{-6}$ Data fitting $a_{0\mu_m}$ In $\mu_m(t)$ 0.1901 Data fitting $a_{1\mu_m}$ In $\mu_m(t)$ -0.0134 Data fitting $a_{2\mu_m}$ In $\mu_m(t)$ 2.739$\times 10^{-4}$ Data fitting $a_{0\gamma_m}$ In $\gamma_m(t)$ 5$\times 10^{4/3}$ Data fitting $a_{1\gamma_m}$ In $\gamma_m(t)$ 0.0768 Data fitting $\beta_{mh}$ In $\beta_m(t)$ 18.9871 Data fitting $N_{\beta_m}$ In $\beta_m(t)$ 7 Data fitting $\epsilon_{\beta m}$ In $\beta_m(t)$ $10^{-6}$ Data fitting $a_{\beta m}$ In $\beta_m(t)$ 0 Data fitting $a_{0\beta_h}$ In $\beta_h(t)$ 1.044$\times 10^{-3}$ Data fitting $a_{1\beta_h}$ In $\beta_h(t)$ 12.286 Data fitting $a_{2\beta_h}$ In $\beta_h(t)$ 32.461 Data fitting $\epsilon_{\beta h}$ In $\beta_h(t)$ $10^{-6}$ Data fitting
 Parameter Description Value Reference $k$ Fraction of female larvae from eggs 0.5 (0-1) [18,26] $b$ Per capita biting rate 0.1 [6,26] $\mu_h$ Natural death rate of humans 4.22$\times 10^{-5}$ d$^{-1}$ Calculated, [16] $1/\gamma_h$ Intrinsic period 10 days [6,16,18,26] $\alpha_h$ Human recovery rate 0.1 d$^{-1}$ [18,26] $\delta_m$ In $\delta(t)$ 9.531 Data fitting $\delta_h$ In $\delta(t)$ 22.55 Data fitting $N_{\delta}$ In $\delta(t)$ 7.084 Data fitting $a_{\delta}$ In $\delta(t)$ 0 Data fitting $\epsilon_{\delta}$ In $\delta(t)$ $10^{-6}$ Data fitting $a_{0\mu_a}$ In $\mu_a(t)$ 2.914 Data fitting $a_{1\mu_a}$ In $\mu_a(t)$ -0.4986 Data fitting $a_{2\mu_a}$ In $\mu_a(t)$ 0.03099 Data fitting $a_{3\mu_a}$ In $\mu_a(t)$ -0.0008236 Data fitting $a_{4\mu_a}$ In $\mu_a(t)$ 7.975$\times 10^{-6}$ Data fitting $a_{0\theta}$ In $\theta(t)$ 8.044$\times 10^{-5}$ Data fitting $a_{1\theta}$ In $\theta(t)$ 11.386 Data fitting $a_{2\theta}$ In $\theta(t)$ 40.1461 Data fitting $\epsilon_{\theta}$ In $\theta(t)$ $10^{-6}$ Data fitting $a_{0\mu_m}$ In $\mu_m(t)$ 0.1901 Data fitting $a_{1\mu_m}$ In $\mu_m(t)$ -0.0134 Data fitting $a_{2\mu_m}$ In $\mu_m(t)$ 2.739$\times 10^{-4}$ Data fitting $a_{0\gamma_m}$ In $\gamma_m(t)$ 5$\times 10^{4/3}$ Data fitting $a_{1\gamma_m}$ In $\gamma_m(t)$ 0.0768 Data fitting $\beta_{mh}$ In $\beta_m(t)$ 18.9871 Data fitting $N_{\beta_m}$ In $\beta_m(t)$ 7 Data fitting $\epsilon_{\beta m}$ In $\beta_m(t)$ $10^{-6}$ Data fitting $a_{\beta m}$ In $\beta_m(t)$ 0 Data fitting $a_{0\beta_h}$ In $\beta_h(t)$ 1.044$\times 10^{-3}$ Data fitting $a_{1\beta_h}$ In $\beta_h(t)$ 12.286 Data fitting $a_{2\beta_h}$ In $\beta_h(t)$ 32.461 Data fitting $\epsilon_{\beta h}$ In $\beta_h(t)$ $10^{-6}$ Data fitting
 [1] Irena PawŃow, Wojciech M. Zajączkowski. Global regular solutions to three-dimensional thermo-visco-elasticity with nonlinear temperature-dependent specific heat. Communications on Pure & Applied Analysis, 2017, 16 (4) : 1331-1372. doi: 10.3934/cpaa.2017065 [2] Rafael Luís, Sandra Mendonça. A note on global stability in the periodic logistic map. Discrete & Continuous Dynamical Systems - B, 2020, 25 (11) : 4211-4220. doi: 10.3934/dcdsb.2020094 [3] Lakmi Niwanthi Wadippuli, Ivan Gudoshnikov, Oleg Makarenkov. Global asymptotic stability of nonconvex sweeping processes. Discrete & Continuous Dynamical Systems - B, 2020, 25 (3) : 1129-1139. doi: 10.3934/dcdsb.2019212 [4] Michael Grinfeld, Amy Novick-Cohen. Some remarks on stability for a phase field model with memory. Discrete & Continuous Dynamical Systems - A, 2006, 15 (4) : 1089-1117. doi: 10.3934/dcds.2006.15.1089 [5] Carlos Gutierrez, Nguyen Van Chau. A remark on an eigenvalue condition for the global injectivity of differentiable maps of $R^2$. Discrete & Continuous Dynamical Systems - A, 2007, 17 (2) : 397-402. doi: 10.3934/dcds.2007.17.397 [6] Bernold Fiedler, Carlos Rocha, Matthias Wolfrum. Sturm global attractors for $S^1$-equivariant parabolic equations. Networks & Heterogeneous Media, 2012, 7 (4) : 617-659. doi: 10.3934/nhm.2012.7.617 [7] Cicely K. Macnamara, Mark A. J. Chaplain. Spatio-temporal models of synthetic genetic oscillators. Mathematical Biosciences & Engineering, 2017, 14 (1) : 249-262. doi: 10.3934/mbe.2017016 [8] Fernando P. da Costa, João T. Pinto, Rafael Sasportes. On the convergence to critical scaling profiles in submonolayer deposition models. Kinetic & Related Models, 2018, 11 (6) : 1359-1376. doi: 10.3934/krm.2018053 [9] Jian Yang, Bendong Lou. Traveling wave solutions of competitive models with free boundaries. Discrete & Continuous Dynamical Systems - B, 2014, 19 (3) : 817-826. doi: 10.3934/dcdsb.2014.19.817 [10] V. V. Zhikov, S. E. Pastukhova. Korn inequalities on thin periodic structures. Networks & Heterogeneous Media, 2009, 4 (1) : 153-175. doi: 10.3934/nhm.2009.4.153 [11] Giovanni Cimatti. Forced periodic solutions for piezoelectric crystals. Communications on Pure & Applied Analysis, 2005, 4 (2) : 475-485. doi: 10.3934/cpaa.2005.4.475 [12] Rui Hu, Yuan Yuan. Stability, bifurcation analysis in a neural network model with delay and diffusion. Conference Publications, 2009, 2009 (Special) : 367-376. doi: 10.3934/proc.2009.2009.367 [13] Gloria Paoli, Gianpaolo Piscitelli, Rossanno Sannipoli. A stability result for the Steklov Laplacian Eigenvalue Problem with a spherical obstacle. Communications on Pure & Applied Analysis, 2021, 20 (1) : 145-158. doi: 10.3934/cpaa.2020261 [14] Shanshan Chen, Junping Shi, Guohong Zhang. Spatial pattern formation in activator-inhibitor models with nonlocal dispersal. Discrete & Continuous Dynamical Systems - B, 2021, 26 (4) : 1843-1866. doi: 10.3934/dcdsb.2020042 [15] Alina Chertock, Alexander Kurganov, Mária Lukáčová-Medvi${\rm{\check{d}}}$ová, Șeyma Nur Özcan. An asymptotic preserving scheme for kinetic chemotaxis models in two space dimensions. Kinetic & Related Models, 2019, 12 (1) : 195-216. doi: 10.3934/krm.2019009 [16] Jaume Llibre, Luci Any Roberto. On the periodic solutions of a class of Duffing differential equations. Discrete & Continuous Dynamical Systems - A, 2013, 33 (1) : 277-282. doi: 10.3934/dcds.2013.33.277 [17] Zhi-Min Chen, Philip A. Wilson. Stability of oscillatory gravity wave trains with energy dissipation and Benjamin-Feir instability. Discrete & Continuous Dynamical Systems - B, 2012, 17 (7) : 2329-2341. doi: 10.3934/dcdsb.2012.17.2329 [18] Cécile Carrère, Grégoire Nadin. Influence of mutations in phenotypically-structured populations in time periodic environment. Discrete & Continuous Dynamical Systems - B, 2020, 25 (9) : 3609-3630. doi: 10.3934/dcdsb.2020075 [19] Marcelo Messias. Periodic perturbation of quadratic systems with two infinite heteroclinic cycles. Discrete & Continuous Dynamical Systems - A, 2012, 32 (5) : 1881-1899. doi: 10.3934/dcds.2012.32.1881 [20] Wei-Jian Bo, Guo Lin, Shigui Ruan. Traveling wave solutions for time periodic reaction-diffusion systems. Discrete & Continuous Dynamical Systems - A, 2018, 38 (9) : 4329-4351. doi: 10.3934/dcds.2018189

2019 Impact Factor: 1.27