Bayesian Calibration of Using CO2 Sensors to Assess Ventilation Conditions and Associated COVID-19 Airborne Aerosol Transmission Risk in Schools ================================================================================================================================================ * Danlin Hou * Ali Katal * Liangzhu (Leon) Wang ## Abstract Ventilation rate plays a significant role in preventing the airborne transmission of diseases in indoor spaces. Classrooms are a considerable challenge during the COVID-19 pandemic because of large occupancy density and mainly poor ventilation conditions. The indoor CO2 level may be used as an index for estimating the ventilation rate and airborne infection risk. In this work, we analyzed a one-day measurement of CO2 levels in three schools to estimate the ventilation rate and airborne infection risk. Sensitivity analysis and Bayesian calibration methods were applied to identify uncertainties and calibrate key parameters. The outdoor ventilation rate with a 95% confidence was 1.96 ± 0.31ACH for Room 1 with mechanical ventilation and fully open window, 0.40 ± 0.08 ACH for Rooms 2, and 0.79 ± 0.06 ACH for Room 3 with only windows open. A time-averaged CO2 level < 450 ppm is equivalent to a ventilation rate > 10 ACH in all three rooms. We also defined the probability of the COVID-19 airborne infection risk associated with ventilation uncertainties. The outdoor ventilation threshold to prevent classroom COVID-19 aerosol spreading is between 3 – 8 ACH, and the CO2 threshold is around 500 ppm of a school day (< 8 hr) for the three schools. **Practical Implications** The actual outdoor ventilation rate in a room cannot be easily measured, but it can be calculated by measuring the transient indoor CO2 level. Uncertainty in input parameters can result in uncertainty in the calculated ventilation rate. Our three classrooms study shows that the estimated ventilation rate considering various input parameters’ uncertainties is between ± 8-20 %. As a result, the uncertainty of the ventilation rate contributes to the estimated COVID-19 airborne aerosol infection risk’s uncertainty up to ± 10 %. Other studies can apply the proposed Bayesian and MCMC method to estimating building ventilation rates and airborne aerosol infection risks based on actual measurement data such as CO2 levels with uncertainties and sensitivity of input parameters identified. The outdoor ventilation rate and CO2 threshold values as functions of exposure times could be used as the baseline models to develop correlations to be implemented by cheap/portable sensors to be applied in similar situations to monitor ventilation conditions and airborne risk levels. Keywords * Airborne infection * COVID-19 * Bayesian calibration * Markov Chain Monte Carlo * Ventilation rate * CO2 * Risk ## 1. Introduction Airborne transmission of relatively small aerosol droplets plays a dominant role in spreading SARS-CoV-2 (hereafter as COVID-19), especially in indoor spaces 1,2. School classrooms pose a considerable challenge because of the large occupants’ density, mandatory presence of students, and ventilation conditions of concern 3. For example, in Quebec province, Canada, a partial lockdown was in effect for non-essential business, with many offices closed whereas primary and secondary schools open. The weekly school-related COVID-19 cases in Quebec at the end of August 2020 showed that schools accounted for 20% of the province’s COVID-19 cases, while students and staffs account for about 18% of Quebec’s population 4. Statistical data shows that 1,781 schools had been observed with at least one positive case in Quebec since the beginning of the pandemic5. Therefore, the rate of COVID-19 transmission in schools was higher than the community transmission, and mitigation measures must be implemented in classrooms to reduce the infection risk. Several studies revealed the significant impact of ventilation rate in reducing or preventing the airborne transmission of diseases in indoor environments 6. There are different recommendations for the minimum required ventilation rate in indoor spaces to achieve an acceptable indoor air quality or preventing indoor airborne transmission. The United States Centers for Disease Control and Prevention (CDC) and World Health Organization (WHO) recommended a minimum ventilation rate of 12 ACH to prevent airborne transmission in health-care facilities 7,8. The Harvard-CU Boulder Portable Air Cleaner Calculator 9 suggests a total of five air changes per hour as a good ventilation condition for reducing airborne transmission risk in classrooms. While recommendations are mainly based on the ventilation rate, it has been a challenge to quantify the outdoor air ventilation rate in a room. Indoor air CO2 concentration is often considered a surrogate/indicator for the ventilation rate. For example, the Montreal school board (Centre de services scolaire de Montreal) stated in an open letter on December 14, 2020: “Establishments without a mechanical ventilation system should apply the window opening guidelines to ensure frequent air changes in our premises”; “Always in order to ensure good indoor air quality, we have also started measuring carbon dioxide (CO2) in our establishments since November. In addition to this initiative, there are the CO2 tests that must be carried out by all school service centers in Quebec before December 16. The level of CO2 is a good indicator of the supply of fresh air in a room. Thus, following these tests, corrective measures will be put forward, if necessary.” 13. In the literature, several studies used a transient CO2 mass balance method and measured CO2 levels to calculate the ventilation rate in different indoor environments such as classrooms and university libraries 10–12. Batterman (2017) 12 estimated the CO2 generation rate based on the age and assumed activity level for CO2 calculation in mechanically ventilated classrooms, but the real activity type was unknown, so the results were subject to uncertainties. They used the whole day data to estimate ventilation rate, but they did not validate the model by calculating the CO2 at a different time or day. For the estimation of ventilation rate using the transient CO2 sensor data, it is essential to find dominant parameters that affect the final results, calibrate the model by measurement data, quantify and report the ventilation rates and infection risks with uncertainties. As a result, due to various factors affecting CO2 levels, such as variable occupant numbers and outdoor conditions, and the fact the unknown uncertainties of these factors, questions have been raised about using CO2 sensors to assess COVID-19 transmission risk 14. In summary, therefore, two questions are often raised: 1). The CO2 for Ventilation Assessment Question: when considering variable occupancy and dynamic surrounding environment, how can we relate CO2 concentrations to ventilation rates with uncertainties taken into account? 2). The CO2 for Risk Assessment Question: If CO2 is used as an indicator for the COVID-19 airborne aerosol transmission risk associated with ventilation conditions, what are the CO2 threshold levels to prevent spreading for ventilation-related risks with uncertainties quantified? Both questions center around the uncertainties and associated sensitivity analysis of parameters. Sensitivity Analysis (SA) and calibration methods such as Bayesian Markov Chain Monte Carlo (MCMC) method 15 can be used along with measured indoor CO2 concentrations to find the dominant parameters for the CO2 levels such as ventilation rate and calibrate them. Bayesian MCMC is a calibration technique proposed in the twentieth century owing to the development of MCMC and modern computer. Its application to the computer models’ calibration was systematically illustrated by Kennedy and O’Hagan (2001) 16. From then on, the boom of Bayesian inference and calibration was signified. Now, Bayesian inference calibration has been utilized in various topics, such as environment 17–20, hydrology 21–23, transportation 24, and medicine research 25. In the field of heating, ventilation, and air conditioning (HVAC), one of the early applications was conducted by Heo et al. 26. Building energy models were calibrated using Bayesian inference and further analyzed considering the uncertainties for retrofit decision making. By propagating parameters using probabilistic analysis, Bayesian inference incorporates uncertainties into real systems’ approximations by computer models. Combining multiple sources of information at different scales and with different reliabilities, the inadequacy of a model, revealed by the discrepancy between the predictions and observed data, can be corrected. Compared with conventional deterministic calibration methods, the advantages of Bayesian inference calibration are: 1) When the calibration measurements are qualitatively/quantitatively insufficient, for traditional methods, the estimated model parameters can be far off from their original values; however, for Bayesian calibration, since the uncertainties are considered, the calibration results are more stable and reasonable; 2) For the traditional calibration method, the results are often deterministic. In comparison, for the Bayesian inference calibration method, the results are derived from quantitative stochastic analysis and with possibilities that can be regarded as a degree of belief. Therefore, to answer the “CO2 for Ventilation Assessment” question, in this work, we investigated three classrooms from three schools in Montreal, Canada, to analyze CO2 and ventilation rate, and estimate COVID-19 airborne aerosol infection risk. To estimate ventilation rate using the transient CO2 mass balance model, we obtained the measurement data of transient CO2 concentrations. We conducted a sensitivity analysis to find the dominant parameters for indoor air CO2 concentration. The Bayesian MCMC method was then used with the three classrooms’ measured CO2 data to calibrate the dominant parameters and quantify the uncertainties. Calibrated ventilation rate helps teachers to compare it with the officially recommended value for preventing indoor airborne transmission. Then we used the calibrated model and tried different ventilation rates to estimate CO2 level for different exposure times. The results help occupants to estimate the ventilation rate from the measured CO2 data at each time. On the other hand, calculation and adjusting the ventilation rate may not be enough for mitigating the airborne transmission of COVID-19 in classrooms because no definite ventilation and CO2 thresholds have been agreed upon for COVID-19. Recommended ventilation rate or indoor air CO2 level by some standards and agencies for indoor air quality condition may not be enough to prevent the indoor airborne transmission. For example, ASHRAE Standard 62-2001 (ASHRAE, 2004) recommended a maximum CO2 level of 1000 ppm in classrooms for acceptable indoor air quality, which may still not be enough for preventing airborne transmission of diseases. The National Institute for Occupational Safety and Health (NIOSH) recommended a CO2 level of 600-1500 ppm for schools and workplaces but only considered comfort and working efficiency 27. An unofficial study by some scientists and teachers in 25 classrooms in Montreal revealed that in 75% of the classrooms, the level of CO2 exceeded 800 ppm because of poor ventilation conditions 28. However, they have been confused about the relations between monitored CO2 levels, ventilation conditions, and COVID-19 airborne transmission risks in the classrooms. The Wells-Riley model 29 can be used to calculate the indoor infection risk using calibrated ventilation rate. The effective reproductive ratio, R (ratio between secondary infectious cases and source cases), is often applied and should be less than one to prevent spreading. Du et al. (2020) studied the impact of ventilation improvement on a real tuberculosis (TB) outbreak in underventilated university buildings. Their result showed that improving indoor ventilation to levels correspond with CO2 < 1000 ppm reduced 97% of TB infection risk. Before improving the ventilation, the initial CO2 concentration was 3,204 ppm in one classroom where the index cases attended. They concluded that enhancing the ventilation rate to 14-15 ACH could effectively control the TB outbreak. However, there is no information given for the COVID-19 regarding the CO2 concentration for monitoring to prevent airborne aerosol spread. Rudnick and Milton (2003) derived an equation to estimate the indoor airborne infection transmission using indoor air CO2 concentration. They calculated the critical rebreathed fraction of indoor air below which airborne propagation of typical respiratory infections will not occur. They modeled several hypothetical cases in their work without measurement data. Peng and Jimenez (2020) 31 derived an analytical expression of CO2-based risk proxies for COVID-19 and used it to estimate CO2 level corresponding to an acceptable airborne risk level in different indoor environments. They showed that acceptable CO2 level varies by over two orders of magnitude for various rooms and activities. It also depends on other factors, such as wearing face masks. No measurements and uncertainties were reported. Eykelbosh (2021) 14 reviewed several studies that used indoor CO2 level in assessing the transmission risk and concluded that indoor air CO2 level could only represent the ventilation condition. The infection risk does not depend only on the ventilation rate, and other factors such as wearing a face mask, using portable air cleaner, and exposure time can also affect the infection risk. Therefore, to estimate the required ventilation rate and critical CO2 level to prevent the transmission of COVID-19 aerosols in a classroom, it is crucial to know the actual room condition such as occupancy profile, activity type, and other parameters that affect the estimation of infection risk. Therefore, to answer the “CO2 for Risk Assessment” question, in this work, we used the calibrated ventilation rate and actual room parameters to calculate the COVID-19 airborne transmission in the three classrooms using the modified Wells-Riley equation. We compared the results with infection risk corresponding to the Reproductive number to be one (R = 1), and different ventilation rates, and CO2 threshold levels at various exposure durations. ## 2. Methodology This section presents the models of CO2 concentration, airborne infection risk, and sensitivity analysis, and finally, the Bayesian calibration methods. Two fully-mixed transient mass balance equations are solved to calculate indoor air CO2 and COVID-19 quanta concentration. The schematic of the models is plotted in Figure 1. ![Fig. 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/03/2021.01.29.21250791/F1.medium.gif) [Fig. 1](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/F1) Fig. 1 Schematic of the indoor air CO2 and quanta concentration models ### 2.1 Indoor CO2 concentration model A fully-mixed transient mass balance model is solved for the calculation of CO2 concentration in the room. ![Formula][1] where *V* is the room volume ![Graphic][2] is the indoor air CO2 concentration (*mg*/*m*3); t is the time duration (*s*); *G**s* is the CO2 generation rate by all occupants (*mg*/*s*), which depends on the age and activity level; *λ*1 is the total outdoor air ventilation rate (*m*3/*s*); and *C**oa* is the outdoor air CO2 concentration (*mg*/*m*3). The transient mass balance, Eq. 1, is useful for solving arbitrary occupancy patterns and generation rates such as classrooms because students leave the classroom for break and launch. The solution of Eq. 1 is: ![Formula][3] where ![Graphic][4] is the observed initial CO2 concentration at each occupancy phase. ### 2.2 Airborne aerosol infection risk model The probability of infection (PI) of a susceptible person in the room is calculated using the Wells-Riley formulation32. The method was first used by Jimenez et al. 33 for calculating infection risk in different indoor environments and is recently applied to the City Reduced Probability of Infection (CityRPI) model and used for city-scale infection risk analysis 34,35. PI is a function of the number of quanta *μ* inhaled by the susceptible person (Eq. 3). We assumed that social distancing is maintained between all occupants, and the current study focuses on airborne aerosol transmission only. We used five assumptions for applying this model: i) there is only one infected person in the room who emits SARS-CoV-2 quanta with a constant rate, ii) the initial quanta concentration is zero, iii) the latent period of the disease is longer than the duration students stay in the classroom. Therefore the quanta emission rate remains constant during the day, iv) the indoor environment is well-mixed, and v) the infectious quanta is removed as a first-order process by the ventilation, filtration, deposition on surfaces, and airborne inactivation. The PI in Eq. 1 is based on the attendance of one infected person in the room, so it calculates the probability that COVID-19 aerosols are transmitted from the infected person to a susceptible person in the room; therefore, it is a conditional probability of infection (PIcond). ![Formula][5] The number of quanta inhaled by the susceptible person at the exposure time T is calculated by time-averaged quanta concentration. ![Formula][6] B is the inhalation rate (*m*3/*h*); *C**q,avg* is the time-average quanta concentration (*q*/*m*3); *T* is the exposure time (*h*); *f**m* is the fraction of people in the room who wears the mask, and *M**in* is the inhalation mask efficiency. A well-mixed transient mass balance equation similar to Eq. 1 is solved to calculate the room’s transient quanta concentration. ![Formula][7] where *C**q* is the indoor quanta concentration (*q*/*m*3); *E* is the net quanta emission rate (*h*−1); and *λ* is the first-order loss rate coefficient for quanta (*h*−1). Assuming that the initial quanta concentration is zero at the beginning of the day, Eq. 5 is solved as follows: ![Formula][8] Because of the change in the occupancy pattern during the day, the time-averaged quanta concentration is calculated using the Trapezoidal integration. *E* is calculated based on the number of infected people in the room *N**inf*, the fraction of people in the room with the mask *f**m*, exhalation mask efficiency *M**ex*, and quanta emission rate by one infected individual ER*q*. ![Formula][9] The first-order loss rate coefficient *λ* reflects several mechanisms: outdoor air ventilation *λ*1, filtration *λ*2, deposition on surfaces *λ*3, and airborne inactivation *λ*4. ![Formula][10] *λ*1 is the air change rate of outdoor air per hour (*h*−1) by the HVAC system or opening windows. *λ*2 is the in-room air filtration using portable air purifiers and/or duct filters in HVAC systems. *λ*3 is the removal by gravitational settling. Finally, *λ*4 is the inactivation/decay rate. As already mentioned, the *PI**cond* is the conditional probability of infection, assuming that there is one infected person in the room. The prevalence rate of the disease in the city *P**prev* can be used to estimate the number of possible infected individuals in the room. Therefore, the absolute probability of infection *PI**abs* is calculated using the *PI**cond*, *P**prev*, and number of susceptible people in the room *N**sus* 33: ![Formula][11] *N**sus* is calculated using the total number of people in the room *N**tot*, number of infected persons, and the fraction of immune people in the community *F**immune*. It can be estimated using the total recovered cases in the study region 36,37. ![Formula][12] The prevalence rate of the disease is estimated using the daily COVID-19 statistics reported by official sources: ![Formula][13] where *N**new* is the daily new cases, *P* is the population, *f**unrep* is the fraction of unreported cases. A study on ten diverse geographical sites in the US shows that the estimated number of infections was much greater than the number of reported cases in all sites 36. By doing more daily tests, the current unreported fraction is around 0.8. *D* is the duration of the infectious period of SARS-CoV-2 of 10 days in this work 38. ### 2.3 Sensitivity analysis model Both CO2 and infection risk models include uncertain parameters such as ventilation rates and emission rates. Some of the uncertain parameters may impact the result’s accuracy and should be calibrated by measurement data. Ideally, with sufficient measurements and computer resources, all the uncertain parameters should be included in the calibration parameters. In reality, limited by data quality/quantity or computer resources, only a few parameters could be included. Many parameters and inputs could also manifest different levels of uncertainties and significances on simulation outputs. So it is impracticable and unnecessary to calibrate all parameters, but for dominant parameters only. Identifying these dominant parameters cannot merely rely on arbitrary parameter selections from modelers’ knowledge but should be based on a scientific process, i.e., a sensitivity analysis. To conduct a sensitivity analysis process, prior distributions and ranges of selected unknown parameters should be determined. Then Monte Carlo (MC) simulation is employed to propagate simulations whose model parameters’ values are randomly chosen from the predefined ranges using a specific sampling method to perform simulation runs iteratively. Here, the Latin Hypercube Sampling (LHS) method 39 is applied since it provides good convergence of parameter space with relatively fewer samples. The obtained input-output dataset is then employed to identify the dominant model parameters that strongly affect the outputs. In this study, 440 parametric simulations were conducted for the sensitivity analysis of the CO2 concentration model 40. The importance ranking results may vary with different combinations of sensitivity methods and outputs depending on the variety of fundamental algorithms and conditions of each sensitivity analysis method 41. To avoid the potential inconsistency, Lim and Zhai proposed a new sensitivity analysis method, sensitivity value index (SVI), to account for the differences in sensitivity analysis methods and target outputs 42. Eq. 12 shows how SVI is applied to recognizing and comparing the importance rankings from different sensitivity analysis methods through the normalization and aggregation process. ![Formula][14] where *V**i,j* is the value of a sensitivity analysis method, *i* is a parameter, *n* is the total number of the parameters, *j* is a sensitivity method, *k* is the total number of sensitivity methods, *l* is the target output, and *m* is the total number of target outputs. ### 2.4 Bayesian calibration and Markov Chain Monte Carlo (MCMC) models As the footstone of all Bayesian statistics, Bayes’ theorem was first proposed by Reverend Thomas Bayes in his doctoral dissertation 43 and can be described as: ![Formula][15] The probability of an event is inferred based on the prior knowledge of conditions related to the event. Bayesian inference is one application of Bayes’ theorem and can be written as: ![Formula][16] where *p*(*θ*|*y*) is the posterior distribution of the unknown parameter *θ* based on known observation *y. p*(*y*|*θ*) is the likelihood function of observation conditional on the unknown parameter. *p*(*θ*) is the prior distribution of the unknown parameter which is the marginal probability that means it is irrespective of the outcome of another variable, and *p*(*y*) is the probability of the observation that is marginal as well to normalize *p*(*y*|*θ*). Therefore, the posterior probability is proportional to the product of the prior probability and the likelihood. In reality, it is impractical to apply the Bayesian inference for analytical solutions to all problems because the likelihood’s integrals can be computationally expensive or sometimes impossible to be calculated. MCMC is a versatile approach to solve the parameter estimation problem with two components. One is the well-known Monte Carlo method. It is a computational algorithm to solve statistically challenging problems relying on repeated random samplings and approximate the target value (e.g., mean value) using the independent samples’ results. The other is the Markov Chain method for solving a sequence of possible events. The probability of each event depends only on the state attained in the previous event. By combining MCMC and Bayesian inference, posterior distribution can be estimated efficiently. Different MCMC algorithms can be classified into either a “random walking” group or a gradient-based group according to the acceptance-rejection criterion’s adoption. In this study, Hamiltonian Monte Carlo (HMC) sampling method 44 was used for the MCMC. Five thousand steps of the HMC algorithms on each of four separate chains were explored in this study to make a total of 20,000 samplers. We used one thousand samples during the “warming-up” stage to move chains toward the highest density area and tune sampler hyperparameters. For each room, the first 2/3 of measurements are used for the calibration, with the remaining for the validation. In this study, the Coefficient of Variance of Root Mean Squared Error (CVRMSE) (Eq. 15) and Normalized Mean Bais Error (NMBE) were used as indicators to estimate the calibration and validation performance. ![Formula][17] where *ŷ**i* is a predicted variable value for period *i, y**i* is an observed value for period *i*, ![Graphic][18] is the mean of the observed value, and *n* is the sample size. ![Formula][19] ## 3. Case study In this study, three typical classrooms of three different schools in Montreal, Canada, were selected for calibration and infection risk analysis. Each classroom was monitored during a typical pandemic day, and occupants’ information (students’ age and number), ventilation system status, window status, and transient indoor CO2 concentration are recorded. The summary information is shown in Table 1 and Figure 2. View this table: [Table 1](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/T1) Table 1 Basic information about the measured classrooms ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/03/2021.01.29.21250791/F2/graphic-19.medium.gif) [](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/F2/graphic-19) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/03/2021.01.29.21250791/F2/graphic-20.medium.gif) [](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/F2/graphic-20) Fig. 2 CO2 concentration measurements, occupancy number profile, outdoor air temperature, and outdoor air pressure of a) Room 1, b) Room 2, c) Room 3. Room 1 is equipped with a mechanical ventilation system, and CO2 is between 500 ∼ 1000 ppm, while for Room 2 and Room 3, with no access to a mechanical ventilation system, the CO2 reaches up to 1800 ppm. The CO2 in Room 1 seems to indicate an acceptable level of air quality (<1000 ppm), and for Rooms 2 and 3, it is higher than the acceptable level of air quality requirements. The occupancy pattern of Room 1 and Room 2 was measured, but for Room 3, no occupancy pattern was recorded, so we assumed that number of students is constant as recorded in the morning during the day. The outdoor air temperature and pressure data are provided by Environment and Climate Change Canada 45. Other parameters required for the CO2 calculation, such as outdoor air CO2, generation rate, and outdoor air ventilation rate, are not available; therefore, we calibrate them using the CO2 measurements. ## 4. Results In this section, the sensitivity analysis finds the dominant parameters for the calculation of CO2 concentrations. Then we use the Bayesian MCMC calibration to estimate the daily average ventilation rate using the CO2 measurement data and occupancy patterns. With the calibrated model, we calculate the CO2 concentrations in the three rooms under different ventilation rates. Finally, we used the calibrated ventilation rate to estimate the infection risk in each classroom. We also find the ventilation rate and CO2 level thresholds to avoid the airborne COVID-19 aerosol spread for different exposure times. ### 4.1 CO2 model sensitivity analysis Outdoor/indoor pressure, outdoor/indoor air temperature, occupancy number, room volume, outdoor air ventilation rate, and CO2 exhale rate are input parameters to predict CO2 concentration. For the sensitivity analysis, the range of selected model inputs/parameters were defined according to the references, codes, and standards available. Table 2 shows the parameters with their sensitivity importance rankings: a smaller number indicates a more important/sensitive parameter. View this table: [Table 2](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/T2) Table 2 Sensitivity analysis with importance ranking for indoor CO2 concentration. It is concluded that, for the classroom CO2 levels, the most dominant parameters are outdoor air ventilation rate, CO2 generation rate per person, number of occupants, and outdoor CO2 concentration. Specifically, the outdoor air ventilation rate’s SVI is two times the second important parameter. Some sensitive parameters are often known, such as occupant number, outdoor temperature, and pressure. So they may not need to be calibrated. Therefore, we selected the outdoor air ventilation rate, CO2 generation rate, outdoor CO2 concentration, and indoor air temperature for the next step’s model calibration. Indoor pressure was assumed to be identical to the outdoor pressure. ### 4.3 Calibration and validation performance For the calibration of the CO2 model, the Bayesian inference method was applied. For each occupancy phase (e.g., between every two measurements), we use the new measured CO2 data as the initial condition for Eq. 2. The posterior distributions are shown in Figure 3 and Table 3. In each subplot, the red dash line represents the parameter’s prior distribution in Table 2. View this table: [Table 3](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/T3) Table 3 Calibrated parameters of the CO2 model ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/03/2021.01.29.21250791/F3/graphic-23.medium.gif) [](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/F3/graphic-23) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/03/2021.01.29.21250791/F3/graphic-24.medium.gif) [](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/F3/graphic-24) Fig. 3 Distribution of calibrated parameters of the indoor CO2 model. The posterior distributions with the Bayesian calibration are indicated by blue, green, and orange for Rooms 1, 2, and 3, respectively. The calibrated outdoor ventilation rate is 1.96 ± 0.31 ACH for Room 1, 0.40 ± 0.08 ACH for Room 2, and 0.79 ± 0.06 ACH for Room 3. Here, the ventilation rate is expressed by the calibrated mean value followed by the uncertainty for a 95% confidence interval. Room 1 is mechanically and naturally ventilated (i.e., open windows), so its ventilation rate is significantly higher than Rooms 2 and 3, with the outdoor air only from open windows. The results of Room 2 and Room 3 are closer since both are naturally ventilated with the same room volumes. The span of the posterior distribution of Room 1 is larger because of its wider prior distribution. For Room 3, a constant occupancy was used, so it is relatively easier for the MCMC iteration to find a posterior distribution. For the outdoor CO2 level and indoor air temperature, Room 2 and Room 3 are closer than Room 1 due to different ventilation modes. Using the mean value of the calibration parameters, we compared the simulation results and measurements of CO2 in Figure 4. The simulation results show similar trends well as the measurements. According to the American Society of Heating, Refrigerating, and Air-conditioning Engineers (ASHRAE) guideline 14 49 and FEMP 50, when the CVRMSE and NMBE values are less than 30% and ±10% for transient data, the calibrated accuracy meet the requirements. The (CVRMSE, NMBE) of validation of Rooms 1-3 are (15.3, 7.6), (10.5, 6.1), and (12.5, 12.3), respectively, which shows that the validated accuracy is satisfied. ![Fig. 4](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/03/2021.01.29.21250791/F4.medium.gif) [Fig. 4](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/F4) Fig. 4 Comparison of simulated and measured CO2 levels in three schools. ### 4.4 CO2 level and ventilation rate evaluations The calibrated ventilation rates in all three classrooms are all less than 2 ACH. The recommended ventilation rate for a safe indoor environment by Harvard-CU Boulder Portable Air Cleaner Calculator for Schools is at least 5 ACH. Therefore, the ventilation rate of all three classrooms seems inadequate. To relate the CO2 levels, e.g., measured from a CO2 sensor, and the ventilation rates, by using the calibrated CO2 model, we calculate the outdoor air ventilation rate in ACH and CFM/person as a function of CO2 levels at different exposure times in Figure 5. It helps teachers find the room ventilation rate directly based on the CO2 sensors at different school hours. For example, for Room 1, when CO2 > 600 ppm, OA (outdoor air) < 5 ACH (24 CFM/person); CO2 > 800 ppm, OA < 2 ACH (9 CFM/person) at any time of the day. A CO2 level less than 480 ppm indicates a ventilation rate greater than 10 ACH (48 CFM/person) at all times. For Room 2, the same CO2 levels correspond to a lower ventilation rate than Room 1. For example, when CO2 > 600 ppm, OA < 2 ACH (13 CFM/person); CO2 < 440 ppm indicates that OA > 10 ACH (66 CFM/person). This shows that Room 1 shows a higher ventilation rate (ACH) for the same CO2 level because of its smaller size. The ventilation rate of Room 3 at a specific CO2 level is higher than Room 2 because of its constant occupancy compared to the variable occupancy of Room 2. It seems lunch breaks indeed lower CO2 levels significantly (thus infectious risk in schools). For Room 3, CO2 > 600 ppm indicates OA < 3 ACH (22 CFM/person). These results show that the indoor CO2 could vary significantly in different classrooms even with the same ventilation rate because of different room sizes, occupants number, and occupancy schedule. All three classrooms show that an indoor CO2 lower than 450 ppm indicates a ventilation rate greater than 10 ACH, close to the recommended 12 ACH value to prevent airborne transmission in health-care facilities7,8. Note that here CO2 is the exposure-time-averaged instead of the instantaneous level. ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/03/2021.01.29.21250791/F5/graphic-26.medium.gif) [](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/F5/graphic-26) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/03/2021.01.29.21250791/F5/graphic-27.medium.gif) [](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/F5/graphic-27) Fig. 5 Relation of indoor CO2 levels and ventilation rates for classrooms. ### 4.5 CO2 level and infection risk evaluations To study CO2 and airborne aerosol infectious risk relation, we first calculate the probability of infection risk with the posterior distribution of the ventilation rate obtained in Section 4.3. Note, Here, the Bayesian and MCMC analysis allow us to quantify the uncertainties of the ventilation rates to estimate airborne infectious risk by defining the probability of the infection risk: the probable range of the infection risk estimated in classrooms. Then, we evaluate different ventilation rates to identify the corresponding CO2 threshold level, at which the reproductive number, R < 1, at all exposure times. We estimate the infection risk and indoor air CO2 threshold based on the actual room conditions in this work. The recommended threshold could be used for other rooms under a similar condition. Table 4 shows the input parameters for Eqs. 3-11 to calculate the COVID-19 airborne infection risk in classrooms. Actual room conditions with age and activity levels were used to determine breathing and quanta emission rates 51,52. All students wore a face mask in the classroom, and mask efficiency was selected based on students’ typical mask type. The prevalence rate and immunity fraction were calculated based on the official reported data53. View this table: [Table 4](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/T4) Table 4 Input parameters for calculation of infection risk in three classrooms Figure 6 plots the transient PIcond and PIabs in three classrooms with the calibrated ventilation rates’ posterior distribution. The orange/green line is the baseline of PIcond/PIabs calculated using the mean value of the posterior distribution of ventilation rate. The orange/green area is the estimated infection risk probability with a 95% confidence interval. For Rooms 1, 2, and 3, the mean *PI**cond* at the end of the day is around 14%, 14%, and 20%, respectively. Rooms 1 and 2 show lower *PI**cond* because no students were present during the break, so the quanta generation was zero. The uncertain band of *PI**cond* in Room 1 is wider than other rooms because of the larger posterior range of ventilation rate due to its mechanical ventilation system. In Rooms 1 and 2, when students leave the classroom, the *PI**abs* are zero because there is no susceptible persons in the room. Room 1 shows the lowest *PI**abs* because of the lower *PI**cond* and prevalence rate: all three measurements were conducted on different days. There are more students in Room 2 than Room 3; however, because of the lower *PI**cond*, the absolute *PI* of Room 2 is thus lower. ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/03/2021.01.29.21250791/F6/graphic-29.medium.gif) [](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/F6/graphic-29) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/03/2021.01.29.21250791/F6/graphic-30.medium.gif) [](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/F6/graphic-30) Fig. 6 Probability of airborne infection risks in classrooms. It also shows that after three hours, *PI**cond* of Room 1 is always smaller than two other rooms due to a higher ventilation rate. More frequent leaves of students (three breaks in a day) also contribute: the quanta level is generally lower than the two other rooms. *PI**cond* of Room 3 continues to rise with the exposure time, unlike other classrooms, again due to the constant occupancy assumed. Therefore, to reduce infection risk, more frequent breaks and leaves could be beneficial when increasing the classroom’s ventilation rate is relatively harder to achieve. The conditional *PI* is the ratio of the number of infections to susceptibles *PI**cond* = *D*/*S*. The reproductive number (R), defined by Rudnick and Milton (2003), is the number of secondary infections when a single infected person is introduced in the room, and everyone in the room is susceptible. If R < 1, then the infectious agent cannot spread in the population. For these three classrooms if *PI**cond* is smaller than 5.3%, 5%, and 5.5%, it is expected that the community spread in the classrooms could be stopped. Figure 6 shows that, for Rooms 1, 2, and 3, the conditional PI exceeds the level of R = 1 at around two ∼ three hours, after which mitigation measures should be taken. To identify the relation between the required ventilation rates, corresponding CO2 levels, and the COVID-19 airborne spreading risks, we calculated them for the three classrooms as shown in Figure 7. At all exposure times, the indoor CO2 and *PI* decrease with an increased ventilation rate. The ventilation rate threshold to prevent the spread (R<1) at all exposure times is 8, 3, and 6 ACH for Rooms 1, 2, and 3, respectively. Room 1 requires a higher ventilation rate because of the smaller room size. Room 2 needs a lower ventilation rate than Room 3 because students left the classroom several times. Therefore, the ventilation rate threshold depends on the occupancy schedule and the size of the room. ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/03/2021.01.29.21250791/F7/graphic-31.medium.gif) [](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/F7/graphic-31) ![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/03/2021.01.29.21250791/F7/graphic-32.medium.gif) [](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/F7/graphic-32) Fig. 7 *PI**cond* and time-averaged CO2 for different exposure times and ventilation rates. According to Figure 7, each classroom needs a different ventilation rate at different exposure times to reach the level of R = 1. On the other hand, we can find the indoor CO2 threshold for each exposure time corresponding *R*** = 1, as shown in Figure 8. The CO2 threshold level decreases with exposure time for all three rooms. For Room 1 at the second hour, the acceptable level is around 1000 ppm and later has to be 520 ppm at the end of the day to avoid spreading. For Rooms 2 and 3, it decreases from 660 ppm to 505 ppm and 750 ppm to 500 ppm, respectively. Meanwhile, the threshold levels for 3-7 hours exposure times are close for all three classrooms. In comparison, Figures 8b and 8c illustrate that the ventilation rate threshold (to prevent the spreading) increases with exposure time and is not a constant number because of the three rooms’ different sizes and schedules. ![Fig. 8](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/02/03/2021.01.29.21250791/F8.medium.gif) [Fig. 8](http://medrxiv.org/content/early/2021/02/03/2021.01.29.21250791/F8) Fig. 8 a) Time-averaged indoor CO2 concentration thresholds, and b, c) outdoor ventilation rate thresholds with exposure times to prevent spreading. In summary, the results of all three classrooms show that the ventilation rate threshold to prevent the airborne transmission of COVID-19 depends on several parameters such as room size, student schedules, and exposure time. The indoor CO2 threshold seems to depend on exposure time mostly, and the time-averaged level of 500 ppm seems acceptable for all three classrooms. ## 5. Conclusion The airborne transmission of COVID-19 is a major infection route in indoor spaces, especially with poor ventilation conditions, large occupancy density, and high exposure time, such as school classrooms. There are some recommended ventilation rates for acceptable indoor air quality or preventing airborne transmission in indoor spaces, but it is not easy to measure the actual room’s ventilation rate. Indoor air CO2 level can be used as an indicator for the ventilation rate, whereas it depends on several parameters that must be estimated. In this work, we used the sensitivity analysis, MCMC Bayesian Calibration method, measured data of indoor CO2 and occupancy profile of three classrooms in Montreal, Canada, to find and calibrate the dominant parameters to estimate CO2 levels indoors. The results showed that the outdoor ventilation rate is the most significant parameter. The calibrated ventilation rate with a 95% confidence level is 1.96 ± 0.31 ACH for Room 1 with mechanical ventilation, and 0.40 ± 0.08 ACH and 0.79 ± 0.06 ACH for Rooms 2 and 3 with windows open only. Based on the calibrated model, we created the correlations between the CO2 levels and the ventilation rates, which help teachers to estimate the ventilation rates from the CO2 sensors at any time. A time-averaged CO2 lower than 450 ppm is equivalent to a ventilation rate greater than 10 ACH in all three rooms, close to the recommended 12 ACH value for a safe indoor environment against airborne transmission. This study also proposed the approach to calculating the probability of the infection risk based on the calibrated ventilation rate, which helps quantify the uncertainty of outdoor ventilation rates. Moreover, this study estimated the required ventilation rate threshold and the CO2 threshold values to prevent the airborne aerosol spreading of COVID-19 as a function of exposure time in the classrooms. The ventilation threshold at all hours is 8, 3, and 6 ACH for rooms 1, 2, and 3, respectively, and the CO2 threshold is around 500 ppm at all exposure times (< 8 hr) of a school day for all three classrooms. This threshold is significantly different from the recommended value of 1000 ppm for acceptable indoor air quality conditions. Therefore, it is recommended that the ventilation rate and indoor CO2 concentration thresholds be reconsidered in indoor spaces, especially classrooms, in the current pandemic. ## Data Availability The data that support the findings of this study are available from the corresponding author, LW, upon reasonable request. ## Acknowledgment The research is supported by the Natural Sciences and Engineering Research Council (NSERC) of Canada through the Discovery Grants Program [#RGPIN-2018-06734] and the Advancing Climate Change Science in Canada Program [#ACCPJ 535986-18]. * Received January 29, 2021. * Revision received January 29, 2021. * Accepted February 3, 2021. * © 2021, Posted by Cold Spring Harbor Laboratory This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/) ## References 1. 1.Asadi S, Bouvier N, Wexler AS, Ristenpart WD. The coronavirus pandemic and aerosols: Does COVID-19 transmit via expiratory particles? Aerosol Sci Technol. 2020. 2. 2.Prather KA, Wang CC, Schooley RT. Reducing transmission of SARS-CoV-2. Science (80-). 2020. 3. 3.Zhang J. Integrating IAQ control strategies to reduce the risk of asymptomatic SARS CoV-2 infections in classrooms and open plan offices. Sci Technol Built Environ. 2020. 4. 4.Remiorz R. School COVID-19 cases consistent with community transmission, but Quebec could be in trouble, data show - The Globe and Mail. 2020. [https://www.theglobeandmail.com/canada/article-school-covid-19-cases-consistent-with-community-transmission-but/](https://www.theglobeandmail.com/canada/article-school-covid-19-cases-consistent-with-community-transmission-but/). Accessed January 13, 2021. 5. 5.Number of schools - Covid Écoles Québec. 2020. [https://www.covidecolesquebec.org/liste-alphabtique](https://www.covidecolesquebec.org/liste-alphabtique). Accessed January 13, 2021. 6. 6.Du CR, Wang SC, Yu MC, et al. Effect of ventilation improvement during a tuberculosis outbreak in underventilated university buildings. Indoor Air. 2020. 7. 7.Sehulster L, Chinn RYW. Guidelines for environmental infection control in health-care facilities. Recommendations of CDC and the Healthcare Infection Control Practices Advisory Committee (HICPAC). MMWR Recomm Rep. 2003. 8. 8.Atkinson J, Chartier Y, Lúcia Pessoa-Silva C, et al. Natural ventilation for infection control in health-care settings. World Heal Organ. 2016. 9. 9.Allen J, Cedeno-Laurent J, Miller S. Harvard-CU Boulder Portable Air Cleaner Calculator for Schools.v1.3 - Google Sheets. 2020. [https://docs.google.com/spreadsheets/d/1NEhk1IEdbEi\_b3wa6gI\_zNs8uBJjlSS-86d4b7bW098/edit#gid=1882881703](https://docs.google.com/spreadsheets/d/1NEhk1IEdbEi_b3wa6gI_zNs8uBJjlSS-86d4b7bW098/edit#gid=1882881703). Accessed January 13, 2021. 10. 10.Penman JM. An experimental determination of ventilation rate in occupied rooms using atmospheric carbon dioxide concentration. Build Environ. 1980. 11. 11.Mumovic D, Palmer J, Davies M, et al. Winter indoor air quality, thermal comfort and acoustic performance of newly built secondary schools in England. Build Environ. 2009. 12. 12.Batterman S. Review and extension of CO2-based methods to determine ventilation rates with application to school classrooms. Int J Environ Res Public Health. 2017. 13. 13.COVID-19: Measures in place as the holidays approach. 2020. [http://app.dialoginsight.com/T/OFSYS/SM2/365/2/S/F/1936/16667709/5P9DnIIH/734118.html](http://app.dialoginsight.com/T/OFSYS/SM2/365/2/S/F/1936/16667709/5P9DnIIH/734118.html). Accessed January 28, 2021. 14. 14.Eykelbosh A. Can CO2 Sensors be Used to Assess COVID-19 Transmission Risk? | National Collaborating Centre for Environmental Health | NCCEH - CCSNE. 2021. [https://ncceh.ca/content/blog/can-co2-sensors-be-used-assess-covid-19-transmission-risk](https://ncceh.ca/content/blog/can-co2-sensors-be-used-assess-covid-19-transmission-risk). Accessed January 22, 2021. 15. 15.Strawderman RL, Gamerman D. Markov Chain Monte Carlo: Stochastic Simulation for Bayesian Inference. J Am Stat Assoc. 2000. 16. 16.Kennedy MC, O’Hagan A. Bayesian calibration of computer models. J R Stat Soc Ser B Stat Methodol. 2001. 17. 17.Van Oijen M, Rougier J, Smith R. Bayesian calibration of process-based forest models: Bridging the gap between models and data. In: Tree Physiology. ; 2005. 18. 18.Steel D. Bayesian statistics in radiocarbon calibration. Philos Sci. 2001. 19. 19.Lehuger S, Gabrielle B, Oijen M van, et al. Bayesian calibration of the nitrous oxide emission module of an agro-ecosystem model. Agric Ecosyst Environ. 2009. 20. 20.Arhonditsis GB, Papantou D, Zhang W, Perhar G, Massos E, Shi M. Bayesian calibration of mechanistic aquatic biogeochemical models and benefits for environmental management. J Mar Syst. 2008. 21. 21.Kuczera G, Kavetski D, Renard B, Thyer M. A limited-memory acceleration strategy for MCMC sampling in hierarchical Bayesian calibration of hydrological models. Water Resour Res. 2010. 22. 22.Hall JW, Manning LJ, Hankin RKS. Bayesian calibration of a flood inundation model using spatial data. Water Resour Res. 2011. 23. 23.Xu T, Valocchi AJ. A Bayesian approach to improved calibration and prediction of groundwater models with structural error. Water Resour Res. 2015. 24. 24.Van Hinsbergen CPIJ, Van Lint HWC, Hoogendoorn SP, Van Zuylen HJ. Bayesian calibration of car-following models. In: IFAC Proceedings Volumes (IFAC-PapersOnline). ; 2009. 25. 25.Whyte S, Walsh C, Chilcott J. Bayesian calibration of a natural history model with application to a population model for colorectal cancer. Med Decis Mak. 2011. 26. 26.Heo Y, Choudhary R, Augenbroe GA. Calibration of building energy models for retrofit analysis under uncertainty. Energy Build. 2012. 27. 27.NIOSH. CDC - The National Institute for Occupational Safety and Health (NIOSH). 2020. [https://www.cdc.gov/niosh/](https://www.cdc.gov/niosh/). Accessed January 21, 2021. 28. 28.Wilton K. 12 Montreal teachers secretly tested classroom ventilation. The results are “problematic” | Montreal Gazette. 2020. [https://montrealgazette.com/news/local-news/12-montreal-teachers-secretly-tested-classroom-ventilation-the-results-are-problematic](https://montrealgazette.com/news/local-news/12-montreal-teachers-secretly-tested-classroom-ventilation-the-results-are-problematic). Accessed January 21, 2021. 29. 29.Barnhart S, Sheppard L, Beaudet N, Stover B, Balmes J. Tuberculosis in health care settings and the estimated benefits of engineering controls and respiratory protection. J Occup Environ Med. 1997. 30. 30.Rudnick SN, Milton DK. Risk of indoor airborne infection transmission estimated from carbon dioxide concentration. Indoor Air. 2003. 31. 31.Peng Z, Jimenez JL. Title: Exhaled CO 2 as COVID-19 infection risk proxy for different indoor environments and activities. medRxiv. 2020. 32. 32.Miller SL, Nazaroff WW, Jimenez JL, et al. Transmission of SARS-CoV-2 by inhalation of respiratory aerosol in the Skagit Valley Chorale superspreading event. Indoor Air. 2020. 33. 33.Jimenez JL. 2020\_COVID-19\_Aerosol\_Transmission\_Estimator - Google Sheets. 2020. [https://docs.google.com/spreadsheets/d/16K1OQkLD4BjgBdO8ePj6ytf-RpPMlJ6aXFg3PrIQBbQ/edit#gid=519189277](https://docs.google.com/spreadsheets/d/16K1OQkLD4BjgBdO8ePj6ytf-RpPMlJ6aXFg3PrIQBbQ/edit#gid=519189277). Accessed October 19, 2020. 34. 34.Katal A, Albettar M, Wang L (Leon). City Reduced Probability of Infection (CityRPI) for Indoor Airborne Transmission of SARS-CoV-2 and Urban Building Energy Impacts. medRxiv. 2021. 2021. [https://www.medrxiv.org/content/early/2021/01/20/2021.01.19.21250046](https://www.medrxiv.org/content/early/2021/01/20/2021.01.19.21250046). 35. 35.Wang L (Leon), Katal A, Albettar M. CityRPI for Indoor Airborne Transmission. 2020. [https://concordia-cityrpi.web.app/](https://concordia-cityrpi.web.app/). Accessed December 30, 2020. 36. 36.Havers FP, Reed C, Lim T, et al. Seroprevalence of Antibodies to SARS-CoV-2 in 10 Sites in the United States, March 23-May 12, 2020. JAMA Intern Med. 2020. 37. 37.Gu Y. Estimating True Infections | COVID-19 Projections Using Machine Learning. 2020. [https://covid19-projections.com/estimating-true-infections/](https://covid19-projections.com/estimating-true-infections/). Accessed October 12, 2020. 38. 38.Byrne AW, McEvoy D, Collins AB, et al. Inferred duration of infectious period of SARS-CoV-2: rapid scoping review and analysis of available evidence for asymptomatic and symptomatic COVID-19 cases. BMJ Open. 2020. 39. 39.Li Q, Augenbroe G, Brown J. Assessment of linear emulators in lightweight Bayesian calibration of dynamic building energy models for parameter estimation and performance prediction. Energy Build. 2016;124:194–202. 40. 40.Matala A. Sample Size Requirement for Monte Carlo simulations using Latin Hypercube Sampling. Indep Res Proj Appl Math. 2008:1–24. 41. 41.Menberg K, Heo Y, Choudhary R. Sensitivity analysis methods for building energy models: Comparing computational costs and extractable information. Energy Build. 2016. 42. 42.Lim H, Zhai ZJ. Comprehensive evaluation of the influence of meta-models on Bayesian calibration. Energy Build. 2017;155:66–75. 43. 43.Bayes TR. An Essay Towards Solving a Problem in the Doctrine of Chances. Philos Trans R Soc London. 1763;53:370–418. [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1098/rstl.1763.0053&link_type=DOI) 44. 44.Betancourt M. A Conceptual Introduction to Hamiltonian Monte Carlo. arXiv Prepr arXiv170102434. 2017. 45. 45.ECCC. Historical Data - Climate - Environment and Climate Change Canada. 2020. [https://climate.weather.gc.ca/historical\_data/search\_historic\_data\_e.html](https://climate.weather.gc.ca/historical\_data/search_historic_data_e.html). Accessed October 5, 2020. 46. 46.ANSI/ASHRAE. ANSI/ASHRAE Standard 62.1-2019, Ventilation for Acceptable Indoor Air Quality. ASHRAE. 2019. [https://ashrae.iwrapper.com/ASHRAE\_PREVIEW\_ONLY\_STANDARDS/STD\_62.1\_2019](https://ashrae.iwrapper.com/ASHRAE\_PREVIEW\_ONLY_STANDARDS/STD_62.1_2019). Accessed December 1, 2020. 47. 47.Persily A, de Jonge L. Carbon dioxide generation rates for building occupants. Indoor Air. 2017. 48. 48. Mike McGee. Daily CO2. A Pro Oxygen. 2016. [https://www.co2.earth/daily-co2](https://www.co2.earth/daily-co2). Accessed December 27, 2016. 49. 49.ANSI/ASHRAE. ASHRAE Guideline 14-2002 Measurement of Energy and Demand Savings. Ashrae. 2002. 50. 50.U.S. Department of Energy. M&V Guidelines: Measurement and Verification for Performance-Based Contracts Version 4.0. 2015. 51. 51.Buonanno G, Stabile L, Morawska L. Estimation of airborne viral emission: Quanta emission rate of SARS-CoV-2 for infection risk assessment. Environ Int. 2020. 52. 52.US Environmental Protection Agency. Exposure Factors Handbook: 2011 Edition. US Environ Prot Agency. 2011. 53. 53.Esri Canada. COVID-19 Canada. 2020. [https://resources-covid19canada.hub.arcgis.com/](https://resources-covid19canada.hub.arcgis.com/). Accessed January 5, 2021. [1]: /embed/graphic-2.gif [2]: /embed/inline-graphic-1.gif [3]: /embed/graphic-3.gif [4]: /embed/inline-graphic-2.gif [5]: /embed/graphic-4.gif [6]: /embed/graphic-5.gif [7]: /embed/graphic-6.gif [8]: /embed/graphic-7.gif [9]: /embed/graphic-8.gif [10]: /embed/graphic-9.gif [11]: /embed/graphic-10.gif [12]: /embed/graphic-11.gif [13]: /embed/graphic-12.gif [14]: /embed/graphic-13.gif [15]: /embed/graphic-14.gif [16]: /embed/graphic-15.gif [17]: /embed/graphic-16.gif [18]: /embed/inline-graphic-3.gif [19]: /embed/graphic-17.gif