Precocity of Covid-19 outbreaks
In this section, we test whether the precocity of the outbreaks of the epidemic could be predicted on the basis of the inflows of airline passengers from the hotbed of Covid-19 (mainland China) or through other air travel routes. In Table 1 we focus on the calendar week of occurrence of the first death in our sample, which ranges from 7 (second week of February 2020, in Southern Canto, Japan) to 23 (last week of May 2020, in New Brunswick, Canada). We estimate models that include passenger volume from Chinese airports and from other world airports separately (tables 1A and 1B) to avoid collinearity issues since the two figures are strongly correlated (r=.60). As a sensitivity check, we also repeat the analysis for the week of occurrence of the 10th Covid-19 death (Models 5-8). Negative coefficients indicate an earlier outbreak. Models 1 and 5 in Table 1 are estimated with mobility-related variables only. The volume of incoming passengers from flights originating from China and from other world airports is prima facie significant, while the area’s network centrality in global air traffic is not. Passenger volume is no longer impactful once we take into account population size, which substantially increases the likelihood of having a Covid-19 outbreak earlier, as well as population density and GDP per capita (models 2, 3 and 7). Demographic and health characteristics – in particular, a larger senior population and a higher prevalence of obesity – appear to amplify the likelihood of mortality at the onset of the epidemic (models 4 and 8). These effects are quite similar whether we estimate them in conjunction with the volume of passengers from China or from the rest of the world. Interestingly, pre-pandemic cardiovascular and cancer death rates are in fact inversely related to the precocity of the onset of Covid-19 deaths, possibly because areas with a higher incidence of these diseases are more prepared to tackle medical emergencies. Equally, a larger availability of hospital beds – which is likely to attest more broadly to the strength of the health system – also counters an earlier spread up to the 10th recorded death.
The same analyses were repeated with the occurrence of the 1st and 10th Covid-19 cases – which are known to be less reliable indicators of the outbreak than death counts, due to different testing and recording efforts across countries. We obtained similar results, although with stronger effects for the volume of passengers (see Supplementary information). We also found that these results were robust to the exclusion of different countries. Four of our sub-national areas are outliers (two in Japan, two in the US), as they receive an extremely high number of air passengers (from China or from the rest of the world, respectively) while having a relatively late occurrence of Covid-19 deaths. If we exclude these areas, the effects of incoming air passengers are larger, but still not significant once we control for all other covariates in our models (see Supplementary information).
On balance, we have to nuance H1 (the ‘travel as timer’ hypothesis). The risk of experiencing an earlier outbreak of Covid-19 deaths is not conditional on airport centrality. The volume of incoming passengers (from China and elsewhere) is at best no more than a concomitant cause with a larger, older, and overweight population – and possibly the randomness of the distribution of a few super-spreaders in the global network of long-distance travels.
Table 1. The impact of late 2019 air passenger traffic on the precocity of Covid-19 outbreaks (occurrence of 1st and 10th death) [1A: traffic from Chinese airports; 1B: traffic from other airports]. OLS regressions without and with continent fixed effects
1A. Airline passengers from Chinese airports
|
(1)
|
(2)
|
(3)
|
(4)
|
(5)
|
(6)
|
(7)
|
(8)
|
|
Wk 1st
|
Wk 1st
|
Wk 1st
|
Wk 1st
|
Wk 10th
|
Wk 10th
|
Wk 10th
|
Wk 10th
|
Inbound Chinese passengers pc (end of 2019)
|
-0.542**
|
0.042
|
-0.240
|
-0.043
|
-0.567
|
0.573
|
-0.042
|
-0.409
|
|
(0.175)
|
(0.204)
|
(0.185)
|
(0.181)
|
(0.348)
|
(0.386)
|
(0.367)
|
(0.358)
|
Airport centrality (end of 2019)
|
-0.156
|
-0.097
|
0.092
|
0.040
|
-0.289
|
-0.068
|
0.149
|
0.086
|
|
(0.191)
|
(0.192)
|
(0.168)
|
(0.158)
|
(0.380)
|
(0.363)
|
(0.331)
|
(0.309)
|
Population
|
|
-2.917***
|
-3.183***
|
-2.843***
|
|
-5.143***
|
-4.720***
|
-4.443***
|
|
|
(0.618)
|
(0.579)
|
(0.553)
|
|
(1.185)
|
(1.153)
|
(1.091)
|
Population density (pop per km²)
|
|
-0.373
|
-0.243
|
-0.349
|
|
-0.714
|
-0.710
|
-0.884*
|
|
|
(0.215)
|
(0.187)
|
(0.179)
|
|
(0.410)
|
(0.371)
|
(0.353)
|
Real GDP pc PPP
|
|
-0.420**
|
0.026
|
0.164
|
|
-1.954***
|
-0.599*
|
-0.365
|
|
|
(0.134)
|
(0.129)
|
(0.129)
|
|
(0.289)
|
(0.302)
|
(0.307)
|
Hospital beds per 1000
|
|
|
|
0.117
|
|
|
|
1.482**
|
|
|
|
|
(0.264)
|
|
|
|
(0.518)
|
Share of 65+
|
|
|
|
-0.711***
|
|
|
|
-1.385***
|
|
|
|
|
(0.181)
|
|
|
|
(0.371)
|
Cardiovascular age-standardized death rate
|
|
|
|
0.978***
|
|
|
|
0.584
|
|
|
|
|
(0.260)
|
|
|
|
(0.535)
|
Cancer age-standardized death rate
|
|
|
|
0.490*
|
|
|
|
1.642***
|
|
|
|
|
(0.198)
|
|
|
|
(0.398)
|
Prevalence of adult obesity
|
|
|
|
-0.773
|
|
|
|
-3.572***
|
|
|
|
|
(0.396)
|
|
|
|
(0.795)
|
Observations
|
439
|
439
|
439
|
439
|
416
|
416
|
416
|
416
|
R-squared
|
0.048
|
0.133
|
0.359
|
0.444
|
0.018
|
0.192
|
0.353
|
0.448
|
Continent Fixed Effects
|
No
|
No
|
Yes
|
Yes
|
No
|
No
|
Yes
|
Yes
|
Standard errors in parentheses |
* p < 0.05, ** p < 0.01, *** p < 0.001 |
Source: Sub-National Covid-19 Incidence and Determinants Dataset |
1B. Airline passengers from non-Chinese airports
|
(1)
|
(2)
|
(3)
|
(4)
|
(5)
|
(6)
|
(7)
|
(8)
|
|
Wk 1st
|
Wk 1st
|
Wk 1st
|
Wk 1st
|
Wk 10th
|
Wk 10th
|
Wk 10th
|
Wk 10th
|
Inbound non-Chinese passengers pc (end of 2019)
|
-0.764***
|
-0.064
|
0.032
|
-0.078
|
-1.549***
|
0.093
|
0.095
|
0.647
|
|
(0.134)
|
(0.221)
|
(0.195)
|
(0.193)
|
(0.265)
|
(0.427)
|
(0.389)
|
(0.382)
|
Airport centrality (end of 2019)
|
-0.289
|
-0.078
|
-0.021
|
0.021
|
-0.198
|
0.229
|
0.132
|
-0.095
|
|
(0.151)
|
(0.162)
|
(0.143)
|
(0.136)
|
(0.297)
|
(0.306)
|
(0.281)
|
(0.265)
|
Population
|
|
-2.664***
|
-3.718***
|
-2.685***
|
|
-4.401**
|
-5.064***
|
-7.018***
|
|
|
(0.801)
|
(0.721)
|
(0.730)
|
|
(1.529)
|
(1.424)
|
(1.438)
|
Population density (pop per km²)
|
|
-0.366
|
-0.244
|
-0.339
|
|
-0.736
|
-0.720
|
-0.955**
|
|
|
(0.217)
|
(0.190)
|
(0.181)
|
|
(0.414)
|
(0.374)
|
(0.355)
|
Real GDP pc PPP
|
|
-0.406**
|
0.0113
|
0.178
|
|
-1.988***
|
-0.628
|
-0.552
|
|
|
(0.143)
|
(0.134)
|
(0.133)
|
|
(0.314)
|
(0.319)
|
(0.322)
|
Hospital beds per 1000
|
|
|
|
0.108
|
|
|
|
1.550**
|
|
|
|
|
(0.264)
|
|
|
|
(0.518)
|
Share of 65+
|
|
|
|
-0.713***
|
|
|
|
-1.383***
|
|
|
|
|
(0.181)
|
|
|
|
(0.370)
|
Cardiovascular age-standardized death rate
|
|
|
|
0.981***
|
|
|
|
0.614
|
|
|
|
|
(0.259)
|
|
|
|
(0.532)
|
Cancer age-standardized death rate
|
|
|
|
0.484*
|
|
|
|
1.609***
|
|
|
|
|
(0.197)
|
|
|
|
(0.395)
|
Prevalence of adult obesity
|
|
|
|
-0.724
|
|
|
|
-3.664***
|
|
|
|
|
(0.399)
|
|
|
|
(0.795)
|
Observations
|
439
|
439
|
439
|
439
|
416
|
416
|
416
|
416
|
R-squared
|
0.094
|
0.133
|
0.357
|
0.444
|
0.088
|
0.188
|
0.354
|
0.451
|
Continent Fixed Effects
|
No
|
No
|
Yes
|
Yes
|
No
|
No
|
Yes
|
Yes
|
Standard errors in parentheses |
* p < 0.05, ** p < 0.01, *** p < 0.001 |
Source: Sub-National Covid-19 Incidence and Determinants Dataset |
Modeling the severity of Covid-19 deaths in April-July 2020
Although travel limitations were pervasive as a first defense against the spread of the new virus, air traffic was not completely stopped in the wake of the pandemic. Our dataset includes information on the number of incoming passengers along all flight routes. The bulk of them diminished drastically after March 2020. Radical drops – with the full cancellation of routes – affected Asia particularly (Figures 1A [April 2019] and 1B [April 2020]). Globally, the drop in the volume of passengers was dramatic (Figure 1C). At its lowest, in May 2020, the number of air passengers was 12 percent of the corresponding 2019 figure and until July 2020 it had recovered to no more than 34 percent. In May 2020, only 3 percent of routes had not experienced any decline in air passengers compared to the same month one year earlier (Figure 1D). However, air traffic was not entirely cancelled, as Figure 1B shows. Some passengers continued to move between most parts of the globe. Potentially, this might have led to higher incidence of the disease where larger numbers of incoming travelers continued to arrive, especially if their origin was from highly infected zones. Therefore, we not only model the impact of the number of incoming air passengers (per capita) on the receiving areas’ incidence of Covid-19 between April and July 2020, but also the interaction between this variable and the average incidence of the disease in the sending areas (Tables 2, 3 and 4).
Together with long-distance mobility, other factors are known to boost heterogeneity in the lethality of the pandemic across locations. Previous research on between-country variations in Covid-19 severity has consistently highlighted the effect of several pre-existing comorbidities including a senior population structure, a higher GDP per capita and a lower number of hospital beds25,26. The impact of population density, obesity and air pollution is empirically less clear, with some within-country ecological studies disproving their independent impact on Covid-19 mortality27,28. We include these potential predictors in the following analyses as Structural predispositions, with the only exception of air pollution, for which we have 20 percent missing cases (analyses including this variable are in the Supplementary information). There is also strong evidence that everyday human interactions are highly conducive to the spread of Covid-19 – and their reduction conversely reduces contagion5,29,30,31. We operationalized such Population mixing through the ‘Oxford-tracker modified stringency index’ of NPIs, lagged by three weeks32. Given that basically all forms of NPIs aim for a reduction of human interactions, we deem this single indicator to be the most parsimonious way of tapping changes in population mixing over time. This is also in light of research based on mobile positioning data that indicates overall compliance with such measures in the first semester of 202033. The index is modified by removing its international travel limitation component, as this is included among our covariates separately. Our core variables revolve around Long-distance mobility: air passenger volume per capita (lagged by one month), the average incidence in the sending regions of these passengers (lagged by one month), the centrality of the subnational area in the global network of commercial flights (lagged by one month), and the air travel restriction policy that was in place three weeks earlier. Air travel restriction policies are operationalized in three categories: total travel ban as a baseline, travel bans for some routes, and less disruptive limitations (like temperature-screening, or quarantines for incoming passengers, or no limitations at all). In a separate analysis, we replace travel restrictions with actual passenger drop (through the dummy variable used in Figure 1D: decline of more than 10 percent of travelers compared to one year earlier). Such a specification simulates a sort of counterfactual scenario in which travel continued in spite of the pandemic (see Supplementary information). Lags of three weeks are introduced assuming, from existing literature, an average Covid-19 infection-to-death delay of 20-23 days.34,35 We extended the lag to one month for our travel indicators given their monthly nature. Finally, we introduced a variable that accounts for Recursive effects – namely the calendar week of the first Covid-19 death (which was the dependent variable in Table 1). This factor incorporates what has been called the ‘surprise effect’ of Covid-19, as a sooner outbreak could bring about a less efficient public health response36.
Weekly incidence: A generalized mixed model specification
In this section, we test H2 (the ‘travel as fuel’ hypothesis) on the impact of airline mobility by estimating a generalized mixed model (or multilevel model) which allows us to pool together all weekly information for the April-July 2020 period. Mixed models adjust for the nested data structure of our sub-national areas within countries, also taking into account the fact that we have longitudinal measures. In other words, each week for which we have data on deaths (level 1) is nested within each of the sub-national areas (level 2), which themselves are nested within countries (level 3) and within continents (level 4). By adjusting for time-based change, we can factor out heterogeneities in individual sub-national areas, as well as other existing heterogeneities between countries and continents37. Our specification also includes a time predictor (i.e., calendar week), as well as time-varying and time-invariant predictors. We fit a generalized linear mixed-effects model (GLMM) for the negative binomial family to our data with a random intercept and slope for the time indicator (calendar week). We assume a known overdispersion parameter of 1.2 (results do not vary much for other choices of the parameter; available upon request). One advantage of this approach is that it models changes over time at the local level simultaneously (e.g., how fast the pandemic evolves) and differences in time-wise change across different areas (e.g., differences in the peaks and valleys of the pandemic).
Table 2 reports the interaction between each of the dynamic variables (i.e., related to long-distance mobility and population mixing) and calendar week (for a stepwise version of the models, see Supplementary information). When interacting with time, the estimated risk ratios for long-distance mobility variables are close to zero. In fact, the effects of minimal travel bans (Model 1, Table 2) or no drops in passenger inflows (Model 2, Table 2) as opposed to full travel bans are above one. This attests to some impact on the severity of the pandemic, but not to a statistically significant level. Structural predispositions to the virus incidence – particularly, population size, GDP per capita, the age structure of the population, as well as the paucity of hospital beds, cancer-risk and obesity (albeit these are not significant statistically) – have a definitively greater bearing on the incidence of the disease. As expected, NPIs other than travel limitations mitigate the pandemic significantly. The precocity of the outbreaks is confirmed to be a strong and significant predictor of its subsequent severity. Finally, the estimated random intercepts and slopes from both models suggest that the higher levels (i.e., countries and continents) explain the total variance in weekly deaths more than the lower sub-national areas.
Table 2
The impact of air passenger traffic on the severity of Covid-19 mortality in April-July 2020, controlling for population mixing (NPI), structural predispositions and recursive effects. Generalized linear mixed-effects models
|
Model 1
|
Model 2
|
|
Risk
Ratio
|
Confidence intervals
|
Risk
ratio
|
Confidence intervals
|
Long-distance mobility (interacted by calendar week)
|
|
|
|
|
Incidence in sending region x inbound passengers
|
1.001
|
[0.986 - 1.015]
|
0.999
|
[0.984 - 1.014]
|
Airport centrality (lag 1 month)
|
1.002
|
[0.991 - 1.012]
|
1.001
|
[0.991 - 1.012]
|
Travel restrictions (ref: total travel ban)
|
|
|
|
|
At most screening and quarantine (lag 3 weeks)
|
1.052
|
[0.993 - 1.115]
|
--
|
--
|
Partial travel ban (lag 3 weeks)
|
1.027
|
[0.986 - 1.070]
|
--
|
--
|
No drop in passengers compared to 2019
|
--
|
--
|
1.016
|
[0.980 - 1.055]
|
Population mixing (interacted by calendar week)
|
|
|
|
|
Modified Stringency Index (lag 3 weeks)
|
0.977***
|
[0.974 - 0.987]
|
0.975***
|
[0.968 - 0.981]
|
Structural predispositions
|
|
|
|
|
Population size (residents)
|
1.878***
|
[1.627-2.168]
|
1.882***
|
[1.631 - 2.171]
|
Real GDP pc PPP
|
1.329***
|
[1.135 - 1.556]
|
1.313***
|
[1.122 - 1.538]
|
Population density (pop per km²)
|
0.991
|
[0.832 - 1.18]
|
0.985
|
[0.828 - 1.172]
|
Hospital beds per 1000 residents
|
0.637
|
[0.348 - 1.166]
|
0.628
|
[0.344 - 1.146]
|
Share of population aged 65+
|
1.744***
|
[1.346 - 2.26]
|
1.674**
|
[1.292 - 2.170]
|
Cardiovascular age-standardized death rate
|
0.765
|
[0.431 - 1.357]
|
0.946
|
[0.54 - 1.657]
|
Cancer age-standardized death rate
|
1.176
|
[0.760 - 1.820]
|
1.177
|
[0.752 - 1.842]
|
Prevalence of adult obesity
|
1.388
|
[0.656 - 2.938]
|
1.412
|
[0.644 - 3.100]
|
Recursive effects
|
|
|
|
|
Week of first COVID death
|
0.588 ***
|
[0.531 - 0.651]
|
0.587***
|
[0.644 - 0.650]
|
Observations
|
|
4,970
|
|
4,970
|
Random Intercepts: Variance (std)
|
|
|
|
|
Nuts / Country / Continent
|
|
2.923 (1.710)
|
|
2.770 (1.664)
|
Country / Continent
|
|
3.102 (1.761)
|
|
3.888 (1.971)
|
Continent
|
|
6.315 (2.513)
|
|
1.278 (11.303)
|
Random Slopes: Variance (std)
|
|
|
|
|
Nuts / Country / Continent
|
|
0.007 (0.086)
|
|
0.007 (0.085)
|
Country / Continent
|
|
0.021 (0.146)
|
|
0.021 (0.145)
|
Continent
|
|
0.022 (0.151)
|
|
0.722 (0.850)
|
Log Likelihood
|
|
-16,719.050
|
|
-16,743.310
|
Akaike Inf. Crit.
|
|
33,508.090
|
|
33,551.620
|
Bayesian Inf. Crit.
|
|
33,753.990
|
|
33,767.490
|
* p < 0.05,*** p < 0.01,**** p < 0.001. Note: the model controls for the main effects of all listed variables. Source: Sub-National Covid-19 Incidence and Determinants Dataset |
Weekly incidence: Modeling months separately
This section presents an alternative modeling strategy. Our dependent variable is still weekly incidence (Covid-19 deaths per capita), but we concentrate on the first weeks of April, May, June and July 2020. With this analysis we aim to check the robustness of the multilevel models discussed in the previous section, and at the same time dig deeper into the time-specific impact of different factors in the first four months of the pandemic, thus testing H3 (the ‘travel as spark’ hypothesis).
Table 3 reports the full models for the first week of each month. The severity of Covid-19 mortality in receiving areas was not significantly impacted by the number of incoming passengers, nor the risk associated with them coming from areas with higher incidence (measured through an interaction factor), nor airport centrality. Travel restrictions, however, did matter. Up until June, areas which imposed at best mild restrictions (such as screening passengers or quarantining them) fared worse than areas that took a stricter approach. Partial travel bans proved as effective as full bans in the early months (or even more in May), but significantly less in June and July, possibly because route triangulation bypassed them over time.
As for controls, their effect is not always as expected or significant. In particular, the lagged Oxford stringency index (modified in order to exclude travel policies, which we considered separately, like in the prior section) was found to associate positively with the severity of the pandemic. A possible explanation is reverse causality: stricter measures against population mixing had to be taken in areas with rising Covid-19 mortality. This is still the case with the 3-week lagged version of this variable, which is included in the model. More hospital beds are associated with a lower Covid-19 mortality at all stages, while GDP per capita and population size boost incidence. The precocity of Covid-19 outbreaks is a constantly significant predictor of higher mortality at all points of our series. Interestingly, with these controls, areas with a larger senior population have slightly (but not significantly) lower incidence in the observed months. Population obesity is a strongly significant predictor in April, but its effect declines and even reverts at the end of our period of observation. In contrast, areas with higher levels of cardiovascular mortality face a progressively increasing incidence.
Table 3
The impact of air passenger traffic on the severity of Covid-19 mortality in April-July 2020, controlling for population mixing (NPI), structural predispositions and recursive effects. Negative binomial regressions with continent fixed effects. Risk ratios and confidence intervals (in parenthesis)
|
(1)
|
(2)
|
(3)
|
(4)
|
|
April
|
May
|
June
|
July
|
Long-distance mobility
|
|
|
|
|
Inbound passengers per capita (lag 1 month)
|
0.831
|
0.863
|
0.767
|
0.764
|
|
[0.506-1.364]
|
[0.529-1.407]
|
[0.514-1.145]
|
[0.523-1.115]
|
Incidence in sending regions (lag 1 month)
|
1.388
|
1.194
|
1.472**
|
0.849
|
|
[0.715-2.693]
|
[0.885-1.612]
|
[1.161-1.866]
|
[0.685-1.052]
|
Incidence in sending x inbound passengers
|
0.964
|
1.089
|
0.884
|
1.031
|
|
[0.361-2.576]
|
[0.686-1.727]
|
[0.643-1.215]
|
[0.824-1.290]
|
Airport centrality (lag 1 month)
|
1.032
|
0.950
|
0.949
|
0.894
|
|
[0.897-1.188]
|
[0.836-1.079]
|
[0.809-1.114]
|
[0.762-1.050]
|
Air travel restrictions (ref: total travel ban)
|
|
|
|
|
At most screening and quarantining (lag 3 weeks)
|
1.997
|
1.631
|
2.792**
|
0.646
|
|
[0.860-4.639]
|
[0.712-3.735]
|
[1.318-5.914]
|
[0.229-1.825]
|
Partial travel ban (lag 3 weeks)
|
0.921
|
0.616**
|
1.612*
|
2.249**
|
|
[0.625-1.356]
|
[0.427-0.887]
|
[1.035-2.509]
|
[1.376-3.674]
|
Population mixing (NPI)
|
|
|
|
|
Modified stringency index (lag 3 weeks)
|
2.125***
|
1.329**
|
1.217
|
2.262***
|
|
[1.600-2.822]
|
[1.118-1.581]
|
[0.968-1.532]
|
[1.757-2.913]
|
Structural predispositions
|
|
|
|
|
Population size (residents)
|
27.79***
|
36.55***
|
48.72***
|
22.81***
|
|
[13.05-59.15]
|
[17.64-75.72]
|
[19.52-121.6]
|
[9.987-52.08]
|
Real GDP pc PPP
|
1.657***
|
1.761***
|
1.448*
|
1.437**
|
|
[1.265-2.171]
|
[1.354-2.290]
|
[1.076-1.948]
|
[1.108-1.864]
|
Population density (pop per km²)
|
1.062
|
1.043
|
1.556***
|
1.160
|
|
[0.869-1.299]
|
[0.861-1.263]
|
[1.217-1.989]
|
[0.933-1.442]
|
Hospital beds per 1000
|
0.554***
|
0.474***
|
0.438***
|
0.411***
|
|
[0.399-0.769]
|
[0.347-0.648]
|
[0.301-0.638]
|
[0.275-0.615]
|
Share of 65+
|
0.780
|
1.267
|
0.864
|
0.936
|
|
[0.567-1.073]
|
[0.986-1.628]
|
[0.624-1.195]
|
[0.714-1.227]
|
Cardiovascular age-standardized death rate
|
0.532***
|
0.839
|
1.681*
|
3.969***
|
|
[0.379-0.746]
|
[0.605-1.164]
|
[1.066-2.652]
|
[2.590-6.081]
|
Cancer age-standardized death rate
|
0.905
|
0.966
|
1.077
|
1.152
|
|
[0.669-1.224]
|
[0.733-1.274]
|
[0.742-1.564]
|
[0.780-1.703]
|
Prevalence of adult obesity
|
3.307***
|
1.126
|
0.786
|
0.400*
|
|
[1.989-5.498]
|
[0.672-1.887]
|
[0.414-1.491]
|
[0.197-0.809]
|
Recursive effects
|
|
|
|
|
Week of first COVID death
|
0.135***
|
0.482***
|
0.405***
|
0.715***
|
|
[0.0752-0.242]
|
[0.379-0.613]
|
[0.307-0.533]
|
[0.611-0.836]
|
Lnalpha
|
1.253**
|
1.214*
|
1.708***
|
1.376***
|
|
[1.066-1.472]
|
[1.045-1.409]
|
[1.452-2.010]
|
[1.152-1.644]
|
Observations
|
361
|
361
|
361
|
361
|
Pseudo R-squared
|
0.177
|
0.131
|
0.107
|
0.138
|
Continent Fixed Effects
|
Yes
|
Yes
|
Yes
|
Yes
|
* p < 0.05, ** p < 0.01, *** p < 0.001 |
Source: Sub-National Covid-19 Incidence and Determinants Dataset
|
Weekly incidence: modeling by epidemiological weeks
The previous analyses treated Covid-19 mortality by calendar weeks of the year 2020. In this section we repeat them aligning data by epidemiological week (i.e., the week since the first Covid-19 death) instead, thus following the evolution of the pandemic (Table 4), regardless of the different timing of the spread of the virus across the world in the first semester of 2020. We focus on the eighteen weeks immediately after the outbreak, with a view to further test H3 on the changing impact of long-distance travel over time.
Table 4 reports results for epidemiological weeks 6, 10, 14 and 18. In line with previous studies38, our findings still show an impact of total travel bans, which peaks at week 14 only to decline at week 18. Partial travel bans (mostly from China) do equally buffer incidence at the very early stage but end up being almost as ineffective as milder passenger monitoring measures. As weeks pass by, it is total travel bans that more clearly reduce mortality. The interaction term between ‘passenger volume’ and ‘average incidence in the area of origin of passengers’ never reaches statistical significance. Notably, the effect does not disappear with the passing of time, as posited by H3, and as we might expect given the progressively increasing controls on travelers’ health at departure. Its magnitude remains relatively modest, while other factors – notably, the precocity of the outbreak – are way more effective in explaining differences in the mortality of Covid-19 at all stages. The effects of structural factors diverge somehow from what was shown in the analysis based on calendar weeks. Richer areas are significantly more at risk in the early weeks, but this effect is reverted over time, as these areas can deploy more resources to contrast the disease. With time, population density becomes an increasingly impactful risk factor. The prevalence of obesity and a larger share of the senior population have a declining effect. This is quite surprising but it may reflect the fact that areas with a higher share of overweight individuals and people ages 65 or more were disproportionately hit in the earlier stages of the pandemic. The opposite trajectory is found for areas with higher cardiovascular death rates, which are subject to the rising severity of Covid-19. The number of hospital beds – a proxy of the solidity of the health system – is in fact a constant and rising predictor of lower mortality.
Table 4
The impact of air passenger traffic on the severity of Covid-19 mortality in the first 18 epidemiological weeks (week1=first Covid-19 death) controlling for population mixing (NPI), structural predispositions and recursive effects. Negative binomial regressions with continent fixed effects. Risk ratios and confidence intervals (in parenthesis)
|
(1)
|
(2)
|
(3)
|
(4)
|
|
Epi wk 6
|
Epi wk 10
|
Epi wk 14
|
Epi wk 18
|
Long-distance mobility
|
|
|
|
|
Inbound passengers per capita (lag 1 month)
|
0.982
|
0.748
|
0.807
|
0.579
|
|
[0.777-1.242]
|
[0.511-1.094]
|
[0.606-1.075]
|
[0.335-1.000]
|
Incidence in sending regions (lag 1 month)
|
1.269
|
0.951
|
1.119
|
1.033
|
|
[0.946-1.701]
|
[0.771-1.172]
|
[0.885-1.415]
|
[0.765-1.395]
|
Incidence in sending x inbound passengers
|
0.915
|
1.061
|
1.044
|
1.426
|
|
[0.662-1.264]
|
[0.781-1.441]
|
[0.888-1.227]
|
[0.972-2.094]
|
Airport centrality (lag 1 month)
|
1.105
|
1.043
|
0.898
|
0.984
|
|
[0.873-1.398]
|
[0.868-1.253]
|
[0.775-1.041]
|
[0.872-1.111]
|
Air travel restrictions (ref: total travel ban)
|
|
|
|
|
At most screening and quarantining (lag 3 weeks)
|
2.280
|
1.502
|
4.556***
|
1.256
|
|
[0.955-5.443]
|
[0.633-3.563]
|
[2.010-10.33]
|
[0.589-2.678]
|
Partial travel ban (lag 3 weeks)
|
0.525**
|
1.076
|
2.988***
|
1.621
|
|
[0.352-0.782]
|
[0.746-1.552]
|
[1.720-5.192]
|
[0.968-2.714]
|
Population mixing (NPI)
|
|
|
|
|
Modified stringency index (lag 3 weeks)
|
0.921
|
1.289*
|
1.038
|
1.016
|
|
[0.771-1.100]
|
[1.052-1.578]
|
[0.981-1.098]
|
[0.988-1.045]
|
Structural predispositions
|
|
|
|
|
Population size (residents)
|
12.89***
|
23.24***
|
18.48***
|
17.78***
|
|
[5.822-28.53]
|
[10.74-50.25]
|
[7.821-43.67]
|
[7.093-44.56]
|
Real GDP pc PPP
|
1.289*
|
1.237
|
0.987
|
0.734**
|
|
[1.023-1.624]
|
[0.999-1.530]
|
[0.791-1.232]
|
[0.585-0.921]
|
Population density (pop per km²)
|
1.066
|
1.304*
|
1.391**
|
1.389**
|
|
[0.854-1.332]
|
[1.044-1.627]
|
[1.094-1.769]
|
[1.084-1.779]
|
Hospital beds per 1000
|
0.567**
|
0.433***
|
0.375***
|
0.313***
|
|
[0.400-0.802]
|
[0.308-0.609]
|
[0.243-0.578]
|
[0.196-0.500]
|
Share of 65+
|
0.835
|
0.864
|
0.553***
|
0.649**
|
|
[0.627-1.111]
|
[0.660-1.131]
|
[0.418-0.730]
|
[0.485-0.869]
|
Cardiovascular age-standardized death rate
|
1.425
|
1.680**
|
4.899***
|
2.938***
|
|
[0.928-2.186]
|
[1.150-2.454]
|
[3.111-7.714]
|
[1.995-4.326]
|
Cancer age-standardized death rate
|
0.859
|
0.850
|
1.310
|
1.027
|
|
[0.618-1.195]
|
[0.611-1.183]
|
[0.855-2.008]
|
[0.643-1.643]
|
Prevalence of adult obesity
|
1.538
|
1.085
|
0.389**
|
0.844
|
|
[0.884-2.676]
|
[0.629-1.870]
|
[0.199-0.760]
|
[0.399-1.783]
|
Recursive effects
|
|
|
|
|
Week of first COVID death
|
0.517***
|
0.675***
|
0.617***
|
0.720***
|
|
[0.445-0.601]
|
[0.590-0.771]
|
[0.537-0.709]
|
[0.618-0.839]
|
Lnalpha
|
1.409***
|
1.250**
|
1.502***
|
1.526***
|
|
[1.211-1.639]
|
[1.067-1.465]
|
[1.265-1.784]
|
[1.274-1.828]
|
Observations
|
357
|
357
|
357
|
357
|
Pseudo R-squared
|
0.125
|
0.116
|
0.122
|
0.123
|
Continent Fixed Effects
|
Yes
|
Yes
|
Yes
|
Yes
|
* p < 0.05, ** p < 0.01, *** p < 0.001. Source: Sub-National Covid-19 Incidence and Determinants Dataset |