Skip to main content

Cross-border mobility responses to COVID-19 in Europe: new evidence from facebook data

Abstract

Background

Assessing the impact of government responses to Covid-19 is crucial to contain the pandemic and improve preparedness for future crises. We investigate here the impact of non-pharmaceutical interventions (NPIs) and infection threats on the daily evolution of cross-border movements of people during the Covid-19 pandemic. We use a unique database on Facebook users’ mobility, and rely on regression and machine learning models to identify the role of infection threats and containment policies. Permutation techniques allow us to compare the impact and predictive power of these two categories of variables.

Results

In contrast with studies on within-border mobility, our models point to a stronger importance of containment policies in explaining changes in cross-border traffic as compared with international travel bans and fears of being infected. The latter are proxied by the numbers of Covid-19 cases and deaths at destination. Although the ranking among coercive policies varies across modelling techniques, containment measures in the destination country (such as cancelling of events, restrictions on internal movements and public gatherings), and school closures in the origin country (influencing parental leaves) have the strongest impacts on cross-border movements.

Conclusion

While descriptive in nature, our findings have policy-relevant implications. Cross-border movements of people predominantly consist of labor commuting flows and business travels. These economic and essential flows are marginally influenced by the fear of infection and international travel bans. They are mostly governed by the stringency of internal containment policies and the ability to travel.

Background

Covid-19 is a disease induced by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). The first known cases of Covid-19 were registered in December 2019 in Wuhan (China), but the virus rapidly spread into many countries, leading to a global pandemic that has held the world captive for many months. As until recently no effective vaccine or medication was available, government policies were mainly targeted towards tracing and disrupting infection chains. Many countries introduced coercive measures and disincentives to limit within- and cross-border mobility of people in the hope of reducing the virus propagation. Whether these policies were effective or not remains an open question, and there is much to be gained from better understanding the evolution and determinants of people’s mobility during Covid times.Footnote 1

There is strong evidence that within-border (or internal) people’s mobility declined during the Covid-19 crisis. Existing literature relies on big data provided by private cellular phone companies,Footnote 2 and documenting spatial movements in real time. [2] uses Vodafone data for Italy, Portugal and Spain, and finds that women and younger people show the largest drop in mobility. [3] combines telco data with household surveys to highlight a sharp decline in short-distance mobility, as proxied by daily time spent at parks, retail and recreation, grocery, transit locations, and workplaces. Using mobility data from the analytics company SafeGraph, [4] finds that the decline in mobility in New York and in four other U.S. cities is mostly driven by the fear of infection, rather than by legal restrictions. Although it also finds a significant impact of non-pharmaceutical interventions (NPIs) such as closing nonessential businesses, sheltering in place, and school closures, the dominant role of infection threats is confirmed by [5], who relies on Google mobility data.

By contrast, there is scant evidence of the impact of Covid-19 on cross-border (or international) mobility, which is due to the absence of high-frequency data on border crossings.Footnote 3 In this paper we aim to fill this research gap by addressing the following research questions: (i) How has the Covid-19 impacted cross-border movements of people? (ii) Are these changes due to coercive measures (such as containment policies or international travel bans) or by the fear of contracting the virus?

Understanding the determinants of cross-border mobility responses to Covid-19 is important for economic and epidemiological reasons. Cross-border movements of people predominantly consists of labor commuting flows and business travels. Economically speaking, labor mobility is a key ingredient for growth and competitiveness in normal times. And in a pandemic context, restrictions placed on how workers move around can slow down economic recovery prospects, by making it more difficult for businesses to hire productive workers. They can also induce severe economic impacts on cross-border workers and their families. Epidemiologically speaking, the role that mobility is playing in the spread of the disease is still unclear. Using SafeGraph data for New York city and for other U.S. cities, [6] find that (internal) mobility increased the spread of the disease in the early stage of the pandemic. [7] also finds a strong association between internal movements and mortality using Google mobility data for the UK. In the same vein, [8] shows that (internal and international) travel bans enacted during the Chinese Lunar New Year holiday helped reduce the spread of the virus, and [9] argue that an appropriate coordination would considerably improve the likelihood of eliminating community transmission throughout Europe. By contrast, others expressed skepticism about the epidemiological consequences of travel bans, arguing that the impacts of these restrictions are not well understood [10] or poorly effective [1115] once patient zero has already spread the virus across regions.

Without taking any position on the fact that cross-border mobility should be limited or encouraged, we use a unique database on daily mobility of European Facebook users to shed light on the evolution of cross-border movements of people during an entire pandemic year, and to compare the effects of coercive policies with those related to the fear of infection. Our results contrast with those obtained for internal mobility. The following sections successively describe our data sources, methods and findings.

Data

Border crossings

Data on daily cross-border mobility is obtained from Facebook (denoted by FB, henceforth) for the period from the 29th of February 2020 to the 28th of February 2021 [16]. The database documents cross-border flows of FB users with location services enabled, who travel from an origin to a destination country by any means of transportation (car, train, air, etc.) during each 24-hour time periodFootnote 4. Only daily flows with a minimum of 1,000 movers are reported in the dataset. All flows below 1,000 are set to zero in order to minimize re-identification risk of FB users. This means that our outcome variable is left-censored. To limit the impact of censoring, we focus on 45 country pairs (involving 30 contiguous European countries) characterized by at least 25% of uncensored values of daily traffic during the period of observation. This selection limits the ability to generalize our results but is necessary to limit the impact of censoring, and allow smooth estimation with Machine Learning methods. Further limiting the impact of censoring, we use the 7-day centered rolling average of daily flows.

Although FB data has high coverage, FB users are not a random sample of the population. This raises concerns about representativeness. Figure A.1 in Appendix mitigates these concerns by showing a strong association between daily movements of FB users and the (estimated) number of daily border crossings in the pre-Covid-19 period, which are presented in Appendix Table A.1.Footnote 5 However [19] shows that FB users are over-represented in the population aged 20-40, with a high level of education and an above-average income level. It means that groups of individuals under 20 or over 65 and those with lower income/education levels are under-represented in FB data. Although people under 20 and over 65 form a tiny minority of the population under investigation (i.e., commuters and business travelers), this is a limitation of our work. FB users are more likely to belong to the richest and healthiest parts of the movers’ population.

Let us denote by \(M_{i\rightarrow jt}^{F}\) the count of FB movers from country i to country j at day t. When focusing on contiguous countries (i.e., the pairs of countries that exhibit the largest numbers of daily cross-border movements by far, and that are the least affected by censoring rules), the number of movers from i to j is almost identical to the number of movers from j to i at each day t (i.e., \(M_{i\rightarrow jt}^{F} \simeq M_{j\rightarrow it}^{F}\)). The reason is that border crossings predominantly consist of back-and-forth movements of commuting workers and business travelers, who move for short periods and for economic reasonsFootnote 6. This is also the case in the summer vacation period when considering a 7-day centered rolling average of daily flows. This means that \(M_{i\rightarrow jt}^{F}\) and \(M_{j\rightarrow it}^{F}\) are reflecting the same reality, and say nothing about the primary direction of the flows. For this reason, we define the level of bilateral traffic of FB users between countries i and j as:

$$ T_{ijt}^{F} \equiv \text{Max}\left[ M_{i\rightarrow jt}^{F},M_{j\rightarrow it}^{F}\right], $$
(1)

and see it as a proxy for the scaled sum of the two unobserved unidirectional flows between the two countries, ϕ(Mijt+Mjit), where the scale factor ϕ denotes the fraction of FB users in the actual number of movers (denoted by Mijt and Mjit). Using the maximum as weighting scheme has the advantage to limit the impact of censoring in the case movers in one of the two directions are below the threshold. As \(T_{ijt}^{F}=T_{jit}^{F}\), we can get rid of the dyadic dimension of the data, treat each country pair as a one-dimensional observation, and divide the size of the sample by two. In the methodological section, however, we explain how priors about the primary direction of the flows can be used to improve the quality of fit of our models.

To avoid dealing with re-scaling issues, we express traffic counts as relative deviations from their initial or pre-Covid-19 levels – in our case, the levels observed at the outset of the pandemic (denoted by day 0). We thus use the relative deviation in bilateral traffic between day 0 and day t, \(\tau _{ijt} \equiv \frac {T_{ijt}-T_{ij0}}{T_{ij0}}\), as a variable of interest instead of focusing on the level of traffic Tijt itself. Modelling relative deviations is also helpful to avoid over-fitting large corridors at the expense of small corridors, and mitigates representativeness issues even if the scale factor (ϕ) varies across country pairs.

Figure 1 portrays these relative deviations in the aggregate level of traffic between all country pairs included in our sample. The curve largely mirrors the three phases of the pandemic, depicting a stark drop in traffic in March 2020, a recovery during the spring and summer periods, and a new contraction in the post-summer period. Between end of February and early April 2020, the aggregate traffic level decreased from 720,000 to 130,000, implying a 82% drop. Aggregate traffic never fully recovered to the February levels in our period of observation. This also holds true during the summer vacation period when international travels were largely liberalized. The pace and strength of these changes vary across the three phases of the pandemic. The drop in March 2020 was strong and sudden, while the summer peak and the post-summer contraction were more gradual.

Fig. 1
figure 1

Aggregate Traffic Deviations from Pre-Covid Levels. Source: Facebook data on daily border crossings. Notes: Y-axis represents the average of τijt, the percentage change (times 100) in the 7-day moving average traffic compared to t=0 over all destinations. The weights are the traffic levels observed in pre-Covid-19 period (i.e., t=0)

Aggregate fluctuations mask large differences across country pairs. Bilateral traffic returned to its pre-Covid level in a minority of cases. For the majority of corridors, however, the traffic level has not fully recovered. This is illustrated in Fig. 2, which depicts the evolution of people’s traffic in corridors involving four open countries, namely Luxembourg, Switzerland, Italy, and Serbia. Luxembourg is the country with the highest share of cross-border workers in Europe. Given the economy’s high reliance on cross-border workers, the Luxembourg government has never implemented international travel restrictions during the pandemic. Luxembourg experienced a significant drop in traffic in March 2020, whatever the partner country. After one month of lockdown, traffic levels recovered pretty quickly until reaching a plateau at about -25% since June 2020. Switzerland is the country with the largest number of cross-border commuters in Europe. This country experienced a larger drop during the first lockdown, and a slower recovery. Furthermore, the variability across corridors is considerably greater than in Luxembourg.

Fig. 2
figure 2

Traffic Deviations from Pre-Covid Levels for Selected Corridors. Source: Facebook data on daily border crossings. Notes: Y-axis represents τijt, i.e. percentage change (times 100) in corridor 7-day moving average traffic compared to t=0 in corridor ij

Italy has been severely impacted by the pandemic, and responded with national and international travel bans. We observe similar patterns of contraction and recovery during the first two quarters of 2020, followed by a substantial increase in traffic during the holiday summer period, and a second lockdown-type contraction in the post-Summer period. Finally, the patterns observed in Serbia are less conclusive as they are more severely affected by censoring rules. Serbia is an important origin and transit country for migrants and refugees entering the EU. Overall these patterns illustrate the need to account for corridor-specific heterogeneity when analyzing the determinants of bilateral traffic. Variations are likely to be influenced by seasonal effects, epidemiological risks, and policy measures implemented in the countries. We now turn to the description of the data sources used to proxy epidemiological conditions and the stringency of national policies.

Explanatory features

We link variations in cross-border traffic during the pandemic to daily changes in epidemiological conditions and containment policies in the countries involved. We proxy the severity of the pandemic with the daily numbers of new Covid-19 cases and new Covid-19 related deaths in each country using data from [20]. With regard to containment measures, we use data on daily policy responses from the Oxford Covid-19 Government Response Tracker (OxCGRT) [21]. The latter database consists of 18 ordinal indicators capturing the levels of nonpharmaceutical interventions (NPIs). Based on our priors as to which policies likely affect mobility, we choose to include the eight mobility-related measures that form the “containment and closure policies” block (denoted by C1-C8 in the database) as well as proxies for the intensity of testing and contact tracing (denoted by H2 and H3). We rescale all sets of predictors between 0 and 1, and align them with the definition of the outcome variable using the centered 7-day rolling average at each day.

For all features and days, Fig. 3 displays the cross-country mean level of each containment index over time. Values close to one represent higher Covid-19 cases/deaths or more stringent responses. As containment policies were implemented in most countries during the second half of March 2020, maximum values are observed during this period. Testing and contract tracing were implemented more heterogeneously across countries and peaked in the summer of 2020. Reported number of new Covid-19 related cases/deaths were much greater during the second wave and peaked at the end of the year 2020.Footnote 7

Fig. 3
figure 3

Evolution of the average value of Covid cases and government policy measures. Source: Oxford Covid-19 Government Response Tracker (OxCGRT). Note: The values of each indicator is scaled between 0 and 1, and the average is computed using the 30 countries included in our sample

Figure 4 shows the cross-country correlations between each of the explanatory variables and the relative deviations in average traffic (i.e. deviation of the country-specific mean level of traffic with all potential partner countries in the sample). Containment policies are positively correlated with each other, and moderately correlated with epidemiological conditions, which allows us to include both sets of variables jointly in our regression and machine learning models. However, the fact that containment policies are correlated with each other raises concerns of multicollinearity, and motivates the usage of a limited number of synthetic policy indices. These indices are obtained by conducting a Principal Component Analysis (PCA) of all policy measures (C1-C8, and H2-H3) over the entire sample of observations, and by extracting the first two components. The first component mainly represents the C1-C8 measures which are strongly correlated with each other while the second component corresponds to a higher variance for the H2-H3 measures. We will compare the results obtained when using a comprehensive specification, including all (collinear) explanatory variables, with those obtained when using a parsimonious specification, including the two synthetic PCA components as predictors.

Fig. 4
figure 4

Correlation Matrix of the different variables. Source: Own computations. Notes: Unilateral traffic growth for each country i is the relative deviation in aggregate traffic involving country i, \(\sum _{j=1}^{n} T_{ijt}\), as compared to the pre-Covid-19 period (t=0)

Turning to the main correlations of interest (i.e., correlations between the relative deviations in cross-border traffic and each predictor), the figure displays a moderate negative correlation with containment policies, as well as low and positive correlations with testing and tracing policies. The strongest associations are obtained for cancellations of public events, restrictions on public gatherings, requirements to stay at home, and restrictions on internal movements. Two interesting observations arise from these partial correlations. First, the negative correlation between measures of epidemiological intensity (Covid-19 cases and deaths) and changes in total traffic are rather low. This might suggest that, in contrast with studies on within-border mobility, the fear of being infected might play a less important role in explaining changes in cross-border movements. Second, among containment policies, the implementation of international travel bans is far from being the most strongly correlated covariate. This suggests that travel bans might be effective in limiting non-essential travels, but less effective to limit labor commuting flows and business travels, which represent the overwhelming majority of daily border crossings between European countries (see Table A.1). By contrast, cancellations of events and constraints on internal movements and gatherings are highly correlated with traffic variations, possibly because such constraints better proxy changes in economic activity and incentives to move. In the following section, we describe the regression and machine learning methods used to test these hypotheses using bilateral traffic growth as a dependent variable.

Methods

Mobility patterns identified in the previous section might result from various factors such as travel bans, sanitary measures influencing economic costs and incentives to move for work and business (e.g. sectoral lockdown, work-from-home practices) or for leisure (e.g. shops, restaurant and bar closings), or the fear of the virus itself. Our goal here is to identify the determinants of the relative deviation in daily traffic of people between country i and country j (τijt), considering all NPIs and epidemiological daily indicators (xit and xjt) during the Covid-19 crisis. Our models are also used to predict the effects of epidemiological restrictions, NPIs and mobility restrictions on traffic counts.

We combine two analytical methods, Econometric Modelling (EM) and Machine Learning (ML). EM and ML techniques are generally used for different purposes. In EM, gravity models are used to explain human mobility flows between two countries.Footnote 8 EM models require imposing one analytical specification for the response function, which governs the derivatives of the dependent with respect to covariates. In the gravity specification, a large set of fixed effects are included to capture relevant confounders, and allow for identification of potential causal effects based on the so-called within variation. ML techniques are at the other extreme of the bias-variance tradeoff. They do not require strong analytical assumptions and allow, by design, to explore a larger set of regression functions including linear or polynomial combinations of the covariates. This increase of the so-called model capacity comes with two drawbacks. First, the models are more complex and are usually not easy to interpret. Contrary to EM, ML techniques might computationally suffer from the inclusion of large sets of control fixed effects. Second, there is always a risk of overfitting the training data and the identification of causation links is usually not an objective per se. Despite these difficulties, we use both EM and ML to identify converging messages, and the triangulation of results can serve to strengthen our evidence base.

We use four models exploring a broad range of learning techniques: (i) A gravity model based on the linear regression method [23]; (ii) A K-nearest neighbors method (KNN), which predicts the dependent variable by interpolation of its nearest observation neighbors in the training set [24]; (iii) A Gradient Boosting method (GBoost), whose predictions are based on a set of decision tree models [25]; (iv) A Multi-layer Perceptron (MLP), which is a classic neural network approach [26, 27]. The last three models rely on different ML regressors, each based on a distinct type of technique. We assess the predictive performance of each model using the very same (and standard) cross-validation ML methodology. The goal of the study is not to design a forecasting model, but rather to identify the main determinants of mobility, and to investigate whether these different approaches generate converging findings. Therefore, instead of the validating our model on a particular sub-period (as is usually done to evaluate a time-series model), the observations composing the cross-validation folds are randomly chosen within the full sample. All models are implemented in Python via the Skicit-learn library [28].

Approaches with or without directional priors

Ideally, mobility models aim to characterize the evolution of the unidirectional flow of people (Mijt) from an origin country i to a destination country j at day t, or of their relative deviation from the initial reference period \(\left (\mu _{i\rightarrow jt} \equiv \frac {M_{i\rightarrow jt}-M_{i \rightarrow j0}}{M_{i\rightarrow j0}}\right)\), based on a set of features available for the same time period. Without loss of generality, the general functional form fM of such a model can be written as:

$$ \mu_{i\rightarrow jt} = f_{M}(x_{it},x_{jt}, d_{ij}, d_{t}) + \eta_{i\rightarrow jt} $$
(2)

where xit represents the set of origin-specific determinant, xjt is a set of destination-specific determinants, dij is a set of bilateral dummies capturing time-invariant bilateral resistance (including initial Mij0, distance, language proximity, cultural proximity, etc.), dt a set of day dummies capturing weekdays and seasonal trends (e.g. holiday season, general feeling of risk when traveling, etc.), and ηijt an error term. In our case, the vectors of explanatory variables xit and xjt capture the set of NPIs and epidemiological variables, and should also be interpreted as variations from period 0 since xi0 and xj0 are equal to zero in the pre-Covid-19 period.

With FB data, the primary direction of the cross-border flows is unknown, which implies that Mijt and Mjit cannot be distinguished a priori. Instead, we observe the relative deviation in bilateral traffic, τijt, and we have to estimate the function fT linking bilateral traffic to the set of explanatory features without being able to distinguish between origin- and destination-specific determinants. A learning approach without directional priors writes as:

$$ \tau_{ijt} = f_{T}(x_{it},x_{jt}, d_{ij}, d_{t}) + \eta_{ijt} $$
(3)

It is possible, however, to discipline the model with priors about the direction of the flows. As bilateral traffic is a proxy for the sum of unidirectional flows (Tijtϕ(Mijt+Mjit)), relative deviations in Tijt can be expressed as a weighted sum of the relative deviations in unidirectional flows: τijt=ωij,0×μijt+ωji,0×μjit, where ωij,0=1−ωji,0 is the pre-Covid-19 share of unidirectional cross-border flows from country i to country j in total traffic between the two countries. Estimates for ωij,0 are constructed using pre-Covid-19 data on commuters, air travels and international migration, and then used as priors to discipline the model (these shares are depicted in Figure A.2 in Appendix).

We can thus create two sets of weighted features, namely \(X_{ijt}^{o}\) for origin-specific effects, and \(X_{ijt}^{d}\) for destination-specific effects, defined as follows:

$$\begin{array}{*{20}l} X_{ijt}^{o} = \omega_{i\rightarrow j0}x_{it}+\omega_{j\rightarrow i0}x_{jt}\\ X_{ijt}^{d} = \omega_{i\rightarrow j0}x_{jt}+\omega_{j\rightarrow i0}x_{it}. \end{array} $$

The model with directional priors is obtained after replacing (xit,xjt) in Eq. (3) by \(\left (X_{ijt}^{o},X_{ijt}^{d}\right)\). It writes as:

$$ \tau_{ijt} = f_{T}\left(X_{ijt}^{o},X_{ijt}^{d},d_{ij}, d_{t}\right) + \eta^{od}_{ijt}. $$
(4)

If the true model for unidirectional flow (fM(.)) was linear, plugging weighted covariates in the estimated model for bilateral traffic (fT(.)) would allow retrieving the true origin- and destination-specific coefficients of interest accurately. Although this is not the case when fM(.) is non linear, using weighted covariates might improve the quality of fit or facilitate the interpretation of the results. The rationale is that the effect of policies depends on where they are implemented and on the primary direction of the flows. Suppose ωij01 (i.e., flows mostly go from i to j), then an increase in restrictions/stringency at destination (resp. at origin) makes \(X_{ijt}^{d}\) positive (resp. \(X_{ijt}^{o}\) positive) and is more (resp. less) likely to reduce the flow of cross-border movements. Using directional priors allows approximating origin- and destination-specific effects without observing the direction of the flows during the pandemic year.

Permutation feature importance

To identify the main causes of daily mobility variations, the importance of each feature is computed for the different regression approaches with permutation feature importance [29]. It has a the advantage of working similarly for all regression models considering them as a black-box models [30]. It is defined to be the decrease in the regression score when a single feature value is randomly shuffled across observations. More exactly, for a given model f, it first calculates a baseline score Sf provided by f when it is fitted, and then evaluated with a certain metric on the whole sample. Then for each possible feature x the modified score \(S_{f,x}^{*}\) is computed by evaluating f on the transformed data set where the values of feature x are randomly permuted across all observations. The mean importance of the feature x for the model f is computed as:

$$ I_{f,x} = \frac{1}{K} \sum_{k=1}^{K} \frac{S_{f}-S_{f,x,k}^{*}}{S_{f}} $$
(5)

where K is the number of random permutations realized for each feature, and \(S_{f,x,k}^{*}\) is the \(S_{f,x}^{*}\) score for the kth permutation.

In order to compare the importance values of the different models in an equivalent manner, the values If,x are scaled between 0% and 100% separately for each model f. The mean features importance If,x are computed over 10 permutations using the negative mean absolute error (MAE) and the Root Mean Squared Error (RMSE). This means that the bigger If,x, the more permutations of the feature x degrades the quality of predictions for the model f and the feature is considered as more importantly associated with the target variable.

Results

We present our results in two steps. First, we assess the predictive power of the various models. This implies comparing learning methods with or without directional priors and with or without day/corridor control dummies. We compare their predictive performance by using out-of-sample predictions and computing the MAE and RMSE. Second, we use the estimated models to rank the importance of different features relying on permutation techniques. Using multiple models allows assessing the robustness of our findings.

Validation of models

We first investigate whether adding priors about the direction of the flows and/or adding a full bunch of day and corridor dummies improves the performance of our learning models. Models without directional priors are described in Eq. (3), while models with priors use weighted regressors, as described in Eq. (4). A 10-fold cross-validation over the whole data set (from the 29th February 2020 until the 28th February 2021) is realized for each model to assess its performance. Table 1 reports the MAE, RMSE and their standard error across cross-validation folds obtained with different learning methods.

Table 1 MAE comparison of the different models with or without directional priors and dummies

It shows that directional priors (ωij,0 and ωji,0) do not bring significant additional predictive power under most learning approaches when day and corridor dummies are not factored in (see Panels A and B). The only exception is the G-Boost method. On the contrary, when day and corridor dummies are included (Panels C and D), adding directional priors slightly improve the quality of fit with virtually all learning techniques (except with KNN). In addition, we show below that distinguishing between origin- and destination-specific effects as in Eq. (4) makes the interpretation of the results much easier. Therefore, the model with directional priors will be prioritized in the rest of the analysis.

Second, we investigate whether the inclusion of day-specific effects – i.e., 366 time dummies, dt, that are common to all corridors and capture unobserved variations such as seasonal changes, synchronized fears of infection, etc. – and corridor-specific effects – i.e., 45 corridor dummies, dij, that are time invariant and capture unobserved variations such as the skill level of the cross-border workforce, linguistic and cultural proximity between countries, etc. – improves the predictive power of our models. Again, we perform another 10-fold cross-validation on different versions of each model. Panels C and D in Table 1 includes both sets of dummies jointly, with or without directional priors. In the absence of directional priors, the inclusion of day and corridor dummies reduces the MAE and RMSE by 20 to 30% whatever the learning technique used. When directional priors are factored in, the dummies improve the performance of the linear and MLP models, whereas they deteriorate the quality of fit under the KNN and G-Boost models. This is because adding day and corridors dummies drastically increases the number of parameters to be estimated, and some ML approaches (like KNN) are known to suffer from the curse of dimensionality.

To further explore this issue, Table A.2 in Appendix considers the model with directional priors and adds one set of dummies at a time. The inclusion of 45 corridor-specific dummies always improves the quality of fit. On the contrary, the inclusion of 366 day-specific dummies deteriorates the performance of KNN and G-Boost methods. This confirms that the gains from adding information about unobserved common time trends, which might already be captured by the relatively well synchronized trends in observed epidemiological conditions and containment measures, is outbalanced by the costs linked to the inflated dimensionality of the computation problem.

Third, ML techniques always outperform the linear EM model. This result was also expected given that ML is based on more complex prediction methods that allow for non-linear relationships between variables, and account for non-stationary variations contained in the matrices of \(X_{ijt}^{o}\) and \(X_{ijt}^{d}\). The KNN always produces the best quality of fit. The error of this approach is minimal when the number of neighbors k used to estimate the relative deviations in traffic is low (say, 2 or 3). Its impressive performance in 10-fold cross-validation can be explained by the fact that the model finds a small number of observations for which the relative deviations in traffic are similar to those that must be predicted. In general, the closest neighbors are observations of the days preceding or following the daily level of traffic observed in the same corridor.

Main sources of variations in cross-border mobility in COVID times

In order to identify the government policy indicators having the greatest impact on relative deviations in traffic, the importance of each feature is computed for the different approaches involving directional priors and dummies. Directional priors allow us to distinguish the effects of origin-specific features from those of destination-specific features. Table 2 presents the results from this exercise. Features are ranked by decreasing order on the basis of the average predictive power across the four learning techniques. The column ‘Avg.’ gives the mean value of error metric averaged over the four models. Panel A provides the results obtained with the saturated models including the large set of corridor and day dummies, which is the first-best model when using linear and MLP learning techniques. Panel B gives the results obtained with corridor dummies only, which is the first-best model when using the KNN and G-Boost techniques. Results of Panel B will be discussed in the next section. In the top part of the table, the models use all individual features depicted in Fig. 3, despite the high level of correlation between some of them. In the bottom part of the table, the models use synthetic containment features derived from a PCA analysis.

Table 2 Feature ranking by origin and destination

When considering all individual features, the ranking based on their predictive power varies across models. We identify, however, several common and interesting findings. First, school closures in the origin country has the largest average impact on the variation in daily traffic. Remember that cross-border traffic predominantly consists of labor commuting flows and business travels. School closures at origin imply that many parent workers are forced to take parental leave and cannot commute to work. In the same vein, school closures in the destination country are also paralyzing economic activity in the destination country and reduce incentives to move. This result is in line with [31], who finds that school closures were among the most important predictors of internal mobility in March 2020. Second, variations in traffic are mostly impacted by containment measures. Based on the average predictive power (col. ‘Avg.’), ten out of the twelve most predictive features involve C-type containment measures implemented in the origin or destination countries. Third, the fear of being infected in the destination country, as proxied by the destination-specific number of Covid-19 deaths and cases (appearing in bold characters in Table 2) are among the least predictive features. Fourth, international travel bans in origin and destination countries (appearing in italics in Table 2) also have a low predictive power. We thus conclude that cross-border daily flows are marginally influenced by the fear of infection and international travel bans. They are mostly governed by the stringency of internal containment policies and by family constraints. It is worth noticing that models without directional priors deliver very similar results, as illustrated in Table A.3 in the Appendix.

In addition to school closures, the top panel of Table 2 suggests that the most important containment measures are the cancellation of public events, restrictions on gatherings, restrictions on internal movements and stay-home requirements. In addition, the twelve most predictive features include 5 destination-specific and 7 origin-specific measures. However, as illustrated in Fig. 4, these measures are highly correlated at the national level. Hence, instead of feeding the model with correlated features, the bottom part of the table uses synthetic indices of containment and sanitary measures. We use a PCA analysis to reduce the dimensionality of the origin- and destination-specific containment measures and we extract the first two components of the PCA.

Remember that the first PCA component can be interpreted as an average index of stringency of containment measures (i.e. C1-C8 indices); the second component captures testing and tracing policies (i.e. H2 and H3). The results clearly reveal that the stringency of containment measures in the destination country has, by far, the greatest predictive power. The average stringency of containment measures at origin is the second most predictive power, with an average importance equal to 40% of that of the destination country. This comforts the idea that cross-border daily flows of people mostly involve economic/essential movements which can only be influenced by changes in incentives to move or coercive mobility constraints. In line with the top part of the table, variables influencing the fear of infection have negligible impacts on border crossings.

Discussion

The EM and ML techniques used in this paper allow highlighting a strong association between the evolution of bilateral traffic between contiguous countries and containment policies in the destination country as well as school closures in the origin country. Association does not imply causation. It could be argued that this statistical association is governed by the influence of unobserved characteristics affecting both policy changes and mobility simultaneously, or that a reverse causation mechanism operates (i.e. cross-border mobility influences policy reforms). The fact that our results are robust to the inclusion of a large set of day and corridor dummies capturing unobserved time- and corridor-specific characteristics strongly mitigates the first misspecification concern.

With regard to reverse causation, concerns are mitigated by the use of high-frequency data. We cannot reject the possibility of a mobility-driven propagation of the virus requiring new containment measures. However, such a mechanism takes time to operate. Mobility shocks at day t do not generate immediate and visible epidemiological consequences, and policy responses are also implemented with a certain delay. By contrast, our estimates suggest that containment policies are immediately associated with changes in cross-border mobility. This prudently supports the existence of a causation effect of containment measures on mobility.

An opposite argument that goes against the reverse causation issue is that it also takes time for information about epidemiological conditions to be assimilated by potential movers. Hence, the fear of infection could be better proxied by the lagged numbers of Covid-19 cases and deaths. Our results are preserved and even reinforced when using lagged proxies for the fear of infection. More precisely, traffic at day t is very badly predicted by the number Covid-19 cases and deaths observed one or two weeks before the date (see Table A.4 in the Appendix). Hence, we can reasonably rule out that the low impact of infection fears at destination is driven by a misspecification problem.

In the same vein, it could be argued that fears are strongly synchronized across countries and captured by the day dummies. These concerns are mitigated by the fact that removing the day-specific dummies does not alter our conclusions. In Panel B of Table 2, day dummies are excluded. The only significant change is that the number of Covid cases in the origin country has a greater predictive power, possibly implying that people in sick leave self-isolate and stop moving. However, the predictive power of international travel bans (in italics) and fears to be infected at destination (epidemiological conditions at destination are in bold characters) remain low when using both individual and synthetic features. Again, the stringency of containment measures in the destination country has the greatest predictive power by far.

Table 3 Estimates of the average count of daily cross-border outflows in 2019
Fig. 5
figure 5

Daily traffic estimates (X-axis) and FB data (Y-axis) by corridor (1 March 2020). Sources: Numbers of daily commuters are extracted from Eurostat data by Nuts2 region in 2019; Numbers of air passengers are extracted from Eurostat monthly statistics on air passenger transport on March 1, 2020; Data on international migrants are extrapolated from [34] for the year 2015, assuming a conservative 50% growth in the flows between 2015 and 2020. FB data are the Facebook data on daily border crossings on March 1, 2020. Note: All variables are expressed in logs

Finally, one may fear that our results can be affected by the quality of data as well as their comparability across countries (see [32] for a discussion of the strengths and weaknesses of data sources on Covid-19). Data on Covid-19 threats and NPIs (such as the number of new cases and deaths) are likely to be subject to country and time biases due to multiple reasons (e.g. changes in testing strategies, seasonal effects, variations in reporting delays or in the classification of Covid-related deaths, etc.). It is worth emphasizing, however, that the Permutation Feature Importance technique that we use is, by construction, robust to noise.Footnote 9 In addition, the systematic use of bilateral and time dummies (dij and dt) allows us to mitigate the biases caused by unobserved variations in the quality of data across country pairs and periods.

Conclusion

Existing literature shows that people within-border mobility has drastically declined in times of Covid-19, primarily because of the fear to be infected in parts of the population. To the best of our knowledge, our study is the first to analyze the effect of Covid-19 and related containment measures on people’s cross-border movements. In line with the findings above, we also document a sharp decline in cross-border mobility in general, especially during the first lockdown and in the second and third waves of the pandemic. However, these variations in cross-border mobility are mostly induced by local containment policies in the destination country, and school closures in both countries. The fear of infection and international travel bans have little influence on cross-border movements.

The likely reason is that cross-border daily flows of people are predominantly made of commuting workers and business travelers who move for economic/essential reasons. These economic flows are observed between contiguous countries, and account for 99% of international movements of people when compared with the flows of migrants and refugees. Their magnitude varies with the economic costs and incentives of moving, which depend on lockdown measures and on the stringency of internal containment policies. In addition, international travel bans do not apply to commuters and businessmen. Although there is no consensus on the fact that these flows contribute to the propagation of the virus, policy-makers must be aware that economic movers hardly adapt their mobility decisions to epidemiological threats. Border crossings can only be controlled with internal coercive policies.

Appendix A: Appendix – figures and tables

Fig. 6
figure 6

Distribution of Weights. Sources: Numbers of daily commuters are extracted from Eurostat data by Nuts2 region in 2019; Numbers of air passengers are extracted from Eurostat monthly statistics on air passenger transport on March 1, 2020; Data on international migrants are extrapolated from [34] for the year 2015, assuming a conservative 50% growth in the flows between 2015 and 2020. Notes: Countries on X-Axis are ordered as countryicountryj, thus linked to ωij. The opposite weight ωji=1−ωij refers to the weight of countryjcountryi

Table 4 MAE comparison of the different models with and without day/corridor dummies
Table 5 Averaged feature ranking across models and specifications without directional priors
Table 6 Feature ranking including lagged epidemiological conditions

A.1 Representativeness of facebook data

A critical question relates to the representativeness of FB data. On average, comparison between the number of FB users and population size by country reveals that the number of FB users with location services enabled represents 5.6% of the population in Europe, on average. However, this fraction is rather uniform across countries, which suggests that FB data might be representative of the population as a whole. This does not mean, however, that data on FB users are representative of the population of cross-border movers, as using internet connection or activating location services abroad can be costly [33]. To address this point, we produce estimates of the number of daily cross-border inflows in the pre-Covid period, and compare them with FB data. We collect data on commuting flows from the Labor Force Survey for 2019,Footnote 10data on air passengers as of March 2019,Footnote 11 and estimates of immigration flows from [34].Footnote 12

Overall, Table A.1 shows that migration flows account for a tiny (say negligible) proportion of daily movements, implying that FB data mostly capture commuting flows and business travels, as well as holiday-maker’s moves in vacation periods. This implies that the largest flows are observed between contiguous countries and business partners. Combining corridors available in our sample and corridors involving non-EU countries, Figure A.1 shows that FB data as of the first week of March 2020 are highly correlated with estimates of actual flows in the year 2019, which suggests that FB data might also be representative of cross-border movements.

A.2 Parameters of the mL methods

We describe in this section the parameters of the different ML models used in this paper. In order to tune the models, a boosting sample corresponding of 70% of the whole dataset is used (from the 29th February 2020 until the 16th November 2020). We use the following parameters:

  • K Nearest Neighbors. The models uses 5 neighbors and it has the distance metric as weight function.

  • Gradient Boosting. The Gradient Boosting model has least square regression as loss function, a learning rate of 0.1 and a number of estimators of 100 (which is equivalent to the number of boosting stages to perform).

  • Multi-Layer Perceptron. The MLP Regressor has 2 hidden layers of 100 neurons each. The activation function for the layer is the rectified linear unit function. It uses an adam solver with minibatches of size 50.

A.3 Additional figures and tables

Some additional figures and tables cited in the core of the text are found below.

Availability of data and materials

The datasets supporting the conclusions of this article are available in the Zenodo repository Cross-border-Mobility-Responses-to-Covid-19-in-Europe, https://doi.org/10.5281/zenodo.4719559

Notes

  1. We understand mobility as “the movement of human beings in space and time” [1]. This broad definition is in line with the mobility measure that we use in our quantitative analysis, as that includes all types of movement.

  2. By big data, we mean extremely large data sets that may be analyzed computationally to reveal patterns, trends and associations, especially those related to human behaviors and interactions, according to the Oxford Languages (https://languages.oup.com/google-dictionary-en/).

  3. The recent OECD migration outlook reveals that issuance of new visas and permits in OECD countries plummeted by 46% in the first half of 2020 (by 72% in the second quarter), as compared with the same period in 2019. However, international migrants and refugees account for a tiny proportion, not a say a negligible proportion, of daily cross-border movements of people between European countries. Daily flows predominantly consist of commuting workers and business travels.

  4. A number of recent studies such as [17, 18] use data extracted from FB advertising API to identify migrants and model migration patterns. Note that our data source is different as we only have access to aggregate bilateral counts, without further details on movers’ characteristics. The advantage of our data source is, however, that it is available on a daily basis.

  5. We also find a strong association between the number of FB users and population size at the regional level. Results are available upon request.

  6. Note that we only observe the total number of movers within a day, and neither the duration of their stay nor the direction of the movement. The vast majority of movers cross the border twice as they return to their origin country within a day or so. This is why we model cross-border traffic and use a trick below to account for the presumed direction of the flows.

  7. Epidemiological conditions are likely to subject to measurement errors. For example, testing and tracing practices played an important role in determining the number of detected cases.

  8. Gravity models aim to predict bilateral flows between two entities based on their size, on their economic characteristics, and on the distance between them [22].

  9. It could be further supplemented with a leave-one-out sensitivity analysis on countries [32].

  10. Labor Force Surveys are provided by Eurostat. We use data on commuting by NUTS2 region, and allocate commuters aged 15 to 64 by destination on the basis of contiguity links with other countries. When a NUTS region shares several contiguity links, we share the number of commuters proportionately to the level of income per capita in the destination country. To identify daily movements, we assume that each commuter spend 4 days per week in the place of work, and multiply LFS number by 4/7.

  11. We use monthly statistics on Air passengers from Eurostat by airport of departure, and allocate them by destination using FB data. To identify daily movements, we divide monthly flows by 30.

  12. The latter are available until the year 2015; we multiply them by 1.5 to proxy the level in 2019. To identify daily movements, we divide annual inflows by 360.

Abbreviations

EM:

Econometric Modelling

EU:

European Union

FB:

Facebook Inc.

G-Boost:

Gradient Boosting

KNN:

K-nearest neighbors method

MAE:

Mean Absolute Error

ML:

Machine Learning

MLP:

Multi-layer Perceptron

NPIs:

non-pharmaceutical interventions

OxCGRT:

Oxford Covid-19 Response Tracker

PCA:

Principal Component Analysis

RSME:

Root Mean Squared Error

References

  1. Zhang A, Wang A. Human Mobility and Smart City. Sch Community Encycl. 2021. https://encyclopedia.pub/entry/6764. Accessed 12 04 2022.

  2. Caselli F, Grigoli F, Sandri D, Spilimbergo A. Mobility under the COVID-19 pandemic: Asymmetric effects across gender and age. IMF Econ Rev. 2020; 70(1):105–38.

    Article  Google Scholar 

  3. Chetty R, Friedman J, Hendren N, Stepner M, Team, The Opportunity Insights. The Economic Impacts of COVID-19: Evidence from a New Public Database Built Using Private Sector Data. Cambridge: National Bureau of Economic Research; 2020, p. w27431. https://doi.org/10.3386/w27431.

    Book  Google Scholar 

  4. Goolsbee A, Syverson C. Fear, lockdown, and diversion: Comparing drivers of pandemic economic decline 2020. J Public Econ. 2021; 193:104311. https://doi.org/10.1016/j.jpubeco.2020.104311.

    Article  PubMed  Google Scholar 

  5. Maloney WF, Taskin T. Determinants of Social Distancing and Economic Activity During COVID-19: A Global View. Rochester: Social Science Research Network; 2020. https://papers.ssrn.com/abstract=3599572.

    Book  Google Scholar 

  6. Glaeser E, Gorback C, Redding S. How Much does COVID-19 Increase with Mobility? Evidence from New York and Four Other U.S. Cities. Cambridge: National Bureau of Economic Research; 2020, p. w27519. https://doi.org/10.3386/w27519.

    Book  Google Scholar 

  7. Basellini U, Alburez-Gutierrez D, Del Fava E, Perrotta D, Bonetti M, Camarda CG, Zagheni E. Linking excess mortality to mobility data during the first wave of covid-19 in england and wales. SSM-Popul Health. 2021; 14:100799.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Chen S, Yang J, Yang W, Wang C, Bärnighausen T. Covid-19 control in china during mass population movements at new year. Lancet. 2020; 395:764–66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Ruktanonchai NW, Floyd JR, Lai S, Ruktanonchai CW, Sadilek A, Rente-Lourenco P, Ben X, Carioli A, Gwinn J, Steele JE, Prosper O, Schneider A, Oplinger A, Eastham P, Tatem AJ. Assessing the impact of coordinated covid-19 exit strategies across europe. Science. 2020; 369(6510):1465–70. https://doi.org/10.1126/science.abc5096. http://arxiv.org/abs/https://science.sciencemag.org/content/369/6510/1465.full.pdf.

    Article  CAS  PubMed  Google Scholar 

  10. Lee K, Worsnop CZ, Grépin KA, Kamradt-Scott A. A global coordination on cross-border travel and trade measures crucial to covid-19 response. Lancet. 2020; 395:1593–95.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Errett NA, Sauer LM, Rutkow L. An integrative review of the limited evidence of international travel bans as an emerging infectious disease disaster control measure. J Emerg Manag. 2020; 18:7–14.

    Article  PubMed  Google Scholar 

  12. Tian H, Liu Y, Li Y, Wu2 C-H, Chen B, Kraemer MUG, Li B, Cai J, Xu B, Yang Q, Wang B, Yang P, Cui Y, Song Y, Zheng P, Wang Q, Bjornstad ON, Yang R, Grenfell BT, Pybus OG, Dye C. An investigation of transmission control measures during the first 50 days of the covid-19 epidemic in china. Science. 2020; 368:638–42.

  13. Kraemer MUG, Yang C-H, Gutierrez B, Wu C-H, Klein B, Pigott DM, Group OC-DW, du Plessis L, Faria NR, Li R, Hanage WP, Brownstein JS, Layan M, Vespignani A, Tian H, Dye C, Pybus OG, Scarpino SV. The effect of human mobility and control measures on the covid-19 epidemic in china. Science. 2020; 497:493–97.

    Article  Google Scholar 

  14. Chinazzi M, Davis JT, Ajelli M, Gioannini C, Litvinova M, Merler S, Piontti APY, Mu K, Ross L, Sun K, Viboud C, Xiong X, Yu H, Halloran ME, Longini Jr IM, Vespignani A. The effect of travel restrictions on the spread of the 2019 novel coronavirus (covid-19) outbreak. Science. 2020; 368:395–400.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Askitas N, Tatsiramos K, Verheyden B. Estimating worldwide effects of non-pharmaceutical interventions on covid-19 incidence and population mobility patterns using a multiple-event study. Sci Rep. 2021; 11(1972):41598–02181442.

    Google Scholar 

  16. Inc. F, Facebook Data for Good. https://dataforgood.fb.com/. Accessed 19 Feb 2021.

  17. Rampazzo F, Bijak J, Vitali A, Weber I, Zagheni E. A Framework for Estimating Migrant Stocks Using Digital Traces and Survey Data: An Application in the United Kingdom. Demography. 2021; 58(6):2193–218. https://doi.org/10.1215/00703370-9578562.

    Article  PubMed  Google Scholar 

  18. Alexander M, Zagheni E, Polimis K. The impact of hurricane maria on out-migration from puerto rico: Evidence from facebook data. Popul Dev Rev. 2019; 45(3):617–30.

    Article  Google Scholar 

  19. Ribeiro FN, Benevenuto F, Zagheni E. How Biased is the Population of Facebook Users? Comparing the Demographics of Facebook Users with Census Data to Generate Correction Factors. In: 12th ACM Conference on Web Science. WebSci ’20. New York: Association for Computing Machinery: 2020. p. 325–34. https://doi.org/10.1145/3394231.3397923. Accessed 01 Dec 2021.

    Google Scholar 

  20. Dong E, Du H, Gardner L. An interactive web-based dashboard to track COVID-19 in real time. Lancet Infect Dis. 2020; 20(5):533–34. https://doi.org/10.1016/S1473-3099(20)30120-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Hale T, Webster S, Petherick A, Phillips T, Kira B. Oxford COVID-19 Government Response Tracker, Blavatnik School of Government. Data retrieved from Oxford Oxford Covid-19 Government Response Tracker (OxCGRT) Github page. 2020. https://github.com/OxCGRT/covid-policy-tracker. Accessed 25 Apr 2021.

  22. Beine M, Bertoli S, Fernández-Huertas Moraga J. A Practitioners’ Guide to Gravity Models of International Migration. World Econ. 2016; 39(4):496–512. https://doi.org/10.1111/twec.12265. _eprint: https://onlinelibrary.wiley.com/doi/pdf/10.1111/twec.12265. Accessed 16 Feb 2022.

  23. Letouz E, Purser M, Rodriguez F, Cummins M. Revisiting the Migration-Development Nexus: A Gravity Model Approach. Hum Dev Rep Off (HDRO), U N Dev Program (UNDP), Hum Dev Res Pap (2009 to present). 2009. https://mpra.ub.uni-muenchen.de/19227/.

  24. Peterson LE. K-nearest neighbor. Scholarpedia. 2009; 4(2):1883. https://doi.org/10.4249/scholarpedia.1883. Accessed 24 Mar 2021.

    Article  Google Scholar 

  25. Friedman JH. Greedy Function Approximation: A Gradient Boosting Machine. Ann Stat. 2001; 29(5):1189–232. Publisher: Institute of Mathematical Statistics. Accessed 24 Mar 2021.

    Article  Google Scholar 

  26. Hinton GE. Connectionist learning procedures In: Kodratoff Y, Michalski RS, editors. Machine Learning. San Francisco: Morgan Kaufmann: 1990. p. 555–610. https://doi.org/10.1016/B978-0-08-051055-2.50029-8.

    Google Scholar 

  27. Glorot X, Bengio Y. Understanding the difficulty of training deep feedforward neural networks. In: Proceedings of the Thirteenth International Conference on Artificial Intelligence And Statistics. Sardinia: JMLR Workshop and Conference Proceedings: 2010. p. 249–56. ISSN: 1938-7228.

    Google Scholar 

  28. Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, Blondel M, Prettenhofer P, Weiss R, Dubourg V, Vanderplas J, Passos A, Cournapeau D, Brucher M, Perrot M, Duchesnay E. Scikit-learn: Machine learning in Python. J Mach Learn Res. 2011; 12:2825–30.

    Google Scholar 

  29. Altmann A, Toloşi L, Sander O, Lengauer T. Permutation importance: a corrected feature importance measure. Bioinformatics. 2010; 26(10):1340–7. https://doi.org/10.1093/bioinformatics/btq134.

    Article  CAS  PubMed  Google Scholar 

  30. Molnar C. Interprtable machine learning: Leanpub; 2018. https://leanpub.com/interpretable-machine-learning.

  31. Vannoni M, McKee M, Semenza JC, Bonell C, Stuckler D. Using volunteered geographic information to assess mobility in the early phases of the covid-19 pandemic: a cross-city time series analysis of 41 cities in 22 countries from march 2nd to 26th 2020. Glob Health. 2020; 16(1):1–9.

    Article  Google Scholar 

  32. Stoto MA, Woolverton A, Kraemer J, Barlow P, Clarke M. COVID-19 data are messy: analytic methods for rigorous impact analyses with imperfect data. Glob Health. 2022; 18(1):2. https://doi.org/10.1186/s12992-021-00795-0. Accessed 09 Feb 2022.

    Article  Google Scholar 

  33. Spyratos S, Vespe M, Natale F, Weber I, Zagheni E, Rango M. Quantifying international human mobility patterns using facebook network data. PloS ONE. 2019; 14(10):0224134.

    Article  Google Scholar 

  34. Abel GJ, Cohen JE. Bilateral international migration flow estimates for 200 countries. Sci Data. 2019; 6(1):1–13.

    Article  Google Scholar 

Download references

Acknowledgments

The second, third and fourth authors acknowledge financial support from the ARC convention on “New approaches to understanding and modeling global migration trends” (convention 18/23-091).

Funding

Felix Stips gratefully acknowledges financial support from Luxembourg National Research Fund (FNR) under the ACROSS Doctoral Training Unit (grant n PRIDE19/14302992).

Author information

Authors and Affiliations

Authors

Contributions

F.D. and P.S. conceived the presented idea. F.D., N.G., P.S. and F.S. contributed to the design and implementation of the research, to the analysis of the results and to the writing of the manuscript. F.D., F.S. and P.S. collected the data. F.S. cleaned the data and produced descriptive graphs. N.G. and F.S. implemented the EM and ML modelling. All authors discussed the results, drafted and critically revised the article. The authors read and approved the final manuscript.

Corresponding author

Correspondence to Pierre Schaus.

Ethics declarations

Ethics approval and consent to participate

No ethics approval or consent required.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Docquier, F., Golenvaux, N., Nijssen, S. et al. Cross-border mobility responses to COVID-19 in Europe: new evidence from facebook data. Global Health 18, 41 (2022). https://doi.org/10.1186/s12992-022-00832-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12992-022-00832-6

Keywords