 Research
 Open access
 Published:
The impact of national and international travel on spatiotemporal transmission of SARSCoV2 in Belgium in 2021
BMC Infectious Diseases volume 23, Article number: 428 (2023)
Abstract
Background
The Severe Acute Respiratory Syndrome Coronavirus 2 (SARSCoV2) has rapidly spread over the world and caused tremendous impacts on global health. Understanding the mechanism responsible for the spread of this pathogen and the impact of specific factors, such as human mobility, will help authorities to tailor interventions for future SARSCoV2 waves or newly emerging airborne infections. In this study, we aim to analyze the spatiotemporal transmission of SARSCoV2 in Belgium at municipality level between January and December 2021 and explore the effect of different levels of human travel on disease incidence through the use of counterfactual scenarios.
Methods
We applied the endemicepidemic modelling framework, in which the disease incidence decomposes into endemic, autoregressive and neighbourhood components. The spatial dependencies among areas are adjusted based on actual connectivity through mobile network data. We also took into account other important factors such as international mobility, vaccination coverage, population size and the stringency of restriction measures.
Results
The results demonstrate the aggravating effect of international travel on the incidence, and simulated counterfactual scenarios further stress the alleviating impact of a reduction in national and international travel on epidemic growth. It is also clear that local transmission contributed the most during 2021, and municipalities with a larger population tended to attract a higher number of cases from neighboring areas.
Conclusions
Although transmission between municipalities was observed, local transmission was dominant. We highlight the positive association between the mobility data and the infection spread over time. Our study provides insight to assist health authorities in decisionmaking, particularly when the disease is airborne and therefore likely influenced by human movement.
Background
The emergence and rapid worldwide spread of the Severe Acute Respiratory Syndrome Coronavirus 2 (SARSCoV2) and its implied Coronavirus Disease 2019 (COVID19) has tremendously impacted global health. Several studies suggested that the transmission dynamics of SARSCoV2 are explained by several spatiotemporal factors, including demographic factors [1], contact networks of individuals [2], containment measures [3, 4], and human mobility [5, 6]. In Belgium, approximately 650,000 COVID19positive cases were confirmed in 2020, and this number increased to 1.5 million confirmed cases over the course of 2021, accounting for 12.9% of the Belgian population. Cases were heterogeneously distributed across municipalities (see, e.g., Fig. 1). Few studies have examined the effects of control strategies and mobility on the spatial spread of COVID19 in the country, and these studies have focused on waves in 2020 [2, 7, 8]. A study incorporating spatiotemporal heterogeneous factors (e.g., demographic factors, containment measures, human mobility) is therefore advisable to understand the transmission dynamics of SARSCoV2 in Belgium, especially in 2021.
Human mobility has been identified earlier as a quintessential factor in the spread of airborne infections [9,10,11,12]. Hence, to mitigate SARSCoV2 spread in Belgium, one of the early responses was the installment of severe restrictions on travel, both at the national and international level. Specifically, international travel was banned during the first wave, MarchMay 2020, and nonessential foreign travel to and from Belgium was prohibited between January 27 and April 19, 2021. However, being a small, centrally located European country, the relatively large influx of travelers  a fraction of which hosting the infection  especially in periods with little or no travel restrictions (e.g., from September 2021), potentially has an important impact on viral spread. The human mobility data in Belgium is available at several spacetime resolutions. In our study, two main sources are used: summary data on national mobility from the mobile network and the number of incoming travelers to the country. As these data reflect changes in mobility patterns at different time periods, they can be used to describe disease spread. More specifically, the key question arises whether we use the available travel data to quantify the impact of travel on the spread of COVID19, in order to use them as an important factor to simulate and predict the disease incidence in different time period.
In response to the pandemic, the nonpharmaceutical interventions (NPIs) issued by the Belgian government have changed with time, leading to variations in the stringency of restriction measures. In 2021, relatively stringent measures were still in place, but these no longer included lockdown and full school closure, although there were still longer periods with mask mandates (e.g., in public transport) and closure of culture and event sectors. A precise estimation of their effects on the evolution of infections is challenging in the face of sparse data availability. Fortunately, the Oxford Coronavirus Government Tracker (OxCGRT) provides a proxy metric, the socalled stringency index, to assess the level of strictness of social policies imposed at a national scale [13]. The stringency index has shown to be useful in recent spatiotemporal studies [14, 15]. Furthermore, vaccination campaigns have largely contributed to an important reduction of the risk of infection and of more severe outcomes such as hospitalization and death [16,17,18,19]. During 2021, mass vaccination campaigns were widely implemented in the Belgian population. However, vaccination coverages were heterogeneous across Belgian geographical units (see, e.g., Fig. 1D, and Faes et al. (2022) [20]). For pathogens such as measles, earlier findings indicated that the heterogeneity in vaccination coverage was associated with disease incidence [21, 22]. Notwithstanding, the impact of vaccination and the relation between vaccine uptake in the population and disease incidence may be different between different vaccinepreventable diseases. Besides, given the importance of demographic characteristics, which are known to vary by location and usually available in spatial resolution, it could not be neglected in studying the transmission process of COVID19. Therefore, we expect that the combination of spacetime available data will help to better unravel the complex dynamics of SARSCoV2, in this paper with a focus on Belgium.
Spatiotemporal modelling of infectious disease spread based on routine incidence data has received great interest in recent years. In the context of COVID19, various modeling frameworks have been employed. For example, compartmental models (e.g., SEIR, metapopulation models) study the interaction between individuals in communities at different spatial scales [2, 23]. Spatial and temporal point processes have also been used to study the disease’s dynamics. Additional information like demographics, human mobility, and control measures have been considered as possible factors influencing transmission [24,25,26]. However, these models can be complex, require highresolution data that may not always be available (e.g., susceptible proportion in communities), and can be timeconsuming to compute. An adequate method is offered by a multivariate time series model in a socalled endemicepidemic (EE) statistical modelling framework introduced by Held et al. (2005) [27]. An important aspect is that it employs a common regression framework, and the model’s parameters can be conveniently estimated using standard optimization techniques, such as maximum likelihood. In essence, the timespace dependence of infections is determined by both background risk independent of the epidemic and autoregression on past counts of within and betweengeographical units. The flexibility of the EE framework enables us to include potential variables and model their effects on disease dynamics, in both time and space. Properly tailored models have been used to investigate and predict the evolution of the COVID19 pandemic [3, 14, 15, 28].
In this paper, we examine both temporal and spatial patterns in disease spread and the way their dynamics are influenced by human mobility, the stringency of NPIs, and the effectiveness thereof. We apply the EE modelling framework with the aim to (i) analyze the spatiotemporal transmission of SARSCoV2 in Belgium between January and December 2021 and (ii) explore the effect of different levels of human travel on disease incidence through the use of counterfactual scenarios. We leverage the available data on the daily time series for case incidence, national mobility between Belgian municipalities, the number of incoming international travelers to the country, the stringency of containment policies, vaccination coverage, and demographic variables. This study aims to provide public health authorities with essential information to make informed policy decision for mitigating the impact of COVID19 as well as that of future pandemics, especially in relation to travel restrictions.
Methods
In this section, we first describe the data sources used in our analyses. A summary of the data used can be found in Table 1. Then, we introduce a spatial dynamic model to assess the impact of travel rate on COVID19 spread. This impact is further illustrated through counterfactual scenarios with different travel levels. Lastly, sensitivity analyses to assess the stability of model estimates are described.
Data
Case data
The daily number of confirmed COVID19 cases in each of the 581 municipalities in Belgium was retrieved for the period between January 4, 2021 and December 5, 2021 from Sciensano, the Scientific Institute for Public Health in Belgium. Figure 1A shows the evolution of the Belgian COVID19 epidemic during 2021. More specifically, the total number of reported cases over time shows two waves: (1) a wave between January and June, 2021, representing the third COVID19 wave in Belgium, and (2) a wave starting from September 2021 onwards. The daily number of cases started to increase again from July onwards and fluctuated around 2,500 confirmed cases before dramatically soaring to approximately 25,000 cases per day at the beginning of December. Furthermore, the lower left panel (Fig. 1B) shows the distribution of confirmed cases by Belgian municipality using heat colors. For areas that are darker, the number of confirmed cases is higher.
National and international mobility data
To assess the amount of travel between municipalities in Belgium, summary information from mobile network data is used, available between September 1, 2020 and June 6, 2021 [29]. These data provide information about the average proportion of time that customers spent in their residing municipality as well as in other municipalities (if any). Based on this, we can present the connectivity among municipalities through a mobility matrix (Fig. 2), calculated as the logtransformed mean percentage time spent in a municipality between September 2020 and June 2021, and averaged this temporal information to obtain a static network. The monthly mobility matrices over this time period are graphically depicted in Fig. A6 in the Appendix.
The international travel rate is available at the provincial level as the number of incoming travelers per 100 inhabitants. This information is based on the COVID19 weekly reports by Sciensano. Figure 1C shows the time trend over the study period. The rate for all provinces was at the minimum during February—April 2021, when a travel ban on nonessential travel was in place (from January 27 until April 19, 2021). After this period, travel rates started to increase, to reach a peak around August and November 2021. Brussels remained the region with the highest travel rate during 2021.
Pharmaceutical and nonpharmaceutical interventions
The number of individuals that received at least one dose of a COVID19 vaccine in 2021 was obtained from Sciensano [30]). Note that these data are available at municipality level and in weekly resolution. The proportion of vaccinated individuals in each Belgian province is depicted in Fig. 1D. Differences in vaccine coverage among provinces are clearly visible. WestFlanders had the highest proportion of individuals receiving at least one vaccine dose, while the lowest level was observed in Brussels.
The Stringency Index ranges from 0 to 100 and reflects the level of severity of nonpharmaceutical mitigation measures imposed at the national level (Oxford Coronavirus Government Tracker 2022) [13]. It is a composite measure of nine response metrics, including school and workplace closures, cancellation of public events, limitations on public gatherings, closure of or restrictions on public transport, stayathome mandates, public awareness campaigns, restrictions on internal mobility, and international movement control. Higher values of the Stringency Index imply that stricter measures were in place. Figure 1E presents the evolution of this index for Belgium throughout 2021. The highest values were observed during February–April 2021, with a peak in April 2021, at the time of the third wave.
Population data
The population data of each Belgian municipality in 2020 was retrieved from the Belgian statistical office [31]. The total population size was approximately 11.5 million individuals and was assumed to be constant during 2021.
Spatial dynamic model
Let \(Y_{it}\) denote the number of confirmed COVID19 cases in municipality i (\(i = 1, \ldots , 581\)) on day t (\(t = 1, \ldots , T\)), with T the length of the time window considered in this analysis. Covariate information is represented by \(\varvec{x}_{it}\). Assuming that the number of new confirmed cases \(Y_{it}\) at day t depends on the (series of) past observations \(Y_{i(td)}\), \(d = 1, \ldots , D\), with up to \(D = 7\) days being of importance to determine \(Y_{it}\), the spatiotemporal model is formulated as [27, 32, 33]:
with mean incidence \(\mu _{it}\) and variance \(\mu _{it}(1+\psi \mu _{it})\). Here \(\psi \ge 0\) is the dispersion parameter to be estimated. Note that if \(\psi = 0\), mean and variance are equal and \(Y_{it}\) follows a Poisson distribution. The conditional mean \(\mu _{it}\) is modelled as:
The natural logarithm of the first term, referred to as the endemic component \(\epsilon _{it}\), and interpreted as the background risk, equals
with \(N_i^* = {N_i}/{100,000}\) defining an offset. The second term is the local epidemic or autoregressive component, describing the local epidemic spread in time. The autoregressive parameter \(\lambda _{it}\), is modeled as
Finally, the third term of mean \(\mu _{it}\) represents the global epidemic or neighbourhood component. Its corresponding parameter \(\phi _{it}\) is modeled as
Two types of weights are included in the model. First, the normalized Poisson weights \(u_d\) in (1) represent the probability for a serial interval (i.e., the average time (in days) between symptom onset in an infectious individual and symptoms emerging in a newly infected individual when both were in close contact) of length d days with d taking values up to \(D = 7\) days [33]. Second, the weights \(w_{ji}^*\) in the neighbourhood component express the strength of connectivity between pairs of municipalities. The definition of the mobility weights is based on the estimated connectivity over time inferred from mobility data, and a comparison is made from defining the connectivity merely on the distance between two municipalities.
In the model, we include componentspecific intercepts (\(\alpha ^{(\epsilon )}\), \(\alpha ^{(\lambda )}\), \(\alpha ^{(\phi )}\)). We account for various covariates, including a weekend effect (WE), population size (N) (in the local and global epidemic components), proportion of nonvaccinated individuals as a proxy for susceptibility (sus), travel rate (travel), stringency index (SI), and calendar time (period). Specifically, two time periods, i.e., from June to September 2021 and from October to December 2021, are considered in the analysis in order to distinguish between different phases in the epidemic. The effect of the proportion of susceptible individuals (sus) is accommodated in the model by considering the natural logarithm of the proportion of residents that are nonvaccinated (i.e., the fraction of individuals without any of the COVID19 vaccine doses) [21, 22]. We add a quadratic term of the travel rate which was centered prior to inclusion in the model.
Inference and model comparison
Model comparison is done based on the Akaike Information Criteria (AIC) [34]. The model with the lowest AIC is selected. Parameter inference is based on maximum likelihood estimation in R package surveillance [35] version 1.20.1 and hhh4addon version 0.0.0.0.9014 [33] under R version 4.0.5 [36].
Counterfactual simulation of case counts
To investigate the impact of travel, we consider three scenarios with different levels of travel and simulate the expected number of cases. In the first scenario, the national mobility and international travel rate is the same as observed. In the second one, the international travel rate is set at the lowest level during the study period, while national mobility is kept the same as observed. In the third case, we show what the situation would be if there would be no mobility within the country.
Sensitivity analyses
In order to check the stability of our model estimates, we performed four sensitivity analyses. First, the travel rate is estimated based on the number of incoming travellers from the Passenger Locator Form (PLF) data (Appendix Fig. A1). This form was required to be filled in by all passengers entering Belgium during the study period. Second, we use serial intervals based on a literature review. This includes a gamma distribution with shape and rate parameters equal to 2.29 and 0.36, respectively [37], and a lognormal distribution with parameters 1.38 and 0.563 [38]. Third, the Stringency Index was adjusted to exclude travel restrictions (i.e., restrictions on internal movement and international travel controls) (Appendix Fig. A2), in order to reduce dependencies between the (original) Stringency Index versus the national mobility and the international travel rate in our model. Finally, we represent municipality connectivity by a power law matrix instead of mobility matrices [39].
Results
Figure 3 depicts the model fit (panel (a)) and the contribution of each component (panel (b)). It can be seen that our model provides a good overall fit to the observed epidemic curve. Nevertheless, the estimated case numbers from around November onward are not as close to the observed data points as in the previous period. This may result as a result of the emergence of the Omicron variant, which started circulating in Belgium from December 2021 onward. The relative contribution of each component to the estimated daily number of cases is illustrated in the right panel of the same figure. The local epidemic spread component usually accounts for the largest part of the circulation, throughout the study period. A sudden increase in the share of the global spread was observed in June  September 2021. Overall, the endemic, local epidemic and global epidemic components accounted for 15.8%, 73.9%, and 10.3% of the estimated cumulative number of cases, respectively.
The exponentiated model parameters, representing the multiplicative effect at the original scale, with their 95% confidence intervals, are presented in Table 2. As expected, the number of confirmed cases in weekends is lower as compared to weekdays due to reporting delays, underreporting, or delays in healthseeking behaviour. Municipalities with a larger population size have a stronger spatial epidemic spread while its effect on withinmunicipality epidemic spread is found to be insignificant. An increase in the proportion of susceptible individuals (as approximated by the fraction of nonvaccinated individuals) tends to lead to an increase in the number of confirmed cases, albeit that the effect is also insignificant. Interestingly, implementation of stricter measures is associated with a higher disease spread even though the effect is minimal. The two periods (i.e., June to September, 2021 and October to December, 2021) show different effects (compared to the reference period January to May, 2021) in each of the three model components. Both periods have a lower mean baseline incidence as compared to the first period (January to May, 2021). In the period of June to September, the speed of local spread is reduced as compared to in the reference period, while in the period October to December it increases again to a level which is higher than January to May, 2021. At the same time, the speed of disease spread between municipalities is elevated in the period of June to September as compared to the period of January to May, and not significantly different between the period of October to December and the period of January to May. Similar results are observed in the sensitivity analyses.
The impact of travel rate on the mean number of cases is visualized in Fig. 4. Note that since the travel rate is centered, a value of zero on the xaxis corresponds to the mean travel rate on the original scale. From Fig. 4, it can be seen that the baseline circulation of COVID19 gets higher if there is more international travel, and in addition that the spatial spread is higher when international travel increases. Indeed, in the endemic and neighbourhood components the baseline mean incidence (assuming a mean centered travel rate of zero) increases multiplicatively with a factor evolving in a nonlinear way in relation to the travel rate, while the opposite (although albeit small) effect holds for the autoregressive component (red line). Furthermore, for large travel rates the increase in multiplicative effect with increasing travel rates stabilizes in the endemic case (green line), thereby deviating from an exponential increase in case of a linear travel rate term at the linear predictor scale.
The comparison between the estimated lag distribution from our model and from literature is illustrated in Fig. 5. This lag distribution reflects the distribution of the serial interval (i.e., the time between symptom onset days of two paired cases). Parameter estimates from using an alternative definition of the mobility matrix, the adjusted Stringency Index, and from fixing the serial intervals are given in Tables A1 and A2 in the Appendix, respectively. Similar conclusions are found in these analyses.
Figure 6 illustrates the impact of travel by means of three simulated scenarios for three periods: 18 January to 25 April (period 1), 26 April to 29 August (period 2) and 30 August to 5 December (period 3). In the first scenario (travel rate as observed), the simulation was able to capture the observed epidemic curve despite a slight underestimation in period 1 and overestimation in period 2 as compared to the observed cases (93.9% and 114%, respectively). In the second scenario (reduced international travel) and third scenario (reduced national travel), the simulation clearly shows that restricting international travel and local mobility would greatly reduce the epidemic growth. No connectivity among municipalities would imply the largest effect on the spread in period 1, whereas in period 2 and 3, minimizing incoming international travelers would be the most effective. Specifically, a travel rate fixed to the minimum level that was observed would have prevented 11.4% of cases in period 1 up to 68.3% in period 3. When there is no connectivity among municipalities, we also observe a decrease in cases with the largest reduction being 43.8% of cases in period 3 as compared to the first scenario (see Table A3).
Discussion
Using the endemicepidemic modelling framework, we performed a detailed analysis concerning the spatiotemporal characteristics of the COVID19 pandemic at municipality level in Belgium. We successfully investigated the impact of different levels of human mobility, both between municipalities as well as by incoming travellers, on the daily incidence during 2021, while considering important covariates, such as vaccination on the one hand, and the Stringency Index as a proxy for the effect of NPIs that are in place.
From the model results, we conclude that a relatively high fraction of cases originated from local transmission within the municipality itself, suggesting that endogenous transmission contributed the most to the COVID19 epidemic, at least during 2021. The spatiotemporal contribution was relatively low in the earlier part of the year. However, it increased during summer, the time when travel related policies became less stringent, suggesting that a larger fraction of cases during this period were imported from other areas. In certain municipalities, the progression of the disease remained local, though an elevated contribution of the neighbourhood component was visible in a part of the Brussels Capital Region, especially in the second half of the year (Appendix Fig. A3). This is plausible since intensive connectivity between municipalities on the one hand and a higher fraction of incoming travelers per capita were observed in the Brussels Capital Region.
It is suggested that differences in levels of travel played a role on daily confirmed cases. This hypothesis was confirmed by constructing three counterfactual scenarios, each with a different level of mobility. Reducing the mobility of individuals in the second and third counterfactual scenario corresponded with an attenuated COVID19 case growth. Previous studies have investigated the relationship between disease data and human movement, but these studies were more focused on the areal mobility matrices, which illustrate the origin–destination connectivity between geographical units [1, 12, 40,41,42]. Our work is one of a few studies that account for mobility patterns outside of the study settings, as a potential source of infection [43, 44]. Nevertheless, regardless of the transmission model applied, the reduction on travel is effective in reducing the spread and size of epidemics. While using aggregated mobility flows can help estimating reductions in incidence at the population level, it would be intriguing to examine the patterns of spread influenced by changes in individuallevel travel behaviors. While not without challenges in data acquisition, one might take into account the finer resolution of mobility data that involves information on individual contacts or clustering patterns. Adopting these mobility and travel behaviors may benefit to identify disease hotspots, adapt appropriate locallevel regulations and control strategies.
From the estimated model coefficients, we found strong evidence of commuterdriven spread when we accounted for the unitspecific population in the neighbourhood component, indicating that regions with a larger population are expected to attract a bigger number of cases from their neighbors. The existence of such an agglomeration effect was also confirmed in previous studies [14, 15, 22, 45]. We identify that vaccination coverage of the first dose in Belgium did not have a significant association with fewer case counts. In endemicepidemic models of areal count time series (socalled hhh4 models within the surveillance package in R), the rate of disease diffusion is allowed to be proportional to the number of susceptible people, which can be inferred from the number of nonvaccinated individuals [21]. In our study, only the local coverage of the first COVID19 vaccine dose was used in the model. The uncertainty about the true immunity levels in the population inferred from onedose coverage should be noted. By including the proportion of unvaccinated individuals in the model, we assumed that the susceptibility level in the population over time is proportionate to the timevarying proportion of unvaccinated individuals. More specifically, immunity in unvaccinated individuals as a result of previous infection with SARSCoV2 would imply an overestimation of the timevarying susceptibility at the population level when solely informed by vaccination coverage. However, when assuming that a constant fraction of unvaccinated individuals has immunity within the study period (and if waning of (vaccineinduced) immunity is limited), the fraction of unvaccinated individuals is proportional to population susceptibility. We acknowledge that this is a strong assumption and an investigation of the sensitivity of our results with regard to this assumption is considered a topic for future research. Besides, for a disease such as measles, the assumption that depletion of susceptibles due to vaccination decreases the transmission risk is a plausible one, but it arguably does not hold for COVID19 given waning of humoral immunity with time since vaccination, the need for a second dose to maintain high protection against severe disease and the fact that vaccination does not completely prevent infection. Given agespecific differences in the fraction of the population that is protected against SARSCoV2 infection due to prior exposure or vaccination [46, 47], an agestratified model would be a better approach to evaluate the impact of vaccination on the disease dynamics. To our surprise, we found that the positive association between stringency of NPIs and the new cases was not significant. Although recent studies indicated that this factor significantly contributed to the lower transmission of cases [14, 15], such a beneficial effect was not observed in this context. One possible reason is that the Stringency Index provides a surrogate picture of control measures by assigning one estimated value for the entire Belgian population. A more granular perspective would provide additional information about the role and impact of NPIs in controlling the epidemic. For example, it would be interesting to replace the Stringency Index with locallevel indices to investigate the impact of NPIs. Furthermore, compliance to intervention measures and its variation across Belgium is not captured by this measure. Stringency should ideally be combined with information on adherence to NPIs and with social contact data. In this way, studying the impact of NPIs on spacetime spread of the virus is a relevant topic of further research, conditional on the availability of sufficient data. Another aspect that may contribute to this controversial result is the fact that we used a lag of 7 days in the model. It might be unlikely that the impact of control measures in reducing transmission will be noticeable within a week, as it may take 912 days or longer to manifest in the trajectory of the epidemics [40]. Additional research could investigate the correlation between the delayed effects of NPIs and changes in disease incidence, for example, by examining the rate of change in reproductive number [48]. Furthermore, we discovered that by fitting a model with the Stringency Index, excluding travel restrictions, we found that dependencies between the Stringency Index versus the national mobility and the international travel rate in our model had minimal influence on model estimates. The results of this sensitivity analysis are shown in Table A2 in the Appendix.
During our study period, the Alpha variant was most prevalent over the largest part of the first half of 2021, while the Delta variant was dominant in the second half of the year. Although the Beta and Gamma variants were concomitantly identified over the first half of 2021, they did not play an important role in subsequent waves of SARSCoV2 variants observed in Belgium [30, 49]. Furthermore, previous studies have shown that both the Alpha and Delta variants have 50%70% higher transmissibility compared to the preexisting variants, and they impacted individuals with similar demographic and comorbidity characteristics [49, 50]. The emergence of variants of concern (VOCs) such as Alpha and Delta variants, which are capable of immune evasion in different molecular mechanisms, may suggest different transmission risks. Adeyinka et al., 2022 suggested that the evaluation of the spread of VOCs should account for other factors such as vaccine uptake, the strength of control measures that are in place, and mobility patterns [51]. However, due to the small size of Belgium as a country, significant geographical variation in the prevalence of VOCs across different provinces is unlikely, as least during our study period  even though there may be some regional differences. Moreover, our model considered different time periods that did not precisely align with the periods when certain VOCs became more dominant. Although it was feasible to adjust the time periods to match, this did not lead to improvement in the model, but rather provoked convergency issue. Therefore, the findings of this study should be interpreted while bearing in mind these VOCs.
A notable strength of our modelling strategy is that we were able to mine the mobile network activity to establish the spatial dependence of COVID19 transmission. Mobility data has gradually shown its potential and importance in (spatial) analyses of infectious diseases [1, 24, 40, 43]. Our study has underscored its usefulness in modelling the connectivity between municipalities. In other words, they allow the modeller to construct betweenmunicipality weights. Although a power law, which describes the shorttime human movement behavior when there is missing information about travel, improved model fit in terms of AIC (Appendix Table A2), it makes more sense to apply the mobile network data to the model as it reflects the “true” activity logs of individuals geographically. In addition, we used the daily case counts and used observations to 7 days back in time. This allowed us to avoid biases in serial interval assumptions when using the aggregated data (in weekly data, for example) [33]. The estimated lag distribution in our study was in line with findings in the studies of Grimée et al. (2022) and Ssentongo et al. (2021) [3, 14]. The sensitivity analyses with different literaturebased serial intervals also showed moderate stability of our results. Furthermore, the covariates used in the model (e.g., travel rate, vaccine coverage, stringency index) were timevarying. This approach is particularly important in the context of COVID19, where policies and interventions tend to be highly variable, potentially over short periods of time.
Several limitations need to be mentioned. First, our mobile network data was not collected during the whole study period. This may lessen the representative quality for the human activity logs in the time period under investigation. Besides, we used a static spatiallydependent network to capture the interaction between localities rather than the timedependent connectivity matrix, which may better describe the commuting behavior of individuals. Second, we assumed a homogeneous mixing scenario, which prevents us from examining the interaction between individuals in disease spread. Endeavors have been made to clarify the transmission dynamics in subgroups (age, gender) using social contact data [52,53,54]. Nowadays, such diseaserelated data and data on human behavior are collected and stored more systematically with high resolution and subject to specific circumstances. Future extensions of the EE framework would ideally address some or all of these limitations.
Conclusions
In conclusion, our research focus centered on capturing the COVID19 diffusion at municipality level in Belgium through a spacetime dynamic model and identifying covariates that serve as contributing factors to the virus spread. Even though transmission between municipalities was observed during the study period, local transmission was dominant. We highlighted the positive association between the mobility data and the infection spread over time. Insight from our study can assist in decisions of public health policies, particularly when the disease is likely influenced by human movement.
Availability of data and materials
The datasets generated and/or analysed during the current study are available from https://github.com/OxCGRT/covidpolicytracker/tree/master/data (stringency index), https://epistat.sciensano.be/covid/ (vaccine coverage) and https://statbel.fgov.be/en (population).
The datasets for case counts, national and international mobility are available from the authors upon reasonable request (contact Minh Hanh Nguyen via minhhanh.nguyen@uhasselt.be or Thi Huyen Trang Nguyen via thihuyentrang.nguyen@uhasselt.be) and with permission of Sciensano.
Abbreviations
 SARSCoV2:

Severe Acute Respiratory Syndrome Coronavirus 2
 COVID19:

Coronavirus Disease 2019
 NPIs:

Nonpharmaceutical interventions
 OxCGRT:

Oxford Coronavirus Government Tracker
 EE:

Endemicepidemic
 AIC:

Akaike Information Criteria
 PLF:

Passenger Locator Form
References
Pullano G, Valdano E, Scarpa N, Rubrichi S, Colizza V. Evaluating the effect of demographic factors, socioe conomic factors, and risk aversion on mobility during the COVID19 epidemic in France under lockdown: A populationbased study. Lancet Digit Health. 2020;2(12):e63849.
Coletti P, Libin P, Petrof O, Willem L, Abrams S, Herzog SA, et al. A datadriven metapopulation model for the Belgian COVID19 epidemic: Assessing the impact of lockdown and exit strategies. BMC Infect Dis. 2021;21(1):1–12.
Grimée M, DunbarMBN, Hofmann F, Held L, et al. Modelling the effect of a border closure between Switzerland and Italy on the spatiotemporal spread of COVID19 in Switzerland. Spat Stat. 2022;49:100552.
Vinceti M, Filippini T, Rothman KJ, Ferrari F, Goffi A, Maffeis G, et al. Lockdown timing and efficacy in controlling COVID19 using mobile phone tracking. EClinicalMedicine. 2020;25:100457.
Cot C, Cacciapaglia G, Sannino F. Mining Google and Apple mobility data: Temporal anatomy for COVID19 social distancing. Sci Rep. 2021;11(1):1–8.
Gatalo O, Tseng K, Hamilton A, Lin G, Klein E. Associations between phone mobility data and COVID19 cases. Lancet Infect Dis. 2021;21(5):e111.
Rollier M, Miranda G, Vergeynst J, Meys J, Alleman T, Van De Vijver E, et al. Mobility and the spatially heterogeneous spread of SARSCoV2 in Belgium. 2022. arXiv preprint arXiv:2202.11528.
Dellicour S, Linard C, Van Goethem N, Da Re D, Artois J, Bihin J, et al. Investigating the drivers of the spatiotemporal heterogeneity in COVID19 hospital incidence—Belgium as a study case. Int J Health Geogr. 2021;20(1):1–11.
Geilhufe M, Held L, Skrovseth SO, Simonsen GS, Godtliebsen F. Power law approximations of movement network data for modeling infectious disease spread. Biom J. 2014;56(3):363–82.
Changruenngam S, Bicout DJ, Modchang C. How the individual human mobility spatiotemporally shapes the disease transmission dynamics. Sci Rep. 2020;10(1):1–13.
Balcan D, Colizza V, Gonçalves B, Hu H, Ramasco JJ, Vespignani A. Multiscale mobility networks and the spatial spreading of infectious diseases. Proc Natl Acad Sci. 2009;106(51):21484–9.
Charu V, Zeger S, Gog J, Bjornstad ON, Kissler S, Simonsen L, et al. Human mobility and the spatial transmission of influenza in the United States. PLoS Comput Biol. 2017;13(2):e1005382.
Hale T, Angrist N, Goldszmidt R, Kira B, Petherick A, Phillips T, et al. A global panel database of pandemic policies (Oxford COVID19 Government Response Tracker). Nat Hum Behav. 2021;5(4):529–38.
Ssentongo P, Fronterre C, Geronimo A, Greybush SJ, Mbabazi PK, Muvawala J, et al. PanAfrican evolution of withinand betweencountry COVID19 dynamics. Proc Natl Acad Sci. 2021;118(28):e2026664118.
Celani A, Giudici P. Endemic–epidemic models to understand COVID19 spatiotemporal evolution. Spat Stat. 2022;49:100528.
Pawlowski C, Lenehan P, Puranik A, Agarwal V, Venkatakrishnan A, Niesen MJ, et al. FDAauthorized mRNA COVID19 vaccines are effective per realworld evidence synthesized across a multistate health system. Med. 2021;2(8):979–92.
Dagan N, Barda N, Kepten E, Miron O, Perchik S, Katz MA, et al. BNT162b2 mRNA Covid19 vaccine in a nationwide mass vaccination setting. N Engl J Med. 2021;384:1412–23.
Braeye T, Cornelissen L, Catteau L, Haarhuis F, Proesmans K, De Ridder K, et al. Vaccine effectiveness against infection and onwards transmission of COVID19: Analysis of Belgian contact tracing data, JanuaryJune 2021. Vaccine. 2021;39(39):5456–60.
Braeye T, Catteau L, Brondeel R, van Loenhout JA, Proesmans K, Cornelissen L, et al. Vaccine effectiveness against onward transmission of SARSCoV2infection by variant of concern and time since vaccination, Belgian contact tracing, 2021. Vaccine. 2022;40(22):3027–37.
Faes C, Molenberghs G, Hens N, Van Bortel L, Vandeboel N, Pellens K, et al. Geographical variation of COVID19 vaccination coverage, ethnic diversity and population composition in Flanders. Vaccine: X. 2022;11:100194.
Herzog S, Paul M, Held L. Heterogeneity in vaccination coverage explains the size and occurrence of measles epidemics in German surveillance data. Epidemiol Infect. 2011;139(4):505–15.
Robert A, Kucharski AJ, Funk S. The impact of local vaccine coverage and recent incidence on measles transmission in France between 2009 and 2018. BMC Med. 2022;20(1):1–20.
Chang S, Pierson E, Koh PW, Gerardin J, Redbird B, Grusky D, et al. Mobility network models of COVID19 explain inequities and inform reopening. Nature. 2021;589(7840):82–7.
Slater JJ, Brown PE, Rosenthal JS, Mateu J. Capturing spatial dependence of COVID19 case counts with cellphone mobility data. Spat Stat. 2022;49:100540. https://www.sciencedirect.com/science/article/pii/S2211675321000506.
Giudici P, Pagnottoni P, Spelta A. Network selfexciting point processes to measure health impacts of COVID19. 2021. SSRN 3892998.
Chiang WH, Liu X, Mohler G. Hawkes process modeling of COVID19 with mobility leading indicators and spatial covariates. Int J Forecast. 2022;38(2):505–20.
Held L, Höhle M, Hofmann M. A statistical framework for the analysis of multivariate infectious disease surveillance counts. Stat Model. 2005;5(3):187–99.
Giuliani D, Dickson MM, Espa G, Santi F. Modelling and predicting the spatiotemporal spread of COVID19 in Italy. BMC Infect Dis. 2020;20(1):1–10.
EnsoyMusoro C, Nguyen MH, Hens N, Molenberghs G, Faes C. Spatiotemporal model to investigate COVID19 spread accounting for the mobility amongst municipalities. Spat Spatiotemporal Epidemiol. 2023;45:100568.
COVID19 Data. 2022. https://epistat.wivisp.be/covid/. Accessed 05 Apr 2022.
The Belgian statistical office. https://statbel.fgov.be/en. Accessed 05 Apr 2022.
Paul M, Held L, Toschke AM. Multivariate modelling of infectious disease surveillance data. Stat Med. 2008;27(29):6250–67.
Bracher J, Held L. Endemicepidemic models with discretetime serial interval distributions for infectious disease prediction. Int J Forecast. 2022;38(3):1221–33.
Akaike H. Information theory and an extention of the maximum likelihood principle. In: 2nd International Symposium on Information Theory, 1973. Akademiai Kiado; 1973. p. 26781.
Meyer S, Held L, Höhle M. Spatiotemporal analysis of epidemic phenomena using the R package surveillance. 2014. arXiv preprint arXiv:1411.0416.
R Core Team. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2022. https://www.Rproject.org/.
Bi Q, Wu Y, Mei S, Ye C, Zou X, Zhang Z, et al. Epidemiology and transmission of COVID19 in 391 cases and 1286 of their close contacts in Shenzhen, China: A retrospective cohort study. Lancet Infect Dis. 2020;20(8):911–9.
Ito K, Piantham C, Nishiura H. Predicted dominance of variant Delta of SARSCoV2 before Tokyo olympic games, Japan, July 2021. Eurosurveillance. 2021;26(27):2100570.
Meyer S, Held L. Powerlaw models for infectious disease spread. Ann Appl Stat. 2014;8(3):1612–39.
Badr HS, Du H, Marshall M, Dong E, Squire MM, Gardner LM. Association between mobility patterns and COVID19 transmission in the USA: A mathematical modelling study. Lancet Infect Dis. 2020;20(11):1247–54.
Arándiga F, Baeza A, CorderoCarrión I, Donat R, Martí MC, Mulet P, et al. A spatialtemporal model for the evolution of the COVID19 pandemic in Spain including mobility. Mathematics. 2020;8(10):1677.
Iacus SM, Santamaria C, Sermi F, Spyratos S, Tarchi D, Vespe M. Human mobility and COVID19 initial dynamics. Nonlinear Dyn. 2020;101(3):1901–19.
Persson J, Parie JF, Feuerriegel S. Monitoring the COVID19 epidemic with nationwide telecommunication data. Proc Natl Acad Sci. 2021;118(26):e2100664118.
Han X, Xu Y, Fan L, Huang Y, Xu M, Gao S. Quantifying COVID19 importation risk in a dynamic network of domestic cities and international countries. Proc Natl Acad Sci. 2021;118(31):e2100201118.
Xia Y, Bjornstad ON, Grenfell BT. Measles metapopulation dynamics: A gravity model for epidemiological coupling and dynamics. Am Nat. 2004;164(2):267–81.
Moghadas SM, Sah P, Shoukat A, Meyers LA, Galvani AP. Population immunity against COVID19 in the United States. Ann Intern Med. 2021;174(11):1586–91.
Herzog SA, De Bie J, Abrams S, Wouters I, Ekinci E, Patteet L, et al. Seroprevalence of IgG antibodies against SARSCoV2a serial prospective crosssectional nationwide study of residual samples, Belgium, March to October 2020. Eurosurveillance. 2022;27(9):2100419.
Rypdal K. The tipping effect of delayed interventions on the evolution of COVID19 incidence. Int J Environ Res Public Health. 2021;18(9):4484.
Seyler L, Van Nedervelde E, De Cock D, Mann C, Pien K, Allard SD, et al. Surfing theWaves: Differences in Hospitalised COVID19 Patients across 4 Variant Waves in a Belgian University Hospital. Viruses. 2023;15(3):618.
MendiolaPastrana IR, LópezOrtiz E, Río de la LozaZamora JG, González J, GómezGarcía A, LópezOrtiz G. SARSCoV2 variants and clinical outcomes: a systematic review. Life. 2022;12(2):170.
Adeyinka DA, Neudorf C, Camillo CA, Marks WN, Muhajarine N. COVID19 Vaccination and Public Health Countermeasures on Variants of Concern in Canada: Evidence From a Spatial Hierarchical Cluster Analysis. JMIR Public Health Surveill. 2022;8(5):e31968.
Meyer S, Held L. Incorporating social contact data in spatiotemporal models for infectious disease spread. Biostatistics. 2017;18(2):338–51.
Ondrikova N, Clough H, Cunliffe N, IturrizaGomaraM, Vivancos R, Harris J. Understanding norovirus reporting patterns in England: A mixed model approach. BMC Public Health. 2021;21(1):19.
Nguyen THT, Faes C, Hens N. Measles epidemic in Southern Vietnam: An agestratified spatiotemporal model for infectious disease counts. Epidemiol Infect. 2022;150:e169.
Acknowledgements
Not applicable.
Funding
NH and CF acknowledge support from the European Union’s Horizon 2020 research and innovation programme  project EpiPose (Grant agreement number 101003688).
Author information
Authors and Affiliations
Contributions
MHN and THTN participated in data aggregation, statistical analysis, interpretation of the study results, and writing the initial draft of the manuscript. GM, SA, NH, and CF participated in study conception, study design and analysis, interpretation of the study results, supervision and critical revision of the manuscript. All authors approved the final version of the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
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.
About this article
Cite this article
Nguyen, M.H., Nguyen, T.H.T., Molenberghs, G. et al. The impact of national and international travel on spatiotemporal transmission of SARSCoV2 in Belgium in 2021. BMC Infect Dis 23, 428 (2023). https://doi.org/10.1186/s12879023083689
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s12879023083689