COVID-19 Flow-Maps an open geographic information system on COVID-19 and human mobility for Spain.

COVID-19 is an infectious disease caused by the SARS-CoV-2 virus, which has spread all over the world leading to a global pandemic. The fast progression of COVID-19 has been mainly related to the high contagion rate of the virus and the worldwide mobility of humans. In the absence of pharmacological therapies, governments from different countries have introduced several non-pharmaceutical interventions to reduce human mobility and social contact. Several studies based on Anonymized Mobile Phone Data have been published analysing the relationship between human mobility and the spread of coronavirus. However, to our knowledge, none of these data-sets integrates cross-referenced geo-localised data on human mobility and COVID-19 cases into one all-inclusive open resource. Herein we present COVID-19 Flow-Maps, a cross-referenced Geographic Information System that integrates regularly updated time-series accounting for population mobility and daily reports of COVID-19 cases in Spain at different scales of time spatial resolution. This integrated and up-to-date data-set can be used to analyse the human dynamics to guide and support the design of more effective non-pharmaceutical interventions.

Background & Summary

COVID-19 is an infectious disease caused by SARS-CoV-2 virus, which has spread all over the world, leading to a global pandemic state. The fast progression of COVID-19 has been mainly related to two crucial factors, the high contagion rate of the virus, as well as, the worldwide mobility of humans. The high contagious rate has represented one of the leading causes of healthcare systems collapse in several countries, such as Italy and Spain, contributing to a large number of deaths. As a result, and in the absence of vaccines or other pharmacological therapies, governments from different countries have introduced several non-pharmaceutical interventions (NPI) to reduce human mobility and social contact[1]. Different NPIs include, but are not limited to, the closing of national borders, the temporal partial or total lockdowns, and social distancing measures with the aim of attenuating the epidemic[2].

Several studies based on Anonymized Mobile (cell) Phone Data (AMPD) have been published analysing the relationship between human mobility and the spread of coronavirus[3]. For instance, in a study conducted in the USA, the authors reported a high correlation between human mobility and COVID-19 transmission rate[4]. A similar study conducted in Italy reports that mobility restrictions and social distancing policies could reduce contagions up to 45%[5]. AMPD has also been used to study how the different NPIs can change the underlying structure of the human mobility networks and how those changes affect the spread of the disease[6,7].

AMPD has also been integrated into epidemiological models to simulate the spread of SARS-CoV-2 in time and space. For instance, a model integrating human mobility was used to predict the first wave of COVID-19 cases in Spain and to call the government to enforce a full lockdown[8,9]. Similarly, models have also been used to study the impact of mobility restrictions, physical distancing, as well as the impact of lockdown on social inequalities[10-12]. Models are critical in the design of NPIs to effectively mitigate the pandemic while reducing the negative impact on the economy, education, and other social activities[13,14]. However, developing realistic models to guide policy-makers decisions requires high quality and up-to-date data[15]. Ideally, data records should be open access, de-identified, periodically updated, and should include data on human mobility, COVID-19 daily cases, and contact tracing, among other sources of information. Furthermore, data records should be properly geo-localized with the highest possible degree of spatio-temporal resolution[16].

Several data-sets on human mobility have been made publicly available to help analyse and model the epidemic dynamics. Pepe et al., have released a data-set based on 170,000 de-identified and aggregated smartphone users that account for movements between Italian provinces[17]. Another notable data-set includes a regularly updated multi-scale human mobility dynamic network across the United States[18]. Furthermore, data-sets on COVID-19 cases across the globe have been also released[19,20]. However, to our knowledge, none of the cited data-sets integrates cross-referenced geo-localized data on human mobility and COVID-19 cases into one all-inclusive open resource.

Herein we present a cross-referenced geographic information system (GIS) named COVID-19 Flow-Maps to manage, retrieve, visualise and analyse regularly updated time-series data on population mobility networks and daily reports of COVID-19 cases in Spain (Fig. 1). Human mobility data is provided as Origin-Destination (OD) matrices at different levels of spatial and temporal resolution (district, province, and autonomous community on a daily and hourly basis). Every mobility data record reports the number of people performing zero, one, two, or more trips, on a daily basis. Health data consisting of daily confirmed COVID-19 cases is also provided on a daily basis and at different levels of spatial resolution. All the data has been gathered from official access points, as detailed in the sec:data-rec section. We offer provenance records that track the origin of the data and information regarding all the processing steps in those cases where the data has been consolidated. All data records are regularly-updated and are accessible through an Application Programming Interface (API) as well as through a dedicated GitHub repository. We also provide an interactive web application for exploring the data. In this work, we present the system and illustrative examples to show how the combination of mobility and health data can help public health directives to effectively mitigate COVID-19 transmission.

Fig. 1

Graphical representation of COVID-19 Flow-Maps Geographic Information Systems. The main data records include geographical layers for different territorial units, COVID-19 daily cases reported at different spatial resolution and phone-based anonymized mobility data in the form of daily origin-destination matrices. All the information is stored in a non-SQL database that can be directly queried through a REST-API, downloaded using provided scripts, and accessed through web-based interactive data dashboards.

Methods

Data acquisition

The COVID-19 data records for Spain are retrieved from endpoints provided by a variety of official sources including national-level reports, as well as those reported by autonomous community governments (see COVID-19 data and Data Records section and Online-only Table 1). As data is obtained from different sources, each record differs in its format, including different fields, field names, and different levels of aggregation (e.g. age, gender). Additionally, data retrieved from different sources are mapped into different geographical layers (e.g. autonomous communities, provinces, municipalities). Therefore, in addition to the COVID-19 reported cases, the spatial information needed for its mapping is also retrieved from the corresponding official sources (see Online-only Table 2). Population mobility records based on AMPD are reported by the Spanish Ministry of Transport, Mobility and Urban Agenda (MITMA, Ministerio de Transportes, Movilidad y Agenda Urbana) in tabular format and it includes two different daily reports, as well as the geographical layer needed to map the data (for further details see Mobility data in Data Records section).

Online-only Table 1

Description of COVID-19 data-sets reported for different geographical layer.

Region	Provider	Initial date	Final date	Update rate	Information provided	Geographic detail	Cases report type	Sex-disaggregated	Age-disaggregated	Event ID	Layer ID	Home page	Link
Provinces	Centro Nacional de Epidemiología	01-01-2020	Present	Daily	Total cases divided by PCR and antigens. Elisa and antibodies tests.	Provinces	Daily	X	X	ES.covid_cpro	cnig_provincias	https://cnecovid.isciii.es/covid19/	https://cnecovid.isciii.es/covid19/resources/casos_diagnostico_provincia.csv
Autonomous communities	Centro Nacional de Epidemiología	01-01-2020	Present	Daily	Total cases divided by PCR, antigens. Elisa and antibodies test.	Autonomous communities	Daily			ES.covid_cca	cnig_ccaa	https://cnecovid.isciii.es/covid19/	https://cnecovid.isciii.es/covid19/resources/casos_diagnostico_ccaa.csv
Cataluña	Generalitat de Catalunya	01-03-2020	Present	Daily	Total cases divided by PCR and antigens. Negative cases. Suspect cases.	Basic Health Areas and Municipalities	Daily	X		09.covid_abs	abs_09	https://analisi.transparenciacatalunya.cat/browse?q=covid&sortBy=relevance	https://analisi.transparenciacatalunya.cat/resource/xuwf-dxjd.json?%24order=data+DESC%2C+abscodi+DESC%2C+absdescripcio+DESC%2C+resultatcoviddescripcio+DESC%2C+sexecodi+DESC&%24limit=181604&%24%24app_token=kuJ8bqzm2w4JgFk62r5wbfYXR
Navarra	Gobierno de Navarra	25-03-2020	Present	Daily	Total cases.	Basic Health Areas and Municipalities	Daily			15.covid_abs	zbs_15		http://www.navarra.es/appsext/DescargarFichero/default.aspx?codigoAcceso=OpenData&fichero=coronavirus\\datos_coronavirus_zonas_basicas.csv http://www.navarra.es/appsext/DescargarFichero/default.aspx?codigoAcceso=OpenData&fichero=coronavirus\\datos_zonas_basicas_salud_poblacion.csv
Castilla y León	Junta de Castilla y León	29-02-2020	Present	Daily	Total cases. Total and positive PCRs. People with symptoms.	Basic Health Areas	Daily			07.covid_abs	zbs_07	https://datosabiertos.jcyl.es/web/es/datos-abiertos-castilla-leon.html	https://analisis.datosabiertos.jcyl.es/explore/dataset/tasa-enfermos-acumulados-por-areas-de-salud/download/?format=json
País Vasco	Gobierno Vasco	21-03-2020	Present	Daily	Total cases.	Basic Health Areas and Municipalities	Daily			16.covid_abs	oe_16		https://opendata.euskadi.eus/contenidos/ds_informes_estudios/covid_19_2020/opendata/aggregated/json/osasun_eremuak-zonas_salud-by_date.json
Madrid	Comunidad de Madrid	27-02-2020	Present	Every two days	Total cases.	Basic Health Areas and Municipalities	Accumulated			13.covid_abs	zon_bas_13	https://datos.comunidad.madrid/catalogo/dataset?tags=COVID-19	https://datos.comunidad.madrid/catalogo/dataset/b3d55e40-8263-4c0b-827d-2bb23b5e7bab/resource/01a7d2e8-67c1-4000-819d-3356bb514d05/download/covid19_tia_zonas_basicas_salud.json https://datos.comunidad.madrid/catalogo/dataset/b3d55e40-8263-4c0b-827d-2bb23b5e7bab/resource/907a2df0-2334-4ca7-aed6-0fa199c893ad/download/covid19_tia_zonas_basicas_salud_s.json
Comunidad Valenciana	Generalitat Valenciana	27-05-2020	Present	Twice per week	Total cases. Deceased.	Municipalities	Accumulated			10.covid_cumun	cnig_municipios	https://dadesobertes.gva.es/es/dataset?tags=COVID-19	https://dadesobertes.gva.es/dataset/15810be9-d797-4bf3-b37c-4c922bee8ef8/resource/609a937e-252c-43eb-948a-5766dbcd0c18/download/covid-19-casos-confirmados-por-pcr-casos-pcr-en-los-ultimos-14-dias-y-personas-fallecidas-por-mu.csv
Cantabria	Servicio Cántabro deSalud	30-03-2020	#########		Total cases. Deceased.	Municipalities	Accumulated			06.covid_cumun	cnig_municipios	https://www.scsalud.es/coronavirus	https://www.scsalud.es/documents/2162705/9255280/2020_covid19_municipal.xlsx
Asturias	Dirección General deSalud Pública deAsturias	29-02-2020	Present	Daily	Total cases divided by PCR and antigens (positive and negative). Deceased. Hospital and ICU occupation, discharge. Traceability of positive cases.	Basic Health Areas and Municipalities	Daily		X	03.covid_cumun	cnig_municipios	https://obsaludasturias.com/obsa/niveles-covid-19/	https://dgspasturias.shinyapps.io/panel_de_indicadores_asturias/DATOS/TABLAS_RESUMEN/CONCEJOS/concejos_resumen.zip

Online-only Table 2

COVID- 19 Flow-Maps Geographical Layer used to georeference data-sets.

Map	Region	Provider	Layer	Home page	Link
World map	World	Geographic Information System of the Cоmmission (GISCO)	world	https://ec.europa.eu/eurostat/web/gisco/geodata/reference-data/administrative-units-statistical-units/countries	http://gisco-services.ec.europa.eu/distribution/v2/countries/download/ref-countries-2020-03m.shp.zip
INE mobility areas	Spain	Instituto Nacional de Estadística (INE)	ine_mov	https://www.ine.es/covid/covid_movilidad.htm	https://www.ine.es/covid/shapefiles_celdas_marzo2020.zip
MITMA mobility areas	Spain	Ministerio de Transporte, Movilidad y Agenda Urbana (MITMA)	mitma_mov	https://www.mitma.gob.es/ministerio/covid-19/evolucion-movilidad-big-data/opendata-movilidad
INE census sections	Spain	INE	ine_sec	https://www.ine.es/ss/Satellite?L=es_ES&c=Page&cid=1259952026632&p=1259952026632&pagename=ProductosYServicios%2FPYSLayout	https://www.ine.es/prodyser/cartografia/seccionado_2020.zip
INE districts	Spain	INE	ine_districts	https://www.ine.es/ss/Satellite?L=es_ES&c=Page&cid=1259952026632&p=1259952026632&pagename=ProductosYServicios%2FPYSLayout	https://www.ine.es/prodyser/cartografia/seccionado_2020.zip
CNIG municipies	Spain	Centro Nacional de Información Geográfica (CNIG)	cnig_municipios	http://centrodedescargas.cnig.es/CentroDescargas/catalogo.do?Serie=LILIM
CNIG provinces	Spain	CNIG	cnig_provincias	http://centrodedescargas.cnig.es/CentroDescargas/catalogo.do?Serie=LILIM
CNIG autonomous communities	Spain	CNIG	cnig_ccaa	http://centrodedescargas.cnig.es/CentroDescargas/catalogo.do?Serie=LILIM
Primary attention centers	Spain	ESRI España	ate_pri	https://comunidadsig.maps.arcgis.com/home/item.html?id=e9e8543a13474b4abef15f6d5a43ab31	https://services1.arcgis.com/nCKYwcSONQTkPA4K/arcgis/rest/services/CentrosAtenci%C3%B3nPrimaria/FeatureServer
Outpatient Urgent Care Centers	Spain	ESRI España	ate_urg_extra	https://comunidadsig.maps.arcgis.com/home/item.html?id=c4ad72ed81594da4955acb2a57fd8533	https://services1.arcgis.com/nCKYwcSONQTkPA4K/arcgis/rest/services/Centros_de_Atenci%C3%B3n_Urgente_Extrahospitalaria/FeatureServer
Hospitals	Spain	ESRI España	hospitales	https://comunidadsig.maps.arcgis.com/home/item.html?id=68745a7fb7a348b6b0d722c8517790af	https://services1.arcgis.com/nCKYwcSONQTkPA4K/arcgis/rest/services/Hospitales/FeatureServer
Andorra state boundary	Andorra	Instituto de Estudios Andorranos (IEA)	andorra	https://www.iea.ad/limit-estatal	https://www.iea.ad/images/stories/Sigma/Mapes_E00/Limit_Estat_Andorra_2012.rar
GenCat Sanitary regions	Cataluña	Generalitat de Cataluña	reg_san_09	https://salutweb.gencat.cat/ca/el_departament/estadistiques_sanitaries/cartografia/	https://salutweb.gencat.cat/web/.content/_departament/estadistiques-sanitaries/cartografia/ABS_2018.zip
GenCat Sanitary sectors	Cataluña	Generalitat de Cataluña	sect_san_09	https://salutweb.gencat.cat/ca/el_departament/estadistiques_sanitaries/cartografia/	https://salutweb.gencat.cat/web/.content/_departament/estadistiques-sanitaries/cartografia/ABS_2018.zip
GenCat healthcare management areas	Cataluña	Generalitat de Cataluña	aga_09	https://salutweb.gencat.cat/ca/el_departament/estadistiques_sanitaries/cartografia/	https://salutweb.gencat.cat/web/.content/_departament/estadistiques-sanitaries/cartografia/ABS_2018.zip
GenCat basic health areas	Cataluña	Generalitat de Cataluña	abs_09	https://salutweb.gencat.cat/ca/el_departament/estadistiques_sanitaries/cartografia/	https://salutweb.gencat.cat/web/.content/_departament/estadistiques-sanitaries/cartografia/ABS_2018.zip
Cataluña shires	Cataluña	Generalitat de Cataluña	shires_09	https://salutweb.gencat.cat/ca/el_departament/estadistiques_sanitaries/cartografia/	http://descarregues.icgc.cat/descarregues2/dldirect.php?p=bm5m&d=shp
Municipies and districts of Madrid	Madrid	Datos abiertos de la Comunidad de Madrid	mun_dis_13	http://origin-datos.comunidad.madrid/catalogo/organization/comunidad-de-madrid	https://datos.comunidad.madrid/catalogo/dataset/7da43feb-8d4d-47e0-abd5-3d022d29d09e/resource/04f71acb-e0c1-47be-84c9-1f57d2e09e40/download/municipios_y_distritos_madrid.zip
Madrid basic zones of Sanity	Madrid	Datos abiertos de la Comunidad de Madrid	zon_bas_13	http://origin-datos.comunidad.madrid/catalogo/dataset/b3d55e40-8263-4c0b-827d-2bb23b5e7bab/resource/f1837bd3-a835-4110-9bbf-fae06c99b56b	https://datos.comunidad.madrid/catalogo/dataset/b3d55e40-8263-4c0b-827d-2bb23b5e7bab/resource/f1837bd3-a835-4110-9bbf-fae06c99b56b/download/zonas_basicas_salud.zip
Madrid health areas	Madrid	Instituto de estadística de la Comunidad de Madrid	area_san_13	http://www.madrid.org/nomecalles/DescargaBDTCorte.icm
Health districts of Madrid	Madrid	Instituto de estadística de la Comunidad de Madrid	dis_san_13	http://www.madrid.org/nomecalles/DescargaBDTCorte.icm
Health zones of Madrid	Madrid	Instituto de estadística de la Comunidad de Madrid	zon_san_13	http://www.madrid.org/nomecalles/DescargaBDTCorte.icm
Basic health zones of Madrid	Madrid	Instituto de estadística de la Comunidad de Madrid	zbs_13	http://www.madrid.org/nomecalles/DescargaBDTCorte.icm
Health zones of Euskadi	País Vasco	Gobierno Vasco	oe_16	https://www.geo.euskadi.eus/zonas-de-salud-de-euskadi/s69-geodir/es/	ftp://ftp.geo.euskadi.eus/cartografia/Salud/Zonas_Salud/OsasunEremuak_ZonasSalud_A_2018.zip
Health centers, outpatient clinics and public hospitals of País Vasco	País Vasco	Gobierno Vasco	cs_hosp_16	https://opendata.euskadi.eus/catalogo/-/centros-de-salud-publicos-en-euskadi/	https://opendata.euskadi.eus/contenidos/ds_localizaciones/centros_salud_en_euskadi/opendata/centros-salud.geojson
Basic health zones of Navara	Navarra	Servicio Navarro de Salud	zbs_15	https://gobiernoabierto.navarra.es/es/open-data/datos/servicio-navarro-salud-delimitacion-zonas-basicas-salud	https://idena.navarra.es/descargas/DOTACI_Pol_SNSZonas.zip
Basic health areas of Navarra	Navarra	Servicio Navarro de Salud	as_15	https://gobiernoabierto.navarra.es/es/open-data/datos/servicio-navarro-salud-delimitacion-areas-salud	https://idena.navarra.es/descargas/DOTACI_Pol_SNSAreas.zip
Basic health zones of Castilla y León	Castilla y León	Junta de Castilla y León	zbs_07	https://datosabiertos.jcyl.es/web/jcyl/set/es/medio-ambiente/zonas-basicas-salud-cyl/1284688135689	https://datosabiertos.jcyl.es/web/jcyl/risp/es/medio-ambiente/zonas-basicas-salud-cyl/1284688135689.shp
Galicia shires	Galicia	Junta de Galicia	comarcas_12	http://mapas.xunta.gal/centro-de-descargas	http://visorgis.cmati.xunta.es/cdix/descargas/visor_basico/Comarcas.zip
Galicia municipies	Galicia	Junta de Galicia	concellos_12	http://mapas.xunta.gal/centro-de-descargas	http://visorgis.cmati.xunta.es/cdix/descargas/visor_basico/Concellos_IGN.zip
Galicia parishes	Galicia	Junta de Galicia	parroquias_12	http://mapas.xunta.gal/centro-de-descargas	http://visorgis.cmati.xunta.es/cdix/descargas/visor_basico/Parroquias.zip
Galicia health areas	Galicia	Servicio Gallego de Salud	as_12	https://www.sergas.es/Saude-publica/GIS-Areas-sanitarias	https://www.sergas.es/Saude-publica/Documents/5916/Areas_Sanitarias.zip
Andalucía sanitary districts	Andalucía	Instituto de Estadística y Cartografía de Andalucía	dis_san_01	https://www.juntadeandalucia.es/institutodeestadisticaycartografia/DERA/g13.htm	https://www.juntadeandalucia.es/institutodeestadisticaycartografia/DERA/ficheros/13_Limites_Administrativos.zip
Andalucía basic health zones	Andalucía	Instituto de Estadística y Cartografía de Andalucía	zbs_01	https://www.juntadeandalucia.es/institutodeestadisticaycartografia/DERA/g13.htm	https://www.juntadeandalucia.es/institutodeestadisticaycartografia/DERA/ficheros/13_Limites_Administrativos.zip
Asturias parishes	Asturias	Gobierno del Principado de Asturias	parroquias_03	https://ideas.asturias.es/es/centro-de-descargas?category=1&subcategory=7&lang=es	https://www.juntadeandalucia.es/institutodeestadisticaycartografia/DERA/ficheros/13_Limites_Administrativos.zip
Asturias councils	Asturias	Gobierno del Principado de Asturias	concejos_03	https://ideas.asturias.es/es/centro-de-descargas?category=1&subcategory=7&lang=es	https://www.juntadeandalucia.es/institutodeestadisticaycartografia/DERA/ficheros/13_Limites_Administrativos.zip
Asturias health areas	Asturias	Gobierno del Principado de Asturias	as_03	https://ideas.asturias.es/es/centro-de-descargas?category=1&subcategory=7&lang=es	https://www.juntadeandalucia.es/institutodeestadisticaycartografia/DERA/ficheros/13_Limites_Administrativos.zip

Data processing and consolidation

To automatically update the COVID-19 Flow-Maps data records we have implemented a workflow that retrieves data from different endpoints, processes and stores new records (see Code availability section). The first step in this workflow is to validate new entries by checking inconsistent and missing values. Secondly, the system controls for duplicated entries ensuring that for any given territorial unit on a given date a single reported value exists. Finally, for all the integrated data sources, the original entries, as well as the registry of any modification, are also stored in the provenance collection (see Data Provenance subsections). For each geographical layer needed to geo-localise each data entry we stored all the geometries, translating the Geographic Coordinates on ETRS89 Datum (EPSG:4258), and assigning them a symbolic and unique layer name (see Geographical Layers in Data Records section).

COVID-19 cases reports

When integrating a source that reports COVID-19 cases, three key attributes are first identified, namely, the date for the reported event, the reported value (e.g. daily/accumulated incidence) and any attribute that can be used to geo-reference the data. Geo-referencing attributes can be geographic coordinates or, more usually, an identifier that matches a geometry of a defined geographical layer, e.g. the unique identifier of a municipality and a reference to that layer. For instance, the government of some autonomous communities reports COVID-19 cases on a daily basis at the level of Basic Health Areas (BHA) whereas others report cases by municipalities (see COVID-19 data and Data Records section). Missing values of the daily incidence for a given date or geographic area are imputed with zero values, while missing values of cumulative incidence are imputed based on the previous available value. After validation, since some data sources either report daily incidence or cumulative incidence, we also calculate one attribute from the other. Moreover, field names are normalised while also keeping the original field names and values to be queried. Additionally, commonly used metrics such as accumulated incidences over one and two weeks are also calculated and stored.

Reconstruction of mobility networks

The reconstruction of mobility networks relies on two main sources of information, namely: recorded events data and mobile phone network topology data. The former corresponds to anonymized data associated with the connection records of the mobile devices with the mobile phone network. These records include both active and passive events. Active events are made up of CDRs (Call Detail Records) that provide a record every time a device interacts with the network (calls, sending text messages, data sessions). These records are joined with passive events data (periodic update of the device position, changes in coverage areas, etc.), providing a very high temporal granularity. Regarding spatial resolution, location information is available at the level of the coverage area of each antenna - which implies a spatial resolution of tens or hundreds of meters in the city and up to several kilometres in rural areas, which provides an idea of the uncertainty that is introduced in the determination of the position according to the areas analysed.

The recorded events data are processed in a secure environment within the mobile operator’s infrastructure where this data is aggregated and anonymized in compliance with the existing European and Spanish legislation, e.g. LOPD-ODD (Ley de Protección de Datos). The phone network data includes the location of the communication towers which are used to obtain a Voronoi Diagram of the cellular coverage map. Additional sources used include land use information from the Spanish Land Use Information System (SIOSE), population data (i.e. the register of inhabitants by districts) and the Spanish transportation network (e.g. airports location, rail network, etc.).

The first step in the data processing workflow consists of the extraction and pseudonymization of the mobile (cell) phone records. The pseudonymization of the records is based on the use of a one-way hash function, that is, a function that allows the calculation of an anonymized identifier from the original identifier (usually the International Mobile Subscriber Identity) in such a way that it is impossible to carry out the reverse process. Furthermore, a perfect hash function is used to avoid collisions, i.e. that two different original identifiers produce the same anonymized identifier. The anonymized phone records are stored in a secure environment within the infrastructure of the mobile operator, where the algorithm used to aggregate the anonymized data will also run.

The processing and analysis of the raw data can be divided into the following steps[21]: (i) pre-processing and cleaning of the data; (ii) construction of the sample; (iii) identification of the place of habitual residence and the place of overnight stay; (iv) extraction and inferences of activities and trips; (v) extrapolation of sample results; and (vi) generation of indicators. First of all, a pre-processing of the mobile phone data is carried out to facilitate its management, ordering and grouping the records in the most convenient way for further analysis. A data integrity check is also carried out to eliminate possible errors in the mobile operator’s data. The process of data cleaning and debugging of errors is essential to ensure the quality of the data, preventing possible source of errors from distorting the results obtained with the algorithms used for the extraction of activity and mobility patterns.

To construct the sample, a selection of valid users (i.e. users who have enough telephone activity to infer their activity and travel diary during the entire analysis day), is first made to provide information on their movements. This selection is made according to different criteria related to their phone activity so that this is sufficient to establish their behaviour patterns with an adequate level of reliability. The unit to measure mobility is the trip. To define the concept of trip we need to introduce two other previous concepts. The first concept is the stay which is the permanence of an individual for more than 20 min in the same mobility area and the concept of activity, which is the reason that motivates the individual to travel there and includes work, home and other. Once defined these two concepts, a trip is defined as the displacement of a person between two consecutive activities. This definition is closely linked to the data source, which are the active and passive events generated in each antenna. Spatially, the detection of the trips is linked to the distance between the towers in which the antennas are located. This distance is variable from some meters in urban areas to some kilometres in rural areas. In some places, the separation among towers may not allow differentiating movements of less than 500 m since only if the terminal changes its connection to a different antenna, the movement can be detected. Because of that, to have homogeneous data, we filter the movements counting as trips for only those of more than 500 meters.

To identify the area of most usual home areas of the users, an analysis of the behavioural patterns is conducted over several weeks. The overnight stay areas are identified in a similarly. To identify activities and trips, we use a combination of criteria based on the duration of stay times, travel itineraries and patterns of behaviour throughout the study period. The information associated with each activity includes the location, the start time and the end time. In particular, intermediate stops subordinate to the journey and made between stages of the journey (e.g. intermediate stops for a transfer at a bus/train station) are filtered to consolidate all the stages as a single trip. In this way, a person travelling from their home to their working place that has to commute from the subway to a bus or train and has to wait in the station is considered a single trip from the origin (home) to the final destination (working place). Thus, by filtering intermediate stops subordinate to the journey and made between stages of the journey (e.g. intermediate stops for a transfer at a bus/train station) the result is a sequence of activities and trips made on the analysed days. As associated information, each trip includes origin (location of the activity immediately before the trip), destination (location of the activity immediately after the trip), the start time of the trip (time of the end of the previous activity) and end time (start time of the next activity), distance, home place, activity in the origin and activity in the destination.

The extrapolation of the sample is carried out by taking as the sampling reference the population residing in the country, according to the data from the Population Register provided by the National Institute of Statistics. Standard sample extrapolation procedures are used (similar to those used, for example, in a household mobility survey), applying expansion factors by area of residence based on the sample/population ratio for each age and gender segment in each district. Finally, the information obtained in the previous steps is presented with the required spatial and temporal resolution and the desired segmentation to generate the origin-destination matrices and the rest of the mobility indicators previously described. The generation of the mobility indicators has been carried out using specialised software developed for this purpose and all these processes are carried out within the mobile operator’s infrastructure, so that the information generated and delivered to our source, the MITMA official site for data release, is already aggregated and anonymized.

The output obtained from the processing steps is then used for generating two mobility indicators: the hourly-based OD matrices referred to as the Maestra 1, and the daily population mobility descriptions referred as the Maestra 2 (see Mobility:data in the Data Records section). Both indicators are geo-referenced to a custom layer referred to as the MITMA mobility layer (see Geographical Layers in the Data Records section). Further details on the analysis and processing of mobility data are provided on the official page of the study (https://www.mitma.gob.es/ministerio/covid-19/evolucion-movilidad-big-data).

Mobility indicators are further processed to obtain more aggregated mobility metrics. Hourly-based OD matrices (Maestra 1) are aggregated to obtain daily mobility matrices. Furthermore, daily population mobility indicators (Maestra 2) that provide the number of people that have done zero, one, two, and more than two trips, are aggregated within each mobility area in a given date to estimate the total population of each mobility area on a daily basis. The total population computed in this way is stored to be used in the calculation of other descriptors (e.g. daily incidence by total population). The population inferred from Maestra 2 is estimated on a daily basis and thus it captures the fluctuations in the population due to net mobility over the year (e.g. mobility during summer vacations).

Data provenance

The database has all the information needed to identify the origin of the data, all the processing carried out, the original files retrieved, and the timestamp of the last update. Furthermore, copies of all the data obtained from the different sources are kept in their original formats and their source URL (if available). All this information can also be queried through the REST-API (see REST-API in the Usage Notes section).

Data projection using geographical layers overlay

In order to combine COVID-19 daily incidence and population mobility data, both data records should be projected onto the same geographical layer. In some cases, one Polygons in one layer must be covered by a single polygon from another layer, with an exclusive overlap (e.g. municipalities are included in provinces). In other cases, the overlap between the two layers is not exclusive, which means that polygons in one layer can be covered by more than one polygon from another layer. For instance, COVID-19 daily incidences and mobility data are geo-referenced into different geographical layers that cannot be combined directly. Thus, to overcome this issue, we have implemented a general approach to project data among different geographical layers. The approach is based on linear interpolation over the overlaying areas between the polygons from the two layers. We call this process “projecting” data from a source layer (e.g. municipalities) into the target layer (e.g. BHAs).

Spatial Overlay Matrix

In general, a geographical layer is composed of several different polygons. For example, in the case of the province layer, each individual province will correspond to a polygon that defines its geographic frontier or border. Therefore, given two layers with polygons and with polygons, we define the Spatial Overlay Matrix where contains the proportion of the polygon that falls into the polygon . In general, and all the rows in the overlay matrix sum up to 1:

Population-based overlay matrix

Using Spatial-based overlays has the implicit assumption that the population is distributed uniformly on the territory. Nonetheless, this assumption may not hold in most cases, and therefore, some authors have proposed using population density grids to improve the projection’s accuracy. The idea is to account for how the population is distributed among the territory and use the population overlays, i.e. the number of people each area of the first layer share with each area of the second layer. Therefore this approach requires information on how the population is distributed with high spatial resolution. Herein we use the population density grid provided by GEOSTAT[22] that accounts for the population distribution of all Europe. This grid has cells of 1 km², and each cell has assigned an estimated population value.

To calculate the population-based overlays between layer A and a second layer B using the population density grid G we first calculate the intersections between A and B obtaining a set of intersections (as in the case of the spatial overlay). Then, for each intersection, we calculate its population by performing a second intersection with the population grid. The population of each intersection is calculated as the sum of the population of the cells that fit completely in the intersection, and the fraction of those cells that overlay but do not fit completely. Once we have the population of each intersection, we build the overlay matrix using population values instead of areas. Finally, as in the case of the Spatial-based overlays, the rows of the overlay matrix are normalised so that each row sums one (see Eq. 1).

Data projection using the overlay matrix

Given a row vector of data on layer A (e.g. COVID-19 incidence value for each polygon of layer A), this data can be projected into the layer B by just multiplying and the overlay matrix:

Figure 2a shows an example of how cases geo-referenced in an origin layer A can be projected into a destination layer B using a pre-calculated overlay matrix. Interestingly, this approach does not depend on how the overlay matrix was constructed, i.e. if it is based on geographic or population-based overlays.

Fig. 2

Toy example to explain the approach for projecting data between layers using Spatial-based overlays. Panel a shows an example of cases projection from layer A to layer B, using the Spatial-based overlays between both layers. Panel b shows an example of trips projection between the same layers.

The same approach is also used to project an OD mobility matrix between geographical layers. Projecting OD matrices to other geographical layer allow the integration of population mobility with data-set reported in a different layer, e.g. one on which COVID-19 cases are reported. In general, given a mobility matrix geo-referenced to the layer A, each entry contains the number of trips from polygon to polygon , (both polygons from layer A) it can be projected to the layer B by multiplying it by the overlay matrix and its transpose:Such calculation can be seen as first projecting the origins and then projecting the destinations (2-b). Overlay matrices for all the combinations of geographical layers have been computed and stored in the database, enabling a fast projection of any data-set between the different geographical layers.

Mobility associated risk

To assess the effect of population mobility on the spreading of COVID-19, we have developed a risk score named Mobility Associated Risk (MAR). The MAR score integrates daily cases with populations flows between different geographical areas (i.e. OD matrices), e.g. provinces or BHAs, to estimate how likely it is that cases spread as a consequence of human mobility (see Fig. 3). Herein, we use the incidence accumulated over two weeks as an estimator of the number of active cases. Then, for a given geographical layer L with n zones , and a given date t, we refer to the cases accumulated over two weeks as . This time window is a proxy for inferring the current number of active cases and it also reduces the potential noise due to delays in the reports. Nonetheless, is worth noting that if cases are reported with a delay larger than the size of the time window used, the estimated number of active cases would be erroneous. The estimator of the total number of cases, , is then normalised to the total population reported in zone j at t, N(t).where is the estimator of the active cases per total number of inhabitants. We then combine the from each zone j together with a daily OD mobility matrix for the date t as follows:where is the number of trips from j to k with both values reported at date t, and is an estimation of the expected number of infected subjects also moving from zone j to k at day t. In general, when the daily COVID-19 cases and the mobility are reported in the same geographical layer L, the risk can be calculated for all pairs of zones j and k by the element-wise, or Hadamard product between the vector of cases densities and the mobility matrix :

Fig. 3

Mobility Associated Risk chart. Panel a is a graphic representation of how the mobility associated risk between zones j and k is calculated from the normalised accumulated incidence and the number of trips. Panels b and c show an example of the total incoming and outgoing risk for the province of Barcelona (highlighted in yellow), respectively. The colour scale on the maps from (c) and (b) indicate the incoming MAR and outgoing MAR, respectively; more intense red/violet indicate a greater incoming/outgoing MAR respect the target zone. Arrows indicate the top ten zones of incoming (b) outgoing (c) MAR.

The matrix is thus a directed weighted network where the nodes correspond to the different n zones from layer L and the flow between any pair of nodes (j, k) is the estimated number of infected subjects moving from the source j to target k at t. The network-based structure of the mobility associated risk allows the definition of the total MAR incoming to zone k by summing the risk over all possible sources (i.e. summing the column of ). In a similar way, that outgoing MAR for a given zone k corresponds to all the weighted edges having j as the source node (i.e. summing the row of ).

The risk network can be represented in the map to analyse the source of incoming MAR to any zone of interest. For instance, the risk network can be calculated between provinces by combining daily incidence reported at province level together with mobility data aggregated at the same level. Figure 3b,c shows a representation of the incoming and outgoing mobility-associated risk for the province of Barcelona for the 10th of October 2020. The incoming risk represents the expected number of imported cases from other provinces (Fig. 3b) whereas the outgoing risk corresponds to the expected number of infected individuals travelling to a different province (Fig. 3c). Mobility Associated Risk (MAR) can be calculated at different scales of spatial resolution by aggregating mobility zones based on different criteria.

Finally, we would like to stress that although we are aware that the MAR score is a raw approximation (e.g. infected cases in quarantine are not expected to travel with the same frequency as healthy people does), and for this reason, the interpretation of the MAR score should be taken with caution. Nonetheless, we find that the MAR score can be used as an approach to evaluate the risk of outbreaks in different zones due to imported cases.

Data visualisation

We have developed a web interactive data dashboard to visualise different metrics as well as results from different analysis pipelines using interactive maps, plots and tables. The COVID-19 Flow-Maps Boards are a set of interactive web dashboards that provide access to different reports of the health situation, the population mobility and its associated risk for the different regions of Spain. Currently, we provide access to three different interactive dashboards for the analysis and visualisation of the evolution of COVID-19 cases in Spain (at three different scales of spatial resolution); the population mobility between different municipalities (or districts, in densely populated areas); and the MAR networks between regions. The COVID-19 Flow-Maps Boards can be accessed through the following link: http://flowmaps.life.bsc.es/.

Data Records

COVID-19 Flow-Maps (http://flowmaps.life.bsc.es/) is a geographic information system that integrates two main sources of information for analysing the COVID-19 pandemic in Spain. In the first place, the system provides daily reports on COVID-19 cases for different regions and at different scales of spatial resolution; secondly, the system provides access to a daily updated data-set on population mobility which includes Origin-Destination (OD) matrices and the total number of trips per person per day. The system also provides access to the different geographical layers to which the different data-sets are geo-referenced. All data collections are stored in a distributed MongoDB database that can be queried through a REST-API and a command-line utility (see Usage Notes section). Furthermore, the data records introduced in this section can be downloaded from our periodically updated GitHub data repository (https://github.com/bsc-flowmaps).

Geographical layers

This data record includes all geographical layers on which the different data records are geo-referenced (e.g. mobility, COVID-19 cases). The different layers can be grouped into those that cover the whole territory of Spain (e.g. municipalities) and those that are restricted to a specific region (Table 1). Among those that cover the full territory of Spain, the record accounts for the first four levels of administrative division, that is, autonomous communities (Fig. 4a), provinces (Fig. 4b), municipalities and districts (layers not shown). These layers are retrieved from the National Centre of Geographic Information (see Online-only Table 2 for details on the sources). Additionally, in the case of the mobility data-set, the events are geo-referenced to mobility zones from a custom layer that is provided together with the mobility data-sets. This layer was reconstructed combining cell-phone antennas coverage areas together with districts and municipalities and includes 2850 different mobility areas (see Fig. 4c) that cover almost the entire territory of Spain. In general, each area in the MITMA mobility layer may correspond to a district or group of districts in densely populated areas, and to municipalities or groups of municipalities in regions with low population density. However, there are some rural areas not covered by cell phone antennas and thus not assigned to any mobility zone.

Table 1

Table describing the different geographical layer on which main data records are reported.

Layer name	ID	Category	Polygons
Autonomous communities	cnig_ccaa	Administrative	19
Provinces	cnig_provincias	Administrative	52
Municipies	cnig_municipios	Administrative	8212
Districts	ine_districts	Administrative	10483
Cataluña Basic Health Areas	abs_09	Sanitary	374
Madrid Basic Health Areas	zon_bas_13	Sanitary	286
Castilla y León Basic Health Areas	zbs_07	Sanitary	247
Navarra Basic Health Areas	zbs_15	Sanitary	57
País Vasco Health Areas	oe_16	Sanitary	135
Mobility Areas	mitma_mov	Custom	2850

Fig. 4

The figure represents the different geographical layers included in the database. In (a) and (b) the coloured polygons correspond to the different autonomous communities and provinces of Spain, respectively. For simplicity, the Canarias islands are represented in the bottom left box. Panel (c) represents the MITMA mobility layer and the coloured polygons correspond to individual mobility zones that match district in high-density populated areas and municipalities or groups of municipalities in less populated areas. Panel (d) represents the layers for which some autonomous communities report COVID-19 cases at a higher spatial resolution than the province level. From top to bottom and left to right the layers: Madrid’s BHAs, Cataluña’s BHAs, Valencia’s Municipalities, Cantabria’s municipalities, Castilla y León BHAs, Navarra’s BHAs, País Vasco BHAs, Asturias’ municipalities. In all the plots colours are only used for visualisation purposes.

The set of layers for specific regions corresponds to the Basic Health Areas (BHA) of different autonomous communities (see Fig. 4d). According to the national health administration of Spain, each autonomous community is divided in a set of BHAs, where each individual area corresponds to the geographic delimitation for the operating of a primary care unit. Because many autonomous communities report COVID-19 cases by BHA, we have included the corresponding layer. Table 1 only includes the geographical layer of the BHA for Cataluña, Castilla y León, Madrid, Navarra and País Vasco, because those are the regions that report COVID-19 cases at the levels of BHAs. The data record also includes other geographical layers (see COVID-19 data in the Data Records section). Geographical layers are distributed in GeoJSON format in one layer per file. In addition to the coordinates of polygons, each layer file also provides the information needed for it crosse-reference with the other data-records, including unique polygon identifiers, province and autonomous communities unique codes, as well as other useful information such as pre-calculated centroids. The released data can be retrieved from Zenodo[23].

COVID-19 data

The COVID-19 record includes daily cases for all Spain reported at the level of autonomous communities as well as provinces (see Fig. 4a,b). Each record has an associated date, the corresponding identifier of the layer and code of the region and a set of COVID-19 related fields, which include the number of cases (daily incidence) as well as the number of cases segregated by test type, i.e. PCR, antigen, antibody, ELISA or unknown. Additionally, information on the total daily hospitalisations and admissions into intensive healthcare units, reported by provinces, are also provided.

This data record also includes COVID-19 cases at a higher spatial resolution (e.g. municipalities) reported by several autonomous communities. Currently, eight out of the nineteen autonomous communities publish reports with local daily COVID-19 cases at the level of municipalities or BHAs. On the one hand, Castilla y León, Madrid, Navarra and País Vasco report COVID-19 cases at the level of local BHAs; on the other hand Asturias, Navarra and Valencia report daily access at the level of municipalities; Cataluña local government reports COVID-19 daily cases at the level of BHA as well as municipalities. Table 2 shows the different COVID-19 data-sets that are integrated into this record, together with the associated geographical layer in which the data is reported. In this way, each entry reporting COVID-19 cases includes the reporting date, the geographical layer and identifier of the specific polygon within the layer and the number of cases reported, reported as daily and cumulative incidence. Additionally, each entry also includes useful metrics, such as the rolling average of the daily incidence over one and two weeks, the population in that area and the number of cases per 100,000 inhabitants. More detailed information regarding the specific data fields reported by each different source is provided in Online-only Table 1. The released data can be retrieved from Zenodo[24].

Table 2

Table describing the different sources of COVID-19 cases and the layer in which the data is reported.

Region	Layer	ID	First record	Extra info.
Spain	Autonomous communities	ES.covid_cca	2020–01–01	✓
Spain	Provinces	ES.covid_cpro	2020–01–01	✓
Cataluña	BHA	09.covid_abs	2020–01–20	✓
Navarra	BHA	15.covid_abs	2020–03–25
Castilla y León	BHA	07.covid_abs	2020–02–29	✓
País Vasco	BHA	16.covid_abs	2020–03–21
Madrid	BHA	13.covid_abs	2020–02–27
Comunidad Valenciana	Municipalities	10.covid_cumun	2020–05–27	✓
Cantabria	Municipalities	06.covid_cumun	2020–03–30	✓
Asturias	Municipalities	03.covid_cumun	2020–02–29	✓

Many sources provide more information besides the incidence.

Mobility data

Mobility data records come from a study conducted by MITMA, which analyses the mobility and distribution of the population in Spain from February 14th 2020 to May 9th 2021. (https://www.mitma.gob.es/ministerio/covid-19/evolucion-movilidad-big-data). The study is based on a sample of more than 13 million anonymized mobile-phone lines provided by a single mobile operator whose subscribers are evenly distributed.

The data record includes two different and complementary indicators, referred to as Maestra 1 and Maestra 2, whose records are geo-referenced to the MITMA mobility layer (see Table 1). The MITMA layer is composed of 2850 zones and is described in Geographical Layers within the Data Records section. The unit to measure mobility is the trip. For any person, it is considered as a trip any movement of more than 500 meters that lasts more than 20 min (see Data Processing and Consolidation). Is worth noting that, the origin and destination zone might be the same one (e.g. a person that goes to the store nearby his/her residence).

Maestra 1: Origin-Destination matrix for the mobility layer, with hourly resolution

Each entry has a date and time period (the range between two consecutive hours), the origin and destination zones and the number of trips from origin to destination. Origin and destination zones correspond to geometries from the MITMA mobility layer and internal trips (same layer of origin and destination) are also reported. Moreover, each entry also reports the putative activity of the grouped trips (from home to work, for instance), ranked travelled distances (in kilometres) as well as the total distance covered by the aggregated trips. The activity at the origin and at the destination for the reported trips in a given entry is classified into one of the following options: home, work, or other; therefore, each data entry groups trips by origin-destination and by all the matrix segmentation (place of residence, distance, travel time, activity at the origin, activity at the destination, etc).

Maestra 2: Trips per person matrix on each mobility area on a daily basis

This indicator reports population-based daily mobility behaviour. For each date and zone from the MITMA mobility layer, the indicator reports how many persons have performed 0, 1, 2 or more than 2 trips. While the indicator does not provide the destination of the trips, it accounts for the fractions of people performing at least one trip or none, as well as the estimated total population in that zone for the given date (considering as population those persons who stay overnight in the zone on that date).

OD matrices have been projected into different geographical layers and at different spatial resolution and are provided to direct use. Table 3 describes the different origin/destination layers into which OD matrices have been projected and are provided. The released data can be retrieved from Zenodo[25]. Mobility data with hourly resolution can be retrieved from the COVID-19 Flow-Maps system using the REST-API. Regarding population data, the released data can be retrieved from Zenodo[26].

Table 3

Table describing the different layer for which OD matrices have been projected into.

data-set ID	Origin layer	Target layer	Projection Type	Frequency
mitma_mov_mitma_mov	mitma_mov	mitma_mov	None	Hourly
mitma_mov_mitma_mov	mitma_mov	mitma_mov	None	Daily
abs_09_abs_09	abs_09	abs_09	Population-based	Daily
abs_09_cnig_provincias	abs_09	cnig_provincias	Population-based	Daily
cnig_ccaa_cnig_ccaa	cnig_ccaa	cnig_ccaa	Population-based	Daily
cnig_provincias_abs_09	cnig_provincias	abs_09	Population-based	Daily
cnig_provincias_cnig_provincias	cnig_provincias	cnig_provincias	Population-based	Daily
cnig_provincias_oe_16	cnig_provincias	oe_16	Population-based	Daily
cnig_provincias_zbs_07	cnig_provincias	zbs_07	Population-based	Daily
cnig_provincias_zbs_15	cnig_provincias	zbs_15	Population-based	Daily
cnig_provincias_zon_bas_13	cnig_provincias	zon_bas_13	Population-based	Daily
oe_16_cnig_provincias	oe_16	cnig_provincias	Population-based	Daily
oe_16_oe_16	oe_16	oe_16	Population-based	Daily
zbs_07_cnig_provincias	zbs_07	cnig_provincias	Population-based	Daily
zbs_07_zbs_07	zbs_07	zbs_07	Population-based	Daily
zbs_15_cnig_provincias	zbs_15	cnig_provincias	Population-based	Daily
zbs_15_zbs_15	zbs_15	zbs_15	Population-based	Daily
zon_bas_13_cnig_provincias	zon_bas_13	cnig_provincias	Population-based	Daily
zon_bas_13_zon_bas_13	zon_bas_13	zon_bas_13	Population-based	Daily

Technical Validation

Validation of estimated population

Population values used for different calculations (e.g. daily cases per total population, mobility associated risk) are estimated using the indicator Maestra 2 from the mobility data record where the population in each MITMA mobility zones is daily reported. We use this estimation because it accounts for the population fluctuations in the different regions over the year. Herein, we compare our estimated population values to those reported in the Spanish census of 2019[27]. The population values from the Spanish census are reported by municipalities whereas our estimations are reported by MITMA mobility areas. To allow the comparison of both data-sets, we aggregate population to province-level and compare the values at this geographical level. Moreover, we compare the population reported in the census with the population estimated from Maestra 2 indicator for four different dates. We calculate the Person correlation and found coefficients greater than 0.99 in the four cases. We also calculate the mean relative difference and in the four cases we find values lower than 3.5% (see Online-only Table 3). The difference may be due to the fact that some people do not live in the city where they are registered, among other things. Furthermore, the results show that the estimated values that are more similar to those reported by the census, correspond to the population estimated for the 2020–04–15, which is just after the beginning of the lockdown. We also find that population estimates that exhibit larger deviations from the reference values are for those of the autonomous cities of Ceuta and Melilla that have a very particular population structure (see Online-only Table 3). In summary, we find that the population estimated from Maestra 2 indicator is in good agreement with those values reported in the Spanish census of 2019.

Online-only Table 3

Validation of estimated population using population values reported in the Spanish census of 2019.

	References	Date: 2020-02-15		Date: 2020-04-15		Date: 2020-06-15		Date: 2020-10-15
Province	Censal Population	Estimated Population	% Error	Estimated Population	% Error	Estimated Population	% Error	Estimated Population	% Error		Mean % Error
Araba/Álava	331549	317029	4.38	327315	1.28	322741	2.66	322939	2.597		2.728
Albacete	388167	397959	2.52	392723	1.17	395635	1.92	392701	1.168		1.697
Alacant/Alicante	1858683	1832457	1.41	1831673	1.45	1834165	1.32	1849220	0.509		1.173
Almería	716820	705128	1.63	707381	1.32	711416	0.75	710269	0.914		1.154
Ávila	157640	164975	4.65	156236	0.89	156242	0.89	157424	0.137		1.642
Badajoz	673559	685995	1.85	689806	2.41	693668	2.99	686803	1.966		2.303
Illes Balears	1149460	1104996	3.87	1112458	3.22	1118858	2.66	1111600	3.294		3.261
Barcelona	5664579	5436723	4.02	5477000	3.31	5479651	3.26	5477995	3.294		3.473
Burgos	356958	359166	0.62	357773	0.23	359512	0.72	360395	0.963		0.631
Cáceres	394151	410989	4.27	403215	2.30	406715	3.19	408629	3.673		3.358
Cádiz	1240155	1231467	0.70	1225441	1.19	1249262	0.73	1233835	0.510		0.783
Castelló/Castellón	579962	586514	1.13	579960	0.00	578093	0.32	585760	1.000		0.613
Ciudad Real	495761	520341	4.96	517203	4.33	517286	4.34	512662	3.409		4.258
Córdoba	782979	798418	1.97	799684	2.13	797259	1.82	791930	1.143		1.768
A Coruña	1119596	1098370	1.90	1103919	1.40	1101147	1.65	1107137	1.113		1.514
Cuenca	196329	205014	4.42	196533	0.10	197465	0.58	201143	2.452		1.889
Girona	771044	769282	0.23	748743	2.89	748747	2.89	742029	3.763		2.444
Granada	914678	938694	2.63	904067	1.16	907687	0.76	921672	0.765		1.329
Guadalajara	257762	256592	0.45	258591	0.32	259983	0.86	260957	1.239		0.719
Gipuzkoa	723576	697862	3.55	705870	2.45	708312	2.11	700789	3.149		2.815
Huelva	521870	533749	2.28	526287	0.85	538307	3.15	524479	0.500		1.693
Huesca	220461	252398	14.49	224472	1.82	228402	3.60	228476	3.636		5.886
Jaén	633564	649823	2.57	649113	2.45	651360	2.81	644405	1.711		2.385
León	460001	474579	3.17	471534	2.51	473228	2.88	477217	3.743		3.074
Lleida	434930	440588	1.30	424243	2.46	428136	1.56	426235	1.999		1.830
La Rioja	316798	321892	1.61	314142	0.84	318047	0.39	315726	0.339		0.795
Lugo	329587	338000	2.55	337522	2.41	339728	3.08	340332	3.260		2.824
Madrid	6663394	6328664	5.02	6355007	4.63	6334976	4.93	6318097	5.182		4.941
Málaga	1661785	1636124	1.54	1630264	1.90	1637035	1.49	1640576	1.276		1.552
Murcia	1493898	1456324	2.52	1466631	1.83	1466243	1.85	1463260	2.051		2.061
Navarra	654214	628538	3.92	634912	2.95	633543	3.16	630826	3.575		3.402
Ourense	307651	312854	1.69	314645	2.27	317514	3.21	315573	2.575		2.436
Asturias	1022800	1025721	0.29	1028144	0.52	1027078	0.42	1031971	0.897		0.531
Palencia	160980	166544	3.46	164315	2.07	165735	2.95	166030	3.137		2.905
Las Palmas	1120406	1089664	2.74	1088945	2.81	1087975	2.89	1085099	3.151		2.899
Pontevedra	942665	927479	1.61	933143	1.01	935773	0.73	933345	0.989		1.085
Salamanca	330119	341372	3.41	326383	1.13	326089	1.22	339554	2.858		2.155
Santa Cruz de Tenerife	1032983	989936	4.17	984390	4.70	983590	4.78	981577	4.977		4.657
Cantabria	581078	599096	3.10	585492	0.76	586730	0.97	592967	2.046		1.720
Segovia	153129	153032	0.06	148041	3.32	147562	3.64	145846	4.756		2.945
Sevilla	1942389	1899330	2.22	1930621	0.61	1905801	1.88	1936663	0.295		1.250
Soria	88636	88125	0.58	85181	3.90	85638	3.38	88979	0.387		2.061
Tarragona	804664	797846	0.85	792664	1.49	792295	1.54	785571	2.373		1.562
Teruel	134137	138601	3.33	134968	0.62	137478	2.49	136631	1.860		2.074
Toledo	694844	725706	4.44	722803	4.02	721957	3.90	718190	3.360		3.932
València/Valencia	2565124	2495425	2.72	2514083	1.99	2513028	2.03	2520861	1.726		2.116
Valladolid	519546	512426	1.37	517502	0.39	515630	0.75	515004	0.874		0.848
Bizkaia	1152651	1096502	4.87	1111404	3.58	1116881	3.10	1100122	4.557		4.028
Zamora	172539	184928	7.18	182166	5.58	182903	6.01	185800	7.686		6.613
Zaragoza	964693	938100	2.76	949771	1.55	943340	2.21	949557	1.569		2.021
Ceuta	84777	66610	21.43	69111	18.48	68237	19.51	68374	19.348		19.692
Melilla	86487	72581	16.08	75598	12.59	74353	14.03	73073	15.510		14.552
		Mean	3.47	Mean	2.53	Mean	2.83	Mean	2.870	Mean	2.925
		R²	0.9998	R²	0.9999	R²	0.9998	R²	0.9998

Evaluation of spatial and Population-based data projection approaches

As explained in the Data projection using geographical layers overlay subsection within the Methods section, data projection is a critical step for the integration of data-sets reported in different geographical layers. The projection is a linear interpolation using a defined overlay matrix. We have implemented two different approaches to reconstruct overlay matrices; one approach uses the ratios between areas (Spatial-based) whereas the other introduces a correction based on the distribution of the population (Population-based). To compare the accuracy of both approaches we first project the population reported by MITMA mobility areas into municipalities using both, Spatial and Population-based overlay matrices. Second, we compare the population values resulting from the projection into municipalities from those of the Spanish census of 2019[27]. Figure 5 shows the comparison between the population projected using both overlay matrices with respect to those of the Spanish census. In both cases, we find correlation coefficients close to one. Nonetheless, the population projected the Spatial-based overlays exhibit a higher dispersion. To quantify these differences we calculate the root mean square error (RMSE) between the projected values and those from the census. The RMSE found for the Spatial-based and Population-based overlays are 3200 and 2500, respectively. Thus, although Spatial-based overlays seem to be a good approximation, population-based overlays introduce less distortion and produce more reliable results.

Fig. 5

Comparison of data projection approaches between geographical layers. The figure represents the census population reported by municipalities with respect to those values estimated from MITMA mobility data after its projection from the MITMA mobility layer into the municipalities layer. Panel (a) and (b) show the result of the projection using Spatial and Population-based overlays, respectively.

Comparison of OD mobility data with other studies

To evaluate the quality of the mobility data from MITMA we compare our daily OD matrices to those from a different data-set. For this comparison, we use data from a study conducted by the National Institute of Statistics (INE, Instituto Nacional de Estadística). This study is based on the analysis of the position of more than 80% of mobile phones throughout Spain and it provides OD flows between areas and between mobility areas as long as they involve more than 15 people[28]. The data-set includes daily reported OD matrices from March to June 2020 and the geographical layer in which the population flows are reported is very similar to MITMA mobility layer. Nonetheless, to compare both data-set we first identify the common set of mobility zones and find 2680 mobility zones common to both data-set. Using the set of common mobility areas, we compare the number of trips for every origin-destination combination in the data. We also aggregate the mobility network at provinces and autonomous communities levels and compare the resulting aggregated OD matrices. Moreover, we compare both mobility data-sets for four different dates. Figure 6 shows the results of the different comparisons. These comparisons show high correlations in all cases. Nonetheless, while the comparison of the aggregated OD matrices at both, province and autonomous community levels, show a very high correlation (), the comparison at the mobility areas level shows more dispersion (). As expected, the discrepancies increase when lower fluxes are compared and decrease when higher ones are considered. Strikingly, the slope of the linear regression is in all the cases, indicating the MITMA data-sets measure more trips than those reported in the INE study. This can be due to the fact that in the INE data-set the definition of trip differs from that used in the MITMA study; specifically, in the INE study, the destination areas are defined as those where the phones are more frequently detected during the 10 am to 4 pm period of a given day, as long as the phone has been in that area for at least two hours[28]. Beside these systemic differences, the comparison shows that there is a good agreement between both studies. Nonetheless MITMA mobility data record, here presented, offers trips reported on an hourly basis with a more accurate trip definition and the study, starting in February 2020 is expected to continue reporting data.

Fig. 6

Comparison of mobility data-sets. The figure shows the comparison between the number of trips reported by the MITMA and the INE data-sets at different levels of aggregation. From left to right, each panel presents the comparison at mobility areas level, province level, and autonomous communities level. From (a) to (d) the plots represent the comparison between different dates indicated on the top of each panel. The Pearson correlation coefficient is indicated on the bottom right corner of each panel.

Usage Notes

In this section, we present simple examples to illustrate the potential uses of this data source, as well as how to access the different data records. In addition, we also provide some suggestions on how to integrate and analyse the different data records. Data records can be accessed in three different ways. We have implemented a REST-API to allow querying the whole COVID-19 Flow-Maps system (see REST-API in Code availability section); on top the REST-API we developed command-line download toolkit to download the data records in a more simple way (see REST-API Code availability section); in addition, all data records are periodically uploaded to our data repository (https://github.com/bsc-flowmaps). Finally, we have also developed a set of web-based interactive data dashboards for exploratory analysis of the data records (https://flowmaps.life.bsc.es/flowboard/).

COVID-19 Incidence

The COVID-19 data record includes cases reported at different levels of spatial resolution (Fig. 7). Further, exploratory analysis shows that data-sets on COVID-19 cases reported by different sources may have different reporting biases. For instance, the data-sets that account for all of Spain, i.e. ES.covid_cca and ES.covid_pro do not report cases on weekends and these cases are reported on Mondays. To avoid this kind of biases we suggest users to apply rolling means on the time series previous to and analysis. Alternatively, users can also use weekly accumulated incidence for any given date. Weekly accumulated incidence can also be used as a proxy to estimate the number of infected individuals in a given date; a useful metric for epidemiological modelling. In general, we have found that the criteria used for reporting can vary between the different sources included in this data record. Furthermore, we find that in order to evaluate the health state of a particular region in absolute numbers, weekly accumulated or averaged incidence are good estimators. However, to compare the incidence between different regions we suggest the use of cases per 100,000 inhabitants to normalise by populations size. In Fig. 7 we use maps to show cases reported in different regions; the maps on the left-side panel show daily incidence whereas those on the right side show cases per 100,000 inhabitants. To facilitate the analysis of the data the weekly accumulated incidence as well as cases per 100,000 inhabitants have been pre-calculated and included in all the provided COVID-19 data-sets. The COVID-19 data record of Flow-Maps also provides cases reported at higher spatial resolution. For instance, Fig. 7 shows two incidences metrics at the level of provinces for all of Spain and at the level of the BHA of Cataluña. The resulting plot shows both the net incidence and the incidence by 100,000 inhabitants. We find maps are powerful visualisation tools and for that reason, we also provide geographical layers in GeoJSON format (see Geographical Layers in Methods sections). All COVID-19 data-set are cross-referenced to the corresponding layer in which the data is reported which facilitate the use of maps for visualisation proposes. Herein, maps were generated using Plotly library (see Code availability for further details on the tools used). To facilitate the exploratory analysis of COVID-19 data-sets our COVID-19 Flow-Board web application provides an integrated dashboard for visualising incidence plots, including interactive maps and time series curves. This analysis tool also allows the users to explore data on cases for different rates and at different scales of spatial resolution and can be accessed through the following link: https://flowmaps.life.bsc.es/flowboard/board_incidence/.

Fig. 7

COVID-19 cases reported at different spatial resolutions. Panel (a) and (b show daily incidence and accumulated incidence in a week per 100.000 people reported at the level of province and Cataluña BHA, respectively.

Mobility data analysis

To evaluate population mobility in general terms we use as an index the percentage of people taking at least one daily trip. This mobility index is more accurate than the total number of trips because a single person can perform several trips in a given date. Figure 8a shows the result of a coarse-grained analysis of population mobility from 14th February of 2020 to 15th January of 2021 in Spain, using the mobility index. The results show a seasonal behaviour with a drop on the weekends. For this reason, when comparing mobility for different dates, we suggest comparing dates that correspond to the same day of the week. Moreover, the mobility index shows an abrupt change at the beginning of the lockdown that lasted from 14th March until the 21st Jun of 2020 (Fig. 8a). The analysis also shows that the number of people taking no trips doubled, while the number of people taking more than two trips was reduced to half (data not shown). Another change in the mobility patterns can be noted at the beginning of January when the Iberian Peninsula was hit by an extraordinary snowstorm, the largest in Spain since 1971[29]. As a consequence, population mobility was severely reduced due to the state of the roads and highways.

Fig. 8

Population mobility patterns in Spain during 2020. Panel (a) shows the population percentage that performs one or more trips in a given day. Panel (b) and (c) show the total number of daily trips between, and within autonomous communities, respectively. Grey shaded area corresponds to the strict lockdown, whereas dashed lines indicate the date of the mobility networks represented in panel (d). Panel (d) represents the mobility networks between MITMA zones for three selected dates that are annotated on the left bottom corner and indicated with a dashed vertical line in panels (a–c).

On the other hand, OD matrices should be used when analyzing population mobility flows between different regions. Figure 8b,c show the number of trips within and between autonomous communities. While the analysis of the flows also show a very similar pattern to the one exhibited by the mobility index (Fig. 8a), the analysis by flows also indicate change during summer. For instance, an increase in the number of trips between autonomous communities is observed in summer whereas the opposite trend is found for internal trips. Finally, to illustrate the granularity of the mobility data Fig. 8d shows a representation of the origin-destination network at the level of MITMA mobility areas for three different states: before the pandemic, during the lockdown, and after the summer. To facilitate the exploratory analysis of the mobility data-sets the COVID-19 Flow-Board provides an integrated dashboard for visualising map-based representation of the mobility networks at different scales of spatial resolution. This analysis tool basic data analytics on population mobility and can be accessed through the following link: https://flowmaps.life.bsc.es/flowboard/board_mobility/.

In order to integrate mobility and COVID-19 cases, we have introduced a risk score named Mobility Associated Risk. The MAR score combines mobility flows and incidence to estimate the potential risk of importing/exporting cases (see Mobility Associated Risk in Methods section). Thus, the MAR score could be of potential use to assist the design of targeted NPI, e.g. confining a specific region, or to inform citizens about the current situation. Moreover, the MAR score can be used to compare the different source of incoming risk for a given region. For instance, Fig. 9 shows the evolution of the MAR score incoming to Asturias during 2020. We have selected the top four sources of risk. To illustrate the idea underlying the MAR score, the Figure shows the evolution of cases in each source together with the number of trips from the source to Asturias (Fig. 9a,b). Combining the former indicators the MAR score is calculated (see Fig. 9c). The Figure allows seeing that during summer, the main source of incoming MAR for Asturias was Madrid. Interestingly, for those dates, the number of trips coming from Madrid is considerably lower than the other three sources. Nonetheless, the number of cases in Madrid was growing. Altogether, this example highlights the importance of combining mobility and cases to have more broad pictures of the epidemiological process. In the presented example, MAR was calculate by aggregating mobility and cases to the level of provinces. Nevertheless, the MAR score can be calculated at different levels of aggregation. Calculated MAR at different levels of spatial resolution can be downloaded from our systems and it is also included in daily reported individual file in the data repository. The COVID-19 Flow-Board also provides an integrated dashboard for visualising MAR scores between regions at different scales of spatial resolution. The MAR data dashboard can be acceded through the following link: https://flowmaps.life.bsc.es/flowboard-test/board_risk/.

Fig. 9

Mobility Associated Risk incoming to Asturias during 2020. The figure shows the time evolution of the COVID-19 cases, the trips and the MAR score for the top four main sources of incoming MAR related to Asturias. Panel (a) shows the evolution of the number of COVID-19 cases per total number of inhabitants. Panel (b) shows the evolution of the trips from the top four sources to Asturias. Panel (c) shows the evolution of MAR incoming to Asturias from the top four main sources of risk.

REST-API

We provide a REST API that gives access to the data stored in the MongoDB database. It allows consumers to retrieve multiple documents using query strings, allowing for filtering and sorting. Filters can be provided using MongoDB syntax or Python conditional expressions. It has been created using the open-source python-eve framework (version 1.1.2). Among other fetch-able collections, the API allows downloading the time-dependent population mobility networks across Spain, and daily reports of COVID-19 cases in Spain, at different levels of spatial resolution. The API can be accessed using the endpoint: https://flowmaps.life.bsc.es/api/. Complete documentation with examples can be found at: https://flowmaps.life.bsc.es/api/docs.

Command-line utility and python library for downloading the data

Although all the data records presented in this work is accessible through a REST API, to enable easy access to most relevant collections, we have implemented a lightweight library and a series of scripts in python language to ease the task of downloading the data-sets. The package allows fetching and ingesting the time-dependent population mobility networks across Spain (provided by MITMA), and daily reports of COVID-19 cases in Spain, at different levels of spatial resolution, provided by the National Centre of Epidemiology (CNE) and the different Autonomous Communities. In addition, the package allows downloading the multiple geographical layers for Spain, in GeoJSON format, at the same levels of spatial resolution as the mobility and COVID-19 data. All this code has been packed and published as the standalone python package flowmaps-data. The flowmaps-data toolkit can be installed from pypi repositories with the command pip install flowmaps-data. The package comes with a command line utility that can be used to list, describe and download the data-sets in several formats (CSV, JSON). A complete list of examples and installation instructions can be found here: https://pypi.org/project/flowmaps-data/. Altogether, this toolkit allows the user to query and retrieve different data records including, geographical layer, COVID-19 cases and population mobility, updated on daily basis for Spain.

Measurement(s)	COVID-19 cases • mobility • population
Technology Type(s)	Real Time PCR • mobile phone
Sample Characteristic - Organism	Homo sapiens • Severe acute respiratory syndrome-related coronavirus
Sample Characteristic - Location	Kingdom of Spain