FormalPara Highlights
  • Focusing on the period June 2020–February 2021, the paper estimates the probability that a randomly-chosen passenger uninfected by Covid-19 contracted the disease during an average (two-hour) domestic jet flight in the United States

  • The analysis uses numerous data sets to estimate the probability that a passenger boarding a US domestic flight over the observation period carried contagious Covid-19. That probability varied considerably over the nine months considered.

  • In calibrating its probabilistic model, the paper uses data from all available peer-reviewed papers that specified the seating location(s) of passengers who harbored Covid-19 when they boarded a flight, and also the seating locations of other passengers who did and did not contract the coronavirus in flight.

  • The point estimate for the probability of contracting Covid-19 on board an average domestic flight was about 1 in 2000 for the nine-month study period. That estimate is subject to sizable uncertainty.

  • The approach to modeling advanced in the paper can be adapted to make risk estimates in other circumstances, in connection with both the coronavirus and other infectious diseases.

1 Introduction

This paper estimates the probability that, early during the Covid-19 pandemic, a US domestic air traveler would contract the coronavirus while in flight. The focus is on the nine-month period from June 2020 to February 2021, a time frame that excludes both (i) the first months of the pandemic, when US air travel came practically to a halt, and (ii) the period starting in March 2021, when the use of vaccines was accelerating. There was no consensus on the magnitude of in-flight infection risk during those nine months. The nationally-renowned expert Dr. Anthony Fauci stated that he would not fly until the pandemic subsided, while the US Centers for Disease Control and Prevention (CDC) warned against air travel throughout the nine-month period. Yet United Airlines declared that Covid-19 transmission risk was “nearly nonexistent” even on its full flights, while Southwest Airlines characterized the risk as “virtually nonexistent.”

It is not easy to estimate the risk that an uninfected passenger on a US domestic flight would contract Covid-19. Here we attempt to do so with a probabilistic model which assumes that three things must occur for an in-flight infection to arise: there must be at least one contagious passenger on board, universal mask wearing must fail to prevent transmission, and the uninfected traveler must be seated fairly close to a contagious one. (As we will explain, we treat other means of on-board transmission as second-order effects.) The model estimates the joint probability of these three events.

The primary focus here is on economy-class passengers who took two-hour US domestic flights on either Boeing 737 or Airbus 320 jets (two hours is the average duration of domestic flights). We make estimates of Covid-19 infection risk on a month-by-month basis during the observation period, using the actual distribution of seat-occupancy levels in each month. For full flights, our approximate distribution for the probability that a randomly chosen uninfected passenger would contract Covid-19 in flight has a median of 1 in 2250 and a mean of 1 in 1450. The 10th percentile of the estimate is about 1 in 14,000, while the 90th percentile is about 1 in 600.

Below we start with a literature review, followed by development of the general risk model for Covid-19 transmission, then the estimation of model parameters, and then the presentation of risk estimates.

2 Some relevant literature

There is a voluminous literature about the transmission of viral infections in closed spaces. Berry et al. [1] summarized much of this work and assessed its general implications about limiting the spread of Covid-19. Noakes et al. [2] modelled the transmission of airborne infections in indoor environments, based on such factors as ventilation rate and room occupancy. Bazant and Bush [3] proposed upper limits on “cumulative exposure time” for uninfected people in an indoor space, given an occupant with contagious Covid-19 who exudes respiratory droplets. . Chen and Liao [4] offered probabilistic analyses about indoor transmission of various influenza viruses, as a function of viral kinetics, aerial transmission potential, and population dynamics . A meta-analysis performed by Chu et al. [5] considered 216 studies about numerous aspects of viral transmission.

Long before the Covid-19 pandemic, the transmission of infectious diseases during passenger flights was recognized as a serious health hazard. Evidence of in-flight infection arose in connection with SARS and the H1N1p flu pandemic [6, 7], and public health agencies including the World Health Organization and US Centers of Disease Control (CDC) concluded that an uninfected traveler is at greatest risk if seated within two rows of a contagious passenger [8, 9].

Modeling viral transmission in airplanes in extremely difficult, but efforts have been made to do so approximately. Infections can spread via respiratory droplets or aerosols, and models of such spread in aircraft have typically focused on one of these modes of transmission. For aerosols alone, several simulations of viral diffusion have been performed based on the Navier-Stokes equation [10,11,12,13]. For droplets, the studies include [14,15,16,17], but only the first of these is specific to passenger aviation.

Several papers have discussed actual experiences with viral transmissions on airplanes.

[18,19,20,21,22,23,24,25,26,27] In this article, we have attempted to locate the full set of peer-reviewed papers with data about transmission of Covid-19 in flight. We then focused on all nine of those papers which reported which seat(s) were occupied by travelers with contagious Covid-19, which nearby seats were occupied, and which nearby travelers contracted Covid-19 during the flight.

Furthermore, much recent research has appeared concerning the course of Covid-19 among individuals infected with it. This research estimates what fraction of carriers are eventually symptomatic [28], the period of contagiousness both before and after symptoms occur [29], the distribution of symptomatic cases by severity [30], and the duration of contagiousness of asymptomatic carriers [31, 32]. These various patterns are germane to estimating the likelihood that a passenger boarding a plane has a contagious case of Covid-19.

3 The general model

We focus on Boeing 737 or Airbus 320 jet aircraft in an all-coach configuration with 174 seats, in 29 rows with six seats apiece. In each row, there are two sets of three seats on opposite sides of a center aisle (seats ABC to the left of the aisle, DEF to the right). We consider nonstop US domestic flights of two-hour duration over the nine-month period June 2020–February 2021. The basic equation for estimating the probability P, that a randomly-chosen uninfected traveler contracts Covid-19 in flight is

$$P=\sum_{N,X}{Q}_{N,X}{P}_{N,X}$$
(1)

where

P:

P(randomly-chosen uninfected passenger contracts Covid-19 in flight)

PN, X:

P(randomly − chosen uninfected passenger boards flight with N passengers,  X of whom board harboring Covid-19)

QN, X:

P(randomly chosen uninfected passenger contracts Covid − 19 on board | N, X)

That P depends on X, the number of contagious passengers on board the flight, is obvious. The total number of passengers (N) matters because, the greater the number on board, the larger in general is the fraction of travelers seated in close proximity to contagious ones.

Below we first estimate PN, X and then the conditional probability QN, X.

3.1 The distribution of N

Suppose that M is the number of nonstop flights on 737/A320 aircraft (with an estimated 174 seats) in a particular month, and let qZ be the proportion of those flights with exactly z passengers on board. Then the total number of passengers carried on flights with z passengers would be Mzqz, and the fraction Qz of all travelers that month aboard such flights would follow:

$${Q}_z=P\left(N=z\right)=\kern0.5em Mz{q}_z/\sum_{j=0}^{174} Mj{q}_j=\frac{z{q}_z}{\sum_{j=0}^{174}j{q}_j}$$
(2)

In practice, we work with data for small ranges of z rather than individual values. As we will discuss, we had access to data about seat-occupancy rates by month on US aircraft between June 2020 and February 2021.

3.2 The distribution of X

For a particular value of N and a particular day, let pCONTAG be the probability that a randomly-chosen passenger boarding one of the flights under study carries contagious Covid-19. We treat different passengers as independent in terms of contagiousness status, meaning that X, the number of contagious travelers on board, would follow a binomial distribution with parameters N and pCONTAG . Estimating pCONTAG is therefore the key to approximating the distribution of X.

Our general strategy for estimating pCONTAG is:

  • Estimate the fraction of Americans who carry Covid-19 at the time of the flight

  • Adjust that fraction to reflect the fact that air travelers are more affluent than average Americans, and the Covid-19 infection rate is negatively correlated with income

  • Further adjust that fraction to reflect the likelihood that carrying Covid-19 reduces the likelihood of actually boarding a plane. Take into account that this statement is least true for those who are asymptomatic/pre-symptomatic, and most true for those with serious/critical cases.

3.3 The fraction γ of potential air passengers who carry Covid-19

To estimate pCONTAG, we first estimate the fraction of potential US air travelers who carry Covid-19 at a given time. The general equation for γ is:

$$\gamma =\left(\frac{C_7}{POP}\right)\rho \beta .$$
(3)

Where

C7:

Number of confirmed cases of Covid − 19 in the US over the past seven days

POP:

Population of the US

ρ:

a multiplier to adjust for the (sizable)extent to which confirmed Covid − 19 cases in the US underestimate the actual number of cases

β:

a"healthy passenger" factor to reflect the fact that US air travelers come disproportionately from communities where Covid − 19 is less common than in the population at large.

The quantity C7 enters (3) because seven days is approximately the duration of contagiousness for a carrier of Covid-19 [29, 32]. Thus, based on confirmed cases, \(\frac{C_7}{POP}\) approximates the fraction of Americans who are contagious at a particular time. Because we are considering US domestic flights as a whole, it seems reasonable to work with national statistics about new infections.

However, a substantial share of new cases of Covid-19 do not get confirmed in the US, meaning that actual cases over a seven-day period considerably exceed confirmed cases. There is no official guidance in the US about the ratio of actual-to-confirmed new cases at any time.. Here the quantity ρ is an estimate of that ratio, which we approximate on a weekly basis over the nine-month observation period.

3.4 Underestimation of US Covid-19 cases (the multiplier ρ for confirmed cases)

To estimate ρ, we follow a process of “reverse engineering,” using death counts at the particular times to estimate numbers of new cases at previous times. Suppose that η is the fraction of Covid-19 cases end in death. Thus, having estimated how many cases that started at a given time ultimately caused deaths, we could multiply by 1/η to approximate the total number of new cases at that time.

We made the approximation based on literature that those who died on a given date contracted their fatal cases of Covid-19 about uniformly over a four-week period that ended fourteen days earlier [30]. Thus, for new cases on day t that ultimately ended in death, these deaths occur in about equal numbers on each day from t + 14 to t + 41. However, on (say) day.

t + 18, there would also be deaths resulting from new cases on days t – 23 to t + 4. We estimate the number of new infections on day t that turned fatal under the rule:

$${N}_D(t)\kern0.5em \approx \left(\frac{1}{28}\right)\kern0.5em \sum_{j=14}^{41}D\left(t+j\right)$$
(4)

Where

ND(t):

Number of new US Covid − 19 cases on day t that ended in death

D(t + j):

recorded number of US Covid − 19 deaths on day t + j

In essence, (4) says that because day t was one of 28 contributors to the death toll on each day from t + 14 to t + 41, it contributed about 1/28 of the total deaths of that period. On average, that rule seems plausible.

Once ND(t) is at hand, it would be scaled up by the aforementioned factor 1/η to approximate the total number of new cases on day t. η would be based on estimates in the literature on the fraction of US Covid-19 cases over June 2020–February 2021 that ended in death. ρ would then be approximated by the ratio of the number of new cases based on death counts to the number of confirmed cases (C7):

$$\rho \approx \left(\frac{1}{\eta}\right)\frac{N_D(t)}{C_7}$$
(5)

3.5 “Healthy traveler” factor β

We consider a region that is divided into various communities, and contrast the overall Covid-19 rate for the region with a rate that weights individual communities by their shares of air travelers. The estimation uses four quantities for each community for 2020: its population, its median income, its per-capita rates of confirmed Covid-19 infections, and the per-capita annual rates at which its residents took air trips (before the Covid-19 pandemic).

If the communities are then combined into n ranges based their median incomes, the chance that a randomly-selected citizen has experienced Covid-19 can be approximated by:

$${R}_{CV}=\sum_{i=1}^n{R}_i{I}_i$$

where

RCV:

probability of having had Covid-19.

Ri:

probability of having had Covid − 19 for citizens of communities in median −income range i.

Ii:

fraction of all citizens living in communities in median − income range i

(n):

number of median − income ranges

Focusing on air travelers rather than randomly-chosen citizens, suppose that Ti = annual air trips per capita in communities in median − income range i

Then \({T}_i{I}_i/\sum_{j=1}^n{T}_j{I}_j\) is the fraction of pre-pandemic air trips taken by people from communities in median-income range i. If we weight each range’s per-capita Covid-19 rates by shares of pre-pandemic air trips, we reach an adjusted risk level RCVT that is given by:

$${R}_{CVT}=\sum_{i=1}^n{R}_i{T}_i{I}_i/\sum_{j=1}^n{T}_j{I}_j$$

R CVT is the probability of having had Covid-19, giving greater weight to communities in the income ranges with disproportionate shares of air travelers.

An estimate of the “healthy traveler” factor β would be:

$$\beta \approx {R}_{CV T}/{R}_{CV}$$
(6)

As we will discuss, we calculated a β−estimate with data from the Commonwealth of Massachusetts, and made another estimate using data from England and Wales.

This approximation for β is not perfect. It treats the pre-Covid distribution of air trips by income as a good proxy for the distribution during the Covid-19 pandemic, and bases trips per capita in a given community on its median income rather than its distribution of income. Avoiding such assumptions is not feasible given the limitations of available data.

\(\mathrm{As}\ \mathrm{noted}\ \mathrm{earlier},\gamma =\left(\frac{C_7}{POP}\right)\rho \beta \kern0.75em\) However, γ is not in itself the probability pCONTAG that we seek, as we explain below.

3.6 The probability p CONTAG that a boarding air passenger carries Covid-19

Suppose that α is the probability that a potential air traveler not suffering from Covid-19 would fly on a given day during the pandemic. (This quantity takes account of the general reluctance to travel because of the pandemic.) Among those who carry Covid-19—the number of which depends on γ−we posit that their chance of flying is α times a filtering parameter f that is less than one. The presumption is that some disease carriers refrain from flying because they suspect or know that they are infected with Covid-19 or simply feel too ill to travel.

We treat f as a function of several other quantities. An individual who carries Covid-19 who might potentially fly would fall into one of four categories:

Category

P(someone carrying Covid-19 is in that Category)

Asymptomatic carrier

PASYM

Pre-symptomatic carrier

PPRESYM

Carrier with mild/moderate symptoms

PMS

Carrier with severe/critical symptoms

PSC

In Section 4, we will estimate these probabilities based on medical literature.

We define three quantities related to potential travelers who harbor Covid-19

ψ NS :

the probability that someone with asymptomatic/pre-symptomatic Covid-19 who would otherwise fly will not do so

ψ MS :

the corresponding probability for someone showing mild/moderate symptoms of Covid-19

ψ CS :

the corresponding probability for someone showing severe/critical symptoms of Covid-19

The quantities ψNS, ψMS, and ψCS  ψNS are “desistance” probabilities, which reflect an unwillingness to fly among Covid-19 carriers beyond the general level of unwillingness signified by the parameter α. ψNS is relevant because someone with asymptomatic/pre-symptomatic Covid-19 might not fly because he knows or suspects that he carries the virus. ψMS is the chance that someone with mild/moderate symptoms foregoes the flight, and ψCS is the chance for severe/critical carriers. One can reasonably postulate that ψNS < ψMSψCS, but assigning probability distributions to ψNS, ψMS, and ψCS is not easy.

If αf is the probability a potential air traveler who carries Covid-19 nonetheless flies, the filter f would follow:

$$f=\left({P}_{ASYM}+{P}_{PRESYM}\ \right)\left(1-{\psi}_{NS}\right)+\kern0.75em {P}_{MS}\left(1-{\psi}_{MS}\right)+{P}_{SC}\left(1-{\psi}_{CS}\right)$$
(7)

The conditional probability pCONTAG that a passenger who boards a flight carries Covid-19 would then follow:

$${\boldsymbol{p}}_{\boldsymbol{CONTAG}}=\frac{\alpha f\gamma}{\alpha f\gamma +\alpha \left(1-\gamma \right)}=\frac{f\gamma}{f\gamma +\left(1-\gamma \right)}$$
(8)

where γ and f arise from (3) and (7)

Combining the various elements above, we reach:

$${P}_{N=j,X=k}\approx \kern0.5em P\left(N=j\right)\ast P\left(X=k\ \right|\ X\ is\ B\left(N,{\boldsymbol{p}}_{\boldsymbol{CONTAG}}\right)\Big)$$
(9)

where P(N = j) arises from (2), and B(N, q)is binomial with parameters N and q)

3.7 Estimation of QN,X

Given that there are N passengers on board, X of whom carry Covid-19, QN,X is the conditional probability that a random-chosen passenger from among the uninfected N – X will contract Covid-19 during the flight. We consider first the case in which X = 1, with a contagious passenger who is seated in seat S1.

For an uninfected in another seat S2, we posit that the probability of an on-board infection depends on:

  • The proximity of S2 to S1, based on both distance and the presence of any transmission barriers between the seats

  • The duration of the flight

  • The failure rate of universal masking

We further assume that short interactions among passengers during boarding, deplaning, or en route to the lavatory present second-order risks relative to those of sitting for hours close to a contagious person. We do so because CDC advised that someone exposed to a person with confirmed Covid-19 need not go into quarantine if the interaction was shorter than 15 minutes If, however, one believes that these hazards are not negligible, they would add to the risks we are estimating here.

Directly estimating the Covid-19 transmission probability on aircraft is very difficult. It depends on the contagious passenger’s emissions of the virus via a mixture of breathing, speaking, coughing or sneezing, as well as on the movement of droplets and aerosols given the geometry of the jet airplane and its powerful HEPA air-purification systems. In modeling QN,X, we lean heavily on two influential journal articles:

  • a dynamic-network simulation of viral transmission on US jet aircraft which appeared in the Proceedings the National Academy of Sciences [14] It assumes a fixed risk per minute that contagious passenger infects a given noninfected passenger within one meter (with a point estimate of 1.8% per minute)

  • a meta-analysis which appeared in The Lancet [5], and estimated that, absent barriers and masks, risk of viral infection (per minute) drops off exponentially with distance, at roughly a factor of two per meter.

Our general equation for estimating risk per minute given S1 and S2 takes the form:

$${V}_{S_{1,}{S}_2}={\pi}_{0}{e}^{-\omega d}{\left(1-\lambda \right)}^R\left(1-{p}_{masks}\right)$$
(10)

Where

d:

grid distance between contagious and given uninfected passenger

ω:

rate of exponential decay of risk with distance assuming no barriers

λ:

P(an individual seatback blocks transmission from a contagious passenger)

π0:

level of risk at distance zero from contagious person

R:

number of rows between contagious and uninfected passengers

pmasks:

P(universal masking blocks the transmission of Covid-19 from a contagious passenger to an uninfected one)

We amplify on (10) below.

In measuring d, we use not the Euclidean distance between S1 and S2 but instead the “grid distance” (i.e., d((x1, y1) to (x2, y2)) = |x2 – x1| + |y2 – y1|). This choice reflects the assumption that emissions from a contagious traveler in 16A that reach the breathing space of a passenger in 15B largely travel there via the breathing space for the passenger in 15A.

In including the quantity λ in (10), we recognize that the seatback directly ahead of a disease sufferer can block some forward transmission of aerosols and droplets, while that person’s own seatback can somewhat protect passengers one row behind. There is apparently no literature about the health benefits of airplane seatbacks, though an article in a medical journal informed potential air travelers that “seat backs provide a partial physical barrier [33].) We define λ as the probability that a single row of seatbacks prevents transmission, in which case 1 - λ is the probability that it fails to do so. If there are R rows separating S1 from S2, we assume that the overall failure risk is (1 − λ)R because R rows of seatbacks all must fail for transmission to occur. R = 0 if both S1 and S2 are in the same row.

The model in (9) follows [14] in assuming a constant infection risk per minute assuming no prior transmission, and follows [5] in assuming exponential decay of risk with distance. In the next section, we estimate the parameters π0, ω, and pmasks based on the literature, assigning them probability distributions as well as point estimates.

For a flight of duration T in minutes, the probability \({Z}_{S_1,{S}_2}\) of Covid-19 transmission conditioned on S1 and S2 follows:

$${Z}_{S_1,{S}_2}=1-{\left(1-{V}_{S_{1,}{S}_2}\ \right)}^T$$
(11)

If N passengers occupy N seats aboard a flight, we assume that the single contagious passenger is equally likely to be in any of those seats while a randomly-chosen uninfected passenger is equally likely to be in any of the other N – 1 seats. Therefore, U(S1), the probability conditioned on S1 that a random uninfected passenger contracts Covid-19 on board would be given by:

$$U\left({S}_1\right)=\left(\frac{1}{N-1}\right)\sum_{all\ {S}_2\ne {S}_1}\ {Z}_{S_1,{S}_2}$$

The overall risk of infection QN, X = 1 would follow:

$${Q}_{N,X=1}=\left(\frac{1}{N}\right)\sum_{S_1= Seat\ 1}^{Seat\ N}U\left({S}_1\right)$$
(12)

where the seats 1 through N are numbered under some convention

If N = 174 on the 737/A320 jet under consideration, then all seats are full. As N drops below 174, we assume that as many of the empty seats as possible are center seats.

Given that QN, X = 1is small (as we will see), the risk of contracting Covid-19 in flight can be treated as essentially linear in the number of contagious people who board. Thus, based on (12) we can write

$${Q}_{N,X=k}\kern0.75em \approx \kern1.25em \left(\frac{k}{N}\right)\sum_{S_1= Seat\ 1}^{Seat\ N}U\left({S}_1\right)$$
(13)

Given (1), (9), and (13) , we have an estimate of P.

4 Probability distribution for key parameters

Here we pursue the general approaches just outlined to get probability distributions for each model parameter. These distributions will underpin the simulation analysis in the next section.

4.1 Distribution of ρ, (actual Covid-19 cases/confirmed cases)

The “reverse engineering” procedure depicted by eq. (4) generates estimates of the average daily number of new fatal US cases of Covid-19 from June 2020 to February 2021. To reach a total number of new cases, one needs divide new fatal cases by the fraction of infections that lead to death (the case fatality rate). That fraction is not precisely known. Baud et al. estimated in June 2020 that 0.7% of Covid-19 infections in the US ended in death [34]. Improvements in treatment for Covid-19 are believed to have lowered the case fatality rate during 2020 by perhaps 20–33%, to about 0.5% [35]. Ioannidis classified the US among nations with a case fatality rate of 0.57% [36].

We estimate ρ on a weekly basis to avoid day-of-week variations in confirmed cases. For a point estimate of the ratio ρ in a given week, we first estimated total deaths that week using [37]) for each of the seven days. Then we multiplied that quantity by 167 to reach a point estimate of total cases, based on a case fatality rate of 0.6% (about the average of the estimates above). We approximated ρ as the ratio of total to estimated cases. We assumed (approximately) a normal distribution for ρ, treating the point estimate based on a 0.6% fatality rate as its median, the ρ − estimate based on a 0.7% rate as its 10th percentile, and the estimate based on 0.5% as the 90th percentile. Fig. 1 presents graphically the movement over time of the estimated 10th, 50th, and 90th percentiles for ρ.

Fig. 1
figure 1

Estimated distribution of ρ, the ratio of new actual US Covid-19 cases to confirmed cases, June 2020–February 2021. Sources: References [34,35,36,37]

4.1.1 Distribution of \(\beta\ \left({}^{"}{\mathsf{Healthy}\ \mathsf{Traveler}}^{"}\kern0.5em \mathrm{effect}\right)\)

Table 1 presents data about both confirmed Covid-19 cases and estimated travel frequencies across several income ranges in Massachusetts, which can be processed in accordance with eq. (3). Somewhat surprisingly, the point estimate of β arising from the calculations is 0.86, which suggests that the “healthy traveler” effect is not all that large.

Table 1 Average annual airline trips and confirmed Covid-19 cases per capita by annual household income, estimates for commonwealth of Massachusetts

For further perspective about β,we use data from England [40] about Covid-19 death rates as a function of the deprivation index used to classify small geographic areas. Doing so yields an estimate of β for England and Wales of about 0.75.

A reasonable model that incorporates these two estimates and some further uncertainty could posit that β follows a beta distribution with parameters α = 28 and β = 7. This distribution has a mean of 0.8, a slightly higher mode, and a standard deviation of 0.066.

4.2 Distributions for ψNS, ψMS, and ψCS, contributors to the “filtering” effect of carrying Covid-19 on the likelihood of air travel

Assigning probability distributions to the “deterrence” quantities ψNS, ψMS, and ψCS is difficult, because no readily-available data allow for direct estimates. One can, however, reasonably postulate that ψNS < ψMSψCS, and try to work with scant existing “clues” about the magnitudes of these parameters.

The quantity ψNS for carriers without symptoms would presumably be low. Such individuals could well decline to travel if they test positive for Covid-19, or have recently been close to others with known infections. But many US cases of Covid-19 are never confirmed, and people who do get tested may well come disproportionately from those with symptoms of Covid-19. For these reasons, it seems plausible to assign ψNS for carriers without symptoms a normal distribution with a low mean and modest standard deviation. As approximate parameters for that distribution, we advance a mean of 0.3 and standard deviation of 0.05. (Recall that ψNS = 0.3 implies a 70% probability of going ahead with travel.)

Some potential travelers with mildly symptomatic Covid-19 will get themselves tested, and those with positive tests would not be expected to fly. But among potential passengers who are not tested and who have mild symptoms that are hard to distinguish from those of a common cold, a fraction of travelers will go ahead with their flights. We assign ψMS a normal distribution with mean 0.6 and a standard deviation .05.

One would expect that few individuals with severe/critical Covid-19 would board US flights. But “few” is not the same as none. In late 2020, two people are known to have died of Covid-19 on US domestic flights. We assign ψCS a normal distribution with mean 0.9 and a standard deviation of 0.033 (also imposing an upper limit of one).

4.3 Distribution of ω (drop in risk with distance)

As noted, the meta-analysis in [5] estimated that viral transmission risk declines exponentially with greater distance from the contagious person. The paper’s point estimate was that risk fell by a factor of 2.02 (equivalently, risk was multiplied by 1/2.02) for each additional meter of separation, with a 95% confidence interval extending from a factor of 1.08 to one of 3.76. This interval—from roughly half the point estimate to double it—is strongly suggestive of a lognormal distribution. We chose such a distribution, and chose its parameters to match the three statistics above (1.08, 2.02, 3.71) Given that ln(1/2.02) = −.703,

ln(1/1.08) = −.077, and ln(1/3.76) = −1.324, we assign the normally-distributed logarithm of the decay parameter ω a mean of −.703 and a standard deviation of .318. (The latter statistic is (1/1.96) times the average of 1.324–.703 = .621 and .703–.077 = .626.) For those parameters, ω itself is assigned a mean of .521 and a standard deviation of .170.

But what of the distance d to which ω would be applied? In economy class in a typical Boeing 737 or Airbus 320 jet, the seats are approximately 18 in. wide while the aisle width is about 30 in.. The seatbacks in consecutive rows are separated by about 31 in.. As noted, we tie viral transmission risk to grid distances between seats. On these aircraft, grid distances on these aircraft in inches from people within three rows of a contagious person in seat 16A are:

 

Seat

Row

A

B

C

D

E

F

13

93

111

129

159

177

195

14

62

80

98

128

146

164

15

31

49

67

97

115

133

16

0

18

36

66

84

102

17

31

49

67

97

115

133

18

62

80

98

128

146

164

19

93

111

129

159

177

195

Similar charts arise when the contagious passenger is in a B or C seat and, by symmetry, in the D, E, or F seats.

4.4 The parameter λ (seatbacks as transmission barriers)

Reflecting considerable uncertainty about the benefits of seat backs, we assign λ--the probability that a row of seatbacks blocks viral transmissions to a preceding or following row of passengers--a normal distribution with a mean of 0.5 and a standard deviation of 0.1. Then 1 - λ is the failure rate of one row of seatbacks to prevent such transmission.

Passengers three or more rows from the contagious person are separated from her by a minimum of 93 in. and three seatbacks. With the distributions we use for ω and λ, transmission risk is so attenuated under three-row separation that we treat it as a second-order effect. Thus, we effectively adopt the two-row rule that, as noted earlier, has been advanced by the World Health Organization and the CDC.

4.5 Distribution of p masks(benefit of masks)

Some passengers on US domestic flights used cloth masks over the study period, while others wear surgical masks. (Powerful N95 or KN95 masks were not generally available during that period.) The meta-analysis in [5] estimates that wearing of cloth or surgical masks reduces by an estimated 67% the chance that a viral infection is successfully transmitted. A study from China [41] about Covid-19 yields the estimate that, when everyone wears such masks, overall infection risk drops 79%. Another paper about viral infections [33] estimates that universal cloth-mask wearing would cut transmissions by 75% while surgical masks would cut transmission by 94% .

Giving by far the greatest weight to the meta-analysis, we approximate the risk multiplier pmasks as approximately normal with a mean of 0.3 (i.e., a 70% reduction) and a standard deviation of 0.075.

4.6 The parameter τ0 (transmission risk per minute)

The estimation of τ0 − the transmission risk per minute of exposure zero distance from the contagious passenger and with no masks--is challenging. One might assume that τ0 can vary from flight to flight, because a contagious passenger who speaks and sneezes throughout the trip might cause more infections than another who sleeps all flight long. Here we rely on several research papers about actual in-flight transmissions of Covid-19 ([19,20,21,22,23,24,25,26,27]), the results of which are summarized in Table 2. While these papers all concern jet flights outside the United States, the flights had functioning HEPA air filters similar to those on US jets. As we discuss in the Appendix, we assign τ0 a beta distribution with parameters α = 1 and β = 520.

Table 2 Some details about ten jet flights on which some boarding passengers carried Covid-19 (Sources: References [19,20,21,22,23,24,25,26,27])

4.7 Estimation of P ASYM, P PRESYM, P MM, and P SC

We follow a different approach in estimating the parameters about the course of contagious Covid-19, namely, PASYM, PPRESYM, PMM, and PSC. Consistent with available literature [28,29,30,31,32], we treat current carriers of Covid-19 as having become infected over the last seven days, and assume that about 30% of them are asymptomatic. Among the 70% infected in the last week who will develop symptoms, we use the estimate that those infected in the last two days (i.e, approximately 2/7 of the group) are pre-symptomatic, while the remaining 5/7 show symptoms of varying severity. Similarly, we work with published estimates that approximately 4/5 of symptomatic carriers will develop mild/moderate symptoms, while 1/5 will have severe/critical ones. For these reasons, we make the approximations:

$${\displaystyle \begin{array}{cc}{P}_{ASYM}\approx 0.3& {P}_{PRESYM}\approx 0.7\ast \left(\frac{2}{7}\right)=0.2\\ {}{P}_{MM}\kern0.5em \approx 0.7\ast \left(\frac{5}{7}\right)\ast 0,8=0.4& {P}_{SC}\approx 0.7\ast \left(\frac{5}{7}\right)\ast 0.2=0.1\end{array}}$$

In the analysis, we worked with these point estimates of Covid-19 infection parameters, rather than with probability distributions. We did so in the belief that adding several more distributions to the modeling would do more to make the formulation unwieldy than to offer greater insight about Covid-19 transmission on airplanes.

5 Results

  1. (i)

    Individual Flights

Given the number of random variables in the expression for PN,X and their differing distributions, the distribution for PN,X does not take any familiar analytic form. It is necessary to use simulation to investigate the behavior of PN,X and we did so, treating different random variables as independent. In most respects, an independence assumption is highly plausible: the uncertainty about the “healthy traveler” parameter seems unrelated to uncertainty in the effectiveness of seatbacks.

We conducted 10,000 simulations for each month from June 2020 to February 2021, assuming an all-economy configuration of the 737 or 320 aircraft with 29 rows and 174 seats. We considered six possible occupancy levels for the seats:

  • All seats full

  • 5/6 of seats full, with half the middle seats empty but all other seats occupied

  • 2/3 of seats full, including all window and aisle seats occupied but none in the middle

  • Half of seats full, consisting of three in each row scattered among the four window and aisle seats

  • 1/3 of seats full, consisting of two in each row, one on each side of the aisle in the window or aisle position

  • 1/6 of seats full, consisting of one window or one aisle seat in each row

Table 3 presents various percentiles for the probability that an uninfected passenger would contract Covid-19 during a two-hour flight on an “average” day over June 2020–February 2021. (i.e., one for which (N7 /N) was at the average for the period). The table reports on three possibilities: the flight was full, 2/3 full, or 1/3 full. In all of the columns, the distribution for risk has high variability, with the 10th percentile about 1/7 of the median and the 90th percentile nearly four times the median. Such variability was unavoidable because several of the “input” distributions in the simulation had high coefficients of variation.

Table 3 Estimated distribution for Covid-19 infection probability for randomly-chosen uninfected travelers on November 1, 2020, for two hour US domestic jet flights on Boeing 737 or Airbus 320, at three levels of seat occupancy

Table 4 presents mean and median risk at all six levels of crowding that were considered. (Rather than designate the mean or the median of the simulation results as the point estimate for infection risk, we give equal prominence to both.) The mean and the median risk estimates diverge somewhat, presumably because several distributions for the parameters had long right tails. But both the mean and median suggest a risk around 1 in 1000 of contracting Covid-19 in flight on a full plane on an average day. The table also indicates that, as the flight becomes less full, infection risk drops more than proportionately. For example, the mean risk estimate on a 55% full flight is not 55% of that for a full flight, but instead 46% of it. If the flight is 55% full, then its expected number of contagious passengers is 55% as high as that for a full flight. But the uninfected passengers are on average further away from the contagious one (especially because the seat next to that traveler is empty). For that reason, risk declines not by 45% but by 54%.

Table 4 Mean and median estimates of Covid-19 infection probability for two hour US domestic jet flights on Boeing 737 or Airbus 320 on November 1, 2020, at six levels of seat occupancy

Tables 5 and 6 elaborate on this last point by offering estimates of the probability of becoming infected during a two-hour flight given a contagious traveler is in seat 16A, as a function of the uninfected traveler’s seat in rows 14–18 (and assuming mask usage). Under the “two row” approximation, passengers in other rows are assigned infection probabilities of zero.

Table 5 Estimated risk of contracting Covid-19 on two-hour US domestic jet flights on Boeing 737 or Airbus 320 for each seat position, given that the flight is full and a single contagious passenger is sitting in Seat 14A
Table 6 Estimated risk of contracting Covid-19 on a two hour US domestic jet flights on Boeing 737 or Airbus 320 for each seat position, given that all of the middle seats on the flight are vacant and a single contagious passenger is sitting in Seat 14A

Fig. 2 shows that the contagious passenger infects more other travelers if she is in an aisle seat than in a middle seat, while the middle seat yields more infections than a window seat. However, the effect is not especially large.

  1. (ii)

    Flights over the Nine-Month Period

Fig. 2
figure 2

Probability of contracting Covid-19 on two hour domestic US jet flight on 11/1/20, based on type of seat of single contagious passenger on fully-loaded Boeing 737 or Airbus 320

So far, we have considered the general question: if pCONTAG is the probability that each boarding passenger has contagious Covid-19 and N is the total number of passengers on the plane, what is the probability that an uninfected passenger will contract Covid-19 in flight? But we have not considered what proportions of actual passengers faced which specific (N, pCONTAG) combinations. Here we try to address this issue by merging the existing “what if?” analysis with data about the ebbs and flows of the Covid-19 pandemic over the nine months of interest, as well as actual seat occupancy rates on US jets during that period.

We are fortunate to have obtained data about the distribution of seat occupancy by for US airlines month from June 2020 to February 2021, broken down by deciles. We can use these data in estimating overall infection risk if we make two assumptions (i) that average occupancy level within a decile is at the middle of that decile (e.g., 55% for the range 50%–60%) and (ii) the distribution of occupancy on 737 and 320 flights is about the same as the distribution over all US flights.

Table 7 offers seat-occupancy data for the month November 2020, along with a tabulation of infection risk for each occupancy range as a multiple of the risk for a full flight.

Table 7 Some details about US domestic jet flights in November 2020, related to the risk of getting infected with Covid-19 on board

To move from this information to infection risk for travelers, we work with eq. (2). Doing so properly reflects the points that a plane with 174 seats that is 80% full carries four times as many passengers as another such plane that is 20% full .

Table 8 presents mean and median estimates of infection risk by month for randomly-chosen travelers of 737 or 320 jets. The risk estimates vary considerably over time, and are highest for January 2021 when rates of new Covid-19 infections in the US were surging. Both the mean and median risk estimates are 1 in 2250 or higher for the full period. Not everyone would agree with Southwest Airlines that this level of risk is “virtually nonexistent.”

Table 8 Mean and median estimates of risk that randomly-chosen uninfected traveler contracts Covid-19 on board US Boeing 737 or Airbus 320, by month from June 2020 to February 2021

It is natural to consider how many US deaths resulted from Covid-19 infections contracted on airplanes. But estimating that number is extremely difficult. Moreover, a substantial fraction of the deaths ultimately caused by on-board Covid-19 infections could befall people who were never on the plane.

One can make rudimentary calculations tied to the estimate that 1 in 2000 uninfected air travelers contracted Covid-19 on US domestic flights over June 2020 to February 2021. Given 204 million passengers in total, the 1 in 2000 rate implies about 100,000 infections in total. If one believes that every 1000 on-board infections ultimately led to one death, then the death toll would be 100. For one death per 500 on-board infections, the corresponding toll would be 200; for one death per 100 such infections, the death toll would reach 1000. While such numbers fall far below the hundreds of thousands of US deaths in the pandemic, even 100 fatalities over a nine-month period compares unfavorably to the average of one passenger death every nine months on US flights over 2010–19.

6 Some other assessments of in-flight Covid-19 transmission risk

As noted, two major studies about on-board transmission of Covid-19 on airplanes did not use data about actual passengers. Here we discuss them briefly.

The Aviation Public Health Initiative of the Harvard School of Public Health extensively studied the risk of Covid-19 transmission during air travel [42]. The study offered many good practical suggestions about flying during the pandemic, but only limited guidance about the magnitude of on-board infection risk. Its overall assessment was the risk is lower than that “in indoor restaurants or grocery stores,” but these were known to be especially hazardous venues during the pandemic. Moreover, HSPH’s assertion that on-board ventilation systems “effectively counter the proximity travelers will be subject to during flights.” was an overstatement given the evidence in Table 3.

The US Department of Defense [43] conducted an imaginative study about viral transmissions of disease on aircraft. It used mannequins wearing surgical masks to simulate coughing passengers and used aerosol tracers to see where the emitted particles went. However, the mannequins in the study faced forward and did not “speak” or move, while possible transmission via droplets was not considered. The researchers suggested that on-board infection risk was low but offered no numerical estimate.

7 Final remarks

We readily acknowledge that this analysis is imperfect, but hope that its strengths push it to the positive side of the ledger. It seemed worthwhile to use various data sets to estimate the prevalence of contagious Covid-19 among Americans who boarded airplanes over June 2020 to February 2021 It seemed sensible to use peer-reviewed reports about Covid-19 transmission on airplanes as guides to the likelihood that on-board-infections will occur. And it seemed desirable to work with actual seat-occupancy data, so that the modeling could assess the experiences of typical passengers, and not just worst-case scenarios like full flights at the height of the pandemic.

The situation presumably changed with the widespread availability of vaccines starting in early 2021, but also with the sharp increase in contagiousness of Covid-19 with the arrival of the Delta and Omicron variants. Modeling like that presented here could help in assessing the changed situation, much as the general approach might help in connection with a future pandemic. And that pandemic might arise sooner than any of us would like.