Spatial Structure: Patch Models

Mathematical modeling at the livestock-wildlife interface: scoping review of drivers of disease transmission between species

Modeling of infectious diseases at the livestock-wildlife interface is a unique subset of mathematical modeling with many innate challenges. To ascertain the characteristics of the models used in these scenarios, a scoping review of the scientific literature was conducted. Fifty-six studies qualified for inclusion. Only 14 diseases at this interface have benefited from the utility of mathematical modeling, despite a far greater number of shared diseases. The most represented species combinations were cattle and badgers (for bovine tuberculosis, 14), and pigs and wild boar [for African (8) and classical (3) swine fever, and foot-and-mouth and disease (1)]. Assessing control strategies was the overwhelming primary research objective (27), with most studies examining control strategies applied to wildlife hosts and the effect on domestic hosts (10) or both wild and domestic hosts (5). In spatially-explicit models, while livestock species can often be represented through explicit and identifiable location data (such as farm, herd, or pasture locations), wildlife locations are often inferred using habitat suitability as a proxy. Though there are innate assumptions that may not be fully accurate when using habitat suitability to represent wildlife presence, especially for wildlife the parsimony principle plays a large role in modeling diseases at this interface, where parameters are difficult to document or require a high level of data for inference. Explaining observed transmission dynamics was another common model objective, though the relative contribution of involved species to epizootic propagation was only ascertained in a few models. More direct evidence of disease spill-over, as can be obtained through genomic approaches based on pathogen sequences, could be a useful complement to further inform such modeling. As computational and programmatic capabilities advance, the resolution of the models and data used in these models will likely be able to increase as well, with a potential goal being the linking of modern complex ecological models with the depth of dynamics responsible for pathogen transmission. Controlling diseases at this interface is a critical step toward improving both livestock and wildlife health, and mechanistic models are becoming increasingly used to explore the strategies needed to confront these diseases.

Spatio-temporal stochastic differential equations for crime incidence modeling

We propose a methodology for the quantitative fitting and forecasting of real spatio-temporal crime data, based on stochastic differential equations. The analysis is focused on the city of Valencia, Spain, for which 90247 robberies and thefts with their latitude-longitude positions are available for a span of eleven years (2010-2020) from records of the 112-emergency phone. The incidents are placed in the 26 zip codes of the city (46001-46026), and monthly time series of crime are built for each of the zip codes. Their annual-trend components are modeled by Itô diffusion, with jointly correlated noises to account for district-level relations. In practice, this study may help simulate spatio-temporal situations and identify risky areas and periods from present and past data.

A linear noise approximation for stochastic epidemic models fit to partially observed incidence counts

Stochastic epidemic models (SEMs) fit to incidence data are critical to elucidating outbreak dynamics, shaping response strategies, and preparing for future epidemics. SEMs typically represent counts of individuals in discrete infection states using Markov jump processes (MJPs), but are computationally challenging as imperfect surveillance, lack of subject-level information, and temporal coarseness of the data obscure the true epidemic. Analytic integration over the latent epidemic process is impossible, and integration via Markov chain Monte Carlo (MCMC) is cumbersome due to the dimensionality and discreteness of the latent state space. Simulation-based computational approaches can address the intractability of the MJP likelihood, but are numerically fragile and prohibitively expensive for complex models. A linear noise approximation (LNA) that approximates the MJP transition density with a Gaussian density has been explored for analyzing prevalence data in large-population settings, but requires modification for analyzing incidence counts without assuming that the data are normally distributed. We demonstrate how to reparameterize SEMs to appropriately analyze incidence data, and fold the LNA into a data augmentation MCMC framework that outperforms deterministic methods, statistically, and simulation-based methods, computationally. Our framework is computationally robust when the model dynamics are complex and applies to a broad class of SEMs. We evaluate our method in simulations that reflect Ebola, influenza, and SARS-CoV-2 dynamics, and apply our method to national surveillance counts from the 2013-2015 West Africa Ebola outbreak.

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.

Methods for Rapid Mobility Estimation to Support Outbreak Response

When pressed for time, outbreak investigators often use homogeneous mixing models to model infectious diseases in data-poor regions. But recent outbreaks such as the 2014 Ebola outbreak in West Africa have shown the limitations of this approach in an era of increasing urbanization and connectivity. Both outbreak detection and predictive modeling depend on realistic estimates of human and disease mobility, but these data are difficult to acquire in a timely manner. This is especially true when dealing with an emerging outbreak in an under-resourced nation. Weighted travel networks with realistic estimates for population flows are often proprietary, expensive, or nonexistent. Here we propose a method for rapidly generating a mobility model from open-source data. As an example, we use road and river network data, along with population estimates, to construct a realistic model of human movement between health zones in the Democratic Republic of the Congo (DRC). Using these mobility data, we then fit an epidemic model to real-world surveillance data from the recent Ebola outbreak in the Nord Kivu region of the DRC to illustrate a potential use of the generated mobility estimation. In addition to providing a way for rapid risk estimation, this approach brings together novel techniques to merge diverse GIS datasets that can then be used to address issues that pertain to public health and global health security.

The SIR dynamic model of infectious disease transmission and its analogy with chemical kinetics

Mathematical models of the dynamics of infectious disease transmission are used to forecast epidemics and assess mitigation strategies. In this article, we highlight the analogy between the dynamics of disease transmission and chemical reaction kinetics while providing an exposition on the classic Susceptible–Infectious–Removed (SIR) epidemic model. Particularly, the SIR model resembles a dynamic model of a batch reactor carrying out an autocatalytic reaction with catalyst deactivation. This analogy between disease transmission and chemical reaction enables the exchange of ideas between epidemic and chemical kinetic modeling communities.

ANALYSIS OF A REACTION–DIFFUSION EPIDEMIC MODEL WITH ASYMPTOMATIC TRANSMISSION

We develop a dynamic model of an evolving epidemic in a spatially inhomogeneous environment. We analyze the model, as a system of reaction–diffusion partial differential equations, to predict the outbreak and spatio-temporal spread of the disease. The model features both an asymptomatic infectious stage and symptomatic infectious stage, with both the asymptomatic and the symptomatic infected populations dispersing through the susceptible population. We prove the existence and uniform boundedness of solutions, and investigate their long-time behavior. We apply the spatially homogeneous version of the model to the current COVID-19 epidemic in Brazil.

The transmission dynamic of different hepatitis B-infected individuals with the effect of hospitalization

We propose an epidemic model for the transmission of hepatitis B virus along with the classification of different infection phases and hospitalized class. We formulate the model and discuss its basic mathematical properties, e.g. existence, positivity, and biological feasibility. Exploiting the next generation matrix approach, we find the basic reproductive number of the model. We perform sensitivity analysis to illustrate the effect of various parameters on the transmission of the disease. We investigate stability of the equilibria of the model in terms of the basic reproduction number. Conditions for the stability of the proposed model are obtained using various approaches. Finally, we perform the numerical simulations to discuss sensitivity analysis and to support our analytical work.

Decision Support for Mitigation of Livestock Disease: Rinderpest as a Case Study

A versatile, interactive model to predict geographically resolved epidemic progression after pathogen introduction into a population is presented. Deterministic simulations incorporating a compartmental disease model run rapidly, facilitating the analysis of mitigations such as vaccination and transmission reduction on epidemic spread and progression. We demonstrate the simulation model using rinderpest infection of cattle, a devastating livestock disease. Rinderpest has been extinguished in the wild, but it is still a threat due to stored virus in some laboratories. Comparison of simulations to historical outbreaks provides some validation of the model. Simulations of potential outbreaks demonstrate potential consequences of rinderpest virus release for a variety of possible disease parameters and mitigations. Our results indicate that a rinderpest outbreak could result in severe social and economic consequences.

Numerical solution of a spatio-temporal gender-structured model for hantavirus infection in rodents

In this article we describe the transmission dynamics of hantavirus in rodents using a spatio-temporal susceptible-exposed-infective-recovered (SEIR) compartmental model that distinguishes between male and female subpopulations [L.J.S. Allen, R.K. McCormack and C.B. Jonsson, Bull. Math. Biol. 68 (2006), 511--524]. Both subpopulations are assumed to differ in their movement with respect to local variations in the densities of their own and the opposite gender group. Three alternative models for the movement of the male individuals are examined. In some cases the movement is not only directed by the gradient of a density (as in the standard diffusive case), but also by a non-local convolution of density values as proposed, in another context, in [R.M. Colombo and E. Rossi, Commun. Math. Sci., 13 (2015), 369--400]. An efficient numerical method for the resulting convection-diffusion-reaction system of partial differential equations is proposed. This method involves techniques of weighted essentially non-oscillatory (WENO) reconstructions in combination with implicit-explicit Runge-Kutta (IMEX-RK) methods for time stepping. The numerical results demonstrate significant differences in the spatio-temporal behavior predicted by the different models, which suggest future research directions.