key: cord-0785187-yo11e55y
authors: Aydin, Mehmet; Evrendilek, Fatih; Aydin, Ismail Erkan; Savas, Seckin Aydin; Evrendilek, Deniz Eren
title: Transport dynamics of SARS-CoV-2 under outdoor conditions
date: 2022-04-02
journal: Air Qual Atmos Health
DOI: 10.1007/s11869-022-01196-x
sha: 22354d9c894a34120e775873fb9f09f205c19327
doc_id: 785187
cord_uid: yo11e55y

This study aimed at estimating the transport dynamics of a single severe acute respiratory syndrome corona-virus 2 (SARS-CoV-2)-laden droplet of 1 to 500 μm in diameter at a wind speed from 1 to 4 m/s. Motion dynamics of SARS-CoV-2-laden respiratory droplets under calm or turbulent air conditions were quantified using a combined model. Dalton’s law was implemented to estimate their evaporation. One-factor-at-a-time procedure was applied for the sensitivity analysis of model of deposition velocity. The transport distance of the single virus ranged from 167 to 1120 m as a function of the droplet size, wind speed, and falling time. The evaporation times of the droplets ≤ 3 and ≤ 14 μm in diameter were shorter than their settling times from 1.7 m in height at midnight and midday, respectively. Such droplets remained in the air for about 5 min as the droplet nuclei with SARS-CoV-2. The minimum transport distance of the respiratory droplets of 1–15 μm varied between 8.99 and 142 m at a wind speed range of 1–4 m/s, based on their deposition velocity. With their short transport distance, the larger droplet (30 to 500 μm) was not suspended in the air even under the windy conditions. The deposition velocity was found most sensitive to the droplet diameter. The droplets < 15 μm in diameter completely evaporated at midday and the droplet nuclei with the single virus can travel a minimum distance of 500 m under a horizontal wind speed of 3 m/s.

The globally rapid spread of severe acute respiratory syndrome corona-virus 2 (SARS-CoV-2) and its predictability are of increasingly great concern in terms of devising preventive and mitigative measures (e.g., social distancing) (Sarhan et al. 2021 ). The first deterministic epidemic model called Susceptible-Infected-Recovered (SIR) was introduced in 1927 and consisted of ordinary differential equations. These deterministic non-linear models of spread and deposition dynamics generate outputs based on parameters and initial condition values (Kermack and McKendrick 1927; Das et al. 2020; Aydin et al. 2021) . The motion dynamics of SARS-CoV-2 were quantified in indoor spaces (Anchordoqui and Chudnovsky 2020; Aydin et al. 2021; Sarhan et al. 2021) , in the calm air (Aydin et al. 2020a) , and under varying outdoor conditions (Das et al. 2020) .

The algorithms to solve the spatiotemporal (3D space + time) spread of SARS-CoV-2 can be non-deterministic polynomial-time (NP)-hard in a dynamic environment (e.g., changes in wind speed and direction, solar radiation, temperature, humidity, and population density) according to computational complexity theory (Aydin et al. 2021) . The nature of its flow (laminar or turbulent motion) depends upon a dimensionless quantity, which is called Reynolds's number (McPherson 1993; Sarhan et al. 2021) . In the presence of air resistance, such light objects may perform Brownian motion and follow a pattern of turbulent convective flows, depending on the particle size. In other words, the velocity of the particles at a given point may behave chaotically with time under the turbulent fluctuations (Anchordoqui and Chudnovsky 2020; Chong et al. 2020) . Pei et al. (2021) reported that buoyancy-driven flow regime led to a longer transmission distance and elevated exposure to viral aerosols than did the mixing airflow. According to Bond et al. (2021) , particle size drives its transport, deposition onto surfaces, and elimination by mitigation measures.

Although recent advances in aerosol instrumentation have enabled numerous deposition measurements (Emerson et al. 2020) , the airborne dynamics of SARS-CoV-2 are challenging to model under outdoor conditions. Therefore, a simple approach still remains to be desired to quantify the transport dynamics of SARS-CoV-2 and devise strategies for alleviating the spread of its respiratory droplets. The objective of this study was to quantify the settling and evaporation durations, dry deposition velocity, 1 and minimum transport distance of single SARS-CoV-2laden respiratory droplets as a function of wind speed (1-4 m/s), vapor pressure deficit (at 22 °C air temperature and 40% relative humidity), and droplet diameter (1-15 and 30-500 μm). Aydin et al. (2021) developed a simple model to estimate the terminal/constant velocity of SARS-CoV-2 under the calm indoor or outdoor air conditions. Although the velocity of a falling object is dependent upon the falling time, the predictions of the physical models show a wide variation for the falling time. The velocity of a micro-and/or nano-scale particle may not be constant due to the different particle size, vertical distance, and possible convection of the air. The velocity may be approximated by holding the time constant as follows (Aydin et al. 2021): where v sm is the effective speed of motion (m/s); m is the particle mass (kg); g is the gravitational acceleration (9.81 m/s 2 ); is 3.1416; r is the particle radius (m); is a

shape-dependent coefficient such as 0.47 for water droplets of a spherical shape which was assumed as 1.99 for a single isolated SARS-CoV-2 due to its spikes; a is the air density (1.2041 kg/m 3 at 20 °C and 1 atm); is the dynamical viscosity of air [1.85 × 10 −5 kg/(m s)]; and e is a mathematical constant (the base of natural logarithm also known as Euler's number: ~ 2.718). The SARS-CoV-2 was assumed as a single particle of a spherical shape with an average diameter of 100 nanometer (nm) and a protein density of 1.35 × 10 3 kg/ m 3 . Thus, the volume and mass of a single virus were estimated at 5.236 × 10 −22 m 3 and 7.07 × 10 −19 kg, respectively (Aydin et al. 2020a (Aydin et al. , 2021 . The cross-surface area of a single virus is assumed to equal the surface area of a hemisphere 2 r 2 where the base is not included, due to the wing area of its spikes. Thus, the expression of √ 2 ∕ a in Eq.

(1) can be parameterized as √ ∕ a for a single SARS-CoV-2. The falling (settling) time (t) of a single SARS-CoV-2 from a given height (z) in the calm indoor air whether it is encapsulated by a droplet or not can be quantified thus (Aydin et al. 2021 ):

Turbulence necessitates detailed physical analysis since the interactions within turbulence create a complex (chaotic) flow pattern that changes with time and space. Under the windy conditions, the interactions among SARS-CoV-2, respiratory droplets, and dust particles may occur. The transport distance of the atmospheric particles depends on flow conditions and particle size. Shao (2008) established an approximate relationship between particle size and transport distance under specific atmospheric conditions: where x d is particle travel distance (m); W s is wind speed (m/s); z is height from the ground surface (m); and v d is dry deposition velocity (m/s). Though generally calculated as a function of particle flux and concentration, the deposition velocity can be obtained as follows (Giardina and Buffa 2018) :

where r z is the total resistance (s/m) to transport as a function of particle diameter (d) and height (Z) (Giardina and Buffa 2018; Aydin et al. 2021 ). In the above equation, we

used v sm as the settling velocity (m/s) of the particle. The total resistance r z can be determined as follows (Giardina and Buffa 2018) :

where r a is an aerodynamic resistance (s/m); and r ql is a quasi-laminar sub-layer (boundary layer) resistance (s/m). The phenomena associated with the boundary layer effects are influenced by Reynolds's number and surface roughness. Re-entrainment can be analyzed considering the drag and frictional forces on the particles on or very close to the solid surfaces (McPherson 1993) . The aerodynamic resistance can be determined using the Monin-Obukhov similarity theory (Giardina and Buffa 2018) :

where z 0 is the surface roughness height above the displacement plane; k is the von Karman constant (0.4); and u * is the friction velocity, the intensity of the atmospheric turbulence. The parameter in Eq. (6) can be calculated according to Brandt et al. (2002) :

and where L is Monin-Obukhov length which was set as 125 m and − 150 m for stable and unstable air conditions, respectively (Sathe et al. 2010 ). To facilitate a comparison, we have selected the same reference height (z) of 1.7 m, an average human height. The boundary-layer resistance can be estimated thus (UK Met Office Hadley Centre 2017):

where S c is the Schmidt number (diffusion vs viscosity); and P r is the Prandtl number of 0.72 for the lower atmosphere (Liu et al. 1979 ; UK Met Office Hadley Centre 2017). At 20 °C in the ambient atmosphere, the Schmidt number was assumed as 1.28 (Wu et al. 1992) . u * was assumed as 0.114 m/s for smooth grounds in Eqs. (6) and (9) (Giardina and Buffa 2018) . The assumed relationship in Eq. (9) is an over-simplification. The deposition speed was modeled for the stable and unstable air conditions, based on a single SARS-CoV-2 particle (protein) density of 1.35 × 10 3 kg/ m 3 , a respiratory water droplet density of 998 kg/m 3 , and a roughness length z 0 of 0.02 m (Giardina and Buffa 2018; Aydin et al. 2020a ).

(5) r z = r a + r ql

The droplets are subject to evaporation once coming out of the respiratory tract. The further the droplets are away from the mouth, the smaller the droplet size is due to evaporation (Anchordoqui and Chudnovsky 2020; Stariolo 2020). Aydin et al. (2021) approximated droplet evaporation as a function of air temperature and relative humidity for indoor conditions thus:

where E di is indoor droplet evaporation (mm/day); T a is daily mean value of air temperature in the range of 2-49 °C; and RH is relative humidity in the range of 10-90%. Solar radiation, wind speed, temperature, and relative humidity are the primary climatic drivers of the outdoor evaporation (Allen et al. 1994; Aydin 2008; Aydin et al. 2020b ). If pan-evaporation data are not available, then outside water evaporation can be computed via the simplest and most common type of Dalton's law as follows (Helfrich et al. 1982) :

where E po is outdoor pan evaporation (mm/day); e s is saturated water vapor pressure (mm of mercury); e a is actual vapor pressure in the air (mm of mercury); and u is a function of wind speed. Saturated vapor pressure (mbar) can be expressed thus (1 mbar = 0.75 mm of mercury) (Allen et al. 1994 ):

e a (mbar) can be estimated as follows (Allen et al. 1994 ):

u can be estimated thus (Helfrich et al. 1982) :

where W s is wind speed (m/s) at a height of 2.0 m (hereafter u 2 ); a and b are constants. To estimate the parameter of u , the daily constants ( a = 0.675 , b = 0.142 , with u 2 expressed in m/s) were derived from the data by Malek (1994) and Srivastava and Jain (2017) . In this study, droplet evaporation was adjusted by considering the effects of the surface area exposed to the atmospheric evaporative demand but by neglecting the collective influence of a typical cloud of droplets, based on Aydin et al. (2021) . Therefore, outdoor droplet evaporation (E do , mm/day) can be derived thus: 

Sensitivity analysis conducted in this study helped to identify the important model parameters, test the model conceptualization, and improve the model structure (Sieber and Uhlenbrook 2005; Aydin and Kececioglu 2010) . Therefore, the sensitivity of the dry deposition velocity v d in Eq. (4) to the droplet size ( d or 2r), the friction velocity (u * ) , Schmidt number S c , and Monin-Obukhov length (L) was determined (see Eq. 1 and from Eqs. 5, 6, 7, 8 and 9). In so doing, one-factor-at-a-time procedure was applied through which a given predictor was changed by a predetermined fractional increase and decrease in several steps, while keeping the others at their nominal values. Thus, their sensitivities were detected by monitoring the corresponding changes in the output relative to the baselines.

The transport distances of an isolated single SARS-CoV-2 are presented as a function of wind speed in Fig. 1 . Its dry deposition velocity was firstly derived from Eq. (4) as deposition is a process by which atmospheric particles are removed by a turbulent transfer and fall onto the earth's surface. Equation (4) may yield irrational results for the airborne nano-scale particles with a high density (such as a protein density of 1.35 × 10 3 kg/m 3 ), while Eq. (3) underestimates the transport distances considerably. Therefore, the transport distance of a single SARS-CoV-2 was estimated by using both v sm and v d for a comparison and ranged from 281 and 167 m at a wind speed of 1 m/s to 1120 and 667 m at a wind speed of 4 m/s, respectively (Fig. 1) . The atmospheric lifetime of the virus was not considered in Fig. 1 .

As a one-way flux, the deposition velocity of a single SARS-CoV-2-laden respiratory droplet is shown in Fig. 2 . Its velocity exhibited an asymptotic exponential growth with its rising droplet diameter. The dry deposition velocity of the respiratory droplet from 1.7 m in height varied between 0.05 m/s for a diameter of 1 μm and 0.19 m/s for a diameter of 15 μm. Evaporation time of a single SARS-CoV-2-laden respiratory droplet is illustrated in Fig. 3 , together with its settling time. On a diurnal average, the falling time of a respiratory droplet > 7 μm from 1.7 m was shorter than its evaporation time. In other words, the droplets larger than 7 μm settled faster than they evaporated, whereas the smaller droplets evaporated faster than they settled down. Similar findings were reported by Bourouiba (2020) and Wells (1934) . However, Anchordoqui and Chudnovsky (2020) reported very short evaporation times for the wide range of droplet sizes of 50 to 2000 µm. According to Aydin et al. (2021) , the atmospheric evaporative demand, aka vapor pressure deficit, increased as the day progressed and peaked ten times faster at midday than midnight based on hourly percentages of daily evaporation data. In this study, the approach suggested by Aydin et al. (2021) was implemented to calculate hourly percentages of daily evaporation. According to their approach, the droplets with the diameters ≤ 3 and ≤ 14 μm were capable of the complete evaporation and remained in the air for about 4.7 min as the droplet nuclei with SARS-CoV-2 at midnight and midday, respectively. The detailed theoretical discussions of the evaporation of aqueous droplets in a rapidly changing outdoor condition were made by Abramzon and Sirignano (1989) . Our simulations were independent of the host-dependent survival duration of SARS-CoV-2 suspended in the air and the possible influence of a typical cloud of droplets. According to Fig. 4 , the minimum transport distance of a single SARS-CoV-2-laden respiratory droplet with a diameter of 1 μm from 1.7 m in height was estimated at 142 m under a wind speed of 4 m/s. However, its evaporation can cause the virus to have a longer transport distance and time in the air since it reduces the droplet diameter (Bathula et al. 2021) . As can be seen in Fig. 5 , the minimum transport distance of the larger droplets (30 to 500 μm) under a wind speed of 3 m/s was very short. Since the settling times of the larger droplets were shorter than their evaporation times, they were not suspended in the air even under the windy conditions. According to Sarhan et al. (2021) , most of the respiratory droplets similarly fell onto the floor due to gravity before travelling 1.5 m in an indoor environment where the airborne transmission of the virus was ignored. Feng et al. (2020) reported that high relative humidity expanded the droplet diameter which in turn increased the deposition fractions, while, with the complex wind and humidity conditions, micro-droplets were influenced by convection and transported from the human coughs/sneezes to other humans in less than 5 s. Likewise, at very low wind speeds, convection was reported to significantly induce the mixing at the reference height, while wind perturbations exerted a more surface stress than did the average wind (Bourassa et al. 1999 ). According to Kayalar et al. (2021) , temperature, wind speed, relative humidity, and precipitation significantly influence the vitality and transmission of biological constituents of aerosols characterized by the suitability and atmospheric lifetime/travel distance of suspended particles. Hadei et al. (2021) concluded that air pollution can exacerbate COVID-19-related incidences and mortality.

The dry deposition velocity was most sensitive to the droplet diameter. The droplet mass and diameter were assumed to be parallel. Friction velocity was rated as the second most effective variable, whereas the contributions of Monin-Obukhov length and Schmidt number to the response were negligible (Table 1 ). According to Aydin et al. (2015) , the sensitivity outcomes may not reflect the site-specific conditions.

Our striking findings are illustrated in Fig. 6 . The droplets ≥ 15 μm settled down at the end of their transport distances. The droplets < 15 μm in diameter completely evaporated at midday and traveled as far as the single virus went. In other words, the droplets < 15 μm evaporated completely within their transport distances, while the droplet nuclei with the single virus traveled at least 500 m under a horizontal wind speed of 3 m/s. These results may be of interest to the environmental and/or public health communities. 

The suggested approach in this study contributes to the better quantification of the transport dynamics of SARS-CoV-2. The transport distance of a single virus showed a wide range of 167 to 1120 m as a function of the droplet size, wind speed, and falling time. The dry deposition velocities of the respiratory droplets increased with their diameters. In other words, the deposition velocity was most sensitive to the droplet diameter and least sensitive to Monin-Obukhov length and Schmidt number. The droplets larger than 7 μm settled down faster than they evaporated based on a daily evaporation rate. The droplets ≤ 14 μm were able to evaporate completely and remained in the air for about 5 min as droplet nuclei with SARS-CoV-2 at midday. The nuclei of these droplets with the single virus traveled at least 500 m under a horizontal wind speed of 3 m/s. Model validations and further experimental observations still remain to be desired to measure and enhance the predictive power.

Author contribution Mehmet Aydin, Fatih Evrendilek, and Ismail Erkan Aydin carried out the conceptualization and design of the study as well as the drafting and editing of the manuscript. Seckin Aydin Savas and Deniz Eren Evrendilek carried out the acquisition, analysis, simulation, and interpretation of data. All the authors reviewed the manuscript.

Data availability Data will be made available on request. Correspondence and requests for materials should be addressed to Mehmet Aydin.

The authors declare no competing interests. 

Transport distance of a single SARS-CoV-2-laden respiratory droplet as a function of its diameter

Droplet vaporization model for spray combustion calculations

An update for the definition of reference evapotranspiration

A physicist view of COVID-19 airborne infection through convective airflow in indoor spaces

A model for evaporation and drainage investigations at ground of ordinary rainfed-areas

Sensitivity analysis of evaporation module of E-DiGOR model

Sensitivity of soil evaporation and reference evapotranspiration to climatic variables in South Korea

Falling dynamics of SARS-CoV-2 as a function of respiratory droplet size and human height

Modeling and computer simulation of Drainage, Evaporation, and Runoff (DEaR) from bare soils

A model for indoor motion dynamics of SARS-CoV-2 as a function of respiratory droplet size and evaporation

Survival of expiratory aerosols in a room: study using a bi-compartment and bi-component indoor air model

Quantifying proximity, confinement, and interventions in disease outbreaks: a decision support framework for air-transported pathogens

A flux parameterization including the effects of capillary waves and sea state

Turbulent gas clouds and respiratory pathogen emissions: potential implications for reducing transmission of COVID-19

Modelling transport and deposition of caesium and iodine from the Chernobyl accident using the DREAM model

Vortices as Brownian particles in turbulent flows

Transmission of airborne virus through sneezed and coughed droplets

Farmer DK (2020) Revisiting particle dry deposition and its role in radiative effect estimates

Influence of wind and relative humidity on the social distancing effectiveness to prevent COVID-19 airborne transmission: A numerical study

A new approach for modeling dry deposition velocity of particles

Effect of short-term exposure to air pollution on COVID-19 mortality and morbidity in Iranian cities

Evaluation of models for predicting evaporative water loss in cooling impoundments

Existence of SARS-CoV-2 RNA on ambient particulate matter samples: a nationwide study in Turkey

A contribution to the mathematical theory of epidemics

Bulk parameterization of air-sea exchanges of heat and water vapor including the molecular constraints at the interface

Calibration of the Penman wind function using the Bowen ratio energy balance method

The Aerodynamics, Sources and Control of Airborne Dust

Human exposure to respiratory aerosols in a ventilated room: Effects of ventilation condition, emission mode, and social distancing

COVID-19 aerodynamic evaluation of social distancing in indoor environments, a numerical study

Estimating the systematic errors in turbulence sensed by wind lidars

Physics and Modelling of Wind Erosion, 2 edn

Sensitivity analyses of a distributed catchment model to verify the model structure

Engineering Hydrology. McGraw-Hill Education

Dry deposition in UKCA -what is it? How does it work? UKCA Training Workshop -Cambridge

On air-borne infection: study II. droplets and droplet nuclei

Dry deposition of atmospheric contaminants: the relative importance of aerodynamic, boundary layer, and surface resistances