American Institute of Mathematical Sciences

August  2018, 15(4): 1033-1054. doi: 10.3934/mbe.2018046

Role of white-tailed deer in geographic spread of the black-legged tick Ixodes scapularis : Analysis of a spatially nonlocal model

 1 Department of Mathematics, University of Surrey, Guildford, Surrey GU2 7XH, UK 2 Institute for Mathematical Sciences, Renmin University of China, Beijing, China 3 Department of Mathematics, College of William and Mary, Williamsburg, VA, USA 4 Key Laboratory of Eco-environments in Three Gorges Reservoir Region, and School of Mathematics and Statistics, Southwest University, Chongqing, China 5 Department of Mathematical Sciences, University of Cincinnati, Cincinnati, OH, USA 6 Department of Applied Mathematics, University of Western Ontario, London, ON, Canada

* Corresponding author: Xingfu Zou

Received  November 15, 2017 Accepted  November 15, 2017 Published  March 2018

Lyme disease is transmitted via blacklegged ticks, the spatial spread of which is believed to be primarily via transport on white-tailed deer. In this paper, we develop a mathematical model to describe the spatial spread of blacklegged ticks due to deer dispersal. The model turns out to be a system of differential equations with a spatially non-local term accounting for the phenomenon that a questing female adult tick that attaches to a deer at one location may later drop to the ground, fully fed, at another location having been transported by the deer. We first justify the well-posedness of the model and analyze the stability of its steady states. We then explore the existence of traveling wave fronts connecting the extinction equilibrium with the positive equilibrium for the system. We derive an algebraic equation that determines a critical value $c^*$ which is at least a lower bound for the wave speed in the sense that, if $c < c^*$, there is no traveling wave front of speed $c$ connecting the extinction steady state to the positive steady state. Numerical simulations of the wave equations suggest that $c^*$ is the minimum wave speed. We also carry out some numerical simulations for the original spatial model system and the results seem to confirm that the actual spread rate of the tick population coincides with $c^*$. We also explore the dependence of $c^*$ on the dispersion rate of the white tailed deer, by which one may evaluate the role of the deer's dispersion in the geographical spread of the ticks.

Citation: Stephen A. Gourley, Xiulan Lai, Junping Shi, Wendi Wang, Yanyu Xiao, Xingfu Zou. Role of white-tailed deer in geographic spread of the black-legged tick Ixodes scapularis : Analysis of a spatially nonlocal model. Mathematical Biosciences & Engineering, 2018, 15 (4) : 1033-1054. doi: 10.3934/mbe.2018046
References:

show all references

References:
The life-stage components of the model: questing larvae ($L$) find a host, feed and moult into questing nymphs ($N$), which then find a new host, feed and moult into questing adults ($A_q$). Adult females that find a deer host ($A_f$) feed, drop to the forest floor, lay $2000$ eggs and then die. Hatching eggs create the next generation of questing larvae. The $r$ parameters are the per-capita transition rates between each compartment
$H_1(\lambda, c)$ and $H_2(\lambda, c)$ for different $c$. (a) $c = 0.55$; (b) $c = c^* = 0.6176844021$, ($\lambda = 0.7081234538$); (c) $c = 0.7$. Here, the model parameters are taken as $b = 3000$, $r_1 = 0.13$, $r_2 = 0.13$, $r_3 = 0.03$, $r_4 = 0.03$, $d_1 = 0.3$, $d_2 = 0.3$, $d_3 = 0.1$, $d_4 = 0.1$, $\tau_1 = 20$, $\tau_2 = 10$, $N_{cap} = 5000$, $h = 100$ and $D = 1$
Dependence of $c^*$ on $b$ and $D$ respectively: (a) with $D = 1$; (b) with $b = 3000$. Other parameters are taken as: $r_1 = 0.13$, $r_2 = 0.13$, $r_3 = 0.03$, $r_4 = 0.03$, $d_1 = 0.3$, $d_2 = 0.3$, $d_3 = 0.1$, $d_4 = 0.1$, $\tau_1 = 20$, $\tau_2 = 10$, $N_{cap} = 5000$ and $h = 100$
There is no biologically relevant traveling wave front solution with speed $c = 0.1<c^* = 0.24$: $\phi_1$ may take negative values
There is a non-negative traveling wave front solution with speed $c = 0.4>c^* = 0.24$
(a): time evolution of $L(x, t)$; (b): time evolution of $N(x, t)$; (c): contours of (a) with region where $L(x, t)>0.1$ shown in grey; (d): contours of (b) with region where $N(x, t)>0.1$ shown in grey
(a): time evolution of $A_q(x, t)$; (b): time evolution of $A_f(x, t)$; (c): contours of (a) with region where $A_q(x, t)>0.1$ shown in grey; (d): contours of (b) with region where $A_f(x, t)>0.1$ shown in grey
Explanation of parameters
 Parameters Meaning Value $b$ Birth rate of tick $3000$ $1/r_1$ average time that a questing larvae needs to feed and moult $1/0.13$ $1/r_2$ average time that a questing nymph needs to feed and moult $1/0.13$ $1/r_3$ average time that a questing adult needs to successfully attach to a deer $1/0.03$ $r_4$ Proportion of fed adults that can lay eggs 0.03 $d_1$ per-capita death rate of larvae 0.3 $d_2$ per-capita death rate of nymphs 0.3 $d_3$ per-capita death rate of questing adults 0.1 $d_4$ per-capita death rate of fed adults 0.1 $\tau_1$ average time between last blood feeding and hatch of laid eggs 20 days $\tau_2$ average time tick is attached to a deer $10$ days
 Parameters Meaning Value $b$ Birth rate of tick $3000$ $1/r_1$ average time that a questing larvae needs to feed and moult $1/0.13$ $1/r_2$ average time that a questing nymph needs to feed and moult $1/0.13$ $1/r_3$ average time that a questing adult needs to successfully attach to a deer $1/0.03$ $r_4$ Proportion of fed adults that can lay eggs 0.03 $d_1$ per-capita death rate of larvae 0.3 $d_2$ per-capita death rate of nymphs 0.3 $d_3$ per-capita death rate of questing adults 0.1 $d_4$ per-capita death rate of fed adults 0.1 $\tau_1$ average time between last blood feeding and hatch of laid eggs 20 days $\tau_2$ average time tick is attached to a deer $10$ days
 [1] Mahir Demir, Suzanne Lenhart. A spatial food chain model for the Black Sea Anchovy, and its optimal fishery. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 155-171. doi: 10.3934/dcdsb.2020373 [2] Yu Jin, Xiang-Qiang Zhao. The spatial dynamics of a Zebra mussel model in river environments. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020362 [3] Ting Liu, Guo-Bao Zhang. Global stability of traveling waves for a spatially discrete diffusion system with time delay. Electronic Research Archive, , () : -. doi: 10.3934/era.2021003 [4] Simone Göttlich, Elisa Iacomini, Thomas Jung. Properties of the LWR model with time delay. Networks & Heterogeneous Media, 2020  doi: 10.3934/nhm.2020032 [5] Ran Zhang, Shengqiang Liu. On the asymptotic behaviour of traveling wave solution for a discrete diffusive epidemic model. Discrete & Continuous Dynamical Systems - B, 2021, 26 (2) : 1197-1204. doi: 10.3934/dcdsb.2020159 [6] Mohammad Ghani, Jingyu Li, Kaijun Zhang. Asymptotic stability of traveling fronts to a chemotaxis model with nonlinear diffusion. Discrete & Continuous Dynamical Systems - B, 2021  doi: 10.3934/dcdsb.2021017 [7] Yoichi Enatsu, Emiko Ishiwata, Takeo Ushijima. Traveling wave solution for a diffusive simple epidemic model with a free boundary. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 835-850. doi: 10.3934/dcdss.2020387 [8] Mugen Huang, Moxun Tang, Jianshe Yu, Bo Zheng. A stage structured model of delay differential equations for Aedes mosquito population suppression. Discrete & Continuous Dynamical Systems - A, 2020, 40 (6) : 3467-3484. doi: 10.3934/dcds.2020042 [9] Guihong Fan, Gail S. K. Wolkowicz. Chaotic dynamics in a simple predator-prey model with discrete delay. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 191-216. doi: 10.3934/dcdsb.2020263 [10] Hai-Feng Huo, Shi-Ke Hu, Hong Xiang. Traveling wave solution for a diffusion SEIR epidemic model with self-protection and treatment. Electronic Research Archive, , () : -. doi: 10.3934/era.2020118 [11] Tin Phan, Bruce Pell, Amy E. Kendig, Elizabeth T. Borer, Yang Kuang. Rich dynamics of a simple delay host-pathogen model of cell-to-cell infection for plant virus. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 515-539. doi: 10.3934/dcdsb.2020261 [12] Thazin Aye, Guanyu Shang, Ying Su. On a stage-structured population model in discrete periodic habitat: III. unimodal growth and delay effect. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2021005 [13] Kengo Nakai, Yoshitaka Saiki. Machine-learning construction of a model for a macroscopic fluid variable using the delay-coordinate of a scalar observable. Discrete & Continuous Dynamical Systems - S, 2021, 14 (3) : 1079-1092. doi: 10.3934/dcdss.2020352 [14] Riccarda Rossi, Ulisse Stefanelli, Marita Thomas. Rate-independent evolution of sets. Discrete & Continuous Dynamical Systems - S, 2021, 14 (1) : 89-119. doi: 10.3934/dcdss.2020304 [15] Shao-Xia Qiao, Li-Jun Du. Propagation dynamics of nonlocal dispersal equations with inhomogeneous bistable nonlinearity. Electronic Research Archive, , () : -. doi: 10.3934/era.2020116 [16] Yong-Jung Kim, Hyowon Seo, Changwook Yoon. Asymmetric dispersal and evolutional selection in two-patch system. Discrete & Continuous Dynamical Systems - A, 2020, 40 (6) : 3571-3593. doi: 10.3934/dcds.2020043 [17] Nabahats Dib-Baghdadli, Rabah Labbas, Tewfik Mahdjoub, Ahmed Medeghri. On some reaction-diffusion equations generated by non-domiciliated triatominae, vectors of Chagas disease. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2021004 [18] Shang Wu, Pengfei Xu, Jianhua Huang, Wei Yan. Ergodicity of stochastic damped Ostrovsky equation driven by white noise. Discrete & Continuous Dynamical Systems - B, 2021, 26 (3) : 1615-1626. doi: 10.3934/dcdsb.2020175 [19] Jian Zhang, Tony T. Lee, Tong Ye, Liang Huang. An approximate mean queue length formula for queueing systems with varying service rate. Journal of Industrial & Management Optimization, 2021, 17 (1) : 185-204. doi: 10.3934/jimo.2019106 [20] Patrick Martinez, Judith Vancostenoble. Lipschitz stability for the growth rate coefficients in a nonlinear Fisher-KPP equation. Discrete & Continuous Dynamical Systems - S, 2021, 14 (2) : 695-721. doi: 10.3934/dcdss.2020362

2018 Impact Factor: 1.313

Tools

Article outline

Figures and Tables