III-1. Basic concepts in (indoor and outdoor) airborne transmission.
Since the dawn of humanity, mankind has suffered from infectious diseases due to a variety of pathogens. In the recent decades, epidemiology has focused more on non-transmissible illnesses, such as heart disease, cancer, or obesity. However, the recent COVID-19 pandemic reminds us that the burden of infectious illnesses, especially respiratory ones, has not been eliminated.
The concepts that are discussed below are partly restricted to respiratory disease in the case where the responsible pathogen is inhaled by a susceptible person. It concerns diseases which affect organs located in the respiratory tract. As discussed in [
12] and references therein, in the case of COVID-19 the route of transmission has been a matter of intense debate but, as stated in the introduction, it is nowadays largely recognized that the major transmission path is through airborne exchange i.e., by inhalation of an aerosol that has been exhaled by an infected person [7-11].
In matters of infectious disease and epidemiology, a key problem is to assess the dose-response relationship i.e., what is the probability of infection resulting from exposure to a pathogen dose i.e. a given level of exposure within a given time [
15]. A dose-response function
relates the dose X to the probability of infection. It is clear that
must be a monotonically increasing function of the dose, starting from zero at zero dose and increasing toward an asymptote
for large values of
. There are several probability laws that can be used for
as discussed in [
15], one of the most widely used being the exponential form:
where Π is a numerical factor which depends on the definition of the dose and its counting unit. In fact, one of the recognized difficulties in the dose-response model is first to define the dose. It is beyond the scope of the present paper to examine this question in detail and the reader is referred to the book of Haas [
16] and to Brouwer et al. [
15]. For an airborne disease,
Wells [
13], using this exponential law,
defined a quantum of contagium as a hypothetical quantity of pathogens that has been inhaled, per susceptible individual, when 63.2% (corresponding to
of these individuals display symptoms of infection. It is linked to a probability of infection which then follows a Poisson law:
The quantum is dimensionless
but is a counting unit (as dozens versus unity, or moles compared to molecules) which is clearly linked to the choice of
in eq. (1). Of course, and as discussed in Brouwer et al. [
15] and Rowe et al. [
6], its value, in terms of the number of pathogens, depends on a variety of mechanisms: inhalation of airborne particles, pathogen inhibition by host defenses or losses by some other processes, before any replication will begin in an infected cell. Obviously, a quantum corresponds statistically to a number of pathogens much greater than one. Hereafter, all discussions will be carried out in terms of quantum concentration; any other choice would not alter the fundamentals of the reasoning.
Considering the concentration of quanta in space (in m
-3 units),
, the inhaled dose
X during a time of exposure
t, can be written as:
being the pulmonary ventilation rate (taken as 0.5 m
3/h in the present investigation). Note that this definition of the dose does not require a homogeneous distribution of quanta in space. Only
at the inhaled location (mouth and nostrils) has to be considered. Note also that due to the extremely low concentration of quanta in air,
is not really a continuous function of
(since a number of viruses is of course an integer) but can be treated as such due to the stochastic character of the problem and the search for a statistical solution. When, the quantum concentration can be considered has constant during the time of exposure, then expression of dose X simplifies as:
In the case of an indoor room or building with well-mixed air, it is possible to write a conservation equation for the quanta, which, together with eq. (3), leads directly to the stationary state dose value and to the well-known Wells-Riley probability [
14]:
where
is the probability of infection for a susceptible person,
a quantum production rate per infector per unit time,
the pulmonary ventilation rate,
the ventilation rate,
the number of infectors, and
the time of exposure
. In the outdoor situation, there are very few models available and below we will present an extension of the model that we recently developed [
12] for the determination of
.
III-2. Box model of outdoor transmission
As discussed in the previous section, the first step in any model, whether indoor or outdoor, is to evaluate the concentration of the virions in inhaled air. The outdoor model that we developed previously [
12] is essentially a “box” model as described in chapter 5 of "Environmental Impact Assessment" [
36], and developed previously by numerous researchers [
37,
38,
39]. Box models are based on mass balance equations and are the simplest atmospheric models that can be used to evaluate the mean concentration of pollutants (molecules or particles), downwind of a source.
Our 2021 model considered mono-sized infectious microdroplets and their airborne behavior. In the following paragraphs, we develop, in some detail, an identical outdoor box model, using the Wells notion of quantum for the counting of virions. We consider an outdoor volume (atmospheric box) as illustrated in
Figure 1, with the wind blowing along the x axis and, as in our previous work; we suppose that there are no quanta escaping the volume above a height
H along z. The evaluation of
H is the most critical part of the model. It is also assumed that the quantum density does not change across the wind:
and hence the quantum concentration is considered only as a function of
along the wind. Although at low values of
,
is a function of height
, assuming that eq. (6) holds everywhere, the height dependency does not change the concentration balance between what is produced in the bulk of the box and what emerges at its downwind border at large
x, where everything has been mixed by the turbulent dispersion. This assumption is inherent to box models [
39]. It can, therefore, safely be concluded that eq. (6) has no influence on the quantum concentration at this border.
Then, assuming stationary state i.e.,
a conservation equation for the quanta can be written as:
where
is a density of infectors per unit surface (assumed homogeneous and therefore constant),
is the wind velocity and
is the virus lifetime defined from the temporal exponential decay of active virions in microparticles, due to natural physicochemical and photochemical processes (see section V for a detailed discussion). Note that with this definition,
is slightly different from the so-called half-life which is the time required to decrease the active virion concentration by a factor of two (since here
corresponds to
).
Note also that the infectors are located at the bottom of the atmospheric box (which can include houses as we will discuss in section IV-2) but this has no influence on the calculation since we assumed an homogeneous dispersion of the viral aerosol in the vertical dimension of the box, as discussed in a previous paragraph.
With
the quantum concentration at
, we can derive the following value for the quantum concentration as a function of distance
x:
In an area where there is no infector,
, eq. (8) leads to a simple downwind exponential decay of the quantum concentration:
On the other hand, and for a virus lifetime much longer than the hydrodynamic time
(i.e.,
, eq. (8) leads to the following value for
:
This is analogous to the equation derived in Rowe et al. [
12] for
, and which expresses the conservation of quanta in the atmospheric box shown in
Figure 1 when there is no decay due to viral inactivation.
Knowing and using equation (2) and (3), it is then possible to calculate the probability of infection at any distance . Note that it is also possible to use a multi-box model (like that discussed further in section IV.2) to take into account more complex situations regarding a non-uniform repartition of infectors.
The areal density of infectors can be taken as:
with
the proportion of infectors and
the areal density of population in the space. The condition
requires that there is no gradient in mean infector density across the wind. If we assume a value of
for the width of the source, then, for
eq. (10) also reads as:
where
is now the total population in the area
In equations (8), (10) and (12),
a key parameter is the value of and therefore , as discussed at length by Rowe et al. [
12] in their supplementary materials. For strong winds i.e
. > 6m/s at 10 m height and at night time or low solar insolation, the atmosphere can be considered as neutral in the so-called Pasquill– Gifford– Turner classification [
40,
41,
42], which means there is no tendency for enhancement of air turbulence (instability) or it being suppressed (stability) through buoyancy effects. In fact, it is known that airborne pollutants emitted locally are transported and dispersed within the so-called atmospheric boundary layer (ABL: the tropospheric bottom layer), whose thickness is usually lower than one thousand meters [
43], except for strongly unstable atmospheres. At any distance from the source, an order of magnitude value of
H versus
x can be estimated by the vertical dispersion length used for Gaussian plumes and shown in
Figure 2 [
44,
45].
The wind itself depends on the altitude but its variations above ten meters are rather small within the ABL [
46]. In the next section therefore, it will be assumed to be independent of altitude and taken as the ten-meter value.
III-3. Possible airborne epidemic triggering by long-range transmission
Let us now consider two strongly populated areas, designated as the “source” and the “target” respectively, separated by an unpopulated area (no man’s land). The source population is considered as infected whereas no sick people are initially present in the target area (hereafter box 3). At the downwind border of the source (hereafter box 1) characterized by a length
L1, it is possible to quantify the quantum concentration
n1,q(L1) following eq. (8) taking into account that at the upstream border of the source
n1,q(x1=0) =
0. The downwind border of the source coincides as well with the upwind border of the "no man's land" section (hereafter box 2). However, as explained further in section IV-2, the dispersive height
H in box 2 (
H2) is higher than in box 1 (
H1). This impacts the initial quantum concentration
n2,q(x2=0) by a factor
H1/H2 at the upstream border of box 2 such as:
Then, the quantum concentration evolution n2,q(x2) in box 2 is ruled only by the virus lifetime according to eq. (9) since in box 2. This leads to a new value n2,q(L2) at the downstream border L2 of box 2. Again, this border coincides with the upstream border of the target area, box 3. At this interface, however, the dispersive height is not modified compared to box 2 (H3= H2) since H2 is already taken as an upper limit of the ABL thickness (see IV-2). Then, the quantum concentration entering box 3 is n3,q(x3=0) = n2,q(L2). The quantum concentration n3,q(x3) can be considered as spatially and temporally constant within box 3 provided that:
the pathogen lifetime is clearly larger than the hydrodynamic time within the target depth which is typically around 10-20 km.
the width of the target is smaller than the width of the source.
the emission source rate and meteorology do not change significantly during the time of exposure.
From this, the calculation of a probability of infection in the target area of population (assumed healthy and therefore susceptible) can be obtained combining eq. (2) and eq. (4).
It follows that the statistical number of contaminated susceptible people
is:
As exemplified below in section IV-2, the value of
will most often be extremely small, which shows
a quasi-zero risk at the individual level as already discussed by Rowe et al. [
12]. However, when the target is composed of a very high number
of individuals, then
a few people (≥ 1) could be infected. Of course, this process alone cannot sustain an epidemic but creates a few infectors (“primary cases”) which may trigger it.
In the next section we examine a hypothetical case study of the creation of COVID-19 primary cases in Northern France from Southern England in wintertime when the strongest winds are most often from west to east in Western Europe.