key: cord-0144628-8cps9c6m authors: Robey, Alison; Fierce, Laura title: Sensitivity of airborne transmission of enveloped viruses to seasonal variation in indoor relative humidity date: 2021-06-28 journal: nan DOI: nan sha: 6d1212f506231782376b8a4b7a43a68af043e4d7 doc_id: 144628 cord_uid: 8cps9c6m In temperate climates, infection rates of enveloped viruses peak during the winter. While these seasonal trends are established in influenza and human coronaviruses, the mechanisms driving the variation remain poorly understood and thus difficult to extend to similar viruses like SARS-CoV-2. In this study, we use the Quadrature-based model of Respiratory Aerosol and Droplets (QuaRAD) to explore the sensitivity of airborne transmission to the seasonal variation in indoor relative humidity across the wide range of relevant conditions, using SARS-CoV-2 as an example. Relative humidity impacts the evaporation rate and equilibrium size of airborne particles, which in turn may impact particle removal rates and virion viability. Across a large ensemble of scenarios, we found that the dry indoor conditions typical of the winter season lead to slower inactivation than in the more humid summer season; in poorly ventilated spaces, this reduction in inactivation rates increases the concentration of active virions, but this effect was important when the susceptible person was farther than 2 m downwind of the infectious person. On the other hand, changes in particle settling velocity with relative humidity did not significantly affect the removal or travel distance of virus-laden scenarios. Late autumn to early spring is often referred to not as the winter season, but as the flu season, reflecting the seasonal variation in influenza. Infection rates of influenza surge in the winter [1] , similar to other airborne viruses sharing its enveloped structure and upper respiratory infection site [2, 3] . One potential cause of this trend is the seasonal variation of environmental conditions. While outdoor temperature and humidity levels vary dramatically throughout the year, humans spend around 90% of their time indoors [4, 5] where temperatures remain stable but relative humidity (RH) displays consistent summer highs and winter lows [6, 7] . Several mechanisms pertinent to airborne transmission of viruses are sensitive to this variability in RH, such as faster inactivation of enveloped virions at extreme RH levels [3, 6] and shorter travel distance of mid-sized particles at high RH levels [8] . However, it is unclear whether the impact of RH on either process is important enough to contribute to virus seasonality. In the context of familiar airborne diseases like the flu or the winter common colds (caused by human coronaviruses (HCoVs) 229E, HKU1, NL63, and OC43) [9, 10] , understanding this relationship between RH and the seasonality of infection rates is a retrospective question; hundreds of years of evidence have clearly established their seasonality [11] , and current research helps clarify why, not if, seasonal cycles have occurred. In the case of emerging viruses like SARS-CoV-2, however, the question is preemptive: will transmission vary seasonally in the future? The first year of an outbreak cannot answer this question [12, 13] , but establishing whether novel viruses are likely to follow seasonal trends is critical to preparing adequate mitigation strategies [11, 14, 15] . Determining the impacts of RH on airborne transmission is vital to progressing our understanding of the seasonal dynamics of any enveloped virus. However, quantifying those impacts requires a model holistic enough to account for the series of processes governing airborne transmission and efficient enough to span the tremendous variability in relevant conditions. Most models complex enough to provide a detailed representation of short-range and long-range airborne transmission [e.g. 16, 17] are too computationally expensive to represent the large uncertainty in indoor conditions, particle properties, and physiological parameters. To obtain a mechanistic understanding of the impact of RH on airborne transmission, while also accounting for the large uncertainties in model parameters, we performed thousands of simulations using the new Quadrature-based model of Respiratory Aerosol and Droplets (QuaRAD) [18] . We use the transmission of SARS-CoV-2 as an example due to its current relevance; however, conclusions are also relevant to other prominent winter viruses that share its enveloped virion structure [19] and ability to transmit through the airborne route [20] . In this paper, we first describe QuARAD (Section 2), emphasizing specific details about the interaction of RH with evaporation (Section 2.1) and viral inactivation (Section 2.2). We then examine the impacts of typical summer and winter RH levels on exposure in a single baseline scenario (Section 3.1) and across a large ensemble of scenarios (Section 3.2). Finally, we discuss the insights gained on the seasonality of airborne enveloped viruses and the utility of different mitigation strategies (Section 4). We assessed the impact of typical indoor RH levels on virion fate and transport using QuaRAD, a broad overview of which is given below. The full model description, including specific parameter distributions, can be found in [18] . The model approximates the size distribution of and viral load within expired particles using numerical quadrature. Such quadrature-based moment methods have been shown to accurately represent integrals over aerosol distributions with only a small number of particles [21, 22] . QuaRAD models the evolution of particles expired by an infectious individual and the subsequent exposure of a susceptible individual. The particle size distribution is represented with three superimposed lognormal modes, based on the measurements of speech emissions detailed in [23] . We estimate the viral loads associated with each particle using these size distribution measurements in combination with Gesundheit II and quantitative RT-PCR measurements of influenza virion emission on fine and coarse particles [24, 25] . We assume that influenza is a reasonable model for viral emissions of SARS-CoV-2 [26] . In QuaRAD, the size distribution of respiratory particles and its evolution is represented using only six weighted particles. For each quadrature point, we simulate the evaporation (Section 2.1) and the subsequent virion inactivation (Section 2.2). Particle dispersion within the expiratory jet of an infectious person is modeled using a Gaussian puff model [27] in combination with a buoyant turbulent jet model [28] , following the approach of [8] . We then sum over the quadrature points to approximate the integral of the size distribution inhaled by the susceptible individual. Input parameters are sampled from the realistic probability distributions detailed in [18] ; the baseline scenario discussed in Section 3.1 was simulated using the median of the probability density function from which the ensemble inputs were sampled. In this study, we analyzed the continuous stream of aerosol particles and droplets expelled when an infectious person speaks and the subsequent exposure of a susceptible individual, which we quantified as the number of virions inhaled per secondṄ inhale .Ṅ inhale was computed using the concentration of virions present at the mouth of the susceptible person, their average number of breaths per minute, and the average volume of air inhaled per breath. We quantify the difference in active virion inhalation between RH cases as described in Section 2.3. Upon expiration, aerosol particles and droplets evaporate to an equilibrium size that depends on their composition and the environmental conditions. Following the approach of [8] and [29] and implemented as detailed in [18] , we model evaporation by solving a coupled pair of ordinary differential equations representing the mass and heat transfer: Equation 1 describes the temporal evolution of the particle mass m p as a function of time t, ambient pressure p, particle diameter D p , molecular weight of water M w , diffusivity of water in air D ∞ , correction factor C T , Sherwood number Sh, gas constant R, surface vapor pressure p v,p , and ambient vapor pressure p v,∞ . Equation 2 describes the evolution of the particle temperature T p as a function of the particle's specific heat C p , the air thermal conductivity k g , the background temperature T ∞ , the Nusselt number Nu, the latent heat of vaporization L v , and Equation 1. These parameters are provided in [18] . The model represents convection-induced evaporation enhancements. We assume that particles are spherical and contain 1%-9% non-water aerosol components by volume, following the measurements of [30] . The vapor pressure over an aqueous droplet is computed from the κ-Köhler model [31] . Seasonal changes in RH modify p v,∞ , such that droplets evaporate more slowly with increasing RH, as shown in Figure 1 . RH is the ratio between p v,∞ and the saturation vapor pressure, which depends on temperature. Although QuaRAD is able to represent the variation in temperature and RH within the turbulent jet of an infectious person, we found that results were insensitive to this variation. Therefore, we assume that expired particles encounter the background temperature and humidity as soon as they are expired. Like other enveloped viruses [6] , SARS-CoV-2 virions exhibit exponential decay [32] with comparable rates on surfaces and while airborne [19, 33] . The inactivation of a virion eliminates its capability of infecting a cell and is analogous to the death of an organism. The trends of inactivation due to temperature and RH are established, though the exact inactivation rates remain poorly constrained; since our focus is on the indoor environment, we do not consider the accelerated inactivation rates caused by UV light [13] . As with influenza [6, 34] and other coronaviruses [35, 36] , SARS-CoV-2 decays more rapidly at high temperatures than low temperatures and shows a U-shaped dependence on RH, as seen in Figure 2 . At room temperature, for example, virion half-lives are longer than 6 hours for RH < 45% or RH > 75% but are around 2.5 hours for RH ≈ 65% [32] . Several mechanisms may cause the temperature dependence of viral inactivation, including thermal denaturation [6, 32, 37] and decreased lipid ordering in viral envelopes at high temperatures [6, 38] . The RH dependence of viral inactivation is likely caused by RH's influence over solute concentrations in infectious particles. Solutes are thought to function as reactants in virion inactivation, and the U-shaped humidity dependence is thus theorized to result from two ways in which the interior solute concentrations are reduced [32, 39] . At low RH levels, evaporation causes solutes in expired particles to effloresce, separating from the aqueous solution as solids on the particle surface. If the RH within the indoor space is above a threshold, known as the efflorescence relative humidity (ERH), solutes within the aqueous particle will remain in the solution. In this regime, the solute concentration decreases with increasing RH because, as long as the RH remains above the ERH, the concentration of solutes will remain more dilute the more water is removed through evaporation. For this reason, inactivation rates are fastest just above the ERH, when solutes are in solution and are also highly concentrated. [32] . We focused on the half-lives of 6.43 h at RH = 40% and 2.41 h at RH = 65% measured at room temperature (shown in purple), corresponding to inactivation rates of 0.108 h −1 and 0.288 h −1 , respectively. In this model, we used the measured distributions of SARS-CoV-2 halflives h in a cell culture medium at three temperatures (283.15 K, 295.15 K, 300.15 K) and RH levels (40%, 65%, 85%) from [32] to estimate the distribution of viral decay rates k given in Figure 2 . Though we recognize that ERH varies with particle size and composition, we followed [32] in assuming an ERH of 45% across particles. These values were chosen due to the explicit variation in environmental conditions but are comparable to those obtained in [19] , [33] , and [40] . Given the focus of this paper on seasonal indoor environments, we analyzed the measurements taken at T = 295.15 K, corresponding to an average indoor temperature of 22 • C, and at RH = 40% or RH = 65%. We quantify the effect of RH on virion exposure as the relative difference between virion inhalation rates,Ṅ inhale , at two RH levels: Positive values of the relative difference ∆ r indicate an increase in the inhalation rate of active virions with a decrease in RH from 65% to 40%. These RH levels match the conditions under which [32] measured virion half-lives and approximately correspond to the median indoor RH of 62% in the summer and 42% in the winter measured by [41] . Susceptible individuals are assumed to inhale at a rate of 14 breaths per minute and a volume of 4.69 × 10 −4 m 3 per breath in the base case and at 12 to 20 breaths per minute and 3.75×10 −4 to 6.25 × 10 −4 m 3 per breath across the ensemble [18] . In this section, we present the sensitivity of virion inhalation rate to RH as a function of downwind distance from an infectious individual. Before describing the effect of RH on airborne transmission across an ensemble of simulations in Section 3.2, we first show how RH impacts key processes in an example scenario in Section 3.1. In the baseline scenario, we found that the sensitivity ofṄ inhale to RH depends on the distance between the infectious and the susceptible person ( Figure 3 ). When directly downwind of an infectious person,Ṅ inhale is insensitive to RH, whereas we find a 10% increase inṄ inhale with a decrease in RH from 65% to 40% (solid line in Figure 3B ). When inactivation is neglected, ∆ r ≈ 0 both near and far from the infectious person (dotted line in Figure 3B ). If a susceptible person is within the expiratory jet of an infectious person, nearly all virions they inhale would have been expelled from an infectious person just seconds earlier. On the other hand, far from an infectious person, a large portion of virions may be contained in particles that were expelled hours earlier. For this reason,Ṅ inhale far from an infectious person depends strongly on the RH-dependent inactivation rates (∆ r ≈ 0.1), whereas nearfield exposure toṄ inhale is insensitive to changes in inactivation rates (∆ r ≈ 0). To quantify the impact of RH across the wide range of conditions expected in indoor spaces, we performed an ensemble of 1000 scenarios using the distributions in input parameters described in [18] . As in the base case, N inhale was insensitive to RH when directly downwind of an infectious person. Beyond distances of approximately 2 m,Ṅ inhale was higher at a typical winter RH of 40% than at a typical summer RH of 65%. The mean ∆ r (Figure 4 ) follows the same trend as the base case. We again found that the sensitivity ofṄ inhale to RH was driven by differences in virion inactivation and not by differences in particle settling. Since RH significantly affects transmission only when the susceptible person is far from the infectious person, we explored the sensitivity ofṄ inhale to conditions controlling far-field virion concentrations: the duration of the encounter and the ventilation rate of the room (expressed here in air changes per hour, or ACH). If the infectious person and the susceptible person are in a poorly ventilated space, the virion concentrations throughout the room increase over time. Under these conditions, virion exposure is highly sensitive to inactivation rate, as shown by the high ∆ r values at low ventilation rates in Figure 5 . Virion concentrations also increase as the duration of the encounter increases. Consequently, longer interactions correspond with higher exposure and a higher ∆ r (comparison between ∆ r in a 1-hour interaction (green) and in a 4-hour interaction (purple) in Figure 5 ). We also found that high ventilation rates reduce the dependence of ∆ r on the length of the encounter. For example, when the ventilation rate is less than 1.0 ACH, the median ∆ r after four hours is nearly twice the median ∆ r after one hour, but when the ventilation rate is greater than 3.0 ACH, the median ∆ r values after four hours and after one hour are approximately equal. Our process analysis revealed that the greater virion exposure at RH = 40% than at RH = 65% was caused by differences in virion inactivation rate. In indoor spaces, SARS-CoV-2 virions have long half-lives (h = 6.4 ± 0.03 hours at RH = 40% and h = 2.4 ± 0.04 hours at RH = 65%), so nearly all virions in freshly expired particles are active, regardless of the background RH. Whereas virion exposure directly downwind of an infectious person is driven by virions in freshly expired particles, virion exposure far from an infectious person may be strongly influenced by particles that remain suspended for hours. Whereas changes in inactivation rates with RH led to seasonal differences in virion survival within respiratory particles, removal rates and travel distances of the respiratory particles themselves were insensitive to seasonal changes in RH. RH does affect the evaporation timescale and equilibrium size of expired particles (see Figure 1 ), but the impact of those differences on travel distance is limited by the expired size distribution. Close to the infectious individual, most expired particles are small enough to travel with the flow of the respiratory jet at their initial diameter, during evaporation, and at their equilibrium diameter; most of the other expired particles are large enough to quickly settle out of the respiratory jet at their initial diameter, during evaporation, and at their equilibrium diameter. Far from the infectious individual, the equilibrium sizes of most particles at RH = 40% and RH = 65% are too similar to cause major differences in particle travel or residence time. The plume does contain mid-sized particles whose spread and settling velocity are sensitive to RH, but these particles make up only a tiny proportion of those expired during speech and consequently have very little effect onṄ inhale . While this finding does not contradict earlier findings of RH sensitivities in the spread of specific particles, it does show that those differences have very little impact on transmission risk when integrating over a population of polydisperse particles. Analysis across the ensemble revealed that at 4 m of distance, RH has a much greater impact on virion exposure over the course of a long encounter than over the course of a short one. Far-field concentrations increase over time, such that inhaling air in the far-field after only a few minutes of emissions results in a much lowerṄ inhale than after an hour of emissions; near-field concentrations do not show this time dependence because, in close proximity,Ṅ inhale is controlled by the freshly expired virions in the expiratory jet instead of the virions building up over time throughout the room. Because longer timescales result in greater virion concentrations than shorter timescales and the decay rate of far-field virion concentrations depend on RH, longer timescales thus allow the decay rate to have a greater impact oṅ N inhale , resulting in an elevated ∆ r . Far-field concentrations are also affected by the ventilation rate of the room. Ventilation removes particles from the air, so when the ACH is low, far-field particle concentrations build up over time. When particles linger for a long time, the difference in inactivation rate between RH scenarios results in a high ∆ r . On the other hand, when the ACH is high, far-field particle concentrations are reduced more quickly. The impact of RH on exposure thus decreases as the ACH increases, resulting in a small ∆ r in well-ventilated spaces. We thus find that the impact of RH on far-field exposure over long timescales can be meaningfully reduced by increasing the ventilation rate. While we found that typical winter and summer indoor RH levels indeed impact virion exposure, it is very likely that other factors not accounted for in this model are also in play. For example, immune response capabilities of susceptible individuals may be impaired by seasonal deficiencies in vitamin D [42] or RH-mediated reductions in mucociliary clearance [43] , each of which may increase the probability of winter infections. Frequency of exposure may also be determined by seasonal differences in human behavior, including the school calendar and the higher frequency of indoor social events in colder seasons [44, 45] . We also note that QuaRAD accounts solely for seasonal differences in airborne transmission; while this is likely the dominant transmission route of SARS-CoV-2 [46] and a prominent one for other discussed viruses, differences in droplet or fomite transmission may independently impact seasonality. In this paper, we showed that a low RH increases airborne virion exposure by slowing the virion inactivation rate, not by increasing the travel distance or residence time of respiratory particles. In poorly ventilated spaces and encounters lasting four hours, we found that lowering the RH from 65% to 40% led to a 20% decrease in the median exposure to active virions when the susceptible individual is at least 2 m from the infectious individual. Close-range exposure was insensitive to this seasonal variation in indoor RH. Whereas previous studies [e.g. 8] have shown that changes in RH can significantly impact the dispersion of particles in specific size ranges, we found that differences in evaporation speed and equilibrium size have a trivial impact on the overall dispersion due to the tiny fraction of the plume those particles make up. The timescale at which virions become inactive, on the other hand, is highly sensitive to RH and thus strongly impacts the inhalation rate of active virions. This finding highlights the importance of measurements that constrain the relationship between virion inactivation and environmental properties [e.g. 32]. Though we focused our analysis on SARS-CoV-2, we expect airborne exposure to other enveloped viruses to follow the same trends. We note that the sensitivity of virion inhalation rate to RH will be more pronounced for viruses with shorter half-lives and, furthermore, that while we exclusively studied the rate of virions inhaled, it is the absolute number of virions inhaled that determines the probability of initial infection. As such, small sensitivities in inhalation rate to RH may, over time, result in larger sensitivities in the actual probability of infection to RH. The observed sensitivity of virion exposure to RH is reduced with higher ventilation and shorter exposure times, suggesting high ventilation and shorter interactions are critical to limiting winter transmission risk even when following safe social distancing guidelines. In settings where short exposure times are impossible or where multiple groups of people transition between the same room, such as classrooms, hospitals, factories, or prisons, employing other protections like effective mask wearing and high ventilation are critical to limiting transmission risks. The QuaRAD source code, input files, and processing script are available for download at: https://github.com/lfierce2/QuaRAD/. Simulation ensembles were created using latin hypercube sampling with pyDOE: https://pythonhosted.org/pyDOE/. Laboratory, a multiprogram national laboratory supported by DOE Contract DE-SC0012704. The authors would like to thank Robert McGraw for helpful discussion and his insights on particle chemistry. Global influenza seasonality: reconciling patterns across temperate and tropical regions Environmental predictors of seasonal influenza epidemics across temperate and tropical climates Dynamics of airborne influenza a viruses indoors and dependence on humidity The national human activity pattern survey (nhaps): a resource for assessing exposure to environmental pollutants Indoor time-microenvironment-activity patterns in seven regions of europe Mechanistic insights into the effect of humidity on airborne influenza virus survival, transmission and incidence Socioeconomic and outdoor meteorological determinants of indoor temperature and humidity in new york city dwellings Enhanced spread of expiratory droplets by turbulence in a cough jet Global seasonality of human coronaviruses: a systematic review Seasonality of coronavirus 229e, hku1, nl63 and oc43 from Seasonality of infectious diseases Grenfell, Susceptible supply limits the role of climate in the early sars-cov-2 pandemic Misconceptions about weather and seasonality must not misguide covid-19 response Forecasting disease risk for increased epidemic preparedness in public health Modelling the trajectory of disease outbreaks works Using large eddy simulation to study particle motions in a room Large-eddy simulation of human-induced contaminant transport in room compartments Simulating near-field enhancement in transmission of airborne viruses with a quadrature-based model Increasing temperature and relative humidity accelerates inactivation of sars-cov-2 on surfaces Airborne transmission of severe acute respiratory syndrome coronavirus 2 (sars-cov-2): What we know Description of aerosol dynamics by the quadrature method of moments Multivariate quadrature for representing cloud condensation nuclei activity of aerosol populations Modality of human expired aerosol size distributions Influenza virus aerosols in human exhaled breath: particle size, culturability, and effect of surgical masks Respiratory virus shedding in exhaled breath and efficacy of face masks Viral load of sars-cov-2 across patients and compared to other respiratory viruses, Microbes and infection Modeling indoor air exposure from short-term point source releases Turbulent jets and plumes: A Lagrangian approach The interdependence of evaporation and settling for airborne freely falling droplets Physico-chemical characteristics of evaporating respiratory fluid droplets A single parameter representation of hygroscopic growth and cloud condensation nucleus activity The effect of temperature and humidity on the stability of sars-cov-2 and other enveloped viruses, bioRxiv Aerosol and surface stability of sars-cov-2 as compared with sars-cov-1 Relationship between humidity and influenza a viability in droplets and implications for influenza's seasonality Effects of air temperature and relative humidity on coronavirus survival on surfaces Influence of relative humidity on the survival of some airborne viruses Thermal inactivation of animal viruses Progressive ordering with decreasing temperature of the phospholipids of influenza virus Humidity-dependent decay of viruses, but not bacteria Stability of sars-cov-2 on surfaces, bioRxiv The relationship between indoor and outdoor temperature, apparent temperature, relative humidity, and absolute humidity Epidemic influenza and vitamin d Criteria for enhancing mucus transport: a systematic scoping review Effects of school closure on incidence of pandemic influenza in alberta, canada Seasonal variation in airborne infection risk in schools due to changes in ventilation inferred from monitored carbon dioxide Ten scientific reasons in support of airborne transmission of sars-cov-2