Multi-patch and multi-group epidemic models: a new framework

Delay epidemic models determined by latency, infection, and immunity duration

We propose new single and two-strain epidemic models represented by systems of delay differential equations and based on the number of newly exposed individuals. Transitions between exposed, infectious, recovered, and back to susceptible compartments are determined by the corresponding time delays. Existence and positiveness of solutions are proved. Reduction of delay differential equations to integral equations allows the analysis of stationary solutions and their stability. In the case of two strains, they compete with each other, and the strain with a larger individual basic reproduction number dominates the other one. However, if the basic reproduction number exceeds some critical values, stationary solution loses its stability resulting in periodic time oscillations. In this case, both strains are present and their dynamics is not completely determined by the basic reproduction numbers but also by other parameters. The results of the work are illustrated by comparison with data on seasonal influenza.

The impacts of spatial-temporal heterogeneity of human-to-human contacts on the extinction probability of infectious disease from branching process model

In this study, we focus on the impacts of spatial-temporal heterogeneity of human-to-human contacts on the spread of infectious diseases and develop a multi-type branching process model by introducing random human-to-human contact mode into a structured population. We provide the general formulas of the generation size, extinction probability, and basic reproduction number of the proposed branching process model. The result shows that the natural temporal heterogeneity (i.e. random contacts over time) can lead to a higher extinction probability while remains the same basic reproduction number and generation size. This is also numerically verified by choosing the real contact distributions from different circumstances of four countries. In addition, we observe a non-monotonic pattern of the differences, against the transmission probability and the mean contact rate, between the extinction probabilities under the constant and random contact patterns. Given the spatial heterogeneity, we show that it can contribute to the increase of basic reproduction number, but also increase the extinction probability of the infectious disease. This study adds novel insights to the course of the impact of heterogeneity on the transmission dynamics and also provides additional evidence for the limited role of reproduction numbers.

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.

Characterization of differential susceptibility and differential infectivity epidemic models

Heterogeneity in susceptibility and infectivity is a central issue in epidemiology. Although the latter has received some attention recently, the former is often neglected in modeling of epidemic systems. Moreover, very few studies consider both of these heterogeneities. This paper is concerned with the characterization of epidemic models with differential susceptibility and differential infectivity under a general setup. Specifically, we investigate the global asymptotic behavior of equilibria of these systems when the network configuration of the Susceptible-Infectious interactions is strongly connected. These results prove two conjectures by Bonzi et al. (J Math Biol 62:39-64, 2011) and Hyman and Li (Math Biosci Eng 3:89-100, 2006). Moreover, we consider the scenario in which the strong connectivity hypothesis is dropped. In this case, the model exhibits a wider range of dynamical behavior, including the rise of boundary and interior equilibria, all based on the topology of network connectivity. Finally, a model with multidirectional transitions between infectious classes is presented and completely analyzed.

On a two-strain epidemic model involving delay equations

We propose an epidemiological model for the interaction of either two viruses or viral strains with cross-immunity, where the individuals infected by the first virus cannot be infected by the second one, and without cross-immunity, where a secondary infection can occur. The model incorporates distributed recovery and death rates and consists of integro-differential equations governing the dynamics of susceptible, infectious, recovered, and dead compartments. Assuming that the recovery and death rates are uniformly distributed in time throughout the duration of the diseases, we can simplify the model to a conventional ordinary differential equation (ODE) model. Another limiting case arises if the recovery and death rates are approximated by the delta-function, thereby resulting in a new point-wise delay model that incorporates two time delays corresponding to the durations of the diseases. We establish the positiveness of solutions for the distributed delay models and determine the basic reproduction number and an estimate for the final size of the epidemic for the delay model. According to the results of the numerical simulations, both strains can coexist in the population if the disease transmission rates for them are close to each other. If the difference between them is sufficiently large, then one of the strains dominates and eliminates the other one.

Inference on a Multi-Patch Epidemic Model with Partial Mobility, Residency, and Demography: Case of the 2020 COVID-19 Outbreak in Hermosillo, Mexico

Most studies modeling population mobility and the spread of infectious diseases, particularly those using meta-population multi-patch models, tend to focus on the theoretical properties and numerical simulation of such models. As such, there is relatively scant literature focused on numerical fit, inference, and uncertainty quantification of epidemic models with population mobility. In this research, we use three estimation techniques to solve an inverse problem and quantify its uncertainty for a human-mobility-based multi-patch epidemic model using mobile phone sensing data and confirmed COVID-19-positive cases in Hermosillo, Mexico. First, we utilize a Brownian bridge model using mobile phone GPS data to estimate the residence and mobility parameters of the epidemic model. In the second step, we estimate the optimal model epidemiological parameters by deterministically inverting the model using a Darwinian-inspired evolutionary algorithm (EA)-that is, a genetic algorithm (GA). The third part of the analysis involves performing inference and uncertainty quantification in the epidemic model using two Bayesian Monte Carlo sampling methods: t-walk and Hamiltonian Monte Carlo (HMC). The results demonstrate that the estimated model parameters and incidence adequately fit the observed daily COVID-19 incidence in Hermosillo. Moreover, the estimated parameters from the HMC method yield large credible intervals, improving their coverage for the observed and predicted daily incidences. Furthermore, we observe that the use of a multi-patch model with mobility yields improved predictions when compared to a single-patch model.

An epidemic model with time delays determined by the infectivity and disease durations

We propose an epidemiological model with distributed recovery and death rates. It represents an integrodifferential system of equations for susceptible, exposed, infectious, recovered and dead compartments. This model can be reduced to the conventional ODE model under the assumption that recovery and death rates are uniformly distributed in time during disease duration. Another limiting case, where recovery and death rates are given by the delta-function, leads to a new point-wise delay model with two time delays corresponding to the infectivity period and disease duration. Existence and positiveness of solutions for the distributed delay model and point-wise delay model are proved. The basic reproduction number and the final size of the epidemic are determined. Both, the ODE model and the delay models are used to describe COVID-19 epidemic progression. The delay model gives a better approximation of the Omicron data than the conventional ODE model from the point of view of parameter estimation.

An age-dependent immuno-epidemiological model with distributed recovery and death rates

The work is devoted to a new immuno-epidemiological model with distributed recovery and death rates considered as functions of time after the infection onset. Disease transmission rate depends on the intra-subject viral load determined from the immunological submodel. The age-dependent model includes the viral load, recovery and death rates as functions of age considered as a continuous variable. Equations for susceptible, infected, recovered and dead compartments are expressed in terms of the number of newly infected cases. The analysis of the model includes the proof of the existence and uniqueness of solution. Furthermore, it is shown how the model can be reduced to age-dependent SIR or delay model under certain assumptions on recovery and death distributions. Basic reproduction number and final size of epidemic are determined for the reduced models. The model is validated with a COVID-19 case data. Modelling results show that proportion of young age groups can influence the epidemic progression since disease transmission rate for them is higher than for other age groups.

The transmission mechanism theory of disease dynamics: Its aims, assumptions and limitations

Most of the progress in the development of single scale mathematical and computational models for the study of infectious disease dynamics which now span over a century is build on a body of knowledge that has been developed to address particular single scale descriptions of infectious disease dynamics based on understanding disease transmission process. Although this single scale understanding of infectious disease dynamics is now founded on a body of knowledge with a long history, dating back to over a century now, that knowledge has not yet been formalized into a scientific theory. In this article, we formalize this accumulated body of knowledge into a scientific theory called the transmission mechanism theory of disease dynamics which states that at every scale of organization of an infectious disease system, disease dynamics is determined by transmission as the main dynamic disease process. Therefore, the transmission mechanism theory of disease dynamics can be seen as formalizing knowledge that has been inherent in the study of infectious disease dynamics using single scale mathematical and computational models for over a century now. The objective of this article is to summarize this existing knowledge about single scale modelling of infectious dynamics by means of a scientific theory called the transmission mechanism theory of disease dynamics and highlight its aims, assumptions and limitations.

Modeling malaria transmission in Nepal: impact of imported cases through cross-border mobility

The cross-border mobility of malaria cases poses an obstacle to malaria elimination programmes in many countries, including Nepal. Here, we develop a novel mathematical model to study how the imported malaria cases through the Nepal-India open-border affect the Nepal government's goal of eliminating malaria by 2026. Mathematical analyses and numerical simulations of our model, validated by malaria case data from Nepal, indicate that eliminating malaria from Nepal is possible if strategies promoting the absence of cross-border mobility, complete protection of transmission abroad, or strict border screening and isolation are implemented. For each strategy, we establish the conditions for the elimination of malaria. We further use our model to identify the control strategies that can help maintain a low endemic level. Our results show that the ideal control strategies should be designed according to the average mosquito biting rates that may depend on the location and season.

A patchy theoretical model for the transmission dynamics of SARS-Cov-2 with optimal control

Short-term human movements play a major part in the transmission and control of COVID-19, within and between countries. Such movements are necessary to be included in mathematical models that aim to assist in understanding the transmission dynamics of COVID-19. A two-patch basic mathematical model for COVID-19 was developed and analyzed, incorporating short-term human mobility. Here, we modeled the human mobility that depended on its epidemiological status, by the Lagrangian approach. A sharp threshold for disease dynamics known as the reproduction number was computed. Particularly, we portrayed that when the disease threshold is less than unity, the disease dies out and the disease persists when the reproduction number is greater than unity. Optimal control theory was also applied to the proposed model, with the aim of investigating the cost-effectiveness strategy. The findings were further investigated through the usage of the results from the cost objective functional, the average cost-effectiveness ratio (ACER), and then the infection averted ratio (IAR).

Analysis of War and Conflict Effect on the Transmission Dynamics of the Tenth Ebola Outbreak in the Democratic Republic of Congo

The tenth Ebola outbreak in the Democratic Republic of Congo (DRC) that occurred from 2018 to 2020 was exacerbated by long-lasting conflicts and war in the region. We propose a deterministic model to investigate the impact of such disruptive events on the transmission dynamics of the Ebola virus disease. It is an extension of the classical susceptible-infectious-recovered model, enriched by an additional class of contaminated environment to account for indirect transmission as well as two classes of hospitalized individuals and patients who escape from the healthcare facility due to violence and attacks perpetrated by armed groups, rebels, etc. The model is formulated using two patches, namely Patch 1 consisting of the three affected eastern provinces in DRC and Patch 2, a war- and conflict-free area consisting of the go-to neighboring provinces for escaped patients. We introduce two key parameters, the escaping rate from hospitals and the destruction of hospitals, in terms of which the effect of war and conflicts is measured. The model is fitted and parameterized using the cumulative mortality data from the region. The basic reproduction number [Formula: see text] is computed and found to have a complex expression due to the high nonlinearity of the model. By using, not a Lyapunov function, but a decomposition theorem in Castillo-Chavez et al.(in Castillo-Chavez et al. (eds) Mathematical approaches for emerging and reemerging infectious diseases: an introduction, vol 126. Springer Science & Business Media, Berlin, 2002), it is shown that the disease-free equilibrium is globally asymptotically stable when [Formula: see text] and unstable when [Formula: see text]. A nonstandard finite difference scheme which replicates the dynamics of the continuous model is designed. In particular, a discrete counterpart of the above-mentioned theorem on the global asymptotic stability of the disease-free equilibrium is investigated. Numerical experiments are presented to support the theoretical results. When [Formula: see text], the numerical simulations suggest that there exists for the full model a unique globally asymptotically stable interior endemic equilibrium point, while it is shown theoretically and computationally that the model possesses at least a one Patch 1 and a one Patch 2 boundary equilibria (i.e., Patch 2 and Patch 1 disease-free equilibrium) points, which are locally asymptotically stable. Some recommendations to tackle Ebola in a conflict zone are stated.

Estimating the effective reproduction number for heterogeneous models using incidence data

The effective reproduction number, R ( t ) , plays a key role in the study of infectious diseases, indicating the current average number of new infections caused by an infected individual in an epidemic process. Estimation methods for the time evolution of R ( t ) , using incidence data, rely on the generation interval distribution, g(τ), which is usually obtained from empirical data or theoretical studies using simple epidemic models. However, for systems that present heterogeneity, either on the host population or in the expression of the disease, there is a lack of data and of a suitable general methodology to obtain g(τ). In this work, we use mathematical models to bridge this gap. We present a general methodology for obtaining explicit expressions of the reproduction numbers and the generation interval distributions, within and between model sub-compartments provided by an arbitrary compartmental model. Additionally, we present the appropriate expressions to evaluate those reproduction numbers using incidence data. To highlight the relevance of such methodology, we apply it to the spread of COVID-19 in municipalities of the state of Rio de Janeiro, Brazil. Using two meta-population models, we estimate the reproduction numbers and the contributions of each municipality in the generation of cases in all others.

Estimating the effective reproduction number for heterogeneous models using incidence data

The effective reproduction number, R ( t ) , plays a key role in the study of infectious diseases, indicating the current average number of new infections caused by an infected individual in an epidemic process. Estimation methods for the time evolution of R ( t ) , using incidence data, rely on the generation interval distribution, g(τ), which is usually obtained from empirical data or theoretical studies using simple epidemic models. However, for systems that present heterogeneity, either on the host population or in the expression of the disease, there is a lack of data and of a suitable general methodology to obtain g(τ). In this work, we use mathematical models to bridge this gap. We present a general methodology for obtaining explicit expressions of the reproduction numbers and the generation interval distributions, within and between model sub-compartments provided by an arbitrary compartmental model. Additionally, we present the appropriate expressions to evaluate those reproduction numbers using incidence data. To highlight the relevance of such methodology, we apply it to the spread of COVID-19 in municipalities of the state of Rio de Janeiro, Brazil. Using two meta-population models, we estimate the reproduction numbers and the contributions of each municipality in the generation of cases in all others.

Estimating the effective reproduction number for heterogeneous models using incidence data

The effective reproduction number, R ( t ) , plays a key role in the study of infectious diseases, indicating the current average number of new infections caused by an infected individual in an epidemic process. Estimation methods for the time evolution of R ( t ) , using incidence data, rely on the generation interval distribution, g(τ), which is usually obtained from empirical data or theoretical studies using simple epidemic models. However, for systems that present heterogeneity, either on the host population or in the expression of the disease, there is a lack of data and of a suitable general methodology to obtain g(τ). In this work, we use mathematical models to bridge this gap. We present a general methodology for obtaining explicit expressions of the reproduction numbers and the generation interval distributions, within and between model sub-compartments provided by an arbitrary compartmental model. Additionally, we present the appropriate expressions to evaluate those reproduction numbers using incidence data. To highlight the relevance of such methodology, we apply it to the spread of COVID-19 in municipalities of the state of Rio de Janeiro, Brazil. Using two meta-population models, we estimate the reproduction numbers and the contributions of each municipality in the generation of cases in all others.

An Epidemic Model with Time Delay Determined by the Disease Duration

Immuno-epidemiological models with distributed recovery and death rates can describe the epidemic progression more precisely than conventional compartmental models. However, the required immunological data to estimate the distributed recovery and death rates are not easily available. An epidemic model with time delay is derived from the previously developed model with distributed recovery and death rates, which does not require precise immunological data. The resulting generic model describes epidemic progression using two parameters, disease transmission rate and disease duration. The disease duration is incorporated as a delay parameter. Various epidemic characteristics of the delay model, namely the basic reproduction number, the maximal number of infected, and the final size of the epidemic are derived. The estimation of disease duration is studied with the help of real data for COVID-19. The delay model gives a good approximation of the COVID-19 data and of the more detailed model with distributed parameters.

Managing public transit during a pandemic: The trade-off between safety and mobility

During a pandemic such as COVID-19, managing public transit effectively becomes a critical policy decision. On the one hand, efficient transportation plays a pivotal role in enabling the movement of essential workers and keeping the economy moving. On the other hand, public transit can be a vector for disease propagation due to travelers' proximity within shared and enclosed spaces. Without strategic preparedness, mass transit facilities are potential hotbeds for spreading infectious diseases. Thus, transportation agencies face a complex trade-off when developing context-specific operating strategies for public transit. This work provides a network-based analysis framework for understanding this trade-off, as well as tools for calculating targeted commute restrictions under different policy constraints, e.g., regarding public health considerations (limiting infection levels) and economic activity (limiting the reduction in travel). The resulting plans ensure that the traffic flow restrictions imposed on each route are adaptive to the time-varying epidemic dynamics. A case study based on the COVID-19 pandemic reveals that a well-planned subway system in New York City can sustain 88% of transit flow while reducing the risk of disease transmission by 50% relative to fully-loaded public transit systems. Transport policy-makers can exploit this optimization-based framework to address safety-and-mobility trade-offs and make proactive transit management plans during an epidemic outbreak.

Geographical network model for COVID-19 spread among dynamic epidemic regions

Pandemic due to SARS-CoV-2 (COVID-19) has affected to world in several aspects: high number of confirmed cases, high number of deaths, low economic growth, among others. Understanding of spatio-temporal dynamics of the virus is helpful and necessary for decision making, for instance to decide where, whether and how, non-pharmaceutical intervention policies are to be applied. This point has not been properly addressed in literature since typical strategies do not consider marked differences on the epidemic spread across country or large territory. Those strategies assume similarities and apply similar interventions instead. This work is focused on posing a methodology where spatio-temporal epidemic dynamics is captured by means of dividing a territory in time-varying epidemic regions, according to geographical closeness and infection level. In addition, a novel Lagrangian-SEIR-based model is posed for describing the dynamic within and between those regions. The capabilities of this methodology for identifying local outbreaks and reproducing the epidemic curve are discussed for the case of COVID-19 epidemic in Jalisco state (Mexico). The contagions from July 31, 2020 to March 31, 2021 are analyzed, with monthly adjustments, and the estimates obtained at the level of the epidemic regions present satisfactory results since Relative Root Mean Squared Error RRMSE is below 15% in most of regions, and at the level of the whole state outstanding with RRMSE below 5%.

A multi-source global-local model for epidemic management

The Effective Reproduction Number Rt provides essential information for the management of an epidemic/pandemic. Projecting Rt into the future could further assist in the management process. This article proposes a methodology based on exposure scenarios to perform such a procedure. The method utilizes a compartmental model and its adequate parametrization; a way to determine suitable parameters for this model in México's case is detailed. In conjunction with the compartmental model, the projection of Rt permits estimating unobserved variables, such as the size of the asymptomatic population, and projecting into the future other relevant variables, like the active hospitalizations, using scenarios. The uses of the proposed methodologies are exemplified by analyzing the pandemic in a Mexican state; the main quantities derived from the compartmental model, such as the active and total cases, are included in the analysis. This article also presents a national summary based on the methodologies to illustrate how these procedures could be further exploited. The supporting information includes an application of the proposed methods to a metropolitan area to show that it also works well at other demographic disaggregation levels. The procedures developed in this article shed light on how to develop an effective surveillance system when information is incomplete and can be applied in cases other than México's.

Dynamical analysis of an SIS epidemic model with migration and residence time

Migration can be divided into temporary and permanent migration, which is related to the residence time of people in the patch, thus we consider an SIS epidemic model with migration and residence time in a patchy environment. If [Formula: see text], the disease-free equilibrium is globally asymptotically stable and the disease dies out. With the same migration rate of susceptible and infectious individuals and without disease-induced death, when [Formula: see text], the endemic equilibrium is unique and globally asymptotically stable. Numerical simulations are carried out to show the effects of residence time and the migration rate on disease prevalence.

Dynamics and asymptotic profiles of endemic equilibrium for SIS epidemic patch models

In this paper, we perform qualitative analysis to two SIS epidemic models in a patchy environment, without and with linear recruitment. The model without linear recruitment was proposed and studied by Allen et al. (SIAM J Appl Math 67(5):1283-1309, 2007). This model possesses a conserved total population number, whereas the model with linear recruitment has a varying total population. However, both models have the same basic reproduction number. For both models, we establish the global stability of endemic equilibrium in a special case, which partially solves an open problem. Then we investigate the asymptotic behavior of endemic equilibrium as the mobility of infected and/or susceptible population tends to zero. Though the basic reproduction number is a well-known critical index, our theoretical results strongly suggest that other factors such as the variation of total population number and individual movement may also play vital roles in disease prediction and control. In particular, our results imply that the variation of total population number can cause infectious disease to become more threatening and difficult to control.