Skip to main content

Estimating the nationwide transmission risk of measles in US schools and impacts of vaccination and supplemental infection control strategies

A Correction to this article was published on 24 July 2020

This article has been updated

Abstract

Background

The spread of airborne infectious diseases such as measles is a critical public health concern. The U.S. was certified measles-free in 2000, but the number of measles cases has increased in recent years breaking the record of the nationwide annual number of cases since 1992. Although the characteristics of schools have made them one of the most vulnerable environments during infection outbreaks, the transmission risk of measles among students is not completely understood. We aimed to evaluate how three factors influence measles transmission in schools: personal (vaccination), social (compartmentalizing), and building systems (ventilation, purification, and filtration).

Methods

We used a combination of a newly developed multi-zone transient Wells-Riley approach, a nationwide representative School Building Archetype (SBA) model, and a Monte-Carlo simulation to estimate measles risk among U.S. students. We compared our risk results with the range of reported transmission rates of measles in school outbreaks to validate the risk model. We also investigated the effectiveness of vaccination and ten supplemental infection control scenarios for reducing the risk of measles transmission among students.

Results

Our best nationwide estimate of measles transmission risk in U.S. schools were 3.5 and 32% among all (both unvaccinated and immunized) and unvaccinated students, respectively. The results showed the transmission risk of measles among unvaccinated students is > 70 times higher than properly immunized ones. We also demonstrated that the transmission risk of measles in primary schools (assuming teacher self-contained classrooms) is less than secondary schools (assuming departmentalized systems). For building-level interventions, schools with ductless-with-air-filter and ductless-without-air-filter systems have the lowest and highest transmission risks of measles, respectively. Finally, our simulation showed that infection control strategies could cut the average number of infected cases among all students in half when a combination of advanced air filtration, ventilation, and purification was adopted in the modeled schools.

Conclusions

Our results highlight the primary importance of vaccination for reducing the risk of measles transmission among students. Yet, additional and significant risk reduction can be achieved through compartmentalizing students and enhancing building ventilation and filtration systems.

Peer Review reports

Background

The spread of airborne infectious diseases is a global public health concern. Measles is an airborne viral respiratory illnesses that remains a significant cause of death worldwide and imposes an extreme economic burden on communities and families despite the availability of a safe and cost-effective vaccine [1]. It is estimated that approximately 110,000 people, mostly children under the age of 5 years, died from measles globally in 2017 [2]. High incidence of measles among healthy children gives rise to considerable morbidity [3], which, in turn, makes the role of school environments crucial in the spread of measles in the community. It is shown that the transmission of respiratory illnesses such as measles in schools leads to large excesses in expenses associated with healthcare, absence from school, and reduction in students’ productivity [4,5,6,7].

Concerns about measles is a public health issue which dates back many centuries. A brief history of measles tells us a Persian physician published one of the first written accounts of measles disease in the ninth century; in 1757 a Scottish doctor, Francis Home, demonstrated that measles is caused by an infectious agent; in 1912 measles became a nationally notifiable disease in the U.S.; and finally in 1963, the measles vaccine became available in the United States [8]. Since then, the annual number of measles cases in the U.S. decreased drastically from ~ 450,000 cases just before the introduction of the measles vaccine to less than 100 cases in 2000, when the U.S. was certified measles-free [9, 10]. The measles vaccine coverage in the U.S. has remained relatively constant since the Vaccines for Children program began in 1994 [11]. However, the number of annual measles cases has increased in recent years. In 2019, 1282 individual cases of measles have been confirmed in 31 states, which was the greatest number of cases reported in the U.S. since 1992 [12]. The majority of the cases have been among unvaccinated individuals. This demonstrates, primarily, the danger that unvaccinated cohorts deliver to the community including to infants younger than 12 months old, pregnant women, and individuals allergic to the vaccine or who haveweakened immune systems and thereforewho should not get vaccinated. Furthermore, this demonstrates the importance of understanding the transmission mechanisms of the measles virus in the built environment and adopting efficient control strategies to reduce the infection risk of measles among susceptible individuals. It is particularly important to investigate measles transmission among students at schools given that educational institutes are considered one of the most vulnerable environments in transmission of airborne infectious disease due to the extensive amount of time that students regularly spend in schools and the high levels of interactions among schoolchildren that occur in these environments [13,14,15].

The airborne transmission of measles in indoor environments such as schools and associated risk of infection presented to susceptible occupants are governed by several complex physical, biological, and epidemiological processes. Mathematical models have long been used to predict the transmission risk of airborne infectious diseases in the built environment. Epidemic modeling approaches such as susceptible-infector-susceptible (SIS) [16], susceptible-infectious-recovered (SIR) [17] competing-risks [18, 19], neural network [20], and Reed-Frost [21, 22] models are used to describe the progression of a disease in a population, although it is shown that these models alone cannot explain the spread of airborne infectious diseases such as measles in indoors environments [22]. Therefore, epidemic models are usually combined with other mathematical approaches to predict the risks associated with indoor spaces such as airplanes, hospitals [3], schools [17], residences [23], and healthcare facilities [24].

Another mathematical approach is the dose-response model, which has been adopted to estimate the airborne infection risks associated with the dose of infectious agents delivered to upper and lower respiratory tracts of a susceptible individual. This model requires the use of an underlying fate and transport model such as Markov chains, computational fluid dynamics (CFD), multi-zone mass balance models, or combinations of these methods to estimate the number of delivered viable pathogens to the infection sites in the respiratory tracts [25,26,27,28,29,30]. The complexity of combination of the dose-response approach and mathematical transport models limits the application of this simulation approach to only well-described indoor environments and diseases such as influenza.

Among all mathematical risk models, the Wells-Riley model is the most common approach [31], which was introduced originally by Wells [32] and Riley et al. [33]. The Wells–Riley model has been extensively used in analyzing ventilation strategy and its association to airborne infections in indoor environments and considered as a valid approach for estimating the transmission risk of airborne infectious diseases [34]. The Wells-Riley model is a relatively simple model that can be deployed with less required information regarding the characteristics of a desired space or disease in comparison to other mathematical risk approaches such as the dose-response model. Despite the validity and popularity of the Wells-Riley model, only a few studies used this approach to estimate the transmission risk of measles in an educational setting and all of them modeled the setting as one well-mixed indoor space [33, 35, 36]. This makes the results of the risk models less accurate, particularly for large and complex indoor environments. A limited number of studies used developed multi-zone versions of the Wells-Riley model for other indoor spaces such as hospitals and multifamily residential buildings [31, 37, 38], which demonstrates the potential capability of a multi-zone version of the model for predicting the infection risk of airborne infectious diseases such as measles in other complex indoor environments such as schools.

There is also a gap in understanding the efficiency of various infection control strategies beyond vaccination in reducing the transmission risk of airborne pathogens in indoor environments and particularly for lowering measles risk in schools. Many studies have demonstrated the most effective control strategy against measles is adequate vaccination [39,40,41] meaning children get two doses of vaccine, starting with the first dose at 12 through 15 months of age, and the second dose at 4 through 6 years of age and teens and adults should also be up to date on their vaccinations [42]. However, despite high vaccination records, explosive measles outbreaks may occur in school environments due to (i) the fact that a portion of adequately immunized individuals remains susceptible to measles viruses, (ii) high interaction and contact rates among students, or (iii) inadequate immunity from vaccinations at younger ages [22]. Heating, ventilation, and air-conditioning (HVAC) systems and air purifiers are shown to have a significant positive impact on reducing the transmission of measles in various indoor environments by enhancing filtration, ventilation and purification rates [43,44,45,46,47]; but, their bio-aerosol removal effectiveness for a typical school setup is understudied.

Moreover, in the best of our knowledge, all of the existing mathematical risk models have been applied to an indoor environment with specific building, HVAC, and occupational characteristics. These model outcomes normally could not be the representative of infection risk in a certain type of environment in a region or nationwide. Consequently, standard developers and policymakers are reluctant to use the existing model results for establishing new guidelines to control the transmission risk of infectious airborne diseases including measles in vulnerable indoor environments such as schools. Developing a nationwide representative School Building Archetype (SBA) model and combining it with a proper mathematical risk approach would form a powerful transmission risk tool that helps to fulfill this shortcoming in the knowledge of airborne infectious disease transmissions.

This research work primarily aims to investigate the transmission risk of measles in U.S. schools. Herein, we combined a newly developed multi-zone transient Wells-Riley model with a nationwide representative SBA model to estimate the transmission risk of measles in primary and secondary educational institutions in the U.S. and evaluate the performance of several control strategies for reducing measles infection risk.

Methods

Developing Wells-Riley model for several microenvironments

The Wells-Riley model was developed originally by Wells [32] and Riley et al. [33] to estimate the probability of airborne transmission of an infectious agent indoors. Rudnick and Milton [35] developed a new derivation of the Wells-Riley model in which the probability of infection transmission (Pinfection) in a well-mixed indoor space can be estimated using Eq. 1.

$$ {P}_{infection}=\frac{Number\ of\ Infected\ Cases}{Number\ of\ Susceptible\ Individuals}=1-{e}^{-\mu } $$
(1)

In this risk model, μ is the number of quantum of infection or ‘quanta’ breathed by a susceptible person. It is important to notice that quanta is not an actual physical unit; rather, it is a hypothetical infectious dose unit, which is typically back-calculated from observational epidemiological studies. Wells suggested the average risk of becoming infected after exposure to one quanta is 63% (i.e. 1 − e−1), which is essentially a 63% infectious dose, ID63 [32].

Herein, we considered three microenvironments or spaces within a typical U.S. school building to estimate the transmission risk of infectious aerosols between students (i.e. droplet nuclei containing measles viruses in this research work) after one index case (infector) enters the school, including:

  1. (i)

    Infector’s classroom: the school classroom, where the index case (infector) spends most of their time

  2. (ii)

    Recirculation spaces: the spaces within a typical school building (e.g. classrooms, labs, and hallways), where generated infectious bio-aerosols would reach there only via HVAC system air recirculation from the infector’s classroom

  3. (iii)

    Common spaces: spaces other than the infector’s classroom, where the index case physically presents for a considerable amount of time and interacts with other students

In this case, the average number of quanta breathed by susceptible students during a typical school day (\( \overline{\mu}\Big) \) can be estimated using Eq. 2.

$$ \overline{\mu}=\frac{1}{N_{total}}\times \overline{p}\times \sum \limits_i{\int}_0^{{\overline{t}}_i}{N}_i\left(\tau \right).{C}_{quanta,i}\left(\tau \right) d\tau $$
(2)

Ntotal: Total number of students in the schools during the infection period.

\( \overline{p} \): Average breathing rate of one student (m3 / hour).

\( {\overline{t}}_i \): Average time that students spend in space i (hour).

Ni(τ): Number of students in space i as a function of time.

Cquanta,i(τ): Concentration of quanta in space i, τ hours after the index case enters the space (quanta / m3).

We made several simplifications for estimating the average number of breathed quanta by susceptible students including:

  1. (i)

    Students stay continuously in each space during an exposure period

  2. (ii)

    The number of students in each space remains constant during an exposure period

  3. (iii)

    Other transmission pathways of measles viruses such as direct contact or fomite are ignored

It is noticeable that the potential impacts of these assumptions on the risk transmission results were taken into the account indirectly by back-calculating the quanta generation rate from actual measles outbreak cases in two U.S. schools (i.e. one elementary and one high school) with different interaction patterns among students. The back-calculation process is described in detail in the “back calculating quanta generation rate” section.

The concentration of quanta in the infector’s classroom and the common space τ hours after presence of the index case, Cqunata,j, can be estimated using Eq. 3, by solving a well-mixed mass balance equation for each of these two spaces as described in Appendix A.

$$ {C}_{quanta,j}\left(\tau \right)=\frac{Iq}{V_j{K}_{total,j}}\left(1-{e}^{-{K}_{total,j}\tau}\right) $$
(3)

I: Number of index cases.

q: Quanta generation rate (quanta / hour).

Vj: Volume of space j – either infector’s classroom or common space – (m3).

Ktotal,j: Total removal rate of measles viruses in space j – either infector’s classroom or common space – (per hour).

In this model, we assumed the recirculated air is the only pathway that the infectious particles can reach the recirculation space from the infector’s classroom. We adopted a discrete time-varying mass balance to estimate the concentration of quanta in the recirculation space. The concentration of quanta in the recirculation space at each time step [Cquanta,recir(τn)] was estimated using Eq. 4. The steps for developing Eq. 4 are shown in Appendix A.

$$ {C}_{quanta, recir}\left({\tau}_n\right)=\Delta \tau \left[-{K}_{total, recir}{C}_{quanta, recir}\left({\tau}_{n-1}\right)\left(\frac{Q_{return, class}{f}_{recir}{f}_{runtime}\left(1-{\eta}_{filter}\right)}{V_{recir}}\times \frac{Q_{supply, recir}}{Q_{supply, total}}\right){C}_{quanta, class}\left({\tau}_{n-1}\right)\right]+{C}_{quanta, recir}\left({\tau}_{n-1}\right) $$
(4)

∆τ: Time step interval, which is considered one minute in this model (hour).

Ktotal,recir: Total removal rate of measles viruses in recirculation space (per hour).

Qreturn,class: Return airflow rate of the infector’s classroom (m3/hour).

frecir: Fraction of recirculated air volume to total airflow capacity of HVAC system.

fruntime: Runtime fraction of HVAC system.

ηfilter: Removal efficiency of HVAC air filter.

Qsupply,recir: Supply airflow rate of the recirculation space (m3/hour).

Qsupply,total: Total supply capacity of HVAC system (m3/hour).

Cquanta,class(τn): Concentration of quanta in infector’s classroom (quanta / m3).

Vrecir: Volume of recirculation space (m3).

The total removal rate of measles viruses in space i (Ktotal,i) – infector’s classroom, recirculation space, and common space – was estimated by summing the rates of five infection removal mechanisms as shown in Eq. 5.

$$ {K}_{total,i}={\lambda}_{infiltration,i}+{K}_{deposition,i}+{K}_{ventilation,i}+{K}_{filtration,i}+{K}_{purification,i} $$
(5)

λinfiltration,i: Natural air ventilation rate or infiltration air exchange rate in space i (per hour).

Kdeposition,i: Deposition rate of measles particles in space i (per hour).

Kventilation,i: Mechanical ventilation rate of HVAC system in space i (per hour).

Kfiltration,i: Infectious particle removal rate due to filtration in space i (per hour).

Kpurification,i: Removal rate of infectious particles by standalone air handling units (AHU) or air purifiers in space i (per hour).

The infiltration air exchange rate (λinfilteration) of a typical school in the U.S. was estimated to be 0.31 per hour ranging between 0.12 and 0.49 per hour based on the U.S. Department of Energy (DOE) commercial reference-building models for educational buildings [48]. The DOE’s commercial reference-building model represents approximately two-thirds of the national building stock in the U.S.

Deposition (Kdeposition) and filtration (Kfiltration) rates of infectious particles in indoor environments depend on the distribution of measles viruses in different bio-aerosol size ranges. We are not aware of any study that reports the size distribution of measles viruses in indoor aerosols. Several studies have reported the size distribution of dry powder measles vaccine for aerosol delivery [49,50,51]; however, there is no evidence that the distribution of measles in the powder vaccine is similar to measles virus size distribution in humans-generated bio-aerosols. In this study, because of the lack of a reliable source, we assumed that the size distribution of measles viruses in indoor bio-aerosols is similar to influenza viruses. This assumption is based on the fact that both diseases are airborne viral respiratory infections with similar virus sizes ranging between 80 and 120 nm [52] and 100–200 nm [53] for influenza and measles viruse, respectively.

As we assumed a similar size distribution for measles and influenza viruses, their estimated deposition rate and air filter removal efficiencies will be similar and equal to the reported values in Azimi and Stephens [54]. Azimi and Stephens estimated the deposition rate of influenza viruses to range between 1 and 2.3 per hour and the average removal efficiency (ηfilter) of air filters for particles containing influenza viruses to be between 10.5 and 99.9% for filters with minimum efficiency reporting value (MERV) of 4 and high-efficiency particulate air (HEPA) filters, respectively, as demonstrated in Table 1.

Table 1 Infectious-particle-size-weighted filtration efficiency for a range of HVAC air filters [54]

The particle removal rate due HVAC filtration (Kfiltration,i) for each space was estimated from Eq. 6.

$$ {K}_{filtration,i}=\frac{f_{runtime}\times {f}_{recir}\times {\eta}_{filter}\times {Q}_{return,i}}{V_i} $$
(6)

Qreturn,i: Return air flow rate of HVAC system in space i (m3/hour).

Vi: Volume of space i (m3).

Similarly, the mechanical ventilation rate of HVAC systems (Kventilation,i) in space i can be estimated from Eq. 7.

$$ {K}_{ventilation,i}=\frac{f_{runtime}\times \left(1-{f}_{recir}\right)\times {Q}_{return,i}}{V_i} $$
(7)

Finally, the removal efficiency of measles viruses by an air purifier in space i (Kpurification,i) was estimated from Eq. 8.

$$ {K}_{purification,i}=\frac{f_{runtime, AP,i}\times CAD{R}_{AP,i}}{V_i} $$
(8)

fruntime,AP: Runtime fraction of air purifier in space i.

CADRAP,i: Clean air delivery rate of air purifier in space i (m3/hour).

The CADR of an air purifier is usually estimated by multiplying the air delivery rate and air filter removal efficiency of the air purifier. As most of air purifiers using HEPA filters with removal efficiencies of more than 99% for all particle sizes and types, we assumed the clean air delivery rate of an air purifier, when challenged with bio-aerosols containing measles viruses, is similar to the factory reported CADR of the air purifier.

Estimating the number of susceptible individuals based on vaccination coverage and age

The transmission risk of measles among occupants of an indoor environment depends on the number of susceptible individuals in that cohort, which is estimated using a variety of approaches before, during, and after a measles outbreaks in existing studies. The simplest approach is to assume that everyone is susceptible [3], which is not a realistic approach, particularly for a population with a high vaccination coverage. The most reliable approach is to measure the level of Immunoglobulin G (IgG) antibodies against measles in the occupants’ blood [55, 56]; however, this approach is costly and time-consuming making it hardly available for every measles outbreak. As another approach, Riley et al. suggested to assume that the total number of susceptible individuals is equal to the number of infected cases at the end of an outbreak assuming the outbreak would stop only after all susceptible individuals were infected [33]. This approach by Riley et al. ignores the fact that the number of infected cases and accordingly their estimate of susceptible individuals would change if a different infection control strategy was deployed during the outbreak.

A proper approach for estimating susceptibility in a cohort is based on the vaccination coverage and age of the individuals in that cohort. Several studies have suggested that measles vaccine efficacy is not flawless, and a small portion (i.e., less than 10%) of vaccinated individuals would always remain susceptible to the disease [57,58,59,60]. Landen et al. assumed a 1 and 5% vaccination failure among students receiving two doses and one dose of vaccine, respectively, and 100% susceptibility for non-vaccinated students during the 1996 measles outbreak in Alaska, U.S. [61]. Choi et al. assumed infants less than 6 months old are immune to measles through maternal antibody, cohorts between 6 months and 14 years old are 100, 10 and 1% susceptible if they had received no vaccine, 1 dose of vaccine, and 2 doses of vaccine, respectively, and 5 and 2% of immunized young students between 14 and 24 years old and adults older than 25 years old are susceptible to measles, respectively [62].

In this study, we estimated the percentage of susceptible individuals in the U.S. schools based on age and vaccination coverage similar to Choi et al. [62]. We assumed elementary students less than 14 years old are 100, 10% or 1%, susceptible to measles viruses if they were not vaccinated or had one dose or two doses of the vaccine, respectively, before the outbreak. We also assumed 5% of secondary students between 14 and 18 years old are susceptible to the measles viruses, if they received one or two doses of measles vaccine before an outbreak. We also compare our assumption for number of susceptible people with the reported number of infected cases during actual measles outbreaks in primary and secondary schools in developed countries to verify our assumptions.

Back calculating quanta generation rate

Quanta generation rate (q) in the Wells-Riley model is a critical parameter back-calculated from observational epidemiological studies indicating the contagiousness of an airborne pathogen. It is important to note that for any new derivation of the Wells-Riley model a new associated quanta generation rate should be back-calculated for a desired pathogen, reflecting the assumptions used in developing the new risk model. For example, Riley et al. estimated the quanta generation rate of measles between 480 and 5589 quanta per hour from an outbreak in an elementary school in New York using a steady state Well-Riley model [33]. Later, Rudnick and Milton reported a quanta generation rate of 570 quanta per hour using a set of new assumptions for the risk model [35] and Chen et al. back-calculated a quanta generation rate of 128 quanta per hour using a transient derivation of the risk model [3] for the exact same measles outbreak in the New York elementary school studied firstly by Riley et al.

In this research work, we relied on two well-known studies that have described measles outbreaks in primary and secondary schools in the U.S. to back-calculate the quanta generation rate of measles for our newly developed risk model [22, 33]. These studies were selected because they reported the detailed characteristics of measles outbreaks in the schools. Both studies include information on the school floor plan, HVAC system operation time and characteristics, index case activity patterns, and vaccination records of students before and during the outbreaks. We selected separate case studies for primary and secondary schools because model parameters, such as students’ interaction, susceptibility, and inhalation rate, are varied among students under 14 years old and between 14 and 18 years old attending primary and secondary schools, respectively.

Similar to the developed risk model, we divided the school environment into three spaces or microenvironments including the infector’s classroom, recirculation spaces, and common spaces. Table 2 demonstrates our primary estimates and ranges for the risk model parameters. Most of the model variables were reported directly in Riley et al. and Chen et al. studies; however, some of the parameters were not reported during the outbreaks. In these cases, we considered a range for the model variables and chose a ‘best’ or ‘primary estimate’ for each model parameter reflecting our finest estimations of that variable as shown in Table 2.

Table 2 Summary of outbreak characteristics in primary and secondary representative schools used in quanta generation rate (q) back-calculation process

In Appendix B, the process for culling the model parameters from the case study articles and other listed references is explained in detail and the results are summarized in Table S1 and Figure S1. It is noticeable that because we back-calculated the quanta generation rate from actual epidemiology studies, any simplification that we deployed during the model development and variable estimation will be considered automatically in the quanta generation estimates. We also determined the boundaries of the quanta generation rates for the primary and secondary schools based on the ranges of the model parameters.

Developing a nationwide representative school model for measles transmission risk

Next, we developed a nationwide school model that represents the majority of educational institutions in the U.S. The combination of the nationwide representative School Building Archetype (SBA) model and the developed multi-zone transient Wells-Riley model was used to estimate the range of measles risk that students are facing in U.S. schools. The SBA model considers two sets of parameters for primary and secondary schools investigating the impacts of students’ age and activity patterns on the risk model results. As the majority of elementary schools in the United States have self-contained classrooms [64, 65], we considered students less than 14 years old in elementary (or primary) schools with teacher self-contained educational formats and students between 14 and 18 years old in departmentalized high (or secondary) schools [66]. The estimates of quanta generation rate back-calculated from the teacher self-contained elementary and departmentalized high school case studies in the previous section were used as the only infection-related inputs of the SBA model. Other parameters were either assumed or culled from other references as demonstrated in Table 3. The detailed methodology for developing the SBA model is presented in Appendix C including Tables S2–S4 and Figure S2.

Table 3 Summary of best estimates and ranges of variables used in the nationwide representative School Building Archetype (SBA) model

One thing to notice is, in Table 3, the best estimate values suggested for each model input could also be treated as an input model values for typical school settings. Herein, we considered six typical school settings including three primary and three secondary schools with difference types of HVAC systems and estimated the transmission risk of measles during heating and cooling seasons in those schools. The average risk of measles transmission among these six example sites was selected as our “best estimate” of measles transmission risk in US schools and shown graphs which can be found in the results section graphs.

Similar to the back-calculation process, we chose a best (primary) estimate and a range for most of SBA model variables. We applied a Monte-Carlo simulation with 10,000 iterations to account for the impacts of the model parameter ranges on the transmission risk results. Each iteration represents the risk of measles transmission in one U.S. school setup. For the Monte-Carlo simulation, we culled the model variables from two decks of primary and secondary inputs with the same proportion as the ratio of primary and secondary educational institutions in the U.S. (i.e., 76.6% of iterations were from primary school inputs and 23.3% were from secondary school inputs) [66]. A similar approach was adopted for the ratio of heating and cooling system types in the SBA model. We divided the heating and cooling systems of the U.S. schools into three categories of central-forced-air systems and ductless HVAC systems with and without air filters based on the U.S. Energy Information Administration, Commercial Buildings Energy Consumption Survey [76, 77]. The number of times that we selected each heating and cooling system type in the Monte-Carlo simulation was based on the ratio of the heating and cooling system type in the U.S. schools as summarized in Table 4.

Table 4 Summary of HVAC system types in U.S. schools based on the U.S. Energy Information Administration, Commercial Buildings Energy Consumption Survey

Evaluating the effects of vaccination and airborne infection control strategies on measles transmission risk

Next, we evaluated the effects of proper vaccination (i.e., ≥2 dose vaccine, the first one after 12 months old) and various control strategies related to HVAC systems on measles transmission risk in the U.S. primary and secondary schools. To evaluate the vaccination effectiveness, we compared the measles transmission risk among unvaccinated and vaccinated cohorts in a variety of infection control scenarios in educational institution setups. To investigate the impacts of infection control strategies on the measles transmission risk, we considered three categories of infectious bio-aerosol removal approaches for schools, including improving air filter efficiencies, increasing ventilation rate, and using air purifiers in classrooms as well as their combinations. The effectiveness of the selected control strategies was evaluated by comparing the measles transmission risk after deploying the strategies with the risk in a basic-infection-control scenario of the SBA model (i.e. the removal efficiency of air filters and the ventilation rate were assumed to be equal to the minimum requirements for schools and no air purifier was used in classrooms). For each control strategy category (i.e., air filtration, ventilation, and purification), we assumed a regular and an advanced risk reduction scenario. The regular risk-reduction scenarios are costly affordable approaches compared to the advanced ones and adopted regularly for indoor environments such as schools. The advanced control scenarios are relatively extreme risk-reduction approaches and less common compared to the regular control strategies; however, they are still feasible techniques for decreasing the risk of airborne pathogens in school environments.

Improving removal efficiency of HVAC air filters

In the SBA model, improving the removal efficiency of HVAC air filters is limited to central-forced-air and ductless with air filter heating and cooling systems. EPA’s “Tool for School” program requires all schools to at least have air filters with MERV8 in all HVAC application [82], while the National Air Filtration Association (NAFA) recommends air filters between MERV 8 and 13 for schools [81]. On the other hand, the best commercially available air filters are called HEPA filters, and have aremoval efficiency of 99% and higher for almost all types of aerosols including droplet nuclei containing viable viral pathogens [54, 83]. Herein, for the SBA model with basic control strategies, we assumed the HVAC systems use MERV8 air filters and evaluated the changes in measles transmission risks after adopting MERV13 and HEPA filters in the heating and cooling systems as regular and advanced control scenarios, respectively.

Increasing ventilation rate

Many studies have highlighted the effects of outdoor air ventilation on reducing the transmission risk of measles [84,85,86]. ASHRAE Standard 62.1–2016 requires a default ventilation rate of 6.7 L/s-person for classrooms with students more than 9 years old changing between 4.0 and 9.5 L/s-person in various types of educational facilities [80]. It also obligates default ventilation rates of 4.7 and 5.1 L/s-person for cafeteria and dining rooms, respectively, which are considered as regular common spaces at U.S. schools in this study as explained in Appendix C [80]. For the SBA model with basic infection control designs, we assumed a minimum required ventilation rate of 6.7 L/s-person for infector’s classroom and the recirculation space, and 4.7 L/s-person in common spaces. We assumed double of the required ventilation rates in classrooms (i.e., 13.4 L/s-person) and cafeteria (i.e., 9.4 L/s-person) as the regular ventilation-related control scenario in the modeled schools. For the advanced ventilation-related control scenario, we assumed double of maximum required ventilation rate in educational facilities for the infector’s classroom and recirculation space (i.e., 19.0 L/s-person), and increased the common space ventilation rates to the double of the required ventilation rates for dining rooms (i.e., 10.2 L/s-person). Fisk (2017) summarized the reported ventilation rates in schools from several studies, where the measurements had conducted during occupancy in 20 or more classrooms [87]. The results show the maximum ventilation rate of 21.7 L/s-person in the classrooms, which demonstrates the feasibility of our advanced ventilation-related infection control scenario in schools (19.0 L/s-person in infector’s classroom and recirculation spaces), although it seems relatively excessive.

Using air purifiers in classrooms

Using air purifiers can reduce the transmission risk of viral airborne disease [88]. Herein, we explored the effectiveness of two air purification scenarios for reducing transmission risk of measles in the school representative model. We did not include utilizing air purifiers in classrooms for the SBA model with basic infection control design. Usually, the capacity of an air purifier for removing all particles of a given size is reported by its CADR in units of air volume per time. Although we are not aware of any standard for the size of air purifiers in indoor environments, a rule of thumb is that for every 23.2 m2 (250 square feet) of space about 0.0472 m3/s (100 cubic feet per minute - CFM) of CADR is desirable. Herein, again, we assumed regular and advanced infection control scenarios for using air purifiers only in the modeled schools’ classrooms as the air-purification-related control strategies. Using the rule of thumb, we assumed air purifiers with 0.189 m3/s (400 CFM) CADR for the modeled classrooms as the regular air purification scenario and doubled the CADR to 0.378 m3/s (800 CFM) in the advanced scenario.

In existing studies, several other control strategies have been deployed to reduce the transmission risk of viral airborne diseases in indoor environments including facial mask protection [89, 90], isolation [91, 92], surface disinfection [93, 94], and ultraviolet germicidal irradiation (UVGI) [44, 45, 95,96,97,98], as well as to increase the immunity of individuals by post-exposure prophylaxis after they are exposed to the infection [99]. Exploring the impacts of these control strategies was out of the scope of this study because:

  1. (i)

    using facial mask is not a feasible control approach in schools particularly before the start of an outbreak and it is not considered as a building-related intervention

  2. (ii)

    isolation strategies would be limited to closing schools in under-vaccinated communities during measles outbreaks or asking students to stay at home if they have the disease symptoms, which is not applicable for the scope of this study as herein we only studied the transmission risk of measles between 2 and 4 days before the appearance of the symptoms in the index case or start of an outbreak

  3. (iii)

    we do not expect surface disinfection at the end of a school day while the students would not come back to the school at least for half a day to reduce the measles transmission risk significantly as the primary transmission pathway of measles is airborne and the measles virus usually does not survive more than a few hours outside of a human’s body

  4. (iv)

    UVGI technology in schools can be deployed either by installing in-duct or by upper room UVGI air disinfection systems in the classrooms [100], while the feasibility both approaches are limited. The in-duct UVGI systems are mostly applicable to classrooms with central forced air systems, and cannot be adopted easily for every existing school (we categorized the portable air handing units with ultraviolet lights in the air purification category). Moreover, using the upper room UVGI air disinfection systems for classrooms, particularly as the primary infection control strategy before the outbreak started, increases the risk of overexposure to UV irradiation among students; although it is demonstrated that careful application of upper-room UVGI can be achieved without an apparent increase in the incidence of the most common side effects of accidental UV overexposure in homeless shelters [101]

  5. (v)

    Post-exposure prophylaxis techniques, such as providing vaccination and immune globulin to people who are at risk for severe illness and complications from measles, are not considered as building-related interventions and increases the immunity of individuals against measles infection instead of reducing the transmission risk of measles among students

Results

Estimates of quanta generation rate

Figure 1 demonstrates the quanta generation rate ranges and primary (best) estimates for elementary and high schools in the U.S.

Fig. 1
figure 1

Best estimates (black line inside the boxes) and ranges of quanta generation rate for typical primary (elementary) teacher self-contained and secondary (high) departmentalized schools in the U.S.

Our best estimates (range) of quanta generation rate were 1925 (1185 – 3345) and 2765 (1430 – 5140), for the elementary and high school case studies, respectively. Figure 1 shows our estimates of quanta generation rate in the high school case study were significantly higher than the elementary school, which could be mostly because the high school had a departmentalized education format in which students switched between classes after each break and consequently, the index case had more interaction with the susceptible students. It is also possible that other factors, such as a longer infection period, changes in the virus infectivity in different stages of the disease (the high school index case attended her classes after rash had appeared) and types of measles virus involved in the outbreaks, caused the higher transmission rate of measles among high school students.

Transmission risk of measles in U.S. schools

Figure 2 demonstrates the distributions, ranges and best estimates of measles transmission risk among U.S. students based on their vaccination status. The median (25th and 75th percentiles) measles transmission risk was estimated 4.2% (2.9 and 6.0%), 39.7% (27.2 and 55.6%), and 0.5% (0.3 and 1.0%) among all, susceptible, and properly vaccinated students, respectively, while our best estimates for the same outputs were 3.5, 31.9, and 0.7%, respectively. The difference between our best estimates and the median transmission risk of measles is largely associated with the difference between the mid-range and the best estimates of the SBA model variables summarizes in Table 3.

Fig. 2
figure 2

Distributions, ranges, and best estimates of measles transmission risk among (a) all students with an average proper vaccination coverage of 91% (changes between 90 and 92%), (b) unvaccinated students, and (c) students with proper measles vaccination assuming 1 and 5% of individuals less than 14 and between 14 and 18 years old remain susceptible, respectively

Figure 2 clearly demonstrates the importance of vaccination for reducing the risk of measles transmission in schools. Looking at the 10,000 infection transmission scenarios (iterations) in the U.S. schools, on average, the infection transmission risk among unvaccinated students is 74 times higher than properly vaccinated students with more than two doses of the vaccine. It also shows that the presence of unvaccinated students in U.S. schools increases the risk of new infection cases significantly from less than 1% if all students were vaccinated to the current nationwide best estimate transmission risk of 3.5%.

Figure 3 compares measles transmission risk among all students (i.e., immunized and unvaccinated) for (i) elementary schools with teacher self-contained education systems and departmentalized high schools, (ii) three types of heating and cooling systems including central forced air and ductless with and without air filter systems, and (iii) heating and cooling seasons. The median (25th and 75th percentiles) transmission risk in elementary and high schools were 3.8% (2.7 and 5.4%) and 5.8% (3.9 and 8.3%) with the best risk estimates of 3.1 and 4.6%, respectively. The effect of HVAC system type on the measles transmission risk were lower than the education system changing the median (25th and 75th percentiles) and best risk estimates from 3.7% (2.6 and 5.2%) and 2.9% in ductless with air filter to 4.0% (2.7 and 5.7%) and 3.5% in central forced air to 6.0% (4.5 and 8.0%) and 5.5% in ductless without air filter heating and cooling systems, respectively. The effects of heating and cooling seasons on the transmission risk of measles were insignificant (using Wilcoxon rank sum test) keeping the median and best estimates of measles transmission risk during both seasons about 4.2 and 3.4%, respectively.

Fig. 3
figure 3

Measles transmission risk among all students in (a) primary teacher self-contained versus secondary departmentalized schools and (b) schools with central forced air and ductless with and without air filter heating and cooling systems, (c) schools during cooling and heating seasons. *Transmission risks were significantly different based on nonparametric Wilcoxon rank sum tests with adjusted p-values for the sample size (i.e., \( P=1-{\left(1-0.05\right)}^{1/\sqrt{N_1{N}_2}} \), where N1 and N2 = number of iterations of compared scenarios)

The best estimate of measles transmission risk was an estimated ~ 1.5 times higher in secondary departmentalized schools comparing to primary teacher self-contained ones. It is mostly because of the higher adopted quanta generation rate due to the high-interaction activity patterns of students in departmentalized schools as well as higher susceptibly rate among students between 14 and 18 years old in comparison to students less than 14 years old. The results also demonstrate using air filters in ductless-with-air-filter and central-forced-air systems decreases our best estimates of measles transmission risk ~ 47% and ~ 37% comparing to ductless-without-air-filter systems, respectively. Moreover, the lower infection transmission risk in schools with ductless-with-air-filter heating and cooling systems comparing to the ones with central-forced-air is due to the air recirculation between the infector’s classroom and other spaces. The insignificant difference between transmission risks during heating and cooling seasons demonstrates the variation in types of HVAC system used for heating and cooling in schools (Table 4) does not influence the model results significantly.

Next, we explored the sensitivity of the results to the changes in the SBA model variables. We divided the SBA model variables into three categories of biological-epidemiological, human interaction, and HVAC-and-building-related variables, as demonstrated in Fig. 4. In each diagram, the x- and y-axis values show the changes in the model variables and the transmission risks associated with those changes, respectively. In the diagrams, each branch shows the relative variation of a model parameter from its primary estimate and the associated relative changes in the transmission risk estimates. The lengths of branches were limited to the ranges of the SBA model variables demonstrated in Table 3. The sensitivity analysis approach is described in detail in Appendix D including Figure S3 and Tables S56.

Fig. 4
figure 4

Sensitivity of the measles transmission model in U.S. schools to changes in (a) biological-epidemiological variables, (b) human-interaction-related parameters and (c) HVAC-building-related variables

Figure 4 results show the higher impacts of biological-epidemiological variables on the measles transmission risk estimates compared to the other two variable categories. On average, the biological-epidemiological, HVAC-and-building-related, and human interaction variables alter the estimates of measles transmission risk 59, 17, and 9% compared to our best infection risk prediction, respectively. Among the biological-epidemiological variables, the model outcomes are most sensitive to the range of number of enrolled students. This drastic change in the transmission risk associated with number of enrolled students is actually driven by the variation in the ratio of infected cases in the infector’s classroom versus the number of enrolled students. Next influential variable is quanta generation rate which changes the transmission risk from 37% less to 54% higher than the primary transmission risk estimate when the lowest and highest estimates of the variable are adopted in the SBA model. The range of infection period is the third most impactful model variable among all parameters changing the estimates of infection risk from 29% less to 26% higher than the primary risk estimation. The occupancy density of classrooms has the highest impact on the measles infection risk among human interaction variables varying the risk from 12% lower to 17% higher risks in comparison to the primary risk estimates when the high and low variable estimates are used in the model, respectively. In the HVAC-and-building-related category, removal efficiency of air filters and recirculation rate of classrooms are the two high impactful variables, shifting the transmission risk from 9 and 10% lower to 25 and 21% higher risk estimates comparing to the primary transmission risk scenario, respectively.

Effects of selected infection control strategies on measles transmission

Next, we explored the effects of vaccination, three categories of common infection control strategies suitable for school environments (i.e., enhancing air filtration, ventilation, and purification) and their combinations on transmission risk of the measles virus. For each infection control category, we evaluated the effectiveness of two scenarios including one regular and one advanced infection removal approach. Moreover, we investigated the infection removal effectiveness of four combination scenarios including (i) regular filtration and ventilation, (ii) regular purification and ventilation, (iii) regular filtration, ventilation, and purification, and (iv) advanced filtration, ventilation, and purification improvements in the infection control designs of the SBA model. It is noticeable that the filtration scenarios were only applied to central-forced-air and ductless with air filter systems in the SBA model. Figure 5 compares the transmission risk of measles for various infection control strategies.

Fig. 5
figure 5

Transmission risk of measles among (a) properly immunized and (b) unvaccinated students and the effects of 10 infection control strategies including regular and advanced filtration, ventilation, and air purification (AP) techniques and their combinations

The first and most important implication from Fig. 5 is the importance of vaccination. The results show regardless of the adopted control strategy, the risk of infection transmission among unvaccinated students is significantly higher than the immunized students. It is shown while the median transmission risk of measles among properly immunized students in schools with basic control strategy is estimated 0.6% and for all other control scenarios remains below 0.5%, the median infection risk among unvaccinated students ranges between 46 and 20% for the basic control designs and the combinations of advanced infection removal strategies, respectively. This demonstrates even with the most advanced infection control mechanisms the measles transmission risk among unvaccinated cohorts remains more than 30 times higher than the risk among properly immunized students staying in school environments with basic infection control designs.

Figure 5-b highlights the role of building designs and particularly effects of adopting regular and advanced control strategies and their combinations on reducing the transmission risk of measles in schools among unvaccinated (susceptible) cohorts. Considering central-forced-air and ductless-with-air-filter heating and cooling systems only, upgrading HVAC air filters from MERV8 in the basic control scenario reduced the median (1st and 3rd quartiles) of infection risk among unvaccinated students from 45% (32 and 60%) to 32% (22 and 45%) and 29% (20 and 42%) using MERV13 and HEPA filters, respectively. Enhancing ventilation rates decreased the median (1st and 3rd quartiles) infection risk from 46% (33 and 62%) among unvaccinated students in the basic control scenario (considering all types of HVAC systems) to 38% (26 and 52%) and 33% (23 and 46%) in the regular and advanced infection control scenarios, respectively. It is also shown that deploying air purifiers in classrooms could be more efficient than the explored ventilation scenarios, which reduced the median (1st and 3rd quartiles) infection risk to 37% (26 and 51%) and 31% (22 and 44%) for air purifiers with CADR of 400 CFM and 800 CFM, respectively.

Figure 5 also demonstrates the effects of adopting more than one control strategy at a time on the transmission risk of measles. Deploying two regular control approaches reduced the median (1st and 3rd quartiles) of measles transmission risks among susceptible students to 28% (19 and 40%) and 31% (22 and 44%) when combinations of regular improvements in filtration-ventilation and ventilation-purification approaches are used in the applicable school environments, respectively. Applying filtration, ventilation, and purification techniques together lowered the median (1st and 3rd quartiles) of measles transmission risk to 24% (16 and 34%) and 19% (13 and 28%) for regular and advanced infection control strategies, respectively.

Figure 6 compares the relative effectiveness of the adopted control strategies, which were estimated by comparing the average number of infected cases among all students in the basic infection control designs with the same numbers after enhancing the removal rates of infectious bio-aerosols in the SBA model using different control scenarios. For the control scenarios related to improving the effectiveness of air filters, we only considered the central-forced-air and ductless-with-air-filter heating and cooling systems.

Fig. 6
figure 6

Relative effectiveness of advanced control strategies on measles transmission risk among all students

Figure 6 results show for the modeled heating and cooling systems, upgrading the air filters to MERV13 and HEPA filters reduce the average number of infected students by approximately 28 and 33%, respectively. Ventilation-related control strategies had average effectiveness of 18 and 28% for regular and advanced scenarios, respectively. Placing air purifiers with regular CADR of 400 cfm in typical school classrooms decreased the number of infected cases by 18%, while doubling the CADR of the air purifiers in the advanced control scenario to 800 cfm increased the effectiveness of the control method to 31%. Using two regular infection control scenarios increased the effectiveness of the control strategies to 37 and 30% when combinations of filtration-ventilation and ventilation-purification scenarios were adopted, respectively, showing the combinations of two regular control approaches can be as effective as adopting one advanced control strategy. Combining all regular and advanced control scenarios reduced the average number of infected cases up to 45 and 56%, respectively, demonstrating the potentially high impacts of building designs on avoiding the new cases of airborne disease infections in school setups.

Discussion

We used a combination of a newly developed multi-zone Wells-Riley approach, a nationwide representative School Building Archetype (SBA) model, and a Monte-Carlo simulation to estimate the transmission risk of measles among U.S. students. In the multi-zone Wells-Riley model, we considered several microenvironments within a typical school building (i.e. infector’s classroom, recirculation space, and common space) and simulated the transmission of measles virus between the zones based on the building and epidemiological characteristics of schools adopted from the SBA model. In the SBA model, we considered two education formats (i.e., teacher self-contained and departmentalized schools) for U.S. schools affecting the interaction of susceptible students with index cases, three categories of HVAC systems (i.e., central-forced-air and ductless with and without air filter systems), and changes in susceptibility of students based on their age and vaccination status. Finally, we investigated the effectiveness of vaccination and ten control strategy scenarios related to enhancing air filtration, ventilation, and purification rates in schools for reducing the risk of measles transmission among primary and secondary students.

As expected, vaccination was shown to be the primary approach for reducing the transmission risk of measles among students. However, the risk still exists when a contagious kid in school encounters vaccinated individuals; therefore, further motivation for our work was understanding the role of factors that influence disease transmission beyond vaccination. Here, we found school educational formats and building and HVAC characteristics play critical roles in measles transmission. Specifically, we found that the transmission risk of measles in primary schools (assuming teacher self-contained classrooms) is less than secondary schools (assuming departmentalized systems) and schools with ductless-with-air-filter and ductless-without-air-filter systems have the lowest and highest transmission risks of measles, respectively.

Comparing the risk model assumptions and results with existing epidemiological studies

In this section, we evaluated our assumptions for the susceptibility of students to measles virus based on their age and vaccination coverage and our best estimates of measles transmission risk among susceptible students as the primary outcome of this study.

We assumed students less than 14 years old are 100, 10 and 1%, susceptible to measles if they were not vaccinated, had one dose of measles vaccine, and had two doses of the vaccine, respectively, and students between 14 and 18 years old are 100 and 5% susceptible to measles viruses if they were not vaccinated and had either one or two doses of the vaccine, respectively. In Tables 5 and 6, we compared our assumptions for susceptibility of students with the reported number of infected cases and measles attack rates during 27 measles outbreaks in primary and secondary schools in developed countries. The selected studies provided information on the total number of enrolled students during the outbreak, final number of infected cases, and vaccination coverage of the students. In all cases, one index case started the outbreak and infected at least one other susceptible student at their school. During the outbreaks if a portion of students were vaccinated or revaccinated, we considered the final vaccination coverage of students in the analysis. We also culled information regarding the number of received vaccination doses, attack rates among unvaccinated (ARU) and immunized (ARI) individuals, and number of infected cases in the first generation of the outbreaks when the data were provided. For validation purposes, we expect the reported measles attack rates among immunized individuals with one dose of vaccination (ARI-1-Dose) and two doses of vaccination (ARI-2-Dose) remains lower than the assumed susceptibility rates for the cohorts and the estimated number of susceptible individuals to be larger than the final number of infected cases during the outbreaks.

Table 5 Characteristics of measles outbreaks in primary schools
Table 6 Characteristics of measles outbreaks in secondary schools

Table 5 summarizes the characteristics of 10 measles outbreaks in the U.S. and 5 outbreaks in other developed countries in primary schools with the average age of students less than 14 years old. Seven of the summarized studies in Table 5 did not provide information about the portion of students who had received one or two doses of measles vaccine. In these cases, we assumed the immunized students had received one dose of the measles vaccine. In all cases, the number of estimated susceptible individuals was higher than the final number of infected cases at the end of the outbreaks.

The maximum ARU in primary schools listed in Table 5 was 84% during an outbreak in New York State, between 1945 and 1946, although the school had deployed a UVGI control system to prevent the spread of measles. The high attack rate of measles during this outbreak and the other two outbreaks in New York State reported in the study by Perkins et al. [44] shows the validity of our assumption for 100% susceptibility to measles among unvaccinated cohorts. The maximum ARI was 6.4% reported in the study by Riley et al. [33], where the ARI-1-Dose and ARI-2-Dose values were not mentioned. Considering the fact that the outbreak happened in 1974 when 2-dose vaccination was not common and assuming immunized students, in this case, had received one dose of vaccination, the maximum ARI in primary schools was in line with our 10% assumption of susceptibility to measles among students under 14 years old who only receive a 1-dose vaccine. The ARI-2-Dose values summarized in Table 5 were also lower than our assumption of 1% susceptibility between primary school students except for one Texas junior high school case, where the attack rate was 1.4%.

Table 6 demonstrates the characteristics of 12 measles outbreaks (9 in the U.S. and 3 in other countries) in secondary schools with the average age of students between 14 and 18 years old. Most of the summarized articles (10 out of 12 studies) reported the proportion of students who had received one or two doses of measles vaccine as well as the measles attack rate during the studied outbreaks. The measles attack rate in all cases remained lower than our assumption of the vaccine failure among students between 14 and 18 years old (i.e. 5%) except in one extreme outbreak in Finland where 29 and 14.3% of students who had received one dose and two doses of measles vaccine, respectively, were infected. It is noticeable that five of the ARI-2-Dose values were between 1 and 5%, demonstrating our assumption for higher vaccination failure among students between 14 and 18 years old receiving two doses of measles vaccine comparing to students less than 14 years old was valid.

We also estimated the number of avoided infected cases by reducing the total number of infected cases from the number of susceptible students. The results show that other control mechanisms than the vaccination such as air filtration, ventilation, and purification and isolation of suspicious infected cases during the outbreaks, on average, could protect ~ 60% (ranged between 11 and 98%) of susceptible students from the infection in the studied outbreaks. This demonstrates the importance of paying attention to other control strategies in addition to the vaccination to reduce the transmission risk of measles in the built environment.

The primary outcome of the developed Wells-Riley model is the infection risk defined as the number of infected cases divided by the number of susceptible individuals during one generation of the infection outbreak as demonstrated in Eq. 1. Therefore, to evaluate our developed model results, we compared our best estimate of nationwide infection risk of measles in U.S. schools with the estimated transmission rate of measles during the first generation of the infection outbreaks in schools in developed countriesreported in existing epidemiological studies. Totally, the number of infected cases during the first generation of measles outbreaks was reported in 24 summarized studies in Tables 5 and 6. The measles transmission rates were estimated by dividing the reported number of infected cases during the first generation of the outbreaks to the estimated number of susceptible students in each school as shown in Fig. 7.

Fig. 7
figure 7

Comparing our best estimate of nationwide measles transmission risk among susceptible students in U.S. schools with estimated transmission rates of measles during first generations of the infection outbreaks in developed countries’ schools among susceptible students reported in existing epidemiological studies

Most of the demonstrated studies in Fig. 7 did not directly report the number of infected cases during the first generation of the outbreaks; instead, they reported the timeline when the infected cases were detected. In these cases, we considered a 14-day (± 7 days) incubation period (i.e., the time elapsed between exposure to a pathogenic organism, and when symptoms and signs are first apparent) for measles, which means we considered the infected cases in a 7-to-21 time period after the index cases had entered the schools as individuals who were infected during the first wave of the outbreaks.

Figure 7 shows an average (± SD) first-generation measles transmission rate of 21% (±18%) for the reported outbreaks in developed countries’ primary and secondary schools, ranging between 3 and 69%, while our best nationwide transmission risk estimate is 32% among susceptible students. We believe the higher estimate of transmission risk derived from the developed Wells-Riley model does not necessarily mean that the developed model overestimates the infection risks as long as our risk estimates are within the range (i.e., average ± SD) of reported first-generation transmission rates, because:

  1. (i)

    The number of selected outbreaks are limited and they are not representative of U.S. schools

  2. (ii)

    Except the two studies used for the quanta generation rate back-calculation process, other selected studies have not provided enough information regarding the student activities, school building properties, and epidemiological characteristics of the outbreaks; therefore, we are unaware of the infection control strategies deployed in most of these schools during the outbreaks

For these reasons, the fact that the primary outcome (i.e. transmission risk among susceptible students) of our risk model is within the range (average ± SD) of the reported measles transmission rates found in the existing literature demonstrates the validity of our models adopted in this study.

We also demonstrated the estimated transmission risk of measles in six typical US school settings among susceptible students (based on our best estimates of school building characteristics presented in Table 3) in Appendix E (Figure S4) and compared them with estimated transmission rates of measles during the first generations of the infection outbreaks in schools from developed countries among susceptible students reported in existing epidemiological studies.

Implications

The results of the current simulation study indicate the primary importance of vaccination for reducing the risk of measles transmission among students at schools. Additionally, our results related to the estimated distribution and range of measles transmission risk in school can aid epidemiologists and risk analyzers to evaluate the chance of a new outbreak in a community. Moreover, the study outcomes shown in Figs. 5 and 6 clearly demonstrate that beyond vaccination, several factors such as increasing filtration, ventilation and air purification rates in indoor environments also were influential in disease transmission. However, none of these interventions were as effective as vaccination and should not be used as a basis or control strategy in place of vaccination. Our findings support their use as a supplemental control strategy that must be combined with vaccination.

In this study, for the first time, we developed a nationwide representative School Building Archetype (SBA) model and a transient multi-zone Wells-Riley model for estimating the transmission risk of an airborne infectious agent (i.e., measles viruses) among U.S. students. We also demonstrated that the combination of the SBA and transient multi-zone Wells-Riley models estimates the nationwide infection risk of measles within the range (i.e., average ± SD) of first-generation transmission rates of measles in schools according to the existing epidemiological studies (Fig. 7). As the only three biological-related variable were involved with the SBA and Wells-Riley models (i.e., quanta generation rate, infection period, and deposition rate), the newly developed models are also capable of estimating the nationwide transmission risk of other airborne infectious disease in school environments as long as the biological-related variables of the desired airborne pathogen are available. Moreover, the methods used in this study can be expanded to other indoor environments such as offices, healthcare facilities, and residences. Therefore, policy makers and standard developers can adopt the methodology used in this study to establish new nationwide and regional policies and requirements for school or other indoor environments to reduce the chance of spread of infectious airborne disease.

We also evaluated the impacts of HVAC system designs, educational format of schools, and several control strategies on the transmission risk of measles in schools in Figs. 3 and 5 and compared the effectiveness of a variety of infection control approaches on reducing the average number of infected cases in the SBA model in Fig. 6. Expanding the framework of this study to other airborne infectious diseases, would provide additional information for building designers and decision-makers to consider before selecting the most appropriate HVAC system types for school buildings. Moreover, the method used in this study is not limited to nationwide estimates of the airborne infection transmission risks, but it can also be deployed by building designers to predict the transmission risk of a variety of infectious diseases in a specific school environment during the design process. For example, as demonstrated in Fig. 3, the transmission risk of measles in schools with ductless HVAC systems and a proper air filter is lower than the other two types of HVAC systems, which most probably would be the same condition for other airborne infectious diseases. As another implication, the comparison between the effectiveness of various control strategies will help school officials to select a financially appropriate infection control approach based on the school’s building and HVAC system characteristics to reduce the transmission risk of measles or other infectious airborne diseases in an existing school building. For example, results shown in Fig. 6 demonstrate adopting HEPA filters instead of MERV13 in HVAC systems would improve the filtration effectiveness less than 5%, while existing studies estimated annual costs of HEPA filters are more than double of MERV13 air filters [54] or if advanced infection control strategies cannot be deployed for a school because of high installation costs or building and HVAC system properties, adopting a combination of two regular infection control scenarios can provide a similar or even higher removal rates of infectious bio-aerosols.

Herein, for the first time we back-calculated quanta generation rate for one infectious disease from multiple studies. This approach helped us to capture a wider potential range for measles quanta generation rate considering different school setups (Fig. 1), which canalso be adopted by researchers in future studies to back-calculate the quanta generation rate of other airborne pathogens. Moreover, the sensitivity analysis (Fig. 4) determined what model variables have the highest impacts on the measles transmission risk estimates in schools. The outcomes of the sensitivity analysis help researchers to identify the critical parameters in the risk model and highlights the most influential research pathways for future studies.

Limitations

One of the most challenging parts of this simulation effort was to find proper ranges for the SBA model variables representing the majority of school building stock in the U.S. The variable ranges were particularly essential for parameters that have high impacts on the risk results including quanta generation rate, infection period, density of students in infector’s classroom, and classrooms’ recirculation and filtration rates. However, the existing knowledge around the ranges of many of the model variables is limited, particularly for quanta generation rate and air recirculation rate in schools. For the quanta generation rate, we only found two studies describing the characteristics of measles outbreaks in schools from 1974 and 1989 [22, 33] in such details that could be used in the back-calculation process; and for recirculation rate, we relied on one study [78] measuring the average air recirculation rate in nine California classrooms. Although the study by Riley et al. is the most well known article used for estimating the range of quanta generation rate, and measurments in the study byPolidori et al. have been used as the representative of U.S. classrooms’ recirculation rate in other peer-reviewed articles [79], more comprehensive studies are required for evaluating the ranges of these variables.

The Wells-Riley approach only considers the airborne transmission pathways of infectious diseases while, the transmission risk of measles through other pathways including fomite and direct contact remains unclear. Although, it is shown that measles primary transmission pathway is airborne, it is necessary to perform more research on the other pathways of measles transmission. This limitation would be more critical for other airborne infectious diseases such as influenza and coronavirus which fomite and direct contact are also shown to have a considerable influence on the risk results.

In developing the transient Wells-Riley model, we made several simplifications such as assuming continuous stay of students in the microenvironments, constant number of students, and a simplified format of student interactions. Although these factors were considered in some levels during the back-calculation process, for individual case studies (not nationwide simulations), where more information on building and human interaction characteristics is available, more advanced and complex derivations of the Wells-Riley approach or other mathematical and statistical models should be deployed.

Conclusion

We used a combination of a newly developed transient multi-zone Wells-Riley approach, a nationwide representative School Building Archetype (SBA) model, and a Monte-Carlo simulation to estimate the transmission risk of measles among students in U.S. schools. We also estimated the number of susceptible students for a school setup based on the age and vaccination record of students and back-calculated quanta generation rate of measles for our newly developed risk model based on two existing epidemiological studies. We considered three microenvironments within school buildings, two education formats, and three types of HVAC systems in the risk model and used the Monte-Carlo simulation with 10,000 iterations to examine the effects of model parameter ranges on the risk results. Our best estimates of nationwide transmission risk of measles in U.S. school were 3.5 and 32% among all and susceptible students, respectively. The results of our study show the transmission risk of measles among unvaccinated students is more than 70 times higher than properly immunized ones. In the back-calculation process, we estimated the quanta generation rate of 2765 and 1925 quanta per hour for primary schools with teacher self-contained classrooms and secondary schools with departmentalized system, respectively, showing that higher student interactions in the departmentalized schools significantly increases the transmission risk of measles. Comparing various types of HVAC systems shows schools with ductless-with-air-filter systems have the lowest transmission risk of measles, while the risk is highest for schools with ductless-without-air-filter systems. Finally, exploring the effectiveness of 10 control scenarios for reducing the transmission risk of measles in schools shows a large difference among the effectiveness of various control strategies and the selected infection control approaches can reduce the average number of infected cases up to 56% when a combination of advanced air filtration, ventilation, and purification approaches was adopted.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.

Change history

  • 24 July 2020

    An amendment to this paper has been published and can be accessed via the original article.

Abbreviations

AP:

Air Purifier

ARI:

Attack Rates among Immunized

ARI-1-Dose:

Attack Rates among Immunized individuals with one Dose of vaccination

ARI-2-Dose:

Attack Rates among Immunized individuals with two Dose of vaccination

ARU:

Attack Rates among Unvaccinated individuals

ASHRAE:

American Society of Heating, Refrigerating and Air-Conditioning Engineers

CADR:

Clean Air Delivery Rate

CDC:

Centers for Disease Control

CFD:

Computational Fluid Dynamics

CFM:

Cubic Feet per Minute

DOE:

Department of Energy

E-M-HS:

Elementary, Middle, and High School

EPA:

Environmental Protection Agency

ES:

Elementary School

HEPA:

High-Efficiency Particulate Air

HS:

High School

HVAC:

Heating, ventilation, and air-conditioning

ID63:

63% Infectious Dose

IgG:

Immunoglobulin G

JHS:

Junior High School

MERV:

Minimum Efficiency Reporting Value

MS:

Middle School

NAFA:

National Air Filtration Association

NCES:

National Center for Education Statistics

NR:

Not Reported

PS:

Primary School

SASS:

Schools and Staffing Survey

SBA:

School Building Archetype

SBS:

Senior Boarding School

SD:

Standard Deviation

SHS:

Senior High School

SIR:

Susceptible-Infectious-Recovered

SIS:

Susceptible-Infector-Susceptible

UVGI:

Ultraviolet Germicidal Irradiation

References

  1. Sugerman DE, Barskey AE, Delea MG, Ortega-Sanchez IR, Bi D, Ralston KJ, et al. Measles outbreak in a Highly vaccinated population, San Diego, 2008: role of the intentionally Undervaccinated. Pediatrics. 2010 Apr 1;125(4):747–55.

    PubMed  Google Scholar 

  2. WHO. Measles . [cited 2019 Oct 1]. Available from: https://www.who.int/news-room/fact-sheets/detail/measles.

  3. Chen S-C, Chang C-F, Liao C-M. Predictive models of control strategies involved in containing indoor airborne infections. Indoor Air. 2006 Dec;16(6):469–81.

    CAS  PubMed  Google Scholar 

  4. Fisk WJ. Health and productivity gains from better indoor environments and their implications for the U.S. Department of Energy. CA (US): Lawrence Berkeley National Lab; 2000. [cited 2019 Aug 29]. Report No.: LBNL-47458. Available from: https://www.osti.gov/biblio/780590.

    Google Scholar 

  5. Azor-Martínez E, Gonzalez-Jimenez Y, Seijas-Vazquez ML, Cobos-Carrascosa E, Santisteban-Martínez J, Martínez-López JM, et al. The impact of common infections on school absenteeism during an academic year. Am J Infect Control. 2014 Jun 1;42(6):632–7.

    PubMed  Google Scholar 

  6. Öhlund LS, Ericsson KB. Elementary school achievement and absence due to illness. J Genet Psychol. 1994;155(4):409–21.

    PubMed  Google Scholar 

  7. Pourabbasi A, Shirvani ME, Khashayar P. Sickness absenteeism rate in Iranian schools during the 2009 epidemic of type a influenza. J Sch Nurs. 2012 Feb 1;28(1):64–9.

    PubMed  Google Scholar 

  8. Measles | History of Measles | CDC [Internet]. 2019 [cited 2019 Jul 30]. Available from: https://www.cdc.gov/measles/about/history.html.

  9. Measles --- United States, 2000 [Internet]. [cited 2019 Oct 2]. Available from: https://www.cdc.gov/mmwr/preview/mmwrhtml/mm5106a2.htm.

  10. MMWR Summary of Notifiable Diseases, United States, 1993. [cited 2019 Oct 2]. Available from: https://www.cdc.gov/mmwr/preview/mmwrhtml/00035381.htm.

    Google Scholar 

  11. CDC. Measles: Answers to Common Questions. Centers for Disease Control and Prevention. 2019 [cited 2019 Aug 8]. Available from: https://www.cdc.gov/measles/about/faqs.html.

    Google Scholar 

  12. CDC. Measles Cases and Outbreaks. Centers for Disease Control and Prevention. 2020 [cited 2020 Jan 22]. Available from: https://www.cdc.gov/measles/cases-outbreaks.html.

    Google Scholar 

  13. Cisse B, Aaby P, Simondon F, Samb B, Soumaré M, Whittle H. Role of schools in the transmission of measles in rural SenegalImplications for measles control in developing countries. Am J Epidemiol. 1999;149(4):295–301.

    CAS  PubMed  Google Scholar 

  14. Orenstein WA, Hinman AR. The immunization system in the United States — the role of school immunization laws. Vaccine. 1999 Oct 29;17:S19–24.

    PubMed  Google Scholar 

  15. van den Hof S, van Spaendonck MAE C, van Steenbergen JE. Measles Epidemic in The Netherlands, 1999–2000. J Infect Dis. 2002;186(10):1483–6.

    PubMed  Google Scholar 

  16. Ferreira SC, Castellano C, Pastor-Satorras R. Epidemic thresholds of the susceptible-infected-susceptible model on networks: a comparison of numerical and theoretical results. Phys Rev E. 2012 Oct 15;86(4):041125.

    Google Scholar 

  17. Liao C-M, Chang C-F, Liang H-M. A probabilistic transmission dynamic model to assess indoor airborne infection risks. Risk Anal. 2005;25(5):1097–107.

    PubMed  Google Scholar 

  18. Brookmeyer R, Johnson E, Bollinger R. Modeling the optimum duration of antibiotic prophylaxis in an anthrax outbreak. Proc Natl Acad Sci. 2003;100(17):10129–32.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. Brookmeyer R, Johnson E, Barry S. Modelling the incubation period of anthrax. Stat Med. 2005;24(4):531–42.

    PubMed  Google Scholar 

  20. Ellner SP, Bailey BA, Bobashev GV, Gallant AR, Grenfell BT, Nychka DW. Noise and nonlinearity in measles epidemics: combining mechanistic and statistical approaches to population modeling. Am Nat. 1998 May 1;151(5):425–40.

    CAS  PubMed  Google Scholar 

  21. Fine PEM, Clarkson JA. Measles in England and Wales—I: an analysis of factors underlying seasonal patterns. Int J Epidemiol. 1982 Mar 1;11(1):5–14.

    CAS  PubMed  Google Scholar 

  22. Chen RT, Goldbaum GM, Wassilak SGF, Markowitz LE, Orenstein WA. AN EXPLOSIVE POINT-SOURCE MEASLES OUTBREAK IN A HIGHLY VACCINATED POPULATIONMODES OF TRANSMISSION AND RISK FACTORS FOR DISEASE. Am J Epidemiol. 1989;129(1):173–82.

    CAS  PubMed  Google Scholar 

  23. Fennelly KP, Davidow AL, Miller SL, Connell N, Ellner JJ. Airborne infection with Bacillus anthracis—from Mills to mail. Emerg Infect Dis. 2004 Jun;10(6):996–1001.

    PubMed  PubMed Central  Google Scholar 

  24. Kubica GP. Exposure risk and prevention of aerial transmission of tuberculosis in health care settings. Indoor Air Hum Health. 1996;141–9.

  25. Nicas M, Jones RM. Relative contributions of four exposure pathways to influenza infection risk. Risk Anal. 2009;29(9):1292–303.

    PubMed  Google Scholar 

  26. Nicas M, Sun G. An integrated model of infection risk in a health-care environment. Risk Anal. 2006;26(4):1085–96.

    PubMed  Google Scholar 

  27. Jones RM, Adida E. Influenza infection risk and predominate exposure route: uncertainty analysis. Risk Anal. 2011;31(10):1622–31.

    PubMed  Google Scholar 

  28. Chen C, Liu W, Lin C-H, Chen Q. A Markov chain model for predicting transient particle transport in enclosed environments. Build Environ. 2015 Aug 1;90:30–6.

    PubMed  PubMed Central  Google Scholar 

  29. King M-F, Noakes CJ, Sleigh PA. Modeling environmental contamination in hospital single- and four-bed rooms. Indoor Air. 2015;25(6):694–707.

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Azimi P. Advancing knowledge of indoor aerosol sources, fate, transport, and control. Illinois Institute of Technology; 2016.

    Google Scholar 

  31. Noakes Catherine J, Andrew SP. Mathematical models for assessing the role of airflow on the risk of airborne infection in hospital wards. J R Soc Interface. 2009;6(suppl_6):S791–800.

    PubMed  PubMed Central  Google Scholar 

  32. Wells MW. VENTILATION IN THE SPREAD OF CHICKENPOX AND MEASLES WITHIN SCHOOL ROOMS. J Am Med Assoc. 1945;129(3):197–200.

    Google Scholar 

  33. Riley EC, Murphy G, Riley RL. AIRBORNE SPREAD OF MEASLES IN A SUBURBAN ELEMENTARY SCHOOL. Am J Epidemiol. 1978;107(5):421–32.

    CAS  PubMed  Google Scholar 

  34. GNS T, Chao CYH. Review and comparison between the Wells–Riley and dose-response approaches to risk assessment of infectious respiratory diseases. Indoor Air. 2010;20(1):2–16.

    Google Scholar 

  35. Rudnick SN, Milton DK. Risk of indoor airborne infection transmission estimated from carbon dioxide concentration. Indoor Air. 2003 Sep 1;13(3):237–45.

    CAS  PubMed  Google Scholar 

  36. Pereira ML, Vilain R, Tribess A, Morawska L. Risk assessment for airborne infectious diseases between natural ventilation and split-system air conditioner in a university classroom. In: Proceedings of the 23rd ABCM International Congress of Mechanical Engineering, COBEM 2015. Associacao Brasileira de Engenharia e Ciencias Mecanicas-ABCM; 2015. p. 1–5.

    Google Scholar 

  37. Noakes CJ, Sleigh PA. Applying the Wells-Riley equation to the risk of airborne infection in hospital environments: The importance of stochastic and proximity effects. In: Indoor Air 2008: The 11th International Conference on Indoor Air Quality and Cl. Leeds; 2008.

    Google Scholar 

  38. Wu Y, Tung TCW, Niu J. On-site measurement of tracer gas transmission between horizontal adjacent flats in residential building and cross-infection risk assessment. Build Environ. 2016 Apr 1;99:13–21.

    PubMed  PubMed Central  Google Scholar 

  39. Watson JC, Hadler SC, Dykewicz CA, Reef S, Phillips L. Measles, Mumps, and Rubella - Vaccine Use and Strategies for Elimination of Measles, Rubella, and Congenital Rubella Syndrome and Control of Mumps: Recommendations of the Advisory Committee on Immunization Practices (ACIP). Vol. 47/No. RR-8. ATLANTA: CENTERS FOR DISEASE CONTROL; 1998. [cited 2019 Sep 23]. Available from: https://apps.dtic.mil/docs/citations/ADA346849.

    Google Scholar 

  40. Cutts FT, Markowitz LE. Successes And Failures In Measles Control. J Infect Dis. 1994;170(Supplement_1):S32–41.

    PubMed  Google Scholar 

  41. Cutts FT, Henderson RH, Clements CJ, Chen RT, Patriarca PA. Principles of measles control. Bull World Health Organ. 1991;69(1):1–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  42. About Measles Vaccination | Vaccination and Immunizations | CDC. 2019 [cited 2019 Oct 4]. Available from: https://www.cdc.gov/vaccines/vpd/measles/index.html.

    Google Scholar 

  43. Remington PL, Hall WN, Davis IH, Herald A, Gunn RA. Airborne transmission of measles in a Physician’s office. JAMA. 1985 Mar 15;253(11):1574–7.

    CAS  PubMed  Google Scholar 

  44. Perkins JE, Bahlke AM, Silverman HF. Effect of ultra-violet irradiation of classrooms on spread of measles in large rural central schools preliminary report. Am J Public Health Nations Health. 1947;37(5):529–37.

    PubMed  PubMed Central  Google Scholar 

  45. Cutchins EC, Dayhuff TR. Photoinactivation of measles virus. Virology. 1962 Jul 1;17(3):420–5.

    CAS  PubMed  Google Scholar 

  46. Biellik RJ, Clements CJ. Strategies for minimizing nosocomial measles transmission. World Health Organ Bull World Health Organ Geneva. 1997;75(4):367–75.

    CAS  Google Scholar 

  47. Davison AM. Pathogen inactivation and filtration efficacy of a new anti-microbial and anti-viral surgical facemask and N 95 against dentistry-associated microorganisms; 2012.

    Google Scholar 

  48. Deru M, Field K, Studer D, Benne K, Griffith B, Torcellini P, et al. U.S. Department of Energy commercial reference building models of the national building stock. Publ E. 2011;1:1–118.

    Google Scholar 

  49. LiCalsi C, Maniaci MJ, Christensen T, Phillips E, Ward GH, Witham C. A powder formulation of measles vaccine for aerosol delivery. Vaccine. 2001;19(17):2629–36.

    CAS  PubMed  Google Scholar 

  50. de Swart RL, LiCalsi C, Quirk AV, van Amerongen G, Nodelman V, Alcock R, et al. Measles vaccination of macaques by dry powder inhalation. Vaccine. 2007 Jan 26;25(7):1183–90.

    PubMed  Google Scholar 

  51. Kisich KO, Higgins MP, Park I, Cape SP, Lindsay L, Bennett DJ, et al. Dry powder measles vaccine: particle deposition, virus replication, and immune response in cotton rats following inhalation. Vaccine. 2011 Jan 29;29(5):905–12.

    CAS  PubMed  Google Scholar 

  52. Stanley WM. THE SIZE OF INFLUENZA VIRUS. J Exp Med. 1944;79(3):267–83.

    CAS  PubMed  PubMed Central  Google Scholar 

  53. Benyesh M, Pollard EC, Opton EM, Black FL, Bellamy WD, Melnick JL. Size and structure of echo, poliomyelitis, and measles viruses determined by ionizing radiation and ultrafiltration. Virology. 1958 Apr 1;5(2):256–74.

    CAS  PubMed  Google Scholar 

  54. Azimi P, Stephens B. HVAC filtration for controlling infectious airborne disease transmission in indoor environments: predicting risk reductions and operational costs. Build Environ. 2013 Dec 1;70:150–60.

    PubMed  PubMed Central  Google Scholar 

  55. Damien B, Huiss S, Schneider F, Muller CP. Estimated susceptibility to asymptomatic secondary immune response against measles in late convalescent and vaccinated persons. J Med Virol. 1998;56(1):85–90.

    CAS  PubMed  Google Scholar 

  56. Sheikh S, Ali A, Zaidi AKM, Agha A, Khowaja A, Allana S, et al. Measles susceptibility in children in Karachi, Pakistan. Vaccine. 2011;29(18):3419–23.

    PubMed  PubMed Central  Google Scholar 

  57. Brunell PA, Weigle K, Murphy MD, Shehab Z, Cobb E. Antibody response following measles-mumps-rubella vaccine under conditions of customary use. JAMA. 1983 Sep 16;250(11):1409–12.

    CAS  PubMed  Google Scholar 

  58. Hayden GF. Clinical review : measles vaccine failure: a survey of causes and means of prevention. Clin Pediatr (Phila). 1979;18(3):155–6.

    CAS  Google Scholar 

  59. Nkowane BM, Bart SW, Orenstein WA, Baltier M. Measles outbreak in a vaccinated school population: epidemiology, chains of transmission and the role of vaccine failures. Am J Public Health. 1987 Apr 1;77(4):434–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  60. Weibel RE, Buynak EB, McLean AA, Hilleman MR. Follow-up surveillance for antibody in human subjects following live attenuated measles, mumps, and rubella virus vaccines. Proc Soc Exp Biol Med. 1979;162(2):328–32.

    CAS  PubMed  Google Scholar 

  61. Landen MG, Beller M, Funk E, Rolka H, Middaugh J. Measles outbreak in Juneau, Alaska, 1996: implications for future outbreak control strategies. Pediatrics. 1998;102(6):E71.

    CAS  PubMed  Google Scholar 

  62. Choi YH, Gay N, Fraser G, Ramsay M. The potential for measles transmission in England. BMC Public Health. 2008 Sep 26;8(1):338.

    PubMed  PubMed Central  Google Scholar 

  63. EPA U. Exposure factors handbook: 2011 edition. Intake Fish Shellfish Natl Cent Environ Assess. 2011;10:66.

    Google Scholar 

  64. Chan TC, Jarman D. Departmentalize elementary schools. Principal. 2004;84(1):70–2.

    Google Scholar 

  65. Strohl A, Schmertzing L, Schmertzing R, Hsiao E. Comparison of self-contained and departmentalized elementary teachers’ perceptions of classroom structure and job satisfaction. J Stud Educ. 2014;4(1):109–27.

    Google Scholar 

  66. Snyder TD, Brey C de, Dillow SA. Digest of Education Statistics 2017 53rd Edition. U.S. Department of Education, National Center for Education Statistics; 2019. Report No.: NCES 2018–070.

    Google Scholar 

  67. Chen C-S. Numbers and Types of Public Elementary and Secondary Schools From the Common Core of Data: School Year 2009–10 First Look. Washington D.C.: U.S. Department of Education, National Center for Education Statistics; 2012. Report No.: NCES 2011-345rev.

  68. Perry RT, Halsey NA. The Clinical Significance of Measles: A Review. J Infect Dis. 2004;189(Supplement_1):S4–16.

    PubMed  Google Scholar 

  69. Moss WJ. Measles. Lancet. 2017 Dec 2;390(10111):2490–502.

    PubMed  Google Scholar 

  70. Moss WJ, Griffin DE. Global measles elimination. Nat Rev Microbiol. 2006 Dec;4(12):900–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  71. Walker TY, Elam-Evans LD, Yankey D, Markowitz LE, Williams CL, Mbaeyi SA, et al. National, regional, state, and selected local area vaccination coverage among adolescents aged 13–17 years — United States, 2017. Morb Mortal Wkly Rep. 2018 Aug 24;67(33):909–17.

    Google Scholar 

  72. Hill HA. Vaccination Coverage Among Children Aged 19–35 Months — United States, 2017. MMWR Morb Mortal Wkly Rep. 2018;67 [cited 2019 Aug 7]. Available from: https://www.cdc.gov/mmwr/volumes/67/wr/mm6740a4.htm.

  73. Schools and Staffing Survey (SASS) [Internet]. [cited 2019 Sep 18]. Available from: https://nces.ed.gov/surveys/sass/tables/sass1112_2013312_s2s_007.asp.

  74. Schools and Staffing Survey (SASS) [Internet]. [cited 2019 Sep 18]. Available from: https://nces.ed.gov/surveys/sass/tables/sass0708_035_s1s.asp.

  75. Conklin MT, Lambert LG. Eating at school: A summary of NFSMI research on time required by students to eat lunch. University: University of Mississippi, National Food Service Management Institute. 2001. Retrieved December 12, 2005.

  76. Winiarski DW, Jiang W, Halverson MA. Review of Pre- and Post-1980 Buildings in CBECS - HVAC Equipment. Richland: Pacific Northwest National Lab. (PNNL); 2006. Dec [cited 2019 Aug 6]. Report No.: PNNL-20346. Available from: https://www.osti.gov/biblio/1013959.

    Google Scholar 

  77. Winiarski DW, Halverson MA, Butzbaugh JB, Cooke AL, Bandyopadhyay GK, Elliott DB. Analysis for Building Envelopes and Mechanical Systems Using 2012 CBECS Data. Richland: Pacific Northwest National Lab.(PNNL); 2018.

    Google Scholar 

  78. Polidori A, Fine PM, White V, Kwon PS. Pilot study of high-performance air filtration for classroom applications. Indoor Air. 2013;23(3):185–95.

    CAS  PubMed  Google Scholar 

  79. Chan WR, Parthasarathy S, Fisk WJ, McKone TE. Estimated effect of ventilation and filtration on chronic health risks in U.S. offices, schools, and retail stores. Indoor Air. 2016;26(2):331–43.

    CAS  PubMed  Google Scholar 

  80. ASHRAE S. Standard 62.1–2016 Ventilation for Acceptable Indoor Air Quality. Am Soc Heat. 2016.

  81. NAFA. Filtration for schools. Madison: National Air Filtration Association; 2012.

    Google Scholar 

  82. EPA. Indoor Air Quality Tools for Schools Coordinator’s Guide A Guide to Implementing an IAQ Program. U.S. Environmental Protection Agency; 2009.

    Google Scholar 

  83. Li H-W, Wu C-Y, Tepper F, Lee J-H, Lee CN. Removal and retention of viral aerosols by a novel alumina nanofiber filter. J Aerosol Sci. 2009 Jan 1;40(1):65–71.

    CAS  PubMed  Google Scholar 

  84. Gao X, Li Y, Leung GM. Ventilation control of indoor transmission of airborne diseases in an Urban Community. Indoor Built Environ. 2009;18(3):205–18.

    Google Scholar 

  85. Noakes CJ, Sleigh PA, Khan A. Appraising healthcare ventilation design from combined infection control and energy perspectives. HVACR Res. 2012;18(4):658–70.

    Google Scholar 

  86. Li Y, Leung GM, Tang JW, Yang X, Chao CY, Lin JZ, et al. Role of ventilation in airborne transmission of infectious agents in the built environment-a multidisciplinary systematic review. Indoor Air. 2007;17(1):2–18.

    CAS  PubMed  Google Scholar 

  87. Fisk WJ. The ventilation problem in schools: literature review. Indoor Air. 2017;27(6):1039–51.

    CAS  PubMed  Google Scholar 

  88. Foarde KK. Methodology to perform clean air delivery rate type determinations with microbiological aerosols. Aerosol Sci Technol. 1999 Feb 1;30(2):235–45.

    CAS  Google Scholar 

  89. Weiss MM, Weiss PD, Weiss DE, Weiss JB. Disrupting the Transmission of Influenza A: Face Masks and Ultraviolet Light as Control Measures. Am J Public Health. 2007;97(Supplement_1):S32–7.

    PubMed  PubMed Central  Google Scholar 

  90. Cowling BJ, Zhou Y, Ip DKM, Leung GM, Aiello AE. Face masks to prevent transmission of influenza virus: a systematic review. Epidemiol Infect. 2010 Apr;138(4):449–56.

    CAS  PubMed  Google Scholar 

  91. Garner JS, Committee HICPA. Guideline for isolation precautions in hospitals. Infect Control Hosp Epidemiol. 1996;17(1):54–80.

    Google Scholar 

  92. Getz WM, Carlson C, Dougherty E, Porco TC, Salter R. An agent-based model of school closing in under-vaccinated communities during measles outbreaks. SIMULATION. 2019;95(5):385–93.

    Google Scholar 

  93. Tuladhar E, Hazeleger WC, Koopmans M, Zwietering MH, Beumer RR, Duizer E. Residual viral and bacterial contamination of surfaces after cleaning and disinfection. Appl Environ Microbiol. 2012 Nov 1;78(21):7769–75.

    CAS  PubMed  PubMed Central  Google Scholar 

  94. Moorer WR. Antiviral activity of alcohol for surface disinfection. Int J Dent Hyg. 2003;1(3):138–42.

    CAS  PubMed  Google Scholar 

  95. Walker CM, Ko G. Effect of ultraviolet germicidal irradiation on viral aerosols. Environ Sci Technol. 2007 Aug 1;41(15):5460–5.

    CAS  PubMed  Google Scholar 

  96. Bahlke AM, Silverman HF, Ingraham HS. Effect of ultra-violet irradiation of classrooms on spread of mumps and Chickenpox in large rural central schools—a Progress report. Am J Public Health Nations Health. 1949;39(10):1321–30.

    CAS  PubMed  PubMed Central  Google Scholar 

  97. Nardell EA. Indoor environmental control of tuberculosis and other airborne infections. Indoor Air. 2016;26(1):79–87.

    CAS  PubMed  Google Scholar 

  98. Ko G, Burge HA, Nardell EA, Thompson KM. Estimation of tuberculosis risk and incidence under upper room ultraviolet germicidal irradiation in a waiting room in a hypothetical scenario. Risk Anal. 2001;21(4):657–74.

    CAS  PubMed  Google Scholar 

  99. Loharikar A, Black S, Ramirez E, Morita J. Measles Post-exposure Prophylaxis and Precautions. Chicago: Department of Public Health; 2015. (Measles Post-exposure Guidelines). Available from: https://www.chicagohan.org/documents/14171/107461/CDPH+HAN%3A+Measles+Post-exposure+Prophylaxis+and+Precautions.pdf/8dabcd10-f54c-4e2e-8a09-987961a20aca.

    Google Scholar 

  100. Kujundzic E, Hernandez M, Miller SL. Ultraviolet germicidal irradiation inactivation of airborne fungal spores and bacteria in upper-room air and HVAC in-duct configurations. J Environ Eng Sci. 2007;6(1):1–9.

    CAS  Google Scholar 

  101. Nardell EA, Bucher SJ, Brickner PW, Wang C, Vincent RL, Becan-McBride K, et al. Safety of upper-room ultraviolet germicidal air disinfection for room occupants: results from the tuberculosis ultraviolet shelter study. Public Health Rep. 2008 Jan 1;123(1):52–60.

    PubMed  PubMed Central  Google Scholar 

  102. Wyll SA, Witte JJ. Measles in previously vaccinated children: an epidemiological study. JAMA. 1971;216(8):1306–10.

    CAS  PubMed  Google Scholar 

  103. Gustafson TL, Lievens AW, Brunell PA, Moellenberg RG, Buttery CMG, Sehulster LM. Measles outbreak in a fully immunized secondary-school population. N Engl J Med. 1987 Mar 26;316(13):771–4.

    CAS  PubMed  Google Scholar 

  104. Morse D, O’shea M, Hamilton G, Soltanpoor N, Leece G, Miller E, et al. Outbreak of measles in a teenage school population: the need to immunize susceptible adolescents. Epidemiol Infect. 1994;113(2):355–65.

    CAS  PubMed  PubMed Central  Google Scholar 

  105. Centers for Disease Control and Prevention (CDC). Measles outbreak among school-aged children--Juneau, Alaska, 1996. MMWR Morb Mortal Wkly Rep. 1996;45(36):777–80.

    Google Scholar 

  106. Mossong J, Muller CP. Estimation of the basic reproduction number of measles during an outbreak in a partially vaccinated population. Epidemiol Infect. 2000 Apr;124(2):273–8.

    CAS  PubMed  PubMed Central  Google Scholar 

  107. Wichmann O, Hellenbrand W, Sagebiel D, Santibanez S, Ahlemeyer G, Vogt G, et al. Large measles outbreak at a German public school, 2006. Pediatr Infect Dis J. 2007 Sep;26(9):782.

    PubMed  Google Scholar 

  108. Centers for Disease Control and Prevention (CDC). Outbreak of measles--San Diego, California, January–February 2008. MMWR Morb Mortal Wkly Rep. 2008;57(8):203–6.

    Google Scholar 

  109. Ma R, Lu L, Zhangzhu J, Chen M, Yu X, Wang F, et al. A measles outbreak in a middle school with high vaccination coverage and evidence of prior immunity among cases, Beijing, P.R. China. Vaccine. 2016 Apr 4;34(15):1853–60.

    PubMed  Google Scholar 

  110. Hutchins SS, Markowitz LE, Mead P, Mixon D, Sheline J, Greenberg N, et al. A school-based measles outbreak: the effect of a selective revaccination policy and risk factors for vaccine failure. Am J Epidemiol. 1990 Jul;132(1):157–68.

    CAS  PubMed  Google Scholar 

  111. Matson DO, Byington C, Canfield M, Albrecht P, Feigin RD. Investigation of a measles outbreak in a fully vaccinated school population including serum studies before and after revaccination. Pediatr Infect Dis J. 1993 Apr;12(4):292–9.

    CAS  PubMed  Google Scholar 

  112. Paunio M, Peltola H, Valle M, Davidkin I, Virtanen M, Heinonen OP. Explosive school-based measles OutbreakIntense exposure may have resulted in high risk, even among Revaccinees. Am J Epidemiol. 1998;148(11):1103–10.

    CAS  PubMed  Google Scholar 

  113. Lyons RA, Jones HI, Salmon RL. Successful control of a school based measles outbreak by immunization. Epidemiol Infect. 1994 Oct;113(2):367–75.

    CAS  PubMed  PubMed Central  Google Scholar 

  114. Lynn TV, Beller M, Funk EA, Middaugh JP, Ritter D, Rota RA, et al. Incremental Effectiveness of 2 Doses of Measles-Containing Vaccine Compared with 1 Dose among High School Students during an Outbreak. J Infect Dis. 2004;189(Supplement_1):S86–90.

    PubMed  Google Scholar 

  115. Yeung LF, Lurie P, Dayan G, Eduardo E, Britz PH, Redd SB, et al. A limited measles outbreak in a Highly vaccinated US boarding school. Pediatrics. 2005 Dec 1;116(6):1287–91.

    PubMed  Google Scholar 

  116. De Serres G, Boulianne N, Defay F, Brousseau N, Benoît M, Lacoursière S, et al. Higher risk of measles when the first dose of a 2-dose schedule of measles vaccine is given at 12–14 months versus 15 months of age. Clin Infect Dis. 2012 Aug 1;55(3):394–402.

    PubMed  Google Scholar 

Download references

Acknowledgements

We are thankful to our colleague Marianne Lahaie Luna from University of Toronto, Dalla Lana School of Public Health, Occupational & Environmental Health Division for her assistance with reviewing the main manuscript and appendices and in that line improved the study significantly.

Funding

Not applicable.

Author information

Authors and Affiliations

Authors

Contributions

PA, JGCL, JGA designed the research. PA and ZK performed the research and wrote the paper. All authors reviewed and approved the final manuscript.

Corresponding author

Correspondence to Parham Azimi.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

Not applicable.

Additional information

Publisher’s Note

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

The original version of this article was revised: the authors identified an error in the author name of dr. Jose Guillermo Cedeno Laurent.

The incorrect author name is: Jose Guillermo Cedeno Cedeno Laurent

The correct author name is: Jose Guillermo Cedeno Laurent

The author group has been updated above and the original article [1] has been corrected.

Supplementary information

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Azimi, P., Keshavarz, Z., Cedeno Laurent, J.G. et al. Estimating the nationwide transmission risk of measles in US schools and impacts of vaccination and supplemental infection control strategies. BMC Infect Dis 20, 497 (2020). https://doi.org/10.1186/s12879-020-05200-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12879-020-05200-6

Keywords