Dynamics of Airborne Influenza A Viruses Indoors and Dependence on Humidity

doi:10.1371/journal.pone.0021481

PLoS One. 2011; 6(6): e21481.

Published online 2011 June 24. doi: 10.1371/journal.pone.0021481

PMCID: PMC3123350

Copyright Yang and Marr. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Dynamics of Airborne Influenza A Viruses Indoors and Dependence on Humidity

Wan Yang and Linsey C. Marr^*

Department of Civil and Environmental Engineering, Virginia Tech, Blacksburg, Virginia, United States of America

Ron A. M. Fouchier, Editor

Erasmus Medical Center, The Netherlands

* E-mail: lmarr/at/vt.edu

Conceived and designed the experiments: WY LCM. Performed the experiments: WY. Analyzed the data: WY LCM. Wrote the paper: WY LCM.

Received March 8, 2011; Accepted May 30, 2011.

This article has been cited by other articles in PMC.

Abstract

There is mounting evidence that the aerosol transmission route plays a significant role in the spread of influenza in temperate regions and that the efficiency of this route depends on humidity. Nevertheless, the precise mechanisms by which humidity might influence transmissibility via the aerosol route have not been elucidated. We hypothesize that airborne concentrations of infectious influenza A viruses (IAVs) vary with humidity through its influence on virus inactivation rate and respiratory droplet size. To gain insight into the mechanisms by which humidity might influence aerosol transmission, we modeled the size distribution and dynamics of IAVs emitted from a cough in typical residential and public settings over a relative humidity (RH) range of 10–90%. The model incorporates the size transformation of virus-containing droplets due to evaporation and then removal by gravitational settling, ventilation, and virus inactivation. The predicted concentration of infectious IAVs in air is 2.4 times higher at 10% RH than at 90% RH after 10 min in a residential setting, and this ratio grows over time. Settling is important for removal of large droplets containing large amounts of IAVs, while ventilation and inactivation are relatively more important for removal of IAVs associated with droplets <5 µm. The inactivation rate increases linearly with RH; at the highest RH, inactivation can remove up to 28% of IAVs in 10 min. Humidity is an important variable in aerosol transmission of IAVs because it both induces droplet size transformation and affects IAV inactivation rates. Our model advances a mechanistic understanding of the aerosol transmission route, and results complement recent studies on the relationship between humidity and influenza's seasonality. Maintaining a high indoor RH and ventilation rate may help reduce chances of IAV infection.

Introduction

Influenza A has a clear seasonal pattern in temperate regions, yet the underlying cause for it remains controversial despite nearly a century of investigation. The literature identifies numerous factors that may influence influenza's seasonality: environmental conditions such as temperature, humidity, and ultraviolet radiation; immune function; school schedules; and human mobility patterns and contact rates [1]. Among these, the leading contenders are humidity and temperature [2], [3], [4], and in indoor environments, where people spend ~90% of their time, humidity is the more variable factor. Particularly in the developed world where heating, ventilation, and air conditioning systems (HVAC) are the norm, indoor temperature tends to fall in a narrower range, and thus its influence is limited. Recent studies using a guinea pig experimental model [3], [4] indicate that low relative humidity (RH) favors aerosol transmission of influenza A viruses (IAVs), in which they are transmitted by small respiratory droplets expelled from infected hosts. Nevertheless, the precise mechanisms by which humidity might influence influenza viability and transmissibility via the aerosol route have not been elucidated.

Humidity may affect airborne IAV transmission via two important variables. The first is droplet size. When released from the respiratory tract (assumed to have 100% RH), droplets experience rapid evaporation and shrinkage upon encountering the unsaturated ambient atmosphere. The ultimate size of a droplet depends on ambient humidity, and size determines aerodynamic behavior and whether the droplet will settle to the ground quickly or remain suspended in the air long enough to possibly cause a secondary infection. Previous studies on evaporation of respiratory droplets usually used water or simple saline solutions (e.g., NaCl) to simulate respiratory fluid [5], [6]. However, respiratory fluid is a complicated combination of water, salts, and various organic compounds [7], [8] that affect the thermodynamics of evaporation, compared to pure water or saline solutions. The equilibrium droplet size is affected by surface curvature and solute effects, the combination of which is described by Köhler theory [9]. While the vapor pressure is enhanced over curved versus flat surfaces, it is reduced by the presence of solutes. These competing effects are magnified at smaller droplet diameters and determine the equilibrium size at a particular RH.

The second variable that is sensitive to humidity is IAV viability [2], [10], [11], [12], [13]. Hemmes et al. [2] linked influenza's seasonality to the seasonal oscillation of RH indoors, based on their experiment on death-rate variation versus RH. Shaman and Kohn [14], on the other hand, concluded that absolute humidity (AH) rather than RH modulates influenza seasonality through constraint on viability, but whether AH or RH is controlling is under debate. For our purposes, differentiating between the two is not possible because this work focuses on a narrow range of typical indoor temperatures. Finally, it is possible that the two variables—final droplet size and viability—are linked, if evaporation and subsequent concentration of solutes in respiratory droplets affects IAV viability in aerosols.

Elucidating the causes of influenza's seasonality will require improved comprehension of transmission mechanisms, especially the aerosol route. To advance a mechanistic understanding of the role of humidity in aerosol transmission, we model the change in size of respiratory droplets and IAV inactivation at RHs ranging from 10% to 90%. Based on these results, we further model the dynamics of droplets emitted from a cough in an indoor environment and illustrate the evolution of infectious IAV concentrations and size distributions, considering removal by gravitational settling, ventilation, and viral inactivation. We are thus able to determine the magnitude by which humidity affects airborne concentrations of infectious IAVs.

Results

Initial size distribution of droplets expelled from a cough

We located seven papers that reported the size distribution of droplets expelled while coughing, sneezing, and speaking [6], [15], [16], [17], [18], [19], [20]. The droplet diameter geometric mean (GM), size range, and droplet number varied greatly among the different studies, as summarized in Table 1. Because ~80% of patients with influenza manifest symptoms of coughing [21], we focus on droplets emitted from coughing to demonstrate the dynamics of airborne IAVs. Droplet diameters as small as 0.3 µm and as large as 2000 µm have been observed from coughing, and the GMs in the studies ranged from 0.25 µm to 96.6 µm, depending on the experimental methods used (Table 1). Four of the studies reported GMs in a much narrower range of 8.4 µm to 16.0 µm. We adopted Duguid's [20] results on the basis of reliability of the methods and care and thoroughness of the experimental design and analysis. Though the study dates to 1946, its results were similar to those of more contemporary work that use modern aerosol characterization equipment [6], [17]. Similar to Nicas et al. [22], we found that droplets between 1–100 µm from coughing follow a log-normal size distribution. Eq. 1 describes the relationship:

A mathematical equation, expression, or formula.
Object name is pone.0021481.e001.jpg

(1)

where D_i is the initial droplet diameter in µm, and z is the corresponding quantile of a normal distribution with the same cumulative probability. According to this relationship, the size distribution of droplets emitted from coughing has a GM of 12.9 µm and a geometric standard deviation (GSD) of 2.3 (Figure 1).

Table 1

Prior studies of respiratory droplet size distributions.

Figure 1

Log-probability plot of droplet size distribution from a cough, adapted from Duguid [20].

Respiratory droplet size transformation

Equilibrium droplet size is attained nearly instantaneously upon release. A 20-µm droplet shrinks to one-half of its original diameter in less than half a second at 20°C [22]. Table 2 shows the equilibrium, or final, diameters (D_eq) of droplets with D_i of 0.1, 1, and 10 µm at 10–90% RH. These were calculated based on a model of droplet transformation that assumes separate solutes and volume additivity (SS-VA) [23]. Due to the Kelvin effect, evaporation of smaller droplets is enhanced, and the equilibrium diameters are smaller. The ratio D_eq/D_i is 0.490 at 90% RH for a respiratory droplet with D_i

0.1 µm, versus 0.516 under the same conditions for a larger one with D_i

10 µm. However, the Kelvin effect is negligible for droplets with D_i>0.1 µm, and D_eq/D_i is independent of RH for droplets larger than 1 µm.

Table 2

Respiratory droplet size transformation.

Efflorescence or crystallization of NaCl, a major component of respiratory droplets, due to loss of water is expected to occur between 40–50% RH [24]. According to our results, between 10–40% RH, D_eq/D_i varies by only 3.7% (0.402–0.417); in comparison, between 50–90% RH, the ratio varies by 21.7% (0.424–0.516). Respiratory droplets lose almost all their water at low RHs. For comparison, we also calculated the D_eq/D_i ratios based on volume additivity using experimental data on dehydration of droplets containing NaCl [25] and on hydration of those containing a glycoprotein [26]. Differences between modeled and experimental results are less than 4% (Table 2).

Inactivation of airborne IAVs

The viability of IAVs decreases over time and is affected by environmental variables such as temperature, humidity, and UV radiation [12], [27]. The inactivation rate (k), derived from experimental data on airborne IAVs [12], is linearly correlated with RH (Figure 2), following the relationship

A mathematical equation, expression, or formula.
Object name is pone.0021481.e002.jpg

(2)

with an r² of 0.977 and p-values for the model, intercept, and slope of 0.0015, 0.059, and 0.0015, respectively.

Figure 2

IAV inactivation rate versus RH.

Evolution of infectious IAV distribution after a cough

The well-mixed indoor air model we developed for infectious IAVs accounts for removal by gravitational settling, ventilation, and viral inactivation. The concentration of infectious IAVs associated with droplets of a specific diameter D_eq in a room at time t is

A mathematical equation, expression, or formula.
Object name is pone.0021481.e003.jpg

(3)

where C₀ is the initial concentration of infectious IAVs associated with droplets of size D_eq in the room, v is the settling velocity, H is the height of the room, λ is the air exchange rate (AER) and assumes no recirculation, and k is the inactivation rate. The inactivation rate k depends on RH, according to Eq. 2, and the settling velocity v depends on D_eq, which also depends on RH. Eq. 3 can be integrated over all droplet sizes to obtain the total concentration of infectious IAVs in a room. We calculate results for AERs of 1 air change per hour (ACH) and 10 ACH, typical of residential and public settings, respectively [28], [29]. For simplicity, we assumed that room heights are the same in residences and public settings such as offices, classrooms, and hospitals, where people aggregate and thus have a higher risk of infection. The equation is applicable for a single rectangular room, does not depend on the volume of the room, and does not account for air exchange between multiple rooms.

Figure 3 shows the evolution of IAV concentrations in time in terms of both the total number of infectious viruses (Figure 3A and 3B) and size distribution (Figure 3C and 3D). Only emitted droplets with D_i≤100 µm are considered, as larger ones will be removed by gravitational settling within seconds. Following a single cough in a well-mixed room, the concentration of infectious IAVs is initially 1.8×10³ # m⁻³ under the assumptions of this study. Figures 3A and 3B show that the total number of infectious IAVs falls rapidly with time and that the loss is greater at higher RH and in public versus residential settings. If one infected person is continuously shedding viruses by coughing 15 times per hour [30], then the concentration of IAVs will be ~2×10³ # m⁻³ in a public setting. This concentration is similar in magnitude to those measured in hospitals, medical clinics, day care facilities, and airplanes [31], [32], [33]. Under conditions of higher RH, removal by settling is more effective because droplets shrink less, and inactivation is more rapid. Removal of 99.9% of the IAVs emitted requires much greater time in a residential versus public setting, indicating that ventilation is an important removal mechanism and that airborne IAVs can persist for longer times in settings with lower AERs.

Figure 3

Evolution of infectious airborne IAV concentrations and size distributions.

Figure 3C and 3D show that, at 50% RH, emitted droplets shrink to about half of their original diameters due to evaporation. Evaporation happens almost instantaneously [22], so while the initial size distribution at 0 min extends out to 100 µm, all subsequent ones end at 42 µm (D_eq/D_i

0.42). This process greatly increases the fraction of IAVs that are associated with smaller droplets, since the virus concentration within a droplet increases by a factor of 8 (2³) as its diameter shrinks by half. For instance, compared to the initial droplets emitted at 0 min, which quickly reach their equilibrium diameters, the number of IAVs associated with equilibrium droplet sizes smaller than 25 µm increases by a factor of 5.2 at 50% RH. Due to more rapid settling, IAVs associated with larger droplets are lost faster than are those associated with smaller ones. Consequently, IAVs associated with smaller droplets become more dominant over time, as indicated by the shifting of the peak of ΔC/ΔD_eq to the left in Figure 3C and 3D. The diameter of droplets containing the most IAVs (i.e., the mode of the distribution) shifts from ~50 µm upon release to ~16 µm at 1 min, ~10 µm at 10 min, and ~5 µm at 60 min in residential settings; a similar trend is shown in Figure 3D for public settings.

Humidity dependency and removal mechanisms

Figure 4 shows the effect of RH on size distributions of infectious IAVs, 10 min after a cough in residential and public settings. For both cases, IAV concentrations decrease with increasing RH across all sizes, but the modes of the distributions remain around 9–10 µm. The total IAV concentration (i.e., the area under each curve) decreases with increased RH. As a result, it takes twice as long to remove 99.9% of IAVs emitted at 10% RH than that at 90% RH in residential settings (>100 min at 10% RH versus <50 min at 90% RH, as shown in Figure 3A).

Figure 4

IAV size distributions.

Humidity affects both settling, because of its dependence on size transformation, and inactivation of IAVs. The relative importance of these two effects can be illustrated by comparing the ratios of virus concentrations at 10% RH versus those at 90% RH at varying times. The ratios increase approximately exponentially with time: 2.4 at 10 min, 5.4 at 30 min, and 16.1 at 60 min. If only inactivation were considered, these factors would instead be 1.4, 2.7, and 7.3, respectively; and if only settling were considered, the corresponding factors would be 1.7, 2.0, and 2.2, respectively. These ratios are independent of the ventilation rate. The much narrower range of factors for settling than for inactivation (i.e., 1.7–2.2 versus 1.4–7.3) indicates that RH has a greater impact on inactivation, especially over long periods (>30 min).

Figure 5 shows the effectiveness of each removal mechanism—settling, ventilation, and inactivation—independently as a function of RH (Figure 5A and 5B) and droplet size (Figure 5C and 5D), 10 min following a cough. Figure 5A shows that gravitational settling is the dominant removal mechanism in residential settings. Settling alone removes over 80% of airborne IAVs within 10 min, and its removal efficiency increases slightly with RH, from 87% to 92% across the range of RHs. In contrast, ventilation only removes 15% of total IAVs, regardless of RH. Removal efficiency by inactivation increases with RH, accounting for up to 28% at the highest RH. Figure 5B shows that ventilation and gravitational settling are both important in removing airborne IAVs from public settings with higher AERs. At an AER of 10 ACH with no recirculation, ventilation removes 81% of airborne IAVs. Settling and inactivation are independent of ventilation rate and remove the same amounts of IAVs as in residential settings.

Figure 5

IAV removal mechanisms.

Removal efficiencies for IAVs vary as a function of droplet size for settling but not ventilation or inactivation. Figure 5C and 5D show removal efficiency versus equilibrium droplet diameter at 50% RH, 10 min following a cough, for each mechanism individually and all three together. Because settling velocity scales with diameter squared, removal efficiencies due to gravitational settling range from only 0.7% for droplets with D_eq

1 µm to 51.2% for those with D_eq

10 µm to >98.8% for those with D_eq>25 µm. Ventilation is equally effective for all sizes, with removal efficiencies of 15% in residential settings (Figure 5C) and 81% in public settings (Figure 5D), depending on the AER.

Overall, gravitational settling is the main removal mechanism in both residential and public settings (Figure 5A and 5B). It removes a disproportionately large fraction of IAVs because it favors larger droplets, which contain far more IAVs, as their numbers are proportional to the initial droplet volume, or D_i³. However, settling is ineffective at removing droplets <5 µm, as shown in Figure 5C and 5D. Ventilation is important in public settings and particularly so for removal of smaller droplets (<5 µm) for which settling is inefficient. It accounts for ~50% of total removal of IAVs associated with droplets <5 µm in residential settings (Figure 5C) and ~80% in public settings (Figure 5D). Inactivation increases with RH and is maximal at 28% at 90% RH, 10 min following a cough (Figure 5A and 5B). Although the removal efficiency by inactivation is relatively low, it is important when removal by ventilation and settling are both minor. For instance, in residential settings (Figure 5C), inactivation accounts for ~50% of total removal of IAVs associated with droplets <5 µm.

Discussion

Removal of infectious IAVs

Higher RH favors removal of infectious IAVs. Since larger droplets have greater settling velocities, higher RHs, at which D_eq/D_i is larger, thereby will accelerate the removal rate. Additionally, the inactivation rate of IAVs increases with increasing RH (Eq. 2). According to our model, the concentration of airborne IAVs resulting from a cough would be reduced by 10% if the RH increases from 35%, the mean indoor RH in heating season [30], to 50%, 10 min following the cough, and by 40% after 1 h in residential settings. These estimates agree in magnitude with those reported by Myatt et al. [30], whose model suggests that influenza virus survival decreases by 17.5–31.6% when indoor RH increases by 11–19% over 15 h. Hence, maintaining a reasonably high indoor RH (e.g., 50%) may accelerate the removal of infectious IAVs and help prevent or reduce influenza infection.

The relative importance of the two mechanisms—droplet size transformation and inactivation—as a function of humidity is of interest. Shaman and Kohn [14] concluded that AH modulates influenza transmission by influencing the virus' survival rate, rather than by enhancing production of airborne droplet nuclei in low humidity conditions. We found that respiratory droplets would shrink to one-half of their original diameters at 90% RH, and to around two-fifths at 10% RH. It thus appears that changes in droplet size are dramatic at unsaturated RHs and that variations due to differences in RH are relatively trivial. Our analysis shows that removal by inactivation is more variable with RH than is removal by settling. This may explain why Shaman and Kohn [14] could find a statistically significant relationship between AH and influenza survival but not transmission. However, this does not suggest that droplet shrinkage in response to unsaturated RHs is not important for influenza transmission, only that it is not as obvious as the induced change in viability.

We have demonstrated the relative importance of the three removal mechanisms. Settling can remove over 80% of droplets emitted from a cough within 10 min; however, it is effective only for larger droplets and allows the smaller ones (<5 µm) to remain suspended. In contrast, ventilation is able to remove all droplets regardless of size simply by air exchange. Therefore, higher AERs will facilitate the elimination of virus-containing droplets from indoor environments, especially to compensate for the inefficacy of settling in removing the small ones. This observation also justifies the requirement to maintain a high AER in public places (e.g., 12 ACH in hospital waiting areas [34]). Removal efficiencies due to virus inactivation are relatively small (i.e., 0–28% in 10 min, if only inactivation were considered). However these estimates are based on experimental data reported by Harper [12], which indicated lower inactivation rates of 0.0031–0.028 min⁻¹ at 20–81% RH, compared to 0.0073±0.0031 min⁻¹ at 15–40% RH and 0.091±0.024 min⁻¹ at 50–90% RH, as reported by Hemmes [2]. If estimated Hemmes' data [2], the corresponding removal efficiencies would be larger: 7.0% at 15–40% RH and 59.8% at 50–90% RH. Virus inactivation may thereby play a more significant role depending on the actual inactivation rate.

IAV viability, seasonality, and humidity dependency

Experimental and/or theoretical models have been constructed to predict the viability of airborne IAVs as a function of humidity [14], [35], but a widely accepted mechanistic explanation for the relationship is still lacking. Studies on the effect of humidity agree that IAVs survive better at lower RHs. However, Hemmes [2], [11] and Harper [12] found higher inactivation rates at both medium and high RHs, in contrast to Shechmeister [10] and Schaffer et al. [13], who found higher inactivation rates at medium but not high RHs. This disparity may stem from the different compositions of media used in each experiment. All media contained salts (approximately 0.5–3%); however, those used in the former two experiments contained far more proteins than did those in the latter two. High concentrations of salts are found to be detrimental to avian IAVs [36]. As water in the droplets evaporates, solute concentrations increase and may consequently accelerate IAV inactivation. However, an NaCl droplet can suddenly lose all of its water and crystallize at the point of efflorescence (45–48% RH) [25], thus eliminating the negative effect of dissolved salts at low RHs. This effect is perceivable in Table 2, which shows that D_eq/D_i varies little when RH≤40%. The combination of increasing salt concentrations followed by efflorescence as RH decreases may explain the trend observed by Shechmeister [10] and Schaffer et al. [13].

Additionally, a study on aerosol transmission between guinea pigs [3] indicated that transmission was inversely related to RH at 5°C, although experiments at 20°C showed a lower transmission rate at 50% RH than at 65% RH. As 0% transmission was observed at 80% RH, the inconsistent result at 50% or 65% RH may be due to the stochastic nature of infection. If higher transmission rates are due to higher viabilities, at least in part, these results appear to agree with the trend reported by Hemmes [2], [11] and Harper [12]. Given the similar constitution of droplets emitted from infected human and guinea pigs (i.e., salts plus proteins), it seems reasonable to believe that IAVs associated with droplets expelled from humans will be subject to higher inactivation at higher RHs.

The relationship between IAV viability and RH may be due to interactions among components of respiratory droplets (i.e., glycoproteins, salts, and water) and the virus that are sensitive to concentration, which depends on the extent of evaporation, which depends in turn on ambient humidity. Proteins may complicate the effect of salt ions on IAVs by interacting with the salt ions and counteracting their adverse effects. Studies have shown that IAVs remain infectious much longer in the presence of respiratory mucus [37], [38]. Investigation into such interactions and the possible complexes formed in respiratory droplets in response to humidity variation at a molecular level is needed.

We speculate that the seasonality of influenza with its wintertime peak in temperate regions is stimulated by more vigorous evaporation of droplets at low RHs leading to higher suspended concentrations of IAVs, combined with the sensitivity of aerosolized IAVs' viability to RH. When RH is <90%, droplets shrink approximately in half, leaving associated IAVs that can remain suspended long enough to cause secondary infections. Our recent measurements of size-resolved airborne IAV concentrations support this assertion: 64% of the IAV genomes detected in a daycare center, a health center, and airplanes were associated with fine particles <2.5 µm (15% in the 0.25–0.5 µm fraction, 10% in the 0.5–1.0 µm fraction, and 28% in the 1.0–2.5 µm fraction) [33]. These particles can remain suspended for hours to days. Because of the many factors involved in infection, it is still not clear which size of droplets is most likely to transmit influenza, nor is it clear which region in human airways is most susceptible to influenza infection. However, if we simply consider deposition efficiency in human airways, the droplet size with the highest deposition efficiency (~95%) in all regions of the airways combined is ~5 µm. For such droplets, deposition efficiency is ~10% in the tracheobrochial and alveolar regions; the majority of the droplets deposit in the nasopharyngeal region. The droplet size with the highest deposition efficiency (~17%) in the tracheobronchial and alveolar regions is ~2.5 µm [39]. Thus these smaller droplets have greater potential both to remain suspended and to deposit deeper into human airways.

At extremely high RHs, for example, close to 100% in tropical regions during the rainy season, the droplets do not shrink as much (D_eq/D_i

0.927 at 99% RH, and 0.755 at 98% RH, according to our calculations). Droplets thus settle more quickly, rendering the aerosol route relatively less important. However, due to less evaporation, salts and glycoproteins remain at concentrations closer to those found in the respiratory tract, and these concentrations are not detrimental to the virus. As suggested by Lowen et al. [40], [41], other transmission routes (e.g., contact) may dominate in the tropics. They also proposed that the airborne route's sensitivity to RH and temperature contributes to seasonality in temperate regions while the contact route's insensitivity to the two variables contributes to year-round influenza in tropical regions. Our analysis supports this hypothesis.

Model limitations

There are several limitations of our model. First, although the model used to predict equilibrium droplet sizes has been confirmed with experiments using NaCl-bovine serum albumin (BSA) particles [23], further verification with actual respiratory fluid is needed due to its complex composition. Furthermore, the composition of respiratory fluid depends on the emission site (nose or mouth) and source (upper or lower respiratory tract), as well as the stage of infection. Inflamed airways secrete larger amounts of mucus which consequently increase the dry mass of respiratory fluid [42]. Therefore, the equilibrium size of emitted droplets may be larger than presented here based on composition under healthy conditions. On the other hand, saliva has much lower concentrations of salts and glycoproteins [43], [44], due to dilution by which droplets emitted from coughing may have lower dry mass.

Second, the model is based on limited data obtained from laboratory experiments. Not only are Harper [12] and other studies of IAV viability in aerosolized droplets [10], [11], [13] decades old, but none investigated inactivation rates as a function of droplet size. More accurate measurements concerning the influence of respiratory droplet size on IAV viability are needed to better predict the fate of airborne IAVs.

Third, we calculated the IAV concentration based on a well-mixed room model with no recirculation. This model assumes that droplets are instantaneously, continuously, and evenly distributed throughout the room. However, according to Lai and Cheng [45], it takes at least 270 s for 10-µm droplets to mix thoroughly at 5 ACH. It may take even longer for the system to become well-mixed at lower AERs. More accurate calculations may be achieved by the use of computational fluid dynamics. Additionally, if recirculation accounts for a large fraction of the AER and if viruses are not removed in the HVAC system, then ventilation will play a relatively smaller role in virus removal compared to settling and inactivation.

Finally, this research demonstrates the evolution of IAV concentrations induced by a cough but not other activities, such as normal respiration, talking, and sneezing. Fabian et al. [46] determined that IAV RNA (i.e., potentially infectious IAVs) was emitted in exhaled breath from infected patients at a rate of <3.2 to 20 RNA particles min⁻¹. These droplets were smaller than those associated with a cough; over 87% of them were <1 µm in diameter. Therefore, for IAVs exhaled during normal breathing, because of the smaller droplet size, airborne IAV concentrations would be lower and removal would rely more on ventilation and inactivation than gravitational settling. IAVs may also be expelled during talking, although to our knowledge, no detailed experimental data on this phenomenon are currently available in the literature. The droplet size distribution for talking is similar to that for coughing (Table 1). Therefore, model results are expected to be similar for IAVs generated by talking in terms of removal efficiencies by different mechanisms. Sneezing is a less common clinical manifestation of influenza than is coughing, which is manifested in ~80% of patients [21], [46], [47], [48]. The main difference between sneezing and coughing is that the former generates far more droplets, especially smaller ones. Thus, IAV concentrations would be higher initially, and ventilation would play a larger role in removal.

Methods

Equations for generating the initial respiratory droplet size distributions

From data on the size distribution of droplets expelled in a cough [20], we considered counts of droplets with D_i≤100 µm. The standard normal distribution z-value with the same cumulative probability as that for droplets with a diameter D_i was computed by the NORMSINV function in Excel 2007. The equation of the least-squares linear regression between z and lnD_i is shown in Eq. 1 and Figure 1.

The initial size distribution of droplets (≤100 µm) from coughing was then generated by the NORMSDIST function in Excel 2007:

A mathematical equation, expression, or formula.
Object name is pone.0021481.e004.jpg

(4)

where n_i is the droplet count in the ith size bin (5-µm step in this study), N is the total number of droplets ≤100 µm (i.e., 4775 according to Duguid [20]), and z_{i, upper} and z_{i, lower} are the upper and lower z values of the ith size bin.

Model for calculating equilibrium respiratory droplet size

The equilibrium droplet sizes resulting from evaporation were estimated based on Köhler theory taking into account the two major constituents of respiratory fluid: inorganic salts and glycoproteins. Effros et al. [7] determined concentrations (mean ± standard error) of the major electrolytes to be, respectively, 91±8 (Na), 60±11 (K), and 102±17 (Cl) mM, of glycoproteins to be 76.3±18.2 g L⁻¹, and of lactate to be 44±17 mM. We thereby assume respiratory fluid contains 150 mM (8.8 g L⁻¹) NaCl to represent the inorganic components and 76 g L⁻¹ of total proteins (TP) to approximate the organic components, as done by Nicas et al. [22].

The SS-VA model derived by Mikhailov et al. [23] is based on the physiochemical properties (practical osmotic coefficients, molecular weights, and densities of the component solutes, etc.) of the droplet and the Kelvin effect. Their modeling results for particles with 90% BSA (dry mass fraction) fitted well with experimental data for dehydration of mixed NaCl-BSA particles. Given the similar composition of respiratory fluid (89.6% TP in dry mass) to their NaCl-BSA particles, we applied their SS-VA model to compute the equilibrium size for respiratory droplets.

The SS-VA model predicts the equilibrium RH with a specific droplet diameter (D_eq) to be:

A mathematical equation, expression, or formula.
Object name is pone.0021481.e005.jpg

(5)

where, σ is the surface tension (approximated by that of water as done in Mikhailov et al. [23], i.e., 0.072 N m⁻¹); M is the molar mass, the subscripts w and y refer to water and component y (either NaCl or TP), respectively, and M_w

18 g mol⁻¹, M_NaC_l

58.4 g mol⁻¹, M_TP≈M_BSA

66.5×10³ g mol⁻¹; ρ, ρ_w, ρ_y are the densities of the entire droplet, water, and component y (ρ_NaCl

2165 kg m⁻³, ρ_TP≈ρ_BSA

1362 kg m⁻³), respectively; R is the ideal gas constant; T is the absolute temperature (298 K in this study); D_eq is the equilibrium diameter of a droplet residue at a given RH; D_m,s is the mass equivalent diameter of a particle consisting of the dry solutes; υ_y is the stoichiometric dissociation number of component y, υ_NaC_l

2, and υ_TP

υ_BSA

1; Φ_y is the molal or practical osmotic coefficient of component y describing the non-ideality of the solution; and x_s,y is the mass fraction of component y (x_NaCl

0.104 and x_TP

0.896 in this study). Given D_i, the mass equivalent diameter, D_m,s, can be calculated and used as an input to further calculate D_eq with its equilibrium RH.

Virus inactivation rate

Harper [12] performed a detailed study on the viability of airborne IAVs over a wide range of both RH and temperature. In the experiment, droplets containing IAVs were generated with an atomizer and stored in a drum turning at 3 rpm, and results were corrected for physical loss by settling and other deposition mechanisms. We used his viability data at 20–24.5°C, typical indoor temperatures, at RHs ranging from 20% to 81% to calculate inactivation rates as a function of RH (Table 3). Because the residence time of air indoors is typically 1–2 h at most, we considered viability data from the first 1 h of the experiment only.

Table 3

Inactivation of airborne IAVs at 20–24.5°C over 1 h.

We quantified viability by assuming that airborne IAVs undergo first-order inactivation upon emission, such that,

A mathematical equation, expression, or formula.
Object name is pone.0021481.e007.jpg

(6)

where N is the number of IAVs emitted, t is time, and the inactivation rate (k) is

A mathematical equation, expression, or formula.
Object name is pone.0021481.e008.jpg

(7)

where N₀ and N_t are the numbers of IAVs at t

0 and time t, and S_t is the survival rate, or viability (%) at time t. Accordingly, we computed k for each RH from the S_t data reported by Harper [12] using the SLOPE function in Excel 2007 (Table 3). The equation of the least-squares linear regression between k and RH is shown in Eq. 2.

Concentration of infectious IAVs indoors

The model for estimating the concentration of infectious IAVs assumes that they are emitted from a cough and instantaneously well-mixed within the whole indoor space such that IAV concentrations in the room and outlet air are the same. The IAVs are subjected to removal by ventilation, inactivation, and gravitational settling. Droplet size transformation is assumed complete at time zero, and D_eq was used for the calculation. Assuming the inlet air contains no IAVs, the change of IAV concentration with time is modeled as:

A mathematical equation, expression, or formula.
Object name is pone.0021481.e009.jpg

(8)

where C is the infectious IAV concentration in the room and outflow (# m⁻³); v is the gravitational settling velocity calculated by Stokes law based on D_eq; H is assumed to be 2.5 m in this study; λ is the air exchange rate assuming no recirculation; and k is the inactivation rate given by Eq. 2.

At time zero, IAVs are released from a cough, and the initial concentration of IAVs associated with droplets in the ith size bin, C_i,0 (# m ⁻³) is

A mathematical equation, expression, or formula.
Object name is pone.0021481.e010.jpg

(9)

where 0.778 is the initial survival rate (i.e., the average survival rate at 1 s according to the results of Harper [12]); 6.3×10⁻³ is the IAV concentration in respiratory fluid (# µm⁻³) obtained by assuming that the respiratory fluid contains the same concentration of IAV as in nasal washes of infected persons (6.3±0.3×10⁶ median tissue culture infectious dose (TCID₅₀) mL⁻¹ [49]) and that 1 TCID₅₀ equals 1000 virus particles [50]; the term in the first set of parentheses is the mean droplet volume for the bin; D_i,upper and D_i,lower are, respectively, the upper and lower diameters of the ith size bin (µm); n_i is the droplet count for the ith size bin given in Eq. 4, and V is the room volume (assumed to be 50 m³). The solution for C is given in Eq. 3.

Removal efficiency of settling, ventilation, and inactivation

Removal efficiency in this study refers to the percentage of IAVs removed by a certain mechanism (i.e., settling, ventilation, inactivation, or a combination of these three) at a given time and RH. In Figure 5, removal efficiencies of settling (E_settling), ventilation (E_vent), and inactivation (E_inactivation), and total removal efficiency (E_total) are calculated by Eq. 10–13:

A mathematical equation, expression, or formula.
Object name is pone.0021481.e011.jpg

(10)

A mathematical equation, expression, or formula.
Object name is pone.0021481.e012.jpg

(11)

A mathematical equation, expression, or formula.
Object name is pone.0021481.e013.jpg

(12)

A mathematical equation, expression, or formula.
Object name is pone.0021481.e014.jpg

(13)

Acknowledgments

We thank Dr. E. Subbiah for reviewing this work and offering insightful suggestions to improve it.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: This research was supported by the National Science Foundation under grant number CBET-0547107. No additional external funding was received for this study. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. Tamerius J, Nelson M, Zhou S, Viboud C, Miller M, et al. Global influenza seasonality: Reconciling patterns across temperate and tropical regions. Environ Health Perspect 2010 [PMC free article] [PubMed]

2. Hemmes JH, Winkler KC, Kool SM. Virus survival as a seasonal factor in influenza and poliomyelitis. Nature. 1960;188:430–431. [PubMed]

3. Lowen AC, Mubareka S, Steel J, Palese P. Influenza virus transmission is dependent on relative humidity and temperature. PLoS Pathog. 2007;3:e151. [PMC free article] [PubMed]

4. Steel J, Palese P, Lowen AC. Transmission of a 2009 pandemic influenza virus shows a sensitivity to temperature and humidity similar to that of an H3N2 seasonal strain. J Virol. 2011;85:1400–1402. [PMC free article] [PubMed]

5. Xie X, Li Y, Chwang ATY, Ho PL, Seto WH. How far droplets can move in indoor environments–revisiting the wells evaporation–falling curve. Indoor Air. 2007;17:211–225. [PubMed]

6. Chao CYH, Wan MP, Morawska L, Johnson GR, Ristovski ZD, et al. Characterization of expiration air jets and droplet size distributions immediately at the mouth opening. J Aerosol Sci. 2009;40:122–133.

7. Effros RM, Hoagland KW, Bosbous M, Castillo D, Foss B, et al. Dilution of respiratory solutes in exhaled condensates. Am J Respir Crit Care Med. 2002;165:663–669. [PubMed]

8. Raphael GD, Jeney EV, Baraniuk JN, Kim I, Meredith SD, et al. Pathophysiology of rhinitis. Lactoferrin and lysozyme in nasal secretions. J Clin Invest. 1989;84:1528–1535. [PMC free article] [PubMed]

9. Seinfeld JH, Pandis SN. Atmospheric chemistry and physics - from air pollution to climate change. New York: John Wiley & Sons; 2006.

10. Shechmeister IL. Studies on the experimental epidemiology of respiratory infections: III. Certain aspects of the behavior of type A influenza virus as an air-borne cloud. J Infect Dis. 1950;87:128–132. [PubMed]

11. Hemmes J, Winkler K, Kool S. Virus survival as a seasonal factor in influenza and poliomyelitis. Antonie Van Leeuwenhoek. 1962;28:221–233. [PubMed]

12. Harper GJ. Airborne micro-organisms: Survival tests with four viruses. J Hyg. 1961;59:479–486. [PMC free article] [PubMed]

13. Schaffer FL, Soergel ME, Straube DC. Survival of airborne influenza virus: Effects of propagating host, relative humidity, and composition of spray fluids. Arch Virol. 1976;51:263–273. [PubMed]

14. Shaman J, Kohn M. Absolute humidity modulates influenza survival, transmission, and seasonality. Proc Natl Acad Sci USA. 2009;106:3243–3248. [PubMed]

15. Loudon RG, Roberts RM. Droplet expulsion from the respiratory tract. Am Rev Respir Dis. 1967;95:435–442. [PubMed]

16. Papineni RS, Rosenthal FS. The size distribution of droplets in the exhaled breath of healthy human subjects. J Aerosol Med. 1997;10:105–116. [PubMed]

17. Yang S, Lee GWM, Chen C-M, Wu C-C, Yu K-P. The size and concentration of droplets generated by coughing in human subjects. J Aerosol Med. 2007;20:484–494. [PubMed]

18. Morawska L, Johnson GR, Ristovski ZD, Hargreaves M, Mengersen K, et al. Size distribution and sites of origin of droplets expelled from the human respiratory tract during expiratory activities. J Aerosol Sci. 2009;40:256–269.

19. Xie X, Li Y, Sun H, Liu L. Exhaled droplets due to talking and coughing. J R Soc Interface. 2009;6:S703–S714. [PMC free article] [PubMed]

20. Duguid JP. The size and the duration of air-carriage of respiratory droplets and droplet-nuclei. J Hyg. 1946;44:471–479. [PMC free article] [PubMed]

21. Tang JWT, Tambyah PA, Lai FYL, Lee HK, Lee CK, et al. Differing symptom patterns in early pandemic vs seasonal influenza infections. Arch Intern Med. 2010;170:861–867. [PubMed]

22. Nicas M, Nazaroff WW, Hubbard A. Toward understanding the risk of secondary airborne infection: Emission of respirable pathogens. J Occup Environ Hyg. 2005;2:143–154. [PubMed]

23. Mikhailov E, Vlasenko S, Niessner R, Poschl U. Interaction of aerosol particles composed of protein and salts with water vapor: Hygroscopic growth and microstructural rearrangement. Atmos Chem Phys. 2004;4:323–350.

24. Posada JA, Redrow J, Celik I. A mathematical model for predicting the viability of airborne viruses. J Virol Methods. 2010;164:88–95. [PubMed]

25. Tang IN, Tridico AC, Fung KH. Thermodynamic and optical properties of sea salt aerosols. J Geophys Res. 1997;102:23269–23275.

26. Bagger H, Fuglsang C, Westh P. Hydration of a glycoprotein: Relative water affinity of peptide and glycan moieties. Eur Biophys J. 2006;35:367–371. [PubMed]

27. Jensen MM. Inactivation of airborne viruses by ultraviolet irradiation. Appl Microbiol. 1964;12:418–420. [PMC free article] [PubMed]

28. Yamamoto N, Shendell DG, Winer AM, Zhang J. Residential air exchange rates in three major US metropolitan areas: Results from the relationship among indoor, outdoor, and personal air study 1999–2001. Indoor Air. 2010;20:85–90. [PubMed]

29. Tompkins JA, White JA, Bozer YA, Tanchoco JMA. 2010. Facilities planning: John Wiley & Sons, Inc.

30. Myatt TA, Kaufman MH, Allen JG, MacIntosh DL, Fabian MP, et al. Modeling the airborne survival of influenza virus in a residential setting: The impacts of home humidification. Environ Health. 2010;9:55. [PMC free article] [PubMed]

31. Blachere FM, Lindsley WG, Pearce TA, Anderson SE, Fisher M, et al. Measurement of airborne influenza virus in a hospital emergency department. Clin Infect Dis. 2009;48:438–440. [PubMed]

32. Lindsley WG, Blachere FM, Davis KA, Pearce TA, Fisher MA, et al. Distribution of airborne influenza virus and respiratory syncytial virus in an urgent care medical clinic. Clin Infect Dis. 2010;50:693–698. [PubMed]

33. Yang W, Subbiah E, Marr LC. Concentrations and size distributions of airborne influenza A viruses measured indoors at a health centre, a day-care centre, and on aeroplanes. J R Soc Interface. 2011 In press. [PMC free article] [PubMed]

34. California Building Standards Commission. 2007. 2007 California mechanical code (part 4, title 24, California code of regulations)

35. McDevitt J, Rudnick S, First M, Spengler J. Role of absolute humidity in the inactivation of influenza viruses on stainless steel surfaces at elevated temperatures. Appl Environ Microbiol. 2010;76:3943–3947. [PMC free article] [PubMed]

36. Brown JD, Goekjian G, Poulson R, Valeika S, Stallknecht DE. Avian influenza virus in water: Infectivity is dependent on pH, salinity and temperature. Vet Microbiol. 2009;136:20–26. [PubMed]

37. Parker ER, Dunham WB, MacNeal WJ. Resistance of the melbourne strain of influenza virus to desiccation. J Lab Clin Med. 1944;29:37–42.

38. Thomas Y, Vogel G, Wunderli W, Suter P, Witschi M, et al. Survival of influenza virus on banknotes. Appl Environ Microbiol. 2008;74:3002–3007. [PMC free article] [PubMed]

39. Oberdorster G, Oberdorster E, Oberdorster J. Nanotoxicology: An emerging discipline evolving from studies of ultrafine particles. Environ Health Perspect. 2005;113:823–839. [PMC free article] [PubMed]

40. Lowen AC, Steel J, Mubareka S, Palese P. High temperature (30°C) blocks aerosol but not contact transmission of influenza virus. J Virol. 2008;82:5650–5652. [PMC free article] [PubMed]

41. Lowen A, Palese P. Transmission of influenza virus in temperate zones is predominantly by aerosol, in the tropics by contact: A hypothesis. PLoS Curr. 2009;1:RRN1002. [PMC free article] [PubMed]

42. Fischer H, Widdicombe JH. Mechanisms of acid and base secretion by the airway epithelium. J Membr Biol. 2006;211:139–150. [PMC free article] [PubMed]

43. Enberg N, Alho H, Loimaranta V, Lenander-Lumikari M. Saliva flow rate, amylase activity, and protein and electrolyte concentrations in saliva after acute alcohol consumption. Oral Surg Oral Med Oral Pathol Oral Radiol Endod. 2001;92:292–298. [PubMed]

44. Michishige F, Kanno K, Yoshinaga S, Hinode D, Takehisa Y, et al. Effect of saliva collection method on the concentration of protein components in saliva. J Med Invest. 2006;53:140–146. [PubMed]

45. Lai ACK, Cheng YC. Study of expiratory droplet dispersion and transport using a new eulerian modeling approach. Atmos Environ. 2007;41:7473–7484.

46. Fabian P, McDevitt JJ, DeHaan WH, Fung ROP, Cowling BJ, et al. Influenza virus in human exhaled breath: An observational study. PLoS ONE. 2008;3:e2691. [PMC free article] [PubMed]

47. Douglas RG., Jr . Influenza in man. In: Kilbourne ED, editor. The influenza viruses and influenza. New York: Academic Press; 1975. pp. 395–447.

48. Lindsley WG, Blachere FM, Thewlis RE, Vishnu A, Davis KA, et al. Measurements of airborne influenza virus in aerosol particles from human coughs. PLoS ONE. 2010;5:e15100. [PMC free article] [PubMed]

49. Murphy BR, Rennels MB, Douglas RG, Jr, Betts RF, Couch RB, et al. Evaluation of influenza A/Hong Kong/123/77 (H1N1) ts-1A2 and cold-adapted recombinant viruses in seronegative adult volunteers. Infect Immun. 1980;29:348–355. [PMC free article] [PubMed]

50. Ward CL, Dempsey MH, Ring CJA, Kempson RE, Zhang L, et al. Design and performance testing of quantitative real time PCR assays for influenza A and B viral load measurement. J Clin Virol. 2004;29:179–188. [PubMed]

51. Hinds WC. Aerosol technology. 2nd ed. New York: Wiley; 1999. pp. 75–110.

Articles from PLoS ONE are provided here courtesy of
Public Library of Science