key: cord-0858170-2g6jxi5k
authors: Safi, Mohammad A.; Gumel, Abba B.
title: Dynamics of a model with quarantine-adjusted incidence and quarantine of susceptible individuals
date: 2013-03-15
journal: J Math Anal Appl
DOI: 10.1016/j.jmaa.2012.10.015
sha: 0c3b3ae35af5de4b137bc3b40d6bf707c6b43cbc
doc_id: 858170
cord_uid: 2g6jxi5k

A new deterministic model for the spread of a communicable disease that is controllable using mass quarantine is designed. Unlike in the case of the vast majority of prior quarantine models in the literature, the new model includes a quarantine-adjusted incidence function for the infection rate and the quarantine of susceptible individuals suspected of being exposed to the disease (thereby making it more realistic epidemiologically). The earlier quarantine models tend to only explicitly consider individuals who are already infected, but show no clinical symptoms of the disease (i.e., those latently-infected), in the quarantine class (while ignoring the quarantine of susceptible individuals). In reality, however, the vast majority of people in quarantine (during a disease outbreak) are susceptible. Rigorous analysis of the model shows that the assumed imperfect nature of quarantine (in preventing the infection of quarantined susceptible individuals) induces the phenomenon of backward bifurcation when the associated reproduction threshold is less than unity (thereby making effective disease control difficult). For the case when the efficacy of quarantine to prevent infection during quarantine is perfect, the disease-free equilibrium is globally-asymptotically stable when the reproduction threshold is less than unity. Furthermore, the model has a unique endemic equilibrium when the reproduction threshold exceeds unity (and the disease persists in the population in this case).

Quarantine (of individuals feared exposed to a communicable disease) is one of the oldest public health control measures for the spread of communicable diseases in given populations. More recently, this measure was successfully used to combat the spread of some emerging and re-emerging human and animal diseases, such as the severe acute respiratory syndrome (SARS) [1] [2] [3] [4] [5] [6] [7] [8] , foot-and-mouth disease [9] and the 2009 swine influenza pandemic [10] . Mathematical models have been designed and used to gain qualitative insights into the spread of diseases in the presence of quarantine of individuals suspected of being exposed to a disease, and the subsequent isolation or hospitalization of those with clinical symptoms of the disease (see, for instance, [1] [2] [3] [4] [5] [6] [7] [8] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] ). In the aforementioned models, with the exception of the models in [2, 5] , the dynamics of the quarantined susceptible individuals is not explicitly incorporated into the model (the models only count quarantined individuals who are actually latently-infected, and essentially ignore those who remain susceptible at the end of the quarantine period).

In other words, in the aforementioned models (with the exception of the models in [2, 5] ) the class (or compartment) of quarantined individuals contain only those who are latently-infected (and the temporary removal of susceptible individuals from the susceptible pool is not accounted for). The justification given for this assumption (of ignoring the back-andforth dynamics of quarantined susceptible individuals: from susceptible to quarantine and back to susceptible class if they show no disease symptoms) is that, in general, the fraction of infected contacts that can be traced and quarantined at the time of infection is very small; and that the total population is large in comparison to the size of the infected individuals [13] . In practice, however, quarantine involves the removal of all individuals suspected of being exposed to a given disease (regardless of infection status) from the rest of the population (this often involves a lot of people; the term mass quarantine is often used to describe this process). Those who show disease symptoms during quarantine are hospitalized (or isolated), and those who remain susceptible at the end of the quarantine period are returned to the susceptible class. It is, therefore, instructive to study the impact of the aforementioned assumption (of not explicitly accounting for the quarantine of susceptible individuals) on the transmission dynamics of a disease that is controllable using quarantine.

The purpose of this study is to investigate the qualitative impact of explicitly including the dynamics of quarantined susceptible individuals on the spread of a disease that is controllable by quarantine and isolation. To achieve this objective, a new deterministic model, which includes the dynamics of quarantined susceptible individuals and quarantine-adjusted incidence function, is designed and analyzed. It should be mentioned that quarantined-adjusted incidence has also been used in a number of quarantine models, such as those in [5, 14, 15, 20] .

The model to be considered is that for the transmission dynamics of a communicable disease that can be controlled using quarantine and isolation, such as ebola, measles, pandemic influenza and SARS [1] [2] [3] [4] [5] 7, 8, 11, 12, 14, 19, 21] . It is based on splitting the total population at time t, denoted by N(t), into the sub-populations of non-quarantined susceptible (S(t)), quarantined susceptible (S q (t)), non-quarantined latently-infected (i.e., infected but show no clinical symptoms of the disease) (E(t)), quarantined latently-infected (E q (t)), non-quarantined symptomatic (I(t)), quarantined symptomatic (I q (t)), hospitalized (H(t)) and recovered (R(t)) individuals, so that

The model is given by the following system of non-linear differential equations (a flow diagram of the model is depicted in Fig. 1) :

where λ a is the infection rate given by Safi et al. [20] :

and N a (t) is the total actively-mixing population given by (see also [20] )

In (2), β is the effective contact rate, and the modification parameters 0 ≤ η 1 , η h < 1 accounts for the assumed reduction of infectiousness of exposed and hospitalized individuals, respectively, in relation to the symptomatically-infected (infectious) individuals in the I class. Similarly, 0 ≤ η q ≤ 1 accounts for the assumed reduction of infectiousness of infected quarantined individuals in relation to individuals in the I class (the parameter 0 ≤ η 2 ≤ 1 represents the assumed reduction of infectiousness of exposed quarantined individuals, in the E q class, in relation to infectious quarantined individuals in the I q class). In (3), η q is a modification parameter that measures the efficacy of quarantine in preventing individuals in quarantine from having contact with members of the general public. If η q = 0, then quarantine is perfectly-implemented (so that individuals in the quarantine classes are not part of the actively-mixing population, and do not transmit infection).

The non-quarantined susceptible population (S) is increased by the recruitment of individuals into the community, at a rate Π. This population is reduced by infection (at the rate λ a ), quarantine (at a rate γ ) and natural death (at a rate µ; populations in all epidemiological compartments are assumed to suffer natural death at this rate). This population is increased by the return of quarantined susceptible individuals at the end of the quarantine period (at a rate ψ 1 ) and the loss of natural immunity by recovered individuals (at a rate ψ 2 ). The population of quarantined susceptible individuals (S q ) is generated by the quarantine of non-quarantined susceptible individuals (at the rate γ ). It is decreased by infection (at a reduced rate rλ a , with 0 < r < 1 accounting for the efficacy of quarantine in preventing infection of quarantined susceptible individuals) and by the reversion to the non-quarantined susceptible class at the end of the quarantine period (at the rate ψ 1 ). Non-quarantined latently-infected individuals (E) are generated at the rate λ a and decreased by the development of clinical symptoms of the disease (at a rate σ 1 ). Similarly, the population of quarantined latently-infected individuals (E q ) is generated at the rate rλ a and decreased by the development of clinical symptoms (at a rate σ 2 ). Non-quarantined symptomatic individuals (I) are generated at the rate σ 1 . This population is decreased by hospitalization (at a rate α 1 ), natural recovery (at a rate φ 1 ) and disease-induced mortality (at a rate δ 1 ). The population of quarantined symptomatic individuals (I q ) is generated at the rate σ 2 and decreased by hospitalization (at a rate α 2 ), recovery (at a rate φ 2 ) and disease-induced death (at a rate δ 2 ).

The population of hospitalized individuals (H) is generated by the hospitalization of non-quarantined and quarantined symptomatic individuals (at the rate α 1 and α 2 , respectively). It is diminished by recovery (at a rate φ 3 ) and disease-induced death (at a rate δ 3 < δ 2 < δ 1 ). The recovered population (R) is generated at the rates φ i (i = 1, 2, 3) and diminished by loss of natural immunity (at the rate ψ 2 ).

The model (1) is an extension of many of the models for quarantine published in the literature (including those in [1, 3, 4, 6, [11] [12] [13] [14] [15] [16] [17] [18] ), by considering the dynamics of quarantined susceptible individuals (this entails adding the epidemiological compartments S q , E q and I q ), in line with the models in [2, 5] . Furthermore, the model (1) extends the model in [5] by:

(i) including a compartment for hospitalized individuals (H); (ii) allowing for the infection of quarantined susceptible individuals (at the rate rλ a ; that is, unlike in [5] , it is assumed that the efficacy of quarantine to prevent infection of quarantined susceptible individuals is not 100%); (iii) allowing for disease transmission by infected quarantined and hospitalized individuals (i.e., it is assumed, unlike in [5] , that quarantine and hospitalization are not 100% effective in preventing transmission by infected quarantined or hospitalized individuals).

The model (1) also extends the model in [2] by:

(a) incorporating demographic parameters (i.e., using an endemic model as against the epidemic model used in [2] ); (b) using a standard incidence function to model the infection rate (the mass action incidence function was used in [2] by assuming constant total population); (c) allowing for the infection of quarantined susceptible individuals (quarantine was assumed to be perfect against infection in [2] ); (d) allowing for disease transmission by non-quarantined latently-infected individuals.

Furthermore, detailed qualitative analysis of the model will be carried out in this study (this is not done in [2] ). The model (1) extends the quarantine models with quarantined-adjusted incidence in [15, 20] , by including the dynamics of quarantined susceptible individuals (S q ) as well as adding classes for quarantined exposed (E q ) and symptomatic (I q ) individuals.

The model (1) will now be analyzed to gain insight into its dynamical features.

Since the model (1) monitors human populations, all its associated parameters are non-negative. The following basic results can be established using the method in Appendix A of [22] . 

Proof. Adding all the equations of the model (1) gives,

Since dN/dt ≤ Π − µN, it follows that dN/dt ≤ 0 if N ≥ Π/µ. Thus, a standard comparison theorem [23] can be used 

The DFE of the model (1) is given by

The local stability of E 0 will be explored using the next generation operator method [24, 25] . Using the notation in [25] , the non-negative matrix, F , of the new infection terms, and the M-matrix, V , of the transition terms associated with the model (1), are given, respectively, by

where,

It follows that the control reproduction number [26, 27] , denoted by R q = ρ(FV −1 ), where ρ is the spectral radius, is given by

Using Theorem 2 in [25] , the following result is established.

The DFE of the model (1), given by (4), is locally-asymptotically stable (LAS) if R q < 1, and unstable if R q > 1.

The threshold quantity R q measures the average number of new infections generated by a single infectious individual in a population where a certain fraction of the susceptible population is vaccinated. Lemma 2 implies that the disease can be eliminated from the community (when R q < 1) if the initial sizes of the sub-populations of the model are in the basin of attraction of the DFE (E 0 ).

In this section, the existence of endemic equilibria (that is, equilibria where the infected compartments of the model are non-zero) of the model (1) is established. Let,

represents any arbitrary endemic equilibrium point (EEP) of the model (1). Further, define

(the force of infection of the model (1) at steady-state). It follows, by solving the equations in (1) at steady-state, that

Substituting the expressions in (6) into (5) shows that the non-zero equilibria of the model satisfy the following quadratic equation (in terms of λ * * a ):

where,

The endemic equilibria of the model (1) can then be obtained by solving for λ * * a from (7), and substituting the positive values of λ * * a into the expressions in (6) . The quadratic equation (7) can be analyzed for the possibility of multiple endemic equilibria when R q < 1. It should be noted that the coefficient, a 0 , of the quadratic (7) is always positive and a 2 is positive (negative) if R q is less (greater) than unity. Hence, the following result is established.

(ii) a unique endemic equilibrium if (a 1 < 0 and a 2 = 0) or a 2 1 − 4a 0 a 2 = 0; (iii) two endemic equilibria if a 2 > 0, a 1 < 0 and a 2 1 − 4a 0 a 2 > 0;

(iv) no endemic equilibrium otherwise.

Thus, it is clear from Case (i) of Theorem 2 that the model (1) has a unique EEP (of the form E 1 ) whenever R q > 1.

Furthermore, Case (iii) of Theorem 2 indicates the possibility of backward bifurcation, where a LAS DFE co-exists with a LAS endemic equilibrium when the associated reproduction number R q is less than unity (see, for instance, [28] [29] [30] for discussions on backward bifurcation) in the model (1). The epidemiological importance of the phenomenon of backward bifurcation is that the classical requirement of having R q < 1 is, although necessary, not sufficient for disease elimination.

In this case, disease elimination will depend upon the initial sizes of the sub-populations of the model. Thus, the existence of backward bifurcation in the transmission dynamics of a disease makes its effective control difficult. A more rigorous proof of the backward bifurcation property of the model (1), based on using center manifold theory (see, for instance, [25, 29, 31, 32] ), is given in the Appendix.

In this section, the scenario where the backward bifurcation property of the model can be lost is explored. Consider the model (1) with a perfect quarantine efficacy against infection (so that, r = 0). In this case, the coefficients a 0 , a 1 and a 2 of the quadratic equation (7) reduce to a 0 = 0, a 1 > 0 and a 2 ≥ 0 wheneverR q = R q | r=0 ≤ 1. Thus, for this case, the quadratic equation (7) has one solution (λ * * a = −a 2 a 1 ≤ 0.) Therefore, the model (1) with a perfect quarantine has no positive endemic equilibrium wheneverR q < 1. This rules out the possibility of backward bifurcation in this case (since backward bifurcation requires the existence of at least two endemic equilibria wheneverR q ≤ 1 [28] [29] [30] ). Furthermore, it can be shown that, for the case when r = 0, the DFE (E 0 ) of the model (1) is globally-asymptotically stable (GAS) under some conditions, as shown below.

Setting r = 0 in the model (1) gives the following reduced model (it should be noted from (1) that, for the case when r = 0, (E q (t), I q (t)) → (0, 0) as t → ∞; hence, these variables are omitted from the asymptotic analysis of the model for the special case with r = 0):

with the associated force of infection λ a = λ 1 , where

It can be shown that the reproduction number associated with the reduced model (8) , with (9), is given bỹ

Define,

The model (8) has a DFE, given by E 01 = (S * , S * q , 0, 0, 0, 0).

The DFE (E 01 ) of the reduced model (8) , with (9), is GAS in D 1 wheneverR q ≤ ν 1 < 1.

with Lyapunov derivative (where a dot represents differentiation with respect to time) given bẏ

Since all the parameters and variables of the model (1) are non-negative (Theorem 1), it follows thatḞ ≤ 0 forR q ≤ ν 1 (it should be noted that ν 1 = S * N * < 1) withḞ = 0 if and only if E = I = H = 0. Hence, F is a Lyapunov function on D 1 . Thus, it follows, by LaSalle's Invariance Principle [33] , that (E(t), I(t), H(t)) → (0, 0, 0) as t → ∞.

Since lim sup t→∞ I(t) = 0 and lim sup t→∞ H(t) = 0 (from (10)), it follows that, for sufficiently small ϖ * > 0, there exist constants M 1 > 0 and M 2 > 0 such that lim sup t→∞ I(t) ≤ ϖ * for all t > M 1 and lim sup t→∞ H(t) ≤ ϖ * for all t > M 2 . Hence, it follows from the last equation of the model (8) that, for t > max{M 1 , M 2 },

Thus, by comparison theorem [23] ,

so that, by letting ϖ * → 0,

Similarly (by using lim inf t→∞ I(t) = 0 and lim inf t→∞ H(t) = 0), it can be shown that

Thus, it follows from (11) and (12) that

Similarly, it can be shown that

Thus, by combining Eqs. (10), (13) and (14), it follows that every solution of the equations of the model (8), with initial conditions in D 1 , approaches E 0 as t → ∞ (forR q ≤ ν 1 < 1).

The above result shows that, for the case when the efficacy of quarantine in preventing infection is perfect (i.e., r = 0), the disease can be eliminated from the community if the associated reproduction number of the model is less than unity.

Furthermore, this result clearly shows that the backward bifurcation property of the model (1) is caused by the imperfect nature of the efficacy of quarantine to prevent infection (see the Appendix).

Conditions for the persistence of the disease in the population will be investigated below.

The model system (1) is said to be uniformly-persistent if there exists a constant c such that any solution (S(t), S q (t), E(t), E q (t), I(t), I q (t), H(t), R(t)) satisfies [34, 35] :

Proof. The theorem can be proved by using the approach used to prove Proposition 3.3 of [36] , by applying a uniform persistence result in [34] and noting that the DFE of the model (1) is unstable whenever R q > 1 (Lemma 2).

The consequence of this result is that each infected variable (E, E q , I, I q , H) of the model will persist above a certain threshold value at steady-state, so that the disease will persist (become endemic) in the population. It should be mentioned that some other models for quarantine, which did not explicitly include the dynamics of quarantined susceptible individuals, also showed the presence of a unique endemic equilibrium and disease persistence when the associated reproduction number exceeds unity (see, for instance, [13, 17, 18] ).

A new deterministic model for the spread of a disease in a population in the presence of quarantine is designed. A major feature of the model is that it incorporate the dynamics of quarantine-adjusted incidence and the quarantine of susceptible individuals (that is, quarantine is modeled in terms of the temporarily removal of susceptible individuals from the susceptible pool as well as the removal of new infected individuals, detected via the contact tracing of known infectious individuals). Rigorous analysis of the model reveals that it undergoes the phenomenon of backward bifurcation when the associated reproduction number (R q ) is less than unity. The presence of this phenomenon, which does not arise if the quarantine is 100% effective, implies that the effective control of the spread of the disease, using an imperfect quarantine, depends on the initial sizes of the sub-populations of the model (when R q < 1). The model has a unique endemic equilibrium whenever R q > 1. Furthermore, it is shown that the disease will persist in the population whenever R q > 1.

ABG acknowledges, with thanks, the support in part of the Natural Science and Engineering Research Council (NSERC) and Mathematics of Information Technology and Complex Systems (MITACS) of Canada. The authors are grateful to the anonymous reviewers for their constructive comments.

The proof is based on using center manifold theory [31, 32] . In particular, Theorem 4.1 of [32] will be used. It is convenient to make the following change of variables: (1) can be re-written in the form dX

The Jacobian of the system (15), at the associated DFE (E 0 ), is given by

Consider the case when R q = 1. Furthermore, suppose that β is chosen as a bifurcation parameter. Solving for β from R q = 1 gives β * = k 1 k 2 k 3 k 4 k 5

The transformed system (15) , with β = β * , has a hyperbolic equilibrium point (i.e., the linearized system has a simple eigenvalue with zero real part, and all other eigenvalues have negative real part), so that the center manifold theory [31, 32] can be used to analyze the dynamics of (15) near β = β * . It can be shown that the right eigenvector of J(E 0 )| β=β * , denoted by w, is given by w = (w 1 , w 2 , . . . , w 7 , w 8 ) T , where,

Similarly, J(E 0 )| β=β * has a left eigenvector, v given by

Consequently, it follows that the associated bifurcation coefficients, a and b (defined in Theorem 4.1 of [32] ), are given, respectively, by

where a 1 = k 3 k 4 k 5 k 6 w 3 + k 3 k 4 k 5 k 6 η q w 4 , a 2 = k 4 k 5 k 6 σ 1 w 3 + k 3 k 5 k 6 η q σ 2 w 4 , a 3 = k 4 k 6 η h α 1 σ 1 w 3 + k 3 k 6 η h α 2 σ 2 w 4 , a 4 = k 4 k 5 φ 1 σ 1 w 3 + k 3 k 5 φ 2 σ 2 w 4 + k 4 φ 3 α 1 σ 1 w 3 + k 3 φ 3 α 2 σ 2 w 4 .

Since the bifurcation coefficient b is always positive, it follows (from Theorem 4.1 of [32] ) that the system (15) will undergo backward bifurcation if the bifurcation coefficient a is positive. This result is summarized below.

Theorem 5. The transformed model (15) , or equivalently (1), exhibits backward bifurcation at R q = 1 whenever the bifurcation coefficient, a, given by (16) , is positive.

It should be noted from (16) that the bifurcation coefficient a, is positive whenever

Thus, the transformed model (15) , exhibits backward bifurcation at R q = 1 whenever r > r c . Furthermore, it should be noted that for the case when quarantined susceptible individuals do not acquire infection during quarantine (i.e., r = 0), the bifurcation coefficient a becomes a = −η q v 3 w 4  µ + σ 2 µ + ψ 1 + η q γ  − v 3 x * 1 k 3 k 4 k 5 k 6 (a 1 + a 2 + a 3 + a 4 ) < 0 since a i > 0 for i = 1, . . . , 4 (since all the model parameters are non-negative). Thus, since a < 0 in this case, it follows from Theorem 4.1 of [32] that the model (1) will not exhibit backward bifurcation if r = 0. In other words, this study shows that the backward bifurcation property of the model (1) arises due to the infection of susceptible individuals in quarantine. This result is consistent with Theorem 3 (where it was shown that the DFE of the model (1) with r = 0 is globally-asymptotically stable).

Modelling strategies for controlling SARS outbreaks

Transmission dynamics and control of severe acute respiratory syndrome

Curtailing transmission of severe acute respiratory syndrome within a community and its hospital

Sensitivity and uncertainty analyses for a SARS model with time-varying inputs and outputs

A cost-based comparison of quarantine strategies for new emerging diseases

Transmission dynamics of etiological agent of SARS in Hong Kong: the impact of public health interventions

Simulating the SARS outbreak in Beijing with limited data

Optimal and sub-optimal quarantine and isolation control in SARS epidemics

Quarantine-based disease control in domesticated animal herds

When should we intervene to control the 2009 influenza A(H1N1) pandemic? Rapid Communications

The basic reproductive number of ebola and the effects of public health measures: the cases of Congo and Uganda

When is quarantine a useful control strategy for emerging infectious diseases? Amer

Epidemiological models with non-exponentially distributed disease stages and application to disease control

Effects of quarantine in six endemic models for infectious diseases

Dynamics of two-strain influenza with isolation and partial cross-immunity

Mathematical study of the impact of quarantine, isolation and vaccination in curtailing an epidemic

Global asymptotic dynamics of a model for quarantine and isolation

The effect of incidence functions on the dynamics of a quarantine/isolation model with time delay

Critical role of nosocomial transmission in the Toronto SARS outbreak

Threshold dynamics of a non-autonomous SEIRS model with quarantine and isolation

Model parameters and outbreak control for SARS

Mathematics in Population Biology

The Theory of the Chemostat

On the definition and computation of the basic reproduction ratio R 0 in models for infectious disease in heterogeneous population

Reproduction numbers and sub-threshold endemic equilibria for compartmental models of disease transmission

The mathematics of infectious diseases

Theoretical assessment of public health impact of imperfect prophylactic HIV-1 vaccines with therapeutic benefits

Backward bifurcations in dengue transmission dynamics

Role of incidence function in vaccine-induced backward bifurcation in some HIV models

Applications of Centre Manifold Theory

Dynamical models of tuberculosis and their applications

Ordinary Differential Equations

Uniform persistence and flows near a closed positively invariant set

Epidemic and demographic interaction in the spread of potentially fatal diseases in growing populations

Global dynamics of a SEIR model with varying total population size