American Institute of Mathematical Sciences

March  2018, 23(2): 609-628. doi: 10.3934/dcdsb.2018036

Stability of travelling waves in a Wolbachia invasion

 School of Mathematics and Statistics, University of Sydney, NSW 2006, Australia

* Corresponding author: Matthew H. Chan

Received  June 2016 Revised  September 2017 Published  December 2017

Fund Project: PSK was supported by the Australian Research Council, Discovery Project (DP160101597).

Numerous studies have examined the growth dynamics of Wolbachia within populations and the resultant rate of spatial spread. This spread is typically characterised as a travelling wave with bistable local growth dynamics due to a strong Allee effect generated from cytoplasmic incompatibility. While this rate of spread has been calculated from numerical solutions of reaction-diffusion models, none have examined the spectral stability of such travelling wave solutions. In this study we analyse the stability of a travelling wave solution generated by the reaction-diffusion model of Chan & Kim [4] by computing the essential and point spectrum of the linearised operator arising in the model. The point spectrum is computed via an Evans function using the compound matrix method, whereby we find that it has no roots with positive real part. Moreover, the essential spectrum lies strictly in the left half plane. Thus, we find that the travelling wave solution found by Chan & Kim [4] corresponding to competition between Wolbachia-infected and -uninfected mosquitoes is linearly stable. We employ a dimension counting argument to suggest that, under realistic conditions, the wavespeed corresponding to such a solution is unique.

Citation: Matthew H. Chan, Peter S. Kim, Robert Marangell. Stability of travelling waves in a Wolbachia invasion. Discrete & Continuous Dynamical Systems - B, 2018, 23 (2) : 609-628. doi: 10.3934/dcdsb.2018036
References:

show all references

References:
The essential spectrum of $\mathcal{L}$ is given by $\lambda$ in the shaded regions. The blue dashed and solid lines represent $\lambda_-^{1, 2}$ and $\lambda_+^{1, 2}$ respectively. The red line indicates the absolute spectrum given by Eq. (17) and the red dot at the origin represents an eigenvalue. Note that these are not drawn to scale for visualisation purposes
The contour $C$
The image of $C$ under $D(\lambda)$, where $r_s =0.1$ and $r_b =10$
The image of $C$ under $D(\lambda)$, where $r_s =0.001$ and $r_b =500$
Plots $(a)$ and $(b)$ show the change in argument for $D[C]$ corresponding to Figures 3 and 4 respectively
Plot of the Evans function given by Eq. (26). The blue and red solid lines in both plots show $D(\lambda)$ with $\alpha = 1.1$ and $\alpha = 1$ respectively. The dashed lines in plot $(a)$ mark the edge of the absolute spectrum corresponding to each $\alpha$. For $\alpha = 1.1$, the only roots are at $\lambda =0$ and at the edge of the absolute spectrum $\lambda =-0.002607$, whereas for $\alpha = 1$ we were unable to detect a zero at the edge of the absolute spectrum due to its proximity to the origin
Solution to the boundary value problem (6). Figure $(a)$ shows the wave profile of $\hat{u}(z)$ and $\hat{v}(z)$, represented by solid and dashed lines respectively. Figure $(b)$ shows the heteroclinic connection between equilibrium states $\mathbf{e_{-}} = (1-\alpha \mu, 0)$ and $\mathbf{e_{+}} = (0, 1-\frac{\mu}{F})$, where the solid and dashed line represent the solution in $u-u'$ and $v-v'$ space respectively
A diagram showing the uniqueness of $c_*$ by dimension counting, where $\tilde{\mathbf{e}}_-$, $\tilde{\mathbf{e}}_+$ denote the equilibria at the end points of the heteroclinic orbit in (33). In the illustration, $W^{u, s}(\tilde{\mathbf{e}}_{-, +}) \times \mathcal{C}_\varepsilon$ are shown as 2-dimensional manifolds which intersect transversally in 3-dimensional space. The one dimensional intersection corresponds to the heteroclinic connection between $\tilde{\mathbf{e}}_-$ and $\tilde{\mathbf{e}}_+$ at $c_*$
Simulations corresponding to parameter set three, listed in Table 2
Simulations corresponding to parameter set two, listed in Table 3
Simulations corresponding to parameter set three, listed in Table 4
Simulations corresponding to parameter set four, listed in Table 5
Parameter values and definitions
 Symbol Definition Value $F$ Relative fecundity of uninfected to infected females $1.0526$ $s_h$ Probability of embryo death due to CI $0.45$ $\mu$ Mortality rate $0.0162$ $\alpha$ Reduction in lifespan due to infection $1.1$
 Symbol Definition Value $F$ Relative fecundity of uninfected to infected females $1.0526$ $s_h$ Probability of embryo death due to CI $0.45$ $\mu$ Mortality rate $0.0162$ $\alpha$ Reduction in lifespan due to infection $1.1$
Parameter set one
 Symbol Definition Value $F$ Relative fecundity of uninfected to infected females $1.05$ $s_h$ Probability of embryo death due to CI $0.7$ $\mu$ Mortality rate $0.03$ $\alpha$ Reduction in lifespan due to infection $1.2$
 Symbol Definition Value $F$ Relative fecundity of uninfected to infected females $1.05$ $s_h$ Probability of embryo death due to CI $0.7$ $\mu$ Mortality rate $0.03$ $\alpha$ Reduction in lifespan due to infection $1.2$
Parameter set two
 Symbol Definition Value $F$ Relative fecundity of uninfected to infected females $1.1$ $s_h$ Probability of embryo death due to CI $0.9$ $\mu$ Mortality rate $0.02$ $\alpha$ Reduction in lifespan due to infection $1.3$
 Symbol Definition Value $F$ Relative fecundity of uninfected to infected females $1.1$ $s_h$ Probability of embryo death due to CI $0.9$ $\mu$ Mortality rate $0.02$ $\alpha$ Reduction in lifespan due to infection $1.3$
Parameter set three
 Symbol Definition Value $F$ Relative fecundity of uninfected to infected females $1.1$ $s_h$ Probability of embryo death due to CI $0.8$ $\mu$ Mortality rate $0.05$ $\alpha$ Reduction in lifespan due to infection $1.2$
 Symbol Definition Value $F$ Relative fecundity of uninfected to infected females $1.1$ $s_h$ Probability of embryo death due to CI $0.8$ $\mu$ Mortality rate $0.05$ $\alpha$ Reduction in lifespan due to infection $1.2$
Parameter set four
 Symbol Definition Value $F$ Relative fecundity of uninfected to infected females $1.4$ $s_h$ Probability of embryo death due to CI $0.9$ $\mu$ Mortality rate $0.1$ $\alpha$ Reduction in lifespan due to infection $1.2$
 Symbol Definition Value $F$ Relative fecundity of uninfected to infected females $1.4$ $s_h$ Probability of embryo death due to CI $0.9$ $\mu$ Mortality rate $0.1$ $\alpha$ Reduction in lifespan due to infection $1.2$
 [1] A. M. Elaiw, N. H. AlShamrani, A. Abdel-Aty, H. Dutta. Stability analysis of a general HIV dynamics model with multi-stages of infected cells and two routes of infection. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020441 [2] Ebraheem O. Alzahrani, Muhammad Altaf Khan. Androgen driven evolutionary population dynamics in prostate cancer growth. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020426 [3] Laurent Di Menza, Virginie Joanne-Fabre. An age group model for the study of a population of trees. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020464 [4] Yifan Chen, Thomas Y. Hou. Function approximation via the subsampled Poincaré inequality. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 169-199. doi: 10.3934/dcds.2020296 [5] Min Chen, Olivier Goubet, Shenghao Li. Mathematical analysis of bump to bucket problem. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5567-5580. doi: 10.3934/cpaa.2020251 [6] Bahaaeldin Abdalla, Thabet Abdeljawad. Oscillation criteria for kernel function dependent fractional dynamic equations. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020443 [7] Reza Chaharpashlou, Abdon Atangana, Reza Saadati. On the fuzzy stability results for fractional stochastic Volterra integral equation. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020432 [8] Scipio Cuccagna, Masaya Maeda. A survey on asymptotic stability of ground states of nonlinear Schrödinger equations II. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020450 [9] Qianqian Han, Xiao-Song Yang. Qualitative analysis of a generalized Nosé-Hoover oscillator. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020346 [10] Laurence Cherfils, Stefania Gatti, Alain Miranville, Rémy Guillevin. Analysis of a model for tumor growth and lactate exchanges in a glioma. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020457 [11] Vieri Benci, Sunra Mosconi, Marco Squassina. Preface: Applications of mathematical analysis to problems in theoretical physics. Discrete & Continuous Dynamical Systems - S, 2020  doi: 10.3934/dcdss.2020446 [12] Weiwei Liu, Jinliang Wang, Yuming Chen. Threshold dynamics of a delayed nonlocal reaction-diffusion cholera model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020316 [13] Manil T. Mohan. First order necessary conditions of optimality for the two dimensional tidal dynamics system. Mathematical Control & Related Fields, 2020  doi: 10.3934/mcrf.2020045 [14] Cuicui Li, Lin Zhou, Zhidong Teng, Buyu Wen. The threshold dynamics of a discrete-time echinococcosis transmission model. Discrete & Continuous Dynamical Systems - B, 2020  doi: 10.3934/dcdsb.2020339 [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] Lingfeng Li, Shousheng Luo, Xue-Cheng Tai, Jiang Yang. A new variational approach based on level-set function for convex hull problem with outliers. Inverse Problems & Imaging, , () : -. doi: 10.3934/ipi.2020070 [17] Mohammed Abdulrazaq Kahya, Suhaib Abduljabbar Altamir, Zakariya Yahya Algamal. Improving whale optimization algorithm for feature selection with a time-varying transfer function. Numerical Algebra, Control & Optimization, 2021, 11 (1) : 87-98. doi: 10.3934/naco.2020017 [18] Yining Cao, Chuck Jia, Roger Temam, Joseph Tribbia. Mathematical analysis of a cloud resolving model including the ice microphysics. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 131-167. doi: 10.3934/dcds.2020219 [19] Xin Guo, Lei Shi. Preface of the special issue on analysis in data science: Methods and applications. Mathematical Foundations of Computing, 2020, 3 (4) : i-ii. doi: 10.3934/mfc.2020026 [20] Xin-Guang Yang, Lu Li, Xingjie Yan, Ling Ding. The structure and stability of pullback attractors for 3D Brinkman-Forchheimer equation with delay. Electronic Research Archive, 2020, 28 (4) : 1395-1418. doi: 10.3934/era.2020074

2019 Impact Factor: 1.27