Next Article in Journal
Partition Coefficients and Diffusion Lengths of 222Rn in Some Polymers at Different Temperatures
Previous Article in Journal
Intersectional Discrimination Is Associated with Housing Instability among Trans Women Living in the San Francisco Bay Area
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJERPH
Volume 16
Issue 22
10.3390/ijerph16224522
ijerph-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Estimation of Individual Exposure Levels to Air Pollution Using Trajectories Reconstructed from Mobile Phone Data

by
Mingxiao Li
1,2,3,
Song Gao
3,
Feng Lu
1,4,5,
Huan Tong
6 and
Hengcai Zhang
1,4,*
1
State Key Laboratory of Resources and Environmental Information System, Institute of Geographic Sciences and Natural Resources Research, Chinese Academy of Sciences, Beijing 100101, China
2
University of the Chinese Academy of Sciences, Beijing 100049, China
3
Geospatial Data Science Lab, Department of Geography, University of Wisconsin-Madison, Madison, WI 53706, USA
4
The Academy of Digital China, Fuzhou University, Fuzhou 350002, China
5
Jiangsu Center for Collaborative Innovation in Geographical Information Resource Development and Application, Nanjing 210023, China
6
UCL Institute for Environmental Design and Engineering, University College London, London WC1E 6BT, UK
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2019, 16(22), 4522; https://doi.org/10.3390/ijerph16224522
Submission received: 17 October 2019 / Revised: 11 November 2019 / Accepted: 13 November 2019 / Published: 15 November 2019
Browse Figures
Versions Notes

Abstract

:
The spatiotemporal variability in air pollutant concentrations raises challenges in linking air pollution exposure to individual health outcomes. Thus, understanding the spatiotemporal patterns of human mobility plays an important role in air pollution epidemiology and health studies. With the advantages of massive users, wide spatial coverage and passive acquisition capability, mobile phone data have become an emerging data source for compiling exposure estimates. However, compared with air pollution monitoring data, the temporal granularity of mobile phone data is not high enough, which limits the performance of individual exposure estimation. To mitigate this problem, we present a novel method of estimating dynamic individual air pollution exposure levels using trajectories reconstructed from mobile phone data. Using the city of Shanghai as a case study, we compared three different types of exposure estimates using (1) reconstructed mobile phone trajectories, (2) recorded mobile phone trajectories, and (3) residential locations. The results demonstrate the necessity of trajectory reconstruction in exposure and health risk assessment. Additionally, we measure the potential health effects of air pollution from both individual and geographical perspectives. This helped reveal the temporal variations in individual exposures and the spatial distribution of residential areas with high exposure levels. The proposed method allows us to perform large-area and long-term exposure estimations for a large number of residents at a high spatiotemporal resolution, which helps support policy-driven environmental actions and reduce potential health risks.

1. Introduction

With the acceleration of urbanization and industrialization in the past few years, ubiquitous and unavoidable air pollution has become a widespread health problem in many developing countries [1,2]. Air pollution has been proven to be associated with higher health risk of respiratory infection, chronic obstructive pulmonary disease, stroke, heart disease, and lung cancer, among others [3,4,5]. According to the European Society of Cardiology, 8.8 million extra deaths can be linked to air pollution worldwide each year [6]. To solve this problem, a better understanding of the negative health effects of air pollution is a prerequisite, which leads to an increased demand of accurate individual air pollution exposure estimates in public health studies [7,8,9].
In the field of human exposure estimation, the importance of human mobility has long been recognized [7,10,11,12]. To address this problem, more detailed spatiotemporal human movement information is required. Although existing literatures have applied survey data [10,13,14,15,16], social media data [9,17,18,19,20] and mobile phone data [21,22,23,24,25] to record human movement behaviors, several limitations remain. As for the survey data-based methods, the cost of personal monitoring and sampling processes limits the number of samples, thereby reducing the persuasiveness of the estimation results [12,26]. In terms of social media data-based methods, representativeness is still problematic, as the data collected only reflects the individual movement characteristics of people who used the application at a specific time.
With the development of information and communication technology (ICT) and the ubiquity of mobile phones, mobile phone data have become an emerging dataset to measure human mobility. It has the advantages of massive users, wide spatial coverage and passive acquisition capability, which has the potential for long-term and well-represented exposure estimation. However, limited by the cost of data transmission and storage, most mobile phone data are only collected when the individual made a phone call or sent a message [27]. This inevitably limits the temporal granularity and regularity of mobile phone data, which introduces errors in exposure estimation [28,29]. As shown in Figure 1, the difference between real trajectory and recorded trajectory of an individual’s movement locations will inevitably result in incorrect estimation of the duration to which individuals are actually exposed to the areas with different air pollution concentrations. Therefore, improving the temporal granularity of mobile phone data is important for the exposure estimation.
In this paper, a novel method is proposed to estimate individuals’ air pollution spatiotemporal exposure levels using trajectories reconstructed from mobile phone data. The contributions of this paper are outlined as follows:
(1)
We present a novel individual air pollution exposure estimate method. Our method mitigates the gap of spatiotemporal resolution between air pollution monitoring data and mobile phone data, which helps improve the accuracy and reliability of fine-scale air pollution exposure estimation.
(2)
By comparing the three different types of exposure estimates using reconstructed mobile phone trajectories, recorded mobile phone trajectories, and residential locations, we demonstrate the necessity of trajectory reconstruction in exposure estimation.
(3)
Using the city of Shanghai as a case study, we quantitatively analyzed the temporal variations in individual exposures and the spatial distribution of residential areas with high exposure levels using large-scale mobile phone data. It provides a more accurate and comprehensively scientific basis for policy-driven environmental actions and potential health risk reduction.
The rest of this paper is organized as follows. In Section 2, we review the research progress in related fields. The details of the proposed dynamic individual exposure estimation method are then presented in Section 3. The analysis of the case study and results are presented in Section 4; and finally, the conclusions and discussion on this study are presented in the final section.

2. Literature Review

2.1. Air Pollution Exposure Estimates

In the existing studies, researchers have attempted to measure air pollution exposure using various methods. According to the human mobility data that researchers have used, the related studies can be grouped into three categories: fixed-location-based methods, survey data-based methods, and mobile phone data-based methods. The fixed-location-based methods assume that individuals remain at a particular location and thus individual exposure can be calculated by measuring the air pollution concentration at their residence or workplace [30,31,32,33]. However, the spatiotemporal variability in air pollutant concentrations makes the human movement pattern a key factor in air pollution exposure estimation. That is, individual pollutant exposure estimates must be determined based on where the individual stays and how long the individual remains at each location. As neither an individual’s residence nor their workplace fully represents their daily movement pattern, the accuracies of these methods are not satisfactory.
To tackle this problem, researchers began to use survey data to record human movement behavior. The two most common techniques were manual sampling surveys [14,34] and Global Navigation Satellite System (GNSS)-enabled personal monitors [8,10,15,16]. For example, Yoo et al. (2015) used GNSS-equipped monitor data of 43 participants to demonstrate how an individual’s mobility affects personal exposure estimates. Similarly, Park and Kwan (2017) used 80 simulated daily movement trajectories to support the argument that ignoring human mobility patterns may lead to misleading results in exposure assessments. Although these techniques can provide fine-scale human movement data and additional personal information for further research, they have two limitations. First, the time requirements of the sampling process may cause the participants to become bored and result in recall bias [35]. Second, although the sampling method could be well designed, the cost of personal monitors and sampling processes limits the sample size and time period, making the results less persuasive.
The popularization of social media applications and location-based services provided a different way to track individual spatiotemporal activities. Compared with the survey data, the social media data enable us to collect human movement data with larger scale, longer term and broader spatial coverage [36]. Therefore, several researchers used social media data to estimate individual human exposure [9,17,18,19,20]. For instance, Song et al. (2019) combined sparse Weibo (a popular social media application in China) geotagged messages and remote sensing derived PM2.5 concentrations to perform monthly exposure estimations for 13 major cities in China. Yu et al. (2019) demonstrated the potential to estimate individual exposure with Google Maps location data on a minute level. However, the social media data-based method suffered from poor representativeness [37]. That is, the individual exposure inferred from these methods presents the exposure characteristics of active social media users rather than the whole population. The exposure risks of some demographic composition such as elderly and poor people tend to be misestimated since these people use the social media applications less often [38].
With recent advances in positioning techniques and the ubiquity of smart phones, mobile phone data, which can be collected without additional procedure and have wide stratified coverage, have been widely used for air pollution exposure estimation [9,21,23,25]. For example, Dewulf et al. (2016) proposed a method to estimate daily NO2 exposure using the mobility data collected from 5 million mobile phone users in Belgium. Yu et al. (2018) used mobile phone data to measure the influence of human movements on air pollutant exposure estimates. The use of mobile phone data enables us to measure individual exposure for a fairly high proportion of the entire human population over large areas and long periods. However, compared with the air pollution monitoring data, the temporal granularity of mobile phone data is not high enough, which limits the performance of individual exposure estimation [28,29]. Therefore, ways to improve the temporal granularity of mobile phone data has become a potential research hotspot in individual exposure estimation.

2.2. Trajectory Reconstruction from Mobile Phone Data

Trajectory reconstruction aims to approximate the locations of an individual that are missing in a mobile phone dataset. The basic methods of approximation rely on spatiotemporal interpolation [39,40,41]. It is assumed that an individual’s trajectory between two recorded points can be reconstructed using interpolation functions, such as the nearest function or a linear function. Therefore, the missing points can be approximated by the time intervals and distances between their contextual recorded points. These methods have been widely used for the reconstruction of intensive and continuous trajectories. However, the limited resolution of time intervals of mobile phone data and the complex individual movement patterns make these methods unsuitable for reconstructing trajectories from mobile phone data.
As vehicles are the main means of transport in the city [42], researchers have attempted to reconstruct missing points in low-frequency trajectories using map-matching-based methods. These methods assume that individuals’ trajectories are strictly dependent on road networks. Therefore, trajectory points can be snapped to road networks and missing points can be approximated based on movement characteristics along transportation networks such as velocity, acceleration, or travel time [43,44,45,46]. However, the precondition of these methods is questionable, as it ignores the fact that human movement behaviors in urban spaces can occur via the metro or on foot. Thus, the result of these methods is doubtful.
With the development of machine learning and information and communications technology, researchers have focused on pattern mining methods. In these methods, the individual movement patterns are explored using historical trajectories and missing points are approximated by trained machine learning models [47,48,49]. However, limited by the high computational cost and the lack of efficient spatiotemporal proximity analysis methods [50], most pattern mining methods only use one individual’s trajectories or a single trajectory segment to approximate missing points. This leaves only a few amounts of data that can be used for model training, which leads to model overfitting and poor generalization issues. Thus, using trajectories of individuals with similar movement patterns as training data has become an efficient way to enhance reconstruction performance [51].

3. Methodology

The workflow of our method is presented in Figure 2, which includes three parts. First, we introduce a trajectory reconstruction algorithm to mitigate the gap of spatiotemporal resolution between air pollution monitoring data and mobile phone data. Then, we explain our algorithm to estimate air pollution concentrations on a fine scale. Finally, we show how dynamic individual exposure levels are calculated by combining the results of fine-scale air pollution concentration measurements and reconstructed trajectories.

3.1. Anchor-Point Based Trajectory Reconstruction Algorithm

Due to the characteristic of low acquisition cost and wide spatiotemporal coverage, mobile phone data have been widely used in air pollution exposure estimation. However, due to the gap of spatiotemporal resolution between air pollution monitoring data and mobile phone data, trajectory reconstruction is important. Thus, an anchor-point-based trajectory reconstruction algorithm is presented in this subsection to grasp highly dynamic human movement trajectories with corresponding time.

3.1.1. Anchor-Point-Based Clustering

As discussed above, using trajectories of individuals with similar movement patterns as training data has become an efficient way to enhance reconstruction performance. Since anchor points can summarize the key locations of individuals’ movement behaviors well, it becomes an efficient means of understanding the similarities among massive numbers of individuals.
In our research, the anchor point is defined as an area where an individual visit more frequently than a specified number of times [52]. In this process, the number of records on each location in a personal dataset is first calculated. Then, the location with the highest record number is selected and merged with all adjacent locations within a distance threshold α as an examinee [53]. Next, the location with the second highest record number is selected and the same process was repeated until all locations in the personal mobile phone data had been traversed. Finally, the total record number is calculated for each examinee, and the examinees with record number exceeding frequency threshold β of their total record number are detected as anchor points and projected onto a subdistrict. Thus, an individual’s trajectories can be generalized as a collection of anchor points.
After that, based on the respective anchor point collections of individuals, the similarity matrix between the individuals could be calculated with the Jaccard index [54]. Then, the hierarchical clustering method [55] is used to divide the individuals into a series of groups. With this process, the movement patterns of the individuals in each group are similar. Therefore, the model can be trained with all the trajectories in each group rather than using an individual’s trajectories, which helps solve the overfitting problem. The flowchart of the anchor-point-based clustering method is presented in Figure 3.
Note that the choices of distance threshold and frequency threshold affect the clustering result. In terms of the distance threshold, we choose 500 m as the distance threshold for two reasons. First, in our case study area, the average distance between two adjacent locations in the mobile phone datasets is about 240 m. The threshold of 500 m could help reduce signal oscillation problem. Moreover, this threshold has been widely applied in existing studies of anchor point detection and human mobility [53,56]. As for the frequency threshold, it mainly refers to the proportion of the total number of records recorded at key activity locations such as residences and workplaces. Choosing a large threshold will result in the number of user identification anchor points being too small, which makes the user movement characteristics not fully summarized. On the contrary, a small threshold will recognize the unimportant locations as anchor points, which results in additional computational costs. For the data we used in the case study, we found that most of the users have about 20% of their recorded locations at residence or workplace (i.e., the anchor points with the most records during nighttime or daytime hours, respectively), therefore the frequency threshold is set to 20%. It is worth noting that such thresholds might vary in different cities.

3.1.2. Reconstruction of Clustered Trajectories Using a Gradient Boosting Decision Tree Model

After dividing the individuals into a series of clusters, we apply a gradient boosting decision tree (GBDT) model to each cluster to infer the location of a target point. GBDT is an ensemble machine learning model with classification and regression tree (CART) as weak learners [57]. It generates a series of decision trees using the gradient boosting method and the result are determined via the summation of the weak learners [58]. A weak learner indicates a classifier that is only weakly correlated with the true classification [59]. To approximate the location of a missing point, a few characteristics were chosen to describe an individual’s movement patterns. The characteristics are selected based on two aspects. First, the following features were chosen to learn the movement pattern from historical trajectory segments, including the locations and times of contextual points p n − 1 ( x n − 1 , y n − 1 ,   t n − 1 ), p n + 1 ( x n + 1 , y n + 1 , t n − 1 ), and the time of the target point p n ( t n ) ; p n − 1 , p n , and p n + 1 denote three consecutive points along a trajectory that were resampled with different time intervals. Second, the radius of gyration (ROG) and the Shannon entropy (Ent) of each individual were chosen to describe their general movement patterns [24,29], which are defined as follows:
R O G k = 1 o ∑ j = 1 l ( p t j → − p t c → )
E n t k = − ∑ j = 1 o p ( j ) log 2 p ( j )
where o indicates the number of locations recorded in the individual’s dataset, p t j → indicates the j th location, p t c → is the barycenter of all the records ( p t c → = 1 o ∑ j = 1 o p t j → ), and p ( j ) is the visit frequency of the j th location.
Therefore, each training vector record can be represented as [ x n − 1 , y n − 1 , t n − 1 , t n , x n + 1 , y n + 1 , t n + 1 , R O G k , E n t k ] and the label data can be represented as [ x n , y n ]. The architecture of the trajectory reconstruction algorithm is shown in Figure 4. More technical details of the proposed trajectory reconstruction algorithm can be found in the paper [51]. All the mobile phone trajectories were reconstructed with one-hour intervals using the trained model (i.e., 00:00, 01:00, …, 23:00, UTC + 8) to conform with the frequency of air pollution monitoring.

3.2. Estimation of Spatiotemporal Concentrations of Air Pollution

Estimating air pollution concentrations is the second step in quantifying individual air pollution exposure. As a practical technique, extending the ordinary linear regression framework, the geographically weighted regression (GWR) model [60] can examine the spatial variations and nonstationary aspects of a continuous surface of parameters at the local scale and is widely used in air pollution estimation [9,61]. In this study, we applied the GWR model to estimate air pollution concentrations based on the meteorological attributes of the surrounding area and used atmospheric particulate matter with a diameter of fewer than 2.5 micrometers (PM2.5) as a case study. It is a major type of air pollution whose concentration shows seasonal patterns and temporal variability [62] and which can be deposited deep in the lungs through simple respiration, causing an increase in respiratory and cardiovascular diseases [63].
Considering the differences in geographic locations between air pollution monitoring stations and meteorological stations, the datasets need to be made consistent with respect to their spatial domains. Thus, an ordinary Kriging model [64] is first used to estimate meteorological variables with a spatial granularity of 1 km. Then, the meteorological observations within 1 km were averaged and assigned to the corresponding air pollution monitoring station to lessen the spatial mismatch bias. After that, a series of GWR models are developed as shown in Figure 5:
where i and t denote the corresponding location ID and time, P M i , t denotes the PM2.5 concentration (μg/m3), V I S i , t denotes the horizontal visibility (m), W S i , t denotes the wind speed (m/s), T E M i , t denotes the air temperature (°C), and β 0 , i , t , β 1 , i , t , β 2 , i , t , and β 3 , i , t denote the regression coefficients of the corresponding features. In this study, as the temporal granularities of air pollution monitoring data and meteorological observation data are 1 hour, in total, 168 GWR models were trained to approximate the air pollution concentration distribution of each hour during one week. Finally, the finalized GWR models are ascertained based on model performance denoted by fitting the highest coefficient of determination ( R 2 ) and the lowest Akaike information criterion (AIC) value. With the finalized modes, the optimal coefficients are used to approximate the air pollution concentration distribution in the study area with a spatial granularity of 1 km.
It is worth noting that there are two reasons for our choice of a spatial granularity of 1 km. First, since air pollution exposure estimation is a typical study of the environmental influences on individual behaviors, its spatiotemporal resolution is limited by the lower resolution data of the human movement data and the air pollution concentration data. Considering that the positioning errors of the mobile phone data in urban spaces range from 100 m to 1000 m [65,66], we chose a 1-km grid to divide the space. Second, this spatial granularity has also been widely used in previous studies to estimate pollutant concentrations and individual exposures to them [10].

3.3. Dynamic Individual Exposure Calculation

Since individuals’ locations and corresponding air pollution concentrations vary in both space and time, we propose an algorithm to incorporate dynamic individual locations, the spatiotemporal variation in air pollution concentrations, and the microenvironment effect to estimate the dynamic individual exposure as follows:
E x p j = ∑ t = 1 T ∑ n = 1 N A P i , t ∗ M E i , n , t ∗ T P i , n , t
where E x p j denotes the dynamic exposure of individual j , i and t denote the corresponding location ID and time, N denotes the total number of microenvironments experienced by individual j within a specified temporal window (e.g., an hour) and n denotes the n t h microenvironment, A P i , t denotes the outdoor air pollution concentration, M E i , n , t denotes the ratio of the air pollution concentration in the n t h microenvironment to the outdoor air pollution concentration, and T P i , n , t denotes the percentage of time that the individual stayed in the n t h microenvironment.
However, research on the impact of air pollution concentrations on the microenvironment has remained at the stage of qualitative analysis of small sample data [10,14]. The temporal resolution requirements and costly monitors make it impossible to acquire accurate observations on a large scale. Moreover, many factors, such as ventilation, air conditioning, smoking, and cooking, are independent of the outdoor environment, but can influence an individual’s microenvironment [67,68,69]. That is, individual exposure can be quite different at the same time and in the same area [35]. Thus, we define air pollution exposure as the outdoor exposure level in this study and simplify the ideal algorithm in Equation (3) to suit the estimation of large-scale dynamic individual exposure using the following equation:
E x p j ′ = ∑ t = 1 T ∑ n = 1 N A P i , t
where E x p j ′ denotes the cumulative individual exposure from the simplified algorithm by only considering outdoor air pollution concentration and individual exposure duration.

4. Case Study

4.1. Data

4.1.1. Mobile Phone Data

The mobile phone data we used were provided anonymously by a mobile network operator through a joint research cooperation. It records the location trajectories of over a million individuals for seven consecutive days in the city of Shanghai, China. A map of case study area is shown in Figure 6. According to the Shanghai Statistical Yearbook 2018, the operator accounts for about 56% of the city’s residents and is widely distributed among all strata of society (Shanghai municipal statistics bureau 2019). The dataset contained call detail records (CDRs, i.e., phone calls and text message) and actively generated records (i.e., regular updates, periodic updates, and cellular handovers). To decrease the signal oscillation problem, a repetition suppression algorithm [70] was used for data preprocessing. Table 1 shows an instance of one individual’s trajectory data. It is worth noting that none of the personal identifiable information (i.e., name, gender, phone number) were provided to protect the individuals’ privacy. In addition, the locations of the trajectory points were projected to the locations of ambient cell phone towers. In other words, there is still a gap between the user’s actual location and the projected trajectory points, which is about 100 m–1000 m on average in the dataset dependent on the specific area (e.g., downtown vs. suburbs). The distribution of time intervals between two adjacent call detail records is shown in Figure 7.

4.1.2. Environmental Data

Hourly ground-station PM2.5 concentration data (in μm/m3) were collected from the data center of the Ministry of Environmental Protection of the People’s Republic of China (http://datacenter.mee.gov.cn) and the World Air Quality Index project (http://aqicn.org/city/shanghai/). Moreover, ground-station meteorological variables, including horizontal visibility, air temperature, and surface wind speed were collected from the National Meteorological Information Center (http://data.cma.cn/). The instances of ground-station PM2.5 Concentration Data and meteorological Data are shown in Table 2 and Table 3.
In line with the mobile phone data, environmental data were modeled during the corresponding days. To mitigate the estimation biases and improve spatial interpolation accuracy at marginal areas of the study area, the environmental data from two neighboring provinces, Zhejiang and Jiangsu, were also adopted for air pollution concentration estimation to add sufficient surrounding information for the locations at the boundary of study areas. Thus, the data from a total of 178 monitoring stations and 150 meteorological stations were collected for monitoring ambient air quality.

4.2. Spatiotemporal Variability in PM2.5 Concentration

The spatiotemporal PM2.5 concentration variation is an important component of the individual exposure estimation. Figure 8 shows examples of the extracted data (collected on a workday and a weekend) on selected hourly maps (e.g., 04:00, 10:00, 16:00, and 22:00) and the temporal variation curves for PM2.5 concentrations, which were approximated by GWR models mentioned in Section 3.2. It clearly reveals the spatiotemporal variations in PM2.5 concentrations in the study area, suggesting that the results of individual exposure estimates are directly related to individual movement behaviors across locations and time. Therefore, temporal mismatch between human movement data and PM2.5 concentrations will inevitably introduce errors into the individual exposure estimation. Thus, the requirement of consistent recording times of human movement data and PM2.5 concentration data makes trajectory reconstruction a key preprocessing step for more accurate individual exposure estimation.
The fitting results of the GWR models were evaluated using the metrics of R-squared and root mean square error (RMSE). The average R-squared between the predicted and observed PM2.5 concentrations is 0.81, and the average RMSE is 21.18 μg/m3, both of which are acceptable for dynamic air pollution exposure estimation [9].

4.3. Performance Evaluation of the Trajectory Reconstruction Algorithm

In our proposed method, the accuracy of air pollution exposure estimation directly depends on the accuracy of trajectory reconstruction. Thus, we compared the performance of our proposed trajectory reconstruction algorithm with other existing methods. To verify the performance of the proposed trajectory reconstruction algorithm, all the CDR data were extracted and used as “recorded data” (about 34% of the total records) and the actively generated records were used as “missing data” (about 66% of the total records). Thus, the average performance of reconstruction algorithms could be evaluated using the mean absolute error (MAE) between the reconstructed locations and the actively generated record locations, and the stability of the reconstruction algorithms could be evaluated using the standard deviation of the errors. In addition, one artificial neural network-based reconstruction algorithm, ANN-TR [49], and two most widely used algorithms, nearest interpolation and linear interpolation [40], were chosen as baselines. The performance comparison of these algorithms is shown in Figure 9.
As shown in Figure 9, our proposed algorithm shows lower reconstruction error and better robustness (with lower StDev) than baselines, which proves the superiority of our proposed method over these existing ones. More technical details about different trajectory reconstruction evaluation results can be found in [51].

4.4. Comparison with Existing Exposure Estimate Methods

To quantitatively analyze how the gap of spatiotemporal resolution between human movement data and air pollution monitoring data affects individual exposure estimates, three types of individual exposure estimates were obtained by (1) using reconstructed mobile phone trajectories for exposure estimation (hereafter, TR-EE); (2) using recorded mobile phone trajectories for exposure estimation (hereafter, REC-EE); and (3) using static locations (home location) for exposure estimation (hereafter, SL-EE). As the statistical results of individual exposure levels do not conform to a normal distribution, the K-S test [71] was applied to access the differences among these three types of exposure estimates, which are shown in Table 4. The results show that all the estimated pairs have larger K-S statistics than the expected values under null hypothesis and with very low p-values, which indicates that the exposure estimates of TR-EE are significantly different from those of the other two.
Figure 10 shows the difference between TR-EE and the other two types of exposure estimates with a box plot. The lengths of the interquartile range boxes of the TR-EE & SL-EE cases are larger than those of the boxes for the TR-EE & REC-EE cases on both weekdays and weekends. This means that the exposure estimation errors for the both middle half of the individuals and all the individuals increase when the home locations are used for exposure estimation than when the recorded mobile phone trajectories are used. In addition, we can see that the median values (the red lines in each box) of the four cases are approximately zero and the boxes are basically symmetrical around the median line. This indicates that the results of exposure estimation can be either overestimated or underestimated when the recorded mobile phone trajectories or static locations are used. It demonstrates that neither the home location nor the recorded mobile phone data can fully represent individuals’ movement behaviors when estimating dynamic individual air pollution exposure. Therefore, the necessity of trajectory reconstruction when using mobile phone data to estimate individual exposure levels is clarified. Further, due to varying patterns of individual movement behaviors on weekends (with the average Shannon entropy of 1.40) compared to that on weekdays (with the average Shannon entropy of 1.45), the exposure estimation errors on weekends are slightly different from those on workdays.

4.5. Potential Health Effects

Based on the ambient air quality standards proposed by the Chinese National Environmental Protection Agency (GB 3095–2012), we classified the PM2.5 concentration into five categories to show the potential health effects of each category, as shown in Table 5. We measured the potential health effects of air pollution from two perspectives: the individual-oriented and geographical space–oriented.
For the individual-oriented perspective, we selected two individuals ( I a and I b ), whose mobilities were considerably different from each other, to show how individual movement behaviors affect exposure estimates and potential health effects during a day. The space-time paths represent the individuals’ trajectories in 3D space, and the colors of the trajectory segments represent the corresponding potential health effects [8]. As shown in Figure 11, both individuals should take appropriate precautions when spending their time outdoors. Individual I b has longer red-colored and purple-colored segments than does individual I a , which means that individual I b experiences PM2.5 concentrations over moderately polluted levels for longer times. If individual I b belongs to the sensitive crowd, the necessary precautions should be taken or avoid outdoor activities. With our method, we can obtain fine spatiotemporal resolution of exposure-level trajectories for each individual more accurately, so as to enhance the individuals’ awareness of protection and reduce their risk of respiratory or heart diseases during trip planning.
As for the geographical space–oriented perspective, understanding the residences of people with high air pollution exposures can help government organizations prioritize resources and strengthen publicity purposefully to reduce the potential health effects. Therefore, we selected individuals with high exposure levels to air pollutions (whose total exposure levels during the whole week were in the top 20% of all individuals) and speculated their residential locations using their most frequently visited location during the nighttime (from 22:00 to 06:00) [72,73]. These possible residential locations of mobile phone users were aggregated to subdistrict levels, and the results are shown in Figure 12.
We can see that the individuals with high exposure levels mostly lived in the western part of Shanghai, and further concentrated in the Songjiang District, Jiading District, and Qingpu District. For more details, we chose the top five subdistricts in which the highest exposure individuals lived and calculated the percentage of time for which residents were away from their residences and the corresponding percentage of exposure away from residence, as shown in Figure 13. The results show that there is a positive linear relationship between the percentage of residents’ time away from their residence and the corresponding percentage of their exposure, which indicates that people with a large percentage of movement time (due to a huge potential of movements in polluted areas) tend to have large air pollution exposures, which are not only restricted to their residential areas. Table 6 shows that the residents in these areas were exposed to air pollution levels above the lightly polluted level for more than 50% of the time and were exposed to severely polluted levels for about 5% of the time, which indicated that these residents may experience irritation and the sensitive individuals among them may experience serious conditions. These results could help government organizations prioritize resources in terms of air pollution issues geographically and provided technical and theoretical support for policy-driven environmental actions.

5. Conclusions

This study proposed a method for estimating dynamic individual air pollution exposures using trajectories reconstructed from mobile phone data. This method mitigates the gap of spatiotemporal resolution between human movement data and air pollution monitoring data, thereby assisting in the estimation of individuals’ air pollution exposures more accurately and comprehensively at a high spatiotemporal resolution. Using the city of Shanghai as a case study, we compared three different types of exposure estimates obtained via (1) reconstructed mobile phone trajectories, (2) recorded mobile phone trajectories, and (3) residential locations. The results show that exposure estimates using reconstructed mobile phone trajectories are significantly different from the other two types of estimates, and such differences are higher when using recorded mobile phone trajectories than when using home locations for exposure estimates. This demonstrated the necessity of trajectory reconstruction in exposure and health risk assessments. Additionally, we measured the potential health effects of air pollution from both individual and geographical perspectives, which helped reveal the temporal variations in individual exposure levels and the spatial distribution of residential areas with high exposure levels.
By using the reconstructed mobile phone trajectories to measure human behaviors, this method provided a more accurate and comprehensive way of estimating dynamic individual air pollution exposure levels across space over time, which can help support policy-driven environmental actions and reduce potential health risks. In addition to the PM2.5 concentrations shown in the case study, the proposed method can be used to estimate individual exposure to other pollutants, such as NO2, SO2, and noise. Further, by updating individual mobile phone data and ambient pollution data, our method can be applied both to near-real-time estimates for individuals who are exposed to poor ambient air quality and the long-term effects of ambient pollution on human health, which could contribute to many crucial applications, such as disease surveillance and disaster loss assessment.
However, this study has several limitations and requires further exploration. First, the individual exposure estimation method suffered from an uncertain geographic context problem (UGCoP). Wherein, two kinds of contextual factors bring uncertainties to exposure estimates: the spatial configuration of corresponding spatial units and the timing and duration of exposure to those units [26,74,75]. For example, although we improve the spatiotemporal granularity of human movement data to the hourly level, an individual may have multiple activities in one hour, which makes the estimation counter-intuitive. Whereas, the proposed method for air pollution estimation using trajectories reconstructed from mobile phone data provides an alternative way to mitigate the influence of UGCoP. More accurate trajectory reconstruction algorithms and human movement data with better spatiotemporal resolution will be applied in this field to further improve the individual exposure estimation performance. In addition, the accuracy of air pollutant estimation is another important factor affecting the accuracy of individual exposure estimates. In the proposed method, we only adopted a GWR model to estimate air pollution concentrations by incorporating meteorological variables across space and air pollution values recorded by air quality monitoring stations. Considering the spatiotemporal coverage and resolution of meteorological observation stations, the satellite-based air pollutant concentration estimation algorithm might be helpful for improving the exposure estimate accuracy from another point of view. In addition, as the purpose of this study was to estimate large-scale dynamic individual exposures and lacking of quantifying methods for determining the effects on the microenvironment, the influence of the microenvironment (as mentioned in Equation (3)) was ignored in this study. How to consider the microenvironment effect in the proposed method is an interesting topic. We will focus on these issues in future work.

Author Contributions

Conceptualization, M.L.; Funding acquisition, F.L.; Methodology, M.L. and S.G.; Resources, H.T. and H.Z.; Supervision, F.L. and S.G.; Visualization, M.L.; Writing—original draft, M.L.; Writing—review & editing, M.L., S.G., F.L., H.T. and H.Z.

Funding

This work was supported by National Natural Science Foundation of China, grant number 41771436; National Key Research and Development Program, grant number 2016YFB0502104.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Brunekreef, B.; Holgate, S.T. Air pollution and health. Lancet 2002, 360, 1233–1242. [Google Scholar] [CrossRef]
  2. Mannucci, P.M.; Franchini, M. Health Effects of Ambient Air Pollution in Developing Countries. Int. J. Environ. Res. Public Health 2017, 14, 1048. [Google Scholar] [CrossRef] [PubMed]
  3. Dominici, F.; Peng, R.D.; Bell, M.L.; Pham, L.; McDermott, A.; Zeger, S.L.; Samet, J.M. Fine particulate air pollution and hospital admission for cardiovascular and respiratory diseases. JAMA 2006, 295, 1127–1134. [Google Scholar] [CrossRef] [PubMed]
  4. Goldberg, M. A Systematic Review of the Relation between Long-term Exposure to Ambient Air Pollution and Chronic Diseases. Rev. Environ. Health 2011, 23, 243–298. [Google Scholar] [CrossRef] [PubMed]
  5. Cohen, A.J.; Brauer, M.; Burnett, R.; Anderson, H.R.; Frostad, J.; Estep, K.; Balakrishnan, K.; Brunekreef, B.; Dandona, L.; Dandona, R.; et al. Estimates and 25-year trends of the global burden of disease attributable to ambient air pollution: An analysis of data from the Global Burden of Diseases Study 2015. Lancet 2017, 389, 1907–1918. [Google Scholar] [CrossRef]
  6. Lelieveld, J.; Klingmüller, K.; Pozzer, A.; Pöschl, U.; Fnais, M.; Daiber, A.; Münzel, T. Cardiovascular disease burden from ambient air pollution in Europe reassessed using novel hazard ratio functions. Eur. Heart J. 2019, 40, 1590–1596. [Google Scholar] [CrossRef]
  7. Kwan, M.-P.; Liu, D.; Vogliano, J. Assessing Dynamic Exposure to Air Pollution. In Space-Time Integration in Geography and GIScience: Research Frontiers in the US and China; Kwan, M.-P., Richardson, D., Wang, D., Zhou, C., Eds.; Springer: Dordrecht, The Netherlands, 2015; pp. 283–300. ISBN 978-94-017-9205-9. [Google Scholar]
  8. Park, Y.M.; Kwan, M.-P. Individual exposure estimates may be erroneous when spatiotemporal variability of air pollution and human mobility are ignored. Health Place 2017, 43, 85–94. [Google Scholar] [CrossRef]
  9. Chen, B.; Song, Y.; Jiang, T.; Chen, Z.; Huang, B.; Xu, B. Real-Time Estimation of Population Exposure to PM2.5 Using Mobile- and Station-Based Big Data. Int. J. Environ. Res. Public Health 2018, 15, 573. [Google Scholar] [CrossRef]
  10. Yoo, E.; Rudra, C.; Glasgow, M.; Mu, L. Geospatial Estimation of Individual Exposure to Air Pollutants: Moving from Static Monitoring to Activity-Based Dynamic Exposure Assessment. Ann. Assoc. Am. Geogr. 2015, 105, 1–12. [Google Scholar] [CrossRef]
  11. Xie, X.; Semanjski, I.; Gautama, S.; Tsiligianni, E.; Deligiannis, N.; Rajan, R.T.; Pasveer, F.; Philips, W. A Review of Urban Air Pollution Monitoring and Exposure Assessment Methods. ISPRS Int. J. Geo-Inf. 2017, 6, 389. [Google Scholar] [CrossRef]
  12. Mennis, J.; Yoo, E.-H.E. Geographic information science and the analysis of place and health. Trans. GIS 2018, 22, 842–854. [Google Scholar] [CrossRef] [PubMed]
  13. Kwan, M.-P. From place-based to people-based exposure measures. Soc. Sci. Med. 2009, 69, 1311–1313. [Google Scholar] [CrossRef] [PubMed]
  14. Lioy, P.J. Exposure Science: A View of the Past and Milestones for the Future. Environ. Health Perspect. 2010, 118, 1081–1090. [Google Scholar] [CrossRef] [PubMed]
  15. Glasgow, M.L.; Rudra, C.B.; Yoo, E.-H.; Demirbas, M.; Merriman, J.; Nayak, P.; Crabtree-Ide, C.; Szpiro, A.A.; Rudra, A.; Wactawski-Wende, J.; et al. Using smartphones to collect time-activity data for long-term personal-level air pollution exposure assessment. J. Expo. Sci. Environ. Epidemiol. 2016, 26, 356–364. [Google Scholar] [CrossRef]
  16. Mihăiţă, A.S.; Dupont, L.; Chery, O.; Camargo, M.; Cai, C. Evaluating air quality by combining stationary, smart mobile pollution monitoring and data-driven modelling. J. Clean. Prod. 2019, 221, 398–418. [Google Scholar] [CrossRef]
  17. Ford, B.; Burke, M.; Lassman, W.; Pfister, G.; Pierce, J.R. Status update: Is smoke on your mind? Using social media to assess smoke exposure. Atmos. Chem. Phys. Discuss. 2017, 17, 7541–7554. [Google Scholar] [CrossRef]
  18. Song, Y.; Huang, B.; Cai, J.; Chen, B. Dynamic assessments of population exposure to urban greenspace using multi-source big data. Sci. Total Environ. 2018, 634, 1315–1325. [Google Scholar] [CrossRef]
  19. Song, Y.; Huang, B.; He, Q.; Chen, B.; Wei, J.; Mahmood, R. Dynamic assessment of PM2.5 exposure and health risk using remote sensing and geo-spatial big data. Environ. Pollut. 2019, 253, 288–296. [Google Scholar] [CrossRef]
  20. Yu, X.; Stuart, A.L.; Liu, Y.; Ivey, C.E.; Russell, A.G.; Kan, H.; Henneman, L.R.; Sarnat, S.E.; Hasan, S.; Sadmani, A.; et al. On the accuracy and potential of Google Maps location history data to characterize individual mobility for air pollution health studies. Environ. Pollut. 2019, 252, 924–930. [Google Scholar] [CrossRef]
  21. Dewulf, B.; Neutens, T.; Lefebvre, W.; Seynaeve, G.; Vanpoucke, C.; Beckx, C.; Van De Weghe, N. Dynamic assessment of exposure to air pollution using mobile phone data. Int. J. Health Geogr. 2016, 15, 14. [Google Scholar] [CrossRef]
  22. Birenboim, A.; Shoval, N. Mobility Research in the Age of the Smartphone. Ann. Am. Assoc. Geogr. 2016, 106, 1–9. [Google Scholar] [CrossRef]
  23. Yu, H.; Russell, A.; Mulholland, J.; Huang, Z. Using cell phone location to assess misclassification errors in air pollution exposure estimation. Environ. Pollut. 2018, 233, 261–266. [Google Scholar] [CrossRef] [PubMed]
  24. Jiang, J.; Li, Q.; Tu, W.; Shaw, S.-L.; Yue, Y. A simple and direct method to analyse the influences of sampling fractions on modelling intra-city human mobility. Int. J. Geogr. Inf. Sci. 2019, 33, 618–644. [Google Scholar] [CrossRef]
  25. Xu, Y.; Jiang, S.; Li, R.; Zhang, J.; Zhao, J.; Abbar, S.; González, M.C. Unraveling environmental justice in ambient PM2.5 exposure in Beijing: A big data approach. Comput. Environ. Urban Syst. 2019, 75, 12–21. [Google Scholar] [CrossRef]
  26. Kwan, M.-P. The Uncertain Geographic Context Problem. Ann. Assoc. Am. Geogr. 2012, 102, 958–968. [Google Scholar] [CrossRef]
  27. Song, C.; Qu, Z.; Blumm, N.; Barabási, A.-L. Limits of Predictability in Human Mobility. Science 2010, 327, 1018–1021. [Google Scholar] [CrossRef]
  28. Ranjan, G.; Zang, H.; Zhang, Z.-L.; Bolot, J. Are call detail records biased for sampling human mobility? SIGMOBILE Mob. Comput. Commun. Rev. 2012, 16, 33. [Google Scholar] [CrossRef]
  29. Zhao, Z.; Shaw, S.-L.; Xu, Y.; Lu, F.; Chen, J.; Yin, L. Understanding the bias of call detail records in human mobility research. Int. J. Geogr. Inf. Sci. 2016, 30, 1–25. [Google Scholar] [CrossRef]
  30. Ashmore, M.; Dimitroulopoulou, C. Personal exposure of children to air pollution. Atmos. Environ. 2009, 43, 128–141. [Google Scholar] [CrossRef]
  31. Du, X.; Kong, Q.; Ge, W.; Zhang, S.; Fu, L. Characterization of personal exposure concentration of fine particles for adults and children exposed to high ambient concentrations in Beijing, China. J. Environ. Sci. 2010, 22, 1757–1764. [Google Scholar] [CrossRef]
  32. Hu, X.; Waller, L.A.; Lyapustin, A.; Wang, Y.; Liu, Y. 10-year spatial and temporal trends of PM2.5 concentrations in the southeastern US estimated using high-resolution satellite data. Atmos. Chem. Phys. Discuss. 2014, 14, 6301–6314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Shafran-Nathan, R.; Yuval; Levy, I.; Broday, D.M. Exposure estimation errors to nitrogen oxides on a population scale due to daytime activity away from home. Sci. Total Environ. 2017, 580, 1401–1409. [Google Scholar] [CrossRef] [PubMed]
  34. Mitchell, C.S.; Zhang, J.; Sigsgaard, T.; Jantunen, M.; Lioy, P.J.; Samson, R.; Karol, M.H. Current State of the Science: Health Effects and Indoor Environmental Quality. Environ. Health Perspect. 2007, 115, 958–964. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Steinle, S.; Reis, S.; Sabel, C.E. Quantifying human exposure to air pollution—Moving from static monitoring to spatio-temporally resolved personal exposure assessment. Sci. Total Environ. 2013, 443, 184–193. [Google Scholar] [CrossRef] [Green Version]
  36. Hawelka, B.; Sitko, I.; Beinat, E.; Sobolevsky, S.; Kazakopoulos, P.; Ratti, C. Geo-located Twitter as proxy for global mobility patterns. Cartogr. Geogr. Inf. Sci. 2014, 41, 260–271. [Google Scholar] [CrossRef] [Green Version]
  37. Zagheni, E.; Weber, I. Demographic research with non-representative internet data. Int. J. Manpow. 2015, 36, 13–25. [Google Scholar] [CrossRef]
  38. Mellon, J.; Prosser, C. Twitter and Facebook are not representative of the general population: Political attitudes and demographics of British social media users. Res. Polit. 2017, 4, 2053168017720008. [Google Scholar] [CrossRef]
  39. Ficek, M.; Kencl, L. Inter-Call Mobility model: A spatio-temporal refinement of Call Data Records using a Gaussian mixture model. In Proceedings of the 2012 Proceedings IEEE INFOCOM, Orlando, FL, USA, 25–30 March 2012; pp. 469–477. [Google Scholar]
  40. Hoteit, S.; Secci, S.; Sobolevsky, S.; Ratti, C.; Pujolle, G. Estimating human trajectories and hotspots through mobile phone data. Comput. Netw. 2014, 64, 296–307. [Google Scholar] [CrossRef] [Green Version]
  41. Perera, K.; Bhattacharya, T.; Kulik, L.; Bailey, J. Trajectory inference for mobile devices using connected cell towers. In Proceedings of the 23rd SIGSPATIAL International Conference on Advances in Geographic Information Systems, Seattle, WA, USA, 3–6 November 2015; pp. 1–10. [Google Scholar]
  42. Shanghai Urban and Rural Construction and Transportation Development Research Institute. The main results of the fifth comprehensive traffic survey in Shanghai. Traffic Transp. 2015, 182, 23–26. [Google Scholar]
  43. Ni, D.; Wang, H. Trajectory Reconstruction for Travel Time Estimation. J. Intell. Transp. Syst. 2008, 12, 113–125. [Google Scholar] [CrossRef]
  44. Jagadeesh, G.R.; Srikanthan, T. Probabilistic Map Matching of Sparse and Noisy Smartphone Location Data. In Proceedings of the 2015 IEEE 18th International Conference on Intelligent Transportation Systems, Gran Canaria, Spain, 15–18 September 2015; pp. 812–817. [Google Scholar]
  45. Schulze, G.; Horn, C.; Kern, R. Map-Matching Cell Phone Trajectories of Low Spatial and Temporal Accuracy. In Proceedings of the 2015 IEEE 18th International Conference on Intelligent Transportation Systems, Gran Canaria, Spain, 15–18 September 2015; pp. 2707–2714. [Google Scholar]
  46. Algizawy, E.; Ogawa, T.; El-Mahdy, A. Real-Time Large-Scale Map Matching Using Mobile Phone Data. ACM Trans. Knowl. Discov. Data 2017, 11, 1–38. [Google Scholar] [CrossRef]
  47. Fan, Z.; Arai, A.; Song, X.; Witayangkurn, A.; Kanasugi, H.; Shibasaki, R. A Collaborative Filtering Approach to Citywide Human Mobility Completion from Sparse Call Records. In Proceedings of the Twenty-Fifth International Joint Conference on Artificial Intelligence, New York, NY, USA, 9–15 July 2016; pp. 2500–2506. [Google Scholar]
  48. Chen, G.; Viana, A.C.; Sarraute, C. Towards an adaptive completion of sparse Call Detail Records for mobility analysis. In Proceedings of the 2017 IEEE International Conference on Pervasive Computing and Communications Workshops (PerCom Workshops), Kona, HI, USA, 13–17 March 2017; pp. 302–305. [Google Scholar]
  49. Liu, Z.; Ma, T.; Du, Y.; Pei, T.; Yi, J.; Peng, H. Mapping hourly dynamics of urban population using trajectories reconstructed from mobile phone records. Trans. GIS 2018, 22, 494–513. [Google Scholar] [CrossRef]
  50. Yuan, H.; Chen, B.Y.; Li, Q.; Shaw, S.-L.; Lam, W.H.K. Toward space-time buffering for spatiotemporal proximity analysis of movement data. Int. J. Geogr. Inf. Sci. 2018, 32, 1211–1246. [Google Scholar] [CrossRef]
  51. Li, M.; Gao, S.; Lu, F.; Zhang, H. Reconstruction of human movement trajectories from large-scale low-frequency mobile phone data. Comput. Environ. Urban Syst. 2019, 77, 101346. [Google Scholar] [CrossRef]
  52. Xu, Y.; Shaw, S.-L.; Fang, Z.; Yin, L. Estimating Potential Demand of Bicycle Trips from Mobile Phone Data—An Anchor-Point Based Approach. ISPRS Int. J. Geo-Inf. 2016, 5, 131. [Google Scholar] [CrossRef] [Green Version]
  53. Long, Y.; Thill, J.-C. Combining smart card data and household travel survey to analyze jobs–housing relationships in Beijing. Comput. Environ. Urban Syst. 2015, 53, 19–35. [Google Scholar] [CrossRef] [Green Version]
  54. Jaccard, P. The distribution of the flora in the alpine zone. New Phytol. 1912, 11, 37–50. [Google Scholar] [CrossRef]
  55. Gil-Garcia, R.; Badia-Contelles, J.; Pons-Porrata, A. A General Framework for Agglomerative Hierarchical Clustering Algorithms. In Proceedings of the 18th International Conference on Pattern Recognition (ICPR’06), Hong Kong, China, 20–24 August 2006; pp. 569–572. [Google Scholar]
  56. Chen, J.; Pei, T.; Shaw, S.-L.; Lu, F.; Li, M.; Cheng, S.; Liu, X.; Zhang, H. Fine-grained prediction of urban population using mobile phone location data. Int. J. Geogr. Inf. Sci. 2018, 32, 1770–1786. [Google Scholar] [CrossRef]
  57. Breiman, L.; Friedman, J.H.; Olshen, R.A.; Stone, C.J. Classification and Regression Trees; Routledge: London, UK, 2017. [Google Scholar]
  58. Friedman, J.H. Greedy Function Approximation: A Gradient Boosting Machine. Ann. Stat. 2001, 29, 1189–1232. [Google Scholar] [CrossRef]
  59. Kearns, M.; Valiant, L. Cryptographic limitations on learning Boolean formulae and finite automata. J. ACM 1994, 41, 67–95. [Google Scholar] [CrossRef]
  60. Brunsdon, C.; Fotheringham, A.S.; Charlton, M.E. Geographically Weighted Regression: A Method for Exploring Spatial Nonstationarity. Geogr. Anal. 1996, 28, 281–298. [Google Scholar] [CrossRef]
  61. You, W.; Zang, Z.; Zhang, L.; Li, Y.; Pan, X.; Wang, W. National-Scale Estimates of Ground-Level PM2.5 Concentration in China Using Geographically Weighted Regression Based on 3 km Resolution MODIS AOD. Remote Sens. 2016, 8, 184. [Google Scholar] [CrossRef] [Green Version]
  62. Zhao, N.; Liu, Y.; Vanos, J.K.; Cao, G. Day-of-week and seasonal patterns of PM2.5 concentrations over the United States: Time-series analyses using the Prophet procedure. Atmos. Environ. 2018, 192, 116–127. [Google Scholar] [CrossRef]
  63. Xu, Q.; Li, X.; Wang, S.; Wang, C.; Huang, F.; Gao, Q.; Wu, L.; Tao, L.; Guo, J.; Wang, W.; et al. Fine Particulate Air Pollution and Hospital Emergency Room Visits for Respiratory Disease in Urban Areas in Beijing, China, in 2013. PLoS ONE 2016, 11, e0153099. [Google Scholar] [CrossRef] [PubMed]
  64. Oliver, M.A.; Webster, R. Kriging: A method of interpolation for geographical information systems. Int. J. Geogr. Inf. Syst. 1990, 4, 313–332. [Google Scholar] [CrossRef]
  65. Li, M.; Lu, F.; Zhang, H.; Chen, J. Predicting future locations of moving objects with deep fuzzy-LSTM networks. Transp. A Transp. Sci. 2018, 1–18. [Google Scholar] [CrossRef]
  66. Gao, S.; Liu, Y.; Wang, Y.; Ma, X. Discovering Spatial Interaction Communities from Mobile Phone Data. Trans. GIS 2013, 17, 463–481. [Google Scholar] [CrossRef] [Green Version]
  67. Braniš, M. Personal Exposure Measurements. In Human Exposure to Pollutants via Dermal Absorption and Inhalation; Lazaridis, M., Colbeck, I., Eds.; Environmental Pollution; Springer: Dordrecht, The Netherlands, 2010; pp. 97–141. ISBN 978-90-481-8663-1. [Google Scholar]
  68. Lai, H.; Bayer-Oglesby, L.; Colvile, R.; Götschi, T.; Jantunen, M.; Künzli, N.; Kulinskaya, E.; Schweizer, C.; Nieuwenhuijsen, M. Determinants of indoor air concentrations of PM2.5, black smoke and NO2 in six European cities (EXPOLIS study). Atmos. Environ. 2006, 40, 1299–1313. [Google Scholar] [CrossRef]
  69. Franklin, P.J. Indoor air quality and respiratory health of children. Paediatr. Respir. Rev. 2007, 8, 281–286. [Google Scholar] [CrossRef]
  70. Fiadino, P.; Valerio, D.; Ricciato, F.; Hummel, K.A. Steps towards the Extraction of Vehicular Mobility Patterns from 3G Signaling Data. In Traffic Monitoring and Analysis; Springer: Berlin/Heidelberg, Germany, 2012; Volume 7189, pp. 66–80. [Google Scholar]
  71. Smirnov, N. Table for Estimating the Goodness of Fit of Empirical Distributions. Ann. Math. Stat. 1948, 19, 279–281. [Google Scholar] [CrossRef]
  72. Calabrese, F.; Diao, M.; Di Lorenzo, G.; Ferreira, J.; Ratti, C. Understanding individual mobility patterns from urban sensing data: A mobile phone trace example. Transp. Res. Part C Emerg. Technol. 2013, 26, 301–313. [Google Scholar] [CrossRef]
  73. Kung, K.S.; Greco, K.; Sobolevsky, S.; Ratti, C. Exploring Universal Patterns in Human Home-Work Commuting from Mobile Phone Data. PLoS ONE 2014, 9, e96180. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Kwan, M.-P. How GIS can help address the uncertain geographic context problem in social science research. Ann. GIS 2012, 18, 245–255. [Google Scholar] [CrossRef]
  75. Chen, X.; Kwan, M.-P. Contextual Uncertainties, Human Mobility, and Perceived Food Environment: The Uncertain Geographic Context Problem in Food Access Research. Am. J. Public Health 2015, 105, 1734–1737. [Google Scholar] [CrossRef]
Ijerph 16 04522 g001 550
Figure 1. Illustration of the influence of the data sparsity problem in human movement data.
Figure 1. Illustration of the influence of the data sparsity problem in human movement data.
Ijerph 16 04522 g001
Ijerph 16 04522 g002 550
Figure 2. The workflow of the dynamic individual exposure estimation method.
Figure 2. The workflow of the dynamic individual exposure estimation method.
Ijerph 16 04522 g002
Ijerph 16 04522 g003 550
Figure 3. Flowchart of the anchor-point-based clustering method.
Figure 3. Flowchart of the anchor-point-based clustering method.
Ijerph 16 04522 g003
Ijerph 16 04522 g004 550
Figure 4. Architecture of the trajectory reconstruction algorithm.
Figure 4. Architecture of the trajectory reconstruction algorithm.
Ijerph 16 04522 g004
Ijerph 16 04522 g005 550
Figure 5. Illustration of the geographically weighted regression (GWR) model.
Figure 5. Illustration of the geographically weighted regression (GWR) model.
Ijerph 16 04522 g005
Ijerph 16 04522 g006 550
Figure 6. Map of the study area.
Figure 6. Map of the study area.
Ijerph 16 04522 g006
Ijerph 16 04522 g007 550
Figure 7. The distributions of the probability density function (PDF) and the cumulative distribution function (CDF) of time intervals between two adjacent call detail records in the recorded mobile phone data.
Figure 7. The distributions of the probability density function (PDF) and the cumulative distribution function (CDF) of time intervals between two adjacent call detail records in the recorded mobile phone data.
Ijerph 16 04522 g007
Ijerph 16 04522 g008 550
Figure 8. Different facets of PM2.5 concentration. (a) Hourly maps of PM2.5 concentration distribution on a workday and a weekend. (b) Temporal variation of PM2.5 concentration on a workday. (c) Temporal variation of PM2.5 concentration on a weekend.
Figure 8. Different facets of PM2.5 concentration. (a) Hourly maps of PM2.5 concentration distribution on a workday and a weekend. (b) Temporal variation of PM2.5 concentration on a workday. (c) Temporal variation of PM2.5 concentration on a weekend.
Ijerph 16 04522 g008
Ijerph 16 04522 g009 550
Figure 9. Comparison of the proposed method with baseline approaches using the indicators of mean absolute error (MAE) and StDev.
Figure 9. Comparison of the proposed method with baseline approaches using the indicators of mean absolute error (MAE) and StDev.
Ijerph 16 04522 g009
Ijerph 16 04522 g010 550
Figure 10. Box plots of the differences in each pair of exposure estimates on a workday and a weekend.
Figure 10. Box plots of the differences in each pair of exposure estimates on a workday and a weekend.
Ijerph 16 04522 g010
Ijerph 16 04522 g011 550
Figure 11. Illustration of individual exposure-level trajectories.
Figure 11. Illustration of individual exposure-level trajectories.
Ijerph 16 04522 g011
Ijerph 16 04522 g012 550
Figure 12. Residential locations with high exposure individuals.
Figure 12. Residential locations with high exposure individuals.
Ijerph 16 04522 g012
Ijerph 16 04522 g013 550
Figure 13. Relationship between the percentage of residents’ time away from residence and the corresponding percentage of exposure levels.
Figure 13. Relationship between the percentage of residents’ time away from residence and the corresponding percentage of exposure levels.
Ijerph 16 04522 g013
Table
Table 1. Instance of one individual’s trajectory data.
Table 1. Instance of one individual’s trajectory data.
User IDDateTime (t)Longitude (x)Latitude (y)Event Type
1EF53 *****102:14:25121.13 **31.06 **Regular update
1EF53 *****108:15:11121.13 **31.02 **Call (inbound)
1EF53 *****109:17:12121.12 **31.02 **Cellular handover
1EF53 *****…………
1EF53 *****721:13:06121.44 **31.08 **Call (outbound)
Note: Accurate coordinate information and user ID were hidden with ** for privacy concern.
Table
Table 2. Instance of ground-station PM2.5 Concentration Data.
Table 2. Instance of ground-station PM2.5 Concentration Data.
Station IDDayTime (t)Longitude (x)Latitude (y)PM2.5 Concentration
(μm/m3)
1144A100:00121.41 **31.16 **43
1144A101:00121.41 **31.16 **49
1144A102:00121.41 **31.16 **52
……………
1150A723:00121.57 **31.20 **20
Note: Accurate coordinate information and user ID were hidden with ** for privacy concern.
Table
Table 3. Instance of ground-station meteorological Data.
Table 3. Instance of ground-station meteorological Data.
Station IDDayTime (t)Longitude (x)Latitude (y)Wind Speed
(m/s)		Horizontal
Visibility (m)		Air Temperature (°C)
58012100:00116.65 **34.66 **1.5200−0.5
58012101:00116.65 **34.66 **1.5300−0.5
58012102:00116.65 **34.66 **1.7200−0.4
……………
58752723:00120.65 **27.78 **1.745008.8
Note: Accurate coordinate information and user ID were hidden with ** for privacy concern.
Table
Table 4. K-S test results.
Table 4. K-S test results.
Day TypeEstimate PairsK-S Statisticsp-Value
WorkdayTR-EE & REC-EE0.0039p < 0.0001
TR-EE & SL-EE0.0214p < 0.0001
WeekendTR-EE & REC-EE0.0036p = 0.0005
TR-EE & SL-EE0.0233p < 0.0001
Table
Table 5. PM2.5 concentrations and health implications.
Table 5. PM2.5 concentrations and health implications.
CategoryPM2.5Health Implications
Excellent<35Without health implications.
Good35–70Outdoor activities normally.
Lightly Polluted70–115Slight irritations for healthy people and slightly impact on sensitive individuals.
Moderately
Polluted		115–150Serious conditions for sensitive individuals. The hearts and respiratory systems of healthy people may be affected.
Severely Polluted>150Significant impact on sensitive individuals. Healthy people will commonly show symptoms.
Table
Table 6. Details of the exposure risk percentage of residents in the top five subdistricts.
Table 6. Details of the exposure risk percentage of residents in the top five subdistricts.
SubdistrictExcellentGoodLightly
Polluted		Moderately
Polluted		Severely
Polluted
Anting County0.9446.0339.438.565.04
Jiangqiao County2.4546.9538.337.304.97
Xiayang Subdistrict12.1329.7246.656.035.46
Huacao County2.4446.4538.757.295.07
Fangsong Subdistrict13.1530.5345.974.765.58

Share and Cite

MDPI and ACS Style

Li, M.; Gao, S.; Lu, F.; Tong, H.; Zhang, H. Dynamic Estimation of Individual Exposure Levels to Air Pollution Using Trajectories Reconstructed from Mobile Phone Data. Int. J. Environ. Res. Public Health 2019, 16, 4522. https://doi.org/10.3390/ijerph16224522

AMA Style

Li M, Gao S, Lu F, Tong H, Zhang H. Dynamic Estimation of Individual Exposure Levels to Air Pollution Using Trajectories Reconstructed from Mobile Phone Data. International Journal of Environmental Research and Public Health. 2019; 16(22):4522. https://doi.org/10.3390/ijerph16224522

Chicago/Turabian Style

Li, Mingxiao, Song Gao, Feng Lu, Huan Tong, and Hengcai Zhang. 2019. "Dynamic Estimation of Individual Exposure Levels to Air Pollution Using Trajectories Reconstructed from Mobile Phone Data" International Journal of Environmental Research and Public Health 16, no. 22: 4522. https://doi.org/10.3390/ijerph16224522

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop