Abstract
The unprecedented pandemic of COVID-19 has created worldwide shortages of personal protective equipment, in particular respiratory protection such as N95 respirators1. SARS-CoV-2 transmission is frequently occurring in hospital settings, with numerous reported cases of nosocomial transmission highlighting the vulnerability of healthcare workers2–4. In general, N95 respirators are designed for single use prior to disposal. Several groups have addressed the potential for re-use of N95 respirators from a mechanical or from a decontamination perspective (for a full literature overview see Supplementary Appendix).
Here, we analyzed four different decontamination methods – UV radiation (260 – 285 nm), 70°C heat, 70% ethanol and vaporized hydrogen peroxide (VHP) – for their ability to reduce contamination with infectious SARS-CoV-2 and their effect on N95 respirator function. For each of the decontamination methods, we compared the inactivation rate of SARS-CoV-2 on N95 filter fabric to that on stainless steel, and we used quantitative fit testing to measure the filtration performance of the N95 respirators after each decontamination run and 2 hours of wear, for three consecutive decontamination and wear sessions (see Appendix). Vaporized hydrogen peroxide and ethanol yielded extremely rapid inactivation both on N95 and on stainless steel (Figure 1A). UV inactivated SARS-CoV-2 rapidly from steel but more slowly on N95 fabric, likely due its porous nature. Heat caused more rapid inactivation on N95 than on steel; inactivation rates on N95 were comparable to UV.
Decontamination of SARS-CoV-2 by four different methods. A) SARS-CoV-2 on N95 fabric and stainless steel surface was exposed to UV, 70 °C dry heat, 70% ethanol and vaporized hydrogen peroxide (VHP). 50 µl of 105 TCID50/mL of SARS-CoV was applied on N95 and stainless steel (SS). Samples were collected at indicted timepoints post exposure to the decontamination method for UV, heat and ethanol and after 10 minutes for VHP. Viable virus titer is shown in TCID50/mL media on a logarithmic scale. All samples were quantified by end-point titration on Vero E6 cells. Plots show estimated mean across three replicates (dots and bars show the posterior median estimate of this mean and the posterior inter-quartile range, or IQR). Lines show predicted decay of virus titer over time (lines; 50 random draws per replicate from the joint posterior distribution of the exponential decay rate, i.e. negative of the slope, and intercept, i.e. initial virus titer). Black dashed line shows maximum likelihood estimate titer at the Limit of Detection (LOD) of the assay: 100.5 TCID50/mL media. B) Mask integrity. Quantitative fit testing results for all the decontamination methods after decontamination and 2 hours of wear, for three consecutive runs. Data from six individual replicates (small dots) for each treatment are shown in addition to the predicted median and IQR (large dots and bars respectively) fit factor. Fit factors are a measure of filtration performance: the ratio of the concentration of particles outside the mask to the concentration inside. The measurement machine reports value up to 200. A minimal fit factor of 100 (red dashed line) is required for a mask to pass a fit test. C) SARS-CoV-2 decontamination performance. Kill rate (y-axis), versus mask integrity after decontamination (x-axis; bar length represents IQR). The three panels report mask integrity after one, two or three decontamination cycles.
Quantitative fit tests showed that the filtration performance of the N95 respirator was not markedly reduced after a single decontamination for any of the four decontamination methods (Figure 1B). Subsequent rounds of decontamination caused sharp drops in filtration performance of the ethanol-treated masks, and to a slightly lesser degree, the heat-treated masks. The VHP- and UV-treated masks retained comparable filtration performance to the control group after two rounds of decontamination, and maintained acceptable performance after three rounds.
Taken together, our findings show that VHP treatment exhibits the best combination of rapid inactivation of SARS-CoV-2 and preservation of N95 respirator integrity, under the experimental conditions used here (Figure 1C). UV radiation kills the virus more slowly and preserves comparable respirator function. 70°C dry heat kills with similar speed and is likely to maintain acceptable fit scores for two rounds of decontamination. Ethanol decontamination is not recommended due to loss of N95 integrity, echoing earlier findings5.
All treatments, particularly UV and dry heat, should be conducted for long enough to ensure that a sufficient reduction in virus concentration has been achieved. The degree of required reduction will depend upon the degree of initial virus contamination. Policymakers can use our estimated decay rates together with estimates of degree of real-world contamination to choose appropriate treatment durations (see Appendix).
Our results indicate that N95 respirators can be decontaminated and re-used in times of shortage for up to three times for UV and HPV, and up to two times for dry heat. However, utmost care should be given to ensure the proper functioning of the N95 respirator after each decontamination using readily available qualitative fit testing tools and to ensure that treatments are carried out for sufficient time to achieve desired risk-reduction.
Data Availability
Code and data to reproduce the Bayesian estimation results and produce corresponding figures are archived online at OSF: and available on Github:
Supplemental methods
Short literature review
The COVID-19 pandemic has highlighted the necessity for large-scale decontamination procedures for PPE, in particular N95 respirator masks1. SARS-CoV-2 has frequently been detected on PPE of healthcare workers2. The environmental stability of SARS-CoV-2 underscores the need for rapid and effective decontamination methods3. Extensive literature is available for decontamination procedures for N95 respirators, using either bacterial spore inactivation tests, bacteria or respiratory viruses (e.g. influenza A virus)4-11. Effective inactivation methods for these pathogens and surrogates include UV, ethylene oxide, vaporized hydrogen peroxide, gamma irradiation, ozone and dry heat4,6,8,10-13. The filtration efficiency and N95 respirator fit has typically been less well explored, but suggest that both filtration efficiency and N95 respirator fit can be affected by the decontamination method used12,14. It will therefore be critical that FDA, CDC and OSHA guidelines with regards to fit testing, seal check and respirator re-use are followed4,15-18.
Laboratory experiments
Viruses and titration
HCoV-19 nCoV-WA1-2020 (MN985325.1) was the SARS-CoV-2 strain used in our comparison19. Virus was quantified by end-point titration on Vero E6 cells as described previously20. Virus titrations were performed by end-point titration in Vero E6 cells. Cells were inoculated with 10-fold serial dilutions in four-fold of samples taken from N95 mask and stainless steel surfaces (see below). One hour after inoculation of cells, the inoculum was removed and replaced with 100 µl (virus titration) DMEM (Sigma-Aldrich) supplemented with 2% fetal bovine serum, 1 mM L-glutamine, 50 U/ml penicillin and 50 µg/ml streptomycin. Six days after inoculation, cytopathogenic effect was scored and the TCID50 was calculated (see below). Wells presenting cytopathogenic effects due to media toxicity (e.g., due to the presence of ethanol or hydrogen peroxide) rather than viral infection were removed from the titer inference procedure.
N95 and stainless steel surface
N95 material discs were made by punching 9/16” (15 mm) fabric discs from N95 respirators, AOSafety N9504C respirators (Aearo Company Southbridge, MA). The stainless steel 304 alloy discs were purchased from Metal Remnants (https://metalremnants.com/) as described previously. 50 µL of SARS-CoV-2 was spotted onto each disc. A 0 time-point measurement was taken prior to exposing the discs to the disinfection treatment. At each sampling time-point, discs were rinsed 5 times by passing the medium over the stainless steel or through the N95 disc. The medium was transferred to a vial and frozen at −80°C until titration. All experimental conditions were performed in triplicate.
Decontamination methods
Ultraviolet light
Plates with fabric and steel discs were placed under an LED high power UV germicidal lamp (effective UV wavelength 260-285nm) without the titanium mesh plate (LEDi2, Houston, Tx) 50 cm from the UV source. At 50 cm the UVAB power was measured at 5 µW/cm2 using a General UVAB digital light meter (General Tools and Instruments New York, NY). Plates were removed at 10, 30 and 60 minutes and 1 mL of cell culture medium added.
Heat treatment
Plates with fabric and steel discs were placed in a 70°C oven. Plates were removed at 10, 20, 30 and 60 minutes and 1 mL of cell culture medium added.
70% ethanol
Fabric and steel discs were placed into the wells of one 24 well plate per time-point and sprayed with 70% ethanol to saturation. The plate was tipped to near vertical and 5 passes of ethanol were sprayed onto the discs from approximately 10 cm. After 10 minutes,, 1 mL of cell culture medium was added.
Vaporized hydrogen peroxide (VHP)
Plates with fabric and steel discs were placed into a Panasonic MCO-19AIC-PT (PHC Corp. of North America Wood Dale, IL) incubator with VHP generation capabilities and exposed to hydrogen peroxide (approximately 1000 ppm). The exposure to VHP was 10 minutes, after the inactivation of the hydrogen peroxide, the plate was removed and 1 mL of cell culture medium was added.
Control
Plates with fabric and steel discs and steel plates were maintained at 21-23°C and 40% relative humidity for up to four days. After the designated time-points, 1 mL of cell culture medium was added.
N95 mask integrity testing
N95 Mask (3M™ Aura™ Particulate Respirator 9211+/37193) integrity testing after 2 hours of wear and decontamination, for three consecutive rounds, was performed for a total of 6 times for each decontamination condition and control condition. Masks were worn by subjects and integrity was quantitatively determined using the Portacount Respirator fit tester (TSI, 8038) with the N95 companion component, following the modified ambient aerosol condensation nuclei counter quantitative fit test protocol approved by the OSHA (Occupational Safety and Health Administration, 2012). Subjects were asked to bend over for 40 seconds, talk for 50 seconds, move head from side-to-side for 50 seconds, and move head up-and-down for 50 seconds whilst aerosols on inside and outside of mask were measured. By convention, this fit test is passed when the final score is ≥100. For the N95 integrity testing, a Honeywell Mistmate humidifier (cat#HUL520B) was used for particle generation.
Statistical analyses
In the model notation that follows, the symbol ∼ denotes that a random variable is distributed according to the given distribution. Normal distributions are parametrized as Normal(mean, standard deviation). Positive-constrained normal distributions (“Half-Normal”) are parametrized as Half-Normal(mode, standard deviation). Normal distributions truncated to the interval [0, 1] are parameterized as TruncNormal(mode, standard deviation).
We use <Distribution Name>CDF(x | parameters) and <Distribution Name>CCDF to denote the cumulative distribution function and complementary cumulative distribution functions of a probability distribution, respectively. So for example NormalCDF(5 | 0, 1) is the value of the Normal(0, 1) cumulative distribution function at 5.
We use logit(x) and invlogit(x) to denote the logit and inverse logit functions, respectively:
Mean titer inference
We inferred mean titers across sets of replicates using a Bayesian model. The log10 titers vijk (the titer\ for the sample from replicate k of timepoint j of experiment i) were assumed to be normally distributed about a mean µij with a standard deviation σ. We placed a very weakly informative normal prior on log10 titers µij:
We placed a weakly informative normal prior on the standard deviation:
We then modeled individual positive and negative wells for sample ijk according to a Poisson single-hit model21. That is, the number of virions that successfully infect cells in a given well is Poisson distributed with mean:
where v is the log10 virus titer in TCID50, where v is the log10 virus titer in TCID50, and the well is infected if at least one virion successfully infects a cell. The value of the mean derives from the fact that our units are TCID50; the probability of infection at v = 0, i.e. 1 TCID50, is equal to 1 – e-ln(2) × 1 = 0.5.
Let Yijkdl be a binary variable indicating whether the lth well of dilution factor d (expressed as log10 dilution factor) of sample ijk was positive (so Yijkdl = 1 if the well was positive and 0 otherwise), which will occur as long as at least one virion successfully infects a cell.
It follows from (5) that the conditional probability of observing Yijkdl = 1 given a true underlying titer log10 titer vijk is given by:
Where
is the expected concentration, measured in log10 TCID50, in the dilute sample. This is simply the probability that a Poisson random variable with mean (– ln(2) × 10x) is greater than 0. Similarly, the conditional probability of observing Yijkdl = 0 given a true underlying titer log10 titer vijk is given by:
which is the probability that the Poisson random variable is 0.
This gives us our likelihood function, assuming independence of outcomes across wells.
Virus inactivation regression
The durations of detectability depend on the decontamination treatment but also initial inoculum and sampling method, as expected. We therefore estimated the decay rates of viable virus titers using a Bayesian regression analogous to that used in van Doremalen et al., 20203. This modeling approach allowed us to account for differences in initial inoculum levels across replicates as well as other sources of experimental noise. The model yields estimates of posterior distributions of viral decay rates and half-lives in the various experimental conditions – that is, estimates of the range of plausible values for these parameters given our data, with an estimate of the overall uncertainty22.
Our data consist of 10 experimental conditions: 2 materials (N95 masks and stainless steel) by 5 treatments (no treatment, ethanol, heat, UV and VHP). Each has three replicates, and multiple time-points for each replicate. We analyze the two materials separately. For each, we denote by Yijkdl the positive or negative status (see above) for well l which has dilution d for the titer vijk from experimental condition i during replicate j at time-point k.
We model each replicate j for experimental condition i as starting with some true initial log10 titer vij(0) = vij0. We assume that viruses in experimental condition i decay exponentially at a rate λi over time t. It follows that:
We use the direct-from-well data likelihood function described above, except that now instead of estimating titer distribution about a shared mean µij we estimate λi under the assumptions that our observed well data Yijkdl reflect the titers vij(t).
Regression prior distributions
We place a weakly informative Normal prior distribution on the initial log10 titers vij0 to rule out implausibly large or small values (e.g. in this case undetectable log10 titers or log10 titers much higher than the deposited concentration), while allowing the data to determine estimates within plausible ranges:
We placed a weakly informative Half-Normal prior on the exponential decay rates λi:
Our plated samples were of volume 0.1 mL, so inferred titers were incremented by 1 to convert to units of log10 TCID50/mL.
Mask integrity estimation
To quantify the decay of mask integrity after repeated decontamination, we used a logit-linear spline Bayesian regression to estimate the rate of degradation of mask fit factors over time, accounting for the fact that fit factors are interval-censored ratios. Fit factors are defined as the ratio of exterior concentration to interior concentration of a test aerosol. They are reported to the nearest integer, up to a maximum readout of 200, but arbitrarily large true fit factors are possible as the mask performance approaches perfect filtration.
We had 6 replicate masks j for each of 5 treatments i (no decontamination, ethanol, heat, UV and VHP). Each mask j was assessed for fit factor at 4 time-points k: before decontamination, and then after 1, 2, and 3 decontamination cycles. We label the control treatment i = 0. So we denote by Fijk the fit factor for the jth mask from the ith treatment after k decontaminations (with k = 0 for the initial value).
We first converted fit factors Fijk to the equivalent observed filtration rate Yijk by:
Observation model and likelihood function
We modeled the censored observation process as follows. logit(Yijk) values are observed with Gaussian error about the true filtration logit(pijk), with an unknown standard deviation σo, and then converted to fit factors, which are then censored:
Because our reported fit factors are known to be within integer values and right-censored at 200, for Fijk ≥ 200 we have a conditional probability of observing the data given the parameters of
That is, we calculate the probability of observing a value of F greater than or equal to 200 (equivalent a value of Y greater than or equal to 1 – 1/200), given our parameters.
For 1.5 ≤ Fijk < 200, we first calculate the upper and lower bounds of our observation Y+ijk = 1 – 1 / (Fijk – 0.5) and Y–ijk = 1 – 1 / (Fijk – 0.5). Then:
That is, we calculate the probability of observing a value between Y+ and Y–, given our parameters.
Decay model
We assumed that each mask had some true initial filtration rate pij0. We assumed that these were logit-normally distributed about some unknown mean mask initial filtration rate pavg with a standard deviation σp, that is:
We then assumed that the logit of the filtration rate, logit(pijk), decreased after each decontamination by a quantity d0k + dik, where d0k is natural degradation during the kth trial in the absence of decontamination (i.e. the degradation rate in the control treatment, i = 0), and dik is the additional degrading effect of the kth decontamination treatment of type i > 0). So for k = 1, 2, 3 and i > 0:
where εijk is a normally-distributed error term with an inferred standard deviation σε:
And for the control i = 0:
Model prior distributions
We placed a weakly informative Half-Normal prior on the control degradation rate d0:
We placed a weakly informative Half-Normal prior on the non-control degradation rates di, i > 0:
reflecting the conservative assumption that decontamination should degrade the mask at least somewhat.
We placed a Truncated Normal prior on the mean initial filtration pavg:
The mode of 0.995 corresponds to the maximum measurable fit factor of 200. The standard deviation of 0.02 leaves it plausible that some masks could start near or below the minimum acceptable threshold fit factor of 100, which corresponds to a p of 0.99.
We placed weakly informative Half-Normal priors on the logit-space standard deviations σp, σε, and σo. σp reflects variation in individual masks’ initial filtration about pavg. σε reflects variation in mask’s true degree of degradation between decontaminations about the expected decay, and σo reflects noise in the observation process.
We chose a standard deviation of 0.5 for the priors because a standard deviation of 1.5 (i.e. 3 σ in the prior) in logit space corresponds to probability values being uniformly distributed between 0 and 1; we therefore wish to tell our model not to use larger standard deviations, as these squash all pijk to one of two modes, one at 0 and one at 123.
Markov Chain Monte Carlo Methods
For all Bayesian models, we drew posterior samples using Stan (Stan Core Team 2018), which implements a No-U-Turn Sampler (a form of Markov Chain Monte Carlo), via its R interface RStan. We ran four replicate chains from random initial conditions for 2000 iterations, with the first 1000 iterations as a warmup/adaptation period. We saved the final 1000 iterations from each chain, giving us a total of 4000 posterior samples. We assessed convergence by inspecting trace plots and examining R□ and effective sample size (neff) statistics.
Supplemental table
Code and data availability
Code and data to reproduce the Bayesian estimation results and produce corresponding figures are archived online at OSF: and available on Github:
Acknowledgements
We would like to thank Madison Hebner, Julia Port, Kimberly Meade-White, Irene Offei Owusu, Victoria Avanzato and Lizzette Perez-Perez for excellent technical assistance. This research was supported by the Intramural Research Program of the National Institute of Allergy and Infectious Diseases (NIAID), National Institutes of Health (NIH). JOL-S and AG were supported by the Defense Advanced Research Projects Agency DARPA PREEMPT # D18AC00031 and the UCLA AIDS Institute and Charity Treks, and JOL-S was supported by the U.S. National Science Foundation (DEB-1557022), the Strategic Environmental Research and Development Program (SERDP, RCL2635) of the U.S. Department of Defense. Names of specific vendors, manufacturers, or products are included for public health and informational purposes; inclusion does not imply endorsement of the vendors, manufacturers, or products by the US Department of Health and Human Services.
