## Abstract

In the event of an influenza pandemic, the most probable way in which the virus would be introduced to an isolated geographical area is by an infected traveller. We use a mathematical model, structured on the location at which infection occurs and based on published parameters for influenza, to describe an epidemic in a community of one million people. The model is then modified to reflect a variety of control strategies based on social distancing measures, targeted antiviral treatment and antiviral prophylaxis and home quarantine, and the effectiveness of the strategies is compared. The results suggest that the only single strategy that would be successful in preventing an epidemic (with _{0}=2.0) is targeted antiviral treatment and prophylaxis, and that closing schools combined with either closing work places or home quarantine would only prevent such an epidemic if these strategies were combined with a modest level of antiviral coverage.

## 1. Introduction

When a pandemic of an emerging infectious disease threatens, there are considerable advantages to isolated geographical regions. New Zealand is one example of such isolation, and this study addresses the threat to the country from pandemic influenza. Travel times from any other major land mass mean that the importation of infection through surface contamination is unlikely, and those arriving by sea would have had time to express symptoms en route. Hence, if pandemic influenza were to reach the country, the most probable scenario would be that an epidemic had already started elsewhere, and the virus would be introduced to New Zealand by an infected traveller. A response plan for dealing with such an event has been formulated (Jennings 2005; Coker & Mournier-Jack 2006). The major point of entry to the country is Auckland, which is also by far the biggest city and the most probable region for an epidemic to be initiated.

We presupposed a single unintentional introduction of a virus by an airline passenger, not detected by screening at the border. We used a structured Kermack–McKendrick integral equation model to describe the spread of a novel influenza-like virus in an isolated population. We used model parameters consistent with other published models of influenza, and investigated the ability of a range of control scenarios to contain the outbreak. It is anticipated that if and when an influenza pandemic begins, observations from other countries and regions will enable us to revise our parameters, and all our calculations will be repeated. In the meantime, it is believed that the model and parameters are based on reasonable estimates of what could happen, and these results are assisting the planning of a response.

The model is summarized in figure 1, and described in appendix A. Similar models have previously been used to evaluate control policies for severe acute respiratory syndrome (SARS; Roberts 2004, 2006) and smallpox (Aldis & Roberts 2005).

A feature of the model is that susceptibles are treated as a homogeneous group, but the incidence of infection is structured according to the location where the infection takes place. The locations at which infection occurs are within the household, at school, within the workplace and within the community. Those that were infected at location ℓ infect others at location *k* at a rate proportional to the component *W*_{kℓ} of the mixing matrix ** W**, which is scaled so that ‖

**‖=1. It was proved by Roberts (submitted) that the final size equation for this model is the same as that for the unstructured model (Kermack & McKendrick 1927; Diekmann & Heesterbeek 2000 and references therein),(1.1)where**

*W**P*is the proportion of the population that becomes infected over the course of the epidemic and

_{0}is the basic reproduction number, defined as the expected number of secondary cases that would arise from a typical primary case introduced into a fully susceptible population (Diekmann

*et al*. 1990; Diekmann & Heesterbeek 2000). It was also proved by Roberts (submitted) that the proportion of the population infected at each location is proportional to the components of the eigenvector of

**that corresponds to the unit eigenvalue. These results are used here to assess the performance of control measures designed to reduce the impact of an influenza epidemic, in terms of the reduction in the total number of infections that result.**

*W*## 2. Method

We solved the model numerically for values of the basic reproduction number _{0}=2.0, 1.1 and 3.0, to determine the potential size and time-course of the epidemic in a region of one million people. The assumed parameter values are summarized in table 1. These are mostly taken from other published modelling exercises, and would be revised in the event of a pandemic as data become available. See the papers cited and their citations for details of the assumptions that led to these estimates. The exception is the assumed 95% uptake of antivirals. It was observed that 81% of passengers on a cruise ship took antiviral prophylaxis when offered during an outbreak of seasonal influenza (Miller *et al*. 2000). As pandemic influenza is a more serious threat to health, a higher uptake was assumed.

We investigated the potential effects of control measures based on social distancing by determining how their implementation would modify the mixing matrix ** W** to form a new matrix

*W*_{c}. The basic reproduction number under control is then

_{c}=‖

*W*_{c}‖

_{0}. Hence, we determined the ‘threshold’ control effort, which reduces

_{c}to 1. For example, reducing transmission in schools would multiply the second row of the matrix

**by a factor**

*W**c*to obtain

*W*_{c}. The threshold value determined was that value of

*c*which led to ‖

*W*_{c}‖=1/

_{0}and hence

_{c}=1.

Finally, we postulated 12 control scenarios that combined social distancing, targeted antiviral treatment and antiviral prophylaxis (TATP) and home quarantine. For each, we calculated the value of _{c} and the proportion of the population infected in an epidemic under the application of selected control scenarios. The scenarios are summarized as follows.

No intervention.

Close schools. Transmission at school ceases.

Close workplaces (70% compliance). Transmission at work is reduced to 30% of the no intervention level.

Close schools and workplaces. A combination of ii and iii above.

TATP. Antiviral treatment of cases prevents 45.6% of all transmission, and when combined with antiviral prophylaxis of home contacts, this leads to an overall reduction of 65.2% of within-household transmission. This is calculated as follows: 80% of transmission is due to symptomatic cases (see table 1) and 95% uptake of antivirals combined with 60% efficacy prevents 57% of this. Hence, out of the transmission that would occur without intervention, 20% is from asymptomatic cases, 45.6% is prevented by antiviral treatment and 34.4% continues from symptomatic cases. Antiviral prophylaxis of home contacts prevents a further 57% of the 34.4% of transmission due to symptomatic cases within the same household.

Home quarantine (70% compliance). Health authorities recommend that members of a household with a symptomatic case remain at home for 6 days and 70% comply. Hence, 56% (equal to 70% of 80%) of all transmission from those infected within their own household is prevented.

Home quarantine (50% compliance). Hence, 40% of transmission from household contacts is prevented.

Close schools and home quarantine. A combination of 2 and 4 above.

Close schools and TATP. A combination of 2 and 5 above.

TATP and home quarantine. A combination of 5 and 6 above. Hence, 66.9% (equal to 70% of the 80% of transmission that would result from symptomatic cases who comply with quarantine, plus 46.6% of the other 30% (of the 80%) prevented by antivirals) of transmission from those infected at home, 65.2% of within-household transmission and 45.6% of all other transmission (see v above) are prevented.

Close schools, TATP and quarantine. A combination of 2, 5 and 6 above.

Close schools and workplaces, TATP and quarantine. A combination of 4–6 above.

For all scenarios, it was assumed that if transmission were reduced at one location, there would be no compensatory increase in transmission elsewhere.

## 3. Results

### 3.1 The uncontrolled epidemic

The solution of the final size equation (1.1) is shown in figure 2. When _{0}=2.0, *P*=0.7968, which implies that 79.7% of the population will be infected in an epidemic. The solutions for _{0}=1.1 and _{0}=3.0 imply that 17.6% or 94.0%, respectively, would be infected. The results from some recent simulation exercises for pandemic influenza are also presented in figure 2 for comparison.

Numerically solving the structured epidemic model without any control intervention for _{0}=2.0 resulted in an epidemic lasting approximately 100 days (figure 3). The other numerical results obtained (not shown) resulted in incidence curves of a similar shape, except that for _{0}=3.0, the epidemic was completed in approximately 80 days and for _{0}=1.1, the epidemic lasted over 600 days. It was confirmed that the final size equation (1.1) estimates the proportion of the population infected in an epidemic, and that the distribution of cases among infection locations (see table 1) is determined by the eigenvector of ** W** that corresponds to the unit eigenvalue.

### 3.2 The controlled epidemic

The threshold values for the level of control required to prevent epidemics are presented in table 2. These are the values at which _{c}=1.0. These show, for example, that if _{0} were 1.1 then transmission in schools would need to be reduced to 63% of its previous (no control) level to prevent epidemics, but this could also be achieved by reducing transmission to 45% in workplaces or to 80% in both schools and workplaces. Note that these values are the percentage of the pre-control value, hence the figure of 80% translates to a 20% reduction in transmission.

The values of _{c}, the basic reproduction number under control, calculated for 12 specific control interventions and _{0}=2.0 are presented in figure 4. If _{c}>1.0 for an intervention, then the final size equation (1.1) was used to compute *P*, the proportion of the population that would be infected if those control measures were in place throughout the epidemic. These are also shown in figure 4. The equivalent results for _{0}=1.1 and _{0}=3.0 are shown in figure 5.

## 4. Discussion

We have constructed and analysed a model for the transmission of a virus, using only the parameter values in table 1. These are based on those in the literature, but they are still only educated guesses. Until a form of pandemic influenza transmitted between humans emerges, we can only base the parameters on previous influenza epidemics in different social conditions and due to different strains of the virus. However, the model has been constructed so that the parameters may be changed and the results revised should a pandemic eventuate. As the results are sensitive to the parameter values, they should only be taken to indicate what could happen should these parameters prove to be approximately correct.

The final size of the epidemic, as computed, shows that should the basic reproduction number (_{0}) be equal to 2.0 we could expect nearly 80% of the population to become infected. The population here is a homogeneously mixing community of one million people, motivated by the City of Auckland (population approx. 1.2 million). The epidemic would peak at approximately 65 days, with peak incidence of approximately 43 000 new cases per day. Increasing _{0} to 3.0 results in a faster epidemic infecting more than 90% of the population. For different size populations the proportion infected is determined only by _{0}, but a larger population size results in an epidemic that takes longer to complete.

The proportion of the population that is infected in an epidemic, as determined by equation (1.1), is plotted in figure 2, where it is compared with the results obtained by other authors from simulation models for outbreaks of pandemic influenza. The results suggest that the final size equation (1.1) overestimates the proportion infected in an epidemic as determined by the more complex models, except in the case of Wu *et al*. (2006). This could be due, in part, to several reasons, apart from the obvious one that the structure of all these models is more complicated and the equation does not necessarily apply. For example, with a stochastic, spatially structured, individual-based model, local exhaustion of susceptible contacts will occur. In addition, the analytic final size equation is for the proportion infected over an infinite time-scale, whereas numerical simulations often cease when the incidence becomes small and the epidemic ‘is largely over’ (Ferguson *et al*. 2005). For complicated models, it is not always clear how the value of _{0} relates to the parameters in the model, and the supplementary information to Germann *et al*. (2006) provides a discussion of this point. Hence, there are reasons for uncertainty about the points plotted in figure 2, in both the vertical and horizontal directions.

We have modelled control measures in two ways. We represented social distancing by reducing the amount of transmission occurring at one or more locations. If this were applied uniformly across all locations and at a level sufficient to prevent a major outbreak, then transmission would need to be multiplied by 1/_{0} to achieve _{c}=1. The results in table 2 also show the effects of social distancing applied to specific locations, for example if _{0} were equal to 2.0, it would be necessary to eliminate 85% of transmission that takes place in schools and workplaces to prevent a major epidemic. This assumes that there is no compensatory increase in transmission, for example within the household. One can envisage that closing schools could result in increased household transmission, as the children would spend more time at home, and possibly the parents too. However, there could be a considerable educational component associated with taking this action, leading to improved hygiene and decreased transmission within the household. Hence, we have taken a neutral position, assuming no change in the model. If _{0} were less than 2.0, then the prospects for control would be more optimistic.

Some specific control strategies have been investigated, with the results presented in figures 4 and 5. All of the strategies successful in preventing an epidemic, with _{0}=2.0 (_{c}<1), involved TATP, with the assumption of 95% uptake and 60% efficacy. However, the strategy of closing schools (100% compliance) combined with either closing workplaces or home quarantine (70% compliance; no. 4 and 8 in figure 4) resulted in values of _{c} just over 1. These strategies, combined with a modest amount of antiviral coverage would be sufficient to reduce _{c} below 1. Under the low basic reproduction number assumption (_{0}=1.1), all strategies considered apart from closing schools were sufficient to contain epidemics. Under the high basic reproduction number assumption (_{0}=3.0), the only strategies sufficient to contain epidemics involved the use of antivirals, school closure and home quarantine. These results are derived assuming that there is no delay from the start of the epidemic until putting the strategy in place. As the scenario is that a pandemic has started elsewhere and is subsequently imported, it is reasonable to assume that the authorities have a response plan in place and the delay is minimal. Hence, the proportion infected will be approximately as calculated. This caveat does not affect the basic reproduction number under control (_{c}) for each strategy.

Experiments with different contact structures in the model have shown the overall results to be quite robust, depending only on _{0} for total numbers infected, and on the latent and infectious periods, and therefore the generation time, for the incidence curve for the epidemic. The number of infections that take place within the household, within school or the workplace, or within the wider community is determined by the assumed percentages at the beginning of an epidemic (table 1). As these weights also determine to a great extent the success or otherwise of control measures, they should be a prime target for scrutiny should a pandemic eventuate.

## Acknowledgments

This project was commissioned and funded by the New Zealand Ministry of Health, and the material was contained in an unpublished report to the Ministry. The authors are very grateful for the information and comments provided by staff in the Public Health Directorate, Ministry of Health (including Alyson Baker, Martin Bonné, John Boyd, Andreya Brown, Andrea Forde, Alison Roberts, David Sinclair and Paul White). They are also grateful for the input of the external peer reviewers, Angela McLean and Patricia Priest. The views expressed are those of the authors, and do not necessarily reflect the views of the New Zealand Ministry of Health.

## Footnotes

- Received October 3, 2006.
- Accepted October 22, 2006.

- © 2006 The Royal Society