No magic bullet: limiting in-school transmission in the face of variable SARS-CoV-2 viral loads =============================================================================================== * Debra Van Egeren * Madison Stoddard * Abir Malakar * Debayan Ghosh * Antu Acharya * Sk Mainuddin * Biswajit Majumdar * Deborah Luo * Ryan Nolan * Diane Joseph-McCarthy * Laura F. White * Natasha S. Hochberg * Saikat Basu * Arijit Chakravarty ## Abstract In the face of a long-running pandemic, understanding the drivers of ongoing SARS-CoV-2 transmission is crucial for the rational management of COVID-19 disease burden. Keeping schools open has emerged as a vital societal imperative during the pandemic, but in-school transmission of SARS-CoV-2 can contribute to further prolonging the pandemic. In this context, the role of schools in driving SARS-CoV-2 transmission acquires critical importance. Here we model in-school transmission from first principles to investigate the effectiveness of layered mitigation strategies on limiting in-school spread. We examine the effect of masks and air quality (ventilation, filtration and ionizers) on steady-state viral load in classrooms, as well as on the number of particles inhaled by an uninfected person. The effectiveness of these measures in limiting viral transmission is assessed for variants with different levels of mean viral load (Wuhan, Delta, Omicron). Our results suggest that a layered mitigation strategy can be used effectively to limit in-school transmission, with certain limitations. First, poorly designed strategies (insufficient ventilation, no masks, staying open under high levels of community transmission) will permit in-school spread even if some level of mitigation is ostensibly present. Second, for viral variants that are sufficiently contagious, it may be difficult to construct any set of interventions capable of blocking transmission once an infected individual is present, underscoring the importance of other measures. Our findings provide several practical recommendations: the use of a layered mitigation strategy that is designed to limit transmission, with other measures such as frequent surveillance testing and smaller class sizes (such as by offering remote schooling options to those who prefer it) as needed. ## Introduction Continued high levels of SARS-CoV-2 transmission strain healthcare systems and accelerate viral evolution, which undermines vaccinal efficacy1,2, and generates new variants with unpredictable epidemiological characteristics. For example, the transmissibility of SARS-CoV-2 has increased over time, with the Delta variant being around twice as transmissible3 as the original Wuhan strain, which had an *R* (reproduction number) of between 3.34 and 5.75. The incubation time of the disease has also demonstrated evolutionary change, going from 5-21 days for the original Wuhan strain6,7 to 2-4 days for the Omicron variant8. In order to slow viral evolution by limiting transmission, it is thus important to understand the role of schools in facilitating SARS-CoV-2 spread. In most countries, a significant fraction of the population consists of K-12 students, staff and first-degree household contacts of students and staff. In the US, 40% of all households have a child at home under 18 years of age9, and 23% of the US population is enrolled in school10. Thus, in-school transmission of SARS-CoV-2 will substantially impact transmission dynamics in the whole population. A number of studies during the early part of the pandemic led to the perception that SARS-CoV-2 did not spread in schools, based on the similarity in case counts between schools and their surrounding communities and a lack of observed transmission chains among children in schools. However, the methodological validity of these conclusions is debatable, as the metrics being used to infer a lack of spread are themselves vulnerable to a “false negative” problem (absence of evidence is not evidence of absence) (see Supplementary Text S2, and 11 for a detailed critique on the limitations of current research arguing that SARS-CoV-2 does not spread in schools). In fact, there is now a robust body of evidence supporting the contention that SARS-CoV-2 spreads efficiently in schools that lack adequate infection control measures. Empirical analyses using county-level panel data in the United States have demonstrated that counties with fully open K– 12 schools with in-person learning had a 5% increase in the growth rate of COVID-19 cases during April-December 2020 (a period of time when US schools were largely closed for the first five months)12. Consistent with this finding, COVID-19 symptom reporting was more common in areas where schools were open compared to areas with remote learning, an effect that was attenuated in communities using multiple mitigation measures13. In-school transmission is apparent when systematic surveillance testing methods are used14–17, and dramatic increases in case detection rates have also been observed in studies that relied on surveillance as opposed to symptomatic testing18. A layered mitigation strategy is one way to limit transmission in the school setting under conditions of widespread community transmission. There are multiple potential interventions available at this point: vaccines, masks, air quality improvements, surveillance testing, contact tracing/ isolation and podding. For reasons of cost and practicality, it is rational to seek a minimal set of infection-control measures. The challenge in this regard is that ongoing viral evolution can yield further changes to the characteristics of the virus, and a set of infection control measures that works well for one viral variant may readily be defeated by the next. An important open question then is: “What is the design of a minimal set of infection-control measures in schools that is robust to variant-to-variant differences in viral load?” In this work, we use mathematical modeling of the steps involved in viral transmission to understand the impact of infection-control measures in a range of different scenarios corresponding to variants with differing viral loads. Our intent was to address both the feasibility and robustness of strategies for limiting SARS-CoV-2 transmission in schools. ## Results To study SARS-CoV-2 transmission and the impact of control measures in schools, we created a multistep mathematical model of viral transmission in indoor settings (Methods) occurring under the assumption of indoor aerosol spread of infectious virus (see Supplementary Text S3 for justification of this assumption). First, we estimated the concentration of virions in the air over time in a room with an infected individual present using a differential equations model. This model assumes that infected individuals emit virions into the air at a constant rate into a room that is modeled as a well-mixed container (see Supplementary Text S4 for details about this assumption). Emitted virions can be inactivated over time or filtered out via air exchanges. These viral concentration estimates were then used as an input to calculate the probability that uninfected individuals in the room will become infected with SARS-CoV-2. ### Viral concentrations in a room with an infected individual present reach steady state quickly Under this model of SARS-CoV-2 emission, the viral concentration in a room reaches a steady-state concentration after the infected individual has been present for a certain time. For a closed room without air exchange or filtration, this concentration *C*SS depends on the emission rate *α*, viral inactivation rate *δ*, and the volume of the room *V* and is estimated as *C**SS* = *α*/(*δV*) (Methods). To understand the impact of evolutionary changes in viral load, we picked three variants of the virus whose viral loads have been reported in the literature to be widely different from one another (see Supplementary Text S5 for details about viral loads for the Wuhan strain and the Delta and Omicron variants). For a typical classroom size (20 feet x 20 feet x 10 feet), the steady-state concentration with an individual infected with the original Wuhan virus is approximately 4 virions per liter, which scales linearly for higher viral load SARS-CoV-2 strains (Fig. 1). We found that this steady-state concentration is reached after an infected individual is present in a room for 1-2 hours, regardless of the rate at which the individual emits viral particles. Masking of the emitter decreases the rate at which infected individuals release virions, lowering the final steady-state concentration, but does not affect the time at which the steady state is reached (see Supplementary Text S6 for details about parameter estimates for various interventions in this and subsequent sections). We further explore the effects of different types of masks on the steady-state viral concentration later in this section (Fig. 6). After the infected individual leaves the room, the viral concentration returns to approximately zero over a 5-hour period in a closed room (Fig. S1). ![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/31/2022.03.25.22272956/F1.medium.gif) [Figure 1:](http://medrxiv.org/content/early/2022/03/31/2022.03.25.22272956/F1) Figure 1: Viral concentrations in a room with an infected individual present reach steady state quickly. **A:** Infected individual with a low viral load strain of virus (similar to the initial Wuhan strain). **B:** Individual infected with an intermediate viral load variant (similar to Omicron). **C:** Individual infected with a high viral load variant (similar to Delta). Room dimensions are 10’ x 20’ x 20’. No air filtration or other mitigation methods were used. In all panels, the blue curve shows the concentration when the infected individual has no mask on, the orange curve shows the concentration when the infected individual has on a mask which filters out 50% of exhaled particles (typical cloth mask), and the green curve shows the concentration when the infected individual has on a mask which filters out 90% of exhaled particles (typical surgical mask). ### Steady-state concentrations of virions are reduced by high-volume air filtration Filtering the air in the room decreases the steady-state concentration of virions by increasing the rate at which infectious virions are eliminated from the room (Methods). For infected individuals with a low viral emission rate (e.g., infected with the original Wuhan SARS-CoV-2 virus), air filtration and masking of the infectious individual can reduce the steady-state concentration of virus in the room to less than one virion per liter of air (Fig. 2). However, for SARS-CoV-2 variants and individuals with higher viral loads and therefore higher viral emission rates (e.g., with Omicron-like viral loads, represented here as the ‘medium’ rate, or with Delta-like viral loads, represented here as the ‘high’ emission rate), the steady-state viral concentration in the room might remain high even with high filtration efficiency and filtration rate (Fig. 2). Notably, for schools using high efficiency filtration systems (e.g., HEPA filters), increasing the rate of air exchange across the filter is important for minimizing transmission. ![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/31/2022.03.25.22272956/F2.medium.gif) [Figure 2:](http://medrxiv.org/content/early/2022/03/31/2022.03.25.22272956/F2) Figure 2: Steady-state concentrations of virions are strongly impacted by air exchange rate across the filter and by the fraction of particles removed when air is passed over the filter. Virion concentrations at steady state, in a room 10’ x 20’ x 20’. **A-C:** Individual infected with a low viral load (A), intermediate viral load (B), or high viral load (C) strain, no mask. **D-F:** Individual infected with a low viral load (D), intermediate viral load (E), or high viral load (F) strain, wearing a 90% effective mask. Typical air purification systems can achieve a filtration rate of approximately 5-6 exchanges/hr (Supplementary Text S6). ### Steady-state concentrations of virions in the air are decreased by ionizers Another strategy to lower the viral concentration is to use ionizers, which inactivate viral particles and remove them from the air19. These devices produce small ions by the corona discharge principle, according to which negatively charged ions transfer their charge to suspended particles upon collision. These charged particles then agglutinate, becoming larger until they fall out of the air under the effect of gravity20. Ionizers generating negatively charged ions have been shown to be efficient at removing bacteria, molds, and viruses from indoor air21–23. The efficiency of particle removal is dependent on the emission rate of ions within an enclosed space, as well as room volume24. Studies conducted with smoke particles in an enclosed room suggest a high efficiency of ionizers in removing particles from the air, that varies between 80 and 100%20,24,25. Although older ionizer technologies generate ozone, which is an undesirable byproduct, newer ionizers do not have this potential liability associated with them26. As with the other control measures we simulated, we found that ionizers can lower viral concentrations in a typical classroom to below one virion per liter in situations where the infected individual is emitting viral particles at a relatively low rate (masked, infected with a low viral load strain) (Fig. 3). However, if the individual is infected with a high viral load SARS-CoV-2 variant (e.g., Delta), the viral concentration in the room will be very high even when ionizers are being used (Fig. 3). ![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/31/2022.03.25.22272956/F3.medium.gif) [Figure 3:](http://medrxiv.org/content/early/2022/03/31/2022.03.25.22272956/F3) Figure 3: Steady-state concentrations are decreased by ionizers (which increase the particle degradation rate). Virion concentrations at steady state, in a room with the same dimensions as Figure 1. **A-C:** Individual infected with a low viral load (A), intermediate viral load (B), or high viral load (C) strain, no mask. **D-F:** Individual infected with a low viral load (D), intermediate viral load (E), or high viral load (F) strain, wearing a 90% effective mask. The degradation rate of airborne SARS-CoV-2 without the use of an ionizer is approximately 0.01/min. ### Only a fraction of inhaled viruses is deposited in the nasopharynx We used the viral concentration estimates calculated in the previous section, estimates of aerosolized particle deposition derived from computational fluid dynamics (CFD) modeling, and the minimum infectious dose of SARS-CoV-2 to estimate the probability that uninfected individuals in the room will become infected with SARS-CoV-2. Infection of a new host requires the virus to be inhaled and deposited on the airway mucosa, where it can replicate. To study transmission dynamics in indoor settings with airborne SARS-CoV-2, we used a published computational fluid dynamics model of airflow in the nasopharynx to estimate the number of inhaled viral particles that are deposited in the airway mucosa27. This computational model of airflow in the human nasopharynx estimates the probability that an inhaled virion will reach the airway mucosa, given the size of the liquid droplet in which it is suspended. We assumed individuals breathe in virions suspended in liquid droplets with the measured steady-state size distribution of expelled respiratory droplets28. We used the computational fluid dynamics modeling results to compute the overall probability that an inhaled virion will hit the nasopharynx, marginalized over the empirical droplet size distribution (Methods, Table S1), and found that approximately 0.6% of virions that are inhaled during each breath are deposited in the mucosa. ### Masks and air filtration reduce the expected number of virions transmitted to uninfected individuals We used this nasopharyngeal deposition probability to estimate the number of virions that reach the airway mucosa per hour in an uninfected student in a classroom with a given concentration of SARS-CoV-2 in the air (Methods). For a classroom with an individual emitting virions at a low rate, interventions such as masking and air filtration can lower the transmission rate to an uninfected individual to less than 1 virion per hour, preventing transmission for short periods in the classroom (Fig. 4). However, exposure over several hours or an entire school day may transmit enough virions to the nasopharynx to cause infection, so limiting contact and masking of uninfected students is still important for controlling viral spread in this situation. If the virus-emitting individual is infected with a high viral load variant (similar to Delta), the number of virions inhaled can increase to hundreds or thousands per hour, possibly causing infection even after only a short exposure period (Fig. 4). In fact, while several hours are required for transmission of a low viral load (Wuhan-like) variant in a classroom, an individual infected with a high viral load (Delta-like) variant can spread the virus within a matter of minutes (Fig. 5). Thus, in this scenario, with multiple environmental control measures (masking and air filtration) in place, transmission can still occur over an 8-hour school day. Thus, in a future scenario where we are faced with a Delta-like variant with high viral loads, schools may have to rely on additional measures (e.g., rapid and widespread surveillance testing, targeted school closures) to limit in-school spread. These measures would need to supplement and not replace in-school mitigation measures in order to be effective. ![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/31/2022.03.25.22272956/F4.medium.gif) [Figure 4:](http://medrxiv.org/content/early/2022/03/31/2022.03.25.22272956/F4) Figure 4: Even with masks and air filtration, a person exposed to a high viral load emitter will be exposed to thousands of virions per day. Virions inhaled by an uninfected individual per hour in the presence of an infected individual, in a room with the same dimensions as Figure 1. The probability that an inhaled virion hits nasopharynx was estimated at 0.6% (Methods). **A-C:** Individual infected with a low viral load (A), intermediate viral load (B), or high viral load (C) strain, no mask. **D-F:** Individual infected with a low viral load (D), intermediate viral load (E), or high viral load (F) strain, with the emitter wearing a 90% effective mask (one-way masking). ![Figure 5:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/31/2022.03.25.22272956/F5.medium.gif) [Figure 5:](http://medrxiv.org/content/early/2022/03/31/2022.03.25.22272956/F5) Figure 5: The presence of robust mitigation measures substantially increases the time to infection for low viral load strains. Plots show the time until an uninfected individual receives 500 virions in a room occupied by an infected individual. **A-C:** Individual infected with a low viral load (A), intermediate viral load (B), or high viral load (C) strain, no mask. **D-F:** Individual infected with a low viral load (D), intermediate viral load (E), or high viral load (F) strain, with the emitter wearing a 90% effective mask (one-way masking). ### Two-way masking can substantially reduce the number of virions transmitted High-quality masks, used correctly, have been demonstrated to reduce infection risk even in high-risk settings29–31. Therefore, such masks can provide an additional level of mitigation in the event that schools are faced with a high viral load (Delta-like) variant. Masks provide a double benefit, as they reduce transmission by filtering out virions emitted by infected individuals and by reducing the number of ambient virions inhaled by uninfected individuals. Masks that filter out >95% of virions increase the time to transmission by approximately 10-fold when worn on either the infected or uninfected individual (one-way masking; Fig. 6). If both the infected and uninfected individuals are masked (two-way masking), less-effective masks (e.g., well-fitted surgical masks) can achieve the same level of protection against transmission. Combining universal N95 masking with excellent ventilation can increase the time to transmission of even high viral load strains to longer than a typical school day (Fig. 6D), suggesting that layered mitigation strategies featuring well-fitted and high-quality masks are critical for the control of high viral load strains in the classroom. ![Figure 6:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/31/2022.03.25.22272956/F6.medium.gif) [Figure 6:](http://medrxiv.org/content/early/2022/03/31/2022.03.25.22272956/F6) Figure 6: Masking both infected and uninfected individuals reduces the rate of viral transmission. **A:** Steady-state airborne viral concentration in an unventilated classroom with a single infected individual wearing a mask with the filtration efficiency given on the x-axis. **B:** Time until infection of uninfected individuals present in an unventilated classroom with a single infected individual wearing a mask with the filtration efficiency given on the x-axis. Uninfected individuals either are not wearing masks (one-way masking; dashed lines) or are wearing a mask with the same filtration efficiency as the infected individual (two-way masking; solid lines). **C:** Steady-state airborne viral concentration in a well-ventilated classroom (6 complete air exchanges per hour) with a single infected individual wearing a mask with the filtration efficiency given on the x-axis. **D:** Time until infection of uninfected individuals present in a well-ventilated classroom with a single infected individual wearing a mask with the filtration efficiency given on the x-axis. Uninfected individuals either are not wearing masks (one-way masking; dashed lines) or are wearing a mask with the same filtration efficiency as the infected individual (two-way masking; solid lines). In all panels, shaded regions denote typical filtration efficiencies for cloth, surgical, and N95 masks (Supplementary Text S6) and curve colors denote viral load of the emitter. ### In a well-ventilated room, risk is strongly dependent on seating position relative to the infected individual Up to this point, we have considered the classroom to be a well-mixed container. This simple modeling approach allows us to identify settings where the risk of in-school transmission is high. The well-mixed container assumption is justified, particularly in settings with limited ventilation (see Supplementary Text S4 for details). However, in some settings, the assumption may not hold, particularly in well-ventilated rooms, which are thought to have a low risk of transmission overall. To better understand the risk of transmission in a well-ventilated setting, we used CFD simulations of air transport inside a classroom. To simulate a best-case scenario for ventilation, we chose to model airflow behavior with the dimensions of a large auditorium-style classroom in a tropical setting, with high ceilings and several windows that are all open (Fig. 7A-C). We placed a single infected individual (emitter) in the room and varied their position to understand the impact of the airflow in the room. ![Figure 7:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/31/2022.03.25.22272956/F7.medium.gif) [Figure 7:](http://medrxiv.org/content/early/2022/03/31/2022.03.25.22272956/F7) Figure 7: In a well-ventilated room, infection risk is strongly dependent on the location of the infected individual. **A:** Digitized geometry of a classroom (based on direct measurements of an auditorium-style classroom in a tropical setting). The position of the instructor’s head is marked by *p*. **B:** Meshed space interior to the bounds of the classroom. **C:** Typical head geometry for the emitter subtracted from the classroom mesh, with realistic nostril diameters. **D**: Mixing of velocity streamlines in the classroom with the windows as the inlet and the door as the outlet. The panel includes the layout and plan dimensions of the ventilated classroom. **E-H:** Flowlines of exhaled virions showing the impact of emitter location (blue dot: p represents the instructor’s head location, *s**1-3* represent the students’ head locations) on local concentrations of virions. Note that *s**1* is in the front seat, *s**2* is in the middle of the classroom, and *s**3* is in the rear seat. Local concentrations are dependent on fluid dynamics within the classroom, suggesting that even if the average concentration of virions in a room is below the infectious threshold, individuals may become infected over time based on their location downstream of the emitter. Figure 7D visually depicts the airflow mixing trends inside the room and the virion-bearing streamline patterns for emitters located at different parts of the auditorium-style classroom. As the room is well-ventilated, the aerosolized virions emitted by the teacher (Fig. 7E) standing on the podium at the front of the room (with a door located proximally on the side), as well as those from the student seated in the front row (Fig. 7F), would escape through the door quite readily. The situation is, however, different if the infected individual were seated in the middle of the room (Fig. 7G), and even worse if the infected individual were seated at the rear (Fig. 7H). In these two situations, the infected individual would be efficiently spreading aerosolized pathogens, via exhaled respiratory ejecta, through the entire room. Thus, for a well-ventilated room (where the well-mixed container assumption cannot be expected to apply) total viral load in the room depends strongly on the position of the infected individual in the room. Additionally, local virion concentrations may be sufficiently high to enable efficient viral spread in the absence of other countermeasures (such as masking). The chaotic airflow patterns (invisible to the naked eye) underscore the unpredictable downside of infection risk in a closed setting (Fig. 7H). Thus, while it is possible to use model-based approaches to identify settings with a high risk of transmission, model-based approaches that rely on the well-mixed container assumption cannot definitively identify indoor settings with a low risk of transmission. This finding further underscores the need for multiple layers of intervention, and a robust ability to detect outbreaks before they spread. ### Measures to reduce the likelihood of infected individuals being present in the classroom setting are crucial Environmental control measures (masking, air filtration, ionizers) can all have an impact on limiting transmission in the in-school setting. However, our work suggests that these measures-both individually and in concert-are all vulnerable to defeat by a sufficiently high burden of virion emission. Thus, reducing the probability of infected individuals being present in the classroom at all is crucial to limiting SARS-CoV-2 transmission in schools. To that end, three measures are critical for reducing the expected number of infected individuals in a classroom: testing, capacity limits, and targeted school closures. Regular screening tests, with infected individuals being identified and isolated before they enter the classroom, can reduce the number of expected infected individuals arriving at school each day. Reducing class sizes by running the school day in shifts or by offering a remote school option to students that prefer it can also reduce the expected number of infected individuals in schools (Fig. 8). Podding, which limits students’ exposure to others outside of their cohort, such as in the lunchroom, can interrupt transmission chains to prevent spread throughout the school. Similarly, keeping schools closed for periods of time when local transmission is high would also have a proportional impact on reducing transmission. Notably, widespread vaccine coverage -- which is essential to limit the mortality and morbidity burden of the ongoing COVID-19 pandemic--may not play a significant role in limiting transmission at present (see Supplementary Text S8 for a summary of vaccinal efficacy against transmission). ![Figure 8:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2022/03/31/2022.03.25.22272956/F8.medium.gif) [Figure 8:](http://medrxiv.org/content/early/2022/03/31/2022.03.25.22272956/F8) Figure 8: Reducing class sizes reduces the risk of infection in school settings. **A:** Probability that at least one student in a classroom is infected, given values for the population-level infection prevalence (y-axis) and the number of students in the classroom (x-axis). **B:** Expected number of students exposed to an infected individual in a 500-person school when the school is broken up into different size classes over the entire school day (see Methods for additional details), given different values for the initial overall SARS-CoV-2 infection prevalence in the student population (blue curve: 1% of students infected, orange curve: 5% infected, green curve: 10% infected). ## Discussion In this study, we have used mathematical modeling to demonstrate the strengths and limitations of a layered mitigation strategy in limiting in-school transmission while keeping schools open. We examined the impact of risk mitigation measures (masking, ionizers, ventilation, filtration) on limiting spread within a classroom when an infected person is present. Our findings underscore the critical importance of layered mitigation strategies in limiting in-school transmission. With that said, all of the examined measures can be readily defeated by sufficiently high viral loads, a biological change that has already been observed during the pandemic (for example between the Wuhan strain and the Delta variant). This is a crucial point: minimal effective measures for the disease as it is at present may have an increased risk of failure in the face of new variants of SARS-CoV-2. Our findings also indicate that the risk of transmission in schools may be hard to predict in certain settings (such as in the turbulent airflow patterns of a well-ventilated room). As a corollary, our work points to the central importance of relying on measures to limit the likelihood of having an infected person in the classroom: testing and isolation, limiting class sizes, and targeted closures when community transmission is high. The ongoing COVID-19 pandemic shows no signs of permitting a return to pre-pandemic life, with high rates of transmission leading to rapid viral evolution. The course of the pandemic under these conditions is likely to be unpredictable, with the potential for catastrophic levels of mortality and morbidity should new variants of the virus emerge that have greater immune evasion potential or virulence. Thus, slowing viral evolution by limiting transmission is now a vital societal imperative. There are many examples in human society where the routine operation of vital services provided by complex systems brings some measure of risk to those involved: road transportation, aviation, medical care, law enforcement, agriculture and power generation, to name a few. In each of these cases, risks are managed by using a systems approach32. The premise is that accidents in complex systems occur through the accumulation of failures. This “Swiss Cheese model”33 of risk mitigation is well known in the epidemiological community, and a number of epidemiologists have advocated for its use from the beginning of this pandemic34,35. The strength of the Swiss Cheese model in risk management lies in building a system that is robust to human error. In the specific case of an evolving virus, the Swiss Cheese model also provides multiple orthogonal selection pressures, making escape more difficult for the virus. To keep schools open, designing guidelines to limit transmission based on the Swiss Cheese model, and then validating those guidelines to ensure robustness to epidemiological changes driven by viral evolution, will be critical. In the school setting, the application of the Swiss Cheese model has been slow for a number of reasons: the risk of in-school transmission has been under-estimated or compared to the wrong outcome (health effects on children), and false dichotomies in strategic thinking have led to an inappropriate focus on some layers of the “Swiss Cheese” in preference to others (for example, arguing for unmasking because children are vaccinated, or arguing for eliminating testing and contact tracing because children are masked). At this point, there is clear evidence supporting the contention that SARS-CoV-2 transmission can occur in schools11,36. The public-health consequences of that transmission are not borne by children alone. In addition to the first-order effect of household transmission of school-acquired COVID-19, transmission in the school setting facilitates viral evolution. An evolving virus benefits from a narrow focus on individual measures-the more focused the selection pressure, the easier it is for the virus to escape it. The mistake of over-reliance on one layer of the Swiss Cheese Strategy (made with the vaccines) can easily be repeated with other measures (such as testing and contact tracing). This work has a number of assumptions and key limitations. We have assumed for most of the work that the classroom is a well-mixed container, and then demonstrated (Figure 7) that the failure of this assumption leads to higher risk than could be estimated from a well-mixed container assumption. We also assume that children are equally susceptible and infectious as adults (see Supplementary Text S1 for an in-depth discussion of this assumption). It assumes perfect compliance with mask-wearing, which is not likely to be true in practice37,38. We also have not considered the effect of vaccination on transmission or risk of infection-primarily because the role of vaccines in limiting SARS-CoV-2 infection (and transmission) has now been demonstrated to be highly time-sensitive and vulnerable to immune evasion (See Supplemental Materials S8 for an in-depth justification of this assumption). There are a number of very thoughtful modeling analyses on this topic that have been published throughout the course of this pandemic. Several other groups have used model-based analyses to demonstrate that in-school transmission of SARS-CoV-2 is likely to be significant12,13,39. The point that reopening schools without robust COVID-19 mitigation could lead to an acceleration of the pandemic has been made in several modeling studies40–42. At the same time, a number of groups have published models focused on a limited set of infection-control measures, to show how schools can be reopened without risking in-school transmission (for example by using portable air purifiers43, limiting class sizes44, or mandating vaccination45,46). Our work adds to the discussion by pointing out the need to take evolutionary-driven epidemiological changes (for example due to viral load from one variant to the next) into account. In his book *The Black Swan* (now a classic in the risk-management community), author Nasim Nicholas Taleb argues that the key to risk management lies not so much in predicting the worst thing that could happen, but in making plans that are robust to that outcome. Our modeling demonstrates that a layered mitigation strategy, implemented properly, can curtail viral transmission under many circumstances. With that said, there are ways to implement infection-control measures that are ineffective, and measures that are effective in the presence of one viral variant can be readily rendered ineffective in the presence of another. Because many of the interventions have a nonlinear effect on risk mitigation, cutting corners on risk mitigation steps can degrade their utility very quickly, turning them into “hygiene theater”. For example, we found that two-way masking with N95 masks increases the time until SARS-CoV-2 transmission by multiple orders of magnitude, as compared to one-way masking with cloth or surgical masks (Fig. 6). As a corollary, it is crucial for schools to have controls in place to ensure that measures taken for mitigation are working as intended. For example, air quality can be monitored using carbon dioxide monitors, and mandatory (as opposed to opt-in) testing can be used to monitor the functional outcomes of in-school mitigation. Mitigation strategies should be pressure-tested using simple mathematical modeling approaches such as the one described in this paper. Thresholds for the acceptable performance of mitigation measures (for example, air filters or ionizers) should be updated periodically to reflect changes in viral epidemiology. Thus, only by iterative optimization of control measures can we expect to have effective suppression of in-school transmission of SARS-CoV-2. Taking a risk-management mindset to the problem of developing layered measures for infection control in the school setting is crucial at this stage. Optimizing such approaches to ensure robustness in the face of viral evolution may allow us to escape the false dichotomy of keeping schools open versus bringing the pandemic to an end. ## Methods To estimate the risk of viral transmission we built a multistep model comprising emission of viral particles into the air and subsequent inhalation and deposition into the nasopharynx of uninfected individuals (Fig. S2). ### Estimating viral steady-state concentrations in a room with an infected individual We designed a mathematical model of viral concentrations in a room in which an infected individual emits airborne SARS-CoV-2 particles at a constant rate. We assumed that the airborne virions are distributed uniformly throughout the room (i.e., the air in the room is well-mixed). Virions are inactivated at a constant rate and have an average lifetime of approximately 1.6 hours indoors at 73°F, 55% humidity47. For a closed room with no air exchange or filtration, the change in virion concentration over time is given by the differential equation ![Formula][1] where *α* is the viral emission rate, *δ* is the viral inactivation rate, *V* is the volume of the room, and *t* is the amount of time elapsed since the infected individual entered the room. Therefore, the viral concentration in the air over time *C*(*t*) is ![Formula][2] Air filtration removes virions in the room by filtering out a certain fraction of viruses that pass through the filter. Assuming the room is well-mixed, this adds an additional first-order virus removal process with rate *εβ*/*V* where *ε* is the fraction of virions that are eliminated while passing through the filter (between 0 and 1), *β* is the rate at which room air is passed through the filter (in units of air exchanges per hour), and *V* is the room volume. With filtration, the concentration in the air is ![Formula][3] Ionizers simply increase the virus inactivation rate *δ*. For the original Wuhan virus strain, we estimated the rate at which infected individuals emit viral particles into the air using published data which show that rooms with infected individuals have approximately 2-6 viral copies per L of air48. Assuming these measurements were taken when the virus was at a steady-state concentration in the room and that these individuals were in a reasonably sized room (4000 cubic feet of air), we estimated that the viral emission rate would be approximately 2319-6937 virions/h using our equation for the viral steady-state contribution shown in the Results. Therefore, we used a value within this range (5000 virions/h) for the emission rate for the low viral load simulations in this study. See Supplemental Materials for details about parameter estimates for various interventions. ### Estimating the risk of infection based on the steady-state viral concentration in a room To cause a new infection, virions must be inhaled by an uninfected host and be deposited in the airway mucosa, where they can replicate and cause disease. The rate at which virions are inhaled is the product of the respiratory tidal volume (volume of air inhaled during each breath) and the respiratory rate. However, computational fluid dynamics suggests that a minority of inhaled virions hit the nasopharyngeal mucosa, and this probability depends on the size of the liquid droplets containing the virions27. We used published estimates of the respiratory droplet size distribution28 and the probability of hitting the nasopharynx for different inhaled droplet sizes27 to estimate the overall fraction of inhaled virions that are deposited in the airway mucosa, marginalized over the empirical respiratory droplet size distribution, as ![Formula][4] where *s**i* is the fraction of exhaled droplets of size *i, d**i* is the probability of mucosal deposition for droplets of size *i*, and *v**i* is the fraction of aerosolized virions that are contained in droplets of size *i*. This fraction *v**i* is calculated as ![Formula][5] where *r**i* is the radius of the droplets of size *i*. To calculate the time until infection of an uninfected host, we conservatively assumed that 500 virions were required for infection (see Supplementary Text S7 for details about the minimum infectious dose). ### Numerical simulations of air mixing and exhalation patterns inside a realistic ventilated classroom Flow physics play a vital role in the distribution of airborne pathogens inside a confined space, e.g., a classroom. To model that, we have implemented state-of-the-art Computational Fluid Dynamics (CFD) simulations of: (a) air transport inside the room owing to incoming flux through the open windows; and (b) exhaled streamline patterns from human emitters situated at different points of the room. The digital reconstruction of the room is based on the measurements from an auditorium-style classroom. The room is sized at 13 m × 8 m × 5 m, with three windows (each 1 m × 1.25 m) and a door (1 m × 2.25 m). A head-sized sphere is used to stand in for the emitter’s head in each simulated flow model and is subtracted from the corresponding interior mesh. The sphere is provided with nasal protuberances, bearing nostrils realistically sized at 107.65 mm2 and 125.33 mm2, based on computed tomography reconstructions of human subjects27. For the emitter’s locations, we tested four different scenarios: (i) streamlines emitted by the teacher standing on a podium at the front of the classroom, (ii) a student seated in the front row, (iii) a student seated in the middle of the auditorium, and (iv) a student situated at the rear row. The classroom interior is meshed on ICEM CFD 2019 R3 (ANSYS Inc., Canonsburg, Pennsylvania) with tetrahedral elements bearing a maximum size of 0.15 m. The resolution was based on an earlier study49 on meshing requirements (maximum size < 0.2 m) to numerically model natural ventilation. Additionally, five layers of prism elements49 with a height ratio of 1.2 are extruded on the boundary walls. The meshed geometries are then imported to ANSYS Fluent 2019 R3 to track the air transport using the SST k-*ω*model50 to account for the turbulent flow scales. It is noted that when wind encounters a blocking effect on its path owing to a building, the velocity pressure is converted into static pressure. Consequently, on the windward side, the pressure would increase, with consequent pressure reduction on the leeward side. The static pressure gradient generated by wind on the building surface can be estimated by: ![Formula][6] where Δ*P**w* is the static pressure difference, *C**p* is the wind pressure coefficient, *ρ**o* is the density of outside air, and *U* is the mean wind velocity. The *C**p* value, which is independent of wind speed but depends upon wind direction and incident angle, can be calculated for low storied buildings (for up to 3 stories)51. The other physical parameters, namely *ρ**o* and *U* are taken as 1.139 Kg/*m*2 and 2.7 m/s, from earlier studies52 and meteorological data53, respectively. Air dynamic viscosity was assumed to be 1.9065 × 10−5 kg/m. However, there will be another pressure difference – one between the classroom and the outside corridor (next to the door). From established building code standards54, we can conclude that pressure differential Δ*P**indoor* between the classroom and corridor will be 5 – 20 Pa, resulting in the net pressure gradient from window-to-door (i.e., main inlet to main outlet) to be ![Formula][7] From the above assessments, the individual pressure gradients approximate to Δ*P**w* = 3.5 Pa and Δ*P**indoor* = 5 Pa. In the simulations: windows were taken as pressure inlets with 0-gauge pressure, nostrils were taken as velocity inlets with volumetric flux at 15 L/min27,55–58 to replicate gentle steady breathing, and at the door, we imposed a pressure outlet condition with negative gauge pressure of Δ*P**total*, i.e., -8.5 Pa. Wall boundary condition was mimicked as a stationary wall with no slip condition. ### Modeling impact of class size on SARS-CoV-2 exposure The probability that at least one student in a classroom arrives infected with SARS-CoV-2 was estimated from the overall prevalence of the infection using Poisson statistics. This probability *P* was estimated as ![Formula][8] where *n* is the number of students in the classroom and *f* is the fraction of infected individuals in the overall population. To estimate the expected number of exposed individuals in a school with *N* total students, the entire student population was broken into classes with *i* students each. The probability that at least one student in each class *P**i* was estimated using the above equation, and the expected number of students exposed to an infected individual was then calculated as *NP**i*. ### Data and materials availability Scripts implementing the viral concentration and transmission models are available on GitHub at [https://github.com/dvanegeren/covid-indoor-transmission](https://github.com/dvanegeren/covid-indoor-transmission). Input data used to estimate the fraction of inhaled droplets that hit the nasopharynx are also available in the same GitHub repository. ## Supporting information Supplementary Text S1-S8, Figs. S1 and S2, Table S1 [[supplements/272956_file02.pdf]](pending:yes) ## Data Availability Airborne viral concentration simulation data are reproducible using the code available at [https://github.com/dvanegeren/covid-indoor-transmission](https://github.com/dvanegeren/covid-indoor-transmission). All additional data are available from the authors upon reasonable request. [https://github.com/dvanegeren/covid-indoor-transmission](https://github.com/dvanegeren/covid-indoor-transmission) * Received March 25, 2022. * Revision received March 25, 2022. * Accepted March 31, 2022. * © 2022, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NoDerivs 4.0 International), CC BY-ND 4.0, as described at [http://creativecommons.org/licenses/by-nd/4.0/](http://creativecommons.org/licenses/by-nd/4.0/) ## References 1. 1.Van Egeren, D. et al. Risk of rapid evolutionary escape from biomedical interventions targeting SARS-CoV-2 spike protein. PLOS ONE 16, e0250780 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0250780&link_type=DOI) 2. 2.Hirabara, S. M. et al. SARS-COV-2 Variants: Differences and Potential of Immune Evasion. Front Cell Infect Microbiol 11, 781429 (2021). 3. 3.Liu, Y. & Rocklöv, J. The reproductive number of the Delta variant of SARS-CoV-2 is far higher compared to the ancestral SARS-CoV-2 virus. J Travel Med 28, taab124 (2021). 4. 4.Alimohamadi, Y., Taghdir, M. & Sepandi, M. Estimate of the Basic Reproduction Number for COVID-19: A Systematic Review and Meta-analysis. J Prev Med Public Health 53, 151– 157 (2020). 5. 5.Sanche, S. et al. High Contagiousness and Rapid Spread of Severe Acute Respiratory Syndrome Coronavirus 2. Emerg Infect Dis 26, 1470–1477 (2020). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F31%2F2022.03.25.22272956.atom) 6. 6.Grant, R. et al. Impact of SARS-CoV-2 Delta variant on incubation, transmission settings and vaccine effectiveness: Results from a nationwide case-control study in France. The Lancet Regional Health – Europe 13, (2022). 7. 7.Lauer, S. A. et al. The Incubation Period of Coronavirus Disease 2019 (COVID-19) From Publicly Reported Confirmed Cases: Estimation and Application. Ann Intern Med 172, 577– 582 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7326/M20-0504&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32150748&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F31%2F2022.03.25.22272956.atom) 8. 8.Jansen, L. et al. Investigation of a SARS-CoV-2 B.1.1.529 (Omicron) Variant Cluster — Nebraska, November–December 2021. MMWR Morb. Mortal. Wkly. Rep. 70, 1782–1784 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.15585/MMWR.MM705152E3&link_type=DOI) 9. 9.Statista Research Department. Share of U.S. family households with children, by type 1970-2020. Statista [https://www.statista.com/statistics/242074/percentages-of-us-family-households-with-children-by-type/](https://www.statista.com/statistics/242074/percentages-of-us-family-households-with-children-by-type/) (2021). 10. 10.US Census Bureau. School Enrollment in the United States: October 2020 - Detailed Tables. Census.gov [https://www.census.gov/data/tables/2020/demo/school-enrollment/2020-cps.html](https://www.census.gov/data/tables/2020/demo/school-enrollment/2020-cps.html) (2021). 11. 11.White, L. F., Murray, E. J. & Chakravarty, A. The role of schools in driving SARS-CoV-2 transmission: not just an open-and-shut case. (2022) doi:10.31219/osf.io/2eh5m. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.31219/osf.io/2eh5m&link_type=DOI) 12. 12.Chernozhukov, V., Kasahara, H. & Schrimpf, P. The association of opening K–12 schools with the spread of COVID-19 in the United States: County-level panel data analysis. Proceedings of the National Academy of Sciences 118, e2103420118 (2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxODoiMTE4LzQyL2UyMTAzNDIwMTE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDMvMzEvMjAyMi4wMy4yNS4yMjI3Mjk1Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 13. 13.Lessler, J. et al. Household COVID-19 risk and in-person schooling. Science 372, 1092– 1097 (2021). [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEzOiIzNzIvNjU0Ni8xMDkyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDMvMzEvMjAyMi4wMy4yNS4yMjI3Mjk1Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 14. 14.Stein-Zamir, C. et al. A large COVID-19 outbreak in a high school 10 days after schools’ reopening, Israel, May 2020. Euro Surveill 25, 2001352 (2020). 15. 15.Torres, J. P. et al. Severe Acute Respiratory Syndrome Coronavirus 2 Antibody Prevalence in Blood in a Large School Community Subject to a Coronavirus Disease 2019 Outbreak: A Cross-sectional Study. Clinical Infectious Diseases 73, e458–e465 (2021). 16. 16.Gold, J. A. W. et al. Clusters of SARS-CoV-2 Infection Among Elementary School Educators and Students in One School District — Georgia, December 2020–January 2021. MMWR Morb. Mortal. Wkly. Rep. 70, 289–292 (2021). 17. 17.Fontanet, A. et al. SARS-CoV-2 infection in schools in a northern French city: a retrospective serological cohort study in an area of high transmission, France, January to April 2020. Eurosurveillance 26, 2001695 (2021). 18. 18.Crowe, J. et al. Assessment of a Program for SARS-CoV-2 Screening and Environmental Monitoring in an Urban Public School District. JAMA Network Open 4, e2126447 (2021). 19. 19.Jiang, S.-Y., Ma, A. & Ramachandran, S. Negative Air Ions and Their Effects on Human Health and Air Quality Improvement. Int J Mol Sci 19, 2966 (2018). 20. 20.Uk Lee, B., Yermakov, M. & Grinshpun, S. A. Removal of fine and ultrafine particles from indoor air environments by the unipolar ion emission. Atmospheric Environment 38, 4815– 4823 (2004). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.atmosenv.2004.06.010&link_type=DOI) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000223646300002&link_type=ISI) 21. 21.Hyun, J., Lee, S.-G. & Hwang, J. Application of corona discharge-generated air ions for filtration of aerosolized virus and inactivation of filtered virus. Journal of Aerosol Science 107, 31–40 (2017). 22. 22.Wang, C., Lu, S. & Zhang, Z. Inactivation of airborne bacteria using different UV sources: Performance modeling, energy utilization, and endotoxin degradation. Science of The Total Environment 655, 787–795 (2019). 23. 23.Hagbom, M. et al. Ionizing air affects influenza virus infectivity and prevents airborne-transmission. Sci Rep 5, 11431 (2015). 24. 24.Sawant, V. S., Meena, G. S. & Jadhav, D. B. Effect of Negative Air Ions on Fog and Smoke. Aerosol Air Qual. Res. 12, 1007–1015 (2012). 25. 25.Grinshpun, S. A. et al. Evaluation of ionic air purifiers for reducing aerosol exposure in confined indoor spaces. Indoor Air 15, 235–245 (2005). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/j.1600-0668.2005.00364.x&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15982270&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F31%2F2022.03.25.22272956.atom) [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000230048500004&link_type=ISI) 26. 26.Sung, J.-H. et al. Ultrafine particle cleaning performance of an ion spray electrostatic air cleaner emitting zero ozone with diffusion charging by carbon fiber. Building and Environment 166, 106422 (2019). 27. 27.Basu, S. Computational characterization of inhaled droplet transport to the nasopharynx. Sci Rep 11, 6652 (2021). 28. 28.Xie, X., Li, Y., Sun, H. & Liu, L. Exhaled droplets due to talking and coughing. J R Soc Interface 6 Suppl 6, S703–714 (2009). 29. 29.Wang, X., Ferro, E. G., Zhou, G., Hashimoto, D. & Bhatt, D. L. Association Between Universal Masking in a Health Care System and SARS-CoV-2 Positivity Among Health Care Workers. JAMA 324, 703–704 (2020). [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F31%2F2022.03.25.22272956.atom) 30. 30.Jehn, M. et al. Association Between K–12 School Mask Policies and School-Associated COVID-19 Outbreaks — Maricopa and Pima Counties, Arizona, July–August 2021. MMWR Morb. Mortal. Wkly. Rep. 70, 1372–1373 (2021). 31. 31.CDC. Science Brief: Community Use of Masks to Control the Spread of SARS-CoV-2. Centers for Disease Control and Prevention [https://www.cdc.gov/coronavirus/2019-ncov/science/science-briefs/masking-science-sars-cov2.html](https://www.cdc.gov/coronavirus/2019-ncov/science/science-briefs/masking-science-sars-cov2.html) (2021). 32. 32.Reason, J. Human error: models and management. BMJ 320, 768–770 (2000). [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiYm1qIjtzOjU6InJlc2lkIjtzOjEyOiIzMjAvNzIzNy83NjgiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMi8wMy8zMS8yMDIyLjAzLjI1LjIyMjcyOTU2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 33. 33.Reason, J. The contribution of latent human failures to the breakdown of complex systems. Philosophical Transactions of the Royal Society of London. B, Biological Sciences 327, (1990). 34. 34.Noh, J. Y. et al. Safe hospital preparedness in the era of COVID-19: The Swiss cheese model. Int J Infect Dis 98, 294–296 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ijid.2020.06.094&link_type=DOI) 35. 35.Roberts, S. The Swiss Cheese Model of Pandemic Defense. The New York Times (2020). 36. 36.Hyde, Z. COVID-19, children and schools: overlooked and at risk. Medical Journal of Australia 213, 444-446.e1 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.5694/mja2.50823&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F31%2F2022.03.25.22272956.atom) 37. 37.Fischer, C. B. et al. Mask adherence and rate of COVID-19 across the United States. PLOS ONE 16, e0249891 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0249891&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33852626&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F31%2F2022.03.25.22272956.atom) 38. 38.Margraf, J., Brailovskaia, J. & Schneider, S. Adherence to behavioral Covid-19 mitigation measures strongly predicts mortality. PLOS ONE 16, e0249392 (2021). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0249392&link_type=DOI) 39. 39.Buja, A. et al. Opening Schools and Trends in SARS-CoV-2 Transmission in European Countries. Int J Public Health 66, 1604076 (2021). 40. 40.Gurdasani, D. et al. School reopening without robust COVID-19 mitigation risks accelerating the pandemic. The Lancet 397, 1177–1178 (2021). 41. 41.Zhang, Y. et al. COVID-19 Projections for K12 Schools in Fall 2021: Significant Transmission without Interventions. 2021.08.10.21261726 (2021) doi:10.1101/2021.08.10.21261726. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMS4wOC4xMC4yMTI2MTcyNnYyIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDMvMzEvMjAyMi4wMy4yNS4yMjI3Mjk1Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 42. 42.Panovska-Griffiths, J. et al. Determining the optimal strategy for reopening schools, the impact of test and trace interventions, and the risk of occurrence of a second COVID-19 epidemic wave in the UK: a modelling study. The Lancet Child & Adolescent Health 4, 817– 827 (2020). 43. 43.Narayanan, S. R. & Yang, S. Airborne transmission of virus-laden aerosols inside a music classroom: Effects of portable purifiers and aerosol injection rates. Physics of Fluids 33, 033307 (2021). 44. 44.Baxter, A., Oruc, B. E., Keskinocak, P., Asplund, J. & Serban, N. Evaluating Scenarios for School Reopening under COVID19. 2020.07.22.20160036 (2020) doi:10.1101/2020.07.22.20160036. [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4wNy4yMi4yMDE2MDAzNnYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjIvMDMvMzEvMjAyMi4wMy4yNS4yMjI3Mjk1Ni5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 45. 45.Panovska-Griffiths, J. et al. Modelling the impact of reopening schools in the UK in early 2021 in the presence of the alpha variant and with roll-out of vaccination against SARS-CoV-2. J Math Anal Appl 126050 (2022) doi:10.1016/j.jmaa.2022.126050. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jmaa.2022.126050&link_type=DOI) 46. 46.Head, J. R., Andrejko, K. L. & Remais, J. V. Model-based assessment of SARS-CoV-2 Delta variant transmission dynamics within partially vaccinated K-12 school populations. Lancet Reg Health Am 5, 100133 (2022). 47. 47.Science and Technology Directorate. Estimated Airborne Decay of SARS-CoV-2. US Department of Homeland Security [https://www.dhs.gov/science-and-technology/sars-airborne-calculator](https://www.dhs.gov/science-and-technology/sars-airborne-calculator) (2022). 48. 48.Santarpia, J. L. et al. Aerosol and surface contamination of SARS-CoV-2 observed in quarantine and isolation care. Sci Rep 10, 12732 (2020). [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41598-020-69286-3&link_type=DOI) [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32728118&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2022%2F03%2F31%2F2022.03.25.22272956.atom) 49. 49.Basu, S., Witten, N. & Kimbell, J. S. Influence of localized mesh refinement on numerical simulations of post-surgical sinonasal airflow. Journal of Aerosol Medicine and Pulmonary Drug Delivery 30, A–14 (2017). 50. 50.Kahsay, M. T., Bitsuamlak, G. T. & Tariku, F. CFD simulation of external CHTC on a high-rise building with and without façade appurtenances. Building and Environment 165, 106350 (2019). 51. 51.Li, Y. & Peterson, F. The modelling of multi-room air flow and its application to exhaust ventilation analysis. Climate and Buildings 2, 15–28 (1990). 52. 52.Rao, K. R. & Balkundi, H. V. Study of Variation in Upper Air Density Over India. Defence Science Journal 30, 31–44 (1980). 53. 53.India Meteorological Department. Climate of West Bengal. (2008). 54. 54.International Organization for Standardization. Cleanrooms and associated controlled environments — Part 4: Design, construction and start-up. (2001). 55. 55.Basu, S. et al. Numerical evaluation of spray position for improved nasal drug delivery. Sci Rep 10, 10568 (2020). 56. 56.Basu, S., Frank-Ito, D. O. & Kimbell, J. S. On computational fluid dynamics models for sinonasal drug transport: Relevance of nozzle subtraction and nasal vestibular dilation. Int J Numer Method Biomed Eng 34, e2946 (2018). 57. 57.Farzal, Z. et al. Comparative Study of Simulated Nebulized and Spray Particle Deposition in Chronic Rhinosinusitis Patients. Int Forum Allergy Rhinol 9, 746–758 (2019). 58. 58.Brandon, B. M. et al. Nasal Airflow Changes With Bioabsorbable Implant, Butterfly, and Spreader Grafts. The Laryngoscope 130, E817–E823 (2020). [1]: /embed/graphic-9.gif [2]: /embed/graphic-10.gif [3]: /embed/graphic-11.gif [4]: /embed/graphic-12.gif [5]: /embed/graphic-13.gif [6]: /embed/graphic-14.gif [7]: /embed/graphic-15.gif [8]: /embed/graphic-16.gif