Skip to main content

The vaccination threshold for SARS-CoV-2 depends on the indoor setting and room ventilation

Abstract

Background

Effective vaccines are now available for SARS-CoV-2 in the 2nd year of the COVID-19 pandemic, but there remains significant uncertainty surrounding the necessary vaccination rate to safely lift occupancy controls in public buildings and return to pre-pandemic norms. The aim of this paper is to estimate setting-specific vaccination thresholds for SARS-CoV-2 to prevent sustained community transmission using classical principles of airborne contagion modeling. We calculated the airborne infection risk in three settings, a classroom, prison cell block, and restaurant, at typical ventilation rates, and then the expected number of infections resulting from this risk at varying percentages of occupant immunity.

Results

We estimate the setting-specific immunity threshold for control of wild-type SARS-CoV-2 to range from a low of 40% for a mechanically ventilation classroom to a high of 85% for a naturally ventilated restaurant.

Conclusions

If vaccination rates are limited to a theoretical minimum of approximately two-thirds of the population, enhanced ventilation above minimum standards for acceptable air quality is needed to reduce the frequency and severity of SARS-CoV-2 superspreading events in high-risk indoor environments.

Peer Review reports

Background

Control of infectious disease is achieved when outbreaks cannot be sustained, and transmission becomes sporadic in nature. For airborne contagion in shared indoor atmospheres, Wells [1] established that the rate of transmission is inversely proportional to the ventilation rate per susceptible occupant. It then follows that to control airborne contagion we can either increase ventilation, or its equivalent through air filtration or disinfection, or decrease the number of susceptible occupants through vaccination [2]. During the COVID-19 pandemic, lockdowns and occupancy controls have been widely applied to reduce transmission of SARS-CoV-2. These are blunt but effective methods of increasing the ventilation rate per susceptible occupant of indoor spaces. As SARS-CoV-2 vaccines become increasingly available worldwide, the question becomes: at what point is the number of susceptibles in public spaces low enough so that occupancy limitations are no longer necessary to control the spread of the virus?

To address this question, we must consider the primary settings of SARS-CoV-2 transmission. As with other agents of airborne contagion such as Mycobacterium tuberculosis, SARS-CoV-2 thrives in congregate living and working spaces with shared air, such as prisons, schools, restaurants, abattoirs, and care homes. The COVID-19 pandemic is also fueled by superspreading events in crowded indoor environments where people vocalize and cannot reliably wear masks. For example, Chang et al. [3] modeled full-service restaurants to produce by far the largest increase in infections upon reopening after lockdown of any non-residential location that people visit. Estimates of necessary vaccination rates for these high-risk community settings should be protective in other microenvironments, and therefore approximate a vaccination threshold to control SARS-CoV-2.

The aim of this paper is to estimate setting-specific vaccination thresholds for the wild-type, original strain of SARS-CoV-2 using classical principles of airborne contagion modeling. We included modeling scenarios for a prison cell block and a full-service restaurant, two settings known to be high risk for SARS-CoV-2 transmission. To compare the vaccination threshold for wild-type SARS-CoV-2 to historical estimates for measles virus, we also included a classroom scenario in our analysis. A secondary aim is to quantify how vaccination and ventilation together reduce the pool of potential infectors in each of the settings by estimating the minimum viral emission rate needed to reproduce infection at varying levels of susceptibility.

Materials and methods

Approach and definitions

To develop our estimates, we defined a representative exposure scenario for each of the three settings (classroom, prison, restaurant) involving one infectious occupant in a room of typical geometry. We used an established airborne infection risk model to calculate the individual risk of infection (R) for each susceptible occupant, and the event reproduction number (Revent) at varying ventilation rates and number of susceptibles. Revent is the expected number of new infections arising from a single infectious occupant at an event [4]. This is distinct from the basic reproduction number (R0), defined as the average number of new infections resulting from the introduction of a single infectious individual into a fully (100%) susceptible host population [5]. For modeling purposes, we quantified the number of susceptibles as the percent of the total occupants who are susceptible to infection (i.e., not immune from vaccination or prior infection). We use the term area concentration of susceptibles to represent the area of indoor space (square meters [m2]) per susceptible occupant. The threshold number of susceptibles and the threshold area concentration of susceptibles occur at a calculated Revent of one, above which an average of at least one new case is expected. For each setting we calculated these two threshold values at a mechanical ventilation rate based on American National Standards Institute (ANSI)/American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) 62.1 standards for acceptable air quality [6], and at a natural ventilation rate when windows cannot be opened and air exchange results solely from infiltration through the building envelope. We then determined the vaccination threshold as the complement of the threshold number of susceptibles assuming no immunity from prior infection and that vaccination confers complete protection. For comparative purposes, we also calculated the threshold values in each setting at a ventilation rate of 15 L per second per person (L s−1 p−1), a typical goal for high indoor air quality consistent with EN 15251 Category I criteria for a non low polluting building [7].

Calculation of the event reproduction number (R event )

We used the Gammaitoni and Nucci [8] equation coupled with a Poisson dose–response model to calculate Revent for SARS-CoV-2 in a prototypical classroom, prison cell block, and full-service restaurant. The first step is calculating the probability of infection (PI) resulting from each exposure through Eqs. (13):

$$n\left( {t,ER_{q} } \right) = \frac{{ER_{q} \cdot I}}{IVRR \cdot V} \cdot \left( {1 - e^{ - IVRR \cdot t} } \right)\quad \left( {{\text{quanta m}}^{{ - {3}}} } \right)$$
(1)
$$D_{q} \left( {ER_{q} } \right) = IR\int_{0}^{T} {n\left( t \right)dt} \quad \left( {{\text{quanta}}} \right)$$
(2)
$$P_{I} = 1 - e^{{ - D_{q} }} \quad \left( \% \right)$$
(3)

where n represents the quanta (infectious dose for 63% of susceptible occupants by droplet nuclei inhalation) concentration in air at time t, ERq is the quanta emission rate (quanta h−1), I is the number of infectious occupants (assumed to be only one), V is the volume of the indoor environment considered (m3), IVRR (h−1) represents the infectious virus removal rate in the space investigated, Dq is the dose of quanta inhaled by susceptible occupants, T is the total time of the exposure (h), and PI is the probability of infection of a susceptible occupant. The infectious virus removal rate is the sum of the air exchange rate (AER) via ventilation in units of air changes per hour, the particle deposition on surfaces (kd, e.g. via gravitational settling), and the viral inactivation in ambient air (λ).

With all other parameters held constant, the probability of infection calculated in Eq. (3) assumes different values based on ERq. To evaluate the individual risk (R) of infection of an exposed susceptible occupant for a given exposure scenario, we then quantify the probability of infection as a function of ERq (PI[ERq]) and the probability of occurrence of each ERq value (PERq) which can be defined by the probability density function (pdfERq) of ERq assuming a lognormal distribution. Since the probability of infection (PI[ERq]) and the probability of occurrence PERq are independent events, R for a given ERq, R(ERq), can be evaluated as the product of the two terms:

$$R\left( {ER_{q} } \right) = P_{I} \left( {ER_{q} } \right) \cdot P_{ERq} \quad \left( \% \right)$$
(4)

where PI(ERq) is the conditional probability of the infection, given a certain ERq, and PERq represents the relative frequency of the specific ERq value. The individual risk (R) of an exposed susceptible occupant is then calculated by integrating the pdfR for all possible ERq values, i.e. summing up the R(ERq) values calculated in Eq. (5):

$$R = \int_{{ER_{q} }} {R\left( {ER_{q} } \right)dER_{q} } = \int_{{ER_{q} }} {\left( {P_{I} \left( {ER_{q} } \right) \cdot P_{ERq} } \right)dER_{q} } \quad \left( \% \right)$$
(5)

Equation (5) represents a numerical solution approximately equaling the average PI that would result from a Monte Carlo simulation randomly sampling ERq from its lognormal distribution. The individual risk R also represents the ratio between the number of new infections and the number of exposed susceptible occupants (S) for a given exposure scenario and considering all possible ERq values from its lognormal distribution for the infectious occupant under investigation. For a single exposure event involving a single infectious occupant, Revent is calculated as the product of R and S as in Eq. (6):

$$R_{event} = R \cdot S\quad \left( {{\text{infections}}} \right)$$
(6)

For a specific event, the threshold number of susceptibles occurs at the value of S where Revent equals one (Sthreshold = 1/R) and is calculated by dividing Sthreshold by the total room occupancy less the infected occupant. The threshold area concentration of susceptibles is calculated by dividing Sthreshold by the room area.

Modeling scenario input parameters

Input parameters for the classroom, prison, and restaurant scenarios are summarized in Table 1. Geometry and default occupancy for the classroom are based on the rooms studied by Wells [1] with an exposure time of 5.5 h representing a single day. The restaurant model encompasses the dining room geometry of the US Department of Energy building prototype for a full-service restaurant, with an exposure time of 1.5 h [9]. The prison model is based on the largest cell block size studied by Hoge et al. [10], which was overcrowded with a median living area of 3.2 m2 per inmate. The exposure time for the prison scenario is likely highly variable, but we assume it to be 36 h since inmates share the same airspace for extended time periods and peak infectiousness has been estimated to occur at 2 days before to 1 day after symptom onset [11]. Thus, a 36-h period where infectiousness is at or near peak but without symptoms that would prompt quarantine can be reasonably expected.

Table 1 Modeling input and ventilation reference parameters

The distributions for the quanta emission rate were modified from Buonanno et al. for wild-type SARS-CoV-2 ([12]; see Additional file 1) for standing and speaking for the classroom (assuming the class instructor is the emitting subject), resting and loudly speaking for the restaurant, and resting and oral breathing for the prison, with the log10 average ERq values indicated in Table 1 and a log10 standard deviation for all distributions of 1.2. All susceptible occupants were assumed to be at rest with an inhalation rate of 0.49 m3 h−1. We used a deposition rate, kd, of 0.24 h−1 based on the ratio between the settling velocity of super-micrometer particles (roughly 1.0 × 10–4 m s−1 [13]) and the height of the emission source (1.5 m). For the SARS-CoV-2 inactivation rate in air, we used a value of 0.63 h−1 based on the measurements reported by van Doremalen et al. [14]. For each scenario, we varied the AER from zero to a maximum of six air changes per hour to calculate R and Revent at a number of susceptibles ranging from 0 to 100%.

Results

The results of our modeling analysis are summarized in Fig. 1 and Table 2 for each setting at an assumed natural ventilation rate of 0.5 air changes per hour, and at a mechanical ventilation rate corresponding to the applicable standard for acceptable air quality based on ANSI/ASHRAE 62.1 and shown in Table 1 [6]. The naturally ventilated restaurant (Fig. 1A) has the lowest threshold number of susceptibles of 15%, and the mechanically ventilated classroom (Fig. 1B) has the highest threshold number of susceptibles of 60%. The threshold number of susceptibles for the prison cell block (Fig. 1C) exhibits the smallest difference between the natural ventilation (23%) and mechanical ventilation (31%) scenarios.

Fig. 1
figure 1

Surface graphs of Revent for wild-type SARS-CoV-2 as a function of the number of susceptibles and air exchange rate (AER) for the restaurant (A), classroom (B) and prison cell block (C) modeling scenarios. Contour lines connect equal Revent values. The black- and white-filled points along the Revent = 1.0 contour line identify the threshold number of susceptibles for natural ventilation and mechanical ventilation scenarios, respectively, at the intersection of the dashed horizontal and vertical lines. The threshold values are labeled in parenthesis in terms of both the percent susceptible and m2 susceptible−1

Table 2 Modeling results

The average threshold number of susceptibles for all three settings calculated under the natural and mechanical ventilation rates is 32%. In the absence of immunity from prior infections and assuming vaccination confers complete protection, these results suggest an average vaccination threshold of 68% with a range of 40–85%. The naturally ventilated prison and restaurant have the highest threshold area concentration of susceptibles at 14 m2 susceptible−1, while the mechanically ventilated classroom has the lowest at approximately 5.0 m2 susceptible−1. The overall average threshold area concentration of susceptibles for mechanical and natural ventilation is approximately 10 m2 susceptible−1.

Increasing the ventilation rate to the high air quality metric of 15 L s−1 p−1 increases the threshold number of susceptibles to 95% in the classroom, 60% in the prison, and 44% in the restaurant. The average threshold number of susceptibles for all three settings becomes 66%, more than twice the average of the natural and mechanical ventilation scenarios. To maintain an Revent of one in a fully susceptible population, the estimated ventilation requirements are 43 L s−1 p−1 (24 air changes per hour), 30 L s−1 p−1 (9.5 air changes per hour) and 17 L s−1 p−1 (7.0 air changes per hour) for the restaurant, prison, and classroom, respectively. Such high air exchange rates are impracticable in most settings, suggesting a role for ultraviolet air disinfection [15, 16].

Increasing ventilation and/or decreasing the number of susceptibles has the effect of increasing the minimum ERq necessary to produce an Revent of one, thereby reducing the number of infected occupants capable of infecting others on average. This is illustrated in Fig. 2 for the prison cell block model. For the naturally ventilated cell block in a fully susceptible population, the minimum ERq is just below 1.0 quanta h−1, occurring at the 58th percentile value of the resting, oral breathing distribution. At a number of susceptibles of 23%, the minimum ERq becomes approximately 4.3 quanta h−1 at the 78th percentile value. Increasing ventilation to 15 L s−1 p−1 further decreases the pool of potential infectors, raising the minimum ERq to approximately 17 quanta h−1 at the 90th percentile value, indicating only a 10% chance of a secondary infection.

Fig. 2
figure 2

Minimum quanta emission rates (ERq) for Revent ≥ 1.0 for the prison scenario under natural ventilation, mechanical ventilation, and high air quality ventilation conditions as a function of the number of susceptibles. Points #1 and #2 identify the minimum emission rates for high air quality ventilation and natural ventilation at their respective threshold number of susceptibles from Fig. 1C. Point #3 identifies the minimum emission rate for high air quality ventilation at the natural ventilation threshold number of susceptibles, representing both high ventilation and high vaccination. The minimum emission values are labeled in parenthesis, denoting the emission in quanta h−1 and its corresponding percentile in the resting, oral breathing ERq distribution

Discussion

The overall average threshold number of susceptibles calculated for the natural and mechanical ventilation scenarios is 32%. This is similar to the threshold inferred from a R0 value of approximately 3, as was estimated for the initial wild-type SARS-CoV-2 outbreaks in Wuhan, China and Northern Italy [17]. Our analysis is also consistent with the overdispersed epidemiological nature of SARS-CoV-2 [18], with a minority of cases accounting for most secondary transmissions. In the naturally ventilated prison, we calculate that emissions approximately below the 60th percentile value will fail to reproduce infection, on average, indicating the median emission is not a significant source of transmission (Fig. 2). Furthermore, application of Eq. (5) for the prison scenario shows that emissions above the 80th percentile value account for at least 85% of the total individual risk, suggesting a dispersion parameter (k) between 0.10 and 0.16. This derivation is provided in the Additional file 1 and enables quantification of the probability of SARS-CoV-2 superspreading and outbreak extinction as defined by Lloyd-Smith et al. [19]. Due to this overdispersion, vaccinating 77% of inmates in a naturally ventilated cell block still leaves the remaining susceptible population vulnerable to emitters above the 78th percentile. As a result, explosive but comparatively rare superspreading events may continue in crowded, poorly ventilated settings, a phenomenon that challenges the eradication of measles virus [20]. For reference, exceedance probability plots for the modeled probability of infection results are provided in the Additional file 1 to allow further assessment of this overdispersion.

Applying both high vaccination and high ventilation raises both the threshold number of susceptibles and the minimum emission rate needed to reproduce infection, decreasing the dispersion parameter and increasing the probability of outbreak extinction. Uniformly increasing ventilation to a high air indoor air quality metric of 15 L s−1 p−1 approximately doubles the average threshold number of susceptibles and therefore halves vaccination requirements for equivalent prevention of infection. Thus, while a ventilation rate of 15 L s−1 p−1 is unlikely to prevent all secondary infections when a high-emitting index case is introduced into a fully susceptible, indoor population [21], it can provide a substantial downstream epidemiological benefit relative to a poorly ventilated baseline condition. This effect is important for pathogens where transmission is overdispersed, with M. tuberculosis being another example [22], as superspreading events (SSEs) facilitate infection of the high-emitting minority that continues the chain of contagion. For our prison cell block model, we estimate that increasing the natural ventilation rate to the high air quality ventilation rate decreases the SSE probability from 16 to 6.6% (see Additional file 1). This is an important finding, as prisons and jails are clear hot spots for SARS-CoV-2 transmission. For example, by March 2021, five California State Prisons (Chuckawalla Valley, California Rehabilitation Center, Avenal, San Quentin, and California Men’s Colony) reported total confirmed COVID-19 case rates above 800 per 1000 inmates [23]. Such high case rates imply a low threshold number of susceptibles, with inadequate ventilation a likely factor. Indeed, during an investigation of the San Quentin State Prison in June 2020, McCoy et al. [24] noted cell blocks with windows that were welded shut and with fan systems that appeared to have been inactive for years.

Historical examples for measles virus illustrating the relationship between ventilation and the threshold number of susceptibles in classrooms are provided by Wells [1, 25] and Thomas [26]. In classic experiments using upper-room air irradiation in primary and upper school classrooms during the 1941 outbreak of measles in suburban Philadelphia, USA, Wells estimated a threshold number of susceptibles of approximately 20% in unirradiated rooms at a then-standard ventilation rate of 14 L s−1 p−1. Irradiated classrooms supported a much higher threshold number of susceptibles of approximately 57% because the weekly probability of infection in the irradiated rooms was approximately four to five times lower than in the unirradiated rooms [1, 25]. The findings of Wells are similar to those of Thomas [26] who studied the spread of measles in primary schools in the Woolwich district of London in 1904. Thomas concluded that outbreaks of measles tend to occur when the number of susceptibles exceeds approximately 33% and generally continue until the proportion is reduced to 18%. However, the spread of measles in the Woolwich classrooms below the 33% threshold was highly heterogeneous, with many experiencing significant outbreaks infecting a majority of susceptible occupants. The three classes with a number of susceptibles below 10% experienced zero cases of measles, and two temporary schools with crowding and poor ventilation had explosive outbreaks that nearly exhausted the population of susceptibles, with a median probability of infection of 87% for the five classes in the two schools. Thomas measured a carbon dioxide (CO2) concentration of 3000 parts per million in one of the temporary schools [26], indicating a steady-state ventilation rate below 2 L s−1 p−1 and comparable to our natural ventilation scenario. The higher contagiousness of measles as compared to SARS-CoV-2 is illustrated by the historical reported threshold number of susceptibles of 20–33% as compared to our classroom estimate of 37–60% despite the lower ventilation standards of present day. This difference is also reflected by the median classroom measles probability of infection of 87% for the poorly ventilated temporary schools studied by Thomas [26] as compared to the individual risk (R) of approximately 14% we calculated for wild-type SARS-CoV-2 (Table 2). A ventilation rate of 14 L s−1 p−1 appears sufficient, on average, to prevent sustained airborne transmission of wild-type SARS-CoV-2 in a classroom with a number of susceptibles up to approximately 90%.

The concept of a threshold number of susceptibles for airborne contagion was also evaluated by Kelker [27] (further described by Gorham [28]), based on induced outbreaks of canine distemper virus (CDV) in small ferret populations. CDV is a highly contagious morbillivirus affecting a range of carnivores, and infection of ferrets with CDV is a common animal model for measles pathogenesis [29]. Kelker [27] determined that unless ferrets were densely packed together, with 70% of the population immune (30% number of susceptibles), and with the immune animals spread uniformly throughout the population, outbreaks caused by the introduction of a single CDV-infected ferret had high extinction probability and faded after only a small fraction of susceptibles were infected. In other words, it was difficult to initiate and sustain an outbreak in a population less than 30% susceptible to CDV [28].

Limitations

A limitation of our infection risk modeling approach is the assumption of a homogeneous concentration of droplet nuclei within the room, with viral emissions being instantaneously and completely mixed. In reality, the risk presented by inhalable virus-laden particles will be substantially higher for a susceptible person directly breathing the respiratory jet in close proximity (generally < 2 m) to an infected host [30, 31], and air speed, local flow direction, and the velocity of the emission itself (for example, a cough versus a normal exhalation) all affect the fate of the jet and the resulting zone of elevated risk relative to the completely mixed condition. While this is a significant limitation, the use of a homogeneous exposure-point concentration is standard practice in environmental risk assessment [32], and it is generally not possible to predict exact locations of infected hosts and their proximity to susceptible occupants of a shared indoor atmosphere, and their relative positions change in time. Similarly, local airflow conditions are complex, variable, and cannot be accurately represented without computational fluid dynamics (CFD) modeling or site-specific tracer testing. Thus, the most defensible approach, also used in environmental risk assessment [32], is to assume susceptible persons are equally exposed to a uniform, average concentration within the exposure unit (the room), fully knowing that some individuals inhale a larger dose than others. From this perspective, the complete and instantaneous mixing of the viral emission into a room is a quite reasonable approach to generate an average, time-variable exposure point concentration. Supporting this, a recent comparison of this box-modeling approach with CFD simulations for a classroom environment indicates relatively minor errors for natural (6%) and forced mechanical (29%) ventilation scenarios [33]. The uncertainty in the emission rate, based on viral loads that vary several orders of magnitude between individuals and over time [34], is likely much more significant than that caused by incomplete mixing at the small scale of our models. For example, the existence of superemitters who generate and emit substantially more particles than the average person, independent of respiratory tract infection, has been demonstrated for breathing [35], the cough [36], and speech [37]. While undoubtedly important, high viral load is not a sine qua non for contagiousness, just as sputum smear status is not for M. tuberculosis [38]. As such, further improvements to the emission rate distributions are needed that incorporate variation in droplet volume concentrations [34, 39], such that a more complete stochastic emission model can be implemented.

An additional limitation is our estimation of vaccination thresholds using singular, setting-specific events, without considering cumulative exposure effects that may result from an infectious person attending class in 2 successive days, for example. The importance of singular SSEs on SARS-CoV-2 transmission is well established, and such events likely occur during a narrow 1–2 day window of peak infectivity [40]. As such we do not expect cumulative exposures to be a significant factor outside of co-habitation environments, which is why our prison scenario used a 36-h duration. Our approach also does not account for extreme examples such as someone visiting multiple similar restaurants for similar durations on the same evening (thus increasing the number of exposed susceptibles to a similar infectious dose), or for a bartender or other vocalizing restaurant employee who may be present for much longer than 1.5 h. Indeed, there are numerous other scenarios, such as choirs or high-intensity exercise rooms, where higher vaccination thresholds are likely, reinforcing the need for high levels of both vaccination and ventilation also considering that vaccines are not 100% protective.

The emergence of the more transmissible Delta variant in 2021 led to a worldwide resurgence of SARS-CoV-2 and revealed waning efficacy of the first-generation vaccines at preventing clinical infection. A case study [41] of an outbreak in a highly vaccinated (79%) prison illustrates the challenges presented by the Delta variant, particularly in a highly contagious environment such as a prison. Attack rates during the outbreak were as follows:

  • 100% for unvaccinated persons without a prior SARS-CoV-2 infection (35 of 35);

  • 78% for vaccinated persons without a prior SARS-CoV-2 infection (128 of 164);

  • 57% for unvaccinated persons with a prior SARS-CoV-2 infection (4 of 7); and

  • 5% for vaccinated persons with a prior SARS-CoV-2 infection (1 of 21).

For this outbreak, a spectrum of susceptibility was observed, with by far the least susceptible group being inmates who were fully vaccinated and had a prior infection. This spectrum is further complicated by different attack rates observed depending on the type of vaccine received by the inmates [41]. The threshold number of susceptibles for this outbreak was extremely low, as all fully naïve inmates were infected. A high attack rate for the Delta variant (50%) was also observed in a classroom setting similar to the scenario modeled herein, with the class instructor as the infected host [42]. When considering our model does not include close proximity interactions, the more transmissible Delta variant, and the observed variable vaccine efficacy and resulting spectrum of susceptibility, it is clear that our estimates reflect a lower bound. Furthermore, in light of the extreme contagiousness of the Delta variant, more emphasis should be placed on improving air flow distribution to minimize the formation of stagnation zones, instead of simply increasing airflow in its current configuration [43].

Conclusions

Our fully prospective airborne infection modeling results are consistent with the transmission dynamics of SARS-CoV-2 and illustrate the challenges presented by substantial heterogeneity in the settings of contagion and a skewed viral emission rate distribution. To support pre-pandemic levels of occupancy, required vaccination rates are much higher for a naturally ventilated restaurant (85%) than for a mechanically ventilated classroom (40%). These estimates reflect a lower bound as they are based on wild-type SARS-CoV-2 and do not consider the more transmissible variants of concern, such as Delta. Regardless, as vaccination campaigns progress around the world, it follows that occupancy limitations should be relaxed for classrooms before full-service indoor restaurants. Maintaining focus on enhanced ventilation together with vaccination is especially important considering the more transmissible Delta variant which is associated with increasing possibility of second infections or vaccine breakthrough infections. Avoidance of overcrowding remains a critical strategy to minimize airborne transmission, as our calculations suggest ensuring an average of 10 m2 per susceptible occupant of an indoor space is approximately equivalent to achieving a number of susceptibles of 32% of normal occupancy. This is because the ventilation rate per susceptible occupant is more than tripled relative to the baseline average occupant loading of 2.7 m2 per susceptible occupant for the three settings evaluated herein.

Availability of data and materials

All data generated or analysed during this study are included in this published article and its Additional file 1.

References

  1. Wells WF. Airborne contagion and air hygiene. Cambridge: Harvard University Press; 1955.

    Google Scholar 

  2. Riley RL. Prevention and control of airborne infection in the community. In: Annals of the New York Academy of Sciences. Airborne contagion, vol. 353. New York: The New York Academy of Sciences; 1980. p. 331–9.

    Google Scholar 

  3. Chang S, Pierson E, Koh PW, et al. Mobility network models of COVID-19 explain inequities and inform reopening. Nature. 2021;589:82–7.

    CAS  Article  Google Scholar 

  4. Tupper P, Boury H, Yerlanov M, Colijn C. Event-specific interventions to minimize COVID-19 transmission. Proc Natl Acad Sci USA. 2020;117:32038–45.

    CAS  Article  Google Scholar 

  5. Anderson RM, May RM. Infectious diseases of humans: dynamics and control. New York: Oxford University Press; 1991.

    Google Scholar 

  6. ASHRAE, ANSI/ASHRAE Standard 62.1-2019. Ventilation for acceptable indoor air quality. Atlanta: American Society of Heating, Refrigerating, and Air-Conditioning Engineers, Inc.; 2019. p. 92.

  7. EN 15251. Indoor environmental input parameters for design and assessment of energy performance of buildings addressing indoor air quality, thermal environment, lighting and acoustics. CEN 2007, Brussels.

  8. Gammaitoni L, Nucci MC. Using a mathematical model to evaluate the efficacy of TB control measures. Emerg Infect Dis. 1997;3(3):335–42.

    CAS  Article  Google Scholar 

  9. Thornton BA, Rosenberg MI, Richman EE, et al. Achieving the 30% goal: energy and cost savings analysis of ASHRAE Standard 90.1-2010. Richland: Pacific Northwest National Laboratory (PNNL); 2011. PNNL-20405.

  10. Hoge CW, Reichler MR, Dominguez EA, et al. An epidemic of pneumococcal disease in an overcrowded, inadequately ventilated jail. N Engl J Med. 1994;331(10):643–8.

    CAS  Article  Google Scholar 

  11. He X, Lau EHY, Wu P, et al. Author correction: temporal dynamics in viral shedding and transmissibility of COVID-19. Nat Med. 2020;26:1491–3.

    CAS  Article  Google Scholar 

  12. Buonanno G, Morawska L, Stabile L. Quantitative assessment of the risk of airborne transmission of SARS-CoV-2 infection: prospective and retrospective applications. Environ Int. 2020;145:106112.

    CAS  Article  Google Scholar 

  13. Chatoutsidou SE, Lazaridis M. Assessment of the impact of particulate dry deposition on soiling of indoor cultural heritage objects found in churches and museums/libraries. J Cult Herit. 2019;39:221–8.

    Article  Google Scholar 

  14. van Doremalen N, Bushmaker T, Morris DH, et al. Aerosol and surface stability of SARS-CoV-2 as compared with SARS-CoV-1. N Engl J Med. 2020;382(16):1564–7.

    Article  Google Scholar 

  15. Riley RL, Nardell EA. Clearing the air. The theory and application of ultraviolet air disinfection. Am Rev Respir Dis. 1989;139(5):1286–94.

    CAS  Article  Google Scholar 

  16. Nardell EA, Nathavitharana RR. Airborne spread of SARS-CoV-2 and a potential role for air disinfection. JAMA. 2020;324(2):141–2.

    CAS  Article  Google Scholar 

  17. D’Arienzo M, Coniglio A. Assessment of the SARS-CoV-2 basic reproduction number, R0, based on the early phase of COVID-19 outbreak in Italy. Biosaf Health. 2020;2(2):57–9.

    Article  Google Scholar 

  18. Adam DC, Wu P, Wong JY, et al. Clustering and superspreading potential of SARS-CoV-2 infections in Hong Kong. Nat Med. 2020;26:1714–9.

    CAS  Article  Google Scholar 

  19. Lloyd-Smith JO, Schreiber SJ, Kopp PE, et al. Superspreading and the effect of individual variation on disease emergence. Nature. 2005;438(7066):355–9.

    CAS  Article  Google Scholar 

  20. Langmuir AD. Changing concepts of airborne infection of acute contagious diseases: a reconsideration of classic epidemiologic theories. In: Annals of the New York Academy of Sciences. Airborne contagion, vol. 353. New York: The New York Academy of Sciences; 1980. p. 35–44.

    Google Scholar 

  21. Nardell EA, Keegan J, Cheney SA, Etkind SC. Airborne infection. Theoretical limits of protection achievable by building ventilation. Am Rev Respir Dis. 1991;144(2):302–6.

    CAS  Article  Google Scholar 

  22. Melsew YA, Gambhir M, Cheng AC, et al. The role of super-spreading events in Mycobacterium tuberculosis transmission: evidence from contact tracing. BMC Infect Dis. 2019;19:244.

    Article  Google Scholar 

  23. California Department of Corrections and Rehabilitation. Population COVID-19 tracking. https://www.cdcr.ca.gov/covid19/population-status-tracking/. Accessed 2 Mar 2021.

  24. McCoy S, Bertozzi SM, Sears D, et al. Urgent memo: COVID-19 outbreak: San Quentin Prison. 2020. Available at https://amend.us/wp-content/uploads/2020/06/COVID19-Outbreak-SQ-Prison-6.15.2020.pdf.

  25. Wells WF. Air disinfection in day schools. Am J Public Health. 1943;33:1436–43.

    CAS  Article  Google Scholar 

  26. Thomas CJ. Measles in the Woolwich district. London County Council, appendix to report of the Medical Officer of Health for London County Council for the year 1904. 1905;46–60. Available at https://wellcomelibrary.org/item/b18252539.

  27. Kelker D. The effect of immunes on the spread of distemper in small ferret populations. Comput Biol Med. 1980;10(1):53–60.

    CAS  Article  Google Scholar 

  28. Gorham JR. Some experiments and field observations of distemper in mink and ferrets. Adv Vet Med. 1999;41:557–70.

    CAS  Article  Google Scholar 

  29. de Vries RD, Ludlow M, de Jong A, Rennick LJ, Verburgh RJ, van Amerongen G, van Riel D, van Run P, Herfst S, Kuiken T, et al. Delineating morbillivirus entry, dissemination and airborne transmission by studying in vivo competition of multicolor canine distemper viruses in ferrets. PLoS Pathog. 2017. https://doi.org/10.1371/journal.ppat.1006371.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Cortellessa G, Stabile L, Arpino F, Faleiros DE, van den Bos W, Morawska L, Buonanno G. Close proximity risk assessment for SARS-CoV-2 infection. Sci Total Environ. 2021;794:148749.

    CAS  Article  Google Scholar 

  31. Fierce L, Robey AJ, Hamilton C. Simulating near-field enhancement in transmission of airborne viruses with a quadrature-based model. Indoor Air. 2021;31:1843–59.

    CAS  Article  Google Scholar 

  32. United States Environmental Protection Agency (USEPA). Calculating upper confidence limits for exposure point concentrations at hazardous waste sites. Washington DC: Office of Emergency and Remedial Response. 2002;OSWER 9285.6-10. Downloaded from https://nepis.epa.gov/Exe/ZyPURL.cgi?Dockey=P100CYCE.txt.

  33. Foster A, Kinzel M. Estimating COVID-19 exposure in a classroom setting: a comparison between mathematical and numerical models. Phys Fluids. 2021;33:021904.

    CAS  Article  Google Scholar 

  34. Chen PZ, Bobrovitz N, Premji Z, Koopmans M, Fisman DN, Gu FX. Heterogeneity in transmissibility and shedding SARS-CoV-2 via droplets and aerosols. eLife. 2021;10:e65774.

    CAS  Article  Google Scholar 

  35. Edwards DA, Man JC, Brand P, Katstra JP, Sommerer K, Stone HA, Nardell E, Scheuch G. Inhaling to mitigate exhaled bioaerosols. Proc Natl Acad Sci USA. 2004;101(50):17383–8.

    CAS  Article  Google Scholar 

  36. Zayas G, Chiang MC, Wong E, MacDonald F, Lange CF, Senthilselvan A, King M. Cough aerosol in healthy participants: fundamental knowledge to optimize droplet-spread infectious respiratory disease management. BMC Pulm Med. 2012;12:11.

    Article  Google Scholar 

  37. Asadi S, Wexler AS, Cappa CD, Barreda S, Bouvier NM, Ristenpart WD. Aerosol emission and superemission during human speech increase with voice loudness. Sci Rep. 2019;9(1):2348.

    Article  Google Scholar 

  38. Fennelly KP, Jones-Lopez EC, Ayakaka I, Kim S, Menyha H, Kirenga B, Muchwa C, Joloba M, Dryden-Peterson S, Reilly N, et al. Variability of infectious aerosols produced during coughing by patients with pulmonary tuberculosis. Am J Respir Crit Care Med. 2012;186(5):450–7.

    Article  Google Scholar 

  39. Edwards DA, Ausiello D, Salzman J, et al. Exhaled aerosol increases with COVID-19 infection, age, and obesity. Proc Natl Acad Sci USA. 2021;118:e2021830118.

    CAS  Article  Google Scholar 

  40. Goyal A, Reeves DB, Cardozo-Ojeda EF, Schiffer JT, Mayer BT. Viral load and contact heterogeneity predict SARS-CoV-2 transmission and super-spreading events. eLife. 2021;10:e63537.

    CAS  Article  Google Scholar 

  41. Hagan LM, McCormick DW, Lee C, Sleweon S, Nicolae L, Dixon T, Banta R, Ogle I, Young C, Dusseau C et al. Outbreak of SARS-CoV-2 B.1.617.2 (Delta) variant infections among incarcerated persons in a federal prison—Texas, July–August 2021. Morb Mortal Wkly Rep. 2021;70(38):1349–54.

  42. Lam-Hine T, McCurdy SA, Santora L, Duncan L, Corbett-Detig R, Kapusinszky B, Willis M. Outbreak associated with SARS-CoV-2 B.1.617.2 (Delta) variant in an elementary school—Marin County, California, May–June 2021. Morb Mortal Wkly Rep. 2021;70(35):1214–9.

    CAS  Article  Google Scholar 

  43. Khankari K. Analysis of spread of airborne contaminants and risk of infection. ASHRAE J. 2021;63(7):14–16,18–20.

Download references

Acknowledgements

The authors thanks Chantal Labbé at QUT ILAQH for her invaluable research support.

Funding

This manuscript was prepared without external funding.

Author information

Authors and Affiliations

Authors

Contributions

AM, LS, GB, LM designed research; AM performed research; AM, LS, GB, LM analyzed data; AM wrote the paper; and LS, GB, LM reviewed the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to L. Morawska.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors have no conflicts or competing interests to disclose.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1:

Supplementary Material.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Mikszewski, A., Stabile, L., Buonanno, G. et al. The vaccination threshold for SARS-CoV-2 depends on the indoor setting and room ventilation. BMC Infect Dis 21, 1193 (2021). https://doi.org/10.1186/s12879-021-06884-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-021-06884-0

Keywords

  • SARS-CoV-2
  • Airborne transmission
  • Ventilation
  • Vaccination