Data-driven control of airborne infection risk and energy use in buildings
==========================================================================

* Michael J. Risbeck
* Alexander E. Cohen
* Jonathan D. Douglas
* Zhanhong Jiang
* Carlo Fanone
* Karen Bowes
* Jim Doughty
* Martin Turnbull
* Louis DiBerardinis
* Young M. Lee
* Martin Z. Bazant

## Abstract

The global devastation of the COVID-19 pandemic has led to calls for a revolution in heating, ventilation, and air conditioning (HVAC) systems to improve indoor air quality (IAQ) [1, 2, 3], due to the dominant role of airborne transmission in disease spread [4, 5, 6]. While simple guidelines have recently been suggested to improve IAQ mainly by increasing ventilation and filtration [7, 8], this goal must be achieved in an energy-efficient and economical manner and include all air cleaning mechanisms. Here, we develop a simple protocol to directly, quantitatively, and optimally control transmission risk while minimizing energy cost. We collect a large dataset of HVAC and IAQ measurements in buildings and show how models of infectious aerosol dynamics and HVAC operation can be combined with sensor data to predict transmission risk and energy consumption. Using this data, we also verify that a simple safety guideline is able to limit transmission risk in full data-driven simulations and thus may be used to guide public health policy. Our results provide a comprehensive framework for quantitative control of transmission risk using all available air cleaning mechanisms in an indoor space while minimizing energy costs to aid in the design and automated operation of healthy, energy-efficient buildings.

## Introduction

The COVID-19 pandemic disrupted the global economy and caused the worldwide shut-down of many public and private buildings essential for daily life, including schools, gyms, religious centers, and offices [9]. At first, public health guidance focused on limiting transmission from fomites and exhaled large droplets, via surface disinfection and social distancing (such as the 6 foot rule), respectively [10]. As the pandemic continued, however, it was recognized that a dominant mode of transmission of COVID-19 is through virus-laden exhaled aerosol droplets, which are small enough to remain suspended in the air for minutes to hours and become well-mixed across indoor rooms [4, 5, 6, 11, 12, 13], so the recommended mitigation strategies shifted from social distancing to masking and improved ventilation and filtration of indoor air [14, 15, 16]. Notably, a collection of leading experts in respiratory disease transmission and building science called for a “paradigm shift” in the design and operation of indoor air control systems, analogous to historical efforts to reduce pathogen transmission through food and water sources [2], as healthy indoor air is becoming recognized as a fundamental human need [1].

In this work, we propose a physics-based, data-driven strategy to achieve this paradigm shift in healthy buildings, which integrates mathematical models of airborne disease transmission with available building data streams, and apply it to data collected from multiple buildings and indoor space types. By combining airflow measurements from building management systems (BMS) with CO2 concentration data and other measures of indoor air quality (IAQ) obtained from portable sensors, the underlying physics-based models are calibrated and used to simulate the transmission risk and energy consumption for each indoor space, as operated. Using our workflow, public health officials can quantitatively assess mitigation strategies for indoor airborne disease transmission and recommend novel building control protocols to limit transmission while minimizing energy costs. Such quantitative analysis and automated building controls is necessary for preventing the next pandemic and minimizing widespread shutdowns during future pandemics.

### Experiment

We collected data from diverse indoor spaces on the Massachusetts Institute of Technology (MIT) campus for several weeks in April 2022. The monitored spaces are all served by heating, ventilation, and air-conditioning (HVAC) systems, which include sensors to measure and record total supply air flow, outdoor air flow (ventilation), and supply air temperature. Temporary in-room Kaiterra sensors were also deployed to collect additional measurements relevant to IAQ, including temperature, relative humidity, CO2 concentration, total volatile organic compounds (TVOC), and size-resolved particle concentrations. Only the first three measurements are used in this study, but the entire dataset is publicly available (SI 8). In addition, we collected information about each room, including floor area, ceiling height, use case, design occupancy, and HVAC configuration (SI 2).

The overall workflow for this study is illustrated in Fig. 1. The key idea is that all data streams are integrated with physics-based models to predict and control transmission risk and energy consumption. This data fusion ensures that all space-specific variables in the model are known or accurately estimated, thus allowing different spaces to be compared with confidence.

![Figure 1:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/20/2023.03.19.23287460/F1.medium.gif)

[Figure 1:](http://medrxiv.org/content/early/2023/03/20/2023.03.19.23287460/F1)

Figure 1: 
Diagram of overall workflow. Portable indoor air quality (IAQ) sensors are placed in each monitored room, and a more durable sensor is placed on the roof to record outdoor conditions. These sensor measurements are then combined with time-series data from the heating, ventilation, and air conditioning (HVAC) system and basic properties of each room (area, ceiling height, design occupancy, etc.) for analysis.

Example time series data for four of the monitored rooms are shown in Fig. 2 and exhibit variations on the scales of days, hours, and minutes. The CO2 and humidity measurements come from the in-zone IAQ sensors, while outdoor-air flow measurements come from the BMS or are estimated from CO2 measurements. Of these measurements, CO2 concentration is the most strongly varying, as it is driven primarily by room occupancy. The peak values for Classroom 3 are significantly higher than for the other rooms, as it is does not have a forced supply of outdoor air provided by the HVAC system and is thus only naturally ventilated.

![Figure 2:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/20/2023.03.19.23287460/F2.medium.gif)

[Figure 2:](http://medrxiv.org/content/early/2023/03/20/2023.03.19.23287460/F2)

Figure 2: 
Plot of sample time series data from the study period at multiple time scales (days, hours, and minutes). The shaded region in the first (or second) column shows the time range covered in the second (or third) column. Missing data points have been filled in via interpolation.

### Theory

Simple mass-balance models [17, 18] have been used for decades to study airborne transmission via aerosols and have been successfully applied to explain transmission in prior diseases [19, 20] and COVID-19 [21, 22, 23, 14, 24]. Each room is assumed to be well-mixed, such that the particle concentration can be treated as uniform throughout the room [17, 25]. This assumption is shown to produce high-quality occupancy estimates, indicating sufficient accuracy for our purposes (SI 1). The time-evolution of infectious pathogen concentration per droplet size can thus be expressed as a partial differential equation, which can be solved numerically (SI 3).

Following [14], the model can also be approximated analytically to derive a “safety guideline” that bounds the indoor reproductive number, ℛ*in*, defined as the expected number of transmissions if an infector were present for a time *τ* in a given room: ![Formula][1]</img>  where *C**q* represents the infectious quanta concentration in exhaled air (SI 5); *Q**b* is the occupants’ breathing rate; *p**m* is a mask penetration factor of aerosols; *V* is the volume of the room; *N**s* is the number of susceptible occupants; and *ϵ* is the desired risk tolerance. *λ*EOA, the “equivalent outdoor air” supply rate, lumps together all infectious particle removal mechanisms including ventilation, filtration, sedimentation, and disinfection [26]. EOA quantifies each removal mechanism in terms of volumetric flow of outdoor-air ventilation that would lead to an equivalent removal rate of infectious particles, thus facilitating comparisons among disparate processes. Droplet size dependencies are integrated out by defining an effective droplet size ![Graphic][2]</img> (SI 4). Typical ranges for these parameters and the values used in this study are provided in Supporting Information (SI 9).

Simply put, ℛ*in* is proportional to the product of susceptible occupants and the time spent, divided by the EOA provided to the room. Thus, any holistic risk assessment and transmission mitigation strategy must consider all three dimensions. For example, simply mandating a maximum occupancy in an indoor space may not adequately reduce transmissions if the occupants spend a large amount of time in the space. Alternatively, occupancy limits may be unnecessary, if an appropriate amount of EOA is delivered to the space.

CO2 is often used as an indicator of transmission risk [27, 28, 29, 8], but we stress that it is only a partial proxy that requires care to interpret, especially when comparing across spaces with significantly different HVAC systems. To derive a CO2-based guideline, a pseudo-steady analysis of the dynamical model for CO2 concentration can be combined with Equation (1) [28], which incorporates the differences between infectious particle dynamics and CO2 dynamics stemming from sources of EOA beyond ventilation (SI 5).

Short-range respiratory flows also contribute to the risk of disease transmission in indoor spaces [30, 31, 32] and should be compared with the risk of long-range airborne transmission in any safety guideline [14]. Here, we conservatively estimate that for most of the spaces monitored, short-range effects account for less than 5% of the expected transmissions caused by long-range mixing (Fig. S25), using an experimentally validated model [33] of respiratory turbulent plumes of warm exhaled air rising by natural convection [34]. The only exception is for rooms where close-range face-to-face contact between occupants is common, where short-range risk may reach up to 30% of the long-range risk. However, the flows responsible for short-range transmission can be eliminated by requiring occupants to wear masks [34].

## Results

### Infection risk

The workflow can be assessed by calculating the transmission risk in the different indoor spaces using the collected data along with the full dynamical model and the pseudo-steady approximation. A key factor affecting the transmission rate is the time-varying occupancy in each space, which is then used to estimate the numbers of susceptible, *N**s*, and infectious, *N**i*, occupants. There are various occupant-counting technologies available, which utilize images, sound, or IAQ data, such as CO2 concentration, temperature, and humidity [35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48]. We estimate occupancy by modifying previous methods that solve the full inverse problem of the CO2 concentration dynamical model and data, with a novel extension to simultaneously estimate ventilation rates (SI 1) if not measured. This method protects occupant privacy and requires no training data, so it can be easily implemented in new spaces.

The transmission rate values predicted by the full model and pseudo-steady approximation are in excellent agreement, and where they differ, the pseudo-steady model produces more conservative estimates (Fig. 3). We assess these models in two ways. In Fig. 3, we extract random segments of data from Classroom 2 (during nominally occupied hours) and plot points for each segment on a plane with axes for average occupancy and time. The color of each point corresponds to the event reproductive number calculated from the full dynamical model. Since EOA delivery is essentially constant for this space, the pseudo-steady model predicts that ℛ*in* ∝ Occupancy × Time, which is the same trend exhibited by the color of the points in Fig. 3. This validates the use of the safety guideline, Eq. (1), to limit occupancy and time such that ℛ*in* is below a given tolerance, as indicated by the dashed curve, which can be shifted by altering the amount of EOA provided in the room.

![Figure 3:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/20/2023.03.19.23287460/F3.medium.gif)

[Figure 3:](http://medrxiv.org/content/early/2023/03/20/2023.03.19.23287460/F3)

Figure 3: 
Transmission predictions for selected rooms in the study. Left: Scatter plot of the reproductive number for randomly chosen time periods from the full dynamical simulations, found to be in good agreement with the safety guideline [14] from the pseudo-steady formula, Eq. (1). Right: Distributions of transmission rates throughout the study period computed using full model simulation and pseudo-steady approximation, again showing good agreement. The distributions are weighted by occupancy and thus predict the expected number of transmissions if one occupant were to be infectious for one hour. Green curves are the distributions from the full dynamic model, with green lines showing the minimum, median, and maximum values. Gray curves are the distributions from the pseudo-steady model, with grey lines showing the same three statistics. Black dots show expected values assuming minimum ventilation rates and occupant density per ASHRAE standard 62.1 [49].

To assess other spaces, Fig. 3 shows distributions of transmission rates across all the monitored spaces. These points are based on a 5-minute sample rate for both models, and the distributions are weighted by the number of occupants within each point. The black dots indicate the corresponding steady-state transmission rate for a space of that type with baseline ventilation rates and occupant density, per ASHRAE standards. We see that, in almost all spaces, the worst-case transmission rate is below the baseline value as expected, since MIT buildings were deliberately operated with extra ventilation during the monitoring period to limit COVID-19 transmission. Again, the distribution of transmission rates calculated from the pseudo-steady model (gray curves) closely matches the distribution of transmission rates calculated from the full dynamic model (green curves). The median values of these distributions are generally within 1% agreement, and the maximum values differ by less that 10%.

The primary outlier from this trend is Classroom 3*, which has no mechanical ventilation, resulting in significantly smaller EOA delivery than in the other spaces. As a result, when a large number of people enter this room, the infectious particle concentration takes longer to approach the pseudo-steady values. As such, the pseudosteady model predicts a more conservative, higher transmission rate. When a large number of people leave those spaces, the opposite transient effect occurs, but since the occupancy is generally lower as people are leaving the space, those events are weighted less in the final distribution. The net result is that the pseudo-steady model predictions are slightly conservative.

The results shown for Classrooms 3 and 3* are for the same space, occupancy profile, and ventilation rate, but different levels of filtration. We assume an active filter delivering 4.5 ACH of HEPA filtration in Classroom 3 and inactive filter in Classrom 3*. When the filter is inactive, there is roughly a threefold increase in median transmission rate. We highlight this distinction because measured CO2 concentrations would be exactly the same for the two scenarios, since filtration is a form of EOA that does not impact CO2, which thus cannot be used by itself to assess transmission risk. Instead, safety guidelines [14, 28], which incorporate the differences between total EOA and ventilation (SI 5), based on a fixed ℛ*in* tolerance, are more appropriate to account for differences in indoor spaces.

### Energy and control analysis

While addressing public health concerns, there are still many opportunities to reduce energy consumption in buildings, which account for 18% of total energy use in the United States [50]. Long term building operation must balance airborne transmission risk and IAQ with energy consumption. Specifically, given a desired level of total transmission risk, buildings should operate to achieve that risk as efficiently as possible, taking advantage of all available mechanisms of infectious particle mitigation.

Formulating all removal processes in terms of EOA provides a common basis to compare various technologies in terms of cost per ACH of EOA. Under this lens, filtration, which can either be provided by in-room air cleaners or recirculated supply air, is often a much more energy-efficient source of EOA than ventilation. As an alternative, ultraviolet (UV) light can provide significant EOA by eradicating any infectious material within particles [51, 52, 53, 54] and can be installed in an upper-room configuration with shielding or as “far-UV” that is not harmful to occupants [55]. If properly installed, such systems can deliver EOA even more efficiently than filtration-based sources [56]. Thus, to optimize energy efficiency, it is necessary to consider all of these options.

Unfortunately, the primary source of EOA for many rooms is in fact ventilation. Considering the large variance in transmission risk for spaces with strongly time-varying occupancy, however, it is possible to significantly reduce energy costs by limiting extra ventilation to periods of high occupancy. Here, we analyze and compare demand-controlled ventilation (DCV) and transmission-controlled ventilation (TCV) operation modes. DCV is a feedback control mechanism implemented in most modern HVAC systems, which adjusts ventilation rates in real time to maintain a setpoint of CO2 concentration. This method thus does not consider any other EOA sources, and thus the mapping from CO2 setpoint to transmission risk can vary strongly from space to space. By contrast, we propose TCV as a novel operating mode to maintain a transmission-rate setpoint by interfacing with HVAC and accounting for other sources of EOA, all of which impact airborne disease transmission.

To quantitatively assess the inherent tradeoff between energy consumption and transmission risk, we estimate the daily energy cost for each room as operated and under various hypothetical scenarios using standard thermodynamic and equipment modeling procedures from our previous work [57, 58]. For the monitored rooms, the primary cost driver is the energy required to heat the outdoor air up to its supply temperature.

We refer to the actual operation during the monitoring period as the “Baseline” scenario. The hypothetical scenarios considered for each room are as follows:

*   Curtailed: ventilation is supplied at the same rate as observed in the data only during nominal occupied hours, assumed to be 8 am through 10 pm.

*   In-Zone Filtration: in addition to the schedule change in the “Curtailed” scenario, 2 ACH of in-room filtration is provided via standalone air cleaners active during occupied hours.

*   In-Zone Far UV: in addition to the schedule change in the “Curtailed” scenario, upper-room far UVC lamps are activated during occupied hours so as to provide 5 ACH of EOA.

*   ASHRAE Minimum: ventilation follows the “Curtailed” schedule and is further adjusted to provide the minimum amount of ventilation required in each space per ASHRAE standard 62.1 [49].

*   Demand Controlled (DCV): ventilation is provided by a standard demand control algorithm for a given CO2 concentration setpoint.

*   Transmission Controlled (TCV): ventilation is provided by a modified algorithm that provides enough ventilation to operate below a given transmission risk as calculated by the pseudo-steady model.

After estimating the time-varying ventilation that would be provided by each hypothetical strategy, transmission risk and energy cost can be calculated using the modeling approach discussed previously. These results are shown for three representative spaces in Fig. 4. Similar plots for other spaces are provided in the Supporting Information.

![Figure 4:](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2023/03/20/2023.03.19.23287460/F4.medium.gif)

[Figure 4:](http://medrxiv.org/content/early/2023/03/20/2023.03.19.23287460/F4)

Figure 4: 
Summary of energy versus transmission rate tradeoffs for hypothetical ventilation scenarios. Points and shaded regions show means and joint standard deviations for daily values within the study period.

In all spaces in Fig. 4, we see there are significant opportunities to reduce energy consumption without large changes in the average and spread of the transmission rate. Simply curtailing ventilation during nighttime unoccupied hours cuts energy consumption roughly in half, with only a slight increase in transmission rate due to a small number of after-hours gatherings in that room. Adding in-room filtration can reduce transmission below baseline values while still providing significant reduction in energy cost. Based on current experimental data [55], far UV disinfection is even cheaper than in-room filtration. The ASHRAE minimum protocol further reduces energy consumption but increases the mean and spread of the transmission rate. However, by applying some of the more advanced control algorithms, transmission rate can be maintained near or below a desired threshold while maximizing energy savings. In particular, DCV at 800 ppm achieves minimum energy cost with transmission risk, while the novel TCV strategy at 0.05 per infector·h forgoes some of the energy savings to achieve further reduction in observed transmission rate. Note that the specific setpoints of these two strategies could be adjusted up or down to further tune the tradeoff. Finally, combining TCV with in-room far UV can achieve the same average and lower spread of transmission rate as the conservative baseline schedule with up to a tenfold decrease in energy cost.

Overall, these results illustrate that advanced control strategies and alternative sources of EOA can be employed to provide similar or better expected transmission rate while significantly reducing energy costs compared to constantly operating at high ventilation rates. We see that the TCV strategies deliver Pareto-optimal performance, which demonstrates that the pseudo-steady transmission model is sufficiently accurate to achieve its control objectives while remaining mathematically simple enough to integrate into existing HVAC control logic. Such TCV systems could interface with other sources of EOA, such as filtration and UV disinfection, to simultaneously control transmission rates and minimize energy consumption by prioritizing lower energy sources.

## Conclusion

Our framework successfully combines data-streams from sensors with accurate physical models of infectious aerosols to predict airborne disease transmission rates in diverse indoor spaces. We have demonstrated the possibility of transmission-controlled ventilation by implementing our models in HVAC control logic. We have also provided a design framework for future buildings, which quantifies the tradeoff between different modes of clean air delivery.

Our results can also inform public health guidance, using data from real buildings. We have validated the use of the simple safety guideline, Eq. (1), to limit infection risk in different classes of indoor spaces, rather than strict occupancy limits [14]. We have shown that the underlying pseudo-steady approximation is consistent with full, dynamical simulations, and whenever small discrepancies arise, the guideline always provides a more conservative estimate of the risk. For normal occupancy in the monitored spaces, we also predict that short-range transmission via respiratory flows can be neglected (compared to the long-range airborne transmission) without imposing physical distance limits.

Controlling disease transmission must also be considered within the context of broader societal needs, such as minimizing energy usage, pollution, and carbon emissions. Our framework is able to quantify and optimize these tradeoffs, as the various sources of EOA are all incorporated, including their different energy requirements. Enabling real-time control of infection risk, prioritized against energy consumption, is a critical first step in the paradigm shift toward more healthy, energy efficient buildings [2].

## Supporting information

Supplementary Information [[supplements/287460_file02.pdf]](pending:yes)

## Data Availability

All data produced are available online at [https://github.com/acoh64/MIT-JCI-IAQ-HVAC](https://github.com/acoh64/MIT-JCI-IAQ-HVAC)

[https://github.com/acoh64/MIT-JCI-IAQ-HVAC](https://github.com/acoh64/MIT-JCI-IAQ-HVAC) 

## Methods

### Occupancy Estimation

A key input to the transmission-rate model is the (time-varying) occupancy in each space, which is then used to estimate the number of susceptible, *N**s*, and infectious, *N**i*, occupants. Our proposed approach to estimate occupancy and ventilation rates proceeds via parameter estimation applied to the following ODE model: ![Formula][3]</img>  The state variable is the CO2 concentration ![Graphic][4]</img>, while *V* is the space volume, *Q**b* ≈ 0.6 m3/h is the occupant breathing rate, ![Graphic][5]</img> ppm is the exhaled-breath excess CO2 concentration, ![Graphic][6]</img> ppm is the outdoor-air CO2 concentration, and *k**a*(*t*) = *Q**a*(*t*)/*V* is the (possibly time-varying) ventilation rate. The values to be inferred are the timevarying occupancy *N**t*(*t*) and also the ventilation rate *k**a* if it is not measured. A key benefit of the proposed approach is that it requires no prior training or actual occupancy counts, and the room volume is the only parameter that has to be specified. Exact details on how we solve the inverse problem to infer *N**t*(*t*) and *k**a* are provided in SI 1.

### Full Dynamical Model

The mass balance for infectious particles in a room results in the following partial differential equation model for the time-evolution of infectious pathogen concentration, *C*(*r, t*), per droplet size in a room of volume *V* and area *A* [14]: ![Formula][7]</img>  where *N**i* are the number of infectors present in the room exhaling infectious droplets with rate *P*. Infectious droplets are removed through ventilation (*Q**a*), filtrationin the recirculated airflow (*p**f* (*r*)*Q**r*), sedimentation (*υ**s* (*r*)*A*), deactivation (*λ**υ* (*r*)*V*), and the action of disinfection devices (∑*d* *p**d*(*r*)*Q**d*). All removal mechanisms can be expressed as rates, *λ**a* = *Q**a*/*V, λ**f* (*r*) = *p**f* (*r*)*Q**r* /*V, λ**s* (*r*) = *υ**s* (*r*)*A*/*V*, and *λ**d* (*r*) = ∑*d* *p**d* (*r*)*Q**d* /*V*, and lumped into a single parameter the describes the supply of “equivalent outdoor air” (EOA) delivered to the space, *λ*EOA = *λ**a* + *λ**f* + *λ**s* + *λ**υ* + *λ**d* [26].

The production rate, *P*, an be further expressed as *P*(*r*) = *Q**b**n**d*(*r*) *V**d*(*r*) *p**m*(*r*) *c**υ*(*r*), where *Q**b* is the breathing flow rate of the individuals; *n**d*(*r*) is the number density of pathogens per volume of breath, which is known to vary with factors that include respiratory activity, time since infection, etc.; *V**d*(*r*) is the volume of the aerosol droplets; *p**m* is a mask penetration factor which accounts for the proportion of pathogen that may be filtered out by the mask (where a value of 1 means all pathogen escapes the mask and 0 means all pathogen is filtered by the mask) [59]; and *c**υ*(*r*) is the pathogen concentration in the droplets.

According to this model, the steady-state value of pathogen concentration if one infector is present is ![Graphic][8]</img>.

### CO2-based Safety Guideline

The CO2-based safety guideline can be derived by considering the steady-state CO2 concentration to the dynamical model given in Equation (2). The steady-state solution is ![Formula][9]</img>  We can rearrange Equation (4) and combine with Equation (1) to arrive at the CO2-based safety guideline ![Formula][10]</img>  where we have assumed that *N**t* /*N**s* *N**i* ≈ 1.

### Short-Range Transmission Risk

Estimates of short-range transmission rates can be derived from the theory of turbulent jets [14]. This analysis predicts that the concentration of infectious particles in the jets of infectors’ exhaled breath decays as 1/*x* where *x* is horizontal distance. A key deficiency of this model is that it does not account for the buoyancy of exhaled breath, which causes it to quickly rise out of the breathing zone of a potential susceptible. Thus, rather than use the turbulent jet models directly, we instead employ an empirical model derived from the experimental results of [33]. In that paper, the authors calculate a “susceptibility index” defined as *ϵ* := *C*(*x*)/*C*∞ where *C*(*x*) is the infectious particle concentration at horizontal distance *x* from the mouth of the infector, and *C*∞ is the background room concentration. For our purposes, we assume the concentration within the jet follows the model ![Formula][11]</img>  where *k* is an unknown constant to be determined. Assuming *C*∞ = *C**b**Q**b*/*V λ*EOA follows the pseudo-steady well-mixed model, we can derive the relationship ![Formula][12]</img>  To quantify the short-range transmission rate, we use the model ![Formula][13]</img>  in which the new parameter *p*short represents the probability that a susceptible is directly within the short-range plume exhaled by each infector. More information about how we determine this parameter is provided in SI 6.

### Transmission-Controlled Ventilation

Given that our primary goal is to control the transmission risk in each room, we propose using the reproductive number directly as a controlled variable, rather than using CO2 as in DCV.

To implement this control strategy, we of course first need to evaluate the current transmission rate ![Graphic][14]</img>. The pseudo-steady model gives ![Graphic][15]</img>, in which we have removed some extra factors for brevity. The value of *λ*EOA*V* can be calculated using flow measurements and filtration parameters for the BMS-provided clean air and the humidity measurements and physics-based models for the deposition and deactivation components of EOA. To estimate *N**s*, the CO2 generation rate ![Graphic][16]</img> can be calculated from successive measurements of ![Graphic][17]</img> and *Q**a* = *λ**a**V*. We could then calculate *N**i* := max(*N**t* − 1, 0), although we propose using *N**i* ≈ *N**t* to add a slight degree of robustness for small rooms. We thus arrive at the formula ![Formula][18]</img>  which can be evaluated by the BMS.

To define the action of the controller, we thus take a transmission-rate setpoint ![Graphic][19]</img> and invert the previous formula to find the corresponding EOA setpoint ![Formula][20]</img>  From this value, the BMS can adjust its various setpoints to deliver the required amount of EOA. In cases where the BMS can control multiple sources of EOA (e.g., ventilation, filtration via recirculation, and possibly in-zone disinfection devices), some form of prioritization would be needed, for example selecting in order of increasing energy consumption. More information on the corresponding control logic is provided in SI 7.

## Data availability

All data used in this study is available at [https://github.com/acoh64/MIT-JCI-IAQ-HVAC](https://github.com/acoh64/MIT-JCI-IAQ-HVAC).

## Code availability

Please contact M.J.R. at michael.james.risbeck{at}jci.com for code.

## Author Contributions

M.Z.B. and Y.M.L. conceived the study. M.J.R. wrote the code and performed the simulations. M.J.R., J.D.D., A.E.C., Z.J., C.F., K.B., J.D., M.T., and L.B. set up the sensors and planned experiments. M.J.R. and A.E.C. prepared the manuscript. All authors discussed the results and edited the manuscript.

## Competing interests

The authors declare no competing interests.

## Supplementary Information

Supplementary Information is available for this work.

## Acknowledgments

A.E.C. was supported by the Department of Defense (DoD) through the National Defense Science and Engineering Graduate (NDSEG) Fellowship Program. Partial funding for this study was provided by Johnson Controls, and portions of the research are subject to patents pending.

*   Received March 19, 2023.
*   Revision received March 19, 2023.
*   Accepted March 20, 2023.


*   © 2023, 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].Morawska, L., Marks, G. B. & Monty, J. Healthy indoor air is our fundamental need: the time to act is now. The Medical Journal of Australia 217, 578–581 (2022).
    
    

2.  [2].Morawska, L. et al. A paradigm shift to combat indoor respiratory infection. Science 372, 689–691 (2021).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNzIvNjU0My82ODkiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMy8wMy8yMC8yMDIzLjAzLjE5LjIzMjg3NDYwLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

3.  [3].Nwanaji-Enwerem, J. C., Allen, J. G. & Beamer, P. I. Another invisible enemy indoors: Covid-19, human health, the home, and united states indoor air policy. Journal of Exposure Science & Environmental Epidemiology 30, 773–775 (2020).
    
    

4.  [4].Morawska, L. & Milton, D. K. It Is Time to Address Airborne Transmission of Coronavirus Disease 2019 (COVID-19). Clinical Infectious Diseases 71, 2311–2313 (2020)…
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/cid/ciaa939&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F20%2F2023.03.19.23287460.atom) 

5.  [5].Morawska, L. et al. COVID-19 and Airborne Transmission: Science Rejected, Lives Lost. Can Society Do Better? Clinical Infectious Diseases (2023)…
    
    

6.  [6].Sachs, J. D. et al. The lancet commission on lessons for the future from the covid-19 pandemic. The Lancet 400, 1224–1280 (2022).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(22)01585-9&link_type=DOI) 

7.  [7].Allen, J. G. & Ibrahim, A. M. Indoor air changes and potential implications for sars-cov-2 transmission. Journal of the American Medical Association 325, 2112–2113 (2021).
    
    

8.  [8].Raymenants, J. et al. Indoor air surveillance and factors associated with respiratory pathogen detection in community settings in belgium. Nature Communications 14, 1332 (2023).
    
    

9.  [9].Pak, A. et al. Economic consequences of the covid-19 outbreak: the need for epidemic preparedness. Frontiers in Public Health 8, 241 (2020).
    
    

10. [10].Van Doremalen, N. et al. Aerosol and surface stability of sars-cov-2 as compared with sars-cov-1. New England Journal of Medicine 382, 1564–1567 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1056/NEJMc2004973&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F20%2F2023.03.19.23287460.atom) 

11. [11].Zhang, X. et al. Monitoring sars-cov-2 in air and on surfaces and estimating infection risk in buildings and buses on a university campus. Journal of Exposure Science & Environmental Epidemiology 32, 751–758 (2022).
    
    

12. [12].Chia, P. Y. et al. Detection of air and surface contamination by sars-cov-2 in hospital rooms of infected patients. Nature communications 11, 2800 (2020).
    
    

13. [13].Nissen, K. et al. Long-distance airborne dispersal of sars-cov-2 in covid-19 wards. Scientific reports 10, 1–9 (2020).
    
    

14. [14].Bazant, M. Z. & Bush, J. W. A guideline to limit indoor airborne transmission of covid-19. Proceedings of the National Academy of Sciences 118, e2018995118 (2021).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxODoiMTE4LzE3L2UyMDE4OTk1MTE4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDMvMjAvMjAyMy4wMy4xOS4yMzI4NzQ2MC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

15. [15].Blocken, B. et al. Ventilation and air cleaning to limit aerosol particle concentrations in a gym during the covid-19 pandemic. Building and Environment 193, 107659 (2021).
    
    

16. [16].Lindsley, W. G. et al. Efficacy of portable air cleaners and masking for reducing indoor exposure to simulated exhaled sars-cov-2 aerosols—united states, 2021. Morbidity and Mortality Weekly Report 70, 972 (2021).
    
    

17. [17].Wells, W. F. Airborne contagion and air hygiene: an ecological study of droplet infections (Harvard University Press, 1955).
    
    

18. [18].Riley, E., Murphy, G. & Riley, R. Airborne spread of measles in a suburban elementary school. American journal of epidemiology 107, 421–432 (1978).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/oxfordjournals.aje.a112560&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=665658&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F20%2F2023.03.19.23287460.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1978EZ85600006&link_type=ISI) 

19. [19].Beggs, C., Noakes, C., Sleigh, P., Fletcher, L. & Siddiqi, K. The transmission of tuberculosis in confined spaces: an analytical review of alternative epidemiological models. The International Journal of Tuberculosis and Lung Disease 7, 1015–1026 (2003).
    
    

20. [20].Nicas, M., Nazaroff, W. W. & Hubbard, A. Toward understanding the risk of secondary airborne infection: emission of respirable pathogens. Journal of Occupational and Environmental Hygiene 2, 143–154 (2005).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/15459620590918466&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15764538&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F20%2F2023.03.19.23287460.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000227927400007&link_type=ISI) 

21. [21].Buonanno, G., Morawska, L. & Stabile, L. Quantitative assessment of the risk of airborne transmission of sars-cov-2 infection: prospective and retrospective applications. Environment international 145, 106112 (2020).
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F20%2F2023.03.19.23287460.atom) 

22. [22].Augenbraun, B. L. et al. Assessment and mitigation of aerosol airborne sars-cov-2 transmission in laboratory and office environments. Journal of Occupational and Environmental Hygiene 17, 447–456 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/15459624.2020.1805117&link_type=DOI) 

23. [23].Miller, S. L. et al. Transmission of sars-cov-2 by inhalation of respiratory aerosol in the skagit valley chorale superspreading event. Indoor Air 31, 314–323 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/ina.12751&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F20%2F2023.03.19.23287460.atom) 

24. [24].Wang, C. C. et al. Airborne transmission of respiratory viruses. Science 373 (2021).
    
    

25. [25].Foster, A. & Kinzel, M. Estimating covid-19 exposure in a classroom setting: A comparison between mathematical and numerical models. Physics of Fluids 33, 021904 (2021).
    
    

26. [26].Force, A. E. T. Building readiness. Tech. Rep. (2022). Accessed: 2022-11-07.
    
    

27. [27].Rudnick, S. & Milton, D. K. Risk of indoor airborne infection transmission estimated from carbon dioxide concentration. Indoor air 13, 237–245 (2003).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1034/j.1600-0668.2003.00189.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12950586&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F20%2F2023.03.19.23287460.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000185114700005&link_type=ISI) 

28. [28].Bazant, M. Z. et al. Monitoring carbon dioxide to quantify the risk of indoor airborne transmission of covid-19. Flow 1, E10 (2021).
    
    

29. [29].Peng, Z. & Jimenez, J. L. Exhaled CO2 as a covid-19 infection risk proxy for different indoor environments and activities. Environmental Science & Technology Letters 8, 392–397 (2021).
    
    

30. [30].Yang, F., Pahlavan, A. A., Mendez, S., Abkarian, M. & Stone, H. A. Towards improved social distancing guidelines: Space and time dependence of virus transmission from speech-driven aerosol transport between two individuals. Physical Review Fluids 5, 122501 (2020).
    
    

31. [31].Abkarian, M., Mendez, S., Xue, N., Yang, F. & Stone, H. A. Speech can produce jet-like transport relevant to asymptomatic spreading of virus. Proceedings of the National Academy of Sciences 117, 25237–25245 (2020).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTE3LzQxLzI1MjM3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjMvMDMvMjAvMjAyMy4wMy4xOS4yMzI4NzQ2MC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

32. [32].Bourouiba, L., Dehandschoewercker, E. & Bush, J. W. Violent expiratory events: on coughing and sneezing. Journal of Fluid Mechanics 745, 537–563 (2014).
    
    

33. [33].Zhang, C. et al. The source control effect of personal protection equipment and physical barrier on shortrange airborne transmission. Building and Environment 211, 108751 (2022).
    
    

34. [34].Bhagat, R. K., Wykes, M. D., Dalziel, S. B. & Linden, P. Effects of ventilation on the indoor spread of covid-19. Journal of Fluid Mechanics 903, F1 (2020).
    
    

35. [35].Franco, A. & Leccese, F. Measurement of CO2 concentration for occupancy estimation in educational buildings with energy efficiency purposes. Journal of Building Engineering 32, 101714 (2020).
    
    

36. [36].Zuraimi, M. et al. Predicting occupancy counts using physical and statistical CO2-based modeling methodologies. Building and Environment 123, 517–528 (2017).
    
    

37. [37].Candanedo, L. M. & Feldheim, V. Accurate occupancy detection of an office room from light, temperature, humidity and CO2 measurements using statistical learning models. Energy and Buildings 112, 28–39 (2016).
    
    

38. [38].Szczurek, A., Maciejewska, M. & Pietrucha, T. Occupancy determination based on time series of CO2 concentration, temperature and relative humidity. Energy and Buildings 147, 142–154 (2017).
    
    

39. [39].Caĺi, D., Matthes, P., Huchtemann, K., Streblow, R. & Müller, D. CO2 based occupancy detection algorithm: Experimental analysis and validation for office and residential buildings. Building and Environment 86, 39–49 (2015).
    
    

40. [40].Rahman, H. & Han, H. Bayesian estimation of occupancy distribution in a multi-room office building based on CO2 concentrations. In Building Simulation, vol. 11, 575–583 (Springer, 2018).
    
    

41. [41].Dong, B. et al. An information technology enabled sustainability test-bed (itest) for occupancy detection through an environmental sensing network. Energy and Buildings 42, 1038–1046 (2010).
    
    

42. [42].Wolf, S., Cali, D., Krogstie, J. & Madsen, H. Carbon dioxide-based occupancy estimation using stochastic differential equations. Applied Energy 236, 32–41 (2019).
    
    

43. [43].Jiang, C., Masood, M. K., Soh, Y. C. & Li, H. Indoor occupancy estimation from carbon dioxide concentration. Energy and Buildings 131, 132–141 (2016).
    
    

44. [44].Ryu, S. H. & Moon, H. J. Development of an occupancy prediction model using indoor environmental data based on machine learning techniques. Building and Environment 107, 1–9 (2016).
    
    

45. [45].Pedersen, T. H., Nielsen, K. U. & Petersen, S. Method for room occupancy detection based on trajectory of indoor climate sensor data. Building and Environment 115, 147–156 (2017).
    
    

46. [46].Wang, S. & Jin, X. CO2-based occupancy detection for on-line outdoor air flow control. Indoor and Built Environment 7, 165–181 (1998).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1420326X9800700305&link_type=DOI) 

47. [47].Wang, S., Burnett, J. & Chong, H. Experimental validation of CO2-based occupancy detection for demandcontrolled ventilation. Indoor and Built Environment 8, 377–391 (1999).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1420326X9900800605&link_type=DOI) 

48. [48].Kusuda, T. Control of ventilation to conserve energy while maintaining acceptable indoor air quality (National Bureau of Standards, Department of Commerce, 1976).
    
    

49. [49].ASHRAE. Ventilation for acceptable indoor air quality. ANSI/ASHRAE/IES Standard 62.1-2019 (2019).
    
    

50. [50].Dunn, D. R. Monthly Energy Review, April 2022, Section 2. Tech. Rep. DOE/EIA-0035(2022/4), US Energy Information Administration (2022).
    
    

51. [51].Abboushi, B. et al. Energy implications of using upper room germicidal ultraviolet radiation and hvac strategies to combat SARS-CoV-2. Tech. Rep., Oak Ridge National Lab.(ORNL), Oak Ridge, TN (United States) (2022).
    
    

52. [52].Escombe, A. R. et al. Upper-room ultraviolet light and negative air ionization to prevent tuberculosis transmission. PLoS medicine 6, e1000043 (2009).
    
    

53. [53].First, M., Rudnick, S. N., Banahan, K. F., Vincent, R. L. & Brickner, P. W. Fundamental factors affecting upperroom ultraviolet germicidal irradiation—part i. experimental. Journal of Occupational and Environmental Hygiene 4, 321–331 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/15459620701271693&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17365506&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F20%2F2023.03.19.23287460.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000246338600007&link_type=ISI) 

54. [54].Rudnick, S. N. & First, M. W. Fundamental factors affecting upper-room ultraviolet germicidal irradiation—part ii. predicting effectiveness. Journal of Occupational and Environmental Hygiene 4, 352–362 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1080/15459620701298167&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17454503&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F20%2F2023.03.19.23287460.atom) 

55. [55].Buonanno, M., Welch, D., Shuryak, I. & Brenner, D. J. Far-UVC light (222 nm) efficiently and safely inactivates airborne human coronaviruses. Scientific Reports 10, 1–8 (2020).
    
    

56. [56].Eadie, E. et al. Far-UVC (222 nm) efficiently inactivates an airborne pathogen in a room-sized chamber. Scientific Reports 12, 1–9 (2022).
    
    

57. [57].Risbeck, M. J. et al. Quantifying the tradeoff between energy consumption and the risk of airborne disease transmission for building hvac systems. Science and Technology for the Built Environment 28, 240–254 (2022).
    
    

58. [58].Risbeck, M. J. et al. Modeling and multiobjective optimization of indoor airborne disease transmission risk and associated energy consumption for building hvac systems. Energy and Buildings 253, 111497 (2021).
    
    

59. [59].Leung, N. H. et al. Respiratory virus shedding in exhaled breath and efficacy of face masks. Nature medicine 26, 676–680 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41591-020-0843-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32371934&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2023%2F03%2F20%2F2023.03.19.23287460.atom)

 [1]: /embed/graphic-3.gif
 [2]: /embed/inline-graphic-1.gif
 [3]: /embed/graphic-6.gif
 [4]: /embed/inline-graphic-2.gif
 [5]: /embed/inline-graphic-3.gif
 [6]: /embed/inline-graphic-4.gif
 [7]: /embed/graphic-7.gif
 [8]: /embed/inline-graphic-5.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/inline-graphic-6.gif
 [15]: /embed/inline-graphic-7.gif
 [16]: /embed/inline-graphic-8.gif
 [17]: /embed/inline-graphic-9.gif
 [18]: /embed/graphic-13.gif
 [19]: /embed/inline-graphic-10.gif
 [20]: /embed/graphic-14.gif