key: cord-0032654-er0ryfy4 authors: Wang, Huizi; Luo, Xiao; Liu, Chao; Fu, Qingyan; Yi, Min title: Spatio-Temporal Variation-Induced Group Disparity of Intra-Urban NO(2) Exposure date: 2022-05-12 journal: Int J Environ Res Public Health DOI: 10.3390/ijerph19105872 sha: 1bac9935b48065d610b0f872b01a2efbfe5e0a3d doc_id: 32654 cord_uid: er0ryfy4 Previous studies on exposure disparity have focused more on spatial variation but ignored the temporal variation of air pollution; thus, it is necessary to explore group disparity in terms of spatio-temporal variation to assist policy-making regarding public health. This study employed the dynamic land use regression (LUR) model and mobile phone signal data to illustrate the variation features of group disparity in Shanghai. The results showed that NO(2) exposure followed a bimodal, diurnal variation pattern and remained at a high level on weekdays but decreased on weekends. The most critical at-risk areas were within the central city in areas with a high population density. Moreover, women and the elderly proved to be more exposed to NO(2) pollution in Shanghai. Furthermore, the results of this study showed that it is vital to focus on land-use planning, transportation improvement programs, and population agglomeration to attenuate exposure inequality. Air pollution has become one of the most serious environmental problems globally. Worldwide, approximately 4.2 million premature deaths could be attributed to air pollution in 2016 [1] . In China, 112.7 of every 100,000 deaths were caused by air pollutants [2] . The rapid development of heavy chemical industries, energy consumption, and motor vehicle ownership has resulted in massive emissions of nitrogen dioxide (NO 2 ), which have led to increasing pollution by secondary pollutants such as acid rain and photochemical pollution. Epidemiological investigations have confirmed that NO 2 pollution is associated with considerable health risks, even below the current WHO air quality guidelines [3] . Shortterm exposure to NO 2 pollution may cause airway responsiveness and lung function injury, while long-term exposure may impair the immune function, even potentially increasing the risk and severity of respiratory viral infection [4] [5] [6] . Furthermore, mounting evidence indicates that vulnerable social groups are more likely to be exposed to air pollution [7, 8] . However, most studies on exposure inequality focused on the annual average data, ignoring the impacts of temporal variation of air pollution. Therefore, it is necessary to conduct studies on the group disparity in population exposure assessment in terms of spatiotemporal variation to assist policy-making aimed at attenuating exposure inequality and improving public health. Exposure is defined by both air pollution concentration in places where people spend time and the amount of time spent in each place [9, 10] . It relies on both air pollutant concentration and population distribution. The land use regression (LUR) model is recognized as a common method to simulate spatial-temporal variations in air pollution. Most service (LBS) data [66] [67] [68] to analyze dynamic exposure levels. These approaches provide more detailed spatio-temporal analyses of individual travel behavior, but, for researchers and participants, these methods are often inefficient, expensive, and limited in terms of how many people they can track. Today, increasingly available mobile phone signal data provide a new opportunity to further improve measurements of population exposure. However, to date, limited efforts have been made to use large-scale mobile phone signal data to examine variations in exposure to air pollution [51, 69] . The first study that used mobile phone data to examine the disparity in exposure to air pollution in China was performed by Guo et al. in 2020 [70] . This study utilized data collected on a weekday (23 March 2012) in Shenzhen. Given the limited time frame, medium-and short-term population exposure remains unclear. Furthermore, this study inappropriately applied time-profile population exposure to assess individual-level exposure over a large spatial area instead of the microenvironment such as an individual's residence or workplace [71] . Environmental justice believes that regardless of socioeconomic status, all residents ought to enjoy the equal benefits of public natural resources and shoulder equal adverse health impacts from deteriorating environmental conditions [72] [73] [74] . Increasingly, the literature has been illustrating that several subpopulations may have born a disproportionate air pollution exposure and health burden over the past several decades [75] . This is of particular concern because a better understanding of the impacts of air pollution exposure on environmental justice is a public health policy-making priority. However, there exists a gap in China due to the lack of data and neglect for environmental justice issues. Many environmental justice studies have noted the inequitable distribution of air pollution exposure among marginal social groups (e.g., poor people, Blacks, and children) [7, 8] , in part because these vulnerable populations are more likely to reside near air pollutant sources where cheaper housing and more job opportunities congregate. For example, a national study in the United States has revealed that low-income non-White children and the elderly were more likely to be exposed to NO 2 [76] . A similar conclusion was drawn in Sweden, where children living in poorer conditions were burdened with higher NO 2 exposure [77] . Race or ethnicity is a particularly consequential element when it comes to air pollution exposure disparity. Many studies in the USA have demonstrated that Blacks and Latinos are disproportionately exposed to higher levels of NO 2 than whites [46, [78] [79] [80] . Evidence from nine European metropolitan areas, indicates that higher NO 2 values were observed in areas with higher populations of individuals born outside the European Union [81] . However, some cases contradict the above studies; for example, Toronto reported that racial minority groups tended to have less air pollution exposure, probably due to the immigration policy aimed at highly educated groups [82] . Low-income populations also suffer more adverse impacts of air pollution, even they have contributed little to the pollution [83] [84] [85] . In addition, recent environmental justice research attempts to understand the evolution of environmental inequality through long-term dynamic analysis. These studies found that vulnerable populations benefit the least from air-quality improvement; in other words, exposure decreased more in less polluted areas, which means that group disparity in air pollution increased further [86] [87] [88] . In this section, the study area and research data we used to model aggregate-level NO 2 exposure are introduced in detail, followed by the modeling and validation of the LUR model and the calculation and comparison of population exposure. The research flow chart is shown in Figure 1 . LUR model and the calculation and comparison of population exposure. The research flow chart is shown in Figure 1 . Shanghai is one of the municipalities under direct administration of the Central Government of China and at the core of the world-class city cluster in the Yangtze River Delta area. It is composed of 16 districts with an area of 6340.5 km 2 and a total population of 24.28 million as of the end of 2019. There were 56 days when air pollution was severe during 2019 [89] , and the primary pollutants were fine particles and NO2. The mean value of NO2 concentration was 42 μg/m 3 in Shanghai, 2 μg/m 3 higher than the standard set by the WHO (40 μg/m 3 ) [90] . In general, severe air pollution still hangs over Shanghai. It should be mentioned that the study area in this study includes all districts except the isolated Chongming Island due to the differences in the geographical location and climate conditions from other districts as shown in Figure 2 . Shanghai is one of the municipalities under direct administration of the Central Government of China and at the core of the world-class city cluster in the Yangtze River Delta area. It is composed of 16 districts with an area of 6340.5 km 2 and a total population of 24.28 million as of the end of 2019. There were 56 days when air pollution was severe during 2019 [89] , and the primary pollutants were fine particles and NO 2 . The mean value of NO 2 concentration was 42 µg/m 3 in Shanghai, 2 µg/m 3 higher than the standard set by the WHO (40 µg/m 3 ) [90] . In general, severe air pollution still hangs over Shanghai. It should be mentioned that the study area in this study includes all districts except the isolated Chongming Island due to the differences in the geographical location and climate conditions from other districts as shown in Figure 2 . Shanghai is one of the municipalities under direct administration of the Central Government of China and at the core of the world-class city cluster in the Yangtze River Delta area. It is composed of 16 districts with an area of 6340.5 km 2 and a total population of 24.28 million as of the end of 2019. There were 56 days when air pollution was severe during 2019 [89] , and the primary pollutants were fine particles and NO2. The mean value of NO2 concentration was 42 μg/m 3 in Shanghai, 2 μg/m 3 higher than the standard set by the WHO (40 μg/m 3 ) [90] . In general, severe air pollution still hangs over Shanghai. It should be mentioned that the study area in this study includes all districts except the isolated Chongming Island due to the differences in the geographical location and climate conditions from other districts as shown in Figure 2 . A spatio-temporal LUR model was developed by utilizing hourly dynamic variables to establish several models in different periods. Considering that road transport is one of the primary contributors to the ambient concentration of NO 2 , a day was divided into four periods, namely, morning peak (7 a.m.-9 a.m.), daytime (9 a.m.-5 p.m.), evening peak (5 p.m.-7 p.m.), and nighttime (7 p.m.-7 a.m.), regarding traffic states in Shanghai [91] . The mean values of the hourly concentrations of NO 2 during different periods for 50 monitors in the 16 districts were retrieved from the Shanghai Environmental Monitoring Center and cover a time span from 11 November to 30 November 2019. The independent variables included location, meteorological elements, road network, land use, point of interest, and other pollutants as listed in Table 1 . It should be mentioned that the suburban and urban areas were set as binary variables, with 1 referring to urban monitoring stations and 0 to suburban monitoring stations, for which the boundary was the outer-ring highway in Shanghai. Mean hourly data for meteorological elements were extracted from Inverse distance weight (IDW) interpolation pollution maps which were drawn in ArcGIS based on the meteorological data from China Meteorological Data Service Center. The road network was elaborated as the intensity of highways and local roads, whose maximum buffer distance of the road networks was set at 1000 m, following the principle proposed by Hock [92] . Land-use data contained five types, namely, residential land use, commercial land use, industrial land use, green space, and water body, which were calculated in each buffer zone from monitor stations. Four types of points of interest (POIs) were extracted using the Baidu Open Map Platform (2019) based on categories and keywords. The numbers of restaurants, bus stops, intersections, and gas stations were calculated within each buffer. µg/m 3 − Note: "+" denoted assumed positive correlation. "−" denoted assumed negative correlation, "O" denoted not assigned because no effect could be assumed. More details were described in Section 3.2.3. The explanatory variables were chosen by a supervised forward regression method [93, 94] . First, each potential explanatory variable was assigned to a prior direction to enhance the applicability of the LUR models described in Table 1 . For example, wind speed could improve the dilution and diffusion process of pollutants in the atmosphere; thus, its prior direction was assumed as negative (denoted as "−"). A stable atmospheric configuration under a high-pressure environment was not conducive to the physical diffusion of pollutants, leading to the accumulation of pollutant concentrations. Therefore, its prior direction was assigned as positive (denoted as "+"). As for relative humidity, some studies have shown that higher relative humidity is detrimental to the diffusion process of gaseous pollutants, worsening environmental pollution, but other studies came to the opposite conclusion that higher relative humidity improves the transition of gaseous pollutants to particulate state, reducing the pollution concentration. Consequently, the prior direction of relative humidity was not specified since it's still a controversial correlation and denoted by "O". Secondly, univariate regression was carried out between the dependent variables and each independent variable, and the resulting F-statistic and corresponding p-value were used to determine the significance of the predictor and the order for its entry into the model. The model with the highest adjusted R 2 that was consistent with the prior direction was selected as the initial model of stepwise regression. Thereafter, a supervised forward regression was conducted in IBM SPSS (version 26.0), and the variables were introduced into the model if they satisfied the following criteria: (1) the adjusted R 2 of the model increased by at least 1%; (2) the coefficient of each variable was consistent with the prior direction; (3) the existing variables in the model did not change their effect directions. All potential explanatory variables were introduced into the model one by one until there were no remaining variables that satisfied the above criteria. Then, the variables with p-values greater than 0.1 were excluded from the model. Thirdly, standardized diagnostic tests were applied to the final models to check the multicollinearity between the variables and influential observations. The variables with the highest variance inflation factors (VIFs) of more than 3 were excluded from the final model and the model was recalibrated. Finally, residual analysis and cross-validation were applied to evaluate the performance of the model. The former plays an important role in validating the regression model, consisting of a test for normality, test for equal variance, test for independence of residuals, and test for spatial autocorrelation, which were developed by ArcGIS and SPSS packages. If the error term satisfies the four basic assumptions of the regression model, then the model is considered valid. The latter was used to assess the accuracy of a linear regression model. In this paper, the leave-one-out cross-validation (LOOCV) method was used to evaluate the model's performance. Each site was subsequently left out, and a model was developed from the remaining sites with the variables unchanged. Predicted concentrations were estimated and compared with the actual concentrations at each left-out site. The LOOCV R 2 and root mean squared error (RMSE) between the predicted and measured concentrations were calculated for all monitoring stations to represent the model performance. The mobile phone signal dataset in Shanghai, including the time, location, and users' information for the period from 11 November to 30 November 2019, was provided by JI SMART (http://daas.smartsteps.com/, accessed on 19 January 2021). Initially, the number of mobile phone users included in the data was approximately 5.28 million, and this number can be used to represent Shanghai's population of approximately 23.55 million, excluding the Chongming district. The main detailed information in the dataset included (1) user attributes, which contained an anonymously processed user ID, gender, age, etc.; (2) date, start time, and end time of stay points; (3) location, which was represented by a grid number with a resolution of 250 m instead of longitude and latitude; (4) grid information containing grid ID, coverage area (well-known text), etc. The dataset we acquired for the study was processed. More details on the data processing, including the identification of the homes and workplaces of users, were not available due to the confidentiality agreement. Briefly, the locations of each user's home and workplace were identified according to how long they remained at certain locations. Home locations were recognized based on the location where users remained the longest between 9 p.m. and 8 a.m. during the month. Similarly, users' workplaces were identified according to the location where they spent the most time between 9 a.m. to 5 p.m. Stay points, excluding the home and workplace, were marked as other activity locations. A correlation test and paired samples t-tests were used to verify whether there was a significant difference between the seventh population census data and mobile phone signal data, regarding the population distribution, sex composition, and age structure in the districts of Shanghai. It is worth mentioning that population proportion was chosen for comparison rather than the population size due to the different magnitudes of the two sets of data. First, population distribution was defined as the proportion of each district's residential population to the total population of Shanghai, designed to indicate whether the two sets of data were consistent in the geographical distribution of the population. Pearson's correlation coefficient of population distribution was 0.998, close to 1, indicating that the two samples had a strong linear correlation. The p-value from the paired-samples t-test was 0.999, which was greater than the significance level of α = 0.05; that is, there was no sufficient evidence to claim that the population distribution from the mobile phone signal dataset differs from the one from the Seventh National Census data as shown in Table 2 . Secondly, the sex ratio defined as the number of females per 1000 males was chosen as a social indicator of sex composition as listed in Table 3 . The correlation between the two samples was considered to be strong, because the absolute value of the Pearson's correlation coefficient was 0.794, greater than 0.75. In addition, the paired sample t-test (p-value = 0.374) confirmed that there was no statistically significant difference between the seventh population census data and mobile phone signal data. That is to say, the mobile phone signal data could accurately capture the population's gender structure in the districts of Shanghai. Thirdly, the proportions of the population aged 0-14, 15-64, and over 65 years were used to describe the population age structure in the districts of Shanghai as listed in Table 4 . The correlation test and paired samples t-tests were performed separately for each age group, and the Pearson's correlation coefficients were 0.518, 0.912, and 0.942, respectively, indicating strong positive linear correlations between the three groups of paired samples. Furthermore, the p-values of the t-tests were 0.997, 0.627, and 0.499, which meant that there was no significant difference in the age structure, suggesting that the mobile phone signal dataset could denote the age structure of different areas in Shanghai. In conclusion, it was proved that the mobile phone signal dataset used in this study was a representative sample of the population distribution, sex composition, and age structure in Shanghai. Population-weighted exposure level (PWEL), proposed by Fu and Kan [95] , was chosen as the exposure assessment indicator in this paper. It was calculated as predicted concentrations combined with population dynamic distributions. First, NO 2 pollution during different periods was estimated with the LUR models. Secondly, hourly gridded population data with a spatial resolution of 500 × 500 m was derived from the mobile phone data. Then, using ArcGIS software, the NO 2 concentration layers were overlaid on the gridded population distribution layers and NO 2 exposure level in each grid was calculated based on Equation (1) as follows: where i indicates grids, n indicates the number of grids, E i represents grid i's potential population exposure, P i denotes grid i's population size in a certain period, and C i denotes grid i's concentration of specific air pollutants. Moreover, the mean exposure was used as a baseline to classify each grid's exposure into one of four levels, namely, low risk (<0.5 SD), medium risk (>0.5 and ≤1.5 SD), high risk (>1.5 and ≤2.5 SD), and critical risk (>2.5 SD). Finally, the population-weighted exposure level of NO 2 for Shanghai was assessed using Equation (2), where E represents Shanghai's potential NO 2 exposure: In this section, NO 2 exposure characteristics in Shanghai were explored for spatiotemporal variations and group disparity (gender and age), which is crucial for assisting policy-making aimed at attenuating health disparities. Firstly, the diurnal and weekly variations in population exposure were identified based on weekly and daily averages of NO 2 exposure for the period from 11 to 30 November 2019. Next, NO 2 exposure levels for each grid area were analyzed to determine spatial variations and to highlight risky areas that required additional attention. Finally, NO 2 exposure was compared according to gender and age to quantify exposure disparities between different groups. Different influencing factors finally entered the LUR models in different periods. Overall, land use was the most significant indicator, where green spaces and water bodies proved to effectively absorb and purify air pollutants, and then meteorological elements and traffic-related variables as listed in Tables 5-8 . Approximately 56-73% of the variations in NO 2 pollution were explained by the final models and the RMSEs were 4.554, 4.732, 5.371, and 4.894, respectively, which indicates good performance. The differences between the LOOCV R 2 and the model R 2 were less than 10%, indicating the stability of the LUR model [96] . Moreover, given the z-scores of −0.757, −1.499, −0.842, and −1.295, respectively, the model residuals did not appear to be significantly different from random, fulfilling the regression model assumptions of spatial independence in residuals. After completing the final models, regression equations were applied to a regular 500 × 500 m grid covering the entire study area. A pollution point map was developed using the predicted values. As shown in Figure 3 , an evident spatio-temporal heterogeneity of NO 2 pollution was observed in Shanghai. Overall spatial variation showed a tendency to decline gradually from downtown in all directions, but the spatial pattern of different periods also presented a uniqueness. Higher NO 2 values were observed to be transferred from the urban districts with the intensive road network to suburban districts with clustered logistics and industrial parks, such as the Baoshan and Minhang Districts, due to the traffic restriction policy aimed at trucks in Shanghai that encourages freight transport at night. Peak hours shared a similar distribution pattern, but transportation contributed more to the evening peak. In general, NO 2 concentrations were higher west of than the east of the Huangpu River. In terms of temporal variation, the average values during the daytime, nighttime, morning peak, and evening peak were 29.71, 41.54, 44.07, and 46.52 µg/m 3 , respectively, which exhibited a bimodal diurnal variation, and the figure of the evening peak was higher than the early peak, approximately resembling daily diurnal traffic patterns. NO2 pollution was observed in Shanghai. Overall spatial variation showed a tendency t decline gradually from downtown in all directions, but the spatial pattern of different pe riods also presented a uniqueness. Higher NO2 values were observed to be transferre from the urban districts with the intensive road network to suburban districts with clus tered logistics and industrial parks, such as the Baoshan and Minhang Districts, due to th traffic restriction policy aimed at trucks in Shanghai that encourages freight transport a night. Peak hours shared a similar distribution pattern, but transportation contribute more to the evening peak. In general, NO2 concentrations were higher west of than th east of the Huangpu River. In terms of temporal variation, the average values during th daytime, nighttime, morning peak, and evening peak were 29.71, 41.54, 44.07, and 46.5 μg/m 3 , respectively, which exhibited a bimodal diurnal variation, and the figure of th evening peak was higher than the early peak, approximately resembling daily diurna traffic patterns. In conclusion, the results of the LUR models suggest that the pollution concentratio could be a result in a dynamic game between the air-pollutant emission intensity and en vironmental self-purification capacity. Consequently, air quality improvements not onl require an emphasis on controlling the air-pollutant sources, such as advocating low-car bon travel, prioritizing public transportation development, and advancing emissio In conclusion, the results of the LUR models suggest that the pollution concentration could be a result in a dynamic game between the air-pollutant emission intensity and environmental self-purification capacity. Consequently, air quality improvements not only require an emphasis on controlling the air-pollutant sources, such as advocating low-carbon travel, prioritizing public transportation development, and advancing emission standards, but also on taking advantage of environmental self-purification capacity. For example, it is necessary to take meteorological elements (wind speed, wind direction, precipitation, etc.) into consideration in the selection of locations for logistics and industrial parks or the height of a chimney. In addition, for air quality improvement, it would be wise to increase the rate of green land use, considering the purification capacity of plants. As previously stated, NO 2 exposure was evaluated by population-weighted exposure levels based on the pollution simulation results from the LUR models and dynamic population distribution from the mobile phone signal data. The mean values, minimum values, maximum values, and standard deviations are summarized in Table 9 , presenting bimodal variations consistent with the differences in the pollutant concentration within a day. The daily mean values of NO 2 exposure are shown in Figure 4 , where the gray shaded areas represent weekends. NO 2 exposure remained at a high level on weekdays but decreased on weekends. This fluctuation could be attributed to the weekly variations in traffic volumes; that is, there was more traffic on weekdays than on weekends due to the fact of commuters heading to work. However, this trend was not distinct during the nighttime. This was mainly because freight transport did not have an obvious weekly variation. In conclusion, NO 2 exposure might be influenced by both pollutant concentration and population dynamics. Table 9 , presenting bimodal variations consistent with the differences in the pollutant concentration within a day. The daily mean values of NO2 exposure are shown in Figure 4 , where the gray shaded areas represent weekends. NO2 exposure remained at a high level on weekdays but decreased on weekends. This fluctuation could be attributed to the weekly variations in traffic volumes; that is, there was more traffic on weekdays than on weekends due to the fact of commuters heading to work. However, this trend was not distinct during the nighttime. This was mainly because freight transport did not have an obvious weekly variation. In conclusion, NO2 exposure might be influenced by both pollutant concentration and population dynamics. In this section, the standard deviation classification method was used to place the NO2 exposure in each grid area at four levels to highlight the differences between the dis- In this section, the standard deviation classification method was used to place the NO 2 exposure in each grid area at four levels to highlight the differences between the districts. It should be mentioned that there was no comparison between different periods, because this classification was based on the mean value of all air pollution exposure during a single corresponding period. Higher risk areas were observed in the downtown area with a dense population, as shown in Figure 5 , mainly west of the Huangpu River, including the Hongkou, Yangpu, Huangpu, Xuhui, and Jing'an districts. In other words, the areas to the east of the Huangpu River were healthier with lower air pollutant exposure. Furthermore, there was a difference in the range of NO 2 exposure levels, which were more clustered during the morning peak, indicating exposure inequality was lower during that period. As previously stated, heavy traffic contributed the most to NO 2 exposure during rush hours. When combined with a large number of cross-regional population dynamics due to the commuting demand, NO 2 exposure was largely homogeneous across different regions. We would also like to note that NO 2 exposure inequality during the morning peak was lower, but exposure levels were higher than during other periods. This characteristic was not observed during the evening peak, in part because the flexible off-duty hours led to more scattered commuting. previously stated, heavy traffic contributed the most to NO2 exposure during rush hours. When combined with a large number of cross-regional population dynamics due to the commuting demand, NO2 exposure was largely homogeneous across different regions. We would also like to note that NO2 exposure inequality during the morning peak was lower, but exposure levels were higher than during other periods. This characteristic was not observed during the evening peak, in part because the flexible off-duty hours led to more scattered commuting. To investigate group disparity, user attribute information, including gender and age, was extracted from the mobile phone signal data. The ages were divided into four groups, namely juvenile (≤18), young , middle-aged (45) (46) (47) (48) (49) (50) (51) (52) (53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) (64) , and elderly (≥65). Females were subject to more NO2 exposure ( Figure 6 ). The exposure disparity between males and females was found to reach a peak during the day, especially during the morning peak, and decreased at night, which could be attributed to gender differences in commuting behavior. In general, females took more trips and had more complex trip chains than men, probably since they undertook a large number of non-work-related trips, such as shopping, delivering children to school, or accompanying the elderly to health centers [97] , namely, women had a higher probability of being exposed to traffic-related pollution or other sources of pollution. An inequality in family responsibilities may have evolved into inequality in terms of pollution exposure due to differences in travel behavior. To investigate group disparity, user attribute information, including gender and age, was extracted from the mobile phone signal data. The ages were divided into four groups, namely juvenile (≤18), young , middle-aged (45) (46) (47) (48) (49) (50) (51) (52) (53) (54) (55) (56) (57) (58) (59) (60) (61) (62) (63) (64) , and elderly (≥65). Females were subject to more NO 2 exposure ( Figure 6 ). The exposure disparity between males and females was found to reach a peak during the day, especially during the morning peak, and decreased at night, which could be attributed to gender differences in commuting behavior. In general, females took more trips and had more complex trip chains than men, probably since they undertook a large number of non-work-related trips, such as shopping, delivering children to school, or accompanying the elderly to health centers [97] , namely, women had a higher probability of being exposed to traffic-related pollution or other sources of pollution. An inequality in family responsibilities may have evolved into inequality in terms of pollution exposure due to differences in travel behavior. As the main bearer of daily housework, women's destinations of most non related trips were concentrated in urban districts with a completely public transp system, convenient shopping areas, abundant medical resources, and children's tional institutions. In a word, urban districts were more attractive to women [98 ever, these places were also heavily polluted areas as described above. Kernel analysis was applied in ArcGIS to explore the gender differences in travel distr based on trajectory data extracted from the mobile phone signal dataset, and the s deviation stretch and gamma correction were applied to the mapping to eliminate in trip volume between different groups and increase the contrast in the raster da highlight the range of hot spots. As shown in Figure 7 , the distribution of different groups in the road network was both generally of unbalanced and aggregated ch istics, and the travel distribution of women was more clustered. Male trip volume a continuous trend of decreasing outward from the urban areas but remained some suburban trunk roads. However, there was a very significant cliff-like declin trip volume of females in the urban periphery, with fewer trips on major arterial r the suburban areas compared to men. It was worth noting that the central area of where women gather is precisely the most polluted area, which may lead to the higher exposure of women than men. As the main bearer of daily housework, women's destinations of most non-workrelated trips were concentrated in urban districts with a completely public transportation system, convenient shopping areas, abundant medical resources, and children's educational institutions. In a word, urban districts were more attractive to women [98] . However, these places were also heavily polluted areas as described above. Kernel density analysis was applied in ArcGIS to explore the gender differences in travel distribution based on trajectory data extracted from the mobile phone signal dataset, and the standard deviation stretch and gamma correction were applied to the mapping to eliminate the gap in trip volume between different groups and increase the contrast in the raster dataset to highlight the range of hot spots. As shown in Figure 7 , the distribution of different gender groups in the road network was both generally of unbalanced and aggregated characteristics, and the travel distribution of women was more clustered. Male trip volume showed a continuous trend of decreasing outward from the urban areas but remained high on some suburban trunk roads. However, there was a very significant cliff-like decline in the trip volume of females in the urban periphery, with fewer trips on major arterial roads in the suburban areas compared to men. It was worth noting that the central area of the city where women gather is precisely the most polluted area, which may lead to the overall higher exposure of women than men. As the main bearer of daily housework, women's destinations of most non-workrelated trips were concentrated in urban districts with a completely public transportation system, convenient shopping areas, abundant medical resources, and children's educational institutions. In a word, urban districts were more attractive to women [98] . However, these places were also heavily polluted areas as described above. Kernel density analysis was applied in ArcGIS to explore the gender differences in travel distribution based on trajectory data extracted from the mobile phone signal dataset, and the standard deviation stretch and gamma correction were applied to the mapping to eliminate the gap in trip volume between different groups and increase the contrast in the raster dataset to highlight the range of hot spots. As shown in Figure 7 , the distribution of different gender groups in the road network was both generally of unbalanced and aggregated characteristics, and the travel distribution of women was more clustered. Male trip volume showed a continuous trend of decreasing outward from the urban areas but remained high on some suburban trunk roads. However, there was a very significant cliff-like decline in the trip volume of females in the urban periphery, with fewer trips on major arterial roads in the suburban areas compared to men. It was worth noting that the central area of the city where women gather is precisely the most polluted area, which may lead to the overall higher exposure of women than men. Regarding the age groups, NO2 exposure increased significantly in the elderly as shown in Figure 8 ; that is, the elderly group endured the greatest exposure, followed by the young and the middle-aged, and then the juvenile groups. Regarding the age groups, NO 2 exposure increased significantly in the elderly as shown in Figure 8 ; that is, the elderly group endured the greatest exposure, followed by the young and the middle-aged, and then the juvenile groups. Residential preference and travel behaviors may have contributed to the observed disparity between age groups. Since the 1990s, the policy of suppressing the second industry and developing the third industry has prompted factories in the central area of the city to relocate to the suburbs and gather to form industrial parks. At the same time, the planning and construction of new towns were in full swing under the influence of the Greater Shanghai Plan. The working-age population in the central urban area and from surrounding cities migrated to the suburbs or new towns due to a large number of employment opportunities and moderately priced housing, which contributed to the age structure differentiation between the urban area and the suburban areas of Shanghai. Finally, the elderly population was concentrated in the urban areas while young and middle-aged groups gathered near the suburbs. According to the seventh Shanghai Census Bulletin, Hongkou and Huangpu districts have the densest elderly population, with a density of more than 5000 people per square kilometer. The elderly population in the suburbs is relatively sparse, such as in Jinshan, Qingpu, and Fengxian Districts, where the density of the elderly population was less than 250 people per square kilometer. Among them, the difference between the densest and the sparsest regions was as high as 35 times. In general, the density of the resident elderly population in Shanghai was low around the middle and uneven in the regional distribution. Furthermore, the characteristics of the circle structure were significant. As listed in Table 10 , the population proportion over 60 years in the urban districts of Shanghai exceeded 25%, with an average of 30.2%. However, suburbs were mostly below 20%, with an average of just 19.6%. To this end, the global Moran's I index was used in ArcGIS to further evaluate the spatial autocorrelation of the people over the age of 65 in the districts of Shanghai. Given the z-score of 3.02, much larger than the critical value of 1.65, there was more than a 95% likelihood that Shanghai's elderly population could be the result of the clustered pattern. It meant that the spatial distribution of the elderly population has strong homogeneity, that is, it is uneven in the regional distribution. Residential preference and travel behaviors may have contributed to the observed disparity between age groups. Since the 1990s, the policy of suppressing the second industry and developing the third industry has prompted factories in the central area of the city to relocate to the suburbs and gather to form industrial parks. At the same time, the planning and construction of new towns were in full swing under the influence of the Greater Shanghai Plan. The working-age population in the central urban area and from surrounding cities migrated to the suburbs or new towns due to a large number of employment opportunities and moderately priced housing, which contributed to the age structure differentiation between the urban area and the suburban areas of Shanghai. Finally, the elderly population was concentrated in the urban areas while young and middle-aged groups gathered near the suburbs. According to the seventh Shanghai Census Bulletin, Hongkou and Huangpu districts have the densest elderly population, with a density of more than 5000 people per square kilometer. The elderly population in the suburbs is relatively sparse, such as in Jinshan, Qingpu, and Fengxian Districts, where the density of the elderly population was less than 250 people per square kilometer. Among them, the difference between the densest and the sparsest regions was as high as 35 times. In general, the density of the resident elderly population in Shanghai was low around the middle and uneven in the regional distribution. Furthermore, the characteristics of the circle structure were significant. As listed in Table 10 , the population proportion over 60 years in the urban districts of Shanghai exceeded 25%, with an average of 30.2%. However, suburbs were mostly below 20%, with an average of just 19.6%. To this end, the global Moran's I index was used in ArcGIS to further evaluate the spatial autocorrelation of the people over the age of 65 in the districts of Shanghai. Given the z-score of 3.02, much larger than the critical value of 1.65, there was more than a 95% likelihood that Shanghai's elderly population could be the result of the clustered pattern. It meant that the spatial distribution of the elderly population has strong homogeneity, that is, it is uneven in the regional distribution. Based on the residence distribution of different age groups extracted from the mobile phone signal dataset, the spatial distribution of the population density between age groups was described through the mean-standard deviation method in ArcGIS. As shown in Figure 9 , the spatial distribution of the elderly population had a more significant clustering characteristic than other age groups. It is worth noting that the urban districts, where the elderly population was densely distributed, were the areas with the most critical NO 2 pollution. This means that the elderly population with the lowest pollution contribution is bearing the highest risk of exposure; that is, Shanghai has serious environmental inequities in terms of air quality. Based on the residence distribution of different age groups extracted from the mobile phone signal dataset, the spatial distribution of the population density between age groups was described through the mean-standard deviation method in ArcGIS. As shown in Figure 9 , the spatial distribution of the elderly population had a more significant clustering characteristic than other age groups. It is worth noting that the urban districts, where the elderly population was densely distributed, were the areas with the most critical NO2 pollution. This means that the elderly population with the lowest pollution contribution is bearing the highest risk of exposure; that is, Shanghai has serious environmental inequities in terms of air quality. In summary, this study attempted to combine dynamic LUR models with mobile phone signal data to explain the spatio-temporal variations and group disparities of NO 2 exposure in Shanghai. Our main conclusions can be summarized as follows. The dynamic LUR models revealed an evident temporal variation of NO 2 exposure in Shanghai. Overall, it followed a bimodal diurnal variation and remains at a high level on weekdays but decreases on weekends, consistent with changes in traffic volume; 2. In terms of spatial variation, higher-risk locations included urban areas with dense populations and busy traffic and were concentrated west of the Huangpu River. Also, lower regional inequality of NO 2 exposure was observed during the morning peak due to a large number of cross-regional commutes that led to NO 2 exposure homogenization across different regions; 3. As for group disparity, women and the elderly proved to suffer more exposure to NO 2 pollution, which could be attributed to gender differences in travel behavior and the preference of residence in different age groups. This study aimed to introduce an improved NO 2 exposure assessment by integrating the dynamic LUR model with mobile phone signal data and explored the spatio-temporal variation of the population exposure disparity in the case of Shanghai. Nevertheless, some limitations in this study should be pointed out. First, limited by the availability of pollution concentration and mobile phone signal data, only diurnal and weekly variations in NO 2 exposure were explored. In the future, the data could be strengthened and validated on a monthly, seasonal, or annual scale over a longer period. Secondly, this paper attempted to establish a link between traffic-related air pollution and human risk. However, the lack of quantitative analysis of traffic emissions made it difficult to strongly support this conclusion. Thirdly, although this paper has revealed the exposure disparity in terms of spatio-temporal variation and age/gender group, there is still a long way to go in addressing the issues of environmental and health equity. The subsequent research is supposed to further dig out whether the implementation of financial support policies such as subsidies and tax cuts could attenuate exposure inequality and improve public health. These results suggest that Shanghai is trapped in an awkward predicament with an uncoordinated relationship between urban development and air quality. Currently, it is necessary to create a comprehensive environmental plan and a target-and results-oriented governance program to improve air quality and public health. Based on our analysis, we propose several recommendations. First, air pollution improvement requires regulators to pay close attention to pollution sources. Scientific land-use planning can be an ideal tool to coordinate urban economic development and quality living environments. For instance, it would be prudent to increase the number of green spaces near the pollutant source to prevent air pollutants from spreading elsewhere. In addition, the needs and interests of vulnerable social groups have to be prioritized in the process of the formulation and implementation of urban planning and environmental health policies to ensure environmental justice and sustainable development. Second, transportation improvement programs are also considered an effective way to alleviate the negative impacts of motorization. Today, TOD (transit-oriented development) has become a priority for improving public health and developing a quality living environment, as it encourages enhanced accessibility with a superior and sustainable transit system along with mixed land use and compact urban development. Moreover, the boom in ridesharing and electric vehicles has also created more opportunities to reduce the contribution of the road networks to air pollution. Third, higher risk areas for air pollution exposure are generally downtown in locations with high population density, and population agglomeration aggravates air pollution, as revealed in previous studies [99] . Therefore, the Shanghai municipal government should develop several effective control strategies for population agglomeration to guide the reasonable spatial distribution of the population to reduce the intensity of population exposure inequality. In summary, this study proposed a fairly novel approach that integrated a dynamic LUR model and mobile phone signal data to improve the accuracy of dynamic population exposure. A more detailed analysis of group disparity was investigated in terms of spatiotemporal variation in the case of Shanghai. Although some limitations exist, our study could help researchers better understand the spatio-temporal patterns of NO 2 pollution exposure and assist policy-making aimed at improving public health. The results showed that it is vital to focus on land-use planning, transportation improvement programs, and population agglomeration to attenuate exposure inequality. The data presented in this study are available on request from the corresponding author. The data are not publicly available due to the privacy restrictions. Ambient (Outdoor) Air Quality and Health; World Health Organization Sustainable Development Goals; World Health Organization Short term associations of ambient nitrogen dioxide with daily total, cardiovascular, and respiratory mortality: Multilocation analysis in 398 cities Long-term exposure to low ambient air pollution concentrations and mortality among 28 million people: Results from seven large European cohorts within the ELAPSE project Links between air pollution and COVID-19 in England Long-term trends in urban NO 2 concentrations and associated paediatric asthma incidence: Estimates from global datasets An environmental justice analysis of British air quality Environmental inequalities in New Zealand: A national study of air pollution and environmental justice Concepts of human exposure to air pollution Air Pollution Assessing Total Exposure in the United States Intercity transferability of land use regression models for estimating ambient concentrations of nitrogen dioxide Estimating Long-Term Average Particulate Air Pollution Concentrations: Application of Traffic Indicators and Geographic Information Systems A regression-based method for mapping traffic-related air pollution: Application and testing in four contrasting urban environments Small area variations in ambient NO 2 concentrations in four European areas Assessing spatial variability of ambient nitrogen dioxide in Montréal, Canada, with a land-use regression model Modeling the Intraurban Variability of Ambient Traffic Pollution in Toronto A land use regression model for predicting ambient concentrations of nitrogen dioxide in Hamilton, Ontario Application of Land Use Regression to Estimate Long-Term Concentrations of Traffic-Related Nitrogen Oxides and Fine Particulate Matter Use of GIS and ancillary variables to predict volatile organic compound and nitrogen dioxide levels at unmonitored locations A land use regression application into assessing spatial variation of intra-urban fine particulate matter (PM 2.5 ) and nitrogen dioxide (NO 2 ) concentrations in City of Shanghai Land use regression models coupled with meteorology to model spatial and temporal variability of NO 2 and PM10 in Changsha Development of land use regression models for PM 2.5 , SO 2 , NO 2 and O 3 in Nanjing Exposure inequality assessment for PM 2.5 and the potential association with environmental health in Beijing Modeling temporal and spatial variability of traffic-related air pollution: Hourly land use regression models for black carbon Long-term exposure to traffic-related air pollution and the risk of coronary heart disease hospitalization and mortality Predicting personal exposure of pregnant women to traffic-related air pollutants Land-use regression panel models of NO 2 concentrations in Seoul Development of temporally refined land-use regression models predicting daily household-level air pollution in a panel study of lung function among asthmatic children Land use regression modeling to estimate historic (1962-1991) concentrations of black smoke and sulfur dioxide for Great Britain A land use regression model incorporating data on industrial point source pollution Spatiotemporal Modeling of Ozone Levels in Quebec (Canada): A Comparison of Kriging, Land-Use Regression (LUR), and Combined Bayesian Maximum Entropy-LUR Approaches An hourly regression model for ultrafine particles in a near-highway urban area Temporal stability of land use regression models for traffic-related air pollution Spatial Modeling of PM 10 and NO 2 in the Continental United States Applying land use regression model to estimate spatial variation of PM2.5 in Beijing Ageing and dynamical scaling in the critical Ising spin glass Evidence on the impact of sustained exposure to air pollution on life expectancy from China's Huai River policy Early life exposure to ambient air pollution and childhood asthma in China The association between lung cancer incidence and ambient air pollution in China: A spatiotemporal analysis GIS-based analysis of population exposure to PM2.5 air pollution-A case study of Beijing Cancer risk from exposure to hazardous air pollutants: Spatial and social inequities in Tampa Bay Uneven magnitude of disparities in cancer risks from air toxics A hybrid model for evaluating exposure of the general population in Israel to air pollutants A risk-based model to assess environmental justice and coronary heart disease burden from traffic-related air pollutants Spatial variability in mortality inequalities, socioeconomic deprivation, and air pollution in small areas of the Barcelona Metropolitan Region Environmental inequality in exposures to airborne particulate matter components in the United States Risk Assessment of Exposure to PM 2.5 in Beijing Using Multi-Source Data Study on the exposure intensity of air pollution in Zhengzhou city based on nighttime lighting Health impact assessment of air pollution using a dynamic exposure profile: Implications for exposure and health impact estimates Predicting vehicular emissions in high spatial resolution using pervasively measured transportation data and microscopic emissions model Population dynamics based on mobile phone data to improve air pollution exposure assessments Dynamic assessment of exposure to air pollution using mobile phone data Improving estimates of air pollution exposure through ubiquitous sensing technologies Measurements of personal exposure to NO 2 and modelling using ambient concentrations and activity data Dynamic assessments of population exposure to urban greenspace using multi-source big data Integrating travel behavior with land use regression to estimate dynamic air pollution exposure in Hong Kong PM 2.5 population exposure in New Delhi using a probabilistic simulation framework Spatial variations in estimated chronic exposure to traffic-related air pollution in working populations: A simulation Regional assessment of exposure to traffic-related air pollution: Impacts of individual mobility and transit investment scenarios A dynamic activity-based population modelling approach to evaluate exposure to air pollution: Methods and application to a Dutch urban area Impact of time-activity patterns on personal exposure to black carbon Dynamic assessment of inhaled air pollution using GPS and accelerometer data Traffic exposure near the Los Angeles-Long Beach port complex: Using GPS-enhanced tracking to assess the implications of unreported travel and locations London hybrid exposure model: Improving human exposure estimates to NO 2 and PM 2.5 in an urban setting Unraveling environmental justice in ambient PM 2.5 exposure in Beijing: A big data approach Dynamic assessment of PM 2.5 exposure and health risk using remote sensing and geo-spatial big data Quantifying personal exposure to air pollution from smartphone-based location data Exposure Track"-The Impact of Mobile-Device-Based Mobility Patterns on Quantifying Population Exposure to Air Pollution Mobile phone tracking: In support of modelling traffic-related air pollution contribution to individual exposure and its implications for public health impact assessment Who are more exposed to PM 2.5 pollution: A mobile phone data approach A Geographic Approach to Modelling Human Exposure to Traffic Air Pollution Using GIS Observed inequality in urban greenspace exposure in China Environmental Inequality: The Social Causes and Consequences of Lead Exposure Race, class and environmental justice Socioeconomic disparities and air pollution exposure: A global review National patterns in environmental injustice and inequality: Outdoor NO 2 air pollution in the United States Children's exposure to nitrogen dioxide in Sweden: Investigating environmental injustice in an egalitarian country Health disparities attributable to air pollutant exposure in North Carolina: Influence of residential environmental and social factors The long-term dynamics of racial/ethnic inequality in neighborhood air pollution exposure Environmental inequality: Air pollution exposures in California's South Coast Air Basin Spatial variability in air pollution exposure in relation to socioeconomic indicators in nine European metropolitan areas: A study on environmental inequality Geographies of susceptibility and exposure in the city: Environmental inequity of traffic-related air pollution in Toronto City planning and population health: A global challenge Environmental justice and health: The implications of the socio-spatial distribution of multiple environmental deprivation for health inequalities in the United Kingdom Migration, environmental hazards, and health outcomes in China Who benefits from environmental policy? An environmental justice analysis of air quality change in Britain Temporal trends in air pollution exposure inequality in Massachusetts Inequality in Exposure to Air Pollution in France: Measurement and Impact of a City-Level Public Policy Air Quality Guidelines Global Update 2005: Particulate Matter, Ozone, Nitrogen Dioxide and Sulfur Dioxide Shanghai Urban and Rural Construction and Transportation Development Research Institute. Main results of the Fifth Comprehensive Transportation Survey in Shanghai A review of land-use regression models to assess spatial variation of outdoor air pollution Development of NO 2 and NO x land use regression models for estimating air pollution exposure in 36 study areas in Europe-The ESCAPE project Land use regression models for estimating individual NOx and NO 2 exposures in a metropolis with a high density of traffic roads and population Air Pollution Dispersion Model and Assessment of Population Weighted Exposure Development of land use regression models for PM 2.5 , PM 2.5 absorbance, PM 10 and PM coarse in 20 European study areas; results of the ESCAPE project Proceedings of the 2011 International Transport Forum The spatial-temporal activity pattern of women in the era of network information: A case study of Nanjing Do economic development and population agglomeration inevitably aggravate haze pollution in China? New evidence from spatial econometric analysis The authors wish to express our gratitude to all those who have helped us during the writing of this thesis. The authors declare no conflict of interest, and the funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.