Abstract
Some infectious diseases, including COVID-19, can be transmitted via aerosols that are emitted by an infectious person and inhaled by susceptible individuals. Although physical distancing effectively reduces short-range airborne transmission, many infections have occurred when sharing room air despite maintaining distancing. We propose two simple parameters as indicators of infection risk for this situation. They combine the key factors that control airborne disease transmission indoors: virus-containing aerosol generation rate, breathing flow rate, masking and its quality, ventilation and air cleaning rates, number of occupants, and duration of exposure. COVID-19 outbreaks show a clear trend in relation to these parameters that is consistent with an airborne infection model, supporting the importance of airborne transmission for these outbreaks. The observed trends of outbreak size vs. risk parameters allow us to recommend values of the parameters to minimize COVID-19 indoor infection risk. All of the pre-pandemic spaces are in a regime where they are highly sensitive to mitigation efforts. Measles outbreaks occur at much lower risk parameter values than COVID-19, while tuberculosis outbreaks are observed at much higher risk parameter values. Since both diseases are accepted as airborne, the fact that COVID-19 is less contagious than measles does not rule out airborne transmission. It is important that future outbreak reports include ventilation information, to allow expanding our knowledge of the circumstances conducive to airborne transmission of different diseases.
1. Introduction
Some respiratory infections can be transmitted through the airborne pathway, in which aerosol particles (< 100 μm) are shed by infected individuals and inhaled by others, causing disease in susceptible individuals (1–4). It is widely accepted that measles, tuberculosis, and chickenpox are transmitted in this way (5, 6), and acceptance is growing that this is a major and potentially the dominant transmission mode of COVID-19 (7–13). There is substantial evidence that smallpox (14), influenza (3), SARS (15, 16), MERS (6), and rhinovirus (17), are or were also transmitted via aerosols.
There are three airborne transmission scenarios of interest in which infectious and susceptible people: (a) are in close proximity to each other (< 1-2 m), so-called “short-range airborne transmission,” (18), which is effectively mitigated by physical distancing; (b) share air in the same room, “shared space airborne”; and (c) “longer-distance airborne transmission” in buildings when individuals are not sharing a room (or are in very large rooms), or even between buildings as in the Amoy Gardens SARS outbreak (16). Often (b) and (c) are lumped together under “long-range transmission,” but it is useful to separate them given the substantial differences in the risk of these situations and actions needed to abate the risk of transmission.
Airborne diseases vary widely in transmissibility, but all of them are most easily transmitted at short range due to the higher concentration of pathogen-containing aerosols close to the infected person. For SARS-CoV-2, a pathogen of moderate infectivity, many instances have been reported implicating transmission in shared indoor spaces. Indeed, multiple outbreaks of COVID-19 have been reported in crowded spaces that were relatively poorly ventilated and that were shared by many people for periods of half an hour or longer. Examples include choir rehearsals (11), religious services (19), buses (20), workshop rooms (20), restaurants (4, 21), and gyms (22), among others. There are only a few documented cases of longer-distance transmission of SARS-CoV-2, in buildings (23–26). However, cases of longer-distance transmission are harder to detect as they require contact tracing teams to have sufficient data to connect cases together and rule out infection elsewhere. Historically, it was only possible to prove longer-distance transmission in the complete absence of community transmission (e.g. ref (14)).
Being able to quickly assess the risk of infection for a wide variety of indoor environments is of the utmost importance given the impact of the continuing pandemic (and the risk of future pandemics) on so many aspects of life in almost every country of the world. We urgently need to improve safety of the air we breathe across a range of environments including child-care facilities, kindergartens, schools, colleges, shops, offices, homes, eldercare facilities, factories, public and private transportation, restaurants, gyms, libraries, cinemas, concert halls, places of worship and mass outdoor events and across different climates and socio-economic conditions. There is very little evidence on the actual ventilation rates and the effectiveness of ventilation systems at reducing risks from viral exposure and other indoor pollutants within the majority of buildings. However, data from COVID-19 outbreaks consistently shows that multiple buildings worldwide have very low ventilation rates despite the requirements set in national building standards. A host of policy questions – from how to safely re-open schools to how to prevent transmission in high-risk occupational settings – require accurate quantification of the multiple interacting variables that influence airborne infection risk.
Qualitative guidance to reduce the risk of airborne transmission has been published (27–29). Different mathematical models have been proposed to help manage risk of airborne transmission (30, 31), and several models have been adapted to COVID-19 (32–35). It is important to define quantitative infection risk criteria for different spaces and types of events to more effectively manage the pandemic (36). Such criteria could then be used by authorities and policy makers to assist in deciding which activities are permitted under what conditions, so as to limit infection risk across a society. To our knowledge, no such quantitative criteria have been proposed.
Here, we use a box model to estimate the viral aerosol concentration indoors, and combine it with the Wells-Riley infection model (30). The combined model is used to derive two quantitative risk parameters that allow comparing the relative risk of transmission in different situations when sharing room air. We explore the trends in infections observed in outbreaks of COVID-19 and other diseases as a function of these parameters. Finally, we use the parameters to quantify a graphical display of the relative risk of different situations and mitigation options.
2. Model development
2.1. Box model of infection
The box model considers a single enclosed space, in which virus-containing aerosols are assumed to be rapidly uniformly mixed compared with the time spent by the occupants in the space. This is approximately applicable in many situations, but there are some exceptions such as outbreaks where clear directional flow occurred in a room (19, 21). The mathematical notation used in the paper is summarized in Table SI-1 in the Supp. Info. The mass balance equation is first written in terms of c, the concentration of infectious quanta in the air in the enclosed space (units of quanta m-3). Compared with a model written in terms of aerosol or viral particle concentrations, c has the advantage of implicitly including effects such as the deposition efficiency of the aerosol particles in the lungs of a susceptible person, as well as the efficiency with which such deposited particles may cause infection, the multiplicity of infection (MOI), etc. The balance of quanta in the space can be written as: where Ep is the emission rate of quanta into the indoor air from an infected person present in the space (quanta / h); fe is the penetration efficiency of virus-carrying particles through masks or face coverings for exhalation (which takes into account the impact of whether the infector wears a face covering); V is the volume of the space; λ0 is the first-order rate of removal of quanta by ventilation with outdoor air (h-1); λcle is the removal of quanta by air cleaning devices (e.g. recirculated air with filtering, germicidal UV, portable air cleaners, etc.); λdec is the infectivity decay rate of the virus; λdep is the deposition rate of airborne virus-containing particles onto surfaces. Ep is a critical parameter that depends strongly on the disease, and it can be estimated with a forward model based on aerosol emission rates and pathogen concentration in saliva and respiratory fluid (33, 37) or by fitting a model such as the one used here to real transmission events (11, 33).
This equation can be solved analytically or numerically for specific situations. Given the enormous number of possible situations, and given the prevalence of outbreaks resulting from longer events where air is shared for a substantial period of time, we consider a sufficiently long event so that a steady state condition (dc/dt ∼ 0) is a reasonable approximation. Writing λ = λ0 + λcle + λdec + λdep for simplicity leads to a steady-state infectious quanta concentration of: Under the assumption of no infectious quanta at the beginning of the event, a multiplicative factor, rss, can be applied for events too short to approximately reach steady state (see section SI-1 in the Supp. Info. for detail) to correct the deviation of the quanta concentration averaged over the event (cavg) from that at steady state: Since the goal is to analyze outbreaks, we assume that only a single infectious person is present in the space, which is thought to be applicable to the outbreaks analyzed below. This allows calculation of the probability of infection, conditional to one infectious person being present. The model can also be formulated to calculate the absolute probability of infection, if we assume that the probability of an infectious person being present reflects the prevalence of a disease at a given location and time (e.g. refs (32, 38)).
The dose expressed in infectious quanta (n) inhaled by each of the susceptible persons present in the space (Nsus) is then: where fi is the penetration efficiency of virus-carrying particles through masks or face coverings for inhalation (which takes into account the effect of the fraction of occupants wearing face coverings); B is the breathing volumetric flow rate of susceptible persons; D is the duration of exposure, assumed to be the same for all the susceptible persons. Substituting: The number of expected secondary infections increases monotonically with increasing n. For an individual susceptible person,by definition of an infectious quantum, the probability of infection is (30): For low values of n, the use of the Taylor expansion for an exponential allows approximating P as: And the total number of secondary infections expected is then: Thus, the number of secondary infections increases linearly with n at lower values, and non-linearly at higher values. We retain the simplified form to define and calculate the risk parameters, but use Eq. 6 for fitting the outbreak results in Figure 1b below.
2.2. Risk parameters for airborne infection
We define the hazard parameter, H, for airborne infection in a shared space. The purpose of H is to capture the dependency of Nsi on the parameters that define an event in a given space, in particular those parameters that can be controlled to reduce risk. To better capture the controllable actions, Ep and B were split into factors that can and cannot be controlled. Ep can be expressed as the product Ep0 x rE where Ep0 and rE are, respectively, the quanta shedding rate of an infectious person resting and only orally breathing (no vocalization); and the shedding rate enhancement factor relative to Ep0 for an activity with a certain degree of vocalization and physical intensity (see Table SI-2a for detail). B can be expressed as B0 x rB, where B0 and rB are the volumetric breathing rate of a sedentary susceptible person in the age group of 41-<51 years (numerically also the average for all age groups) and the relative breathing rate enhancement factor (vs. B0) for an activity with a certain physical intensity and for a certain age group (see Table SI-2b for detail). Ep0 is uncertain, likely highly variable across the population, and variable over time during the period of infectiousness (33, 37, 39–41). It may also increase due to new virus variants such as the COVID-19 B1.1.7 variant, that are more contagious, assuming that the increased contagiousness is due to increased viral emission or reduced infectious dose (both of which would increase the quanta emission rate) (42, 43). We note that some variants could in principle also increase transmissibility by lengthening the period of infectiousness for a given person, which by itself would not increase the quanta emission rate in a given situation. B0 is relatively well known, and varies with a susceptible person’s age, sex, and body weight, in addition to physical activity level. rE and rB are less uncertain than Ep0 and are functions of the specific physical and vocalization activities (33, 37, 44). Thus, they are useful in capturing the quantitative impact of specific controllable factors. There could be factors beyond those considered here that lead to variation of viral emissions, such as respiratory effort of patients with breathing disorders such as pneumonia or asthma (45). Such factors can be incorporated into updated Tables in the future.
Then P and Nsi can be expressed as functions of Ep0, B0 and the product of the other controllable factors as: Where H is the hazard parameter: The four terms that make up λ may vary in relative importance for different diseases and conditions. λdec ∼ 1.1 h-1 (46) has been reported for COVID-19. λdec depends on temperature and relative humidity (47, 48). λdep depends on particle size and the geometry and airflow in a given space. Respiratory particle sizes in the range from 1-5 μm are thought to play a role in aerosol transmission of COVID-19, due to a combination of high emission rates by activities such as talking (49) and low deposition rates. λdep for a typical furnished indoor space span 0.2-2 h-1 over this size range, with faster deposition for larger particles (50). λ0 varies from ∼12 h-1 for airborne infection isolation rooms (51), ∼6 h-1 for laboratories, ∼0.5 h-1 for residences (52), and ∼1 h-1 for offices (1, 53). Very little ventilation data is available for many semi-public spaces such as retail, restaurants and bars or transportation. λcle can vary from 0, if such systems are not in use, to several h-1 for adequately sized systems. Ventilation with clean outdoor air will be important in most situations, while virus decay and deposition likely contribute but are more uncertain for COVID-19, based on current information. In particular, the size distribution of aerosols containing infectious viruses is uncertain.
We consider a worst-case scenario where rates of deposition and infectivity decay are small compared with ventilation and air cleaning and can therefore be neglected. This also allows using the same risk parameter to compare different airborne diseases. This yields: H can be recast as: where: Assuming that all of the people present in a space are susceptible to infection, L is equivalent to the ventilation plus air cleaning rate per person present in the space, (typically expressed in liters s-1 person-1 in guidelines such as from refs (54–56). If some fraction of the people present are immune to the disease, then L is larger than the corresponding personal ventilation rate in the guidance. While this recasting will be useful to persons familiar with ventilation guidelines, we keep the form in Eq. 11 for most further analyses, since the number of people allowed in a space is one of the critical variables that can be examined with this risk parameter.
We define another parameter, which does not depend on activity type or face covering choice and only captures the characteristics of susceptible people’s presence in the indoor space, not of their behavior.
For the infection risk of individuals in the presence of one infector, we define the relative risk parameter, a third parameter proportional to the conditional probability of infection, Hr, as follows: When it is unknown whether an infector is present, H is an approximate indicator of the absolute probability of infection (Pa), since the expected value of number of infectors (Ni) is the product of number of occupants (N) and probability of an occupant being infectious, a measure of prevalence of infectious people in local population (ηI). Then, For a situation with multiple potential infectors present (e.g. a COVID ward in a hospital), the risk parameters should be multiplied by the number of infectors.
3. Applications
3.1. Value of the risk parameters for documented outbreaks of COVID-19
An important advantage of the simplified risk parameters is that their values can be calculated for outbreaks that are documented in the scientific literature. Values for documented COVID-19 outbreaks are shown in Table 1 (rE and rB are estimated based on the likely types of activities in each case (33, 37, 44), see Table SI-2 for typical values). Also included are values for outbreaks documented in the literature for tuberculosis and measles, which are widely accepted to transmit through the air, and an influenza outbreak that was clearly due to airborne transmission.
We see that the COVID-19 outbreaks that have been documented span ∼2.5 orders of magnitude range of the risk parameter H ∼ 0.09-30 persons h2 m-3. Numbers of secondary cases of these outbreaks generally increase with H (Figure 1a). Aiming to maintain values far below the threshold of 0.05 persons h2 m-3 should help reduce outbreaks.
Hr correlates well with the attack rates for the outbreaks reported in Table 1 (Figure 1b, Hr can be calculated for more literature outbreaks than H). COVID-19 outbreaks are observed for Hr > 0.001 h2 m-3, and thus indoor activities should be limited to conditions below this value during the pandemic whenever possible.
A trend line was fitted to the attack rate vs. Hr dataset with the Box/Wells-Riley model, with the fitting parameter being Ep0, i.e. the basic quanta shedding rate (when breathing only, no vocalization). An Ep0 of 18.6 quanta h-1 was obtained by fitting (with B0 = 0.288 m3 h-1 assumed for all occupants for simplicity). This value is higher than that suggested by Buonanno et al. (2 quanta h-1) (33, 37), but within the uncertainties provided by those authors. The attack rates estimated according to this trend line have a high correlation with the actual attack rates (r2 = 0.90; Figure SI-2).
The good agreement supports the dominant airborne character of these COVID-19 outbreaks. If the outbreaks had major components of fomite or large droplet transmission, we would expect a dependence on other parameters not considered here and overall much lower correlation with the risk parameters.
3.2. Effect of building parameters vs. human activities
The type of activity performed in each case (captured by rE and rB) contributes substantially to the difference in H between these cases. When human activities are not taken into account, the parameter H’ only spans a narrow range, 0.09-0.56. This is probably due to similar per-person ventilation rates in many public indoor spaces (on the order of a few liter s-1 person-1 (54)) and similar lengths of common events (in hours).
Similar to H, the variation in the values of Hr for the outbreaks in Table 1 is also largely due to rE and rB. If they are not taken into account, Hr for all outbreaks would vary in the narrow range 0.001 to ∼0.01 h2 m-3, as V and λ0 are building characteristics and D, as discussed above, is usually in hours. This implies that, in the presence of a single infector, reducing vocalization and/or physical intensity levels of the indoor activity is a very effective way to lower the infection risk of susceptible individuals. Reducing event length can also help, while reducing occupancy cannot in this case, as shown by Eq. 15.
It is possible that some of the most visible outbreaks are associated with super-emitter individuals, who shed virus particles at higher rates than others (45, 66–68). If that is the case, the actual Hr values for significant transmission in the presence of individuals that are not super emitters may be higher than those determined here. However, if super-emitters are important, so will be their contribution to total spread, and thus one should try to reduce the risk to reduce the probability of such events occurring.
3.3. Values of the risk parameters for outbreaks of other airborne diseases
In Table 1 and Figure 1 we also include a few reported indoor outbreaks of three other diseases with significant airborne transmission, i.e., tuberculosis, influenza, and measles. For outbreaks to have a similar number of secondary cases or attack rate, H or Hr needs to be higher for tuberculosis and influenza and lower for measles than that for COVID-19 (Figure 1). Note that many of the children present in the measles outbreak were vaccinated, but the risk parameter framework can still be applied by considering the number of susceptible children present. This difference is mainly due to differences in Ep0 (lower for tuberculosis and influenza and higher for measles; Figure 1b). A higher Ep0 for measles may indicate a larger amount of airborne measles virus in breath or a steeper dose-response curve for the measles virus than SARS-CoV-2, or both. A novel disease as contagious as measles would make almost any indoor situation prone to superspreading. On the other hand, tuberculosis and influenza are less contagious. Tuberculosis transmission is helped because untreated infected people remain contagious for years (69). The influenza outbreak occurred in an airplane without ventilation with the index case constantly coughing, and represents an extreme for this disease (64). Most influenza patients emit significantly less virus (70).
Given that both the measles and tuberculosis pathogens are widely accepted as airborne, the intermediate risk profile for COVID-19 in Figure 1b is not inconsistent with airborne transmission, contrary to frequently made arguments (71, 72). Contagiousness of a disease does not necessarily indicate the transmission route (73). Airborne diseases can vary in their contagiousness depending on parameters such as the amount of virus shed, the survival of the virus in the air, the dose-response relationship for infection, and other parameters. The only fundamental requirement is that transmission needs to be sufficient for the disease to survive as such, something COVID-19 has had no trouble with so far.
3.4. Graphical representation of relative risks of different situations
When it is not known whether infectors are present in an indoor event, all occupants need to be considered as possible infectors. We will assume that the probability of an occupant being infectious is the same as the fraction of infectious people in the local population (ηI). H indicates the risk of an outbreak. Consequently, the risk also depends on occupancy, besides vocalization level, event duration, ventilation, and mask wearing. Jones et al. (27) estimated the dependency of the infection risk on these factors in a tabulated manner similar as Table 2. However, they only did so qualitatively. Having introduced H as a risk parameter, we can assess the risk more quantitatively based on H values under different conditions (Table 2). Although the actual risk also depends on ηI and the choice of the threshold for high risk (red cells in Table 2) is subjective, the risk parameter (H value) in Table 2 seems to vary in a smaller range than the corresponding table in ref (27). We also show that being outdoors (with much better ventilation than indoors) has a greater expected benefit than Jones et al. estimated.
Note that, although occupancy has no impact on attack rate if an infector is present, occupancy affects the risk in two ways when the presence of infector(s) is unknown, i.e., i) the probability of the presence of an infector in a certain locality and ii) the size (number of secondary cases) of the outbreak if it occurs. Therefore, lowering occupancy has double benefits.
3.5. Risk evaluation for indoor spaces with pre-pandemic and mitigation scenarios
Values of the risk parameter H for some typical public spaces under pre-pandemic conditions are shown in Figure 2 and tabulated in Table SI-4. H in all pre-pandemic settings is on the order of 0.05 persons h2 m-3 or higher, implying significant risk of outbreak during the pandemic. Often, ventilation rates may be lower than official guidance due to e.g. malfunction, lack of maintenance, or attempts to save energy. Substandard ventilation is associated with substantial increases in the risk of outbreak. However, all of the pre-pandemic spaces are in a regime where they are highly sensitive to mitigation efforts. Therefore mitigation measures such as increasing ventilation or air cleaning, reducing voice volume when speaking, reducing occupancy, shortening duration of occupancy and mask wearing are required to reduce the risk of transmission in similar settings. With mitigation measures implemented, H in these settings can be lowered to the order of 0.01 persons h2 m-3, low enough to avoid major outbreaks.
3.6. Calculation of risk parameters for specific situations
The calculation of H, H’, and Hr for specific situations of interest has been implemented in the COVID-19 aerosol transmission estimator, which is freely available online (32). The estimator is a series of spreadsheets that implement the same aerosol transmission model described in ref (11) and in this paper. It allows the user to make a copy into an online Google spreadsheet or download it as a Microsoft Excel file for adaptation to the situations of interest to each user. The model can be used to estimate the risk of specific situations, to explore the reduction in transmission due to different control measures (e.g. increased ventilation, masking etc.), and to understand aerosol transmission modeling for incorporation into more complex models. The model also allows the estimation of the average CO2 concentration during an activity, as an additional indicator of indoor risk, and to facilitate the investigation of the relationship between infection risk and CO2 concentrations (38, 74). A screenshot of the estimator is shown in Figure SI-1.
In addition a sheet allows recalculating Table 2 in this paper for sets of parameters different from those used here (and shown in Table SI-2).
4. Conclusions and recommendations
We have explored the relationship between airborne infection transmission when sharing indoor spaces and the parameters of the space using a box / Wells-Riley model. We have derived an expression for the number of secondary infections, and isolated the controllable terms in this expression in two airborne transmission risk parameters, H and Hr.
We find a consistent relationship, with increasing attack rate in the known COVID-19 outbreaks as the value of Hr increases. This provides some confidence that airborne transmission is important in these outbreaks, and that the models used here capture the key processes important for airborne transmission.
When Hr is on order of 0.001 h2 m-3 and higher, individual occupants start to have significant COVID-19 infection risk. The lowest H for the major COVID-19 outbreaks in indoor settings reported in the literature is ∼0.1 person h2 m-3. H can be orders of magnitude higher for the superspreading events where most attendees were infected (e.g. the Skagit Choir (11)). However, if human activity-dependent factors are not taken into account, H’ for all the outbreaks discussed in this paper is ∼0.1-0.5 person h2 m-3. H’ values for public indoor spaces usually fall in or near this range probably due to similar per-person ventilation rates and public event durations. Substandard ventilation is associated with substantial increases in the risk of outbreak. However, all of the pre-pandemic example spaces analyzed are in a regime where they are highly sensitive to mitigation efforts.
The relative risk of COVID-19 infection falls in between two well-known airborne diseases, the more transmissible measles and the less transmissible tuberculosis. This shows that the fact that COVID-19 is less transmissible than measles does not rule out airborne transmission. These risk parameters can be applied for other airborne diseases, once some outbreaks are characterized in this framework. This approach may be useful in the design and renovation of building systems. For a novel disease that was as transmissible as measles, it would be very difficult to make any indoor activities safe.
In conclusion, our analysis shows that mitigation measures are needed whenever COVID-19 is spreading in a community to limit aerosol transmission risk in most indoor spaces. Among effective measures are reducing vocalization, avoiding intense physical activities, shortening duration, reducing occupancy, wearing high-quality well-fitting masks, increasing ventilation and applying additional virus removal measures (e.g. using HEPA filters). The use of multiple “layers of protection” is needed in many situations, while a single measure (e.g. masking) may not be able to reduce risk to low levels. We have shown that combinations of some or all of these measures are able to lower H close to 0.01 person h2 m-3, so that the expected number of secondary cases is substantially lower than 1 even in the presence of an infectious person, hence would be likely to avoid major outbreaks.
Data Availability
All data is available in the main text or the Supporting Information.
Supplementary Information
S1. Deviation of quanta concentration from steady state due to finite event duration
In case of short events where quanta concentration does not reach steady state, a correction factor, rss, can be introduced to account for the deviation of average quanta concentration (cavg) from that at steady state (c): Under the assumption of no infectious quanta in the air at the beginning of the event, cavg can be easily obtained from the integration of Eq. (1). Details of the derivation can be found elsewhere (11, 38). For a period [0, D], Inserting Eqs. (2) and (S2) into Eq. (S1) yields: The value of rss as a function of λD is shown in Fig. SI-3. rss approaches to λD/2 when λD is very small and to 1 when λD is very large, and reaches 0.6 at λD ∼ 2.
Tables
Acknowledgements
ZP and JLJ were partially supported by NSF AGS-1822664. TG was supported by ESRC ES/V010069/1.