A hybrid Lagrangian-Eulerian model for vector-borne diseases.

In this paper, a multi-patch and multi-group vector-borne disease model is proposed to study the effects of host commuting (Lagrangian approach) and/or vector migration (Eulerian approach) on disease spread. We first define the basic reproduction number of the model, R 0 , which completely determines the global dynamics of the model system. Namely, if R 0 ≤ 1 , then the disease-free equilibrium is globally asymptotically stable, and if R 0 > 1 , then there exists a unique endemic equilibrium which is globally asymptotically stable. Then, we show that the basic reproduction number has lower and upper bounds which are independent of the host residence times matrix and the vector migration matrix. In particular, nonhomogeneous mixing of hosts and vectors in a homogeneous environment generally increases disease persistence and the basic reproduction number of the model attains its minimum when the distributions of hosts and vectors are proportional. Moreover, R 0 can also be estimated by the basic reproduction numbers of disconnected patches if the environment is homogeneous. The optimal vector control strategy is obtained for a special scenario. In the two-patch and two-group case, we numerically analyze the dependence of the basic reproduction number and the total number of infected people on the host residence times matrix and illustrate the optimal vector control strategy in homogeneous and heterogeneous environments.

Vector-borne disease models with Lagrangian approach.

We develop a multi-group and multi-patch model to study the effects of population dispersal on the spatial spread of vector-borne diseases across a heterogeneous environment. The movement of host and/or vector is described by Lagrangian approach in which the origin or identity of each individual stays unchanged regardless of movement. The basic reproduction number [Formula: see text] of the model is defined and the strong connectivity of the host-vector network is succinctly characterized by the residence times matrices of hosts and vectors. Furthermore, the definition and criterion of the strong connectivity of general infectious disease networks are given and applied to establish the global stability of the disease-free equilibrium. The global dynamics of the model system are shown to be entirely determined by its basic reproduction number. We then obtain several biologically meaningful upper and lower bounds on the basic reproduction number which are independent or dependent of the residence times matrices. In particular, the heterogeneous mixing of hosts and vectors in a homogeneous environment always increases the basic reproduction number. There is a substantial difference on the upper bound of [Formula: see text] between Lagrangian and Eulerian modeling approaches. When only host movement between two patches is concerned, the subdivision of hosts (more host groups) can lead to a larger basic reproduction number. In addition, we numerically investigate the dependence of the basic reproduction number and the total number of infected hosts on the residence times matrix of hosts, and compare the impact of different vector control strategies on disease transmission.

Relative prevalence-based dispersal in an epidemic patch model.

In this paper, we propose a two-patch SIRS model with a nonlinear incidence rate: [Formula: see text] and nonconstant dispersal rates, where the dispersal rates of susceptible and recovered individuals depend on the relative disease prevalence in two patches. In an isolated environment, the model admits Bogdanov-Takens bifurcation of codimension 3 (cusp case) and Hopf bifurcation of codimension up to 2 as the parameters vary, and exhibits rich dynamics such as multiple coexistent steady states and periodic orbits, homoclinic orbits and multitype bistability. The long-term dynamics can be classified in terms of the infection rates [Formula: see text] (due to single contact) and [Formula: see text] (due to double exposures). In a connected environment, we establish a threshold [Formula: see text] between disease extinction and uniform persistence under certain conditions. We numerically explore the effect of population dispersal on disease spread when [Formula: see text] and patch 1 has a lower infection rate, our results indicate: (i) [Formula: see text] can be nonmonotonic in dispersal rates and [Formula: see text] ([Formula: see text] is the basic reproduction number of patch i) may fail; (ii) the constant dispersal of susceptible individuals (or infective individuals) between two patches (or from patch 2 to patch 1) will increase (or reduce) the overall disease prevalence; (iii) the relative prevalence-based dispersal may reduce the overall disease prevalence. When [Formula: see text] and the disease outbreaks periodically in each isolated patch, we find that: (a) small unidirectional and constant dispersal can lead to complex periodic patterns like relaxation oscillations or mixed-mode oscillations, whereas large ones can make the disease go extinct in one patch and persist in the form of a positive steady state or a periodic solution in the other patch; (b) relative prevalence-based and unidirectional dispersal can make periodic outbreak earlier.

Evaluation of Effectiveness of Global COVID-19 Vaccination Campaign.

To model estimated deaths averted by COVID-19 vaccines, we used state-of-the-art mathematical modeling, likelihood-based inference, and reported COVID-19 death and vaccination data. We estimated that >1.5 million deaths were averted in 12 countries. Our model can help assess effectiveness of the vaccination program, which is crucial for curbing the COVID-19 pandemic.

Total biomass of a single population in two-patch environments.

For the two-patch logistic model, we study the effect of dispersal intensity and dispersal asymmetry on the total population abundance and its distribution. Two complete classifications of the model parameter space are given: one concerning when dispersal causes smaller or larger total biomass than no dispersal, and the other addressing how the total biomass changes with dispersal intensity and dispersal asymmetry. The dependencies of the population abundance of each individual patch on dispersal intensity and dispersal asymmetry are also fully characterized. In addition, the maximal and minimal total population sizes induced by dispersal are determined for the logistic model with an arbitrary number of patches, and a weak order-preserving result correlated the local population abundances with and without dispersal is established.

Shrinkage in serial intervals across transmission generations of COVID-19.

One of the key epidemiological characteristics that shape the transmission of coronavirus disease 2019 (COVID-19) is the serial interval (SI). Although SI is commonly considered following a probability distribution at a population scale, recent studies reported a slight shrinkage (or contraction) of the mean of effective SI across transmission generations or over time. Here, we develop a likelihood-based statistical inference framework with truncation to explore the change in SI across transmission generations after adjusting the impacts of case isolation. The COVID-19 contact tracing surveillance data in Hong Kong are used for exemplification. We find that for COVID-19, the mean of individual SI is likely to shrink with a factor at 0.72 per generation (95%CI: 0.54, 0.96) as the transmission generation increases, where a threshold may exist as the lower boundary of this shrinking process. We speculate that one of the probable explanations for the shrinkage in SI might be an outcome due to the competition among multiple candidate infectors within the same case cluster. Thus, the nonpharmaceutical interventive strategies are crucially important to block the transmission chains, and mitigate the COVID-19 epidemic.

A Zika Endemic Model for the Contribution of Multiple Transmission Routes.

Zika virus disease is a viral disease primarily transmitted to humans through the bite of infected female mosquitoes. Recent evidence indicates that the virus can also be sexually transmitted in hosts and vertically transmitted in vectors. In this paper, we propose a Zika model with three transmission routes, that is, vector-borne transmission between humans and mosquitoes, sexual transmission within humans and vertical transmission within mosquitoes. The basic reproduction number [Formula: see text] is computed and shown to be a sharp threshold quantity. Namely, the disease-free equilibrium is globally asymptotically stable as [Formula: see text], whereas there exists a unique endemic equilibrium which is globally asymptotically stable as [Formula: see text]. The relative contributions of each transmission route on the reproduction number, and the short- and long-term host infections are analyzed. Numerical simulations confirm that vectorial transmission contributes the most to the initial and subsequent transmission. The role of sexual transmission in the early phase of a Zika outbreak is greater than the long term, while vertical transmission is the opposite. Reducing mosquito bites is the most effective measure in lowering the risk of Zika virus infection.

Mechanistic modelling of the large-scale Lassa fever epidemics in Nigeria from 2016 to 2019.

Lassa fever, also known as Lassa hemorrhagic fever, is a virus that has generated recurrent outbreaks in West Africa. We use mechanistic modelling to study the Lassa fever epidemics in Nigeria from 2016-19. Our model describes the interaction between human and rodent populations with the consideration of quarantine, isolation and hospitalization processes. Our model supports the phenomenon of forward bifurcation where the stability between disease-free equilibrium and endemic equilibrium exchanges. Moreover, our model captures well the incidence curves from surveillance data. In particular, our model is able to reconstruct the periodic rodent and human forces of infection. Furthermore, we suggest that the three major epidemics from 2016-19 can be modelled by properly characterizing the rodent (or human) force of infection while the estimated human force of infection also present similar patterns across outbreaks. Our results suggest that the initial susceptibility likely increased across the three outbreaks from 2016-19. Our results highlight the similarity of the transmission dynamics driving three major Lassa fever outbreaks in the endemic areas.

Mathematical Analysis of the Ross-Macdonald Model with Quarantine.

People infected with malaria may receive less mosquito bites when they are treated in well-equipped hospitals or follow doctors' advice for reducing exposure to mosquitoes at home. This quarantine-like intervention measure is especially feasible in countries and areas approaching malaria elimination. Motivated by mathematical models with quarantine for directly transmitted diseases, we develop a mosquito-borne disease model where imperfect quarantine is considered to mitigate the disease transmission from infected humans to susceptible mosquitoes. The basic reproduction number [Formula: see text] is computed and the model equilibria and their stabilities are analyzed when the incidence rate is standard or bilinear. In particular, the model system may undergo a subcritical (backward) bifurcation at [Formula: see text] when standard incidence is adopted, whereas the disease-free equilibrium is globally asymptotically stable as [Formula: see text] and the unique endemic equilibrium is locally asymptotically stable as [Formula: see text] when the infection incidence is bilinear. Numerical simulations suggest that the quarantine strategy can play an important role in decreasing malaria transmission. The success of quarantine mainly relies on the reduction of bites on quarantined individuals.

Habitat fragmentation promotes malaria persistence.

Based on a Ross-Macdonald type model with a number of identical patches, we study the role of the movement of humans and/or mosquitoes on the persistence of malaria and many other vector-borne diseases. By using a theorem on line-sum symmetric matrices, we establish an eigenvalue inequality on the product of a class of nonnegative matrices and then apply it to prove that the basic reproduction number of the multipatch model is always greater than or equal to that of the single patch model. Biologically, this means that habitat fragmentation or patchiness promotes disease outbreaks and intensifies disease persistence. The risk of infection is minimized when the distribution of mosquitoes is proportional to that of humans. Numerical examples for the two-patch submodel are given to investigate how the multipatch reproduction number varies with human and/or mosquito movement. The reproduction number can surpass any given value whenever an appropriate travel pattern is chosen. Fast human and/or mosquito movement decreases the infection risk, but may increase the total number of infected humans.

Modelling diapause in mosquito population growth.

Diapause, a period of arrested development caused by adverse environmental conditions, serves as a key survival mechanism for insects and other invertebrate organisms in temperate and subtropical areas. In this paper, a novel modelling framework, motivated by mosquito species, is proposed to investigate the effects of diapause on seasonal population growth, where the diapause period is taken as an independent growth process, during which the population dynamics are completely different from that in the normal developmental and post-diapause periods. More specifically, the annual growth period is divided into three intervals, and the population dynamics during each interval are described by different sets of equations. We formulate two models of delay differential equations (DDE) to explicitly describe mosquito population growth with a single diapausing stage, either immature or adult. These two models can be further unified into one DDE model, on which the well-posedness of the solutions and the global stability of the trivial and positive periodic solutions in terms of an index [Formula: see text] are analysed. The seasonal population abundances of two temperate mosquito species with different diapausing stages are simulated to identify the essential role on population persistence that diapause plays and predict that killing mosquitoes during the diapause period can lower but fail to prevent the occurrence of peak abundance in the following season. Instead, culling mosquitoes during the normal growth period is much more efficient to decrease the outbreak size. Our modelling framework may shed light on the diapause-induced variations in spatiotemporal distributions of different mosquito species.

Modeling the 2016-2017 Yemen cholera outbreak with the impact of limited medical resources.

We present a mathematical model to investigate the transmission dynamics of the 2016-2017 Yemen cholera outbreak. Our model describes the interaction between the human hosts and the pathogenic bacteria, under the impact of limited medical resources. We fit our model to Yemen epidemic data published by the World Health Organization, at both the country and regional levels. We find that the Yemen cholera outbreak is shaped by the interplay of environmental, socioeconomic, and climatic factors. Our results suggest that improvement of the public health system and strategic implementation of control measures with respect to time and location are key to future cholera prevention and intervention in Yemen.

Impact of Awareness Programs on Cholera Dynamics: Two Modeling Approaches.

We propose two differential equation-based models to investigate the impact of awareness programs on cholera dynamics. The first model represents the disease transmission rates as decreasing functions of the number of awareness programs, whereas the second model divides the susceptible individuals into two distinct classes depending on their awareness/unawareness of the risk of infection. We study the essential dynamical properties of each model, using both analytical and numerical approaches. We find that the two models, though closely related, exhibit significantly different dynamical behaviors. Namely, the first model follows regular threshold dynamics while rich dynamical behaviors such as backward bifurcation may arise from the second one. Our results highlight the importance of validating key modeling assumptions in the development and selection of mathematical models toward practical application.

A PERIODIC ROSS-MACDONALD MODEL IN A PATCHY ENVIRONMENT.

Based on the classical Ross-Macdonald model, in this paper we propose a periodic malaria model to incorporate the effects of temporal and spatial heterogeneity on disease transmission. The temporal heterogeneity is described by assuming that some model coefficients are time-periodic, while the spatial heterogeneity is modeled by using a multi-patch structure and assuming that individuals travel among patches. We calculate the basic reproduction number [Formula: see text] and show that either the disease-free periodic solution is globally asymptotically stable if [Formula: see text] or the positive periodic solution is globally asymptotically stable if [Formula: see text]. Numerical simulations are conducted to confirm the analytical results and explore the effect of travel control on the disease prevalence.

Modeling the spatial spread of Rift Valley fever in Egypt.

Rift Valley fever (RVF) is a severe viral zoonosis in Africa and the Middle East that harms both human health and livestock production. It is believed that RVF in Egypt has been repeatedly introduced by the importation of infected animals from Sudan. In this paper, we propose a three-patch model for the process by which animals enter Egypt from Sudan, are moved up the Nile, and then consumed at population centers. The basic reproduction number for each patch is introduced and then the threshold dynamics of the model are established. We simulate an interesting scenario showing a possible explanation of the observed phenomenon of the geographic spread of RVF in Egypt.

Effects of behaviour change on HFMD transmission.

We propose a hand, foot and mouth disease (HFMD) transmission model for children with behaviour change and imperfect quarantine. The symptomatic and quarantined states obey constant behaviour change while others follow variable behaviour change depending on the numbers of new and recent infections. The basic reproduction number R0 of the model is defined and shown to be a threshold for disease persistence and eradication. Namely, the disease-free equilibrium is globally asymptotically stable if R0≤1 whereas the disease persists and there is a unique endemic equilibrium otherwise. By fitting the model to weekly HFMD data of Shanghai in 2019, the reproduction number is estimated at 2.41. Sensitivity analysis for R0 shows that avoiding contagious contacts and implementing strict quarantine are essential to lower HFMD persistence. Numerical simulations suggest that strong behaviour change not only reduces the peak size and endemic level dramatically but also impairs the role of asymptomatic transmission.

Towards a theory of ecotone resilience: coastal vegetation on a salinity gradient.

Ecotones represent locations where vegetation change is likely to occur as a result of climate and other environmental changes. Using a model of an ecotone vulnerable to such future changes, we estimated the resilience of the ecotone to disturbances. The specific ecotone is that between two different vegetation types, salinity-tolerant and salinity-intolerant, along a gradient in groundwater salinity. In the case studied, each vegetation type, through soil feedback loops, promoted local soil salinity levels that favor itself in competition with the other type. Bifurcation analysis was used to study the system of equations for the two vegetation types and soil salinity. Alternative stable equilibria, one for salinity-tolerant and one for salinity intolerant vegetation, were shown to exist over a region of the groundwater salinity gradient, bounded by two bifurcation points. This region was shown to depend sensitively on parameters such as the rate of upward infiltration of salinity from groundwater into the soil due to evaporation. We showed also that increasing diffusion rates of vegetation can lead to shrinkage of the range between the two bifurcation points. Sharp ecotones are typical of salt-tolerant vegetation (mangroves) near the coastline and salt-intolerant vegetation inland, even though the underlying elevation and groundwater salinity change very gradually. A disturbance such as an input of salinity to the soil from a storm surge could upset this stable boundary, leading to a regime shift of salinity-tolerant vegetation inland. We showed, however, that, for our model as least, a simple pulse disturbance would not be sufficient; the salinity would have to be held at a high level, as a 'press', for some time. The approach used here should be generalizable to study the resilience of a variety of ecotones to disturbances.

Modelling COVID-19 outbreak on the Diamond Princess ship using the public surveillance data.

The novel coronavirus disease 2019 (COVID-19) outbreak on the Diamond Princess (DP) ship has caused over 634 cases as of February 20, 2020. We model the transmission process on DP ship as a stochastic branching process, and estimate the reproduction number at the innitial phase of 2.9 (95%CrI: 1.7-7.7). The epidemic doubling time is 3.4 days, and thus timely actions on COVID-19 control were crucial. We estimate the COVID-19 transmissibility reduced 34% after the quarantine program on the DP ship which was implemented on February 5. According to the model simulation, relocating the population at risk may sustainably decrease the epidemic size, postpone the timing of epidemic peak, and thus relieve the tensive demands in the healthcare. The lesson learnt on the ship should be considered in other similar settings.

Impact of State-Dependent Dispersal on Disease Prevalence.

Based on a susceptible-infected-susceptible patch model, we study the influence of dispersal on the disease prevalence of an individual patch and all patches at the endemic equilibrium. Specifically, we estimate the disease prevalence of each patch and obtain a weak order-preserving result that correlated the patch reproduction number with the patch disease prevalence. Then we assume that dispersal rates of the susceptible and infected populations are proportional and derive the overall disease prevalence, or equivalently, the total infection size at no dispersal or infinite dispersal as well as the right derivative of the total infection size at no dispersal. Furthermore, for the two-patch submodel, two complete classifications of the model parameter space are given: one addressing when dispersal leads to higher or lower overall disease prevalence than no dispersal, and the other concerning how the overall disease prevalence varies with dispersal rate. Numerical simulations are performed to further investigate the effect of movement on disease prevalence.

Modelling triatomine bug population and Trypanosoma rangeli transmission dynamics: Co-feeding, pathogenic effect and linkage with chagas disease.

Trypanosoma rangeli (T. rangeli), a parasite, is not pathogenic to human but pathogenic to some vector species to induce the behavior changes of infected vectors and subsequently impact the transmission dynamics of other diseases such as Chagas disease which shares the same vector species. Here we develop a mathematical model and conduct qualitative analysis for the transmission dynamics of T. rangeli. We incorporate both systemic and co-feeding transmission routes, and account for the pathogenic effect using infection-induced fecundity and fertility change of the triatomine bugs. We derive two thresholds Rv (the triatomine bug basic reproduction number) and R0 (the T. rangeli basic reproduction number) to delineate the dynamical behaviors of the ecological and epidemiological systems. We show that when Rv>1 and R0>1, a unique parasite positive equilibrium E* appears. We find that E* can be unstable and periodic oscillations can be observed where the pathogenic effect plays a significant role. Implications of the qualitative analysis and numerical simulations suggest the need of an integrative vector-borne disease prevention and control strategy when multiple vector-borne diseases are transmitted by the same set of vector species.