# American Institute of Mathematical Sciences

August  2018, 15(4): 841-862. doi: 10.3934/mbe.2018038

## Continuous and pulsed epidemiological models for onchocerciasis with implications for eradication strategy

 1 Department of Mathematics, University of Nebraska-Lincoln, Lincoln, NE 68588-0130, USA 2 Mathematics Department, Seattle University, 901 12th Ave, Seattle, WA 98122-1090, USA 3 Department of Mathematics, Indiana University of Pennsylvania, Indiana, PA 15705, USA 4 Department of Mathematics, New York City College of Technology -CUNY, Brooklyn, NY 11201, USA

* Corresponding author: Glenn Ledder

Received  March 2017 Accepted  January 03, 2018 Published  March 2018

Fund Project: This work was supported by NSF grant DMS 1239280.

Onchocerciasis is an endemic disease in parts of sub-Saharan Africa. Complex mathematical models are being used to assess the likely efficacy of efforts to eradicate the disease; however, their predictions have not always been borne out in practice. In this paper, we represent the immunological aspects of the disease with a single empirical parameter in order to reduce the model complexity. Asymptotic approximation allows us to reduce the vector-borne epidemiological model to a model of an infectious disease with nonlinear incidence. We then consider two versions, one with continuous treatment and a more realistic one where treatment occurs only at intervals. Thorough mathematical analysis of these models yields equilibrium solutions for the continuous case, periodic solutions for the pulsed case, and conditions for the existence of endemic disease equilibria in both cases, thereby leading to simple model criteria for eradication. The analytical results and numerical experiments show that the continuous treatment version is an excellent approximation for the pulsed version and that the current onchocerciasis eradication strategy is inadequate for regions where the incidence is highest and unacceptably slow even when the long-term behavior is the disease-free state.

Citation: Glenn Ledder, Donna Sylvester, Rachelle R. Bouchat, Johann A. Thiel. Continuous and pulsed epidemiological models for onchocerciasis with implications for eradication strategy. Mathematical Biosciences & Engineering, 2018, 15 (4) : 841-862. doi: 10.3934/mbe.2018038
##### References:

show all references

##### References:
The O. volvulus life cycle. Microfilariae produced in the human host are transmitted to the black fly of the genus Simulium via a bite. Within the black fly, the larvae pass through larval stages L1-L3. At larval stage L3, they are transmitted to a human host via a bite, where they pass through the final larval stages L3-L5 and become adults [8]
The SEIPMS-UV epidemiological model. Fly bites transfer the disease from infected flies ($V$) to susceptible humans ($S$) and from infective humans ($I$, $P$, and $M$) to uninfected flies ($U$), with the transfer from medicated humans ($M$) decreased by a factor $\nu$. Exposed humans ($E$) become infective, with the fraction $p$ counting as premedicated ($P$) while waiting for treatment and the remaining fraction $q = 1-p$ becoming unmedicated infectives ($I$). Premedicated humans become medicated when they receive treatment. All three infective classes can become susceptible by clearing all the adult worms. Birth and death rates for humans are equal, with all births into the susceptible class; similarly, birth and death rates for flies are equal with all births into the uninfected class
Simulation of the introduction of a small population of human infectives into a previously unexposed population, showing $E/N$, $V/F$, $I/N$, and $S/N$ in Figure 3a, from bottom up, and $V/F$, $I/N$ in Figure 3b, from bottom up, using $a = 0.9$, $b = 3.0$, $\eta = 0.1$, $\epsilon = 0$, $p = 0$, with solid for $\delta = 0.01$ and dash-dot for the asymptotic simplification
Periodic solutions for the exposed ($x$, dashed) and total infective ($y = h+i$, solid) classes, with treatment intervals of 2 years (top), 1 year (middle), and 6 months (bottom), using $a = 0.9$, $b = 3.5$, $\eta = 0.1$, $\epsilon = 0$, $\nu p = 0.6$
Time average infective populations, with $a = 0.9$, $b = 3$, $\eta = 0.1$, $\epsilon = 0$, $\nu p = 0.6$. Humans: top 2; Flies: bottom 2; Pulsed: solid; Continuous: dashed
Simulations of various treatment scenarios, with $a = 0.9$, $\eta = 0.1$, $\epsilon = 0$
 [1] Guo-Niu Han, Huan Xiong. Skew doubled shifted plane partitions: Calculus and asymptotics. Electronic Research Archive, 2021, 29 (1) : 1841-1857. doi: 10.3934/era.2020094 [2] Mauricio Achigar. Extensions of expansive dynamical systems. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020399 [3] Serena Dipierro, Benedetta Pellacci, Enrico Valdinoci, Gianmaria Verzini. Time-fractional equations with reaction terms: Fundamental solutions and asymptotics. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 257-275. doi: 10.3934/dcds.2020137 [4] Joel Kübler, Tobias Weth. Spectral asymptotics of radial solutions and nonradial bifurcation for the Hénon equation. Discrete & Continuous Dynamical Systems - A, 2020, 40 (6) : 3629-3656. doi: 10.3934/dcds.2020032 [5] The Editors. The 2019 Michael Brin Prize in Dynamical Systems. Journal of Modern Dynamics, 2020, 16: 349-350. doi: 10.3934/jmd.2020013 [6] Nitha Niralda P C, Sunil Mathew. On properties of similarity boundary of attractors in product dynamical systems. Discrete & Continuous Dynamical Systems - S, 2021  doi: 10.3934/dcdss.2021004 [7] Isabeau Birindelli, Françoise Demengel, Fabiana Leoni. Boundary asymptotics of the ergodic functions associated with fully nonlinear operators through a Liouville type theorem. Discrete & Continuous Dynamical Systems - A, 2020  doi: 10.3934/dcds.2020395 [8] Jiahao Qiu, Jianjie Zhao. Maximal factors of order $d$ of dynamical cubespaces. Discrete & Continuous Dynamical Systems - A, 2021, 41 (2) : 601-620. doi: 10.3934/dcds.2020278 [9] Fanni M. Sélley. A self-consistent dynamical system with multiple absolutely continuous invariant measures. Journal of Computational Dynamics, 2021, 8 (1) : 9-32. doi: 10.3934/jcd.2021002 [10] Toshiko Ogiwara, Danielle Hilhorst, Hiroshi Matano. Convergence and structure theorems for order-preserving dynamical systems with mass conservation. Discrete & Continuous Dynamical Systems - A, 2020, 40 (6) : 3883-3907. doi: 10.3934/dcds.2020129 [11] Peter Giesl, Zachary Langhorne, Carlos Argáez, Sigurdur Hafstein. Computing complete Lyapunov functions for discrete-time dynamical systems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 299-336. doi: 10.3934/dcdsb.2020331 [12] Alessandro Fonda, Rodica Toader. A dynamical approach to lower and upper solutions for planar systems "To the memory of Massimo Tarallo". Discrete & Continuous Dynamical Systems - A, 2021  doi: 10.3934/dcds.2021012 [13] Mikhail I. Belishev, Sergey A. Simonov. A canonical model of the one-dimensional dynamical Dirac system with boundary control. Evolution Equations & Control Theory, 2021  doi: 10.3934/eect.2021003 [14] Stefan Siegmund, Petr Stehlík. Time scale-induced asynchronous discrete dynamical systems. Discrete & Continuous Dynamical Systems - B, 2021, 26 (2) : 1011-1029. doi: 10.3934/dcdsb.2020151 [15] Anna Abbatiello, Eduard Feireisl, Antoní Novotný. Generalized solutions to models of compressible viscous fluids. Discrete & Continuous Dynamical Systems - A, 2021, 41 (1) : 1-28. doi: 10.3934/dcds.2020345 [16] Xin Guo, Lexin Li, Qiang Wu. Modeling interactive components by coordinate kernel polynomial models. Mathematical Foundations of Computing, 2020, 3 (4) : 263-277. doi: 10.3934/mfc.2020010 [17] Urszula Ledzewicz, Heinz Schättler. On the role of pharmacometrics in mathematical models for cancer treatments. Discrete & Continuous Dynamical Systems - B, 2021, 26 (1) : 483-499. doi: 10.3934/dcdsb.2020213 [18] P. K. Jha, R. Lipton. Finite element approximation of nonlocal dynamic fracture models. Discrete & Continuous Dynamical Systems - B, 2021, 26 (3) : 1675-1710. doi: 10.3934/dcdsb.2020178 [19] Wei Feng, Michael Freeze, Xin Lu. On competition models under allee effect: Asymptotic behavior and traveling waves. Communications on Pure & Applied Analysis, 2020, 19 (12) : 5609-5626. doi: 10.3934/cpaa.2020256 [20] Zhimin Li, Tailei Zhang, Xiuqing Li. Threshold dynamics of stochastic models with time delays: A case study for Yunnan, China. Electronic Research Archive, 2021, 29 (1) : 1661-1679. doi: 10.3934/era.2020085

2018 Impact Factor: 1.313