Infodemics: A new challenge for public health.

The COVID-19 information epidemic, or "infodemic," demonstrates how unlimited access to information may confuse and influence behaviors during a health emergency. However, the study of infodemics is relatively new, and little is known about their relationship with epidemics management. Here, we discuss unresolved issues and propose research directions to enhance preparedness for future health crises.

Evaluating the Impact of Meningococcal Vaccines With Synthetic Controls.

Invasive meningococcal disease (IMD) has a low and unpredictable incidence, presenting challenges for real-world evaluations of meningococcal vaccines. Traditionally, meningococcal vaccine impact is evaluated by predicting counterfactuals from pre-immunization IMD incidences, possibly controlling for IMD in unvaccinated age groups, but the selection of controls can influence results. We retrospectively applied a synthetic control (SC) method, previously used for pneumococcal disease, to data from 2 programs for immunization of infants against serogroups B and C IMD in England and Brazil. Time series of infectious/noninfectious diseases in infants and IMD cases in older unvaccinated age groups were used as candidate controls, automatically combined in a SC through Bayesian variable selection. SC closely predicted IMD in absence of vaccination, adjusting for nontrivial changes in IMD incidence. Vaccine impact estimates were in line with previous assessments. IMD cases in unvaccinated age groups were the most frequent SC-selected controls. Similar results were obtained when excluding IMD from control sets and using other diseases only, particularly respiratory diseases and measles. Using non-IMD controls may be important where there are herd immunity effects. SC is a robust and flexible method that addresses uncertainty introduced when equally plausible controls exhibit different post-immunization behaviors, allowing objective comparisons of IMD programs between countries.

Interplay between mobility, multi-seeding and lockdowns shapes COVID-19 local impact.

Assessing the impact of mobility on epidemic spreading is of crucial importance for understanding the effect of policies like mass quarantines and selective re-openings. While many factors affect disease incidence at a local level, making it more or less homogeneous with respect to other areas, the importance of multi-seeding has often been overlooked. Multi-seeding occurs when several independent (non-clustered) infected individuals arrive at a susceptible population. This can lead to independent outbreaks that spark from distinct areas of the local contact (social) network. Such mechanism has the potential to boost incidence, making control efforts and contact tracing less effective. Here, through a modeling approach we show that the effect produced by the number of initial infections is non-linear on the incidence peak and peak time. When case importations are carried by mobility from an already infected area, this effect is further enhanced by the local demography and underlying mixing patterns: the impact of every seed is larger in smaller populations. Finally, both in the model simulations and the analysis, we show that a multi-seeding effect combined with mobility restrictions can explain the observed spatial heterogeneities in the first wave of COVID-19 incidence and mortality in five European countries. Our results allow us for identifying what we have called epidemic epicenter: an area that shapes incidence and mortality peaks in the entire country. The present work further clarifies the nonlinear effects that mobility can have on the evolution of an epidemic and highlight their relevance for epidemic control.

The impact of news exposure on collective attention in the United States during the 2016 Zika epidemic.

In recent years, many studies have drawn attention to the important role of collective awareness and human behaviour during epidemic outbreaks. A number of modelling efforts have investigated the interaction between the disease transmission dynamics and human behaviour change mediated by news coverage and by information spreading in the population. Yet, given the scarcity of data on public awareness during an epidemic, few studies have relied on empirical data. Here, we use fine-grained, geo-referenced data from three online sources-Wikipedia, the GDELT Project and the Internet Archive-to quantify population-scale information seeking about the 2016 Zika virus epidemic in the U.S., explicitly linking such behavioural signal to epidemiological data. Geo-localized Wikipedia pageview data reveal that visiting patterns of Zika-related pages in Wikipedia were highly synchronized across the United States and largely explained by exposure to national television broadcast. Contrary to the assumption of some theoretical epidemic models, news volume and Wikipedia visiting patterns were not significantly correlated with the magnitude or the extent of the epidemic. Attention to Zika, in terms of Zika-related Wikipedia pageviews, was high at the beginning of the outbreak, when public health agencies raised an international alert and triggered media coverage, but subsequently exhibited an activity profile that suggests nonlinear dependencies and memory effects in the relation between information seeking, media pressure, and disease dynamics. This calls for a new and more general modelling framework to describe the interaction between media exposure, public awareness and disease dynamics during epidemic outbreaks.

The Framing of COVID-19 in Italian Media and Its Relationship with Community Mobility: A Mixed-Method Approach.

Media framing of epidemics was found to influence public perceptions and behaviors in experiments, yet no research has been conducted on real-world behaviors during public health crises. We examined the relationship between Italian news media coverage of COVID-19 and compliance with stay-at-home orders, which could impact the spread of epidemics. We used a computational method for framing analysis (ANTMN) and combined it with Google's Community Mobility data. A time-series analysis using vector autoregressive models showed that the Italian media used media frames that were largely congruent with ones used by journalists in other countries: A scientific frame focusing on symptoms and health effects, a containment frame focusing on attempts to ameliorate risks, and a social frame, focusing on political and social impact. The prominence of different media frames over time was associated with changes in Italians' mobility patterns. Specifically, we found that the social frame was associated with increased mobility, whereas the containment frame was associated with decreased mobility. The results demonstrate that the ways the news media discuss epidemics can influence changes in community mobility, above and beyond the effect of the number of deaths per day.

How Search Engine Data Enhance the Understanding of Determinants of Suicide in India and Inform Prevention: Observational Study.

BACKGROUND: India is home to 20% of the world's suicide deaths. Although statistics regarding suicide in India are distressingly high, data and cultural issues likely contribute to a widespread underreporting of the problem. Social stigma and only recent decriminalization of suicide are among the factors hampering official agencies' collection and reporting of suicide rates. OBJECTIVE: As the product of a data collaborative, this paper leverages private-sector search engine data toward gaining a fuller, more accurate picture of the suicide issue among young people in India. By combining official statistics on suicide with data generated through search queries, this paper seeks to: add an additional layer of information to more accurately represent the magnitude of the problem, determine whether search query data can serve as an effective proxy for factors contributing to suicide that are not represented in traditional datasets, and consider how data collaboratives built on search query data could inform future suicide prevention efforts in India and beyond. METHODS: We combined official statistics on demographic information with data generated through search queries from Bing to gain insight into suicide rates per state in India as reported by the National Crimes Record Bureau of India. We extracted English language queries on "suicide," "depression," "hanging," "pesticide," and "poison". We also collected data on demographic information at the state level in India, including urbanization, growth rate, sex ratio, internet penetration, and population. We modeled the suicide rate per state as a function of the queries on each of the 5 topics considered as linear independent variables. A second model was built by integrating the demographic information as additional linear independent variables. RESULTS: Results of the first model fit (R2) when modeling the suicide rates from the fraction of queries in each of the 5 topics, as well as the fraction of all suicide methods, show a correlation of about 0.5. This increases significantly with the removal of 3 outliers and improves slightly when 5 outliers are removed. Results for the second model fit using both query and demographic data show that for all categories, if no outliers are removed, demographic data can model suicide rates better than query data. However, when 3 outliers are removed, query data about pesticides or poisons improves the model over using demographic data. CONCLUSIONS: In this work, we used search data and demographics to model suicide rates. In this way, search data serve as a proxy for unmeasured (hidden) factors corresponding to suicide rates. Moreover, our procedure for outlier rejection serves to single out states where the suicide rates have substantially different correlations with demographic factors and query rates.

Wearable Proximity Sensors for Monitoring a Mass Casualty Incident Exercise: Feasibility Study.

BACKGROUND: Over the past several decades, naturally occurring and man-made mass casualty incidents (MCIs) have increased in frequency and number worldwide. To test the impact of such events on medical resources, simulations can provide a safe, controlled setting while replicating the chaotic environment typical of an actual disaster. A standardized method to collect and analyze data from mass casualty exercises is needed to assess preparedness and performance of the health care staff involved. OBJECTIVE: In this study, we aimed to assess the feasibility of using wearable proximity sensors to measure proximity events during an MCI simulation. In the first instance, our objective was to demonstrate how proximity sensors can collect spatial and temporal information about the interactions between medical staff and patients during an MCI exercise in a quasi-autonomous way. In addition, we assessed how the deployment of this technology could help improve future simulations by analyzing the flow of patients in the hospital. METHODS: Data were obtained and collected through the deployment of wearable proximity sensors during an MCI functional exercise. The scenario included 2 areas: the accident site and the Advanced Medical Post, and the exercise lasted 3 hours. A total of 238 participants were involved in the exercise and classified in categories according to their role: 14 medical doctors, 16 nurses, 134 victims, 47 Emergency Medical Services staff members, and 27 health care assistants and other hospital support staff. Each victim was assigned a score related to the severity of his/her injury. Each participant wore a proximity sensor, and in addition, 30 fixed devices were placed in the field hospital. RESULTS: The contact networks show a heterogeneous distribution of the cumulative time spent in proximity by the participants. We obtained contact matrices based on the cumulative time spent in proximity between the victims and rescuers. Our results showed that the time spent in proximity by the health care teams with the victims is related to the severity of the patient's injury. The analysis of patients' flow showed that the presence of patients in the rooms of the hospital is consistent with the triage code and diagnosis, and no obvious bottlenecks were found. CONCLUSIONS: Our study shows the feasibility of the use of wearable sensors for tracking close contacts among individuals during an MCI simulation. It represents, to our knowledge, the first example of unsupervised data collection-ie, without the need for the involvement of observers, which could compromise the realism of the exercise-of face-to-face contacts during an MCI exercise. Moreover, by permitting detailed data collection about the simulation, such as data related to the flow of patients in the hospital, such deployment provides highly relevant input for the improvement of MCI resource allocation and management.

Fast and accurate dynamic estimation of field effectiveness of meningococcal vaccines.

BACKGROUND: Estimating the effectiveness of meningococcal vaccines with high accuracy and precision can be challenging due to the low incidence of the invasive disease, which ranges between 0.5 and 1 cases per 100,000 in Europe and North America. Vaccine effectiveness (VE) is usually estimated with a screening method that combines in one formula the proportion of meningococcal disease cases that have been vaccinated and the proportion of vaccinated in the overall population. Due to the small number of cases, initial point estimates are affected by large uncertainties and several years may be required to estimate VE with a small confidence interval. METHODS: We used a Monte Carlo maximum likelihood (MCML) approach to estimate the effectiveness of meningococcal vaccines, based on stochastic simulations of a dynamic model for meningococcal transmission and vaccination. We calibrated the model to describe two immunization campaigns: the campaign against MenC in England and the Bexsero campaign that started in the UK in September 2015. First, the MCML method provided estimates for both the direct and indirect effects of the MenC vaccine that were validated against results published in the literature. Then, we assessed the performance of the MCML method in terms of time gain with respect to the screening method under different assumptions of VE for Bexsero. RESULTS: MCML estimates of VE for the MenC immunization campaign are in good agreement with results based on the screening method and carriage studies, yet characterized by smaller confidence intervals and obtained using only incidence data collected within 2 years of scheduled vaccination. Also, we show that the MCML method could provide a fast and accurate estimate of the effectiveness of Bexsero, with a time gain, with respect to the screening method, that could range from 2 to 15 years, depending on the value of VE measured from field data. CONCLUSIONS: Results indicate that inference methods based on dynamic computational models can be successfully used to quantify in near real time the effectiveness of immunization campaigns against Neisseria meningitidis. Such an approach could represent an important tool to complement and support traditional observational studies, in the initial phase of a campaign.

On the use of human mobility proxies for modeling epidemics.

Human mobility is a key component of large-scale spatial-transmission models of infectious diseases. Correctly modeling and quantifying human mobility is critical for improving epidemic control, but may be hindered by data incompleteness or unavailability. Here we explore the opportunity of using proxies for individual mobility to describe commuting flows and predict the diffusion of an influenza-like-illness epidemic. We consider three European countries and the corresponding commuting networks at different resolution scales, obtained from (i) official census surveys, (ii) proxy mobility data extracted from mobile phone call records, and (iii) the radiation model calibrated with census data. Metapopulation models defined on these countries and integrating the different mobility layers are compared in terms of epidemic observables. We show that commuting networks from mobile phone data capture the empirical commuting patterns well, accounting for more than 87% of the total fluxes. The distributions of commuting fluxes per link from mobile phones and census sources are similar and highly correlated, however a systematic overestimation of commuting traffic in the mobile phone data is observed. This leads to epidemics that spread faster than on census commuting networks, once the mobile phone commuting network is considered in the epidemic model, however preserving to a high degree the order of infection of newly affected locations. Proxies' calibration affects the arrival times' agreement across different models, and the observed topological and traffic discrepancies among mobility sources alter the resulting epidemic invasion patterns. Results also suggest that proxies perform differently in approximating commuting patterns for disease spread at different resolution scales, with the radiation model showing higher accuracy than mobile phone data when the seed is central in the network, the opposite being observed for peripheral locations. Proxies should therefore be chosen in light of the desired accuracy for the epidemic situation under study.

Generalized contact matrices allow integrating socioeconomic variables into epidemic models.

Variables related to socioeconomic status (SES), including income, ethnicity, and education, shape contact structures and affect the spread of infectious diseases. However, these factors are often overlooked in epidemic models, which typically stratify social contacts by age and interaction contexts. Here, we introduce and study generalized contact matrices that stratify contacts across multiple dimensions. We demonstrate a lower-bound theorem proving that disregarding additional dimensions, besides age and context, might lead to an underestimation of the basic reproductive number. By using SES variables in both synthetic and empirical data, we illustrate how generalized contact matrices enhance epidemic models, capturing variations in behaviors such as heterogeneous levels of adherence to nonpharmaceutical interventions among demographic groups. Moreover, we highlight the importance of integrating SES traits into epidemic models, as neglecting them might lead to substantial misrepresentation of epidemic outcomes and dynamics. Our research contributes to the efforts aiming at incorporating socioeconomic and other dimensions into epidemic modeling.

Estimating household contact matrices structure from easily collectable metadata.

Contact matrices are a commonly adopted data representation, used to develop compartmental models for epidemic spreading, accounting for the contact heterogeneities across age groups. Their estimation, however, is generally time and effort consuming and model-driven strategies to quantify the contacts are often needed. In this article we focus on household contact matrices, describing the contacts among the members of a family and develop a parametric model to describe them. This model combines demographic and easily quantifiable survey-based data and is tested on high resolution proximity data collected in two sites in South Africa. Given its simplicity and interpretability, we expect our method to be easily applied to other contexts as well and we identify relevant questions that need to be addressed during the data collection procedure.

Human mobility and time spent at destination: impact on spatial epidemic spreading.

Host mobility plays a fundamental role in the spatial spread of infectious diseases. Previous theoretical works based on the integration of network theory into the metapopulation framework have shown that the heterogeneities that characterize real mobility networks favor the propagation of epidemics. Nevertheless, the studies conducted so far assumed the mobility process to be either Markovian (in which the memory of the origin of each traveler is lost) or non-Markovian with a fixed traveling time scale (in which individuals travel to a destination and come back at a constant rate). Available statistics however show that the time spent by travelers at destination is characterized by wide fluctuations, ranging from a single day up to several months. Such varying length of stay crucially affects the chance and duration of mixing events among hosts and may therefore have a strong impact on the spread of an emerging disease. Here, we present an analytical and a computational study of epidemic processes on a complex subpopulation network where travelers have memory of their origin and spend a heterogeneously distributed time interval at their destination. Through analytical calculations and numerical simulations we show that the heterogeneity of the length of stay alters the expression of the threshold between local outbreak and global invasion, and, moreover, it changes the epidemic behavior of the system in case of a global outbreak. Additionally, our theoretical framework allows us to study the effect of changes in the traveling behavior in response to the infection, by considering a scenario in which sick individuals do not leave their home location. Finally, we compare the results of our non-Markovian framework with those obtained with a classic Markovian approach and find relevant differences between the two, in the estimate of the epidemic invasion potential, as well as of the timing and the pattern of its spatial spread. These results highlight the importance of properly accounting for host trip duration in epidemic models and open the path to the inclusion of such an additional layer of complexity to the existing modeling approaches.

The limits of human mobility traces to predict the spread of COVID-19: A transfer entropy approach.

Mobile phone data have been widely used to model the spread of COVID-19; however, quantifying and comparing their predictive value across different settings is challenging. Their quality is affected by various factors and their relationship with epidemiological indicators varies over time. Here, we adopt a model-free approach based on transfer entropy to quantify the relationship between mobile phone-derived mobility metrics and COVID-19 cases and deaths in more than 200 European subnational regions. Using multiple data sources over a one-year period, we found that past knowledge of mobility does not systematically provide statistically significant information on COVID-19 spread. Our approach allows us to determine the best metric for predicting disease incidence in a particular location, at different spatial scales. Additionally, we identify geographic and demographic factors, such as users' coverage and commuting patterns, that explain the (non)observed relationship between mobility and epidemic patterns. Our work provides epidemiologists and public health officials with a general-not limited to COVID-19-framework to evaluate the usefulness of human mobility data in responding to epidemics.

On the forecastability of food insecurity.

Food insecurity, defined as the lack of physical or economic access to safe, nutritious and sufficient food, remains one of the main challenges included in the 2030 Agenda for Sustainable Development. Near real-time data on the food insecurity situation collected by international organizations such as the World Food Programme can be crucial to monitor and forecast time trends of insufficient food consumption levels in countries at risk. Here, using food consumption observations in combination with secondary data on conflict, extreme weather events and economic shocks, we build a forecasting model based on gradient boosted regression trees to create predictions on the evolution of insufficient food consumption trends up to 30 days in to the future in 6 countries (Burkina Faso, Cameroon, Mali, Nigeria, Syria and Yemen). Results show that the number of available historical observations is a key element for the forecasting model performance. Among the 6 countries studied in this work, for those with the longest food insecurity time series, that is Syria and Yemen, the proposed forecasting model allows to forecast the prevalence of people with insufficient food consumption up to 30 days into the future with higher accuracy than a naive approach based on the last measured prevalence only. The framework developed in this work could provide decision makers with a tool to assess how the food insecurity situation will evolve in the near future in countries at risk. Results clearly point to the added value of continuous near real-time data collection at sub-national level.

Crime, inequality and public health: a survey of emerging trends in urban data science.

Urban agglomerations are constantly and rapidly evolving ecosystems, with globalization and increasing urbanization posing new challenges in sustainable urban development well summarized in the United Nations' Sustainable Development Goals (SDGs). The advent of the digital age generated by modern alternative data sources provides new tools to tackle these challenges with spatio-temporal scales that were previously unavailable with census statistics. In this review, we present how new digital data sources are employed to provide data-driven insights to study and track (i) urban crime and public safety; (ii) socioeconomic inequalities and segregation; and (iii) public health, with a particular focus on the city scale.

Association of close-range contact patterns with SARS-CoV-2: a household transmission study.

Background: Households are an important location for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission, especially during periods when travel and work was restricted to essential services. We aimed to assess the association of close-range contact patterns with SARS-CoV-2 transmission. Methods: We deployed proximity sensors for two weeks to measure face-to-face interactions between household members after SARS-CoV-2 was identified in the household, in South Africa, 2020-2021. We calculated the duration, frequency, and average duration of close-range proximity events with SARS-CoV-2 index cases. We assessed the association of contact parameters with SARS-CoV-2 transmission using mixed effects logistic regression accounting for index and household member characteristics. Results: We included 340 individuals (88 SARS-CoV-2 index cases and 252 household members). On multivariable analysis, factors associated with SARS-CoV-2 acquisition were index cases with minimum Ct value <30 (aOR 16.8 95% CI 3.1-93.1) vs >35, and female contacts (aOR 2.5 95% CI 1.3-5.0). No contact parameters were associated with acquisition (aOR 1.0-1.1) for any of the duration, frequency, cumulative time in contact, or average duration parameters. Conclusions: We did not find an association between close-range proximity events and SARS-CoV-2 household transmission. Our findings may be due to study limitations, that droplet-mediated transmission during close-proximity contacts plays a smaller role than airborne transmission of SARS-CoV-2 in the household, or due to high contact rates in households. Funding: Wellcome Trust (Grant number 221003/Z/20/Z) in collaboration with the Foreign, Commonwealth, and Development Office, United Kingdom.

Real-time numerical forecast of global epidemic spreading: case study of 2009 A/H1N1pdm.

BACKGROUND: Mathematical and computational models for infectious diseases are increasingly used to support public-health decisions; however, their reliability is currently under debate. Real-time forecasts of epidemic spread using data-driven models have been hindered by the technical challenges posed by parameter estimation and validation. Data gathered for the 2009 H1N1 influenza crisis represent an unprecedented opportunity to validate real-time model predictions and define the main success criteria for different approaches. METHODS: We used the Global Epidemic and Mobility Model to generate stochastic simulations of epidemic spread worldwide, yielding (among other measures) the incidence and seeding events at a daily resolution for 3,362 subpopulations in 220 countries. Using a Monte Carlo Maximum Likelihood analysis, the model provided an estimate of the seasonal transmission potential during the early phase of the H1N1 pandemic and generated ensemble forecasts for the activity peaks in the northern hemisphere in the fall/winter wave. These results were validated against the real-life surveillance data collected in 48 countries, and their robustness assessed by focusing on 1) the peak timing of the pandemic; 2) the level of spatial resolution allowed by the model; and 3) the clinical attack rate and the effectiveness of the vaccine. In addition, we studied the effect of data incompleteness on the prediction reliability. RESULTS: Real-time predictions of the peak timing are found to be in good agreement with the empirical data, showing strong robustness to data that may not be accessible in real time (such as pre-exposure immunity and adherence to vaccination campaigns), but that affect the predictions for the attack rates. The timing and spatial unfolding of the pandemic are critically sensitive to the level of mobility data integrated into the model. CONCLUSIONS: Our results show that large-scale models can be used to provide valuable real-time forecasts of influenza spreading, but they require high-performance computing. The quality of the forecast depends on the level of data integration, thus stressing the need for high-quality data in population-based models, and of progressive updates of validated available empirical knowledge to inform these models.

Evidence of pandemic fatigue associated with stricter tiered COVID-19 restrictions.

Despite the availability of effective vaccines against SARS-CoV-2, non-pharmaceutical interventions remain an important part of the effort to reduce viral circulation caused by emerging variants with the capability of evading vaccine-induced immunity. With the aim of striking a balance between effective mitigation and long-term sustainability, several governments worldwide have adopted systems of tiered interventions, of increasing stringency, that are calibrated according to periodic risk assessments. A key challenge remains in quantifying temporal changes in adherence to interventions, which can decrease over time due to pandemic fatigue, under such kind of multilevel strategies. Here, we examine whether there was a reduction in adherence to tiered restrictions that were imposed in Italy from November 2020 through May 2021, and in particular we assess whether temporal trends in adherence depended on the intensity of the restrictions adopted. We analyzed daily changes in movements and in residential time, combining mobility data with the restriction tier enforced in the Italian regions. Through mixed-effects regression models, we identified a general trend of reduction in adherence and an additional effect of faster waning associated with the most stringent tier. We estimated both effects being of the same order of magnitude, suggesting that adherence decreased twice as fast during the strictest tier as in the least stringent one. Our results provide a quantitative measure of behavioral responses to tiered interventions-a metric of pandemic fatigue-that can be integrated into mathematical models to evaluate future epidemic scenarios.

Mapping urban socioeconomic inequalities in developing countries through Facebook advertising data.

Ending poverty in all its forms everywhere is the number one Sustainable Development Goal of the UN 2030 Agenda. To monitor the progress toward such an ambitious target, reliable, up-to-date and fine-grained measurements of socioeconomic indicators are necessary. When it comes to socioeconomic development, novel digital traces can provide a complementary data source to overcome the limits of traditional data collection methods, which are often not regularly updated and lack adequate spatial resolution. In this study, we collect publicly available and anonymous advertising audience estimates from Facebook to predict socioeconomic conditions of urban residents, at a fine spatial granularity, in four large urban areas: Atlanta (USA), Bogotá (Colombia), Santiago (Chile), and Casablanca (Morocco). We find that behavioral attributes inferred from the Facebook marketing platform can accurately map the socioeconomic status of residential areas within cities, and that predictive performance is comparable in both high and low-resource settings. Our work provides additional evidence of the value of social advertising media data to measure human development and it also shows the limitations in generalizing the use of these data to make predictions across countries.

Comparing sources of mobility for modelling the epidemic spread of Zika virus in Colombia.

Timely, accurate, and comparative data on human mobility is of paramount importance for epidemic preparedness and response, but generally not available or easily accessible. Mobile phone metadata, typically in the form of Call Detail Records (CDRs), represents a powerful source of information on human movements at an unprecedented scale. In this work, we investigate the potential benefits of harnessing aggregated CDR-derived mobility to predict the 2015-2016 Zika virus (ZIKV) outbreak in Colombia, when compared to other traditional data sources. To simulate the spread of ZIKV at sub-national level in Colombia, we employ a stochastic metapopulation epidemic model for vector-borne diseases. Our model integrates detailed data on the key drivers of ZIKV spread, including the spatial heterogeneity of the mosquito abundance, and the exposure of the population to the virus due to environmental and socio-economic factors. Given the same modelling settings (i.e. initial conditions and epidemiological parameters), we perform in-silico simulations for each mobility network and assess their ability in reproducing the local outbreak as reported by the official surveillance data. We assess the performance of our epidemic modelling approach in capturing the ZIKV outbreak both nationally and sub-nationally. Our model estimates are strongly correlated with the surveillance data at the country level (Pearson's r = 0.92 for the CDR-informed network). Moreover, we found strong performance of the model estimates generated by the CDR-informed mobility networks in reproducing the local outbreak observed at the sub-national level. Compared to the CDR-informed networks, the performance of the other mobility networks is either comparatively similar or substantially lower, with no added value in predicting the local epidemic. This suggests that mobile phone data captures a better picture of human mobility patterns. This work contributes to the ongoing discussion on the value of aggregated mobility estimates from CDRs data that, with appropriate data protection and privacy safeguards, can be used for social impact applications and humanitarian action.