Monitoring carbon dioxide to quantify the risk of indoor airborne transmission of COVID-19

Martin Z. Bazant, Ousmane Kodio, Alexander E. Cohen, Kasim Khan, Zongyu Gu, John W. M. Bush
doi: https://doi.org/10.1101/2021.04.04.21254903

Abstract

A new guideline for mitigating indoor airborne transmission of COVID-19 prescribes a limit on the time spent in a shared space with an infected individual (Bazant and Bush, 2021). Here, we rephrase this safety guideline in terms of occupancy time and mean exhaled carbon dioxide concentration in an indoor space, thereby enabling the use of CO2 monitors in the risk assessment of airborne transmission of respiratory diseases. While CO2 concentration is related to airborne pathogen concentration (Rudnick and Milton, 2003), the guideline developed here accounts for the different physical processes affecting their evolution, such as enhanced pathogen production from vocal activity and pathogen removal via face-mask use, filtration, sedimentation and deactivation. Critically, transmission risk depends on the total infectious dose, so necessarily depends on both the pathogen concentration and exposure time. The transmission risk is also modulated by the fractions of susceptible, infected and immune persons within a population, which evolve as the pandemic runs its course. A mathematical model is developed that enables a prediction of airborne transmission risk from real-time CO2 measurements. Illustrative examples of implementing our guideline are presented using data from CO2 monitoring in university classrooms and office spaces.

Impact Statement

Impact Statement There is mounting scientific evidence that COVID-19 is primarily transmitted through indoor airborne transmission, as arises when a susceptible person inhales virus-laden aerosol droplets exhaled by an infectious person. A safety guideline to limit indoor airborne transmission (Bazant and Bush, 2021) has recently been derived that complements the public health guidelines on surface cleaning and social distancing. We here recast this safety guideline in terms of total inhaled carbon dioxide, as can be readily monitored in most indoor spaces. Our approach paves the way for optimizing air handling systems by balancing health and financial concerns, informs policy for safely re-opening schools and businesses as the pandemic runs its course, and may be applied quite generally in the mitigation of airborne respiratory illnesses, including influenza.

1. Introduction

Coronavirus disease 2019 (COVID-19) has caused a devastating global pandemic since it was first identified in Wuhan, China in December 2019 (Chen et al., 2020; Li et al., 2020). For over a year, public health guidance has focused on disinfecting surfaces in order to limit transmission through fomites (Van Doremalen et al., 2020) and maintaining social distance in order to limit transmission via large drops generated by coughs and sneezes (Bourouiba et al., 2014). The efficacy of these measures has been increasingly called into question, however, since there is scant evidence for fomite transmission (Lewis, 2021) and large-drop transmission is effectively eliminated by masks (Moghadas et al., 2020).

There is now overwhelming evidence that the pathogen responsible for COVID-19, severe-acute-respiratory-syndrome coronavirus 2 (SARS-CoV-2), is transmitted primarily through exhaled aerosol droplets suspended in indoor air (Prather et al., 2020; Morawska and Milton, 2020; Morawska and Cao, 2020; Jayaweera et al., 2020; Zhang et al., 2020b; Bazant and Bush, 2021). Notably, airborne transmission provides the only rational explanation for the so-called “super-spreading events”, which have now been well chronicled and all took place indoors (Miller et al., 2020; Moriarty, 2020; Hamner, 2020; Shen et al., 2020; Nishiura et al., 2020; Kwon et al., 2020; Hwang et al., 2020). The dominance of indoor airborne transmission is further supported by the fact that face-mask directives have been more effective in limiting the spread of COVID-19 than either social distancing directives or lockdowns (Zhang et al., 2020b; Stutt et al., 2020). Indeed, a recent analysis of spreading data from Massachusetts public schools where masking was strictly enforced found no statistically significant effect of social distance restrictions that ranged from 3 feet to 6 feet (van den Berg et al., 2021). Finally, the detection of infectious SARS-CoV-2 virions suspended in hospital room air as far as 18 feet from an infected patient provides direct evidence for the viability of airborne transmission of COVID-19 (Lednicky et al., 2020; Santarpia et al., 2020).

With a view to informing public health policy, we proceed by developing a quantitative approach to mitigating the indoor airborne transmission of COVID-19, an approach that might be similarly applied to other airborne respiratory diseases. The canonical theoretical framework of Wells (1955) and Riley et al. (1978) describes airborne transmission in an indoor space that is well-mixed by ambient air flows, so that infectious aerosols are uniformly dispersed throughout the space (Gammaitoni and Nucci, 1997; Beggs et al., 2003; Nicas et al., 2005; Noakes et al., 2006; Stilianakis and Drossinos, 2010). While exceptions to the well-mixed room are known to arise (Bhagat et al., 2020), supporting evidence for the well-mixed approximation may be found in both theoretical arguments (Bazant and Bush, 2021) and computer simulations of natural and forced convection (Foster and Kinzel, 2021). The Wells-Riley model and its extensions have been applied to a number of super-spreading events and used to assess the risk of COVID-19 transmission in a variety of indoor settings (Miller et al., 2020; Buonanno et al., 2020b,a; Prentiss et al., 2020; Evans, 2020).

A safety guideline for mitigating indoor airborne transmission of COVID-19 has recently been derived that indicates an upper bound on the cumulative exposure time, that is, the product of the number of occupants and the exposure time (Bazant and Bush, 2021). This bound may be simply expressed in terms of the relevant variables, including the room dimensions, ventilation, air filtration, mask efficiency and respiratory activity. The guideline has been calibrated for COVID-19 using epidemiological data from the best characterized super-spreading events and incorporates the measured dependence of expiratory droplet-size distributions on respiratory and vocal activity (Morawska et al., 2009; Asadi et al., 2019, 2020a). An online app has facilitated its widespread use during the pandemic (Khan et al., 2020). The authors also considered the additional risk of turbulent respiratory plumes and jets (Abkarian et al., 2020a,b), as need be considered when masks are not worn. The accuracy of the guideline is necessarily limited by uncertainties in a number of model parameters, which will presumably be reduced as more data is analyzed from indoor spreading events.

Carbon dioxide measurements have been used for decades to quantify airflow and zonal mixing in buildings and so guide the design of heating, ventilation, and air-conditioning (HVAC) systems (Fisk and De Almeida, 1998; Seppänen et al., 1999). Such measurements thus represent a natural source of data for assessing indoor air quality, especially as they rely only on relatively inexpensive, widely available CO2 sensors. Quite generally, high carbon dioxide levels in indoor settings are known to be associated with poor health (Salisbury, 1986; Hung and Derossis, 1989; Seppänen et al., 1999). Statistically significant correlations between CO2 levels and illness-related absenteeism in both the work place (Milton et al., 2000) and classrooms (Mendell et al., 2013; Shendell et al., 2004) have been widely reported (Li et al., 2007). Direct correlations between CO2 levels and concentration of airborne bacteria have been found in schools (Liu et al., 2000). Correlations between outdoor air exchange rates and respiratory infections in dorm rooms have also been reported (Sun et al., 2011; Bueno de Mesquita et al., 2020). Despite the overwhelming evidence of such correlations and the numerous economic analyses that underscore their negative societal impacts (Milton et al., 2000; Fisk, 2000), using CO2 monitors to make quantitative assessments of the risk of indoor disease transmission is a relatively recent notion (Li et al., 2007).

Rudnick and Milton (2003) first proposed the use of Wells-Riley models, in conjunction with measurements of CO2 concentration, to assess airborne transmission risk indoors. Their model treats CO2 concentration as a proxy for infectious aerosols: the two were assumed to be produced proportionally by the exhalation of an infected individual and removed at the same rate by ventilation. The current pandemic has generated considerable interest in using CO2 monitoring as a tool for risk management of COVID-19 (Bhagat et al., 2020; Hartmann and Kriegel, 2020). The Rudnick-Milton model has recently been extended by Peng and Jimenez (2021) through consideration of the different removal rates of CO2 and airborne pathogen. They conclude by predicting safe CO2 levels for COVID-19 transmission in various indoor spaces, which vary by up to two orders of magnitude.

We here develop a safety guideline for limiting indoor airborne transmission of COVID-19 by expressing the safety guideline of Bazant and Bush (2021) in terms of CO2 concentration. Doing so makes clear that one must limit not only the CO2 concentration, but also the occupancy time. Our model accounts for the effects of pathogen filtration, sedimentation and deactivation in addition to the variable aerosol production rates associated with different respiratory and vocal activities, all of which alter the relative concentrations of airborne pathogen and CO2. Our guideline thus quantifies the extent to which safety limits may be extended by mitigation strategies such as mask directives, air filtration and the imposition of ‘quiet spaces’.

In §2, we rephrase the indoor safety guideline of Bazant and Bush (2021) in terms of the room’s carbon dioxide concentration. In §3, we present theoretical descriptions of the evolution of CO2 concentration and infectious aerosol concentration in an indoor space, and highlight the different physical processes influencing the two. We then model the disease transmission dynamics, which allows for the risk of indoor airborne transmission to be assessed from CO2 measurements taken in real time. In §4, we apply our model to a pair of data sets tracking the evolution of CO2 concentration in specific office and classroom settings. These examples illustrate how CO2 monitoring, when coupled with our safety guideline, provides a means of assessing and mitigating the risk of indoor airborne transmission of respiratory pathogens.

2. Safety Guideline for the Time-Averaged Carbon Dioxide Concentration

2.1. Occupancy-based safety guideline

We begin by recalling the safety guideline Bazant and Bush (2021) for limiting indoor airborne disease transmission in a well-mixed space. The guideline would impose an upper bound on the cumulative exposure time: Embedded Image where Nt is the number of possible transmissions (pairs of infected and susceptible persons) and τ is the time in the presence of the infected person(s). The reader is referred to Table 1 for a glossary of symbols and their characteristic values. Qb is the mean breathing flow rate and V the room volume. The risk tolerance ϵ < 1 is the prescribed bound on the probability of at least one transmission, as should be chosen judiciously according to the vulnerability of the population (Garg, 2020); for example, Bazant and Bush (2021) suggested ϵ = 10% for children and 1% for the elderly.

View this table:
Table 1.

Glossary of symbols arising in our theory, their units and characteristics values.

The only epidemiological parameter in the guideline, Cq, is the infectiousness of exhaled air, measured in units of infection quanta per volume for a given aerosolized pathogen. The notion of “infection quantum” introduced by Wells (1955) is widely used in epidemiology to measure the expected rate of disease transmission, which may be seen as a transfer of infection quanta between pairs of infected and susceptible individuals. For airborne transmission, a suitable concentration of infection quanta per volume, Cq, can thus be associated with exhaled air without reference to the microscopic pathogen concentration. Notably, Cq is known to depend on the type of respiratory and vocal activity (resting, exercising, speaking, singing, etc.), being larger for the more vigorous activities (Buonanno et al., 2020b; Bazant and Bush, 2021) The relative susceptibility sr is introduced as a scaling factor for Cq that accounts for differences in the transmissibility of different respiratory pathogens, such as bacteria or viruses (Rudnick and Milton, 2003; Li et al., 2008) with different strains (Volz et al., 2021; Davies et al., 2020), and for differences in the susceptibility of different populations, such as children and adults (Riediker and Morawska, 2020; Zhang et al., 2020a; Zhu et al., 2020).

The mask penetration probability, pm (r), is a function of drop size that is bounded below by 0 (for the ideal limit of perfect mask filtration) and above by 1 (appropriate when no mask is worn). Standard surgical masks at low flow rates allow only 0.04-1.5% of the most infectious (submicron) aerosols to penetrate (Chen and Willeke, 1992), values that should be increased by a factor of 2 − 10 to account for imperfect fit (Oberg and Brosseau, 2008). Cloth masks show much greater variability (Konda et al., 2020). The mask penetration probability may also depend on respiratory activity (Asadi et al., 2020b) and direction of airflow (Pan et al., 2020). Here, for the sake of simplicity, we treat it as being constant over the limited aerosol size range of interest, and evaluate Embedded Image at the effective aerosol radius Embedded Image to be defined below. Bazant and Bush (2021) thus suggested Embedded Image for surgical masks (Li et al., 2008; Oberg and Brosseau, 2008), Embedded Image for hybrid cloth face coverings and Embedded Image for single-layer fabrics (Konda et al., 2020). Notably, even low quality masks can significantly reduce transmission risk since the bound on cumulative exposure time, Eq. (1), scales as Embedded Image.

Finally, we define Embedded Image as an effective relaxation rate of the infectious aerosol-borne pathogen concentration, C (r, t), evaluated at the effective aerosol radius Embedded Image. The size-dependent relaxation rate of the droplet-borne pathogen has four distinct contributions, Embedded Image where λa is the ventilation rate (outdoor air exchanges per time). and λ f (r) = pf (r) λr is the filtration rate, where pf (r) is the droplet removal efficiency for air filtration at a rate λr (recirculated air changes per time). λs (r) = vs (r) A/V is the net sedimentation rate for infectious droplets with the Stokes settling velocity vs (r) sedimenting through a well-mixed ambient to a floor of area A (Corner and Pendlebury, 1951; Martin and Nokes, 1988). Finally, λv (r) is the deactivation rate of the aerosolized pathogen, which depends weakly on humidity and droplet size (Yang and Marr, 2011; Lin and Marr, 2019; Marr et al., 2019), and may be enhanced by other factors such as ultraviolet (UV-C) irradiation (Hitchman, 2021; García de Abajo et al., 2020), chemical disinfectants (Schwartz et al., 2020), or cold plasma release (FilipiΔ et al., 2020; Lai et al., 2016).

Notably, only the first of the four removal rates apparent in Eq. (2) is relevant in the evolution of CO2; thus, the concentrations of CO2 and airborne pathogen are not strictly slaved to one another. Specifically, the proportionality between the two equilibrium concentrations varies in different indoor settings (Peng and Jimenez, 2021), for example in response to room filtration (Hartmann and Kriegel, 2020). Moreover, when transient effects arise, for example, following the arrival of an infectious individual or the opening of a window, the two concentrations adjust at different rates. Finally, we note that there may also be sources of CO2 other than human respiration, such as emissions from animals, stoves, furnaces, fireplaces, or carbonated beverages, as well as sinks of CO2, such as plants, construction materials or pools of water, which we neglect for simplicity. As such, following Rudnick and Milton (2003), we assume that the primary source of excess CO2 is exhalation by the human occupants of the indoor space.

In order to prevent the growth of an epidemic, the safety guideline should bound the indoor reproductive number, ℛin, which is the expected number of transmissions if an infectious person enters a room full of susceptible persons. Indeed, the safety guideline, Eq. (1), corresponds to the bound ℛin < ϵ with the choice Nt = N − 1, and so would limit the risk of an infected person entering the room of occupancy N transmitting to any other during the exposure time τ. If the epidemic is well underway or subsiding, the guideline should take into account the prevalence of infection pi and immunity pim (as achieved by previous exposure or vaccination) in the local population. Assuming a trinomial distribution of N persons who are infected, immune or susceptible, with mutually exclusive probabilities pi, pim and ps = 1 − pi −pim, respectively, the expected number of infected-susceptible pairs is N (N −1) pips. It is natural to switch between these two limits (Nt = N − 1 and Nt = N (N − 1) pips) when one infected person is expected to be in the room, Npi = 1, and thus set Embedded Image One may thus account for the changing infection prevalence pi and increasing immunity pim in the local population as the pandemic evolves.

2.2. CO2-based safety guideline

The total rate of CO2 production by respiration in the room is given by p2 = NQbC2,b, where C2,b is the CO2 concentration of exhaled air, approximately C2,b = 38, 000 ppm, although the net CO2 production rate, QbC2,b varies considerably with body mass and physical activity (Persily and de Jonge, 2017). If the production rate p2 and the ventilation flow rate Q = λaV are constant, then the steady-state value of the excess CO2 concentration, relative to the steady background concentration C0 at zero occupancy, is given by Embedded Image which is simply the ratio of the individual CO2 flow rate, QbC2,b, to the ventilation flow rate per person, Q/N. We note that the outdoor CO2 concentration is typically in the range C0 = 250 − 450 ppm, with higher values in urban environments (Prill et al., 2000). In the absence of other indoor CO2 sources, human occupancy in poorly ventilated spaces can easily lead to CO2 levels of several thousand ppm. People have reported headaches, slight nausea, drowsiness, and decreased decision-making performance for levels above 1000 ppm (Fisk et al., 2013; Krawczyk et al., 2016), while short exposures to much higher levels may go unnoticed. As an example of CO2 limits in industry, the American Conference of Governmental Industrial Hygienists recommends a limit of 5000 ppm for an 8-hour period and 30,000 ppm for 10 minutes. A value of 40, 000 ppm is considered to be immediately life-threatening.

The safety guideline, Eq. (1), was derived on the basis of the conservative assumption that the infectious aerosol concentration has reached its maximum, steady-state value. If we assume, for consistency, that the CO2 concentration has done likewise, and so approached the value expressed in Eq. (4), then the guideline can be recast as a bound on the safe mean excess CO2 concentration, Embedded Image where we replace the steady excess CO2 concentration with its time average, ⟨C2 ⟩ ≈ C2,s, and define the mean quanta emission rate per infected person, λq = QbCq. For the early to middle stages of an epidemic or when pi and pim are not known, we recommend setting N/Nt = 1 < N/ (N − 1) ≈ 1, for a conservative CO2 bound that limits the indoor reproductive number. In the later stages of an epidemic, as the population approaches herd immunity (pi →0, pim →1), the safe CO2 bound diverges, N/Nt → ∞, and so may be supplanted by the limits on carbon dioxide toxicity noted above, that lie in the range 5, 000 − 30, 000 ppm for 8-hour and 10-minute exposures.

Our simple CO2-based safety guideline, Eq. (5), reveals scaling laws for exposure time, filtration, mask use, infection prevalence and immunity, factors that are not accounted for by directives that would simply impose a limit on CO2 concentration. The substantial increase in safe occupancy times, as one proceeds from the peak to the late stages of the pandemic, is evident in the difference between the solid and dashed lines in Figure 1, which were evaluated for the case of a typical classroom in the United States (Bazant and Bush, 2021). This example shows the critical role of exposure time in determining the safe CO2 level, a limit that can be increased dramatically by mask use and to a lesser extent by filtration. When infection prevalence pi falls below 10 per 100,000 (an arbitrarily chosen small value), the chance of transmission is extremely low, allowing for long occupancy times. The risk of transmission at higher levels of prevalence, as may be deduced by interpolating between the solid and dashed lines in Figure 1, could also be rationally managed by monitoring the CO2 concentration and adhering to the guideline.

Figure 1.
Figure 1.

Illustration of the safety guideline, Eq. (5), which bounds the safe excess CO2 (ppm) and exposure time g (hours). Here, we consider the case of a classroom with N = 25 occupants, assumed to be children engaging in normal speech and light activity (λq sr = 30 quanta/h) with moderate risk tolerance (ϵ = 10%). Compared to the most restrictive bound on the indoor reproductive number without any precautions (red line), the safe CO2 level or occupancy time is increased by at least an order of magnitude by the use of face masks (blue line), even with relatively inconsistent use of cloth masks (pm = 30%). The effect of air filtration (green line) is relatively small, shown here for a case of efficient HEPA filtration (pf = 99%) with 17% outdoor air fraction (λf = 5λa). All three bounds are increased by several orders of magnitude (dashed lines) during late pandemic conditions (pi ps = 10 per 100, 000), when it becomes increasingly unlikely to find an infected-susceptible pair in the room. The other parameters satisfy (λ ν + λs)/ λa = 0.5, as could correspond to, for example, λ ν= 0.3/h, λB = 0.2/h and λa = 1/h (1 ACH).

3. Mathematical Model of CO2 Monitoring to Predict Airborne Disease Transmission Risk

3.1. CO2 dynamics

We follow the traditional approach of modeling gas dynamics in a well mixed room (Shair and Heitner, 1974), as a continuous stirred tank reactor (Davis and Davis, 2012). Given the time dependence of occupancy, N(t), mean breathing flow rate, Qb (t), and ventilation flow rate, Qa (t) = λa (t)V, one may express the evolution of the excess CO2 concentration C2(t) in a well-mixed room through Embedded Image where Embedded Image is the exhaled CO2 production rate. The relaxation rate for excess CO2 in response to changes in P2(t) is precisely equal to the ventilation rate, λa (t) = Qa (t)/V. For constant λa, the general solution of Eq. (6) for C2(0) = 0 is given by Embedded Image which can be derived by Laplace transform or using an integrating factor. The time-averaged excess CO2 concentration can be expressed as Embedded Image by switching the order of time integration. If P2 (t) is slowly varying over the ventilation time scale Embedded Image, the time-averaged CO2 concentration may be approximated as Embedded Image where the excess CO2 concentration approaches the ratio of the mean exhaled CO2 production rate to the ventilation flow rate at long times, Embedded Image, as indicated in Eq. (4).

3.2. Infectious aerosol dynamics

Following Bazant and Bush (2021), we assume that the radius-resolved concentration of infectious aerosol-borne pathogen, C(r, t), evolves according to Embedded Image where the mean production rate, Embedded Image depends on the number of infected persons in the room, I (t), and the size distribution nd (r, t) of exhaled droplets of volume Vd (r) containing pathogen (i.e. virions) at microscopic concentration, cν (r). The droplet size distribution is known to depend on expiratory and vocal activity (Morawska et al., 2009; Asadi et al., 2019, 2020c). Quite generally, the aerosols evolve according to a dynamic sorting process (Bazant and Bush, 2021): the drop-size distribution evolves with time until an equilibrium distribution is obtained.

Given the time evolution of excess CO2 concentration, C2(t), one may deduce the radius-resolved pathogen concentration C(r, t) by integrating the coupled differential equation, Embedded Image This integration can be done numerically or analytically via Laplace transform or integrating factors if one assumes that λa, λc, P and P2 all vary slowly over the ventilation (air change) time scale, Embedded Image. In that case, the general solution takes the form Embedded Image where we consider the infectious aerosol build-up from C(r, 0) = 0.

3.3. Disease transmission dynamics

According to Markov’s inequality, the probability of at least one transmission taking place during the exposure time τ is bounded above by the expected number of airborne transmissions, Ta (τ), and the two become equal in the (typical) limit of rare transmissions, Ta (τ) ≪ 1. The expected number of transmissions to S(t) susceptible persons is obtained by integrating the mask-filtered inhalation rate of infection quanta over both droplet radius and time, Embedded Image where Ci (r) is the infectivity of the aerosolized pathogen. The infectivity is measured in units of infection quanta per pathogen and generally depends on droplet size. One might expect pathogens contained in smaller aerosol droplets with r < 5 μm to be more infectious than those in larger drops, as reported by Santarpia et al. (2020) for SARS-CoV-2, on the grounds that smaller drops more easily traverse the respiratory tract, absorb and coalesce onto exposed tissues, and allow pathogens to escape more quickly by diffusion to infect target cells. The mask penetration probability Pm(r) also decreases rapidly with increasing drop size above the aerosol range for most filtration materials (Chen and Willeke, 1992; Oberg and Brosseau, 2008; Konda et al., 2020; Li et al., 2008), so the integration over radius in Eq. (15) gives the most weight to the aerosol size range, roughly r < 5μm, which also coincides with the maxima in exhaled droplet size distributions (Morawska et al., 2009; Asadi et al., 2019, 2020c).

The inverse of the infectivity, Embedded Image, is equal to the “infectious dose” of pathogens from inhaled aerosol droplets that would cause infection with probability 1−(1/e) = 63%. Bazant and Bush (2021) estimated the infectious dose for SARS-CoV-2 to be on the order of ten aerosol-borne virions. Notably, the corresponding infectivity, Ci ∼0.1, is an order of magnitude larger than previous estimates for SARS-CoV (Watanabe et al., 2010; Buonanno et al., 2020b), which is consistent with only COVID-19 reaching pandemic status. The infectivity is known to vary across different age groups and pathogen strains, a variability that is captured by the relative susceptibility, sr. For example, Bazant and Bush (2021) suggest assigning sr = 1 for the elderly (over 65 years old), sr = 0.68 for adults (aged 15-64) and sr = 0.23 for children (aged 0-14) for the original Wuhan strain of SARS-CoV-2, based on a study of transmission in quarantined households in China (Zhang et al., 2020a). The authors further suggested multiplying these values by 1.6 for the more infectious variant of concern of the lineage B.1.1.7 (VOC 202012/01), which recently emerged in the United Kingdom with a reproductive number that was 60% larger than that of the original strain (Volz et al., 2021; Davies et al., 2020).

3.4. Approximate Formula for the Airborne Transmission Risk from CO2 Measurements

Equations (14) and (15) provide an approximate solution to the full model that depends on the exhaled droplet size distribution,nd(r,t), and mean breathing rate, Qb (t), of the population in the room. Since the droplet distributions nd(r) A have only been characterized in certain idealized experimental conditions (Morawska et al., 2009; Asadi et al., 2020a,c), it is useful to integrate over r to obtain a simpler model that can be directly calibrated for different modes of respiration using epidemiological data (Bazant and Bush 2021). Assuming Qb(t), I(t), S(t), N(t) and nd(r,t)vary slowly over the relaxation time,Embedded Image we may substitute Eq. (14) into Eq. (15) and perform the time integral of the second term to obtain Embedded Image where nq(r,t)= nd(r,t)Vd(r)cv(r)ci(r) is the radius-resolved exhaled quanta concentration.

Following Bazant and Bush (2021), we define an effective radius of infectious aerosols Embedded Image such that Embedded Image where Embedded Image is the exhaled quanta concentration, which may vary in time with changes in expiratory activity, for example, following a transition from nose breathing to speaking. In principle, the effective radius Embedded Image can be evaluated, given a complete knowledge of the dependence on drop radius of the mask penetration probability, Pm(r), and of all the factors that determine the exhaled quanta concentration, nq(r,t)and pathogen removal rate, λc(r). While these dependencies are not readily characterized, typical values of Embedded Image are at the scale of several microns, based on the size dependencies of nd(r,t), ci(r),and pm(r) noted above.

Further simplifications allow us to derive a formula relating CO2 measurements to transmission risk. By assuming that Cq(t) varies slowly over the timescale of concentration relaxation, one may approximate the memory integral with the same effective radius Embedded Image. Thus, accounting for immunity and infection prevalence in the population via Embedded Image we obtain a formula for the expected number of airborne transmissions, Embedded Image in terms of the excess CO2 time series, C2(t), where λq(t) = Cq(t)Qb(t)is the mean quanta emission rate. It is also useful to define the expected transmission rate, Embedded Image which allows for direct assessment of airborne transmission risk based on CO2 levels. A pair of examples of such assessments will be presented in §4.

Notably, the mean airborne transmission rate per expected infected-susceptible pair, βa(t), reflects the environment’s memory of the recent past, which persists over the pathogen relaxation time scale,Embedded Image The CO2 concentration in Eq. (9) has a longer memory of past changes in CO2 sources or ventilation, which persists over the air change time scale, Embedded Image since CO2 is unaffected by the filtration, sedimentation and deactivation rates evident in Eq. (2). The time delays between the production of CO2 and infectious aerosols by exhalation and their buildup in the well mixed air of a room shows that CO2 variation and airborne transmission are inherently non-Markovian stochastic processes. As such, any attempt to predict fluctuations in airborne transmission risk would require stochastic generalizations of the differential equations governing the mean variables, Eqs. (6) and (11), and so represent a stochastic formulation of the Wells-Riley model (Noakes and Sleigh, 2009).

3.5. Reduction to the CO2-based Safety Guideline

Finally, we connect the general result, Eq. (19), with the CO2 based safety guideline derived above, Eq. (5). Since C2(t) varies on the ventilation time scale, Embedded Image, which is necessarily longer than the relaxation time scale of the infectious aerosols, Embedded Image, we may assume that λq(t)2C2(t) is slowly varying and evaluate the integral in Eq. (19). We thus arrive at the approximation Embedded Image Since λq(t) C2(t) is slowly varying, the second term in brackets is negligible relative to the first for times longer than the ventilation time, Embedded Image. In this limit, the imposed bound on expected transmissions, Ta (τ) < ϵ, is approximated by Embedded Image This formula reduces to the safety guideline, Eq. (5), in the limit of constant mean quanta emission rate, λq, which confirms the consistency of our assumptions.

4. Examples

We proceed by illustrating the process by which the guideline, Eq. (5), can be coupled to real data obtained from CO2 monitors. Specifically, we consider time series of CO2 concentration gathered in classroom and in office settings at the Massachusetts Institute of Technology using an Atlas Scientific EZO-CO2 Embedded NDIR CO2 Sensor controlled with an Arduino Uno, and an Aranet4, respectively. Social distancing guidelines were adhered to, and masks were worn by all participants. We assume a constant exhaled CO2 concentration of 38, 000 ppm, and use the global minimum of the CO2 series as the background CO2 level C0 from which the excess concentration C2 (t) was deduced. Notably, the relatively small fluctuations in the CO2 measurements recorded in a variety of settings support the notion of a well-mixed room.

From Equations (19) and (20), we calculate the expected number of transmissions and transmission rate, assuming that there is one infected person in the room (Nt = N − 1). In this case, the expected number of transmissions is equal to the indoor reproductive number, Ta (τ) = ℛin(τ), and the transmission rate is Embedded Image. The approximations in the derivation are valid in these examples, so a direct numerical solution of Eq. (13) would yield indistinguishable results. In particular, the “slowly varying” assumptions are satisfied, since we keep N and I constant, and any time-dependence of Qb cancels in the ratio P /P2. The droplet distributions nd (r,t) A, C may vary in time, but no significant changes in mean respiratory activity were observed or measured. Moreover, accounting for any variations in λa and λc would require additional measurements in the same space, so these parameters may also be taken as constants during the relatively short exposure times considered.

We choose realistic values of the parameters that fall within the typical ranges estimated by Bazant and Bush (2021). The mean breathing rate is set to Qb = 0.5 m3 / h for light activity, and Embedded Image for sedimentation with Embedded Image. The viral deactivation rate for SARS-CoV-2 is set to λv = 0.3 / h, an estimate appropriate for 50% relative humidity. We plot Nt βa (τ) and ℛin= (τ) for cases where masks are and are not worn, and choose a mask penetration probability of pm ≈ 0.3, as is roughly appropriate for a cloth mask. The parameters Embedded Image, H and N are chosen according to the specific scenario presented. Exhaled COVID-19 quanta concentrations Cq for various expiratory activities are estimated from Fig. 2 of Bazant and Bush (2021). In order to be conservative, we assume that sr = 1, suitable for the high-risk individuals exposed to the Wuhan strain of SARS-CoV-2.

Figure 2.
Figure 2.

Measured CO2 concentration and calculated transmission rate in a two-person office. a) Black dots represent the concentration of CO2. The solid blue, dashed magenta and dash-dot green curves represents the transmission rate, as calculated from (20) for three different scenarios, two of which were hypothetical: (blue) the pair are not wearing masks and there is no filtration present; (magenta) the pair are not wearing masks and there is filtration present; (green) the pair are wearing masks and there is no filtration present. The orange solid curve denotes the period of exponential relaxation following the exit of the room’s occupants, from which one may infer the room’s ventilation rate, λa = 2.3 /h. b) Corresponding blue, magenta and green curves, deduced by integrating (20), indicate the total risk of transmission over the time of shared occupancy. If the pair were not wearing masks, the safety limit ℛin < 0.1 would be violated after approximately an hour.

Figure 3.
Figure 3.

Measured CO2 concentration and calculated transmission rate for 12 masked students in a university lecture hall. a) Black dots represent the concentration of CO2. The solid blue and dash-dot green curves represent the transmission rate, as calculated from (20), when the occupants are wearing masks, and in the hypothetical case where they are not, respectively. b) The solid blue and dash-dot green curves indicate the total risk of transmission with and without masks, respectively, as deduced by integrating (20) over time. Even had masks not been worn, the safety guideline would not have been violated during the lecture.

4.1. Small office with two workers

Figure 2(a) shows CO2 measurements taken in an office of length L = 4.2 m, width, W= 3 m, and height H = 3 m. Initially, a single worker is present, but at 19:00, a second worker arrives at the office. The workers exit the office at 21:09, return at 21:39, and exit again at 22:30. As the participants were speaking, we use Cq= 72 quanta /m3. We compute the room’s ventilation rate, λa 2.3 / h, from the exponential relaxation that follows the occupants’ exit from the office (as indicated by the orange curve in Figure 2). The office was equipped with moderate ventilation, which we characterize with pf = 0.99 and λr = 6 / h.

As shown in Figure 2(b), if the office mates were not wearing masks, the safety guideline of expected transmissions Ta < 10% would be violated after approximately an hour together. However, office mates wearing cloth masks remain well under the safety guideline during the 2.5 hours spent together. Filtration without masks also extends the occupancy time limits, however ℛin approaches the safety limit after approximately 2 hours. While this example should not be taken as a definitive statement of danger or safety in this setting, it does serves to illustrate how our CO2 guideline can be implemented in a real-life situation.

4.2. University classroom adhering to social distancing guidelines

We next monitor CO2 levels during a university lecture. There were N= 12 participants in a lecture hall of length L = 13 m, width, W= 12 m and height H = 3 m. The lecture started at 13:05 and finished at 13:50. Four people remained in the room for 30 minutes after class. The classroom has mechanical ventilation, which we characterize in terms of λa = 2 / h (30 minute outdoor air change). During the lecture, the professor spoke while the students were quiet and sedentary; thus, we assume Cq = 30 quanta /m3. The maximum CO2 concentration reached in the classroom was ≈550 ppm, the excess level no more than 100 ppm. Thus, the safety limit was never exceeded, and would not have been even if masks had not been worn. We note that this lecture hall was particularly large, well-ventilated and sparsely populated, and so should not be taken as being representative of a classroom setting.

A more complete assessment of safety in schools would require integrating CO2 concentrations over a considerably longer time interval. For example, if students are tested weekly for COVID-19, then one should assess the mean CO2 concentration in class during the course of an entire school week. Nevertheless, this second example further illustrates the manner in which our model may be applied to real-world settings and suggests that precautions such as ventilation, filtration and mask use, can substantially increase safe occupancy times.

5. Conclusion

Mounting evidence suggests that COVID-19 is spread primarily via indoor airborne transmission. Such an inference is no surprise, as such is also the case for many other respiratory illnesses, including influenza, tuberculosis, measles and severe-acute respiratory syndrome (as is caused by a precursor to SARS-CoV-2, the coronavirus SARS-CoV). More than a year into the pandemic, public health guidance continues to emphasize the importance of social distancing and surface cleaning, despite evidence that mask directives are much more effective than either in limiting airborne transmission. We have here illustrated the manner in which CO2 monitoring may be used in conjunction with the safety guideline of Bazant and Bush (2021) in assessing the risk of indoor airborne respiratory disease transmission, including that of COVID-19. We thus hope to inform personal and policy decisions about closing and re-opening indoor spaces, such as schools and businesses.

We have here reformulated the COVID-19 indoor safety guideline of Bazant and Bush (2021), expressing it in terms of cumulative exposure to carbon dioxide, which can be readily monitored in real time for most indoor spaces. In so doing, we have built upon the important early work of Rudnick and Milton (2003), as was recently extended and applied to COVID-19 by Peng and Jimenez (2021). The guideline of Bazant and Bush (2021) makes clear that, since the risk of indoor airborne infection is determined by the total volume of pathogen inhaled, safety limits intended to protect against it must be expressed in terms of occupancy time. Likewise, in the context of CO2 measurements, safety limits cannot be expressed solely in terms of a limit on carbon dioxide levels, but must also depend on occupancy time. As we have demonstrated with our two case studies, the safety guideline (5), when coupled with CO2 monitors, allows for real-time assessment of risk of airborne disease transmission ℛin in indoor spaces. Moreover, this approach has the distinct advantage that one can assess certain key model parameters, including the background concentration of CO2 and the room’s ventilation rate, directly from the CO2 measurements.

Within a well-mixed space, carbon dioxide is effectively a passive scalar that tracks the ambient flow, and is removed only through exchange with outdoor air. Aerosol-borne pathogen is subject to additional removal mechanisms, including filtration (by face masks and internal circulation), sedimentation and deactivation. Thus, the concentration of CO2 cannot be taken as a proxy for that of aerosol-borne pathogen without resolving the proportionality constant between the two that results from these additional removal mechanisms. We stress that the effects of face mask use are dramatic in reducing the ratio of aerosol-borne pathogen to CO2 concentration, and so in reducing the risk of indoor transmission. Finally, we note that the additional removal mechanisms acting on the droplet-borne pathogen alter not only its equilibrium concentrations relative to that of CO2, but their relaxation times in transient situations, as may be treated using the mathematical formalism of Bazant and Bush (2021).

We emphasize that the caveats enumerated by Bazant and Bush (2021) concerning the limitations of their safety guideline apply similarly here. First, there is considerable uncertainty in a number of model parameters, including the critical viral load and relative susceptibility. While our inferences are consistent with available data, we hope that these uncertainties will be reduced as more COVID-19 spreading events are characterized and analyzed. Second, when masks are not worn, there is substantial additional risk of short-range airborne transmission from respiratory jets and plumes, as accompany breathing, speaking (Abkarian et al., 2020a; Abkarian and Stone, 2020), coughing and sneezing (Bourouiba et al., 2014). Third, while the assumption of the well-mixed room is widely applied and represents a reasonable first approximation, it is known to have limitations (Bhagat et al., 2020; Linden et al., 1990; Linden, 1999). Notably, measured fluctuations in CO2 levels provide a direct means of assessing the validity of the well-mixed-room hypothesis, especially when several sensors are used simultaneously at different locations in the same space. Indeed, CO2 monitoring has been used for decades to assess the quality of air handling and zonal mixing in buildings (Seppänen et al., 1999; Fisk and De Almeida, 1998; Hung and Derossis, 1989; Salisbury, 1986; Cheng et al., 2011), and can now be re-purposed to assess the risk of indoor airborne disease transmission.

Our transmission theory and safety guideline provide a quantitative basis for the use of CO2 monitors in assessing the risk of indoor airborne disease transmission. Specifically, the simple guideline, Eq. (5), and mathematical formulae connecting CO2 data to the evolving transmission risk, Eqs. (19)-(20), pave the way for real-time assessment of personal risk in indoor spaces. Moreover, our results should allow for building-scale optimization of public health using CO2 sensors, wherein risk assessment might be weighed against the energy requirements of enhanced ventilation. Finally, our model provides a general framework for using CO2 monitors to mitigate the indoor airborne transmission of other respiratory illnesses, including the seasonal flu.

Extensions of our study would include implementing our guideline in spaces where not only CO2 is monitored, but also the room occupancy along with other relevant parameters appearing in our model. For example, monitoring decibel levels and type of vocalization could serve to inform the infectious aerosol production rate (Asadi et al., 2019, 2020c; Morawska et al., 2009). Changes in occupancy could also be monitored at entrances and exits (while maintaining anonymity). One could further envision feeding all such data into air regulation controls in order to ensure that our CO2-based indoor safety limit is never violated. We note that such a prospect would be most easily achieved in quasi-steady circumstances in which a fixed population is behaving in a predictable fashion over the course of an event of known duration, for example, for students in a lecture hall or passengers on a charter bus. In more complicated situations, our model provides a framework for optimizing sensor-based demand-controlled ventilation (Fisk and De Almeida, 1998) with a view to limiting transmission risk while reducing energy consumption and system costs.

In order to facilitate the application of our safety guideline, in the Supplementary Material we provide a link to an online app that computes the safety guideline in terms of both room occupancy and CO2 levels (Khan et al., 2020). The CO2-based guideline, available in the app’s Advanced Mode, may be used in conjunction with CO2 monitors to formulate safe re-opening policies for indoor spaces in the later stages of the pandemic.

Data Availability

The data presented in the article is available upon request.

https://indoor-covid-safety.herokuapp.com

Funding Statement

The authors received no funding for this work.

Declaration of Interests

The authors declare no conflict of interest.

Author Contributions

J.W.M.B. and M.Z.B. conceived the study. M.Z.B. developed the mathematical model. O.K. performed the experiments and analyzed the data with A.E.C. K.K. developed the online app with M.Z.B and J.W.M.B. All authors contributed to writing the manuscript.

Data Availability Statement

Raw data are available from M.Z.B. and O.K.

Ethical Standards

The research meets all ethical guidelines, including adherence to the legal requirements of the study country.

Supplementary Material

An online app is provided to compute the COVID-19 Indoor Safety Guideline, where the safe CO2 level versus time can be displayed in Advanced Mode (Khan et al., 2020).

Acknowledgements

We thank Laura Champion for valuable discussions and for sharing preliminary data.

Footnotes

  • Expanded discussion of the assumptions and results.

References

  1. ↵
    Abkarian, M., Mendez, S., Xue, N., Yang, F., and Stone, H. A. (2020a). Puff trains in speaking produce long-range turbulent jet-like transport potentially relevant to asymptomatic spreading of viruses. arXiv preprint arXiv:2006.10671.
  2. Abkarian, M., Mendez, S., Xue, N., Yang, F., and Stone, H. A. (2020b). Speech can produce jet-like transport relevant to asymptomatic spreading of virus. Proceedings of the National Academy of Sciences, 117(41):25237.
    OpenUrlAbstract/FREE Full Text
  3. ↵
    Abkarian, M. and Stone, H. (2020). Stretching and break-up of saliva filaments during speech: A route for pathogen aerosolization and its potential mitigation. Physical Review Fluids, 5(10):102301.
    OpenUrlCrossRef
  4. ↵
    Asadi, S., Bouvier, N., Wexler, A. S., and Ristenpart, W. D. (2020a). The coronavirus pandemic and aerosols: Does COVID-19 transmit via expiratory particles? Aerosol Science and Technology, 54(6):635–638.
    OpenUrlCrossRef
  5. ↵
    Asadi, S., Cappa, C. D., Barreda, S., Wexler, A. S., Bouvier, N. M., and Ristenpart, W. D. (2020b). Efficacy of masks and face coverings in controlling outward aerosol particle emission from expiratory activities. Scientific reports, 10(1):1–13.
    OpenUrl
  6. ↵
    Asadi, S., Wexler, A. S., Cappa, C. D., Barreda, S., Bouvier, N. M., and Ristenpart, W. D. (2019). Aerosol emission and superemission during human speech increase with voice loudness. Scientific Reports, 9(1):1–10.
    OpenUrl
  7. ↵
    Asadi, S., Wexler, A. S., Cappa, C. D., Barreda, S., Bouvier, N. M., and Ristenpart, W. D. (2020c). Effect of voicing and articulation manner on aerosol particle emission during human speech. PloS One, 15(1):e0227699.
    OpenUrlCrossRefPubMed
  8. ↵
    Bazant, M. Z. and Bush, J. W. M. (2021). A guideline to limit indoor airborne transmission of COVID-19. PNAS, in press. Preprint: Beyond six feet, medRxiv (2020), https://doi.org/10.1101/2020.08.26.20182824.
  9. ↵
    Beggs, C., Noakes, C., Sleigh, P., Fletcher, L., and Siddiqi, K. (2003). The transmission of tuberculosis in confined spaces: an analytical review of alternative epidemiological models. Int. J. Tuberculosis and Lung Disease, 7(11):1015–1026.
    OpenUrl
  10. ↵
    Bhagat, R. K., Wykes, M. D., Dalziel, S. B., and Linden, P. (2020). Effects of ventilation on the indoor spread of COVID-19. Journal of Fluid Mechanics, 903.
  11. ↵
    Bourouiba, L., Dehandschoewercker, E., and Bush, J. W. M. (2014). Violent expiratory events: on coughing and sneezing. Journal of Fluid Mechanics, 745:537–563.
    OpenUrl
  12. ↵
    Bueno de Mesquita, P. J., Noakes, C. J., and Milton, D. K. (2020). Quantitative aerobiologic analysis of an influenza human challenge-transmission trial. Indoor air, 30(6):1189–1198.
    OpenUrl
  13. Buonanno, G., Morawska, L., and Stabile, L. (2020a). Quantitative assessment of the risk of airborne transmission of SARS-CoV-2 infection. Environment International, 145:106112.
    OpenUrlPubMed
  14. ↵
    Buonanno, G., Stabile, L., and Morawska, L. (2020b). Estimation of airborne viral emission: quanta emission rate of SARS-CoV-2 for infection risk assessment. Environment International, 141:105794.
    OpenUrlCrossRefPubMed
  15. ↵
    Chen, C.-C. and Willeke, K. (1992). Aerosol penetration through surgical masks. American journal of infection control, 20(4):177–184.
    OpenUrlCrossRefPubMedWeb of Science
  16. ↵
    Chen, N., Zhou, M., Dong, X., Qu, J., Gong, F., Han, Y., Qiu, Y., Wang, J., Liu, Y., Wei, Y., et al. (2020). Epidemiological and clinical characteristics of 99 cases of 2019 novel coronavirus pneumonia in Wuhan, China: a descriptive study. The Lancet, 395(10223):507–513.
    OpenUrlCrossRef
  17. ↵
    Cheng, K.-C., Acevedo-Bolton, V., Jiang, R.-T., Klepeis, N. E., Ott, W. R., Fringer, O. B., and Hildemann, L. M. (2011). Modeling exposure close to air pollution sources in naturally ventilated residences: Association of turbulent diffusion coefficient with air change rate. Environmental Science & Technology, 45(9):4016–4022.
    OpenUrl
  18. ↵
    Corner, J. and Pendlebury, E. (1951). The coagulation and deposition of a stirred aerosol. Proceedings of the Physical Society. Section B, 64(8):645.
    OpenUrl
  19. ↵
    Davies, N. G., Barnard, R. C., Jarvis, C. I., Kucharski, A. J., Munday, J., Pearson, C. A., Russell, T. W., Tully, D. C., Abbott, S., Gimma, A., et al. (2020). Estimated transmissibility and severity of novel SARS-CoV-2 variant of concern 202012/01 in england. medRxiv.
  20. ↵
    Davis, M. E. and Davis, R. J. (2012). Fundamentals of Chemical Reaction Engineering. ourier Corporation.
  21. ↵
    Evans, M. (2020). Avoiding COVID-19: Aerosol guidelines. arXiv preprint arXiv:2005.10988.
  22. Filipić, A., Gutierrez-Aguirre, I., Primc, G., MozetiÃ, M., and Dobnik, D. (2020). Cold plasma, a new hope in the field of virus inactivation. Trends in Biotechnology, 38:1278–1291.
    OpenUrl
  23. ↵
    Fisk, W. J. (2000). Health and productivity gains from better indoor environments and their relationship with building energy efficiency. Annual review of energy and the environment, 25(1):537–566.
    OpenUrl
  24. ↵
    Fisk, W. J. and De Almeida, A. T. (1998). Sensor-based demand-controlled ventilation: a review. Energy and buildings, 29(1):35–45.
    OpenUrl
  25. ↵
    Fisk, W. J., Satish, U., Mendell, M. J., Hotchi, T., and Sullivan, D. (2013). Is co2 an indoor pollutant? higher levels of co2 may diminish decision making performance. ASHRAE JOURNAL, 55(LBNL-6148E).
  26. ↵
    Foster, A. and Kinzel, M. (2021). Estimating covid-19 exposure in a classroom setting: A comparison between mathematical and numerical models. Physics of Fluids, 33(2):021904.
    OpenUrl
  27. ↵
    Gammaitoni, L. and Nucci, M. C. (1997). Using a mathematical model to evaluate the efficacy of TB control measures. Emerging Infectious Diseases, 3(3):335.
    OpenUrlCrossRefPubMedWeb of Science
  28. ↵
    García de Abajo, F. J., Hernández, R. J., Kaminer, I., Meyerhans, A., Rosell-Llompart, J., and Sanchez-Elsner, T. (2020). Back to normal: An old physics route to reduce SARS-CoV-2 transmission in indoor spaces. ACS Nano, 14(7):7704–7713.
    OpenUrl
  29. ↵
    Garg, S. (2020). Hospitalization rates and characteristics of patients hospitalized with laboratory-confirmed coronavirus disease 2019: COVID-NET, 14 States, March 1–30, 2020. MMWR. Morbidity and Mortality Weekly Report, 69.
  30. ↵
    Hamner, L. (2020). High SARS-CoV-2 attack rate following exposure at a choir practice, Skagit County, Washington, March 2020. MMWR. Morbidity and Mortality Weekly Report, 69.
  31. ↵
    Hartmann, A. and Kriegel, M. (2020). Risk assessment of aerosols loaded with virus based on co2-concentration. TU Berlin preprint, accessed April 3, 2021.
  32. ↵
    Hitchman, M. L. (2021). A new perspective of the chemistry and kinetics of inactivation of covid-19 coronavirus aerosols. Future Virology, 15(12):823–835.
    OpenUrl
  33. ↵
    Hung, I.-F. and Derossis, P. (1989). Carbon dioxide concentration as indicator of indoor air quality. Journal of Environmental Science and Health. Part A: Environmental Science and Engineering, 24(4):379–388.
    OpenUrl
  34. ↵
    Hwang, S. E., Chang, J. H., Bumjo, O., and Heo, J. (2020). Possible aerosol transmission of COVID-19 associated with an outbreak in an apartment in Seoul, South Korea. International Journal of Infectious Diseases.
  35. ↵
    Jayaweera, M., Perera, H., Gunawardana, B., and Manatunge, J. (2020). Transmission of COVID-19 virus by droplets and aerosols. Environ Res., 188(109819).
  36. ↵
    Khan, K., Bazant, M. Z., and Bush, J. W. M. (2020). COVID-19 indoor safety guideline. Online app, https://indoor-covid-safety.herokuapp.com.
  37. ↵
    Konda, A., Prakash, A., Moss, G. A., Schmoldt, M., Grant, G. D., and Guha, S. (2020). Response to letters to the editor on aerosol filtration efficiency of common fabrics used in respiratory cloth masks: Revised and expanded results. ACS Nano, 14(9):10764–10770.
    OpenUrl
  38. ↵
    Krawczyk, D., Rodero, A., Gladyszewska-Fiedoruk, K., and Gajewski, A. (2016). Co2 concentration in naturally ventilated classrooms located in different climates—measurements and simulations. Energy and Buildings, 129:491–498.
    OpenUrl
  39. ↵
    Kwon, K.-S., Park, J.-I., Park, Y. J., Jung, D.-M., Ryu, K.-W., and Lee, J.-H. (2020). Evidence of long-distance droplet transmission of SARS-CoV-2 by direct air flow in a restaurant in Korea. J. Korean Med. Sci., 35((46):e415).
    OpenUrl
  40. ↵
    Lai, A., Cheung, A., Wong, M., and Li, W. (2016). Evaluation of cold plasma inactivation efficacy against different airborne bacteria in ventilation duct flow. Building and Environment, 98:39–46.
    OpenUrl
  41. ↵
    Lednicky, J. A., Lauzard, M., Fan, Z. H., Jutla, A., Tilly, T. B., Gangwar, M., Usmani, M., Shankar, S. N., Mohamed, K., Eiguren-Fernandez, A., et al. (2020). Viable SARS-CoV-2 in the air of a hospital room with COVID-19 patients. International Journal of Infectious Diseases, 100:476–482.
    OpenUrlCrossRef
  42. ↵
    Lewis, D. (2021). Covid-19 rarely spreads through surfaces. so why are we still deep cleaning. Nature, 590(7844):26–28.
    OpenUrlCrossRef
  43. ↵
    Li, Q., Guan, X., Wu, P., Wang, X., Zhou, L., Tong, Y., Ren, R., Leung, K. S., Lau, E. H., Wong, J. Y., et al. (2020). Early transmission dynamics in Wuhan, China, of novel coronavirus–infected pneumonia. New England Journal of Medicine, 382:1199–1207.
    OpenUrlCrossRefPubMed
  44. ↵
    Li, Y., Guo, Y. P., Wong, K. C. T., Chung, W. Y. J., Gohel, M. D. I., and Leung, H. M. P. (2008). Transmission of communicable respiratory infections and facemasks. Journal of multidisciplinary healthcare, 1:17.
    OpenUrlCrossRef
  45. ↵
    Li, Y., Leung, G. M., Tang, J., Yang, X., Chao, C., Lin, J. Z., Lu, J., Nielsen, P. V., Niu, J., Qian, H., et al. (2007). Role of ventilation in airborne transmission of infectious agents in the built environment-a multidisciplinary systematic review. Indoor air, 17(1):2–18.
    OpenUrlCrossRefPubMedWeb of Science
  46. ↵
    Lin, K. and Marr, L. C. (2019). Humidity-dependent decay of viruses, but not bacteria, in aerosols and droplets follows disinfection kinetics. Environmental Science & Technology, 54(2):1024–1032.
    OpenUrl
  47. ↵
    Linden, P., Lane-Serff, G., and Smeed, D. (1990). Emptying filling boxes: the fluid mechanics of natural ventilation. Journal of Fluid Mechanics, 212:309–335.
    OpenUrl
  48. ↵
    Linden, P. F. (1999). The fluid mechanics of natural ventilation. Annual Review of Fluid Mechanics, 31(1):201–238.
    OpenUrlCrossRefWeb of Science
  49. ↵
    Liu, L.-J. S., Krahmer, M., Fox, A., Feigley, C. E., Featherstone, A., Saraf, A., and Larsson, L. (2000). Investigation of the concentration of bacteria and their cell envelope components in indoor air in two elementary schools. Journal of the Air & Waste Management Association, 50(11):1957–1967.
    OpenUrl
  50. ↵
    Marr, L. C., Tang, J. W., Van Mullekom, J., and Lakdawala, S. S. (2019). Mechanistic insights into the effect of humidity on airborne influenza virus survival, transmission and incidence. Journal of the Royal Society Interface, 16(150):20180298.
    OpenUrl
  51. ↵
    Martin, D. and Nokes, R. (1988). Crystal settling in a vigorously converting magma chamber. Nature, 332(6164):534–536.
    OpenUrlCrossRefGeoRef
  52. ↵
    Mendell, M. J., Eliseeva, E. A., Davies, M. M., Spears, M., Lobscheid, A., Fisk, W. J., and Apte, M. G. (2013). Association of classroom ventilation with reduced illness absence: a prospective study in california elementary schools. Indoor air, 23(6):515–528.
    OpenUrlCrossRefPubMedWeb of Science
  53. ↵
    Miller, S. L., Nazaroff, W. W., Jimenez, J. L., Boerstra, A., Buonanno, G., Dancer, S. J., Kurnitski, J., Marr, L. C., Morawska, L., and Noakes, C. (2020). Transmission of SARS-CoV-2 by inhalation of respiratory aerosol in the Skagit Valley Chorale superspreading event. Indoor Air, page in press.
  54. ↵
    Milton, D. K., Glencross, P. M., and Walters, M. D. (2000). Risk of sick leave associated with outdoor air supply rate, humidification, and occupant complaints. Indoor air, 10(4):212–221.
    OpenUrlCrossRefPubMedWeb of Science
  55. ↵
    Moghadas, S. M., Fitzpatrick, M. C., Sah, P., Pandey, A., Shoukat, A., Singer, B. H., and Galvani, A. P. (2020). The implications of silent transmission for the control of covid-19 outbreaks. Proceedings of the National Academy of Sciences, 117(30):17513–17515.
    OpenUrlAbstract/FREE Full Text
  56. ↵
    Morawska, L. and Cao, J. (2020). Airborne transmission of SARS-CoV-2: The world should face the reality. Environment International, 139:105730.
    OpenUrlCrossRefPubMed
  57. ↵
    Morawska, L., Johnson, G., Ristovski, Z., Hargreaves, M., Mengersen, K., Corbett, S., Chao, C., Li, Y., and Katoshevski, D. (2009). Size distribution and sites of origin of droplets expelled from the human respiratory tract during expiratory activities. Journal of Aerosol Science, 40:256–269.
    OpenUrlCrossRefWeb of Science
  58. ↵
    Morawska, L. and Milton, D. K. (2020). It is time to address airborne transmission of COVID-19. Clinical Infectious Diseases, 71:2311–2313.
    OpenUrlCrossRefPubMed
  59. ↵
    Moriarty, L. F. (2020). Public health responses to COVID-19 outbreaks on cruise ships worldwide, February–March 2020. MMWR. Morbidity and Mortality Weekly Report, 69.
  60. ↵
    Nicas, M., Nazaroff, W. W., and Hubbard, A. (2005). Toward understanding the risk of secondary airborne infection: emission of respirable pathogens. J. Occ. Env. Hygiene, 2(3):143–154.
    OpenUrl
  61. ↵
    Nishiura, H., Oshitani, H., Kobayashi, T., Saito, T., Sunagawa, T., Matsui, T., Wakita, T., Covid, M., and Suzuki, M. (2020). Closed environments facilitate secondary transmission of coronavirus disease 2019 (COVID-19). medRxiv.
  62. ↵
    Noakes, C., Beggs, C., Sleigh, P., and Kerr, K. (2006). Modelling the transmission of airborne infections in enclosed spaces. Epidemiology & Infection, 134(5):1082–1091.
    OpenUrlCrossRefPubMedWeb of Science
  63. ↵
    Noakes, C. J. and Sleigh, P. A. (2009). Mathematical models for assessing the role of airflow on the risk of airborne infection in hospital wards. Journal of the Royal Society Interface, 6(Suppl_6):S791–S800.
    OpenUrl
  64. ↵
    Oberg, T. and Brosseau, L. M. (2008). Surgical mask filter and fit performance. American journal of infection control, 36(4):276–282.
    OpenUrlCrossRefPubMedWeb of Science
  65. ↵
    Pan, J., Harb, C., Leng, W., and Marr, L. C. (2020). Inward and outward effectiveness of cloth masks, a surgical mask, and a face shield. medRxiv.
  66. ↵
    Peng, Z. and Jimenez, J. L. (2021). Exhaled co2 as covid-19 infection risk proxy for different indoor environments and activities. Environmental Science and Technology Letters.
  67. ↵
    Persily, A. and de Jonge, L. (2017). Carbon dioxide generation rates for building occupants. Indoor air, 27(5):868–879.
    OpenUrlCrossRef
  68. ↵
    Prather, K. A., Marr, L. C., Schooley, R. T., McDiarmid, M. A., Wilson, M. E., and Milton, D. K. (2020). Airborne transmission of sars-cov-2. Science, 370(6514):303–304.
    OpenUrlFREE Full Text
  69. ↵
    Prentiss, M. G., Chu, A., and Berggren, K. K. (2020). Superspreading events without superspreaders: Using high attack rate events to estimate n<sub>o</sub>for airborne transmission of COVID-19. medRxiv.
  70. Prill, R. et al. (2000). Why measure carbon dioxide inside buildings. Published by Washington State University Extension Energy Program WSUEEP07, 3.
    OpenUrl
  71. ↵
    Riediker, M. and Morawska, L. (2020). Low exhaled breath droplet formation may explain why children are poor SARS-CoV-2 transmitters. Aerosol and Air Quality Research, 20(7):1513–1515.
    OpenUrl
  72. ↵
    Riley, E. C., Murphy, G., and Riley, R. L. (1978). Airborne spread of measles in a suburban elementary school. American Journal of Epidemiology, 107(5):421–432.
    OpenUrlCrossRefPubMedWeb of Science
  73. ↵
    Rudnick, S. and Milton, D. (2003). Risk of indoor airborne infection transmission estimated from carbon dioxide concentration. Indoor air, 13(3):237–245.
    OpenUrlCrossRefPubMedWeb of Science
  74. ↵
    Salisbury, S. (1986). Measuring carbon dioxide levels as an indicator of poor building ventilation: A case study. In Proc IAQ’86, pages 78–82.
    OpenUrl
  75. ↵
    Santarpia, J. L., Herrera, V. L., Rivera, D. N., Ratnesar-Shumate, S., Reid, S., Denton, P. W., Martens, J., Fang, Y., Conoan, N., Callahan, M. V., Lawler, J. V., Brett-Major, D. M., and Lowe, J. J. (2020). The infectious nature of patient-generated SARS-CoV-2 aerosol. medRxiv.
  76. ↵
    Schwartz, A., Stiegel, M., Greeson, N., Vogel, A., Thomann, W., Brown, M., Sempowski, G. D., Alderman, T. S., Condreay, J. P., Burch, J., et al. (2020). Decontamination and reuse of N95 respirators with hydrogen peroxide vapor to address worldwide personal protective equipment shortages during the SARS-CoV-2 (COVID-19) pandemic. Applied Biosafety, 25(2):67–70.
    OpenUrl
  77. ↵
    Seppänen, O., Fisk, W., and Mendell, M. (1999). Association of ventilation rates and co2 concentrations with health andother responses in commercial and institutional buildings. Indoor air, 9(4):226–252.
    OpenUrlCrossRefPubMedWeb of Science
  78. ↵
    Shair, F. H. and Heitner, K. L. (1974). Theoretical model for relating indoor pollutant concentrations to those outside. Environmental Science & Technology, 8(5):444–451.
    OpenUrl
  79. ↵
    Shen, Y., Li, C., Dong, H., Wang, Z., Martinez, L., Sun, Z., Handel, A., Chen, Z., Chen, E., Ebell, M. H., et al. (2020). Community outbreak investigation of SARS-CoV-2 transmission among bus riders in eastern China. JAMA Internal Medicine.
  80. ↵
    Shendell, D. G., Prill, R., Fisk, W. J., Apte, M. G., Blake, D., and Faulkner, D. (2004). Associations between classroom co2 concentrations and student attendance in washington and idaho. Indoor Air, 14(5):333–341.
    OpenUrlCrossRefPubMedWeb of Science
  81. ↵
    Stilianakis, N. I. and Drossinos, Y. (2010). Dynamics of infectious disease transmission by inhalable respiratory droplets. Journal of the Royal Society Interface, 7(50):1355–1366.
    OpenUrl
  82. ↵
    Stutt, R. O. J. H., Retkute, R., Bradley, M., Gilligan, C. A., and Colvin, J. (2020). A modelling framework to assess the likely effectiveness of facemasks in combination with lock-down in managing the COVID-19 pandemic. Proc. R. Soc. A., 476:20200376.
    OpenUrlCrossRef
  83. ↵
    Sun, Y., Wang, Z., Zhang, Y., and Sundell, J. (2011). In china, students in crowded dormitories with a low ventilation rate have more common colds: evidence for airborne transmission. PloS One, 6(11):e27140.
    OpenUrlPubMed
  84. ↵
    van den Berg, P., Schechter-Perkins, E. M., Jack, R. S., Epshtein, I., Nelson, R., Oster, E., and Branch-Elliman, W. (2021). Effectiveness of three versus six feet of physical distancing for controlling spread of covid-19 among primary and secondary students and staff: A retrospective, state-wide cohort study. Clinical Infectious Diseases.
  85. ↵
    Van Doremalen, N., Bushmaker, T., Morris, D. H., Holbrook, M. G., Gamble, A., Williamson, B. N., Tamin, A., Harcourt, J. L., Thornburg, N. J., Gerber, S. I., et al. (2020). Aerosol and surface stability of SARS-CoV-2 as compared with SARS-CoV-1. New England Journal of Medicine, 382(16):1564–1567.
    OpenUrlCrossRefPubMed
  86. ↵
    Volz, E., Mishra, S., Chand, M., Barrett, J. C., Johnson, R., Geidelberg, L., Hinsley, W. R., Laydon, D. J., Dabrera, G., O’Toole, Á., Amato, R., Ragonnet-Cronin, M., Harrison, I., Jackson, B., Ariani, C. V., Boyd, O., Loman, N., McCrone, J. T., Gonçalves, S., Jorgensen, D., Myers, R., Hill, V., Jackson, D. K., Gaythorpe, K., Groves, N., Sillitoe, J., Kwiatkowski, D. P., Flaxman, S., Ratman, O., Bhatt, S., Hopkins, S., Gandy, A., Rambaut, A., and Ferguson, N. M. (2021). Transmission of SARS-CoV-2 lineage B.1.1.7 in England: Insights from linking epidemiological and genetic data. medRxiv.
  87. ↵
    Watanabe, T., Bartrand, T. A., Weir, M. H., Omura, T., and Haas, C. N. (2010). Development of a dose-response model for SARS coronavirus. Risk Analysis, 30(7):1129–1138.
    OpenUrlCrossRefPubMed
  88. ↵
    Wells, W. F. (1955). Airborne Contagion and Air Hygiene: An Ecological Study of Droplet Infections. Harvard University Press.
  89. ↵
    Yang, W. and Marr, L. C. (2011). Dynamics of airborne influenza a viruses indoors and dependence on humidity. PLOS ONE, 6(6):1–10.
    OpenUrlCrossRefPubMed
  90. ↵
    Zhang, J., Litvinova, M., Liang, Y., Wang, Y., Wang, W., Zhao, S., Wu, Q., Merler, S., Viboud, C., Vespignani, A., et al. (2020a). Changes in contact patterns shape the dynamics of the COVID-19 outbreak in china. Science, 368:1481–1486.
    OpenUrlAbstract/FREE Full Text
  91. ↵
    Zhang, R., Li, Y., Zhang, A. L., Wang, Y., and Molina, M. J. (2020b). Identifying airborne transmission as the dominant route for the spread of COVID-19. Proc. Nat. Acad. Sci., 117(26):14857–14863.
    OpenUrlAbstract/FREE Full Text
  92. ↵
    Zhu, Y., Bloxham, C. J., Hulme, K. D., Sinclair, J. E., Tong, Z. W. M., Steele, L. E., Noye, E. C., Lu, J., Xia, Y., Chew, K. Y., Pickering, J., Gilks, C., Bowen, A. C., and Short, K. R. (2020). A meta-analysis on the role of children in SARS-CoV-2 in household transmission clusters. Clinical Infectious Diseases.
Back to top
Posted April 16, 2021.
Download PDF
Data/Code
Email
Monitoring carbon dioxide to quantify the risk of indoor airborne transmission of COVID-19
Martin Z. Bazant, Ousmane Kodio, Alexander E. Cohen, Kasim Khan, Zongyu Gu, John W. M. Bush
medRxiv 2021.04.04.21254903; doi: https://doi.org/10.1101/2021.04.04.21254903
Twitter logo Facebook logo LinkedIn logo Mendeley logo
Citation Tools

Subject Area