Building Ventilation as an Effective Disease Intervention Strategy in a Dense Indoor Contact Network in an Ideal City

Xiaolei Gao; Jianjian Wei; Hao Lei; Pengcheng Xu; Benjamin J. Cowling; Yuguo Li

doi:10.1371/journal.pone.0162481

Abstract

Emerging diseases may spread rapidly through dense and large urban contact networks, especially they are transmitted by the airborne route, before new vaccines can be made available. Airborne diseases may spread rapidly as people visit different indoor environments and are in frequent contact with others. We constructed a simple indoor contact model for an ideal city with 7 million people and 3 million indoor spaces, and estimated the probability and duration of contact between any two individuals during one day. To do this, we used data from actual censuses, social behavior surveys, building surveys, and ventilation measurements in Hong Kong to define eight population groups and seven indoor location groups. Our indoor contact model was integrated with an existing epidemiological Susceptible, Exposed, Infectious, and Recovered (SEIR) model to estimate disease spread and with the Wells-Riley equation to calculate local infection risks, resulting in an integrated indoor transmission network model. This model was used to estimate the probability of an infected individual infecting others in the city and to study the disease transmission dynamics. We predicted the infection probability of each sub-population under different ventilation systems in each location type in the case of a hypothetical airborne disease outbreak, which is assumed to have the same natural history and infectiousness as smallpox. We compared the effectiveness of controlling ventilation in each location type with other intervention strategies. We conclude that increasing building ventilation rates using methods such as natural ventilation in classrooms, offices, and homes is a relatively effective strategy for airborne diseases in a large city.

Citation: Gao X, Wei J, Lei H, Xu P, Cowling BJ, Li Y (2016) Building Ventilation as an Effective Disease Intervention Strategy in a Dense Indoor Contact Network in an Ideal City. PLoS ONE 11(9): e0162481. https://doi.org/10.1371/journal.pone.0162481

Editor: Jeffrey Shaman, Columbia University, UNITED STATES

Received: March 4, 2016; Accepted: August 23, 2016; Published: September 9, 2016

Copyright: © 2016 Gao et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the paper and its Supporting Information file.

Funding: LY received a General Research Fund grant from The Research Grants Council (17205014) and a grant from National Natural Science Foundation of China (51278440). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: BJC has received research funding from MedImmune Inc. and Sanofi Pasteur, and consults for Crucell NV. BJC is a Section Editor of PLOS ONE. This does not alter the authors’ adherence to all the PLOS ONE policies on sharing data and materials. The authors report no other potential conflicts of interest.

Introduction

Infectious disease epidemics such as the 2003 SARS epidemic, the 2009 H1N1 pandemic, and the 2015 MERS epidemic are a threat to public health. More than 10% of annual deaths worldwide, i.e., 6 million deaths, are caused by viral or bacterial respiratory infections [1]. Due to increasing population density and a more connected world, new or non-vaccinable respiratory infectious diseases [2], especially airborne diseases [3–5], can be widely and rapidly transmitted.

Many interventions for controlling airborne disease transmission have been studied, such as masks, hand hygiene, isolation, quarantine, vaccination, school closure, etc. [6, 7]. Building ventilation is known to be effective for reducing the spread of airborne diseases such as tuberculosis, SARS, smallpox, chicken pox, and influenza in single indoor environments [8–12], but has not been studied at the community level. Ventilation intervention decreases transmission probability by directing the flow of airborne infectious agents away from susceptible persons and/or by removing infectious agents from room air. Because it relies less on individual compliance, ventilation has an advantage over other non-pharmaceutical interventions (e.g., hand washing or mask use). But questions remain. Is ventilation as effective at the community level as other interventions?

Airborne infection, if any, occurs predominantly in indoor environments, where people in modern society spend more than 90% of their time [13]. An individual infected in one building may infect others in the building(s) that he/she subsequently visits. People move from one indoor environment to another, and are in frequent contact with others, creating an indoor contact network in which an airborne disease can spread. To study ventilation as a community-level intervention, the disease spread dynamics must be ideally modeled at the individual person and location level for an entire urban social network [14]. Many widely used social contact network models, such as a random network [15] or a small world network [16, 17], may be too simplistic to describe the complicated dynamic between individuals and indoor environments. Environment parameters such as ventilation cannot be considered in these models, as in these models location is not refined to the individual indoor-space level. More recently, Eubank and colleagues [18, 19] developed a location-based social network model. However, the resources for building a realistic urban network in a large city are generally not available.

One of the significant limitations on the simulation of an indoor-contact network at the community level is computational power. In a large city such as Hong Kong, there are around 7 million people (N_P) and an estimated 3 million indoor locations (N_L, as a building can have many indoor locations). Ideally, we need to know the probability and duration of any two people meeting in any one indoor space in Hong Kong at any time. That is, we need to know “Who meets whom, where, when, and for how long.” Clearly, to model a large city there is need for a fast algorithm and huge computer storage capacity. The computer storage requirement can be formidably large; it is proportional to the cubic of the population (), assuming that the indoor space population is at the same order of population as in Hong Kong. It is thus unrealistic to trace every individual on an hourly scale in a large city, i.e., when and where each person enters and leaves a room.

Consequently, we constructed a probability indoor contact model for an ideal city using real Hong Kong data with a linear memory requirement. Our model used the available data on population, buildings, and indoor spaces; the probability of people visiting each location type; ventilation rates; etc. Then, we developed a probability urban airborne transmission model by integrating the indoor contact network model with both the SEIR model for the spread of community diseases and the Wells-Riley equation for local infection risks. The transmission dynamics of a hypothetical airborne disease outbreak was studied and the effectiveness of different interventions such as isolation, home isolation, and increasing ventilation were evaluated and compared. In this feasibility study, we assumed that this hypothetical airborne disease has the same natural history and infectiousness as smallpox with no pre-existing immunity in the population.

Methods

Building an indoor contact network for an ideal large city

Indoor contact was defined as two or more individuals visiting the same indoor location at the same time. In other studies, contacts are often defined as close, conversational, or physical contact [20, 21]. Our indoor contact network model starts with the so-called individual/location approach, which has three stages.

Stage 1.

Generating both people and indoor space populations. A synthetic population was first generated that includes every individual in Hong Kong. At the same time, a synthetic indoor space (location) was generated that includes every indoor space in Hong Kong. 6,857,100 people were divided into eight population groups (home stayers, office workers, classroom attendees, food service workers, shop workers, drivers, public space workers, and others), and 2,923,035 locations into seven location groups (home, office, classroom, restaurant, shop, transport vehicles, and other public locations). Each individual was assigned a unique person ID and an occupation (population group), and each location a unique location ID, a function (location group), a size, an occupant density, and a ventilation rate. The inclusion of ventilation rate as a characteristic of locations allows us to study for the first time the effect of ventilation on the spread of disease.

Stage 2.

Defining the location visiting schedule and behavior of each individual. This is a crucial stage for tackling the formidably cubic computer storage requirement. To construct a simple probability indoor contact model with a linear memory requirement, we proposed an ideal city with the following three simplifications.

First, each population group has a fixed daily schedule. Individuals in each population group have a fixed and identical location visiting schedule. All of the individuals in the entire population spend the same duration at any location in a location group. The location visiting schedule for every individual in a population group is defined, e.g., for office workers it is “home → transportation → office → lunch → office→ transportation → dinner → home” (see Fig 1). We denote g = 1,2,…, G as the sequence of location groups of individuals’ location visiting behavior during a simulation day. As we assume that all of the visitors enter and leave a location at the same time, we can simply define the meeting duration in a location group g, as τ_g. We refer to this assumption as the “parallel session” approach, as in typical academic conferences.

Download:

Fig 1. An example of hourly indoor activity sequences for a family of three individuals with different occupations.

They meet other people in different types of locations (represented with different colors and shapes) during a typical workday. There are 6,857,100 people and 2,923,035 locations in Hong Kong.

https://doi.org/10.1371/journal.pone.0162481.g001

Second, we assumed that any individual in a population group has the same number of favorite locations in each location group. For each location type, individuals choose a fixed number of favorite locations such as favorite shops, restaurants, etc., although offices and homes are unique to each individual. Individuals only visit their favorite locations. We use l = 1,2,…L_g to represent an individual’s collection of possible choices in each location group.

Third, the individual location visiting probability for locations in each location group was assumed the same for all individuals from the same population group. Hence, the probability of visiting a favorite location can be determined for each individual in each location group.

The first simplification deals with “when and how long,” and the second with “where.” These two simplifications allow us to use an array η_j,l (l = 1,2,…L_g) to represent the probability of individual j visiting location l. As there are only g = 1,2,…, G location groups, the total storage requirement is reduced to a linear requirement N_P × L_g × G. This simplification allows us to handle a population of 7 million visiting 3 million locations on our desktop computers.

The third simplification implies that the location visiting behavior of individuals is determined solely by their occupations. McPherson et al. [22] demonstrated that social contact behavior is strongly related to the occupation and education level of individuals. Edmunds et al. [20] and Read et al. [23] also suggested that occupation is a primary influence on face-to-face contacts, especially during weekdays. The contact patterns and contact time are correlated with contact locations [23].

Hence, a pre-defined 24-hour activity sequence can be easily assigned to each synthetic individual, according to their profession and social duties. The time step for the activity sequence was chosen to be 30 minutes. Finer time steps can also be used, but are not necessary as the incubation period is relatively long.

Note that we do not associate each individual with other demographic attributes such as income level or age. In our simple model, each population group has the same location visiting schedule, and each location group has the same occupant density. Visiting probabilities for individuals in a population group to a type location were generated randomly assuming the population distribution follows a normal distribution (the mean is the location-visiting proportion and the standard deviation is10% of the mean value, see S1 File for detailed description). We used the available statistical data on population, buildings, and the probability of people visiting different indoor environments. For example, the probability that students will have lunch in a restaurant, or that a teacher will take a bus to go to work, were estimated according to the available data drawn from social behavior studies and household expenditure surveys.

Stage 3.

Constructing the indoor contact network (people-location probability network, ). For our approach, a probability network is preferred, given the uncertainty of individual daily routes. In this people-location probability network, all of the individuals are connected through their favorite locations, see Fig 2. Both individuals and locations are represented by nodes. An edge between an individual node and a location node denotes a connection between the individual and location (i.e., a location visiting). The weight of the edge represents the probability that the individual visits the location in a day. In implementation, location nodes were selected to connect each individual with their favorite locations, and the edges were weighted according to the probability that the individual will visit the location group. A recursion algorithm was used to ensure that the occupants of all of the locations satisfy their limits.

Download:

Fig 2. Revised Eubank’s bipartite graph of a small indoor contact network.

(a) A bipartite graph of a probability network with two types of vertex and four people (represented with a p or a smiling face), and four indoor spaces (represented by l or an open square). P1 has a favorite location. Each person p can visit each location l with a given probability, so there is an edge in the graph between each p and each l. People vertices are labeled with occupation and favorite locations. Location vertices are labeled with type, floor area, and ventilation rate. Edges are labeled with arrival and departure times, and visiting probability. (b) A bipartite graph of a fixed network, which is otherwise the same as (a), but the edges are labeled with arrival and departure times. (c, d) The two disconnected graphs Gp and Gl induced by connecting vertices that were separated by exactly two edges in Gpl. (e, f). The static projections resulting from ignoring time labels in Gp and Gl.

https://doi.org/10.1371/journal.pone.0162481.g002

Detailed descriptions of the data collection process and the indoor contact network construction methodologies are given in S1 File. In reality, each individual in the network might connect with many other locations with different edge weights. A real network is also dynamic with some edges changing periodically (e.g., weekly or seasonal changes) or permanently. However, it is in practice impossible to describe mathematically these complicated features of the network. Here, we neglected weekly or seasonal changes when calculating .

Integrating the indoor contact network model with the Wells-Riley equation

The airborne transmission risk is estimated using the Wells-Riley equation [24]. In an enclosed indoor environment, infectious airborne droplet nuclei are assumed to be evenly dispersed. Let the probability of susceptible individual i being infected by the only infector j in indoor environment l be p_i,j,l, which can be expressed by the following equation [24]: (1) where Q is the average number of infectious quanta generated by a infector (quanta/h); p is the pulmonary ventilation rate, and here we use 0.38 m³/h [25]; τ_l is the duration of exposure to infection (h); and q_l is the ventilation rate of the location (m³/h). The infection bears a Poisson relation to the dose of inhaled infective agents (droplets or droplet nuclei), according to Wells [24].

One possible solution for integration is to generate the social contact network from the people-location probability network, and simulate the disease’s dynamic with the Markov chain Monte Carlo (MCMC) method. However, this method is too time-consuming for a large population. Hence, we simulated the epidemic as a Markov chain with each simulation day as a time step. We defined the infection states of individuals as probabilities. The calculated local infection risk of an individual in a contaminated location using the Wells-Riley equation becomes a joint probability.

Let be the probability of individual j being an infector on day t, and η_j,l represent the probability of individual j visiting location l. The visiting probability of an individual to a location is assumed to be independent of other individuals’ choices. In the meantime, the visiting probabilities within the set of favorite locations are given, as the indoor contact network is built. Therefore, the location visiting probability of individual i to location l is independent of the other probabilities, including the probability of infection. This means that and η_j,l are independent, and hence the probability of individual i being infected by individual j in location l in the g^th location group in day t, , can be calculated as (2)

The visiting time of any individual to a given location is identical, according to the following equation: (3)

The probabilities of individual i being infected from different infectors in location l can be assumed to be independent. This is not difficult to understand, as a quanta is defined as a dose of infective agents that, when inhaled, can cause an infection. Assuming the concentration of each quanta is constant and well-mixed, the event of breathing in a quanta became a Poisson process [23]. Hence, the probability of individual i being infected in location l can be calculated as (4)

Here, N_l is the number of possible visitors to location l. Assuming there are favorite locations in the g^th location group, we mark them as . It is assumed that the infection risks from the favorite locations that an individual may choose in a visiting step are non-correlated. As individual i can only be in one location at a time and has a given probability of being infected in the selected location, the infection probabilities of individual i for different locations in the same location group are mutually exclusive. Thus, we can calculate the infection probability for individual i in the location group , as (5)

is the probability of individual i being susceptible to infection on simulation day t − 1. As the daily schedule of an individual’s location visiting behavior is given, for each step of the individual’s movement, the infection stage of the individual can be assumed to correlate only with the previous step. Hence, we use to demonstrate the susceptibility of an individual after visiting group g as follows: (6)

Then, (7)

According to Eqs (3), (4), (5) and (7), (8)

In the sketch of the natural history of the airborne disease (the same as smallpox), shown in Fig 3, each state of an individual’s infection is represented by the probability of this event. A susceptible individual i, after exposure, will pass through a latent period () before infection, followed by a pre-symptomatic infection period (), eventually reaching an infectious state. Patients with severe symptoms will be hospitalized and isolated (), and people with symptoms will comply with home isolation (). The death rate is not considered in this case.

Download:

Fig 3. The natural history of the airborne disease infection.

We assume that in the natural history of airborne diseases a susceptible individual will pass through a latent period (, 11.6 days) after being exposed, followed by a symptomatic state (, 2.49 days), and will eventually reach an infectious state during which a rash appears (, 18.3 days). Both the fever onset and rash states are considered infectious [26]. Infected individuals will either recover or die (). We assume that in the baseline outbreak scenario, 50% of the patients with a rash will be hospitalized and isolated (), 95% of the rest will be isolated at home (); during the latter period the individual will stay home and not visit any other locations.

https://doi.org/10.1371/journal.pone.0162481.g003

A simple case of three individuals visiting four locations in two location groups is shown in Fig 4.

Download:

Fig 4. Calculation of transmission probability based on different assumptions.

A. A simple example of the location visiting behavior of A, B, and C. Individual A is susceptible to infection and individuals B and C both have the probability of being an infector. Each of the three individuals has two location visiting steps in their daily schedule. In each simulation day they can choose to visit either location J or K from group 1 in the first step, and location H or M from group 2 in the second step. The location visiting probability of individual A to location J is defined as η_A,J. B. The diagrammatic flow chart of the method for calculating infection risks from individuals B and C to individual A (Solid arrow) and the correlations of variables (dashed arrow) given by Assumptions 1 to 4.

https://doi.org/10.1371/journal.pone.0162481.g004

The epidemic was modeled discretely with a stated time step of one day. The state of each individual in the infection-recovery process in the current time step depends on their state in the previous time step. Hence, the transmission process becomes a Markov chain. The Markov chain transmission probability for a single individual can be represented by the following equations: (9) where , , , , , , and are the probabilities that individual i on day t is, respectively, susceptible, exposed (latent), pre-symptomatic, symptomatic infectious, hospitalized, home-isolated, and recovered or removed from the social network at day t. 1 / θ, 1 / ϕ, and 1 / γ are the latent, pre-symptomatic infectious, and total infectious period. ϖ and ϑ are the hospitalization rate and home isolation rate. According to Eq (8) we have (10) where is the location visiting probability after home isolation. If l is the home of j, , otherwise .

Eq (10) is referred to as the integrated indoor contact network and disease transmission model. It uses the probabilities of individuals i and j visiting location l, that is η_i,l, and η_j,l, from the probability indoor contact network model.

Input parameters

To test the indoor contact network model, statistical data describing population groups and location groups were collected from the Hong Kong Census and Statistics Department and other related studies. The ventilation rates of the locations were estimated by reviewing relevant field measurements conducted in Hong Kong (see S1 File). The baseline ventilation rates were set as: 0.7 air change per hour (ACH) in homes, 2 ACH in classrooms, 1 ACH in offices, 1 ACH in restaurants, 1 ACH in shops, 4 ACH in transportations, and 1.4 ACH in other public locations. Unlike the agent-based approach of Eubank et al. [18], the input data for our simple network model are minimal. We are not interested in the exact choices of activity, travel mode, departure and estimation (and or time) choice, car following, etc.

One major difficulty in using the Wells-Riley equation for our purpose is obtaining a meaningful quanta generation rate. It may be estimated from reported outbreaks, but for some diseases the estimated values of different outbreaks are significantly different [10, 12, 25, 27–29]. Such significant variations may be due to differences in the sub-strains of the viruses, respiratory activities (number and size distribution of droplets or droplet nuclei generated from the index patient) [30, 31], the state of infection of the source patients [32], possible existence of super-spreaders [33], environment conditions such as temperature or humidity that relate to virus or bacterial survival rate [34], and significant uncertainties in the estimated ventilation condition [28]. Therefore, the estimated quanta generation rate from reported outbreaks may not be suitable for analyzing population transmission dynamics.

In intervention studies, the basic reproductive number, R_o, is commonly used as an indicator of the transmission ability of the infection [35]. The transmission (diffusion) of diseases through a social network is generally estimated using the effective reproductive number. The correlation of quanta generation rate and R_o can be calculated for a single room outbreak by integrating the Wells-Riley equation and the compartment model [28]. This inspired us to propose an alternative approach for obtaining the quanta generation rate. We assumed that infectors have a constant quanta generation rate in the disease outbreak. We modeled such an outbreak in Hong Kong by assuming that the disease is transmitted by the airborne route only. The simulations were carried out with a range of quanta generation rates of 1, 2 or 3 quanta/h. The basic reproductive number was calculated for each simulation, with a respectively value of 2.22, 4.03 and 5.83 A correlation was then established between the basic reproductive numbers and the quanta generation rates Q.

The range of the estimated basic reproductive numbers agrees with the commonly adopted values of smallpox transmission, which are in the range of 3.5 to 10 [26, 35, 36]. Hence, we adopted three cases representing low, medium, and high transmission rate scenarios.

Results

Structure of the indoor contact network

Our model produced a bipartite graph with two types of vertices, people and location. In our case, all of the locations are indoor spaces. People are connected by location (indoor spaces), and indoor spaces are connected by people. To understand the basic feature of the network structure, we generated a fixed network G_I,I, based on a people-location probability network , using random numbers to generate fixed edges between connected individual nodes and location nodes in network . The degree distribution of G_I,L is plotted in Fig 5B. The location degree distribution of graph G_I,L is better fitted to the power-law distribution than the population degrees distribution of graph G_I,I, which agrees with the findings of Eubank et al. [18]. The location graph is scale-free, meaning that a few locations (hubs) have many “visitors,” but the majority of locations have relatively few “visitors.” Eubank et al. argued that removing locations with more than 100 degrees (visitors) may improve the results. Such a network structure can be revealing in terms of infection control [18], e.g., locations of over 100 degrees need better ventilation.

Download:

Fig 5. Degree distributions of people-location graph.

A. Degree distribution of individuals with different numbers of daily visiting locations. B. Degree distribution of locations with different numbers of visitors (per day) on a log-log scale. C. Degree distributions of people-people graph on a log-log scale. There are 6,857,100 people and 2,923,035 locations in Hong Kong.

https://doi.org/10.1371/journal.pone.0162481.g005

The people-people graph G_I,I can then be generated by connecting individuals who are linked by common locations based on G_I,L with degree distribution; the result is shown in Fig 5C. The people contact graph is like a small-world, highly connected (not scale-free) network. Eubank et al. [18] suggested that vaccinating long-distance travelers is crucial for infection control. Our study cannot validate this conclusion, as we do not model travel distance. However, as Hong Kong has effective public transport and a small land area (less than 1000 km²), most of the locations visited by an individual are within one-hour travel. The people-people graph fits reasonably well with the predictions in [18]. The noise of two degree distributions might have been introduced by the compartmentalization of population and locations.

A hypothetical airborne disease outbreak under the baseline scenario

In this study, we assume the complete susceptibility of the whole population. In the baseline outbreak scenario, the air change rates of different types of locations are set to be consistent with data from investigations of ventilation rates (see S1 File). The baseline outbreak case is simulated based on the indoor contact model without any interventions. We present the epi-curve of the outbreak, infection of individuals in different population groups, and infection risk in different types of locations in Fig 6.

Download:

Fig 6. Outbreak dynamics of a hypothetical airborne disease outbreak in Hong Kong under the baseline scenario.

A. The expected percentages of susceptible, exposed, infectious, and recovered individuals. B. The expected percentages of infectors in each population group. C. Average infection risks exposed in locations in different location groups.

https://doi.org/10.1371/journal.pone.0162481.g006

The changes in the daily numbers of individuals in each infection state are plotted in Fig 6A. The peak of infection in the community occurs 154 days after the first infector is introduced into the community. Most of the individuals in the community are infected and the outbreak lasts for about 500 days. According to Fig 6B, the rate of infection among students rises first, but the number of infectors is highest in offices workers and lowest in home stayers.

Among the different locations, the highest infection risks are in homes, followed by classrooms and offices. The infection risk for individuals working in restaurants, shops, public places, and transportation systems are not listed here, as the data in Fig 6C were calculated using the exposure time of the general population rather than those who work in these locations.

Effects of increasing ventilation rates under the medium transmission scenario

The effects under different control methods are plotted in Fig 7. The attack rate, δ, time of peak infection, T_p (day), and percentage of peak infection, λ_p, under Policies A to P are listed in Table 1. The control effects of increasing ventilation rates in one type of location (Policies A to H) under the medium transmission condition are shown in Fig 7A. Increasing the ventilation rate to 5 ACH (Policy A) or 10 ACH (Policy G) in all homes is an effective method for reducing the epidemic, and it reduced the attack rate substantially (around 50%). However, the differences between the two policies are not significant. Increasing the ventilation rate of all offices from 1 ACH to 5 ACH delays the arrival of peak infection, but it does not reduce much of the overall attack rate.

Download:

Fig 7. Control effects of increasing ventilation rates, isolation, and vaccination under a quanta generation rate of 2 quanta/h.

A. Daily percentages of infectors under the baseline condition and control policies of increasing the ventilation rate in one type of location (Policy A to H). B. Baseline condition and increasing ventilation rate in all locations or a few types of locations (Policy I to P). The curves for Policy K and P are flat. C. Different ventilation, isolation, and vaccination control methods. The curves for Policy VAC3, VEN2, VEN3, and COM2 are all nearly flat (not shown).

https://doi.org/10.1371/journal.pone.0162481.g007

Download:

Table 1. Different ventilation control methods and their effects under a quanta generation rate of 2 quanta/h.

https://doi.org/10.1371/journal.pone.0162481.t001

Increasing ventilation rates in all locations, demonstrated in Fig 7B, is more effective in delaying the peak infection than ventilation control in any one type of location, as it reduces the rapid transmission of the disease in large hubs at the very beginning of the outbreak. The policies of doubling ventilation rates (Policy I), applying 8 liters per second per person in all locations (Policy L), and applying 12 liters per second per person in all locations (Policy M) are all more effective than increasing the household ventilation rate to 5 ACH (about 85 liter per second per person; see S1 File). Increasing the ventilation rate in homes and offices largely delays and reduces the peak infection and has a control effect similar to that of a policy of applying 3 ACH in all locations (Policy J). Applying 5 ACH in all locations (Policy K) will diminish the outbreak.

Comparison of effects of increasing ventilation rates with isolation and mass-vaccination

When the quanta generation rate is 2 quanta/h, the epi-curves of the daily percentages of infectors under different control policies (Table 2) are as plotted in Fig 7C. The control effects of all of the methods under different quanta generation rates are given in Table 2.

Download:

Table 2. Effects of isolation, vaccination, and ventilation on attack rate (δ), day of peak infection (T_p), and percentage of infectors on day T_p (λ_p) under different quanta generation rates.

https://doi.org/10.1371/journal.pone.0162481.t002

The ability to control the epidemic with isolation is limited in the medium and high transmission scenarios, even though 90% of the infectors with rashes are isolated (13% of the population at the peak state, 2 quant/h) and 50% of the infectors with pre-symptomatic transmission are isolated in their households (5% of the population at the peak state, 2 quant/h). Also, the effect of mass-vaccination depends on the ability to distribute vaccinations. Vaccinating 100,000 susceptible individuals per day can diminish the infection in all three cases, whereas giving 10,000 vaccinations per day can only diminish the outbreak under the low transmission scenario; it leads to a reduction of less than 20% in the attack rates in the medium and high transmission scenarios.

Table 2 shows that higher ventilation rates can significantly reduce the disease transmission; in particular, it can delay and reduce the peak infection. Increasing the household ventilation rate (Policy VEN1) can mostly eliminate an outbreak under the low transmission rate, and is more effective than isolation (Policy ISO) or vaccination (Policy VAC1) under the medium and high transmission rate conditions. Applying isolation and increasing the household ventilation rates together (Policy COM1) largely restricts the outbreak in all three cases, and the number of isolated individuals is also 50 times smaller (less than 0.4% of population) than when isolation is applied alone in the medium transmission scenario.

Discussion

Ventilation is an effective intervention, particularly when the transmission rate is high

In this study, when the disease is completely airborne transmitted, increasing ventilation rates has demonstrated similar or even better control effects than isolation and vaccination. Furthermore, ventilation rates in homes, classrooms, and offices can reach up to 5 ACH simply by opening windows [37] or through the proper design of ventilation systems [38, 39], which is easier, more convenient, and more economical than isolating 12% of the population or mass-vaccination. Additionally, in the low transmission rate condition, both isolation and vaccination have acceptable control effects.

However, when the transmission rate is high, infection spreads rapidly before vaccinations can be distributed and the existence of the pre-symptomatic transmission limits the effect of isolation. In this case, combining increased ventilation rates with isolation has the best control effect and decreases the load of isolation requirements.

Increasing ventilation in critical locations is important

Our comparison of increasing ventilation rate policies suggests that it is essential to apply targeted control strategies to the high infection risk locations. Homes are the most important locations for controlling infection throughout the community, due to longer exposure periods, household isolation of infectors, and most importantly the role of homes in connecting the infection transmission network. The ventilation rate can be readily increased to 5 ACH taking use of natural ventilation [40, 41].

In contrast, classrooms and offices play a role in connecting small hubs (homes) with medium and large hubs in the network. Infection in classrooms develops in the early stages, as shown in Fig 6, perhaps because of the indoor social contact networks and the ventilation rate settings in classrooms. In all classrooms, if an infector is introduced to the classroom, the rest of the students will develop infection rapidly (due to the long exposure time and low ventilation rate) and the transmission conditions in different classrooms are similar (similar ventilation rates per person). Natural ventilation was suggested by Lee and Chang [42] to be efficient in increasing the ventilation rates of classrooms. As offices are extremely large hubs in the network, reducing infection risks in offices will delay the peak infection. Hence, increasing ventilation rates in homes is essential to protect the public from infection from airborne diseases, whereas applying targeted control methods such as increasing ventilation rates to large location hubs may delay the rise of the epi-curve and provide time for the public to prepare for the outbreak.

Advantages and limitations of our approaches

The advantage of our indoor contact network model is that only a small amount of data is required to build a model of the indoor contact network of a large city and to simultaneously capture the heterogeneous population contact pattern and uncertainty of the network structure. The indoor contact network in Hong Kong has a larger number of vertices with a higher degree of contacts than the people-people contact network in the city of Portland [18]. The highly connected graph developed in this study is potentially related to the common habits of eating out and using public transportation among inhabitants of Hong Kong.

Our limited input data (i.e. number locations, location visiting frequencies) and our assumptions (i.e. number of favorite locations), may influence the accuracy of our network topology. Further studies can address this problem by using more detailed data or fitting the network in a small-scale, real contact network investigation. This study also did not consider the geographical contact distribution, which may have a significant influence on the dynamics of infectious disease transmission [43, 44]; the local outbreaks cannot be predicted in the current model. Future studies are necessary to evaluate the geographical aspects of the social contact network and the role of possible transmission hubs that are not included in this study such as hospitals.

Contacts between friends are not simulated in our contact network model due to the difficulty of identifying the number of friends an individual may have within and outside the individual’s subgroup and the lack of social investigation data. The gathering activities are also hard to define in terms of frequency and location. Similarly, the individual diversity is not considered. Age is not integrated into the indoor contact network model due to the lack of data on age distribution within each population subgroup. In addition, the individual vaccine histories and the presence of super infectors have important effects on predicting the outbreak. A more complex model can be built when more detailed data is available.

Although only a limited number of intervention methods were compared, as increasing ventilation rate is the major concern of in this study, the effects of other intervention methods such as school closure, contact tracing, ultraviolet germicidal irradiation, and wearing masks could be evaluated in the same social contact network.

One significant limitation of this study is that we assumed that the disease (similar with smallpox) is transmitted only via the airborne route. In reality, respiratory infections may have multiple possible transmission routes [45, 46]. The transmission dynamics and the control effects of different methods might shift with the dominant transmission route or routes. However, the aim of the study is to demonstrate the effectiveness of ventilation in combination with other intervention methods in controlling a possible airborne disease outbreak on a social contact network, as well as our ability to simulate different scenarios on the theoretical basis. In another study by the authors [47], the role of ventilation in controlling the partially airborne influenza was investigated; it was found that increasing the ventilation rate significantly impacts on the transmission dynamics, even when the airborne route only contributes 20% to the total infection risk. Besides, as the methodology of constructing an indoor contact network provided by this study contains information about both individual-based social contact behavior and local environment transmission risks, a similar methodology could be adopted to model the multiple-routes-transmitted respiratory infections in large cities.

When applying the Wells-Riley model to estimate the local transmission rate, the selection of quanta generation rate is still a crucial problem. Most studies have applied a quanta generation rate calculated from historical outbreaks [10, 25, 28, 48, 49]. This method might lead to an over-estimation of the infectivity of the diseases due to the super-spreader effect [33]. Furthermore, the number of outbreak cases with sufficient information to estimate the quanta generation rate is very small, and a significant bias will be introduced by this small number of statistical data. However, other infectious diseases outbreak parameters, such as the basic reproductive number, are well documented for the past infectious disease outbreaks. In this study, we provide a method for investigating the correlations between basic reproductive number and quanta generation rate. The quanta generation in the previous outbreak can be estimated from the basic reproductive number by rebuilding the social contact network during the outbreak period.

Supporting Information

S1 File. Supplementary information.

The following information is described in detail, e.g. construction of an indoor contact network, ventilation rate profile of indoor environments in Hong Kong, and determining the quanta generation rate.

https://doi.org/10.1371/journal.pone.0162481.s001

(PDF)

Acknowledgments

This study was funded by a General Research Fund grant from RGC (17205014) and an NSFC grant (51278440).

Author Contributions

Conceptualization: YL XG.
Data curation: XG.
Formal analysis: XG JW YL.
Funding acquisition: YL.
Investigation: XG YL JW.
Methodology: XG PX YL.
Project administration: YL.
Resources: YL.
Software: XG.
Supervision: YL BJC.
Validation: XG.
Visualization: XG JW YL.
Writing – original draft: XG YL.
Writing – review & editing: XG JW HL BJC YL.

References

1. WHO. The Global Burden of Disease: 2004 Available: http://www.who.int/healthinfo/global_burden_disease/2004_report_update/en/index.html.
2. Feldmann H. Truly emerging—A new disease caused by a novel virus. New England Journal Medicine. 2011; 364: 1561–1563
- View Article
- Google Scholar
3. Yu IT, Li Y, Wong TW, Tam W, Chan AT, Lee JHW, et al. Evidence of airborne transmission of the severe acute respiratory syndrome virus. New England Journal Medicine. 2004; 350: 1731–1739.
- View Article
- Google Scholar
4. Tellier R. Review of aerosol transmission of influenza A virus. Emerging Infectious Diseases. 2006; 12: 1657–1662. pmid:17283614
- View Article
- PubMed/NCBI
- Google Scholar
5. Tellier R. Aerosol transmission of influenza A virus: A review of new studies. Journal of the Royal Society Interface. 2009; vol. 6 (Suppl 6): S783–S790.
- View Article
- Google Scholar
6. Jefferson T, Del Mar CB, Dooley L, Ferroni E, Al-Ansary LA, Bawazeer GA, et al. Physical interventions to interrupt or reduce the spread of respiratory viruses: Systematic review. BMJ. 2009; 339: b3675. pmid:19773323
- View Article
- PubMed/NCBI
- Google Scholar
7. Cowling BJ, Zhou Y, Ip DKM, Leung GM, Aiello AE. Face masks to prevent transmission of influenza virus: A systematic review. Epidemiology and Infection. 2010; 138: 449–456. pmid:20092668
- View Article
- PubMed/NCBI
- Google Scholar
8. Li Y, Leung GM, Tang JW, Yang X, Chao CYH, Lin JZ, et al. Role of ventilation in airborne transmission of infectious agents in the built environment—a multidisciplinary systematic review. Indoor Air. 2007; 17: 2–18. pmid:17257148
- View Article
- PubMed/NCBI
- Google Scholar
9. Nardell EA, Keegan J, Cheney SA, Etkind SC. Airborne infection—Theoretical limits of protection achievable by building ventilation. American Review of Respiratory Disease. 1991; 144: 302–306. pmid:1907115
- View Article
- PubMed/NCBI
- Google Scholar
10. Beggs CB, Noakes CJ, Sleigh PA, Fletcher LA, Siddiqi K. The transmission of tuberculosis in confined spaces: An analytical review of alternative epidemiological models. The International Journal of Tuberculosis and Lung Disease. 2003; 7: 1015–1026. pmid:14598959
- View Article
- PubMed/NCBI
- Google Scholar
11. Noakes CJ, Sleigh PA. Mathematical models for assessing the role of airflow on the risk of airborne infection in hospital wards. Journal of the Royal Society Interface. 2009; 6(Supp 6): S791–S800.
- View Article
- Google Scholar
12. Gao X, Li Y, Leung GM. Ventilation control of indoor transmission of airborne diseases in an urban community. Indoor and Built Environment. 2009; 18: 205–218.
- View Article
- Google Scholar
13. Leech JA, Nelson WC, Burnett RT, Aaron S, Raizenne ME. It’s about time: a comparison of Canadian and American timeactivity patterns. J Expo Anal Environ Epidemiol. 2002; 12: 427–432. pmid:12415491
- View Article
- PubMed/NCBI
- Google Scholar
14. Keeling MJ, Eames KTD. Networks and epidemic models. Journal of the Royal Society Interface. 2005; 2: 295–307.
- View Article
- Google Scholar
15. Newman MEJ. Spread of epidemic disease on networks. Physical Review E. 2002; 66: No. 016128
- View Article
- Google Scholar
16. Lin G, Jia X, Ouyang Q. Predict SARS infection with the small world network model, Beijing Da Xue Xue Bao. 2003; 35: 66–69. pmid:12914222
- View Article
- PubMed/NCBI
- Google Scholar
17. Moore C, Newman MEJ. Epidemics and percolation in small-world networks. Physical Review E. 2000; 61: 5678–5682.
- View Article
- Google Scholar
18. Eubank S, Guclu H, Kumar VS, Marathe MV, Srinivasan A, Toroczkai Z, et al. Modelling disease outbreaks in realistic urban social networks. Nature. 2004; 429: 180–184. pmid:15141212
- View Article
- PubMed/NCBI
- Google Scholar
19. Eubank S, Kumar VS, Marathe MV, Srinivasan A, Wang N. Structure of social contact networks and their impact on epidemics. AMS-DIMACS Special Volume on Epidemiology. 2006; 70: 181–213.
- View Article
- Google Scholar
20. Edmunds WJ, O’Callaghan CJ, Nokes DJ. Who mixes with whom? A method to determine the contact patterns of adults that may lead to the spread of airborne infections. Proceedings of the Royal Society B: Biological Sciences. 1997; 264: 949–957. pmid:9263464
- View Article
- PubMed/NCBI
- Google Scholar
21. Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, et al. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008; 5: e74. pmid:18366252
- View Article
- PubMed/NCBI
- Google Scholar
22. McPherson M, Smith-Lovin L, Cook JM. Birds of a feather: Homophily in social networks. Annual Review of Sociology. 2001; 27: 415–444.
- View Article
- Google Scholar
23. Read JM, Eames KTD, Edmunds WJ. Dynamic social networks and the implications for the spread of infectious disease. Journal of the Royal Society Interface. 2008; 5: 1001–1007.
- View Article
- Google Scholar
24. Wells WF. Airborne contagion and air hygiene: An ecological study of droplet infections. Harvard Univ. Press, Cambridge, MA. 1955. pp423.
25. Chen SC, Chang CF, Liao CM. Predictive models of control strategies involved in containing indoor airborne infections. Indoor Air. 2006; 16: 469–481. pmid:17100668
- View Article
- PubMed/NCBI
- Google Scholar
26. Riley S, Ferguson NM. Smallpox Transmission and Control: Spatial Dynamics in Great Britain. Proc Natl Acad Sci USA. 2006; 103: 12637–12642 pmid:16894173
- View Article
- PubMed/NCBI
- Google Scholar
27. Qian H, Li Y, Nielsen PV, Huang XH. Spatial distribution of infection risk of SARS transmission in a hospital ward. Building and Environment. 2009; 44: 1651–1658.
- View Article
- Google Scholar
28. Rudnick SN, Milton DK. Risk of indoor airborne infection transmission estimated from carbon dioxide concentration. Indoor Air. 2003; 13: 237–245. pmid:12950586
- View Article
- PubMed/NCBI
- Google Scholar
29. Gammaitoni L, Nucci MC. Using a mathematical model to evaluate the efficacy of TB control measures. Emerging Infectious Diseases. 1997; 3: 335–342. pmid:9284378
- View Article
- PubMed/NCBI
- Google Scholar
30. Xie X. Evaporation and Movement of Respiratory Droplets in Indoor Environments. PhD thesis, The University of Hong Kong. 2008. pp. 196.
31. Morawska L. Droplet fate in indoor environments, or can we prevent the spread of infection? Indoor Air. 2006; 16: 335–347. pmid:16948710
- View Article
- PubMed/NCBI
- Google Scholar
32. Ferguson NM, Cummings DAT, Cauchemez S, Fraser C, Riley S, Meeyai A, et al. Strategies for containing an emerging influenza pandemic in Southeast Asia. Nature. 2005; 437: 209–214. pmid:16079797
- View Article
- PubMed/NCBI
- Google Scholar
33. Li Y, Yu ITS, Xu P, Lee JHW, Wong TW, Ooi PL, et al. Predicting super spreading events during the 2003 severe acute respiratory syndrome epidemics in Hong Kong and Singapore. American Journal of Epidemiology. 2004; 160: 719–728. pmid:15466494
- View Article
- PubMed/NCBI
- Google Scholar
34. Lowen AC, Mubareka S, Steel J, Palese P. Influenza virus transmission is dependent on relative humidity and temperature. PLoS Pathogens. 2007; 3: e151.
- View Article
- Google Scholar
35. Gani R, Leach S. Transmission potential of smallpox in contemporary populations. Nature. 2001; 414: 748–751 pmid:11742399
- View Article
- PubMed/NCBI
- Google Scholar
36. Ferguson NM, Keeling MJ, Edmunds WJ, Gani R, Grenfell BT, Anderson RM, et al. Planning for smallpox outbreaks. Nature. 2003; 425: 681–685. pmid:14562094
- View Article
- PubMed/NCBI
- Google Scholar
37. Escombe AR, Oeser CC, Gilman RH, Navincopa M, Ticona E, Pan W, et al. Natural ventilation for the prevention of airborne contagion. PLoS Med. 2007; 4: e68. pmid:17326709
- View Article
- PubMed/NCBI
- Google Scholar
38. Allard F, editor. Natural ventilation in buildings: A design handbook. London: James & James. 1998.
39. Atkinson J, Chartier Y, Pessoa-Silva C, Jensen P, Li Y, Seto W, editors. Natural Ventilation for Infection Control in Health-Care Settings- WHO Guidelines: World Health organization. 2009.
40. Chao CYH. Comparison between indoor and outdoor air contaminant levels in residential buildings from passive sampler study. Building and Environment. 2001; 36, 999–1007.
- View Article
- Google Scholar
41. Chao CYH, Tung TC. An empirical model for outdoor contaminant transmission into residential buildings and experimental verification, Atmospheric Environment.2001; 35, 1585–1596.
- View Article
- Google Scholar
42. Lee SC, Chang M. Indoor and outdoor air quality investigation at schools in Hong Kong, Chemosphere. 2000; 41, 109–113. pmid:10819186
- View Article
- PubMed/NCBI
- Google Scholar
43. Anderson RM, May RM. Infectious Diseases of Humans: Dynamics and Control. Oxford, New York, Oxford University Press. 1991. pp. 688.
44. Riley S. Large-scale spatial-transmission models of infectious disease. Science. 2007; 316: 1298–1301. pmid:17540894
- View Article
- PubMed/NCBI
- Google Scholar
45. Roy CJ, Milton DK. Airborne transmission of communicable infection—the elusive pathway. New England Journal Medicine. 2004; 350: 1710–1712.
- View Article
- Google Scholar
46. Tang JW, Li Y, Eames I, Chan PK, Ridgway GL. Factors involved in the aerosol transmission of infection and control of ventilation in healthcare premises. Journal of Hospital Infection. 2006; 64: 100–114. pmid:16916564
- View Article
- PubMed/NCBI
- Google Scholar
47. Gao X, Wei J, Cowling BJ, Li Y. Potential impact of a ventilation intervention for influenza in the context of a dense indoor contact network in Hong Kong. Science of the Total Environment. 2016; 569–570: 373–381 pmid:27351145
- View Article
- PubMed/NCBI
- Google Scholar
48. Liao CM, Chen SC, Chang CF. Modelling respiratory infection control measure effects. Epidemiology and Infection. 2008; 136: 299–308. pmid:17475088
- View Article
- PubMed/NCBI
- Google Scholar
49. Riley EC, Murphy G, Riley RL. Airborne spread of measles in a suburban elementary school. American Journal of Epidemiology. 1978; 107: 421–432. pmid:665658
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. WHO. The Global Burden of Disease: 2004 Available: http://www.who.int/healthinfo/global_burden_disease/2004_report_update/en/index.html.

[ref2] 2. Feldmann H. Truly emerging—A new disease caused by a novel virus. New England Journal Medicine. 2011; 364: 1561–1563
View Article
Google Scholar

[3] View Article

[4] Google Scholar

[ref3] 3. Yu IT, Li Y, Wong TW, Tam W, Chan AT, Lee JHW, et al. Evidence of airborne transmission of the severe acute respiratory syndrome virus. New England Journal Medicine. 2004; 350: 1731–1739.
View Article
Google Scholar

[6] View Article

[7] Google Scholar

[ref4] 4. Tellier R. Review of aerosol transmission of influenza A virus. Emerging Infectious Diseases. 2006; 12: 1657–1662. pmid:17283614
View Article
PubMed/NCBI
Google Scholar

[9] View Article

[10] PubMed/NCBI

[11] Google Scholar

[ref5] 5. Tellier R. Aerosol transmission of influenza A virus: A review of new studies. Journal of the Royal Society Interface. 2009; vol. 6 (Suppl 6): S783–S790.
View Article
Google Scholar

[13] View Article

[14] Google Scholar

[ref6] 6. Jefferson T, Del Mar CB, Dooley L, Ferroni E, Al-Ansary LA, Bawazeer GA, et al. Physical interventions to interrupt or reduce the spread of respiratory viruses: Systematic review. BMJ. 2009; 339: b3675. pmid:19773323
View Article
PubMed/NCBI
Google Scholar

[16] View Article

[17] PubMed/NCBI

[18] Google Scholar

[ref7] 7. Cowling BJ, Zhou Y, Ip DKM, Leung GM, Aiello AE. Face masks to prevent transmission of influenza virus: A systematic review. Epidemiology and Infection. 2010; 138: 449–456. pmid:20092668
View Article
PubMed/NCBI
Google Scholar

[20] View Article

[21] PubMed/NCBI

[22] Google Scholar

[ref8] 8. Li Y, Leung GM, Tang JW, Yang X, Chao CYH, Lin JZ, et al. Role of ventilation in airborne transmission of infectious agents in the built environment—a multidisciplinary systematic review. Indoor Air. 2007; 17: 2–18. pmid:17257148
View Article
PubMed/NCBI
Google Scholar

[24] View Article

[25] PubMed/NCBI

[26] Google Scholar

[ref9] 9. Nardell EA, Keegan J, Cheney SA, Etkind SC. Airborne infection—Theoretical limits of protection achievable by building ventilation. American Review of Respiratory Disease. 1991; 144: 302–306. pmid:1907115
View Article
PubMed/NCBI
Google Scholar

[28] View Article

[29] PubMed/NCBI

[30] Google Scholar

[ref10] 10. Beggs CB, Noakes CJ, Sleigh PA, Fletcher LA, Siddiqi K. The transmission of tuberculosis in confined spaces: An analytical review of alternative epidemiological models. The International Journal of Tuberculosis and Lung Disease. 2003; 7: 1015–1026. pmid:14598959
View Article
PubMed/NCBI
Google Scholar

[32] View Article

[33] PubMed/NCBI

[34] Google Scholar

[ref11] 11. Noakes CJ, Sleigh PA. Mathematical models for assessing the role of airflow on the risk of airborne infection in hospital wards. Journal of the Royal Society Interface. 2009; 6(Supp 6): S791–S800.
View Article
Google Scholar

[36] View Article

[37] Google Scholar

[ref12] 12. Gao X, Li Y, Leung GM. Ventilation control of indoor transmission of airborne diseases in an urban community. Indoor and Built Environment. 2009; 18: 205–218.
View Article
Google Scholar

[39] View Article

[40] Google Scholar

[ref13] 13. Leech JA, Nelson WC, Burnett RT, Aaron S, Raizenne ME. It’s about time: a comparison of Canadian and American timeactivity patterns. J Expo Anal Environ Epidemiol. 2002; 12: 427–432. pmid:12415491
View Article
PubMed/NCBI
Google Scholar

[42] View Article

[43] PubMed/NCBI

[44] Google Scholar

[ref14] 14. Keeling MJ, Eames KTD. Networks and epidemic models. Journal of the Royal Society Interface. 2005; 2: 295–307.
View Article
Google Scholar

[46] View Article

[47] Google Scholar

[ref15] 15. Newman MEJ. Spread of epidemic disease on networks. Physical Review E. 2002; 66: No. 016128
View Article
Google Scholar

[49] View Article

[50] Google Scholar

[ref16] 16. Lin G, Jia X, Ouyang Q. Predict SARS infection with the small world network model, Beijing Da Xue Xue Bao. 2003; 35: 66–69. pmid:12914222
View Article
PubMed/NCBI
Google Scholar

[52] View Article

[53] PubMed/NCBI

[54] Google Scholar

[ref17] 17. Moore C, Newman MEJ. Epidemics and percolation in small-world networks. Physical Review E. 2000; 61: 5678–5682.
View Article
Google Scholar

[56] View Article

[57] Google Scholar

[ref18] 18. Eubank S, Guclu H, Kumar VS, Marathe MV, Srinivasan A, Toroczkai Z, et al. Modelling disease outbreaks in realistic urban social networks. Nature. 2004; 429: 180–184. pmid:15141212
View Article
PubMed/NCBI
Google Scholar

[59] View Article

[60] PubMed/NCBI

[61] Google Scholar

[ref19] 19. Eubank S, Kumar VS, Marathe MV, Srinivasan A, Wang N. Structure of social contact networks and their impact on epidemics. AMS-DIMACS Special Volume on Epidemiology. 2006; 70: 181–213.
View Article
Google Scholar

[63] View Article

[64] Google Scholar

[ref20] 20. Edmunds WJ, O’Callaghan CJ, Nokes DJ. Who mixes with whom? A method to determine the contact patterns of adults that may lead to the spread of airborne infections. Proceedings of the Royal Society B: Biological Sciences. 1997; 264: 949–957. pmid:9263464
View Article
PubMed/NCBI
Google Scholar

[66] View Article

[67] PubMed/NCBI

[68] Google Scholar

[ref21] 21. Mossong J, Hens N, Jit M, Beutels P, Auranen K, Mikolajczyk R, et al. Social contacts and mixing patterns relevant to the spread of infectious diseases. PLoS Med. 2008; 5: e74. pmid:18366252
View Article
PubMed/NCBI
Google Scholar

[70] View Article

[71] PubMed/NCBI

[72] Google Scholar

[ref22] 22. McPherson M, Smith-Lovin L, Cook JM. Birds of a feather: Homophily in social networks. Annual Review of Sociology. 2001; 27: 415–444.
View Article
Google Scholar

[74] View Article

[75] Google Scholar

[ref23] 23. Read JM, Eames KTD, Edmunds WJ. Dynamic social networks and the implications for the spread of infectious disease. Journal of the Royal Society Interface. 2008; 5: 1001–1007.
View Article
Google Scholar

[77] View Article

[78] Google Scholar

[ref24] 24. Wells WF. Airborne contagion and air hygiene: An ecological study of droplet infections. Harvard Univ. Press, Cambridge, MA. 1955. pp423.

[ref25] 25. Chen SC, Chang CF, Liao CM. Predictive models of control strategies involved in containing indoor airborne infections. Indoor Air. 2006; 16: 469–481. pmid:17100668
View Article
PubMed/NCBI
Google Scholar

[81] View Article

[82] PubMed/NCBI

[83] Google Scholar

[ref26] 26. Riley S, Ferguson NM. Smallpox Transmission and Control: Spatial Dynamics in Great Britain. Proc Natl Acad Sci USA. 2006; 103: 12637–12642 pmid:16894173
View Article
PubMed/NCBI
Google Scholar

[85] View Article

[86] PubMed/NCBI

[87] Google Scholar

[ref27] 27. Qian H, Li Y, Nielsen PV, Huang XH. Spatial distribution of infection risk of SARS transmission in a hospital ward. Building and Environment. 2009; 44: 1651–1658.
View Article
Google Scholar

[89] View Article

[90] Google Scholar

[ref28] 28. Rudnick SN, Milton DK. Risk of indoor airborne infection transmission estimated from carbon dioxide concentration. Indoor Air. 2003; 13: 237–245. pmid:12950586
View Article
PubMed/NCBI
Google Scholar

[92] View Article

[93] PubMed/NCBI

[94] Google Scholar

[ref29] 29. Gammaitoni L, Nucci MC. Using a mathematical model to evaluate the efficacy of TB control measures. Emerging Infectious Diseases. 1997; 3: 335–342. pmid:9284378
View Article
PubMed/NCBI
Google Scholar

[96] View Article

[97] PubMed/NCBI

[98] Google Scholar

[ref30] 30. Xie X. Evaporation and Movement of Respiratory Droplets in Indoor Environments. PhD thesis, The University of Hong Kong. 2008. pp. 196.

[ref31] 31. Morawska L. Droplet fate in indoor environments, or can we prevent the spread of infection? Indoor Air. 2006; 16: 335–347. pmid:16948710
View Article
PubMed/NCBI
Google Scholar

[101] View Article

[102] PubMed/NCBI

[103] Google Scholar

[ref32] 32. Ferguson NM, Cummings DAT, Cauchemez S, Fraser C, Riley S, Meeyai A, et al. Strategies for containing an emerging influenza pandemic in Southeast Asia. Nature. 2005; 437: 209–214. pmid:16079797
View Article
PubMed/NCBI
Google Scholar

[105] View Article

[106] PubMed/NCBI

[107] Google Scholar

[ref33] 33. Li Y, Yu ITS, Xu P, Lee JHW, Wong TW, Ooi PL, et al. Predicting super spreading events during the 2003 severe acute respiratory syndrome epidemics in Hong Kong and Singapore. American Journal of Epidemiology. 2004; 160: 719–728. pmid:15466494
View Article
PubMed/NCBI
Google Scholar

[109] View Article

[110] PubMed/NCBI

[111] Google Scholar

[ref34] 34. Lowen AC, Mubareka S, Steel J, Palese P. Influenza virus transmission is dependent on relative humidity and temperature. PLoS Pathogens. 2007; 3: e151.
View Article
Google Scholar

[113] View Article

[114] Google Scholar

[ref35] 35. Gani R, Leach S. Transmission potential of smallpox in contemporary populations. Nature. 2001; 414: 748–751 pmid:11742399
View Article
PubMed/NCBI
Google Scholar

[116] View Article

[117] PubMed/NCBI

[118] Google Scholar

[ref36] 36. Ferguson NM, Keeling MJ, Edmunds WJ, Gani R, Grenfell BT, Anderson RM, et al. Planning for smallpox outbreaks. Nature. 2003; 425: 681–685. pmid:14562094
View Article
PubMed/NCBI
Google Scholar

[120] View Article

[121] PubMed/NCBI

[122] Google Scholar

[ref37] 37. Escombe AR, Oeser CC, Gilman RH, Navincopa M, Ticona E, Pan W, et al. Natural ventilation for the prevention of airborne contagion. PLoS Med. 2007; 4: e68. pmid:17326709
View Article
PubMed/NCBI
Google Scholar

[124] View Article

[125] PubMed/NCBI

[126] Google Scholar

[ref38] 38. Allard F, editor. Natural ventilation in buildings: A design handbook. London: James & James. 1998.

[ref39] 39. Atkinson J, Chartier Y, Pessoa-Silva C, Jensen P, Li Y, Seto W, editors. Natural Ventilation for Infection Control in Health-Care Settings- WHO Guidelines: World Health organization. 2009.

[ref40] 40. Chao CYH. Comparison between indoor and outdoor air contaminant levels in residential buildings from passive sampler study. Building and Environment. 2001; 36, 999–1007.
View Article
Google Scholar

[130] View Article

[131] Google Scholar

[ref41] 41. Chao CYH, Tung TC. An empirical model for outdoor contaminant transmission into residential buildings and experimental verification, Atmospheric Environment.2001; 35, 1585–1596.
View Article
Google Scholar

[133] View Article

[134] Google Scholar

[ref42] 42. Lee SC, Chang M. Indoor and outdoor air quality investigation at schools in Hong Kong, Chemosphere. 2000; 41, 109–113. pmid:10819186
View Article
PubMed/NCBI
Google Scholar

[136] View Article

[137] PubMed/NCBI

[138] Google Scholar

[ref43] 43. Anderson RM, May RM. Infectious Diseases of Humans: Dynamics and Control. Oxford, New York, Oxford University Press. 1991. pp. 688.

[ref44] 44. Riley S. Large-scale spatial-transmission models of infectious disease. Science. 2007; 316: 1298–1301. pmid:17540894
View Article
PubMed/NCBI
Google Scholar

[141] View Article

[142] PubMed/NCBI

[143] Google Scholar

[ref45] 45. Roy CJ, Milton DK. Airborne transmission of communicable infection—the elusive pathway. New England Journal Medicine. 2004; 350: 1710–1712.
View Article
Google Scholar

[145] View Article

[146] Google Scholar

[ref46] 46. Tang JW, Li Y, Eames I, Chan PK, Ridgway GL. Factors involved in the aerosol transmission of infection and control of ventilation in healthcare premises. Journal of Hospital Infection. 2006; 64: 100–114. pmid:16916564
View Article
PubMed/NCBI
Google Scholar

[148] View Article

[149] PubMed/NCBI

[150] Google Scholar

[ref47] 47. Gao X, Wei J, Cowling BJ, Li Y. Potential impact of a ventilation intervention for influenza in the context of a dense indoor contact network in Hong Kong. Science of the Total Environment. 2016; 569–570: 373–381 pmid:27351145
View Article
PubMed/NCBI
Google Scholar

[152] View Article

[153] PubMed/NCBI

[154] Google Scholar

[ref48] 48. Liao CM, Chen SC, Chang CF. Modelling respiratory infection control measure effects. Epidemiology and Infection. 2008; 136: 299–308. pmid:17475088
View Article
PubMed/NCBI
Google Scholar

[156] View Article

[157] PubMed/NCBI

[158] Google Scholar

[ref49] 49. Riley EC, Murphy G, Riley RL. Airborne spread of measles in a suburban elementary school. American Journal of Epidemiology. 1978; 107: 421–432. pmid:665658
View Article
PubMed/NCBI
Google Scholar

[160] View Article

[161] PubMed/NCBI

[162] Google Scholar

Figures

Abstract

Introduction

Methods

Building an indoor contact network for an ideal large city

Stage 1.

Stage 2.

Stage 3.

Integrating the indoor contact network model with the Wells-Riley equation

Input parameters

Results

Structure of the indoor contact network

A hypothetical airborne disease outbreak under the baseline scenario

Effects of increasing ventilation rates under the medium transmission scenario

Comparison of effects of increasing ventilation rates with isolation and mass-vaccination

Discussion

Ventilation is an effective intervention, particularly when the transmission rate is high

Increasing ventilation in critical locations is important

Advantages and limitations of our approaches

Supporting Information

S1 File. Supplementary information.

Acknowledgments

Author Contributions

References