Abstract
Background This study examined the correlation of classroom ventilation (air exchanges per hour (ACH)) and exposure to CO2 ≥1,000 ppm with the incidence of SARS-CoV-2 over a 20-month period in a specialized school for students with intellectual and developmental disabilities (IDD). These students were at a higher risk of respiratory infection from SARS-CoV-2 due to challenges in tolerating mitigation measures (e.g. masking). One in-school measure proposed to help mitigate the risk of SARS-CoV-2 infection in schools is increased ventilation.
Methods We established a community-engaged research partnership between the University of Rochester and the Mary Cariola Center school for students with IDD. Ambient CO2 levels were measured in 100 school rooms, and air changes per hour (ACH) were calculated. The number of SARS-CoV-2 cases for each room was collected over 20 months.
Results 97% of rooms had an estimated ACH ≤4.0, with 7% having CO2 levels ≥2,000 ppm for up to 3 hours per school day. A statistically significant correlation was found between the time that a room had CO2 levels ≥1,000 ppm and SARS-CoV-2 PCR tests normalized to room occupancy, accounting for 43% of the variance. No statistically significant correlation was found for room ACH and per-room SARS-CoV-2 cases. Rooms with ventilation systems using MERV-13 filters had lower SARS-CoV-2-positive PCR counts. These findings led to ongoing efforts to upgrade the ventilation systems in this community-engaged research project.
Conclusions There was a statistically significant correlation between the total time of room CO2 concentrations ≥1,000 and SARS-CoV-2 cases in an IDD school. Merv-13 filters appear to decrease the incidence of SARS-CoV-2 infection. This research partnership identified areas for improving in-school ventilation.
Introduction
A variety of mitigation measures have been employed to prevent SARS-CoV-2 exposure and infection, including immunologic (e.g. vaccination, monoclonal antibodies), antiviral (e.g. paxlovid), and isolation methods (e.g. masking, physical distancing, stay-at-home policies). A key approach to environmental mitigation has been the use of high volume ventilation of enclosed spaces (i.e. adequate flow of fresh, uncontaminated air), which has been shown to decrease the risk of infection from a variety of respiratory pathogens [1–3], and this effect has been widely presumed to include SARS-CoV-2. Schools were an immediate focus in many areas with respect to mitigation measures for SARS-CoV-2, including potentially costly recommendations for improving ventilation. However, very little data are available assessing the actual impact of ventilation on in-school transmission of SARS-CoV-2.
The SARS-CoV-2 β-coronavirus is spread primarily via airborne droplets and aerosols, which can remain circulating for up to three hours [4]. Ventilation is defined as the process of supplying clean air to the indoor air of a dwelling by natural (e.g., open windows if the outdoor air is sufficiently clean) or mechanical means (e.g. heating, ventilation, and air conditioning; HVAC) often involving recirculating the air after cleaning with a filter [5].
Increasing ventilation is one means of reducing exposure to airborne viral pathogens, and increasing ventilation rates and upgrading inline HVAC particle filters (e.g. from MERV-8 to MERV-13) were widely recommended for public venues including businesses, transportation, and schools [6–8]. For schools specifically, improving ventilation as a risk reduction strategy was the subject of many recommendations [8, 8–11]. However, very little robust data are available assessing the actual impact of improved ventilation on in-school transmission of SARS-CoV-2. One study found 74% lower aggregate rates of COVID-19 in schools with mechanical versus natural ventilation [12]. While there have been reports of outbreaks within schools [13, 14], a number of epidemiological studies have concluded that students and school staff had relatively low risk of SARS-CoV-2 in-school transmission [15–17] and schools were unlikely to be a major factor in the pandemic spread [18].
Despite the larger number of epidemiological studies in public schools, we are not aware of similar studies in schools for students with intellectual and developmental disabilities (IDD), who face greater risks [19–24], and impacts from, infection. The majority of students with IDD are served by IDD-specialized schools that provide educational services, comprehensive clinical services, speech and psychological interventions, nutritional needs, and physical care. Staff-to-student ratios are much higher than in typical K-12 schools, with more frequent, sustained and close physical contact between staff and students [25]. Students with an IDD are at higher risk of exposure to SARS-CoV-2 [25], severe infection, and mortality [24]. The IDD spectrum covers many conditions (e.g. autism spectrum disorder, genetic disorders), which are associated with immune abnormalities that increase risk of SARS-CoV-2 and poor vaccine responses [19–24]. Many students with an IDD have difficulty with conventional infection prevention measures such as masking and social distancing. This also puts staff in IDD settings at higher risk of exposure to SARS-CoV-2 from asymptomatic but infected students. Conversely, the impact on IDD students, and their caregivers, of school closures, isolation and quarantine is heightened [26]. IDD students receive many specialized services (e.g. physical and speech therapy) in school that are not readily available at home. Thus, it is essential to identify interventions that are implementable in the unique environment of IDD-specialized schools that retain students in school while mitigating the risk of SARS-CoV-2 infection. Increased ventilation has been proposed as one such measure.
Time series measurements of CO2 levels have been used as a frequent surrogate for ventilation adequacy in schools [8, 27–32]. As CO2 is exhaled during respiration, it accumulates in indoor spaces depending on the air exchange rates. The dynamic balance between CO2 production and ventilation with air having a lower CO2 concentration determines the time-varying CO2 concentration within a room. Accordingly, time series measurements of room CO2 concentrations can be used to derive a surrogate marker for ventilation adequacy, typically measured as air changes per hour (ACH) [27]. International guidelines recommend a minimum of 4 ACH for school classrooms, and 8-15 ACH for auditoriums, athletic facilities, cafeterias and other high-density or high-activity enclosed spaces [5]. Mismatched ventilation can result in elevated indoor CO2 levels that can affect learning and cognitive performance [33], areas of key focus in schools for students with an IDD. Unfortunately, many school buildings have outdated HVAC systems, resulting in low or very low ACH [6, 34, 35].
The National Institutes of Health (NIH) program “Rapid Acceleration of Diagnostics for Underserved Populations (RADx-UP)” was designed to study ways to support vulnerable populations during the COVID-19 pandemic. A major goal of the program was to study strategies to keep vulnerable populations of students in school while mitigating the risk of exposure and in-school spread of SARS-CoV-2 among students and staff. Primarily designed to study the utility of screening for infection with frequent molecular testing, our RADx-UP project work has also focused on the effect of a variety of measures on mitigating exposure and infection in an IDD-specialized school setting. Here we report the results of an extensive study of ventilation across 3 buildings and 100 rooms in an IDD school for children and young adults (ages 3-21 years) during the SARS-CoV-2 pandemic.
Methods
Human subjects protection
This study was approved by the University of Rochester Medical Center Research Subjects Review Board (RSRB STUDY00005294). Participants provided written informed consent, and parents or guardians provided consent for minors. Recruitment occurred April 1, 2021 and through May 31 2023. All subjects’ information and research data were coded in compliance with the Department of Health and Human Services Regulations for the Protection of Human Subjects (45 CFR 46.101(b) (4)).
Room and building data
The school in this study is composed of 3 different buildings, each with a varying number of classrooms, therapy and activity rooms, gyms, and a variety of offices. We purposefully selected 100 school rooms, from a total of 267, reflecting the different room types in use across the 3 different buildings. Metadata collected included a unique room identifier, building, square footage, ceiling height, total occupancy range, and room type (e.g. classroom, therapy, etc.). In addition, we obtained further data regarding the heating, ventilation, and air conditioning (HVAC) systems used in each building relevant to the study, including type of ventilation (1 versus 2 stage), percentage of outside air mixing, and filtering (Supplemental File S1 File).
SARS-CoV-2 case data
Total counts of positive SARS-CoV-2 PCR tests were tabulated for the period from August 2021 - August 2022 for the rooms where CO2 concentrations were measured. During this period at the school, most individuals spent the majority of the day in a single classroom, and thus positive tests were associated with the primary classroom assignment for any individual. Individuals who tested positive within 30 days of a prior positive PCR were considered single infections.
CO2 measurements
In November 2022, we measured ambient CO2 levels, temperature, and humidity in 100 school rooms using ARANET4 sensors (Aranet, Aurora Colorado) with firmware version 1.2, with one sensor per classroom, twenty classrooms at a time. Monitors were placed on the ceiling near the air return vents in the center of the room to ensure mixing and avoid dead spaces within each room. In a small number of rooms, where ceiling placement was not possible, monitors were placed on side walls near air return vents. The indicator screens of the devices were covered so that classroom occupants would not alter ventilation (e.g. open windows or doors) outside of their usual pattern. Sensors were set to record ambient CO2, temperature, and humidity levels every 2 minutes continuously for 24-72 hours. Data were downloaded to an iPad (Apple Inc., Cupertino CA) using the Aranet manufacturer’s software as a .CSV file. Each record in a .CSV file identified the sensor and time of each measurement. A separate .CSV file identified the school room and sampling start and stop times. Individual .CSV files were then processed with a custom R program that added the unique room ID number to each record based on the sampling start and stop times and the unique sensor ID. Duplicate records were removed. The final data set contained 150,286 individual measurements (Supplemental File S2 File). To assess intra-day variation, some rooms had data continuously recorded for up to 120 hours, both during the normal weekly school day, and over the weekend to assess temperature variations.
Air changes per hour estimation
To estimate air changes per hour (ACH), CO2 time series for each room (Fig 1) were first subjected to Gaussian blurring to scale σ = 10 and then local maxima identified using Mathematica (vers 13.1; FindPeaks; Wolfram Research). Local minima were similarly identified by analyzing the negative time series, resulting in pairs of , which were then filtered to exclude pairs with small declines in CO2 levels where
(2 rooms) [27, 36]. This resulted in 98 rooms with measurements adequate for estimating ACH. We then fitted an exponential decay curve [27] of the form:
to each time series segment between successive
pairs, where i = number of segments identified for an individual room, α is air exchanges per hour, and t is the time from
to
. Fitting was performed using the NonlinearModelFit function in Mathematica with the Levenberg-Marquardt algorithm option (Fig 1). The room airflow needed to achieve 4 air changes per hour was estimated using the volume of the room (cubic feet) divided by 4.
Example of CO2 time series measurement, peak identification and fitting. A: CO2 levels versus time in curve with Gaussian smoothing window of 10, and peak (red triangle) and valley (blue triangle) identification. Segments with a long tail were truncated at the point where subsequent readings were only 2% lower than the mean of the previous 10 readings. B. Eq 1 was fitted to each set of interval data using NonlinearModelFit function in Mathematica, and the Levenberg-Marquardt algorithm option. In those rooms with multiple (peak, valley) segments, the final room ACH value was estimated by averaging the estimates from all segments.
Statistical Analysis
The incidence of positive PCR tests in each room over 12 months was normalized to mean room occupancy from 1 month worth of attendance and occupancy data during a fall semester of the school year in 2022 during the COVID-19 pandemic. Correlation between occupancy normalized counts of positive SARS-CoV-2 tests and air changes per hour (ACH) or exposure time to CO2 ≥ 1,000 ppm were calculated using the non-parametric Spearman’s rank correlation test, with the explained variance reported as the Spearman’s ρ statistic.
Data Availability
Room data is available in Supplemental File S1 File, and the complete CO2 measurements used in the analysis can be found in Supplemental File S2 File.
Results
Room characteristics
We sampled 100 rooms spread across three buildings. Two buildings had single stage HVAC systems, which supplied ventilation only when triggered by the need for heating during the sampling period (November), while the third building had a two stage HVAC system. Set-points for heating were set at 70°F from 6 AM - 6 PM, and 65°F outside those hours. Outside temperatures during the measurement periods were all ≤65°F. For all HVAC systems, air provided to rooms was 20% outside and 80% recirculated air. One building had HVAC systems with inline MERV-13 filters installed to mitigate SARS-CoV-2 transmission from recirculating air. The HVAC systems in two other buildings had MERV-11 filters as the blowers were not able to overcome the higher filter resistance from MERV-13 filters. The distribution of maximum room occupancy is shown in Figure 2A. Room dimensions were measured, with mean room area 752 ± 492 ft2, and mean ceiling heights 8.6 ± 2.0 ft. As expected, peak room occupancy varied linearly with room volume (Fig. 2B). Windows remained closed in all classrooms during the measurements.
A. Distribution of maximum room occupancy in persons for the 100 rooms in the study. B. Room occupancy versus volume in ft3. Most rooms had ceiling heights of 9-10 feet, and the relationship between room volume and area is close to linear. Larger rooms had higher occupancy levels.
CO2 level profiles
Time series measurements of room CO2 for the 98 rooms analyzed showed an expected increase in CO2 concentrations during the day, with a decline overnight to baseline. Several patterns were evident (Fig 3A), including a simple peak-trough, dual peak-troughs, and multiple peak-trough pairs. We then examined if peak CO2 levels in each room could be correlated with room volume (Fig 3B), air changes per hour calculated from each peak-trough pair (Fig 3C), or maximum room occupancy (Fig 3D). No obvious correlations were present, likely due to the variation in ventilation systems and room types among the 3 buildings.
A. Selected patterns of CO2 concentration over the school day. Color coding indicates CO2 concentration stratification by concentration bands commonly used by the American Society of Heating, Refrigerating and Cooling Engineers (ASHRAE) (green, < 1,000 ppm; yellow 1,000 - 2,000 ppm; red, geq 2,000 ppm). Minimum (blue) and maximum (red) CO2 concentration versus B. Room volume, C. Estimated ACH, and D. Maximum room occupancy.
Room CO2 concentrations
CO2 levels ≥2500 ppm have been shown to decrease cognitive task performance in neurotypical adults on standard test instruments [37, 38]. We were therefore interested in the amount of time each room’s CO2 concentrations were within four different ranges of CO2 during an 8 hour school day (7AM - 3PM). Several patterns were evident (Fig 4. Rooms in Group 1 (n=18) remained within the recommended values of CO2 (≤1000 ppm) [5] throughout the entire period. In contrast rooms in Group 2 (n=74) had CO2 levels of 1001-2000 ppm, which have been shown to cause drowsiness and modest decreases in performance of cognitive tasks in neurotypical adults [37, 38] and children [39–42]. Of note, occupants of rooms in Group 3 (n=6; 2001-3000 ppm) and Group 4 (n=2; ≥3000 ppm) had approximately 1-6 hours at moderately high CO2 levels, although still below the OSHA occupational limits of ≤5000 ppm. We confirmed the elevated CO2 levels in Group 3 and 4 rooms with 3 more full day time series measurements.
The distribution of CO2 exposure during an 8 hour school day. Stacked plots indicate the amount of time spent at each CO2 level during the school day. Horizontal grid lines indicate quartiles. Each bar represents a room, with labels coded by building. Rooms are labeled by function: classroom, therapy, treatment, activity, gym, office, kitchen, break, copy, music, and technology.
We also investigated whether the maximum CO2 concentrations observed in each room correlated with functional room type and building. Several different functional room types were present among the 3 buildings: classroom, therapy, treatment, activity, gym, office, kitchen, break, copy, music, and technology. Classrooms generally are occupied during the day by 7-8 students, 1 teacher, and 2-4 aids, with activity levels that produce moderate amounts of CO2. In contrast, music, activity, and gym rooms all have activities that may generate larger emissions of CO2. Each building also has a different type and vintage of HVAC system, with potential variation in ventilation. To assess the distribution of peak CO2 values, we sorted rooms by type and building, as shown in Fig 5. Note that building B had the widest distribution of maximum measured CO2 values, and these were found in classrooms. While a number of the measured CO2 concentrations were above the ASHRAE guidelines for classrooms for varying periods of time [5], none exceeded OSHA’s exposure limits [43].
Each plot displays a density distribution, along with individual measurements, for the room types indicated. Horizontal grid lines indicate ranges of CO2 which are acceptable (green), moderately elevated (yellow), and high (red). Distributions for room types are further categorized by building.
Temperature versus CO2 relationship
The heating HVAC systems for each of the 3 buildings had zone heating and cooling. Two buildings were described as single stage, with mechanical room ventilation occurring only when the room needed to be heated or cooled, while one was described as dual stage, with both a baseline level of ventilation and then increased airflow when heating or cooling was required. CO2 values are a function of both ongoing production and ventilation. We hypothesized that room CO2 measurements would decrease during periods when the HVAC system was activated during room occupancy. Fig 6 shows superimposed plots of CO2 levels and temperature readings over time in a representative sample of 9 out of 100 rooms. Note the oscillations in temperature as the heating system engages (e.g. Figure 6A, E, F, I) in some rooms. Of note, Fig 6 panels C,E,F and I all show instances where the temperatures increased rapidly, indicating active mechanical heating and ventilation, with a corresponding decline in room CO2 levels. These findings suggest that, at ambient temperatures above (in winter) or below (in summer) the thermostat set point, HVAC cycling coupled to temperature regulation may have decreased CO2 clearance, since the blowers circulate air at higher flow rates only when temperatures fell outside of the thermostat set point.
Plots of CO2 levels (left axis) and temperature (right axis) showing representative room patterns. (A-C) show rooms that had CO2 levels ≥2000 ppm, (D-F) had CO2 peak levels 1000-2000 ppm, and (G-I) had CO2 ≤1000 from 7AM-7PM. CO2 are coded by ≤ 1,000 ppm (green), 1,000-2,000 ppm (yellow), and > 2,000 ppm (red). Temperature is represented by the blue line.
Air changes per hour estimates
Using the CO2 level profiles for each room, we derived air changes per hour (ACH) by fitting an exponential decay model to the time series segments where ΔCO2 ≥ 120ppm between peak and trough in a decay segment (Figure S1 Fig). The distribution of ACH values is shown in Fig 7A, stratified by very low (red; ACH≤1.0), low (yellow; ACH 1.0-2.0), and moderately acceptable ACH≥2.0). All were less than the recommended ACH level of ≥4.0 [5]. Fig 7B demonstrates clearance of a 1200 ppm CO2 load at the estimated ACH for each room. Clearance times of 4-9 hours were consistent with measured, overnight, CO2 level declines. Finally, we compared the measured ACH values for each room with the minimum ACH recommended based on maximal room occupancy and measured room area and volume (Fig 7C. The estimate of recommended ACH based on ASHRAE 62.1 Table 6-1, Classrooms (ages 5 to 8) require an air flow 10 ft3 min−1 per person and 0.12 ft3 min−1 per ft2 of area. The air flow was divided by the room volume to obtain recommended ACH. This calculation may underestimate required air flow, and thus recommended ACH, as it does not account for re-mixing of outside (i.e. fresh) air (20%) and recirculated air (80%) used by the HVAC systems.
(A) Distribution of estimated ACH. Only four rooms had an estimated ACH ≥2.0. The recommended ACH for schools is ≥ 4.0 exchanges of fresh air per hour [5, 9]. (B) Each ACH measure was used to model the decrease in CO2 levels over time from a starting concentration of 1200 ppm CO2 (gray lines). Colored reference lines show expected CO2 clearance with 0.1, 0.5, 1.0, 2.0, and 4.0 ACH. (C) Estimated room ACH plotted against recommended ACH for each room based maximum room occupancy as recommended for school rooms American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) standard 61.1 calculated from room area and occupancy [5]. Dashed line is where ACHmeasured = ACHrecommended.
Incidence of SARS-CoV-2-positive PCR tests
We next sought to determine if rooms with lower ACH had higher numbers of SARS-CoV-2 infections (Fig. 8), and whether MERV-13 filters altered this result. During the COVID-19 pandemic, the use of high efficiency inline particle filters (e.g. MERV-13) were recommended to improve filtration of SARS-CoV-2 within HVAC systems [2, 8–10]. Higher filter efficiency, however, often results in lower ventilation rates and greater strain on the HVAC mechanical blowers such that older HVAC systems may not accommodate MERV-13 filters. In this school, two building HVAC systems could not provide the higher mechanical load required to circulate air through MERV-13 filters, and continued to use MERV-11 filters. We first gathered data for a 1-year period for positive SARS-CoV-2 PCR test counts per classroom for those rooms that had measurements of CO2 levels, estimates of ACH, and MERV filter data. Positive PCR tests were normalized to mean classroom occupancy. A comparison of rooms with MERV-11 filters demonstrates that they have a higher occupancy normalized positive PCR tests per room compared to rooms with MERV-13 filters (Fig. 8A). Next, total hours of CO2 levels ≥1,000 ppm were tabulated for each room. A significant correlation was found between hours of classroom CO2 levels ≥ 1,000 ppm and the incidence of SARS-CoV-2 infection (Spearman rank correlation; p=0.01937, ρ=0.2864). Lastly, no statistically significant correlation was found between estimated ACH and the incidence of positive PCR tests normalized to mean classroom occupancy.
(A) Difference in PCR tests per room, normalized by room occupancy, between rooms with HVAC systems using MERV-11 (red) versus MERV-13 filters. Statistical comparison with the Mann-Whitney U-test (p<0.0012). (B) The time in each room where the ambient CO2 ≥ 1,000 ppm was plotted against positive PCR tests in that room normalized by mean room occupancy. Rooms with CO2 ≤ 1,000 ppm for the entire day were excluded (lighter markers). MERV filter status for the building rooms is shown (red, MERV-11; blue, MERV-13). Linear regression (dashed line) had an R2 = 0.46, with an R2 = 0.4537 explaining 45% of the variance. (C) ACH versus positive SARS-CoV-2 PCR tests per room normalized to room occupancy. No statistically significant correlation was found (R2 = 0.0036).
Discussion
Our overall findings revealed only a small number of PCR-detected SARS-CoV-2 cases in this IDD school over the course of one year. Our finding that less than 20% of staff and students had PCR-confirmed SARS-CoV-2 infections is in agreement with other school studies. Investigators in Norway [17], Germany [44], Italy [13], United States [45–47] and Ireland [16] found minimal or no evidence for in-school SARS-CoV-2 transmission. Although there have been reports of small secondary infection clusters [14, 48], larger studies have concluded that, overall, schools were unlikely to be a major factor in the pandemic spread [18]. Although one retrospective cohort study in Italy found that schools with mechanical ventilation had a 78% lower incidence of SARS-CoV-2 cases, their findings may have limited applicability in geographic regions where seasonal heating and cooling is necessary. Additionally, the conditions in IDD schools are likely somewhat different than public schools reported in the above publications. In the school studied in this work, only 30% of the students were able to tolerate masking, while all of the professional staff was masked. Thus, further study of the combined effects of mitigation measures is needed.
Using CO2 levels to estimate ACH, we found that the school ventilation systems did not achieve ACH ≥ 4, which is consistent with other reports [49–52]. However, the use of CO2 measurements to derive ACH comes with several important caveats [53]. To fit exponential decay models to CO2 concentrations, a sufficient decay time is needed after CO2 loading [27, 54]. In this study, the CO2 concentration decay phases occurred at the end of the school day and overnight. However, the HVAC systems had slightly different temperature set-points at night, which may have led to different ventilation rates, underestimating daytime ACH [27]. In addition, the HVAC systems were set primarily to maintain temperature, and thus higher CO2 concentrations may have been due to lower air circulation when the heat generated by room occupancy decreased the need for higher rates of heated air circulation [27]. Ventilation rates can also vary in both naturally- and mechanically-ventilated buildings, depending the indoor-outdoor temperature difference, wind speed, heating and cooling load, building envelope leakage, and other factors that may change over time [27, 54]. Our measurements were performed in the winter season (November) over 3 contiguous days, when heating was necessary due to outside temperatures, which is relevant to both HVAC use and the respiratory virus season in the northern hemisphere. These factors suggest that use of CO2 level ranges normalized to occupancy, with the above caveats, may be a practical method to estimate ventilation adequacy in schools than CO2 derived ACH measures.
Many have suggested that CO2 monitoring could be used as a surrogate for adequate ventilation [27, 54–58]. While the primary goal of this work was to provide an estimate of the correlation between SARS-CoV-2 infections and room ventilation in an IDD school, elevated CO2 concentrations on their own may decrease cognitive task performance [37, 39–42, 59], concentration [60], increase cough and rhinitis [61], and teacher’s reporting of health symptoms (e.g. fatigue, headache) [32]. Especially in schools for IDD students, improving ventilation, and achieving recommended CO2 levels is likely to be important for school attendance [62], learning, and cognitive task performance [63, 64]. Monitoring CO2 levels would provide real-time data in classrooms, and could be used to motivate interventions to improve ventilation and decrease CO2 levels, from simple measures such as opening windows when possible to large-scale upgrades of school HVAC systems. Such measures should, however, be correlated with more rigorous measures of ventilation to account for numerous factors important for interpretation of the results [53]. Finally, it is important to note that we did not measure an individual’s exposure over time to CO2 concentrations, but rather concentrations in rooms. Individuals generally move between classrooms, therapy rooms, the outdoors, and other school areas, although less so during COVID.
This study also highlights the potential trade-offs in terms of improving in-line filtration and ventilation of particles and droplets to prevent infection and the potential expense required to upgrade HVAC systems in many schools [49]. Specifically, the rooms supplied by HVAC systems using MERV-11 inline filters had a higher incidence of SARS-CoV-2 than those with MERV-13 filters. Circulating air through a higher resistance particulate filter may not be possible for many older HVAC systems, requiring an expensive upgrade. In addition, filters may decrease ventilation rates, raising ambient CO2, particulate, and volatile organic compound concentrations. Increasing the percentage of outdoor air brought in creates a trade-off between improved ventilation and the need to heat, cool, and condition the humidity of the incoming fresh air, increasing energy costs [37, 65]. Similarly, for IDD schools, measures such as opening windows create another trade off between improved ventilation and student safety. Given the vintage of most public school buildings and HVAC systems, as well as the budgets available for upgrades and remediation measures, this issue may benefit from more public funding and policy discussions..
There are a few important limitations of this work. As noted above, CO2 concentration level decay may vary with the time of day and the room occupancy, and thus tracer gas methods for estimating ACH may have yielded more accurate results. These methods, however, also have limitations, including difficulty in performing them in actively occupied classrooms. In addition, CO2 levels may vary seasonally with the ability to open windows, changes in classroom occupancy patterns, and air blower activation patterns at summer versus winter temperatures. Our CO2 measurements were performed during the fall season over 48-72 hours, while the incidence per room occupancy of SARS-CoV-2 was calculated over an entire year.
One goal of community engaged research partnerships is to have direct impact on improving outcomes highly valued by the community. The comprehensive sampling in 100 rooms across 3 buildings revealed a wide variation in ventilation rates, as well as number of hours exposed to CO2 levels greater than 1,000 ppm. Only 8 rooms had levels of CO2 ≥ 2,000 ppm for more than 30-130 minutes over the course of a school day, and of those 7 of 8 rooms were in the same building. These findings were subsequently validated by an independent HVAC contractor engaged by the school. Of note, the rooms with ventilation served by HVAC systems with MERV-13 filters had lower rates of PCR test positives for SARS-CoV-2. Overall, these results have sparked discussion about upgrade of the ventilation systems. The school continues to seek state and federal funding in an on-going effort to ensure and exceed the ventilation standards recommended for schools.
Conclusions
There was a statistically significant correlation between rooms with CO2 levels ≥ 1,000 ppm and SARS-CoV-2 cases in an IDD school. Rooms served by HVAC systems with in-line MERV-13 filters had a lower incidence of SARS-CoV-2-positive PCR tests. This research partnership identified areas for improving in-school ventilation and the use of in-line, high efficiency air filtration.
Supporting information
S1 Fig. ACH Estimation. ACH curve fitting results for each room in the study.
S1 File. Room Data. A CSV file containing data with room characteristics, including building, square footage and volume, occupancy, CO2 concentration times, SARS-CoV-2 PCR positive subjects per room.
S2 File. CO2 Time Series Data. A CSV file containing the time series CO2 concentration data from each room used to derive the ACH estimates.
Authorship Contributions
This study was conceived and designed by MZ, KZ, AP, SS, AR, JH, AD, SD, EF, CLS, AC, AMC, and JF. Data were collected by SS, AP, and AR. Data wrangling was performed by KK, AMC, MZ, LA, and JW. Analysis was performed by MZ, CLS, and JF. The manuscript was drafted by MZ and edited by all co-authors.
Acknowledgments
The authors would like to thank the Mary Cariola Center school community for their enthusiastic support and extensive partnership in this work. The authors would also like to thank the technologists of the UR Medicine Central Laboratory for their rapid performance of SARS-CoV-2 PCR tests. We are very grateful for Dr. Philip K. Hopke’s (University of Rochester) expert guidance on the American Society of Heating, Refrigerating and Air-Conditioning Engineers (ASHRAE) standards for schools, CO2 measurements and ACH calculations, and his critical reading of the manuscript. The inspiration for CO2 measurements came from the outstanding social media videos of Dr. Katrine Wallace (University of Illinois; @epidemiologistkat). Research reported in this Rapid Acceleration of Diagnostics – Underserved Populations (RADx-UP) publication was supported by the National Institutes of Health under OT2 HD107553. This work was also supported by the University of Rochester Intellectual and Developmental Disabilities Research Center (P50 HD103536) from the NIH Eunice Kennedy Shriver Institute of Child Health and Human Development, and the University of Rochester Clinical and Translational Science Award (UL1 TR002001) from the National Center for Advancing Translational Sciences of the National Institutes of Health. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
Footnotes
This version of the manuscript has been revised to update the format with the figures displayed within the text, and to add caveats to the study based on reader comments.