Dewulf et al. Int J Health Geogr (2016) 15:14 DOI 10.1186/sl 2942-016-0042-z
International Journal of Health Geographies
RESEARCH
Open Access
Dynamic assessment of exposure to air pollution using mobile phone data
1
CrossMark
Bart Dewulf1,2,3*, Tijs Neutens1, Wouter Lefebvre3, Gerdy Seynaeve4, Charlotte Vanpoucke5, Carolien Beckx3 and Nico Van de Weghe1
Abstract
Background: Exposure to air pollution can have major health impacts, such as respiratory and cardiovascular diseases. Traditionally, only the air pollution concentration at the home location is taken into account in health impact assessments and epidemiological studies. Neglecting individual travel patterns can lead to a bias in air pollution exposure assessments.
Methods: In this work, we present a novel approach to calculate the daily exposure to air pollution using mobile phone data of approximately 5 million mobile phone users living in Belgium. At present, this data is collected and stored by telecom operators mainly for management of the mobile network. Yet it represents a major source of information in the study of human mobility. We calculate the exposure to N02 using two approaches: assuming people stay at home the entire day (traditional static approach), and incorporating individual travel patterns using their location inferred from their use of the mobile phone network (dynamic approach).
Results: The mean exposure to N02 increases with 1.27 ug/m3 (4.3 %) during the week and with 0.12 ug/m3 (0.4 %) during the weekend when incorporating individual travel patterns. During the week, mostly people living in municipalities surrounding larger cities experience the highest increase in N02 exposure when incorporating their travel patterns, probably because most of them work in these larger cities with higher N02 concentrations.
Conclusions: It is relevant for health impact assessments and epidemiological studies to incorporate individual travel patterns in estimating air pollution exposure. Mobile phone data is a promising data source to determine individual travel patterns, because of the advantages (e.g. low costs, large sample size, passive data collection) compared to travel surveys, GPS, and smartphone data (i.e. data captured by applications on smartphones).
Keywords: Air pollution, Dynamic assessment, Exposure, Mobile phone data,Travel patterns
Background
A large body of evidence indicates that exposure to air pollution causes various acute and chronic health effects, such as respiratory and cardiovascular diseases [1-9]. Approximately 2 million deaths worldwide are caused by air pollution annually [10]. Mainly black carbon (BC), particulate matter (PM), and nitrogen dioxide (N02) are identified as culprits of negative health effects.
Current health impact assessments and epidemiological studies examining exposure to air pollution often only
'Correspondence: bartd.dewu lf@ugent.be
1 Department of Geography, Ghent University, Krijgslaan 281, S8,
9000 Ghent, Belgium
Full list of author information is available at the end of the article
take the air pollution concentration at the home location into account [11-17]. Such static approach does not incorporate individual travel patterns and may lead to a bias in exposure and health assessments [18-23].
Detailed information on travel patterns is thus needed to obtain more dynamic estimates of the exposure to air pollution. Previous research showed an increase in exposure to air pollution by incorporating individual travel patterns [24], but the outcome depends on the air pollution concentration at the home location [25]. To assess individual travel patterns, often self-reported household travel surveys are used [26]. Major disadvantages of this approach are the large non-response rate [27], non-representative samples [28], and high costs [26]. Alternatively,
O© 2016 Dewulf et al. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License BioRflCCl CGntfSl      (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/ publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Dewulf ef a/. Int J Health Geogr (2016) 15:14
Page 2 of 14
mathematical models of travel patterns can be used [18, 24, 29]. This approach allows to draw more quantitative conclusions from a larger population size, but results are however only valid for situations similar to those for which their initial parameters were estimated. More recently, Global Positioning Systems (GPS) [25, 30] or smartphone data [31, 32] were used to provide detailed information on people's travel patterns. However, data collection with GPS or smartphone devices is often intensive for both researchers and participants, expensive and only a limited number of people can be tracked.
To overcome the limitations of travel surveys, travel models and GPS/smartphone data, mobile phone data can be used to derive information on individual travel patterns. At present, this kind of data is collected and stored by telecom operators mainly for management of the mobile network. However, it represents a major source of information in the study of human mobility. This data is continuously available, does not need additional costs to collect, and is often available for millions of phone users. With over 6 billion mobile subscriptions globally and a growing awareness of telecom operators of the potential, this data source offers a wide range of applications and research possibilities [33]. However, the number of studies published with such data is limited up to now because of privacy issues and problems accessing the data [34]. Previous studies using this type of data mostly analyse population densities [35-37], tourism [38, 39], and mobility [33-35, 40-44]. To our knowledge, no studies have used mobile phone data to dynamically estimate the exposure to air pollution.
This research will add knowledge to the existing strand of literature by calculating the exposure to air pollution using mobile phone data of more than 5 million people in Belgium. Our main objective is to bring evidence on how this innovative, underused data source can offer more dynamic estimations of the exposure to air pollution. Further, we explore how using daily averaged and hourly air pollution concentrations influences the results.
Methods
Data
Mobile phone data
Mobile phone data (or passive mobile positioning data) is based on signalling information that is exchanged between mobile devices and the mobile network. When using the mobile network, there is a flow of signalling information between the device and the mobile network. The mobile device switches to the antenna with the strongest radio coverage, which is typically the closest one. The signalling messages contain an indication of the antenna in use.
The data used for this study is available from probes installed in the Proximus network, which capture this information. We have data available from more than 4000 antennae sites. On each antenna site there are typically three or four antennae, delivering network coverage in diverged directions. As sites can be equipped with 2G, 3G and 4G technologies in different frequency bands, we make abstraction of the different technologies and group all cells that are co-located on the same antenna site and cover the same sector to considerably reduce complexity. This leads to more than 10,000 macro cells covering the entire country of Belgium. The mobile phone location is thus available at the precision of these macro cells, with each cell having its own, unique geographical coverage area and identity code. Figure 1 shows the antennae with the associated macro cells for the region of Ghent, overlaid on the road network. Because of the higher capacity needs, macro cells are smaller in urban areas and larger in rural areas. Figure 2 shows a histogram of the area of the macro cells. Some macro cells are larger than 10 km2 (with a maximum of 49 km2), but 50 % of them have an area smaller than 2 km2.
Data is collected from more than 5 million users of the Proximus network, which are representative for the Belgian population [45]. In Belgium, Proximus has a market share of about 41 %, which is higher than the other Belgian telecom operators: Mobistar (27 %), Telenet (14 %), Base (11 %), and others (7 %) [46].
The network probing system collects from all active users. Each data point consists of an anonymised user ID, the date and time of the transaction, the cell where the transaction occurred, and the transaction type. The following transaction types are possible: (a) turning on and off the phone; (b) setting up, maintaining and terminating calls; (c) sending and receiving text messages; (d) setting up, maintaining and terminating data sessions; (e) location update (when changing from location area; a location area is a group of cells of which there are approximately 65 in the Proximus network); (f) periodic location update (automatic update every 3 h when there is no activity). For this study, we used mobile phone data for both 1 week and 1 weekend day: Thursday October 8 and Saturday October 11, 2015. Because of regulation terms we had limited access to the data, but these 2 days were chosen to be as representative as possible, in terms of weather conditions for the time of the year, and travel behaviour (e.g. no holidays). No data of the home location was available due to privacy issues. Therefore, we used the location of the users at 4 am as a proxy for their home location (hereafter called reference location), since it is assumed most of the people are at home at that time.
Privacy issues of using mobile phone network data are a major concern of phone owners, telecom operators,
Dewulf efa/./nf J Health Geogr (2016) 15:14
Page 3 of 1
Dewulf ef a/. Int J Health Geogr (2016) 15:14
Page 4 of 14
researchers, and the general public. Because of this, no personal information is linked to the mobile phone data, and IDs that can link directly to individuals are removed. Individual exposure measures were aggregated to postal code level for mapping purposes.
Air pollution data
We focus on N02, an understudied pollutant that can cause an increase in pulmonary morbidity, a worsening of obstructive lung disease, and a higher susceptibility to airway infections [9, 47, 48]. Hourly N02 concentration data (in ug/m3) for Belgium was provided by the coupled RIO-IFDM model [49]. This model couples the land use regression model RIO, the road emissions model MIMOSA4 (taking into account COPERT4 emission functions, vehicle fleet and vehicle counts), and the Gaussian plume model IFDM. The latter is used to incorporate large concentration variations close to the major air pollution sources, such as roads and point sources. The model has been validated extensively for the discussed region [49, 50]. Hourly air quality measurements are provided by the Belgian Interregional Environment Agency [51].
In line with the mobile phone data, N02 concentration patterns in Belgium were modelled during 2 days. N02 concentration levels varied from 3 to 63 ug/m3 on the weekday (Thursday October 8, 2015) and from 5 to 54 ug/m3 on the weekend day (Saturday October 11, 2015). Figure 3 shows the mean N02 concentration for the entire country of Belgium, for both Thursday October 8 and Saturday October 11, 2015.
Data processing
Mobile phone data are collected and stored by the telecom operator, mainly for management of the mobile network and technical operations. Because each user in the mobile network has a different mobile activity, the temporal resolution of the data varies. The last known position (cell) of each user was used at a temporal resolution of 15 min. We assume that, when there is no new data point within 15 min, the user is still at the same location as before. As an example, Fig. 4 shows the user density (number of users per cell divided by the cell area) of Thursday October 8, 2015 at 12 am UTC.
Proximus has 5,574,000 active costumers [52]. Active costumers are costumers who have made or received at least one call, or sent or received at least one message in the last 3 months, or if at least one data connection has been made on the last month. From the initial dataset, users are omitted if:
- they are international users,
- their data relate to machine to machine transactions (e.g. car kits; to avoid duplicate data),
- their travel patterns exceed the borders of Belgium during the selected days (no air pollution concentration available),
- they have no known position from 1:00 until 4:00 in the morning (necessary to derive the reference location).
This results in a dataset of 3,465,917 users on the weekday and 3,495,453 on the weekend day.
The point dataset of the N02 concentration is triangulated to a 50 x 50 m grid, using the SAGA 'gridding triangulation tool' in QGIS. Following, we calculate the average N02 concentration per cell using the SAGA 'grid statistics for polygons tool' in QGIS to combine this with the location data. The mean N02 concentration per cell is 29.36 ug/m3 on the weekday and 27.32 ug/m3 on the weekend day, with a mean standard deviation per cell of respectively 3.62 (ig/m3 and 2.73 (ig/m3.
The location of the users is combined with the air pollution concentration, to calculate the exposure to air pollution. The air pollution data is in Coordinated Universal Time (UTC), and the mobile phone data is in local time (UTC+1). Therefore, we have an overlap of 23 h (92 quarters) per day, and are thus able to combine the datasets from 0 am UTC to 11 pm UTC. The exposure to N02 is calculated using either a static or dynamic location. For the static approach, we use the cell where the user is at 4 am UTC as their reference location. For the dynamic approach, we use the exact cell where the user is, at a temporal resolution of 15 min. Additionally, we use the N02 concentration per cell in two different ways. We either use the hourly concentration or the daily averaged concentration per cell. This results in four possible average air pollution concentrations each user is exposed to during the day.
Data analyses
First, having these four average air pollution exposure values for 2 days, we check what the influence of using hourly air pollution concentrations (hour) is on the exposure to air pollution, compared to using the daily average air pollution concentration (day). Second, we compare the effect of using the reference location of the user (static) with taking into account the actual location of the user (dynamic) on the calculated exposure to air pollution. Third, this comparison is also analysed geographically.
Dewulf ef a/. Int J Health Geogr (2016) 15:14
Page 5 of 14
Mean N02 concentration (ug/m3) I     I 0.13 I     I 5.70 I     I 11.27
□ 16-83 N I I 22.40 I I 27.97 I I 33.54 I I 38.68 I     I 42.96
week
Mean N02 concentration (ug/m3) I     I 0.13 I     I 5.70 I     I 11.27
□ 16-83 ^ I | 22.40 I | 27.97 I I 33.54 I I 38.68 I     I 42.96
weekend
Fig. 3 Map showing the mean N02 concentration for the entire country of Belgium, for both Thursday October 8 and Saturday October 11 201 5
Dewulf ef a/. Int J Health Geogr (2016) 15:14
Page 6 of 14
The data was statistically analysed using R 3.2.2 . To check for significant differences between the approaches, paired-samples t tests (hour vs. daily averaged air pollution concentration, static vs. dynamic approach, week vs. weekend days) were performed. The data does not have to be tested for normality, because of the large sample size [53]. Statistical significance was set at p < 0.05. Geographical analyses were performed in QGIS 2.12™. Averages of individual exposure values were calculated per municipality and visualised using choropleth maps.
Results
To gain insight into the origin of the four average values per user, Fig. 5 shows the exposure to N02 during the weekday for a random user, calculated statically and dynamically, with both hourly and daily averaged N02 concentrations. Using hourly N02 concentrations leads to a higher level of detail of the exposure to air pollution. However, since we calculate the average exposure to air pollution per day, this has limited effects on the results. It is also clear that taking into account actual travel patterns
(instead of assuming the person stays at the reference location) leads to a different exposure to air pollution. In this case, the person spends time in cells with a higher N02 concentration than at his or her reference location.
The mean N02 exposure per person was calculated statically and dynamically, using hourly N02 concentration values for both the week and weekend day, and is presented in Table 1. The significance of the differences was tested using multiple paired-samples t tests (daily averaged vs. hourly N02 concentrations, static vs. dynamic, week vs. weekend).
Comparison of using hourly or daily averaged air pollution concentrations
From Table 1, we observe practically no difference in the mean N02 exposure calculated with hourly and daily averaged N02 concentrations for the static approach, which is expected. For the dynamic approach, we observe a small significant difference (p < 0.001) only during the week. Here, the calculated N02 exposure is 0.13 (J.g/m3 (0.4 %) higher when using hourly values compared to when using daily averaged values.
Dewulf ef a/. Int J Health Geogr (2016) 15:14
Page 7 of 14
N02 concentration (ug/m3)
45 -
40
35
30
25
20
10
Dynamic_day
Static_day
Static_hour Dynamic_hour
our_ /
Time (h)
0    1    2    3    4   5    6    7    8    9    10  11  12 13   14   15 16   17 18   19 20   21 22 23
Fig. 5 Exposure to N02 during the weekday for a random user, using the four different approaches (static_hour, dynamic_hour, static_day, dynamic_day)
Table 1 Mean exposure to N02 per person, calculated statically and dynamically, using hourly and daily averaged N02 concentrations, for both the weekday (n = 3,465,917) and weekend day (n = 3,495,453)
Method	Mean N02 exposure (ug/m3i		
	Weekday	Weekend day	
	N02 per hour N02 per day	N02 per hour	N02 per day
Static Dynamic Difference	29.69 [12.03]     29.69 [12.03] 30.96 [11.26]     30.83 [11.25] 1.27 [5.02]        1.14 [4.43]	27.47 [8.58] 27.59 [7.99] 0.12 [2.82]	27.47 [8.58] 27.57 [8.01] 0.10 [2.41]
Comparison of the static and dynamic calculation of the exposure to air pollution for a week and weekend day
To compare the static with the dynamic approach, we will only consider the values calculated with the hourly N02 concentrations, since this way the highest level of detail is obtained.
Table 1 shows that by incorporating individual travel patterns (dynamic), the mean exposure to N02 increases with 1.27 ug/m3 (4.3 %) on the weekday and with 0.12 ug/m3 (0.4 %) on the weekend day, compared to assuming the person stays at the reference location (static). These values are all significantly different from each other (p < 0.001). Figure 6 combines the static and dynamic approach for the two days in a histogram. It is
clear that during the week, there are more users who experience an increase in exposure to N02. During the weekend, the values are more central and the increase in N02 exposure is less pronounced. During the week, 12.4 % of the users have no change, 54.5 % have an increase, and 33.1 % have a decrease. During the weekend, 20.1 % have no change, 43.3 % have an increase, and 36.6 % have a decrease.
Figure 7 shows a scatterplot of the exposure to N02 calculated statically and dynamically, for both the week and weekend day. During the week, it is clear that users with a low N02 exposure calculated statically (thus with a low average N02 concentration at the reference location) experience a strong increase in N02 exposure when dynamically calculated (indicated in blue), a pattern that is less pronounced during the weekend. This is also true the other way round: users with high N02 concentrations at the reference location experience a decrease in N02 exposure when their travel patterns are considered (indicated in green), which is also observed during the weekend.
Geographically analysing the comparison between the static and dynamic calculation of the exposure to air pollution
Next to these statistical analyses, we also performed geographical analyses on the comparison between the static and dynamic approach to calculate the N02 exposure.
Dewulf ef a/. Int J Health Geogr (2016) 15:14
Page 8 of 14
week
dynamic
St
20
40
N02 concentration (ug/m3)
60
[
80
weekend
dynamic
i-1--1-1-
0 20 40 60
N02 concentration (ug/m3)
Fig. 6 Histogram of the average N02 concentration that the users are exposed to using the static and dynamic approach, for the week (n = 3,465,917) and weekend day (n = 3,495,453), including the mean reference line for both approaches
i
80
Here, again only hourly N02 concentrations were used in the calculations.
Figure 8 shows the average exposure to N02 per municipality, for both the week and weekend day, calculated statically and dynamically with hourly N02 concentrations. Figure 9 shows the difference between the average exposure to N02 calculated with static and dynamic approach (dynamic minus static), using hourly N02 concentrations, per municipality for both the week and weekend day. During the week, there is a large increase mostly in the municipalities surrounding larger cities (Brussels, Antwerp, Ghent) and a decrease in these
larger cities. During the weekend we observe a similar pattern, but with lower increases and more decreases in the difference between the static and dynamic approach.
Discussion
General discussion
First, our study shows that using daily averaged instead of hourly N02 concentrations leads to only a 0.4 % decrease in dynamically calculated exposure to N02 in our analyses. If hourly data is available, it is preferred to use it [18, 31]. If, however, no detailed hourly N02 concentration data is available, the impact will be limited.
Dewulf ef a/. Int J Health Geogr (2016) 15:14
Page 9 of 14
Second, our study supports the findings of several exposure [18, 20-23, 25], an issue sometimes overlooked recent studies stating the importance of incorporat- [11]. Mobile phone data makes it possible to esti-ing individual travel patterns in estimating air pollution    mate individual travel patterns to use in health impact
Dewulf ef a/. Int J Health Geogr (2016) 15:14
Page 10 of 14
Week - static
Exposure to N02 (ug/m3) I    I 0.0 -10.0 I    I 10.1 -20.0
■ 20.1 -30.0
■ 30.1 -40.0 I    I >40.0
I   I No value
Weekend - static
Exposure to N02 (ug/m3) I    I 0.0 -10.0 I    I 10.1 -20.0
■ 20.1 -30.0
■ 30.1 -40.0 I    I >40.0
I   I No value
Weekend - dynamic
Fig. 8 Maps of Belgium, showing the statically and dynamically calculated exposure to N02, for the week (n = 3,465,917) and weekend day (n = 3,495,453)
assessments and epidemiological studies. We observe a mean increase in N02 exposure of 4.3 % during the weekday and 0.4 % during the weekend day when incorporating individual travel patterns, which means that current health impact assessments underestimate the exposure to N02 and the related acute and chronic health effects. These increases were also found in previous research, where integrating time-activity information lead to a 1.2 % increase air pollution exposure than when assuming people are always at their home location [24]. We observed an increase or decrease in N02 exposure for respectively 54.5 and 33.1 % of the users because of their travel patterns during the week. In the weekend, respectively 43.3 and 36.6 % of the users experience an increase or decrease in exposure. Thus, people tend to make more
trips to areas that are less polluted than their reference location in the weekend than during the week. During the week, people living in areas with a low N02 concentration undergo an increase in N02 exposure because of their travel patterns (going to work in a more polluted area) whereas people living in highly polluted areas undergo a decrease in N02 exposure, which is similar to our previous study [25].
Third, concerning the geographical analysis, our study reports that people living near Brussels are most exposed to N02 both during the week and the weekend, because of the highest density of air pollution sources (industry and roads), and Brussels being one of the most congested cities in Europe [54, 55]. People living in the south of Belgium are least exposed to N02. Mainly people from
Dewulf efa/./nf J Health Geogr (2016) 15:14
Page 11 of 14
Dewulf ef a/. Int J Health Geogr (2016) 15:14
Page 12 of 14
municipalities around larger cities experience an increase in exposure to air pollution during the week because of their travel patterns (going to work in these cities). The average difference for people living in these cities is negative since people working in the city do not experience any change and people working outside the city experience a decrease in N02 exposure. During the weekend, we observe lower increases and more decreases in N02 exposure because of individual travel patterns, because during the weekend people tend to visit more areas with lower air pollution concentrations than at the reference location. In more rural areas (e.g. the south of Belgium) there is an increase in N02 exposure during the weekend when incorporating travel patterns, because every trip people living in this area make leads to an increased N02 exposure due to the low air pollution concentration at people's reference location.
Strengths and limitations
This study has several strengths compared to similar studies. First, to our knowledge, this is the first study that combines mobile phone data with air pollution concentration data to dynamically estimate the exposure to air pollution. Using mobile phone data has several advantages above GPS data or questionnaires: no additional costs have to be made to collect the data, a very large number of people can be traced because of the wide adoption of the mobile phone, the data collection is passive so people are not disturbed and does not influence the battery of their mobile devices, and they are tracked without them knowing so they don't change their normal behaviour. Numerous practical applications can be developed based on the presented method, both for individual as on a community level. Policy makers can for example be interested to follow-up the average population exposure indicator or to assess the impact of a policy measure such as on the exposure. However, it is not easy to access mobile phone data. Concerning privacy issues, good agreement with the telecom operator is needed, as well as a clear understanding of the data use [33, 38]. A second strength of the current study involves the type of location data that were used. Previous studies using mobile phone data only collected a location when users made a phone call or sent a text message. In our study, we additionally locate users when turning the phone on or off, during a data session, when changing location area, or when periodically updating the location area by the telecom operator. This approach significantly improves the spatial accuracy of all the users but also includes relevant information on non-frequent callers in the population. Third, using modelled air pollution concentrations instead of personal measurements offers nation-wide data on a detailed geographical scale [18, 24] and is easier
to generate than personal measurements [19]. Fourth, we used both the hourly and daily averaged N02 concentrations in contrast to our previous research [25], making a comparative analysis possible. Using daily averaged values does not alter the results extremely, since we calculated daily exposure values. It is however preferred to use the hourly values when available, to obtain more accurate results.
Apart from these strengths, this study also has some limitations that open up interesting avenues for future work. First, we only used air pollution concentration and individual travel data for two days. Despite the fact that these days were chosen to be as representative as possible, it would be better to use more data, e.g. for an entire week, month or even year. Also, in order to assess the associated health impacts, more data is required. Second, following the privacy issues of mobile phone data, it is difficult to combine this data with personal sociodemographic variables or other semantic information (e.g. transport purpose and trip mode), which limits the analysis possibilities. Future research could try to deduce the sociodemographic characteristics from the most likely living place (possible to determine using long-term location data) to get an idea of the socioeconomic status of the users. These privacy issues could be addressed by using privacy-enhancing technologies [56]. One possible solution is to slightly obfuscate the location of the user, while keeping enough information to perform satisfying analyses [57]. An other possibility is that telecom providers could ask for an opt-in consent from their costumers to make use of their location data for scientific research and try to build a trust relationship with them, and build services where costumers benefit from. Third, since we had no information on the user's home location, we used the location at 4 am as a proxy for their home location (and use it as reference location). A better solution would be to use mobile phone data from a longer period (e.g. 1 month), to make a more accurate estimate of the most likely living location. Fourth, the spatial resolution of mobile phone data is limited to that of the used cells, which is low compared to the spatial accuracy of GPS data. Because of this low spatial resolution, local differences in air pollution concentration (e.g. near roads) may not be taken into account. This might mean that the exposure to air pollution is probably even higher because of people spending time in traffic [21]. On the other hand, the pollution concentration gradients in the rural areas, where the macro cell size is larger, are in general very small. As a result, smaller rural macro cells, if available, would not increase the accuracy of our results. The spatial resolution of the data could be increased by applying triangulation [36, 40].
Dewulf ef a/. Int J Health Geogr (2016) 15:14
Page 13 of 14
Conclusions
Hourly air pollution concentrations are preferably used over daily averages to maximise the level of detail when combining air pollution with individual travel patterns. This study shows that for epidemiological studies and exposure assessments, it is relevant to incorporate individual travel patterns to estimate the exposure to air pollution. The change in exposure to air pollution depends on the air pollution concentration at the reference location and someone's individual travel patterns, but on average we found an increase of 4.3 % in the exposure to N02 during the week and 0.4 % during the weekend. People living in and near large cities are most exposed to N02. However, people from other areas experience a higher increase in N02 exposure when taking their travel patterns into account. Mainly people living in municipalities surrounding larger cities have an increase in N02 exposure because they work in these cities. Aside from privacy issues, we strongly believe that using mobile phone data has several advantages (e.g. low costs, large sample, passive data collection) over travel surveys, GPS, and smartphone data. Especially for air pollution research the applications of using mobile phone data are numerous. Policy makers can use this information to assess the impact of air pollution on the population. Also, they can analyse the impact of a certain policy measure or occurring events (e.g. festivals, strikes) on the individual travel patterns and assess the associated impacts on exposure to air pollution. Mobile phone data is therefore a promising data source for air pollution research.
Authors' contributions
BD carried out the main research and drafted the manuscript. TN and CB participated in the study design and helped to draft the manuscript. NVdW and CB coordinated the study. WL and CV provided the air quality modelling. All authors read and approved the final manuscript.
Author details
1 Department of Geography, Ghent University, Krijgslaan 281, S8, 9000 Ghent, Belgium.2 Research Foundation Flanders, Egmontstraat 5,1000 Brussels, Belgium.3 VITO, Boeretang 200,2400 Mol, Belgium.4 Proximus, Koning Albert ll-laan 27,1030 Brussels, Belgium.5IRCEL, Kunstlaan 10-11,1210 Brussels, Belgium.
Acknowledgements
This research is supported by Proximus.
Competing interests
The authors declare that they have no competing interests. Availability of data and materials
Due to privacy issues, the data used for this research are not available. Ethics approval
All data was anonymously treated to protect the privacy of the costumers. Funding
This research is funded by the Research Foundation Flanders (FWO) B/l3230/01.
Received: 12 January 2016 Accepted:21 March 2016 Published online: 21 April 2016
References
1. Beelen R, Stafoggia M, Raaschou-Nielsen O, Andersen ZJ, Xun WW, Katsouyanni K, Dimakopoulou K, Brunekreef B, Weinmayr G, Hoffmann B, Wolf K, Samoli E, Houthuijs D, Nieuwenhuijsen M, Oudin A, Forsberg B.OIsson D, SalomaaV, LankiTYli-TuomiT Oftedal B, AamodtG, Nafstad P De Faire U, Pedersen NL, Ostenson C-G, Fratiglioni L, Penell J, Korek M, Pyko A, et al. Long-term exposure to air pollution and cardiovascular mortality: an analysis of 22 European cohorts. Epidemiology. 2014;25:368-78.
2. Brook RD, Rajagopalan S, Pope CA, Brook JR, Bhatnagar A, Diez-Roux AV, Holguin F, Hong Y, Luepker RV, Mittleman MA, Peters A, Siscovick D, Smith SC, Whitsel L, Kaufman JD. Particulate matter air pollution and cardiovascular disease: an update to the scientific statement from the American Heart Association. Circulation. 2010;121:2331 -78.
3. Brugge D, Durant JL, Rioux C. Near-highway pollutants in motor vehicle exhaust: a review of epidemiologic evidence of cardiac and pulmonary health risks. Environ Heal. 2007;6:12.
4. Gehring U, Gruzieva O, Agius RM, Beelen R, Custovic A, Cyrys J, Eeftens M, FlexederC, Fuertes E, Heinrich J, Hoffmann B, DeJongsteJC, Kerkhof M. Air pollution exposure and lung function in children: the ESCAPE project. Environ Health Perspect. 2013;121:1357-64.
5. HEI. Traffic-related air pollution: a critical review of the literature on emissions, exposure, and health effects. 2010. http://pubs.healtheffects.org/ getfile.php?u=553. Accessed 13 Apr 2015.
6. Peters A, von Klot S, Heier M,Trentinaglia I, Hdrmann A, Wichmann HE, Ldwel H. Exposure to traffic and the onset of myocardial infarction. N Eng J Med. 2004;351:1 721-30.
7. Pope CA III, Dockery DW. Health effects of fine particulate air pollution: lines that connect. Air Waste Manag Assoc. 2006;56:709-42.
8. Riediker M, Cascio WE, Griggs TR, Herbst MC, Bromberg PA, Neas L, Williams RW, Devlin RB. Particulate matter exposure in cars is associated with cardiovascular effects in healthy young men. Am J Respir Crit Care Med. 2004;169:934-40.
9. WHO. Health aspects of air pollution with particulate matter, ozone
and nitrogen dioxide. 2003. http://www.euro.who.int/_data/assets/
pdf_file/0005/l 12199/E79097.pdf. Accessed 23 Mar 2015.
10. WHO. Review of evidence on health aspects of air pollution—REVI-
HAAP Project. 2013. http://www.euro.who.int/_data/assets/pdf_
file/0004/193108/REVIHAAP-Final-technical-report-final-version.pdf. Accessed 23 Mar 2015.
11. Brunekreef B, Hoek G, Schouten L, Bausch-goldbohm S, Fischer P, Armstrong B, Hughes E, Jerrett M, Brandt PVan Den. Effects of long-term exposure on respiratory and cardiovascular mortality in the Netherlands: the NLCS-AIR study. 2009. http://www.n65.nl/NCLS-AIR-Study-2009.pdf. Accessed 5 June 2015.
12. Jerrett M, Burnett RT, Beckerman BS, Turner MC, Krewski D, Thurston G, Martin RV, van Donkelaar A, Hughes E, Shi Y, Gapstur SM.Thun MJ, Pope CA. Spatial analysis of air pollution and mortality in California. Am J Respir Crit Care Med. 2013; 188:593-9.
13. Tenailleau QM, Mauny F, Joly D, Francois S, Bernard N. Air pollution in moderately polluted urban areas: how does the definition of "neighborhood" impact exposure assessment? Environ Pollut. 2015;206:437-48.
14. Bell ML, Ebisu K, Belanger K. Ambient air pollution and low birth weight in Connecticut and Massachusetts. Child Heal. 2007;115:1118-24.
15. Cesaroni G, Badaloni C, Porta D, Forastiere F, Perucci CA. Comparison between various indices of exposure to traffic-related air pollution and their impact on respiratory health in adults. Occup Environ Med. 2008;65:683-90.
16. Hoek G, Brunekreef B, Goldbohm S, Fischer P, Van Den Brandt PA. Association between mortality and indicators of traffic-related air pollution in the Netherlands: a cohort study. Lancet. 2002;360:1203-9.
17. Huynh M, Woodruff TJ, Parker J D, Schoendorf KC. Relationships between air pollution and preterm birth in California. Pediatr Perinat Epidemiol. 2006;20:454-61.
18. BeckxC, Int Panis L, ArentzeT, Janssens D.Torfs R, Broekx S, Wets G. A dynamic activity-based population modelling approach to evaluate
Dewulf et al. Int J Health Geogr (2016) 15:14
Page 14 of 14
exposure to air pollution: methods and application to a Dutch urban area. Environ Impact Assess Rev. 2009;29:179-85.
19. Dons E, Van Poppel M, Kochan B, Wets G, Int Panis L. Implementation and validation of a modeling framework to assess personal exposure to black carbon. Environ Int. 2014;62:64-71.
20. Setton E.The impact of daily mobility on exposure to traffic-related air pollution and health effect estimates. J Expo Sci Environ Epidemiol. 2011;21:42-8.
21. Dons E, Int Panis L.Van Poppel M.Theunis J, Willems H.Torfs R.WetsG. mpact of time-activity patterns on personal exposure to black carbon. Atmos Environ. 2011;45:3594-602.
22. Steinle S, Reis S, Sabel CE. Quantifying human exposure to air pollution— moving from static monitoring to spatio-temporally resolved persona exposure assessment. Sci Total Environ. 2013;443:184-93.
23. Valero N, Aguilera I, Llop S, Esplugues A, de Nazelle A, Ballester F, Sunyer J. Concentrations and determinants of outdoor, indoor and personal nitrogen dioxide in pregnant women from two Spanish birth cohorts. Environ Int. 2009;35:1196-201.
24. DhondtS, BeckxC, Degraeuwe B, LefebvreW, Kochan B, BellemansT, Int Panis L, Macharis C, Putman K. Integration of population mobility in the evaluation of air guality measures on local and regional scales. Atmos Environ. 2012;59:67-74.
25. Dewulf B, NeutensT, Van Dyck D, De Bourdeaudhuij I, Int Panis L, Beckx C, Van de Weghe N. Dynamic assessment of inhaled air pollution using GPS and accelerometer data. JTransp Heal. 2016. doi:l0.1016/j.jth.2015.10.004.
26. Stopher PR, Greaves SP. Household travel surveys: where are we going? Transp Res Part A Policy Pract. 2007;41:367-81.
27. Wilson J. Measuring personal travel and goods movement: a review of the bureau of transportation statistics'surveys. 2004. http://onlinepubs. trb.org/onlinepubs/sr/sr277.pdf. Accessed 23 Nov 2015.
28. Murakami E: Hard to reach populations (Presentation at the NYMTC survey workshop). 2008. https://www.nymtc.org/project/surveys/WORK-SHOP/Murakami.pdf. Accessed 18 June 2015.
29. Fecht D, Beale L, Briggs D. A GIS-based urban simulation model for environmental health analysis. Environ Model Softw. 2014;58:1-11.
30. Houston D, Ong P Jaimes G, Winer A. Traffic exposure near the Los Angeles-Long Beach port complex: using GPS-enhanced tracking to assess the implications of unreported travel and locations. J Transp Geogr. 2011;19:1399-409.
31. de Nazelle A, Seto E, Donaire-Gonzalez D, Mendez M, Matamala J, Nieu-wenhuijsen MJ, Jerrett M. Improving estimates of air pollution exposure through ubiguitous sensing technologies. Environ Pollut. 2013; 176:92-9.
32. Su JG, Jerrett M, Meng Y-Y, Pickett M, Ritz B. Integrating smart-phone based momentary location tracking with fixed site air guality monitoring for personal exposure assessment. Sci Total Environ. 2015;506-507:518-26.
33. Calabrese F, Ferrari L, Blondel VD. Urban sensing using mobile phone network data: a survey of research. ACM Comput Surv. 2014;47:1-20.
34. Ahas R, Silm S, Jarv O, Saluveer E,Tiru M. Using mobile positioning data to model locations meaningful to users of mobile phones. J Urban Technol. 2010;17:3-27.
35. Jonge E De, Pelt M Van, Roos M.Time patterns, geospatial clustering nd mobility statistics based on mobile phone network data. Stat Netherlands 2012:1-26. http://www.cbs.nl/NR/rdonlyres/4EDB51 ED-927A-4A69-38F3-4DC57A44DDE4/0/Timepatternsgeospatialclusteringandmobilitys-tatistics.pdf. Accessed 29 Oct 2015.
36. Ratti C, Frenchman D, Pulselli RM, Williams S. Mobile landscapes: using location data from cell phones for urban analysis. Environ Plan B Plan Des. 2006;33:727-48.
37. Deville P Linard C, Martin S, Gilbert M, Stevens FR, Gaughan AE, Blondel VD.Tatem AJ. Dynamic population mapping using mobile phone data. Proc Natl Acad Sci. 2014; 111:15888-93.
38. Ahas R, Aasa A, Roose A, Mark Ü, Silm S. Evaluating passive mobile positioning data for tourism surveys: an Estonian case study. Tour Manag. 2008;29:469-86.
39. Kuusik A, Nilbe K, MehineT, Ahas R. Country as a free sample: the ability of tourism events to generate repeat visits. Case study with mobile positioning data in estonia. Procedia Soc Behav Sci. 2014;148:262-70.
40. Alexander L, Jiang S, Murga M, Gönz MC. Validation of origin-destination trips by purpose and time of day inferred from mobile phone data. Transp Res Part C. 2015;58:1-20.
41. Widhalm P, Yang Y Ulm M, Athavale S, Gönz MC. Discovering urban activity patterns in cell phone data.Transportation (Amst). 2015;42:597-623.
42. Calabrese F, Di LG, Liang L, Carlo R. Estimating origin-destination flows using mobile phone location data. Pervasive Comput. 2011;10:36-44.
43. Pappalardo L, Simini F, Rinzivillo S, Pedreschi D, Giannotti F. Returners and explorers dichotomy in human mobility. Nat Commun. 2015;6:8.
44. Chen C, Bian L, Ma J. From traces to trajectories: how well can we guess activity locations from mobile phone traces?Transp Res Part C Emerg Technol. 2014;46:326-37.
45. Proximus. Proximus telecom universe survey. 2015.
46. Smart Business Strategies. Deel 6: telecom and netwerken. ICTjaarboek 2014-2015.2014;141:186-189.
47. Blomberg A, Krishna MT, Bocchino V, Biscione GL, Shute JK, Kelly FJ, Frew AJ, Holgate ST, Sandström T. The inflammatory effects of 2 ppm N02 on the airways of healthy subjects. Am J Respir Crit Care Med. 1997;156:418-24.
48. Cesaroni G, Stafoggia M, Galassi C, Hilding A, Hoffmann B, Houthuijs D, Ostenson C, Overvad K, Pedersen NL. Long term exposure to ambient air pollution and incidence of acute coronary events: prospective cohort study and meta-analysis in 11 European cohorts from the ESCAPE project. BMJ.2014;348:16.
49. LefebvreW, Degrawe B, BeckxC.Vanhulsel M, Kochan B, BellemansT, Janssens D, Wets G. Presentation and evaluation of an integrated model chain to respond to traffic- and health-related policy guestions. Environ Model Softw. 201 3;40:160-70.
50. VMM, VITO. Validation of the IFDM-model for use in urban applications. 2013. http://www.atmosys.eu/faces/doc/ATMOSYS_Deliverable_10_ FDM_Model_Validation.pdf. Accessed 15 Sept 2015.
51. IRCEL. http://www.irceline.be.
52. Leroy D, DufourS, Vandervoort P Van Den Meersche B, Standaert G, Tilmans R, Georgis M, Lybaert D. 2015 Q3 guarterly report. 2015. http:// cdn.proximus.com/proximusprod.artwhere.net/sites/default/files/Docu-ments/lnvestors/Reports/2015/Q3/Quarterly%20Report%20-%20Q3%20 2015%20results_final.pdf.
53. Altman DG, Bland JM.The normal distribution. BMJ. 1995;310:298.
54. OECD. OECD economic surveys: Belgium, vol. 2013. Paris: OECD Publishing; 2013.
55. FOD Mobiliteit en Vervoer. Kilometers Afgelegd Door Belgische Voertuigen in Het Jaar 2011. 2013. https://www.google.be/url?sa= t&rct=j&g=&esrc=s&source=web&cd=3&ved=0ahUKEwilvoCfm dzLAhXhC5oKHTQiCGYQFggpMAI&url=http%3A%2F%2Fmobilit. belgium.be%2Fsites%2Fdefault%2Ffiles%2Fdownloads%2Fpl3060 3an_studie%2520voertuigkilometersx.doc&usg=AFQjCNET4cwbwSwPCt noKLDMCUssS-dtpw&sig2=Sg3Wn42nuzQU14eJvht44w. Accessed 2 Oct 2015.
56. Giannotti F, Pedreschi D. Mobility, data mining and privacy. Berlin: Springer; 2008.
57. Wightman P Coronell W, Jabba D, JimenoM, Labrador M. Evaluation of location obfuscation technigues for privacy in location based information systems. In 2011 IEEE Latin-American conference on communications, LATINCOM 2011.2011.