Introduction
Bioaerosol particles in the stratosphere are of interest from aerobiological and astrobiological perspectives. Some previous studies have focused on bioaerosol particle sampling (Bryan et al., Reference Bryan, Stewart and Granger2014, Reference Bryan, Christner and Guzik2019) with the aim of developing sampling methods to determine the extent of biosphere sampling of bioaerosol particles in the stratosphere. Although bioaerosol particles have been found in the stratosphere, the mechanism by which these particles reach tens of thousands of meters above sea level is not fully understood. Previous studies reported the injection of tropospheric air to the stratosphere (Mote et al., Reference Mote, Rosenlof and McIntyre1996) or vertical diffusion of the air in the stratosphere (Mote et al., Reference Mote, Dunkerton and McIntyre1998). Some meteorological mechanisms, such as quasi-periodic oscillation or Brewer–Dobson circulation were reported to transport mass or momentum upwards in the stratosphere (Plumb, Reference Plumb2002; Butchart, Reference Butchart2014). However, these previous studies on troposphere–stratosphere transport focused on the transport of gas, which requires much less momentum compared with bioaerosol particles. Another study suggested that chemical aerosol particles are distributed in the lower stratosphere because of volcanic eruptions, biomass burning and dust. However, these particles only have the size of accumulation mode (0.1–1.0 μm) at altitudes as low as 12 km. These previous studies cannot fully explain how microorganisms as large as 10 μm, which are believed to be too large to reach the stratospheric level, were sampled in the stratosphere at an altitude of 41 km (Harris et al., Reference Harris, Wickramasinghe and Lloyd2002). Some previous studies concluded that it is highly possible that large microorganism particles in stratosphere originate from space (Wainwright et al., Reference Wainwright, Alharbi and Wickramasinghe2006; Alshammari et al., Reference Alshammari, Wainwright and Alabri2011).
In this study, we simulated vertical particle motion in the stratosphere using parameters that have been suggested to represent the main driving forces. The results revealed the physical characteristics of bioaerosol particles that could potentially stay within the stratospheric level, or be lifted up to the stratospheric level.
Materials and methods
Virtual particles with various physical characteristics of bioaerosol particles were simulated based on equation (19) of Miki (Reference Miki2020) with the addition of the vertical wind velocity, W, and the neglect of the gravitational effect of the sun. In the simulation, the motion of a particle was assumed to be governed by gravity, buoyancy, friction from the air and photophoretic force, which as discussed by Wainwright et al. (Reference Wainwright, Alharbi and Wickramasinghe2006). The photophoretic force results from the temperature difference in the gas next to the particles that results from the difference in solar heating. For micron sized particles exposed to sunlight the photophoretic force can be 20% of the weight of the particle (e.g. Tehranian et al., Reference Tehranian, Giovane and Blum2001). The terminal velocity of a particle (V) is given by:
where D ae is the aerodynamic diameter of the particle; W is the vertical wind speed in the stratosphere; r is the radius of the particle; ρ p is the density of the particle, assumed to be the same as water (1000 kg m−3) based on the aerodynamic diameter; ρ f is the density of air; F Δα and F ΔT are the photophoretic forces, which are driven by the difference in the thermal accommodation coefficient inside a particle and the thermal gradient of the particle, respectively; C D is the drag coefficient; Re p is the Reynolds number of the particle; and η is the viscosity of air.
The viscosity of air (η) was derived from the Sutherland's formula as:
where T 0 is the temperature at sea level and T is the temperature at each altitude, which was calculated as equation (3) using an approximate formula used in the US standard atmosphere (United States Committee on Extension to the Standard Atmosphere, 1976) as first introduced in Miki (Reference Miki2020).
Gryazin and Beresnev (Reference Gryazin and Beresnev2011) reported that in the stratosphere, vertical wind blows at a monthly average speed of approximately 5.0 mm s−1 in upward and downward directions and an annual average speed of approximately 1.0 mm s−1 in upward and downward directions. In this study, wind speeds were set as 5.0 mm s−1 (5.0 × 10−3 m s−1) and 1.0 mm s−1 (±1.0 × 10−3 m s−1) in upward and downward directions, and 0 mm s−1 because the simulation duration significantly varied depending on the variable settings such as particle size or density.
Reynold's number (Re p) was derived using the following equation:
The density of air (ρ f) varied depending on the altitude of the stratosphere, and for each altitude, ρ f was calculated using equations (5) and (6):
where M is the molar mass of the air, R is the gas constant, P is the air pressure at each altitude and H is the scale height of the Earth's atmosphere (8432 m). The approximate formula for the drag coefficient (C D) changed depending on the amplitude of the Reynold's number of the particle (Re p) as shown in equations (7) and (8). These were first introduced by Oseen (Reference Oseen1910).
When Re p < 1:
and when Re p ≥ 1:
The two types of photophoretic forces, F ΔT and F Δα, were calculated using equations (9)–(14), which have been previously introduced (Rohatschek, Reference Rohatschek1996; Wurm and Krauss, Reference Wurm and Krauss2008). It was assumed that F ΔT works downwards and F Δαworks upward (Keith, Reference Keith2010).
where, $F_\ast$ is the maximum force under pressure $p_{\rm \ast }$, J is the asymmetric parameter, C T is the thermal conductivity, κ is the thermal creep parameter, c m is the average speed of the atmospheric molecule, α is the thermal accommodation coefficient and Δα is the difference in the thermal accommodation coefficient inside a particle.
When S is the solar radiation intensity at each altitude,
where S 0 is the solar radiation at wavelength (L) and k is a coefficient.
From equation (15), the relationship between the solar radiation above the atmosphere and that on the Earth's surface is described as follows:
where ΔS is the difference between S 0 and the solar radiation above the atmosphere S earth.
From equation (16), k is described as follows:
When E(h) is given as equation (18), the effective irradiance is derived as equation (19).
Equation (19) can be approximated using the piecewise quadrature as equation (20)
As the actual solar radiation data, ASTM G-173-03, the distributions of power obtained by the National Renewable Energy Laboratory (NREL) were used (https://www.nrel.gov/grid/solar-resource/spectra-am1.5.html 2022/10/27; Fig. 1a). As the radiation data are obtained discretely every 0.5 or 5.0 nm, the effective irradiance was approximately calculated using equation (21):
where L 1 = 280.0, L 2 = 3990.0, $L_1^{\prime} = 280.5$ and $L_2^{\prime} = 3995.0$. The solar irradiance at each altitude was derived as Fig. 1(b).
The asymmetric parameter (J) and thermal creep parameter (κ) were set at the values J = 0.5 and κ = 1.14 under the assumption that particle surfaces perfectly absorb light, as described in a previous study (Wurm and Krauss, Reference Wurm and Krauss2008). The thermal conductivity (C T) and thermal accommodation (α and Δα) depend on the characteristics of the particles and ambient air; these parameters of bioaerosol particles have not yet been studied. We assigned 0.3 and 0.9 values to α to account for low and high thermal accommodations, respectively, and 0 and 0.2 to Δα. The values of thermal conductivity (C T) were set as 0.1 and 10, respectively, to cover a broad range of high and low thermal conductivities. Lastly, in order to analyse bioaerosol particle motion in the stratosphere, the diameters of some representative bioaerosol particles, such as pollen, spores and microorganisms, were selected as aerodynamic diameters in the simulation. For example, an aerodynamic diameter of 10 μm can be applied for small pollen grains because of various primary pollen taxa, such as small Salix pollen. Additionally, the maximum size of possible microorganism found in the stratosphere was found to be 10 μm (Harris et al., Reference Harris, Wickramasinghe and Lloyd2002; Wainwright et al., Reference Wainwright, Alharbi and Wickramasinghe2006), but was 2 μm in another study (Bryan et al., Reference Bryan, Christner and Guzik2019). Spores were represented by a particle with an aerodynamic diameter of 1–3 μm based on a previous morphology study (Reponen et al., Reference Reponen, Grishpun and Conwell2001). Thus, 1, 2, 5 and 10 μm were substituted as the representative aerodynamic diameters in this study.
In summary, every combination of the parameters given below were simulated.
The initial altitude was set to 30 000 or 60 000 m, and the terminal velocity was adjusted every 3600 s based on the equations explained above at each altitude.
Results and discussion
The results of the particle vertical motion simulation calculated based on equation (1) showed that the relationship between Δα and solar radiation plays an important role in determining whether a particle in the stratosphere rises or falls when there is no vertical wind (e.g. Figs. 2–4(a) and (c)). When the particle diameter is 1 or 2 μm, the particle is uplifted and deposited when the wind blows upwards and downwards, respectively (Figs. 2 and 3). When the vertical wind speed is 0 mm s−1, the particle is deposited when Δα is 0 though the particle is uplifted when Δα is 0.2, although if the particle is lifted or deposited is determined by the initial altitude only when the particle diameter is 2 μm, α is 0.9, Δα is 0.2 and C T is 0.1 (Fig. 3(d)).
When the particle diameter is 5 μm, upward F Δα is strong enough to carry the particle upwards regardless of a downward 1 mm s−1 wind (Fig. 4). In addition, the simulation results revealed that thermal conductivity influences particle motion when the particle size is 5 μm and when its magnitude was changed by two digits (Fig. 4(f)). Additionally, the body temperature of the particle, a variable of thermal accommodation and the colour of the particle seems to be largely influenced by the body temperature since the colour determines the solar radiation absorption rate of the particle.
The small pollen-sized particles (10 μm particles) struggled to levitate in the stratosphere unless the upward wind was strong enough and the ratio of the difference in the thermal accommodation coefficient inside a particle to the thermal coefficient of the particle itself was large (Fig. 5), although spore-sized particles always levitated when Δα was not zero. Thus, the results show that the 10 μm microorganisms found in a previous study could be from the Earth, assuming that the particle characteristics satisfy the conditions. Previous research has suggested that the morphology of a particle may influence the photophoretic force (Redding et al., Reference Redding, Hill and Alexson2015). When the spore images available on the PAAA website (https://www.paaa.org/gallery/spores_eh/#) were analysed using ImageJ (Schneider et al., Reference Schneider, Rasband and Eliceiri2012), the circularity of each spore (Epicoocum, Curvularia and Spegazzinia) was found to be:
Epicoocum: 0.843; Curvularia: 0.452; Spegazzinia: 0.669
Here, the circularity was derived using equation (22):
Thus, because some spores are not morphologically point-symmetric, the actual photophoretic force (F Δα) on them could be larger than that in the simulation. Although previous research indicated that a particle with a diameter of 5 μm could possibly be uplifted into the stratosphere when only the uplifting wind is considered (Gryazin and Beresnev, Reference Gryazin and Beresnev2011), the simulation in the present study, in which the photophoretic force was considered, showed that even a particle with a radius of 5 μm can be lifted when there is a difference in the accommodation coefficient of a particle or weak vertical uplift wind. This rise of a particle into the stratosphere appears to be identical to that of a microorganism attached to a dust particle (Barberán et al., Reference Barberán, Ladau and Leff2015; Hu et al., Reference Hu, Murata and Fan2020). This phenomenon has also been addressed in previous research (Wainwright et al., Reference Wainwright, Alharbi and Wickramasinghe2006).
In this study, bioaerosol particles were assumed to have a density of 1000 kg m−3 regardless of the particle type. However, Amaranthus pollen are reported to have a density of approximately 14 000 kg m−3 (Sosnoskie et al., Reference Sosnoskie, Webster and Dales2009). Another previous study indicated that dry ragweed pollen has a density of 840 kg m−3, but 1280 kg m−3 in 100% humidity environments (Harrington and Metzger, Reference Harrington and Metzger1963). A previous study on spore particles assumed that a spore of Gymnosporangium juniperi-virginianae is 1200 kg m−3 (Fischer et al., Reference Fischer, Stolze-Rybczynski and Cui2010) while fungal spores of Lycoperdon perlatum have a density of 770 kg m−3 (Tesmer and Schnittler, Reference Tesmer and Schnittler2007). Concerning dust particles, a dust particle in the Asian dust-storm theoretically has a density of 2600 kg m−3 (Iwasaki et al., Reference Iwasaki, Minoura and Nagaya1983); however, because the chemical compounds vary significantly, the mass density will also vary significantly (Fergusson et al., Reference Fergusson, Forbes and Schroeder1986). Thus, understanding the motion of a specific type of particle in the stratosphere requires conversion of the geometric particle size to the aerodynamic particle size.
A previous study collected Bacillus spp. and Penicillium at an altitude of 20 km and found that cultured Penicillium had characteristic blue/green fruiting bodies with white mycelium (Smith et al., Reference Smith, Griffin and Schuerger2010). In order to understand the relationship between thermal accommodation and the uplifting force in the stratosphere, sampling the bioaerosol particles as they are and performing morphological analysis are important.
Although the results indicate that some types of bioaerosol particles can be lifted up to the high altitude (i.e. into the stratosphere), bioaerosol particles reported to be sampled in the stratosphere could have been contaminated. In addition, this study did not take into account global horizontal circulation of the stratospheric atmosphere and seasonal changes in the stratospheric environment. Previous research has suggested that the driving force by the electric field above thunderstorms can also be an uplifting force in the stratosphere (Dehel et al., Reference Dehel, Lorge and Dickinson2008). Thus, performing more in situ bioaerosol sampling in the stratosphere and more multiple bioaerosol transportation simulations considering stratospheric phenomenon are required for a full understanding of stratospheric airborne ecology.
Conclusions
The simulation of various bioaerosol particles in the stratosphere showed the significance of a difference in thermal accommodation in determining particle motion changes for a diameter size of 5 μm. This shows that microorganisms such as fungus or bacteria attached to dust can possibly levitate in the stratosphere. In addition, 10 μm bioaerosol particles can be levitated when the morphological or physical characteristics are ideal. In contrast to Wainwright et al. (Reference Wainwright, Alharbi and Wickramasinghe2006), these results suggest that 10 μm particles are not necessarily of extra-terrestrial origin and can be lifted from lower in the Earth’ atmosphere.
Conflict of interest
None.