Hot, Humid Air Decontamination of Aircraft Confirmed That High Temperature and High Humidity Are Critical for Inactivation of Infectious, Enveloped Ribonucleic Acid (RNA) Virus

Aims: To develop infectious (live/dead) enveloped virus test indicators and response surface methodology (RSM) models that evaluate survival of an enveloped ribonucleic acid (RNA) virus on contaminated aircraft materials after exposure to hot, humid air (HHA). Methods and Results: Enveloped RNA bacteriophage Phi6 (Φ6) was dried on wiring insulation, aircraft performance coating (APC), polypropylene, and nylon at ≥ 8 log10 plaque-forming units (PFU) test coupon−1. Only 2.4 log10 inactivation was measured on APC at 70°Celsius (°C), 5% relative humidity (RH) after 24 h. In contrast, HHA RSM models showed a 90% probability of a 7 log10 inactivation at ≥63°C, 90% RH after 1 h, and decontamination kinetics were similar across different materials. HHA decontamination of C-130 and C-17 aircraft showed >7 log10 and ≥5.9 log10 inactivation of enveloped virus on 100 and 110 test indicators, respectively, with a 1-h treatment, excluding ramp-up and ramp-down times. Conclusions: Enveloped RNA virus test indicators were successfully developed, lab tested for HHA decontamination, analyzed for RSM, and field-tested in aircraft demonstrations. Significance and Impact of the Study: The utility of HHA decontamination was demonstrated after inactivating enveloped RNA virus on aircraft with a 1-h HHA treatment within aircraft temperature and RH limits.


INTRODUCTION
There is a need to develop aircraft decontamination procedures for enveloped RNA viruses ranging from hemorrhagic fever viruses to coronaviruses (Centers for Disease Control Prevention, 2019; National Business Aviation Association, 2020). A candidate procedure for aircraft decontamination is hot, humid air (HHA), which has been tested for Bacillus spore decontamination (Buhr et al., 2012(Buhr et al., , 2015(Buhr et al., , 2016. HHA decontamination is highly feasible for enveloped virus decontamination since enveloped viruses are considered much easier to kill than spores (Spaulding, 1957). In order to develop HHA decontamination profiles and test in aircraft field demonstrations, a biosafety level 1 (BSL-1) enveloped RNA virus was needed as a test indicator to meet the requirements imposed by biosafety and environmental reviews for field testing. Field test approval is a lengthy process that can require decades of compiled research to approve BSL-1 organisms for field testing, regardless of whether the organism is contained or released during the test (Bishop and Robinson, 2014;Buhr et al., 2016). It is particularly challenging to justify and approve enveloped virus strains for rapid field testing (2 months for biosafety approval at different facilities, enveloped virus test indicator preparation, and 2 aircraft field tests including results) during a pandemic such as the SARS-CoV-2/COVID-19 pandemic. Furthermore, active, infectious virus was needed for a live/dead (infectious/non-infectious for viruses) assay since live/dead assays are the hallmark of decontamination testing.
6 is a BSL-1 enveloped RNA virus originally isolated in a bean field as a lytic virus that infected the plant pathogenic bacterium Pseudomonas syringae pathovar phaseolicola (Vidaver et al., 1973;Van Etten et al., 1976;Mindich, 2004). Early hopes to produce the virus in large quantities to spray on bean fields as an environmentally friendly biocontrol agent were never commercially realized. However, due to its rare combination as a BSL-1, enveloped RNA virus it has been proposed as a general surrogate for a number of different enveloped RNA viruses particularly for field testing (Gallandat and Lantagne, 2017). The 6 envelope structure is similar to many other enveloped viruses as the envelope consists of a glycoprotein/protein-embedded lipid membrane, and the host cell has similar temperature sensitivity to mammalian cells at around 40 • C. This is important since the envelope components are considered a primary target for inactivation by many different decontaminants (McDonnell and Burke, 2011;Wiggington et al., 2012). Here 6 was prepared to develop BSL-1 enveloped virus test indicators at ≥ 7.6 log 10 coupon −1 , test HHA decontamination parameters, develop response surface models for HHA decontamination, and measure enveloped virus inactivation during C-130 and C-17 aircraft field tests.
6 was prepared after inoculating broth cultures of HB10Y. A frozen stock prep of HB10Y was thawed at 22 ± 3 • C. HB10Y was added either directly from a frozen stock or by transferring a single colony from a streaked TSA plate to 100-200 ml of 3% TSB in a 1-l smooth-bottom Erlenmeyer flask with a HEPA cap and incubated at 26 ± 2 • C, 200 rev min −1 overnight. Cells were then diluted and grown to mid-log phase. The host flask was inoculated with 0.5-1 ml of 6 at a concentration of ∼1-2e11 plaque-forming units (PFU) ml −1 at the time of inoculation. The culture was incubated at 26 ± 2 • C, 200 rev min −1 for 24 ± 2 h. The 6 preparation was stored at 4 • C until after titering was completed. After titer determination was completed, typically around 1-2e11 PFU ml −1 , then 1-1.3 ml volumes were aliquoted into 1.5-ml screw-cap tubes with O-rings, inverted, and stored at −80 • C.

Environmental Test Chamber Setup and Validation
Thermotron SM-8-8200 (Thermotron Industries, Holland, MI, USA) environmental chambers were used to control temperature and RH as described (Buhr et al., 2015).

Coupon Materials and Sterilization
Square 2 centimeter (cm) × 2 cm coupons of different test materials or the inside surface of 50-ml Techno Plastic Products AG, Switzerland (TPP R ), polypropylene conical tubes were inoculated with >8 log 10 6 virus inoculum. The inoculated 2 cm × 2 cm coupons were set inside sterile TPP R tubes during testing (Buhr et al., 2012(Buhr et al., , 2015(Buhr et al., , 2016. Materials for inoculation included aluminum 2024-T3 coupons painted with water-based aircraft performance coating (APC), stainless-steel 304 coupons painted with Navy Top Coat (NTC) [Coatings Group at the University of Dayton Research Institute (Dayton, OH, USA)], wiring insulation (Kapton Film Type HN, 1 mil, Reference No. 6197844-00 from Cole Parmer, Vernon Hills, IL, USA), and nylon webbing (nylon) (US Netting, Erie, PA, USA) with the ends of each coupon cauterized to prevent fraying. Polypropylene plastic and wiring insulation represent non-porous materials found on commercial and military aircraft. APC and NTC represent semiporous surfaces found on military aircraft and ships, with APC the primary focus for this testing. Nylon webbing represents porous materials found on military and commercial aircraft since nylon is used to manufacture various carpets and fabrics. Prior to inoculating and testing, coupons were rinsed with 18 mega-Ohm-cm, de-ionized water, placed on absorbent paper in an autoclave-safe container, and autoclaved for 30 min at 121 • C, 100 kilopascals. Autoclaved coupons were stored in sterile containers until used. Microbial recovery from autoclaved APC coupons was variable. There was no visible damage, changes in contact angle measurements, or evidence of surface variability on autoclaved APC coupons after inspection with a scanning electron microscope (data not shown). In order to clean and remove interferents such as dried, residual paint solvent, APC coupons were soaked in 70% ethanol for 20 min, rinsed at least 3 times with sterile 18 mega-ohm-cm water, and dried.

Response Surface Methodology (RSM)
RSM was used for previous work with HHA decontamination of materials contaminated with spores because the RSM Design of Experiments (DOE) was requested by the funding agencies as an industry-accepted method to describe HHA points of failure (Beauregard et al., 1992;Lenth, 2009;Myers et al., 2009;Buhr et al., 2012Buhr et al., , 2015. RSM required three equally spaced test parameters for each variable; time, temperature, and humidity. An extensive amount of trial-and-error testing showed a very rapid and dramatic impact of temperature and RH on virus decontamination. Test temperatures were 55, 60, and 65 • C. Test RH was 50, 70, and 90%. Test times were 1, 5, and 9 h to provide a practical timeframe for decontamination of ≤24 h. RSM DOE test parameters are displayed in Figure 1. The first 13 test runs (first iteration) included the center and corners of the experimental test box. The latter 6 test runs (second iteration) included the face-points at the center of each side of the test box.

Virus Test Method and Design
An overlay procedure used for 6 was adapted from procedures used for Escherichia coli and its bacteriophages using TSA and TSB media instead of Luria-Bertani agar and broth. Plates were incubated at 26 ± 2 • C for 20 ± 2 h before quantifying plaques. Plates were then incubated another 24 ± 2 h at RT and scored a final time.
6 stock preps were removed from −80 • C storage and thawed at RT prior to preparing inoculum by transferring stock 6 into 50-ml conical tubes containing 10 mM HEPES + 10% sucrose pH 7.0 for a final concentration of 1-2e9 PFU ml −1 . The infectious virus concentration was confirmed by titering. Test coupons and conical tubes were inoculated with 0.1 ml of 6 inoculum, held at RT for 24 ± 2 h to dry allowing virus to adhere to the test material, then transferred to TPP R tubes. Following exposure of 6 to various hot humid air decontamination testing parameters, samples were extracted and plated.
For 6 extraction, 5 ml of 10 mM HEPES + 10% sucrose pH 7 were added to each conical tube with a virus-inoculated coupon and vortexed for 2 min. After vortexing, 5 ml of HB10Y log-phase culture was added and allowed to infect at RT for 15 min, followed by 2 min of vortexing. Each sample was serially diluted in 900 µl of 10 mM HEPES + 10% sucrose pH 7 out to the −6 dilution. For each 6 dilution, 1,000 µl for the first dilution and 200 µl for subsequent dilutions were transferred into individual tubes containing 200 µl log-phase HB10Y. Then, 1,200 or 200 µl of those 6/HB10Y mixtures was added to individual TSB overlay tubes, poured onto individual TSA plates, and allowed to solidify for ≥30 min. Solidified plates were then inverted, incubated for 20 ± 2 h at 26 • C, and quantified. Plates were incubated another 24 h at RT and quantified a final time. Quantitation and calculations of survival were performed as previously described (Buhr et al., 2012(Buhr et al., , 2014. The maximum (100%) recoverable virus reference, the inoculum titer, was used to calculate extraction efficiency. Extraction efficiency was the number of PFU removed from non-treated control coupons divided by the number of PFU in the inoculum.

Response Surface Methodology (RSM) Performance Envelopes-Mathematical Model Development
RSM to represent survivability of 6 bacteriophage were developed. For a given DOE point, sample data varied with respect to initial 6 bacteriophage population and a final surviving 6 bacteriophage population. To facilitate model development, the results of the samples were normalized into a Survival Fraction (SF). A SF is simply the surviving population divided by the initial population. SF is assumed to never be greater than one or less than zero.
The smallest non-zero SF that may be observed is limited due to the initial population being a finite number. Specifically, the smallest observable SF is 1 divided by the initial population. However, if the results reflect that of a 1-ml sample taken from a 10-ml sample (as is true for this data), then the smallest observable SF is 10 divided by the initial population. If extraction efficiency is <100%, then the smallest observable SF becomes 10 divided by the initial population divided by the extraction efficiency. If test data shows a surviving population of zero, it is likely because the theoretical SF is smaller than the smallest observable SF. Because of the limitation imposed by a finite initial population, any test result showing a survival population of zero was treated as a SF that is less than or equal to the smallest observable SF.
The geometric mean survival fraction (GMSF) equation below shows how GMSF was calculated for each DOE point. In Equation 1, "N" represents the number of samples for a given DOE point.
All GMSF values from the experiments lie between "0" and "1, " and they differ from each other by orders of magnitude, which can make it difficult to find a good fit to the data. To avoid this problem and to facilitate the data fitting process, it is desirable to transform the widely varying GMSF values such that the transformed values are distributed in a linear fashion. To achieve this, the natural log of the GMSF values was computed. However, this transform still resulted in a very non-linear distribution. A second transform was performed, using Equation 2 below, that flipped the sign of the first transform and then took its natural log. This yielded a more linear distribution of data.
One benefit of this transform is that it constrains all model predictions to a value between "0" and "1" and thereby prevents a result that contradicts the definition of SF. Data fit exploration focused on obtaining a good fit to the data set of transformed GMSF values. During the data fitting process, a constraint was imposed on the fit such that as temperature and/or RH and/or time increase, SF decreases consistently without exception. The final functional form for predicting GMSF is shown in the system of equations below (Equation 3).
Where: T = temperature ( • C), R = relative humidity (%), H = time (hours) c 0 to c 7 = addition and multiplication coefficients n 1 to n 3 = power coefficients T 0 = temperature offset coefficient R 0 = relative humidity offset coefficient H 0 = time offset coefficient GMSF = geometric mean survival fraction. As is evident in the data, a given exposure condition will result in a distribution of various SF values that can vary geometrically. Because of this geometric variation, it is appropriate to assume the natural log of the survival fractions for a given sample follow a normal distribution for a given exposure condition. An estimate of the expected distribution of results is important because it helps a user identify a probability that an expected total kill will occur for a given exposure condition. For example, a user can utilize the expected distribution to identify the exposure condition required to achieve at least a 7 log 10 kill with a 95% probability.
Most of the exposure DOE points consist of only 5 samples, and this small sample size increases the uncertainty in determining the population geometric standard deviation (GSD) from the sample geometric standard deviation. There was no clear and consistent pattern in the sample geometric standard deviations that would allow for the development of a reasonable equation to predict geometric standard deviation as a function of exposure conditions. However, the DOE center point (Temp = 60, RH = 70, H = 5) consists of 25 samples for each medium and consists of enough data to reasonably estimate a geometric standard deviation that is assumed to be representative for the entire DOE factor space. Using the test results from the DOE center point for each medium, the geometric standard deviation was computed using Equation 4 below.
Where: n = number of samples i = indexed reference to a particular sample SF i = survival fraction for sample "i" GMSF = geometric mean survival fraction GSD = geometric standard deviation. The normal distribution of possible Survival Fraction logs for a given exposure condition is described by the natural log of the GMSF predicted using Equation 3 and the natural log of the GSD computed using Equation 4. The GMSF will vary based on exposure conditions and medium, but the GSD will only vary by medium. The equations for the Normal distribution are provided in Equations 5 and 6 below. Where: Frontiers in Bioengineering and Biotechnology | www.frontiersin.org PDF = probability distribution function. Characterizes the relative frequency of occurrence. CDF = cumulative distribution function. The integral of the PDF with results ranging from 0 to 1. x = horizontal axis of possible Survival Fractions µ = natural log of geometric mean survival fraction (GMSF) σ = natural log of geometric standard deviation (GSD) erf = the error function. The cumulative distribution function equation returns the area for a region under the Probability Distribution Function curve by integrating from -∞ to a specified SF. A CDF result of "0.95" indicates that 95% of the SF values observed will be less than the specified SF.
The Solver add-in for Microsoft Excel 2016 was used to explore various fits to the transformed survival fractions. The Solver add-in has the added benefit of allowing fit constraints which was crucial to obtaining a good and reasonable fit. The "Evolutionary" method in Solver proved a valuable tool in optimizing the coefficients shown in Equation 3. Excel was also used to produce the various plots in this report and assess model performance.

C-130 and C-17 Aircraft Test Demonstrations
Enveloped virus test indicators consisted of APC coupons inoculated with 1 of 5 independent virus preparations at >8 log 10 coupon −1 , and the test coupons were placed inside TPP R tubes prior to storage, shipping, and testing. Test indicators were placed in 20-22 locations with 5 test indicators per location within the aircraft. Each of the 5 test indicators was inoculated with an independent virus preparation. HHA tests were conducted inside a C-130 fuselage located at Army Combat Capabilities Developmental Command, Edgewood, MD, USA (Figure 2), or in a C-17 located at Dover Air Force Base, Dover, DE, USA (Figure 3), using a previously tested HHA system (Buhr et al., 2016) ran and managed/monitored by a consortium listed in the acknowledgments. After HHA exposure, the test indicators were shipped to a laboratory for virus extraction. Shipping controls were shipped, but not treated with HHA. Laboratory controls were not shipped. Both control sets were prepared, extracted, and quantified alongside the test indicators. Shipping and laboratory controls were maintained inside 50-ml conical tubes with 0.2µm filter caps that allowed equilibration of temperature and humidity exactly the same as the test samples. The controls were maintained at ambient temperature and humidity and measured using data loggers, while test samples were exposed to HHA.

HHA Inactivation of Enveloped RNA Virus
In order to develop RSM performance envelopes, partial inactivation data was needed especially at the DOE mid-point. A small test at 50, 60, 70, and 80 • C, RH at 5 or 90% for 1 and 24 h, produced critical data showing that high RH was critical for enveloped virus inactivation ( Table 1). The virus inoculum was held at RT while test samples were exposed to temp/RH. There was only 2.4 log 10 enveloped virus inactivation at 70 • C, 5% RH, for 24 h compared to complete inactivation at 70 • C, 90% RH, TABLE 1 | Log 10 PFU recovered (survival) and reduction (L.R.) after drying 6 on APC coupons (enveloped virus test indicator) and then hot air exposure at 5 or 90% RH for 1 or 24 h.

Parameters
Survival L.R. Parameters Survival L.R. for 1 h. In addition, there was no enveloped virus inactivation after drying enveloped virus on NTC and treating at 26.7 • C at 5, 30, 55, and 80% RH for up to 10 d (data not shown). These data drove the test parameters shown in Figure 1. Although 65 • C was greater than the original 60 • C temperature goal, 65 • C was selected for the experimental design because inactivation at <60 • C was too slow to generate RSM that met the <24 h objective. Therefore, 65 • C was incorporated in the DOE and 60 • C, 70% RH, 5 h, was set as the mid-point for the final DOE in order to obtain sufficient partial kill data to develop the RSM performance envelopes.
Overall, there were several thousand trial-and-error data points leading to the HHA decontamination DOE data set in Table 2, which shows data from nineteen test runs with five independent virus preparations per coupon type. There were 475 6 virus test samples and 475 corresponding controls, plus non-inoculated coupon controls. Each test variable (temperature, RH, and time) had an impact on decontamination. Importantly, high RH, high temperature, and longer times increased inactivation rates as expected. Inactivation of enveloped virus dried on to materials was found to be ≥7 log 10 at 60 • C, 90%, for 9 h. This surpassed the original goal of ≥7 log 10 at ≤60 • C, 24 h. Furthermore, the enveloped virus inoculum held in aqueous solution (solution controls) was shown to be much more sensitive to the rise in temperature than virus dried on to materials and subjected to hot, humid air. Therefore, water and/or high humidity were significant parameters for inactivating enveloped virus.
There was one statistical outlier for APC at 60 • C, 50% RH for 5 h ( Table 2). The subsequent RSM modeling showed that decontamination kinetics of all materials was essentially the same so long as this data point was treated as an outlier. A full set of tests was not repeated for 60 • C, 50% RH, for 5 h because there was no information advantage to update this outlier for this large data set. Much trial and error indicated that APC coupons required 70% ethanol cleaning to remove residue that interfered with microbe recovery and ethanol cleaning was specific to APC; hence, the source of this outlier was likely due to insufficiently cleaned APC coupons.

HHA RSM for Enveloped RNA Virus Inactivation
RSM for each material (Figure 4) were constructed from the DOE 6 inactivation data in Table 2. Inactivation of 6 virus was similar on all coupon types indicating that material type had little impact on inactivation kinetics. Since virus inactivation was similar across each material, an additional composite RSM was generated (Figure 4). The models more clearly showed virus inactivation kinetics compared to the table data. At 1 h, there was a 90% probability of at least a 7 log 10 virus inactivation at 65 • C, 90% RH, after virus-contaminated materials were treated with HHA. The composite model after 1 h of virus inactivation was critical for directing HHA parameters during the aircraft field testing. The RSM models were calculated using relative humidity because the laboratory environmental chambers controlled humidity using relative humidity. After extensive testing to develop methods and results to develop the RSM models using relative humidity, the sponsoring agencies then requested that relative humidity be converted to dew points for large aircraft testing. Based on a 4 log 10 objective and the 4 log 10 inactivation parameters in the RSM composite model (Figure 4), relative humidity was converted to dew points. This conversion showed that 4 log 10 inactivation occurred at a dew point humidity of ≥59.4 • C for 1 h, excluding ramp-up and ramp-down times, so long as the temperature was above the dew point. Therefore, the translation made it apparent that absolute humidity or dew point strongly correlated with inactivation. For example, there was a 90% probability of a 4 log 10 kill at ≥61.7 • C and a dew point humidity ≥59.4 • C. The temperature could go up or down with the same inactivation kinetics so long as the absolute humidity/dew point was maintained. Hence, the inactivation kinetics correlated with the total amount of water in the air. Models with dew points and/or absolute humidity are in the queue for future development.

C-130 Aircraft Field Test
A total of 100 enveloped virus test indicators were distributed at 20 sites throughout a C-130 fuselage as depicted in Figure 2. Based on the RSM in Figure 4 and with a >4 log 10 inactivation goal, the HHA goals for the C-130 were ≥61.7 • C and a dew point humidity of ≥59.4 • C for 1 h, excluding ramp-up and ramp-down times. The measured HHA parameters for the C-130 are shown in Figure 5. For clarity, both RH and dew points are shown in Figure 5 so that the parameters can be compared to both the RSM expressed in RH and the large-scale aircraft objective expressed as a dew point. As shown, the temperature and humidity goals were exceeded. In addition, the ramp-up time was 2.5-3 h, which was longer than lab testing with a ramp-up of 1 h; hence, ramp-up conditions on the aircraft likely contributed to enveloped virus inactivation. Ramp-down time was nearly 2 h, which was also much longer than lab conditions. The total HHA time including ramp-up and ramp-down for the C-130 was 5.5 h compared to 2 h in the lab-based RSM. The C-130 temperatures were always above the dew points so there was no condensation on the virusinoculated APC samples (Figure 5). Samples were shipped to FIGURE 4 | Models with a 90% statistical probability of 6 virus inactivation after virus was dried on to nylon, polypropylene, wiring insulation, or aircraft performance coating after incubation in hot humid air. Composite models are the accumulation of data for all 4 materials. Areas are represented as >7 log 10 (red), 6-7 log 10 (orange), 5-6 log 10 (gold), 4-5 log 10 (yellow), 3-4 log 10 (green), 2-3 log 10 (dark green), 1-2 log 10 (blue), and <1 log 10 (dark blue) virus inactivation.
Frontiers in Bioengineering and Biotechnology | www.frontiersin.org Frontiers in Bioengineering and Biotechnology | www.frontiersin.org 9 October 2020 | Volume 8 | Article 592621  (100) 0.0 ± 0.0 8.1 ± 0.1 The detection range was 1.0-8.1 log 10 where the upper limit was the amount extracted from the controls. Since there was no enveloped virus recovery on any test indicator, there was practical confidence of a complete virus kill. The number of test indicators are in parentheses. Not applicable (NA).
a lab for extraction. Table 3 shows 8.1 ± 0.1 log 10 PFU were recovered on 10 non-shipped control samples and 25 groundshipped control samples out of 8.3 ± 0.1 log 10 PFU dried on to the APC coupons. The loss of 0.2 log 10 PFU on the control coupons could have been due to extraction efficiency and/or a small drop in virus viability. In either case, the 8.1 log 10 PFU recovery on the control coupons was high, and if validated there was no significant loss of viability during transportation to and from the field site. The detection range was 1.0-8.1 log 10 . Since the temperature and RH limits were exceeded for a 7 log 10 kill and there was no viable virus recovered from any test indicator, there is practical confidence that a complete >8 log 10 kill was achieved.

C-17 Aircraft Field Test
A total of 110 enveloped virus test indicators were distributed throughout the C-17 as depicted in Figure 3. Based on the RSM in Figure 4, results from the C-130 test, and with a >4 log 10 inactivation goal, the HHA parameters were set to ≥61.7 • C and a dew point humidity of ≥57.2 • C for the C-17, at which point a 1-h "hold time" was started. Another goal was to increase the temperature and humidity during the 1 h "hold time" to 61.7 • C and a dew point humidity of ≥59.4 • C. The measured HHA parameters for the C-17 are shown in Figure 6. For clarity, both RH and dew points are shown in Figure 6 so that the parameters can be compared to both the RSM expressed in RH and the large-scale aircraft objective expressed as a dew point. The temperature reached 65.0 • C and humidity reached a 60.6 • C dew point, surpassing the temperature and humidity goals during the 1 h hold time. In addition, the ramp-up time was 1.5 h, which saved >1 h compared to the C-130 test, but was 0.5 h longer than lab tests. Ramp-down time was 45 min compared to 2 h for the C-130 and only one min for lab tests. The total HHA time for the C-17 was <3.5 h compared to 5.5 h for the C-130 and 2 h for lab tests. The C-17 temperatures were always above the dew points so there was no condensation on the virus-inoculated APC samples (Figure 6). Samples were ground-shipped to a lab for extraction. Table 4 shows 7.6 ± 0.1 log 10 PFU on 10 nonshipped control samples and 25 ground-shipped control samples out of 8.2 ± 0.1 log 10 PFU dried on to the APC coupons. The loss of 0.6 log 10 PFU on the control coupons could have been due to extraction efficiency and/or a drop in virus viability. In either case, the 7.6 log 10 PFU recovery on the control coupons was high and it showed no practical significant loss of virus viability during transportation to and from the field site. The limit of detection was 1.0 log 10 . Since there was viable PFU from at least one of the 110 test indicators, all samples with no viable PFU were assigned 1.0 log 10 . Site 15 on the C-17 showed the lowest log reduction of 5.9 ± 0.6 log 10 PFU. This site had the most air movement. Although there were no TPP R tubes with temperature and RH readings at Site 15, we speculate that the rapid air movement slowed the penetration of moist air across the 0.2 um membrane of the TPP R tubes, hence the lowest level of inactivation. Every other C-17 location had at least a 6 log 10 reduction. Ramp-up/ramp-down times were longer than the lab tests used for RSM, and raising the heat and humidity during the 1 h hold time contributed to the virus inactivation. Hence the goal of >4 log 10 reduction was surpassed as predicted by RSM.

DISCUSSION
6 was selected as a BSL-1, enveloped RNA virus test indicator for both lab and field tests. 6 has been widely used as an enveloped virus surrogate (de Carvalgo et al., 2017). The structural similarity of 6 to many other enveloped viruses including coronaviruses suggest that the 6 structure should be similarly susceptible to general decontaminants. Furthermore, the structural molecules of the virus are produced by host cells with similar temperature sensitivity at around 40 • C, further suggesting that 6 should be similarly susceptible to general decontaminants.
For this study, virus was dried ≥24 h before exposing to the environmental test conditions because (1) the level of virus drying/level of free water was an important consideration for standardizing these tests and ≥24 h drying mitigated drying variability and (2) ≥24 h drying was done in attempt to mimic environments where free water is limited, especially for virus bound in debris. There was no 6 inactivation after the virus was dried on different surfaces for at least 24 h at RT followed by a 10-d exposure to 26.7 • C at 80% RH, and there was only 2.4 log 10 inactivation at 70 • C, 5% RH, for 24 h. This is a higher level of temperature and humidity stability than reported in other enveloped virus studies, but none of the articles we reviewed tested virus that was dried for at least 24 h. A review of articles describing enveloped virus stability to temperature and/or humidity, including coronavirus, showed that most stability tests were with wet virus or virus that was only dried for an hour. Sizun et al. (2000), Rabenau et al. (2005), and Chan et al. (2011) tested similar coronaviruses including dried SARS-CoV, CoV-229E, and CoV-OC-43 at different virus concentrations, virus purities, additives such as fetal calf serum, and fomites, but all the tests were at ambient environmental conditions of 21-25 • C, which contrasts with no virus inactivation at 26.7 • C, 80% RH, 10 d in this study. Complete virus inactivation ranged from 3 h up to 13 d across these reports. Numerous published reports for enveloped virus stability have tested wet virus preparations (e.g., Anonymous-Department of Homeland Security, 2020; Biryukov et al., 2020;Castaño et al., 2020;Harbourt et al., 2020;Matson et al., 2020;Ren et al., 2020) and also showed much higher sensitivity to temperature and RH than reported in the current study. While there are numerous differences in methods and results described across these studies, a key difference in the testing described herein was that virus was dried at least 24 h before exposing to the environmental test conditions and that may be a key variable for virus stability.
There is general agreement in the literature that higher temperature and humidity will increase the inactivation rates for enveloped virus. However, a direct side-by-side study would eliminate the impact of methods to define temperature and humidity profiles. This would also highlight physiochemical similarities/differences inherent to a specific virus or common to virus families (Nims and Plavsic, 2013;Castaño et al., 2020). Side-by-side testing comparing dried 6 to several dried coronaviruses, including 6 SARS-CoV-2, are moving forward in our lab. Once these data are collected, a more standardized comparison can be made between the temperature and humidity susceptibility among these enveloped viruses.
There are numerous mathematical tools that can be used to analyze decontamination data. RSM models were very useful here since they provided simplified graphical data to rapidly adjust to changing goals for log reduction, time, temperature, and humidity. This became very important since the HHA goals (requirements) were frequently redefined by end users (aircraft maintenance authorities at Aircraft Mobility Command) over the past 20 years, and numerous goal changes were directed during the COVID-19 pandemic. The original goals were to show enveloped virus inactivation of ≥7 log 10 out of a ≥8 log 10 challenge. This challenge level was set because measurements with high concentrations of microbes greatly increased the confidence in inactivation and it mitigated the risk of incomplete decontamination (Hamilton et al., 2013), an important consideration for decision makers that assess contamination/decontamination. High challenge levels are important since exposure limits (infectious dosages) are not well-defined for many viruses such as SARS-CoV and SARS-CoV-2. Although HHA does not fall under the EPA's regulatory jurisdiction because it is not classified as a chemical disinfectant, the inactivation goal was reduced from 7 log 10 to 3 log 10 inactivation by the funding agencies during the COVID-19 pandemic to match the EPA N-list for decontaminants and provide decontaminant options to end users. The HHA RSM showed 90% probability levels of inactivation for log reductions of 0 through 7 log 10 , which allowed for rapid adjustments during the aircraft testing.
There were originally no time limitations set for aircraft sterilization (spore decontamination). The HHA time goal was reduced to ≤24 h for enveloped virus to increase HHA throughput and user acceptance. The time goal was further reduced by the end users and funding agencies to "1 h hold time; <6 h total time" after the COVID-19 pandemic had begun.
The HHA temperature goal was originally capped at the accepted temperature limit of C-130 aircraft (80 • C), but that was for sterilization (spore decontamination). The HHA temperature goal was reduced to 60 • C for the enveloped virus to increase HHA acceptance over a range of aircraft since 60 • C is a standard, accepted test temperature for most aircraft materials. The time goal of ≤24 h provided the rationale for measuring HHA over several hours to develop RSM models. However, the end users and funding agencies raised the accepted temperature limit to 68.3 • C during COVID-19 to reduce decontamination time.
Although there was a strong preference to reduce/eliminate humidity from HHA to reduce cost and complexity, all of the experimental data here showed that 6 was stable after drying. Virus was quickly inactivated after it was diluted in aqueous solutions and heated (solution control data). For HHA, the combination of high heat and high humidity, i.e., a high dew point, was critical to virus inactivation on materials in hour-long timeframes. This is consistent with long-known observations that enveloped virus is stabilized after drying on fomites. As an example, the enveloped virus smallpox was stabilized after it peeled/sloughed along with host cells and then dried on porous materials including clothes and blankets. Furthermore, historic data shows that enveloped viruses were stabilized in cold, dry climates (Fenn, 2001).
The end goal of the RSM was simplistic graphs where all decontamination test variables could be compared, understood, and utilized for field testing and field application. This was particularly useful for large hot, humid air decontamination systems with fluctuating temperature and RH, potential power disruptions, and decontamination goals that have fluctuated widely over time. Since the RSM models showed little difference in inactivation kinetics among the different materials, a composite RSM was developed as a primary reference guide. This reference guide served the purpose of guiding HHA field testing in aircraft. The RSM models translated to field demonstrations because there was a ≥5.9 log 10 enveloped virus inactivation after both aircraft field tests, as predicted. The RSM models were also used as a guide for a number of iterative improvements between the first and second aircraft demonstration that shortened the overall HHA decontamination time from 5.5 to <3.5 h.
There were numerous conclusions. The 6 test indicators were useful measurements of inactivation for both lab and field testing. HHA successfully inactivated enveloped virus at high temperature and high humidity. The RSM models were adaptable to changes in user requirements and successfully utilized to predict HHA results for field applications. Much methods development and testing remains to confirm HHA inactivation kinetics across virus species, and to develop HHA models that include HHA ramp-up and ramp-down times. Lastly, HHA decontamination of aircraft was successful because the enveloped virus on the C-17 aircraft was inactivated, the aircraft was flown after the HHA field test, and then the C-17 was certified as an operational aircraft.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/supplementary material.

AUTHOR CONTRIBUTIONS
TB was the principal investigator and the principal for methods development, test design, execution analysis, and writing. AY was a lead tester and contributed to methods development test design, execution, and data assembly. EB-K, NK, ZM, MBo, and EO participated in methods development, test design, test execution, and data assembly. SH, MK, MH, CM, RM, JI, MBe, BG, VY, SL, EH, and BL contributed to test execution and data assembly. HB was the math modeler. All team members participated or contributed to writing, figures, tables, and/or editing.