Stochastic epidemics: major outbreaks and the duration of the endemic period

Using early detection data to estimate the date of emergence of an epidemic outbreak

While the first infection of an emerging disease is often unknown, information on early cases can be used to date it. In the context of the COVID-19 pandemic, previous studies have estimated dates of emergence (e.g., first human SARS-CoV-2 infection, emergence of the Alpha SARS-CoV-2 variant) using mainly genomic data. Another dating attempt used a stochastic population dynamics approach and the date of the first reported case. Here, we extend this approach to use a larger set of early reported cases to estimate the delay from first infection to the Nth case. We first validate our framework by running our model on simulated data. We then apply our model using data on Alpha variant infections in the UK, dating the first Alpha infection at (median) August 21, 2020 (95% interpercentile range across retained simulations (IPR): July 23-September 5, 2020). Next, we apply our model to data on COVID-19 cases with symptom onset before mid-January 2020. We date the first SARS-CoV-2 infection in Wuhan at (median) November 28, 2019 (95% IPR: November 2-December 9, 2019). Our results fall within ranges previously estimated by studies relying on genomic data. Our population dynamics-based modelling framework is generic and flexible, and thus can be applied to estimate the starting time of outbreaks in contexts other than COVID-19.

The probability of epidemic burnout in the stochastic SIR model with vital dynamics

We present an approach to computing the probability of epidemic "burnout," i.e., the probability that a newly emergent pathogen will go extinct after a major epidemic. Our analysis is based on the standard stochastic formulation of the Susceptible-Infectious-Removed (SIR) epidemic model including host demography (births and deaths) and corresponds to the standard SIR ordinary differential equations (ODEs) in the infinite population limit. Exploiting a boundary layer approximation to the ODEs and a birth-death process approximation to the stochastic dynamics within the boundary layer, we derive convenient, fully analytical approximations for the burnout probability. We demonstrate-by comparing with computationally demanding individual-based stochastic simulations and with semi-analytical approximations derived previously-that our fully analytical approximations are highly accurate for biologically plausible parameters. We show that the probability of burnout always decreases with increased mean infectious period. However, for typical biological parameters, there is a relevant local minimum in the probability of persistence as a function of the basic reproduction number [Formula: see text]. For the shortest infectious periods, persistence is least likely if [Formula: see text]; for longer infectious periods, the minimum point decreases to [Formula: see text]. For typical acute immunizing infections in human populations of realistic size, our analysis of the SIR model shows that burnout is almost certain in a well-mixed population, implying that susceptible recruitment through births is insufficient on its own to explain disease persistence.

The effect of screening on the health burden of chlamydia: An evaluation of compartmental models based on person-days of infection

Sexually transmitted diseases (STDs) are detrimental to the health and economic well-being of society. Consequently, predicting outbreaks and identifying effective disease interventions through epidemiological tools, such as compartmental models, is of the utmost importance. Unfortunately, the ordinary differential equation compartmental models attributed to the work of Kermack and McKendrick require a duration of infection that follows the exponential or Erlang distribution, despite the biological invalidity of such assumptions. As these assumptions negatively impact the quality of predictions, alternative approaches are required that capture how the variability in the duration of infection affects the trajectory of disease and the evaluation of disease interventions. So, we apply a new family of ordinary differential equation compartmental models based on the quantity person-days of infection to predict the trajectory of disease. Importantly, this new family of models features non-exponential and non-Erlang duration of infection distributions without requiring more complex integral and integrodifferential equation compartmental model formulations. As proof of concept, we calibrate our model to recent trends of chlamydia incidence in the U.S. and utilize a novel duration of infection distribution that features periodic hazard rates. We then evaluate how increasing STD screening rates alter predictions of incidence and disability adjusted life-years over a five-year horizon. Our findings illustrate that our family of compartmental models provides a better fit to chlamydia incidence trends than traditional compartmental models, based on Akaike information criterion. They also show new asymptomatic and symptomatic infections of chlamydia peak over drastically different time frames and that increasing the annual STD screening rates from 35% to 40%-70% would annually avert 6.1-40.3 incidence while saving 1.68-11.14 disability adjusted life-years per 1000 people. This suggests increasing the STD screening rate in the U.S. would greatly aid in ongoing public health efforts to curtail the rising trends in preventable STDs.

Global asymptotic stability, extinction and ergodic stationary distribution in a stochastic model for dual variants of SARS-CoV-2

Several mathematical models have been developed to investigate the dynamics SARS-CoV-2 and its different variants. Most of the multi-strain SARS-CoV-2 models do not capture an important and more realistic feature of such models known as randomness. As the dynamical behavior of most epidemics, especially SARS-CoV-2, is unarguably influenced by several random factors, it is appropriate to consider a stochastic vaccination co-infection model for two strains of SARS-CoV-2. In this work, a new stochastic model for two variants of SARS-CoV-2 is presented. The conditions of existence and the uniqueness of a unique global solution of the stochastic model are derived. Constructing an appropriate Lyapunov function, the conditions for the stochastic system to fluctuate around endemic equilibrium of the deterministic system are derived. Stationary distribution and ergodicity for the new co-infection model are also studied. Numerical simulations are carried out to validate theoretical results. It is observed that when the white noise intensities are larger than certain thresholds and the associated stochastic reproduction numbers are less than unity, both strains die out and go into extinction with unit probability. More-over, it is observed that, for weak white noise intensities, the solution of the stochastic system fluctuates around the endemic equilibrium (EE) of the deterministic model. Frequency distributions are also studied to show random fluctuations due to stochastic white noise intensities. The results presented herein also reveal the impact of vaccination in reducing the co-circulation of SARS-CoV-2 variants within a given population.

Threshold Dynamics and the Density Function of the Stochastic Coronavirus Epidemic Model

Since November 2019, each country in the world has been affected by COVID-19, which has claimed more than four million lives. As an infectious disease, COVID-19 has a stronger transmission power and faster propagation speed. In fact, environmental noise is an inevitable important factor in the real world. This paper mainly gives a new random infectious disease system under infection rate environmental noise. We give the existence and uniqueness of the solution of the system and discuss the ergodic stationary distribution and the extinction conditions of the system. The probability density function of the stochastic system is studied. Some digital simulations are used to demonstrate the probability density function and the extinction of the system.

A Temperature Conditioned Markov Chain Model for Predicting the Dynamics of Mosquito Vectors of Disease

Understanding and predicting mosquito population dynamics is crucial for gaining insight into the abundance of arthropod disease vectors and for the design of effective vector control strategies. In this work, a climate-conditioned Markov chain (CMC) model was developed and applied for the first time to predict the dynamics of vectors of important medical diseases. Temporal changes in mosquito population profiles were generated to simulate the probabilities of a high population impact. The simulated transition probabilities of the mosquito populations achieved from the trained model are very near to the observed data transitions that have been used to parameterize and validate the model. Thus, the CMC model satisfactorily describes the temporal evolution of the mosquito population process. In general, our numerical results, when temperature is considered as the driver of change, indicate that it is more likely for the population system to move into a state of high population level when the former is a state of a lower population level than the opposite. Field data on frequencies of successive mosquito population levels, which were not used for the data inferred MC modeling, were assembled to obtain an empirical intensity transition matrix and the frequencies observed. Our findings match to a certain degree the empirical results in which the probabilities follow analogous patterns while no significant differences were observed between the transition matrices of the CMC model and the validation data (ChiSq = 14.58013, df = 24, p = 0.9324451). The proposed modeling approach is a valuable eco-epidemiological study. Moreover, compared to traditional Markov chains, the benefit of the current CMC model is that it takes into account the stochastic conditional properties of ecological-related climate variables. The current modeling approach could save costs and time in establishing vector eradication programs and mosquito surveillance programs.

Threshold Dynamics of a Non-Linear Stochastic Viral Model with Time Delay and CTL Responsiveness

This article focuses on a stochastic viral model with distributed delay and CTL responsiveness. It is shown that the viral disease will be extinct if the stochastic reproductive ratio is less than one. However, when the stochastic reproductive ratio is more than one, the viral infection system consists of an ergodic stationary distribution. Furthermore, we obtain the existence and uniqueness of the global positive solution by constructing a suitable Lyapunov function. Finally, we illustrate our results by numerical simulation.

The effect of the definition of 'pandemic' on quantitative assessments of infectious disease outbreak risk

In the early stages of an outbreak, the term 'pandemic' can be used to communicate about infectious disease risk, particularly by those who wish to encourage a large-scale public health response. However, the term lacks a widely accepted quantitative definition. We show that, under alternate quantitative definitions of 'pandemic', an epidemiological metapopulation model produces different estimates of the probability of a pandemic. Critically, we show that using different definitions alters the projected effects of key parameters-such as inter-regional travel rates, degree of pre-existing immunity, and heterogeneity in transmission rates between regions-on the risk of a pandemic. Our analysis provides a foundation for understanding the scientific importance of precise language when discussing pandemic risk, illustrating how alternative definitions affect the conclusions of modelling studies. This serves to highlight that those working on pandemic preparedness must remain alert to the variability in the use of the term 'pandemic', and provide specific quantitative definitions when undertaking one of the types of analysis that we show to be sensitive to the pandemic definition.

Will an outbreak exceed available resources for control? Estimating the risk from invading pathogens using practical definitions of a severe epidemic

Forecasting whether or not initial reports of disease will be followed by a severe epidemic is an important component of disease management. Standard epidemic risk estimates involve assuming that infections occur according to a branching process and correspond to the probability that the outbreak persists beyond the initial stochastic phase. However, an alternative assessment is to predict whether or not initial cases will lead to a severe epidemic in which available control resources are exceeded. We show how this risk can be estimated by considering three practically relevant potential definitions of a severe epidemic; namely, an outbreak in which: (i) a large number of hosts are infected simultaneously; (ii) a large total number of infections occur; and (iii) the pathogen remains in the population for a long period. We show that the probability of a severe epidemic under these definitions often coincides with the standard branching process estimate for the major epidemic probability. However, these practically relevant risk assessments can also be different from the major epidemic probability, as well as from each other. This holds in different epidemiological systems, highlighting that careful consideration of how to classify a severe epidemic is vital for accurate epidemic risk quantification.

Village-scale persistence and elimination of gambiense human African trypanosomiasis

Gambiense human African trypanosomiasis (gHAT) is one of several neglected tropical diseases that is targeted for elimination by the World Health Organization. Recent years have seen a substantial decline in the number of globally reported cases, largely driven by an intensive process of screening and treatment. However, this infection is highly focal, continuing to persist at low prevalence even in small populations. Regional elimination, and ultimately global eradication, rests on understanding the dynamics and persistence of this infection at the local population scale. Here we develop a stochastic model of gHAT dynamics, which is underpinned by screening and reporting data from one of the highest gHAT incidence regions, Kwilu Province, in the Democratic Republic of Congo. We use this model to explore the persistence of gHAT in villages of different population sizes and subject to different patterns of screening. Our models demonstrate that infection is expected to persist for long periods even in relatively small isolated populations. We further use the model to assess the risk of recrudescence following local elimination and consider how failing to detect cases during active screening events informs the probability of elimination. These quantitative results provide insights for public health policy in the region, particularly highlighting the difficulties in achieving and measuring the 2030 elimination goal.

A stochastic SIR model on a graph with epidemiological and population dynamics occurring over the same time scale

We define and study an open stochastic SIR (Susceptible-Infected-Removed) model on a graph in order to describe the spread of an epidemic on a cattle trade network with epidemiological and demographic dynamics occurring over the same time scale. Population transition intensities are assumed to be density-dependent with a constant component, the amplitude of which determines the overall scale of the population process. Standard branching approximation results for the epidemic process are first given, along with a numerical computation method for the probability of a major epidemic outbreak. This procedure is illustrated using real data on trade-related cattle movements from a densely populated livestock farming region in western France (Finistère) and epidemiological parameters corresponding to an infectious epizootic disease. Then we exhibit an exponential lower bound for the extinction time and the total size of the epidemic in the stable endemic case as a scaling parameter goes to infinity using results inspired by the Freidlin-Wentzell theory of large deviations from a dynamical system.

Precise Estimates of Persistence Time for SIS Infections in Heterogeneous Populations

For a susceptible-infectious-susceptible infection model in a heterogeneous population, we derive simple and precise estimates of mean persistence time, from a quasi-stationary endemic state to extinction of infection. Heterogeneity may be in either individuals' levels of infectiousness or of susceptibility, as well as in individuals' infectious period distributions. Infectious periods are allowed to follow arbitrary non-negative distributions. We also obtain a new and accurate approximation to the quasi-stationary distribution of the process, as well as demonstrating the use of our estimates to investigate the effects of different forms of heterogeneity. Our model may alternatively be interpreted as describing an infection spreading through a heterogeneous directed network, under the annealed network approximation.

Stochastic lattice-based modelling of malaria dynamics

BACKGROUND

The transmission of malaria is highly variable and depends on a range of climatic and anthropogenic factors. In addition, the dispersal of Anopheles mosquitoes is a key determinant that affects the persistence and dynamics of malaria. Simple, lumped-population models of malaria prevalence have been insufficient for predicting the complex responses of malaria to environmental changes.

METHODS AND RESULTS

A stochastic lattice-based model that couples a mosquito dispersal and a susceptible-exposed-infected-recovered epidemics model was developed for predicting the dynamics of malaria in heterogeneous environments. The It[Formula: see text] approximation of stochastic integrals with respect to Brownian motion was used to derive a model of stochastic differential equations. The results show that stochastic equations that capture uncertainties in the life cycle of mosquitoes and interactions among vectors, parasites, and hosts provide a mechanism for the disruptions of malaria. Finally, model simulations for a case study in the rural area of Kilifi county, Kenya are presented.

CONCLUSIONS

A stochastic lattice-based integrated malaria model has been developed. The applicability of the model for capturing the climate-driven hydrologic factors and demographic variability on malaria transmission has been demonstrated.

Persistence time of SIS infections in heterogeneous populations and networks

For a susceptible–infectious–susceptible infection model in a heterogeneous population, we present simple formulae giving the leading-order asymptotic (large population) behaviour of the mean persistence time, from an endemic state to extinction of infection. Our model may be interpreted as describing an infection spreading through either (1) a population with heterogeneity in individuals’ susceptibility and/or infectiousness; or (2) a heterogeneous directed network. Using our asymptotic formulae, we show that such heterogeneity can only reduce (to leading order) the mean persistence time compared to a corresponding homogeneous population, and that the greater the degree of heterogeneity, the more quickly infection will die out.

Modeling the within-host dynamics of cholera: bacterial-viral interaction

Novel deterministic and stochastic models are proposed in this paper for the within-host dynamics of cholera, with a focus on the bacterial-viral interaction. The deterministic model is a system of differential equations describing the interaction among the two types of vibrios and the viruses. The stochastic model is a system of Markov jump processes that is derived based on the dynamics of the deterministic model. The multitype branching process approximation is applied to estimate the extinction probability of bacteria and viruses within a human host during the early stage of the bacterial-viral infection. Accordingly, a closed-form expression is derived for the disease extinction probability, and analytic estimates are validated with numerical simulations. The local and global dynamics of the bacterial-viral interaction are analysed using the deterministic model, and the result indicates that there is a sharp disease threshold characterized by the basic reproduction number [Formula: see text]: if [Formula: see text], vibrios ingested from the environment into human body will not cause cholera infection; if [Formula: see text], vibrios will grow with increased toxicity and persist within the host, leading to human cholera. In contrast, the stochastic model indicates, more realistically, that there is always a positive probability of disease extinction within the human host.

An Alternative to Moment Closure

Moment closure methods are widely used to analyze mathematical models. They are specifically geared toward derivation of approximations of moments of stochastic models, and of similar quantities in other models. The methods possess several weaknesses: Conditions for validity of the approximations are not known, magnitudes of approximation errors are not easily evaluated, spurious solutions can be generated that require large efforts to eliminate, and expressions for the approximations are in many cases too complex to be useful. We describe an alternative method that provides improvements in these regards. The new method leads to asymptotic approximations of the first few cumulants that are explicit in the model's parameters. We analyze the univariate stochastic logistic Verhulst model and a bivariate stochastic epidemic SIR model with the new method. Errors that were made in early applications of moment closure to the Verhulst model are explained and corrected.

The quasi-stationary distribution of the closed endemic sis model

The quasi-stationary distribution of the closed stochastic SIS model changes drastically as the basic reproduction ratio R0 passes the deterministic threshold value 1. Approximations are derived that describe these changes. The quasi-stationary distribution is approximated by a geometric distribution (discrete!) for R0 distinctly below 1 and by a normal distribution (continuous!) for R0 distinctly above 1. Uniformity of the approximation with respect to R0 allows one to study the transition between these two extreme distributions. We also study the time to extinction and the invasion and persistence thresholds of the model.

The quasi-stationary distribution of the closed endemic sis model

The quasi-stationary distribution of the closed stochastic SIS model changes drastically as the basic reproduction ratio R 0 passes the deterministic threshold value 1. Approximations are derived that describe these changes. The quasi-stationary distribution is approximated by a geometric distribution (discrete!) for R 0 distinctly below 1 and by a normal distribution (continuous!) for R 0 distinctly above 1. Uniformity of the approximation with respect to R 0 allows one to study the transition between these two extreme distributions. We also study the time to extinction and the invasion and persistence thresholds of the model.

Balls, cups, and quasi-potentials: quantifying stability in stochastic systems

When a system has more than one stable state, how can the stability of these states be compared? This deceptively simple question has important consequences for ecosystems, because systems with alternative stable states can undergo dramatic regime shifts. The probability, frequency, duration, and dynamics of these shifts will all depend on the relative stability of the stable states. Unfortunately, the concept of “stability” in ecology has suffered from substantial confusion and this is particularly problematic for systems where stochastic perturbations can cause shifts between coexisting alternative stable states. A useful way to visualize stable states in stochastic systems is with a ball‐in‐cup‐diagram, in which the state of the system is represented as the position of a ball rolling on a surface, and the random perturbations can push the ball from one basin of attraction to another. The surface is determined by a potential function, which provides a natural stability metric. Systems amenable to this representation, called gradient systems, are quite rare, however. As a result, the potential function is not widely used and other approaches based on linear stability analysis have become standard. Linear stability analysis is designed for local analysis of deterministic systems and, as we show, can produce a highly misleading picture of how the system will behave under continual, stochastic perturbations. In this paper, we show how the potential function can be generalized so that it can be applied broadly, employing a concept from stochastic analysis called the quasi‐potential. Using three classic ecological models, we demonstrate that the quasi‐potential provides a useful way to quantify stability in stochastic systems. We show that the quasi‐potential framework helps clarify long‐standing confusion about stability in stochastic ecological systems, and we argue that ecologists should adopt it as a practical tool for analyzing these systems.

The probability of epidemic fade-out is non-monotonic in transmission rate for the Markovian SIR model with demography

Epidemic fade-out refers to infection elimination in the trough between the first and second waves of an outbreak. The number of infectious individuals drops to a relatively low level between these waves of infection, and if elimination does not occur at this stage, then the disease is likely to become endemic. For this reason, it appears to be an ideal target for control efforts. Despite this obvious public health importance, the probability of epidemic fade-out is not well understood. Here we present new algorithms for approximating the probability of epidemic fade-out for the Markovian SIR model with demography. These algorithms are more accurate than previously published formulae, and one of them scales well to large population sizes. This method allows us to investigate the probability of epidemic fade-out as a function of the effective transmission rate, recovery rate, population turnover rate, and population size. We identify an interesting feature: the probability of epidemic fade-out is very often greatest when the basic reproduction number, R0, is approximately 2 (restricting consideration to cases where a major outbreak is possible, i.e., R0>1). The public health implication is that there may be instances where a non-lethal infection should be allowed to spread, or antiviral usage should be moderated, to maximise the chance of the infection being eliminated before it becomes endemic.

Stochastic tunneling and metastable states during the somatic evolution of cancer

Tumors initiate when a population of proliferating cells accumulates a certain number and type of genetic and/or epigenetic alterations. The population dynamics of such sequential acquisition of (epi)genetic alterations has been the topic of much investigation. The phenomenon of stochastic tunneling, where an intermediate mutant in a sequence does not reach fixation in a population before generating a double mutant, has been studied using a variety of computational and mathematical methods. However, the field still lacks a comprehensive analytical description since theoretical predictions of fixation times are available only for cases in which the second mutant is advantageous. Here, we study stochastic tunneling in a Moran model. Analyzing the deterministic dynamics of large populations we systematically identify the parameter regimes captured by existing approaches. Our analysis also reveals fitness landscapes and mutation rates for which finite populations are found in long-lived metastable states. These are landscapes in which the final mutant is not the most advantageous in the sequence, and resulting metastable states are a consequence of a mutation-selection balance. The escape from these states is driven by intrinsic noise, and their location affects the probability of tunneling. Existing methods no longer apply. In these regimes it is the escape from the metastable states that is the key bottleneck; fixation is no longer limited by the emergence of a successful mutant lineage. We used the so-called Wentzel-Kramers-Brillouin method to compute fixation times in these parameter regimes, successfully validated by stochastic simulations. Our work fills a gap left by previous approaches and provides a more comprehensive description of the acquisition of multiple mutations in populations of somatic cells.

Noise-induced switching and extinction in systems with delay

We consider the rates of noise-induced switching between the stable states of dissipative dynamical systems with delay and also the rates of noise-induced extinction, where such systems model population dynamics. We study a class of systems where the evolution depends on the dynamical variables at a preceding time with a fixed time delay, which we call hard delay. For weak noise, the rates of interattractor switching and extinction are exponentially small. Finding these rates to logarithmic accuracy is reduced to variational problems. The solutions of the variational problems give the most probable paths followed in switching or extinction. We show that the equations for the most probable paths are acausal and formulate the appropriate boundary conditions. Explicit results are obtained for small delay compared to the relaxation rate. We also develop a direct variational method to find the rates. We find that the analytical results agree well with the numerical simulations for both switching and extinction rates.

Two-strain competition in quasineutral stochastic disease dynamics

We develop a perturbation method for studying quasineutral competition in a broad class of stochastic competition models and apply it to the analysis of fixation of competing strains in two epidemic models. The first model is a two-strain generalization of the stochastic susceptible-infected-susceptible (SIS) model. Here we extend previous results due to Parsons and Quince [Theor. Popul. Biol. 72, 468 (2007)], Parsons et al. [Theor. Popul. Biol. 74, 302 (2008)], and Lin, Kim, and Doering [J. Stat. Phys. 148, 646 (2012)]. The second model, a two-strain generalization of the stochastic susceptible-infected-recovered (SIR) model with population turnover, has not been studied previously. In each of the two models, when the basic reproduction numbers of the two strains are identical, a system with an infinite population size approaches a point on the deterministic coexistence line (CL): a straight line of fixed points in the phase space of subpopulation sizes. Shot noise drives one of the strain populations to fixation, and the other to extinction, on a time scale proportional to the total population size. Our perturbation method explicitly tracks the dynamics of the probability distribution of the subpopulations in the vicinity of the CL. We argue that, whereas the slow strain has a competitive advantage for mathematically "typical" initial conditions, it is the fast strain that is more likely to win in the important situation when a few infectives of both strains are introduced into a susceptible population.

Estimating the probability of an extinction or major outbreak for an environmentally transmitted infectious disease

Indirect transmission through the environment, pathogen shedding by infectious hosts, replication of free-living pathogens within the environment, and environmental decontamination are suspected to play important roles in the spread and control of environmentally transmitted infectious diseases. To account for these factors, the classic Susceptible-Infectious-Recovered-Susceptible epidemic model is modified to include a compartment representing the amount of free-living pathogen within the environment. The model accounts for host demography, direct and indirect transmission, replication of free-living pathogens in the environment, and removal of free-living pathogens by natural death or environmental decontamination. Based on the assumptions of the deterministic model, a continuous-time Markov chain model is developed. An estimate for the probability of disease extinction or a major outbreak is obtained by approximating the Markov chain with a multitype branching process. Numerical simulations illustrate important differences between the deterministic and stochastic counterparts, relevant for outbreak prevention, that depend on indirect transmission, pathogen shedding by infectious hosts, replication of free-living pathogens, and environmental decontamination. The probability of a major outbreak is computed for salmonellosis in a herd of dairy cattle as well as cholera in a human population. An explicit expression for the probability of disease extinction or a major outbreak in terms of the model parameters is obtained for systems with no direct transmission or replication of free-living pathogens.

Dynamics of infectious diseases

Modern infectious disease epidemiology has a strong history of using mathematics both for prediction and to gain a deeper understanding. However the study of infectious diseases is a highly interdisciplinary subject requiring insights from multiple disciplines, in particular a biological knowledge of the pathogen, a statistical description of the available data and a mathematical framework for prediction. Here we begin with the basic building blocks of infectious disease epidemiology--the SIS and SIR type models--before considering the progress that has been made over the recent decades and the challenges that lie ahead. Throughout we focus on the understanding that can be developed from relatively simple models, although accurate prediction will inevitably require far greater complexity beyond the scope of this review. In particular, we focus on three critical aspects of infectious disease models that we feel fundamentally shape their dynamics: heterogeneously structured populations, stochasticity and spatial structure. Throughout we relate the mathematical models and their results to a variety of real-world problems.

Integrated Sensor Systems and Data Fusion for Homeland Protection

This chapter addresses the application of data and information fusion to the design of integrated systems in the Homeland Protection (HP) domain. HP is a wide and complex domain: systems in this domain are large (in terms of size and scope) integrated (each subsystem cannot be considered as an isolated system) and different in purpose. Such systems require a multidisciplinary approach for their design and analysis and they are necessarily required to provide data and information fusion in the most general sense. The first data fusion algorithms employed in real systems in the radar field go back to the early seventies; now a days new concepts have been developed and spread to be applied to very complex systems with the aim to achieve the highest level of intelligence as possible and hopefully to support decision. Data fusion is aimed to enhance situation awareness and decision making through the combination of information/data obtained by networks of homogeneous and/or heterogeneous sensors. The aim of this chapter is to give an overview of the several approaches that can be followed to design and analyze systems for homeland protection. Different fusion architectures can be drawn on the basis of the employed algorithms: they are analyzed under several aspects in this chapter. Real study cases applied to real world problems of homeland protection are provided in the chapter.

Multiple extinction routes in stochastic population models

Isolated populations ultimately go extinct because of the intrinsic noise of elementary processes. In multipopulation systems extinction of a population may occur via more than one route. We investigate this generic situation in a simple predator-prey (or infected-susceptible) model. The predator and prey populations may coexist for a long time, but ultimately both go extinct. In the first extinction route the predators go extinct first, whereas the prey thrive for a long time and then also go extinct. In the second route the prey go extinct first, causing a rapid extinction of the predators. Assuming large subpopulation sizes in the coexistence state, we compare the probabilities of each of the two extinction routes and predict the most likely path of the subpopulations to extinction. We also suggest an effective three-state master equation for the probabilities to observe the coexistence state, the predator-free state, and the empty state.

Switching between phenotypes and population extinction

Many types of bacteria can survive under stress by switching stochastically between two different phenotypes: the "normals" who multiply fast, but are vulnerable to stress, and the "persisters" who hardly multiply, but are resilient to stress. Previous theoretical studies of such bacterial populations have focused on the fitness: the asymptotic rate of unbounded growth of the population. Yet for an isolated population of established (and not very large) size, a more relevant measure may be the population extinction risk due to the interplay of adverse extrinsic variations and intrinsic noise of birth, death and switching processes. Applying a WKB approximation to the pertinent master equation of such a two-population system, we quantify the extinction risk, and find the most likely path to extinction under both favorable and adverse conditions. Analytical results are obtained both in the biologically relevant regime when the switching is rare compared with the birth and death processes, and in the opposite regime of frequent switching. We show that rare switches are most beneficial in reducing the extinction risk.

Control of rare events in reaction and population systems by deterministically imposed transitions

We consider control of reaction and population systems by imposing transitions between states with different numbers of particles or individuals. The transitions take place at predetermined instants of time. Even where they are significantly less frequent than spontaneous transitions, they can exponentially strongly modify the rates of rare events, including switching between metastable states or population extinction. We also study optimal control of rare events. Specifically, we are interested in the optimal control of disease extinction for a limited vaccine supply. A comparison is made with control of rare events by modulating the rates of elementary transitions rather than imposing transitions deterministically. It is found that, unexpectedly, for the same mean control parameters, controlling the transitions rates can be more efficient.

Stochastic dynamics of cholera epidemics

We describe the predictions of an analytically tractable stochastic model for cholera epidemics following a single initial outbreak. The exact model relies on a set of assumptions that may restrict the generality of the approach and yet provides a realm of powerful tools and results. Without resorting to the depletion of susceptible individuals, as usually assumed in deterministic susceptible-infected-recovered models, we show that a simple stochastic equation for the number of ill individuals provides a mechanism for the decay of the epidemics occurring on the typical time scale of seasonality. The model is shown to provide a reasonably accurate description of the empirical data of the 2000/2001 cholera epidemic which took place in the Kwa Zulu-Natal Province, South Africa, with possibly notable epidemiological implications.

Speeding up disease extinction with a limited amount of vaccine

We consider optimal vaccination protocol where the vaccine is in short supply. In this case, the endemic state remains dynamically stable; disease extinction happens at random and requires a large fluctuation, which can come from the intrinsic randomness of the population dynamics. We show that vaccination can exponentially increase the disease extinction rate. For a time-periodic vaccination with fixed average rate, the optimal vaccination protocol is model independent and presents a sequence of short pulses. The effect can be resonantly enhanced if the vaccination pulse period coincides with the characteristic period of the disease dynamics or its multiples. This resonant effect is illustrated using a simple epidemic model. The analysis is based on the theory of fluctuation-induced population extinction in periodically modulated systems that we develop. If the system is strongly modulated (for example, by seasonal variations) and vaccination has the same period, the vaccination pulses must be properly synchronized; a wrong vaccination phase can impede disease extinction.

Enhanced vaccine control of epidemics in adaptive networks

We study vaccine control for disease spread on an adaptive network modeling disease avoidance behavior. Control is implemented by adding Poisson-distributed vaccination of susceptibles. We show that vaccine control is much more effective in adaptive networks than in static networks due to feedback interaction between the adaptive network rewiring and the vaccine application. When compared to extinction rates in static social networks, we find that the amount of vaccine resources required to sustain similar rates of extinction are as much as two orders of magnitude lower in adaptive networks.

Extinction rate fragility in population dynamics

Population extinction is of central interest for population dynamics. It may occur from a large rare fluctuation. We find that, in contrast to related large-fluctuation effects like noise-induced interstate switching, quite generally extinction rates in multipopulation systems display fragility, where the height of the effective barrier to be overcome in the fluctuation depends on the system parameters nonanalytically. We show that one of the best-known models of epidemiology, the susceptible-infectious-susceptible model, is fragile to total population fluctuations.

Switching-path distribution in multidimensional systems

We explore the distribution of paths followed in fluctuation-induced switching between coexisting stable states. We introduce a quantitative characteristic of the path distribution in phase space that does not require a priori knowledge of system dynamics. The theory of the distribution is developed and its direct measurement is performed in a micromechanical oscillator driven into parametric resonance. The experimental and theoretical results on the shape and position of the distribution are in excellent agreement, with no adjustable parameters. In addition, the experiment provides the first demonstration of the lack of time-reversal symmetry in switching of systems far from thermal equilibrium. The results open the possibility of efficient control of the switching probability based on the measured narrow path distribution.

Disease extinction in the presence of random vaccination

We investigate disease extinction in an epidemic model described by a birth-death process. We show that, in the absence of vaccination, the effective entropic barrier for extinction displays scaling with the distance to the bifurcation point, with an unusual critical exponent. Even a comparatively weak Poisson-distributed random vaccination leads to an exponential increase in the extinction rate, with the exponent that strongly depends on the vaccination parameters.

Extinction of an infectious disease: a large fluctuation in a nonequilibrium system

We develop a theory of first passage processes in stochastic nonequilibrium systems of birth-death type using two closely related epidemiological models as examples. Our method employs the probability generating function technique in conjunction with the eikonal approximation. In this way the problem is reduced to finding the optimal path to extinction: a heteroclinic trajectory of an effective multidimensional classical Hamiltonian system. We compute this trajectory and mean extinction time of the disease numerically and uncover a nonmonotone, spiral path to extinction of a disease. We also obtain analytical results close to a bifurcation point, where the problem is described by a Hamiltonian previously identified in one-species population models.

Endemic persistence or disease extinction: The effect of separation into sub-communities

Consider an infectious disease which is endemic in a population divided into several large sub-communities that interact. Our aim is to understand how the time to extinction is affected by the level of interaction between communities. We present two approximations of the expected time to extinction in a population consisting of a small number of large sub-communities. These approximations are described for an SIR epidemic model, with focus on diseases with short infectious period in relation to life length, such as childhood diseases. Both approximations are based on Markov jump processes. Simulations indicate that the time to extinction is increasing in the degree of interaction between communities. This behaviour can also be seen in our approximations in relevant regions of the parameter space.

Spatial heterogeneity and the persistence of infectious diseases

The endemic persistence of infectious diseases can often not be understood without taking into account the relevant heterogeneities of host mixing. Here, we consider spatial heterogeneity, defined as 'patchiness' of the host population. After briefly reviewing how disease persistence is influenced by population size, reproduction number and infectious period, we explore its dependence on the level of spatial heterogeneity. Analysis and simulation of disease transmission in a symmetric meta-population suggest that disease persistence typically becomes worse as spatial heterogeneity increases, although local persistence optima can occur for infections with oscillatory population dynamics. We obtain insight into the dynamics that underlie the observed persistence patterns by studying the infection prevalence correlation between patches and by comparing full-model simulations to results obtained using simplified patch-level descriptions of the interplay between local extinctions and between-patch transmissions. The observed patterns are interpreted in terms of rescue effects for strong spatial heterogeneity and in terms of between-patch coherence and synchronization effects at intermediate and weak levels of heterogeneity.

Measures of durability of resistance

ABSTRACT Conventional models for the durability of resistant cultivars focus on the dynamics of the frequency of resistance genes. This leads to a definition of the durability of resistance as the time from introduction of the cultivar to the time when the frequency of the virulence gene reaches a preset threshold. It is questionable whether this is the most appropriate way to measure durability. Here we use a simple epidemiological model to link population dynamics and population genetics to compare three measures of durability: (i) the expected time until invasion of the virulent genotype, by mutation or immigration, and subsequent establishment of a population (T(invasion)); (ii) the virulence frequency related measure of the time for the virulent genotype to take-over the pathogen population ( T(take-over)); and (iii) the additional yield, measured by the additional number of uninfected host growth days (T(additional)). Specifically, we show how the measures of durability are affected by deployment and epidemiological parameters. We use a combination of numerical solution and analytical approximation of a model for the population dynamics of avirulent and virulent genotypes of a pathogen growing in dynamically changing populations of resistant and susceptible cultivars. The three measures of durability are compared. Some consequences of the results for durable resistance in multilines and mixtures and the regional deployment of resistant cultivars are discussed.

Realistic distributions of infectious periods in epidemic models: changing patterns of persistence and dynamics

Most mathematical models used to study the epidemiology of childhood viral diseases, such as measles, describe the period of infectiousness by an exponential distribution. The effects of including more realistic descriptions of the infectious period within SIR (susceptible/infectious/recovered) models are studied. Less dispersed distributions are seen to have two important epidemiological consequences. First, less stable behaviour is seen within the model: incidence patterns become more complex. Second, disease persistence is diminished: in models with a finite population, the minimum population size needed to allow disease persistence increases. The assumption made concerning the infectious period distribution is of a kind routinely made in the formulation of mathematical models in population biology. Since it has a major effect on the central issues of population persistence and dynamics, the results of this study have broad implications for mathematical modellers of a wide range of biological systems.

Sustained oscillations in stochastic systems

Many non-linear deterministic models for interacting populations present damped oscillations towards the corresponding equilibrium values. However, simulations produced with related stochastic models usually present sustained oscillations which preserve the natural frequency of the damped oscillations of the deterministic model but showing non-vanishing amplitudes. The relation between the amplitude of the stochastic oscillations and the values of the equilibrium populations is not intuitive in general but scales with the square root of the populations when the ratio between different populations is kept fixed. In this work, we explain such phenomena for the case of a general epidemic model. We estimate the stochastic fluctuations of the populations around the equilibrium point in the epidemiological model showing their (approximated) relation with the mean values.

Stochastic epidemics: the expected duration of the endemic period in higher dimensional models

A method is presented to approximate the long-term stochastic dynamics of an epidemic modelled by state variables denoting the various classes of the population such as in SIR and SEIR model. The modelling includes epidemics in populations at different locations with migration between these populations. A logistic stochastic process for the total infectious population is formulated; it fits the long-term stochastic behaviour of the total infectious population in the full model. A good approximation is obtained if only the dynamics near the equilibria is fit.

Implications derived from a mathematical model for eradication of pseudorabies virus

Simple mathematical models based on experimental and observational data were applied to evaluate the feasibility of eradicating pseudorabies virus (PRV) regionally by vaccination and to determine which factors can jeopardise eradication. As much as possible, the models were uncomplicated and our conclusions were based on mathematical analysis. For complicated situations, Monte-Carlo simulation was used to support the conclusions. For eradication, it is sufficient that the reproduction ratio R (the number of units infected by one infectious unit) is < 1. However, R can be determined at different scales: at one end the region with the herds as units and at the other end compartments with the pigs as units. Results from modelling within herds showed that contacts between groups within a herd is important whenever R between individuals (R(ind)) is > 1 in one or more groups. This is the case within finishing herds. In addition, if the R(ind) is more than 1 within a herd, the size of the herd determines whether PRV can persist in the herd and determines the duration of persistence. Moreover, when reactivation of PRV in well-vaccinated sows is taken into account, R(ind) in sow herds is still less than 1. In sow herds with group-housing systems, it is possible that in those groups R(ind) is > 1. Results from modelling between herds showed that whether or not Rherd is < 1 in a particular region is determined by two factors: (1) the transmission of infection between nucleus herds and rearing herds through transfer of animals and (2) contacts among finishing herds and among rearing herds. The transmission between herds can be reduced by reduction of the contact rate between herds. reduction of the herd size, and reduction of the transmission within herds.