Skip to main content
  • Research article
  • Open access
  • Published:

Spatiotemporal patterns of particulate matter (PM) and associations between PM and mortality in Shenzhen, China

Abstract

Background

Most studies on air pollution exposure and its associations with human health in China have focused on the heavily polluted industrial areas and/or mega-cities, and studies on cities with comparatively low air pollutant concentrations are still rare. Only a few studies have attempted to analyse particulate matter (PM) for the vibrant economic centre Shenzhen in the Pearl River Delta. So far no systematic investigation of PM spatiotemporal patterns in Shenzhen has been undertaken and the understanding of pollution exposure in urban agglomerations with comparatively low pollution is still limited.

Methods

We analyze daily and hourly particulate matter concentrations and all-cause mortality during 2013 in Shenzhen, China. Temporal patterns of PM (PM2.5 and PM10) with aerodynamic diameters of 2.5 (10) μm or less (or less (including particles with a diameter that equals to 2.5 (10) μm) are studied, along with the ratio of PM2.5 to PM10. Spatial distributions of PM10 and PM2.5 are addressed and associations of PM10 or PM2.5 and all-cause mortality are analyzed.

Results

Annual average PM10 and PM2.5 concentrations were 61.3 and 39.6 μg/m3 in 2013. PM2.5 failed to meet the Class 2 annual limit of the National Ambient Air Quality Standard. PM2.5 was the primary air pollutant, with 8.8 % of days having heavy PM2.5 pollution. The daily PM2.5/PM10 ratios were high. Hourly PM2.5 concentrations in the tourist area were lower than downtown throughout the day. PM10 and PM2.5 concentrations were higher in western parts of Shenzhen than in eastern parts. Excess risks in the number of all-cause mortality with a 10 μg/m3 increase of PM were 0.61 % (95 % confidence interval [CI]: 0.50–0.72) for PM10, and 0.69 % (95 % CI: 0.55–0.83) for PM2.5, respectively. The greatest ERs of PM10 and PM2.5 were in 2-day cumulative measures for the all-cause mortality, 2-day lag for females and the young (0–65 years), and L02 for males and the elder (>65 years). PM2.5 had higher risks on all-cause mortality than PM10. Effects of high PM pollution on mortality were stronger in the elder and male.

Conclusions

Our findings provide additional relevant information on air quality monitoring and associations of PM and human health, valuable data for further scientific research in Shenzhen and for the on-going discourse on improving environmental policies.

Peer Review reports

Background

Airborne particulate matter (PM) consistently associated with adverse health effects at current levels of exposure in urban populations [14]. Air pollution has serious direct and indirect effects on public health in China [2, 58]. PM with aerodynamic diameters less than 2.5 μm (PM2.5) has become the fourth prominent threat to the health of Chinese people [9].

The range of adverse health effects of air pollution is broad [2, 10, 11]. Susceptibility to pollution may vary depending on overall health condition and age [5, 6, 1214]. Risk of various effects has been shown to increase with exposure, but there is little evidence to suggest a threshold below which no adverse health effects can be anticipated [15, 16]. The lowest concentration at which such effects begin to manifest is not much greater than the background concentration, which has been estimated at 3–5 μg/m3 for PM2.5 in the United States and western Europe [15]. Most studies on air pollution exposure and its effects on human health in China have focused on heavily polluted cities or mega-cities [8, 1719], whereas studies on cities with relatively low air pollutant concentrations are rare.

Shenzhen is a major coastal city with a population of some 15 million. It is situated within the Pearl River Delta (PRD) and Guangdong Province, immediately north of Hong Kong. Shenzhen has become China’s most crowded city and is the fifth most densely populated city in the world, with a population density of 17,150 per square kilometre. Shenzhen is listed as the fourth most important economic centre among Chinese cities. As China’s first and still one of the most successful special economic zones the city has an important position in the PRD region and the country. Compared with other cities in China, air quality in Shenzhen is high. Nevertheless, the city has been experiencing elevated levels of PM pollution in recent years because of rapid economic development [20]. As one of the first-stage cities implementing the National Ambient Air Quality Standard (GB3095-2012) in 2013, Shenzhen provided real-time hourly monitoring concentrations of air pollutants to the general public since January 1 2013. According to air quality monitoring data from the China National Environmental Monitoring Center (CNEMC), respective annual average concentrations of PM10 and PM2.5 were 61.3 and 39.6 μg/m3 in Shenzhen in 2013. The annual average PM10 concentration was higher than in 2012 (52 μg/m3). However, comprehensive studies on PM in Shenzhen have been rare and there have been no systematic investigations of PM spatiotemporal patterns in Shenzhen.

We carried out a time-series analysis on daily and hourly PM concentrations and daily number of all-cause mortality (excluding accidental deaths) during the first year (2013) of National Ambient Air Quality Standard implementation in Shenzhen. Daily and hourly patterns of PM2.5 and PM10 were summarized and the daily PM2.5/PM10 ratio was calculated. Spatial distributions of PM10 and PM2.5 were investigated. Associations of PM and all-cause mortality were analysed and the susceptibility differentiated according to gender and age were addressed. The objectives were to provide daily/hourly PM10 and PM2.5 monitoring information for Shenzhen during 2013 to the general public and scientific researchers, investigate spatiotemporal characteristics of PM10 and PM2.5, evaluate changes of PM2.5/PM10 ratio, and to discover potential relationships between daily exposure to PM2.5 or PM10 and all-cause mortality.

Rang of health effects of air pollution was broad [2, 10, 11]. Susceptibility to pollution may vary with health or age [5, 6, 1214]. Risk of various effects has been shown to increase with exposure, but there is little evidence to suggest a threshold below which no adverse health effects can be anticipated [15, 16]. The lowest concentration at which such effects begin to manifest is not much greater than the background concentration, which has been estimated at 3–5 μg/m3 for PM2.5 in the United States and western Europe [15]. Most studies on air pollution exposure and its effects on human health in China have focused on heavily polluted cities or mega-cities [8, 1719], whereas studies on cities with relatively low air pollutant concentrations are rare.

Methods

Study area

Shenzhen is in southern China, 113°46–114°37E and 22°27–22°52 N, with an area of 1991.64 km2. There is a subtropical oceanic climate, with warm temperatures and abundant rainfall. Annual average temperature is 22.4 °C. The monthly average temperature in January is 15.4 °C, and 28.9 °C in July.

Data sources

Mortality data

All non-accidental mortality data for calendar year 2013 were obtained from death certificates recorded at the Shenzhen Center for Disease Control and Prevention. In the death registry, causes are coded by the International Classification of Disease revision 10 (ICD10).

Air pollutant monitoring data

Daily air quality monitoring data were provided by the Shenzhen Environmental Monitoring Center and CNEMC. Daily PM10 and PM2.5 concentrations were derived from the average of available hourly data measured at 11 state-controlled monitoring stations across Shenzhen, the locations of which are presented in Fig. 1. According to Technical regulation for ambient air quality assessment (on trial) (HJ 663—2013), when calculating daily means of a city, at least 75 % hourly concentrations from the monitoring stations of the city had to be available in a single day. If more than 25 % of the data in a monitoring station was missing in the whole study period, the entire station would be excluded. According to technical guidelines of the Chinese government, these locations must not be in the immediate vicinity of traffic intersections or major industrial polluters, and should be sufficiently distant from any other emission sources. Thus, the monitoring data reflect the general background urban air pollution level in our study area.

Fig. 1
figure 1

Distribution of 11 state-controlled air quality monitoring sites in Shenzhen

To discern spatiotemporal changes of hourly PM2.5 and PM10 concentrations, we acquired hourly monitoring data from 1 January through 30 November 2013 at two state-controlled monitoring stations. These stations were downtown (Huaqiaocheng, HQC) and in a tourist area (Nan’ao, NA). Hourly PM2.5 and PM10 monitoring data were from the National Real-Time Air Quality Monitoring Data Publishing Platform developed by CNEMC, which is publicly accessible via the website http://113.108.142.147:20035/emcpublish/.

Meteorological data

To control for effects of weather on mortality, meteorological data (temperature, relative humidity, barometric pressure and wind speed) were obtained from the Meteorological Bureau of Shenzhen Municipality. The weather data was monitored at a weather station belonging to that bureau. The monitoring standard is consistent with the international WMO (World Meteorological Organization) standard. There were no missing meteorological data.

Data analysis

Statistical analysis

Spearman correlation coefficients were used to reflect the relationship between PM and meteorological factors during the study period.

Spatial analysis

In the Macroscopic regional scale, spatial distribution of PM2.5 concentration follow the basic assumption of ‘the first law of geography’, namely the regional concentrations in nearby areas are more similar than in the more distant areas. Therefore, inverse distance weighted model (Inverse Distance Weighted, IDW) interpolation analysis was used to analyze spatial distributions of PM2.5 and PM10.

Associations between daily concentration of PM and mortality

Consistent with other time-series studies [21, 22], we used a generalized additive model (GAM) with penalized splines to analyze mortality, PM, and confounding factors (calendar time, day of week, temperature, barometric pressure, wind speed and humidity). Because the daily mortality number was small and typically followed a Poisson distribution [23, 24], the core analysis was via a GAM with log link and Poisson error that accounted for smooth fluctuations of that number.

In preparation for conducting the model analyses, we conducted two steps in the procedure of the model building and model fit: development of the best base model (without a pollutant) and development of the main model (with a pollutant). The latter is achieved by adding the PM to the final cause-specific best base model, assuming a linear relationship between the logarithmic mortality number and PM concentration.

First, we constructed the basic pattern of mortality number excluding PM. We incorporated smoothed spline functions of time and weather conditions, which can include non-linear and non-monotonic links between mortality and time/weather conditions, offering a flexible modelling tool. Day of the week was also included in the basic models.

After we established the basic models, we introduced the PM and analyzed their associations with mortality. To compare the relative quality of the mortality predictions across these non-nested models, Akaike’s Information Criterion (AIC) was used as a measure of how well the model fitted the data. Smaller AIC values indicate the preferred model. Briefly, we fitted the following log-linear generalized additive models to obtain the estimated pollution log-relative rate β in the study district:

$$ \log \left[\mathrm{E}\left(\mathrm{Y}\mathrm{t}\right)\right]=\upalpha +{\displaystyle \sum_{\mathrm{i}=1}^{\mathrm{q}}}\upbeta \mathrm{i}\left(\mathrm{Xi}\right)+{\displaystyle \sum_{\mathrm{j}=1}^{\mathrm{p}}}\mathrm{f}\mathrm{j}\left(\mathrm{Z}\mathrm{j},\mathrm{d}\mathrm{f}\right)+\mathrm{W}\mathrm{t}\left(\mathrm{week}\right) $$

Here E(Yt) represents the expected number of mortality at day t; β represents the log-relative rate of mortality associated with an unit increase of PM; Xi indicates the concentrations of pollutants at day t; Wt(week) is the dummy variable for day of the week. \( {\displaystyle {\sum}_{\mathrm{j}=1}^{\mathrm{p}}}\mathrm{f}\mathrm{j}\left(\mathrm{Z}\mathrm{j},\mathrm{d}\mathrm{f}\right) \) is the non-parametric spline function of calendar time, temperature, barometric pressure, wind speed and humidity. A detailed introduction to GAM is given in Wood [24]. We initialized the df as 7 df/year for time, 3 df for temperature, barometric pressure, wind speed and humidity [25].

Results were expressed as excess risk (ER) in mortality number per 10 μg/m3 increases in PM concentrations (ER = (eβxΔC-1) × 100, where ΔC is the incremental PM amount, which was 10 μg/m3 here for comparison with similar studies in other locations of China).

Values of p < 0.05 were considered statistically significant.

We also examined PM effects with different lag (L) structures of single-day (distributed lag; L0–L3) and multi-day (moving average lag; L01–L03) lags. Here, a lag of 0 day (L0) corresponds to current-day pollution and a lag of 1 day to the previous-day concentration. In multi-day lag models, L03 corresponds to a 4-day moving average pollutant concentration of the current and previous 3 days [26, 27]. Meteorological factors used in the lag models (distributed and moving average) were from current-day data. While running the models we also considered lags of more than three days for each of the pollutants, but very few associations were identified and these results have been excluded from further analyses.

Software used

Temporal changes of PM2.5 and PM10 were summarized by Origin 9.0 software, and their spatial differences were presented by ArcGIS 10.2 using Inverse distance weighted (IDW). Other statistical analyses were conducted in R3.1.0, and MGCV package in R3.1.0 was used for the GAM analysis.

Results

Descriptive results

Table 1 summarizes annual means and percentages of daily mortality number, PM10 and PM2.5 concentrations, and meteorological factors for Shenzhen in 2013.

Table 1 Statistical characteristics of air pollutants, meteorological factors, and daily mortality number

During the study period, mean daily temperature and humidity were 23.1 °C and 74.8 %, respectively. Mean daily temperature was 9.8–31.2 °C and mean daily humidity ranged from 24–100 %, reflecting the subtropical oceanic climate of Shenzhen.

Annual average concentrations were 61.3 μg/m3 for PM10 and 39.6 μg/m3 for PM2.5. Averages were higher than median values for the two air pollutants. PM2.5/PM10 ratios ranged from 37.0 to 88.3 %, with a mean of 62.4 %.

A total mortality of 11,919 people from all causes was observed in 2013. Among these, 7494 were male and 4425 female. There were 6421 people in the 0–65 year age group and 5498 in the 65+ group. The daily number of all-cause mortality was between 9 and 51.

Spearman correlation coefficients for PM and meteorological factors are presented in Table 2. Significant positive correlations were found between PM and barometric pressure. Temperature was negatively and significantly correlated with PM. Similar patterns of correlations were also found for humidity and wind speed.

Table 2 Spearman correlation coefficients between air pollutants and meteorological factors

Temporal changes

Daily concentrations of PM10 and PM2.5

Figure 2 shows temporal characteristics of PM10 and PM2.5. Daily PM10 and PM2.5 concentrations showed significant similar temporal trends, with relatively high levels during October–December and low levels for May–September.

Fig. 2
figure 2

Daily concentrations of PM10 and PM2.5 and daily number of mortality

Daily PM10 concentrations were 10–179 μg/m3 with an average of 61.3 μg/m3, and PM2.5 concentrations were 9–135 μg/m3 with an average of 39.6 μg/m3. PM10 and PM2.5 concentration ranges were wide, and their maxima were twice the Class 2 limits of National Ambient Air Quality Standard.

Ratios of PM2.5 to PM10

Temporal characteristics of daily PM2.5/PM10 ratios are shown in Fig. 2. The ratios peaked during December–February and April–May, with low values from June to August, which mean high fine particulate ratio in December–February and April–May. The ratios ranged from 37.0 to 88.3 %, with an average of 62.4 %.

Temporal trends of daily mortality number

Daily trends on number of mortality for all-cause, male, female, young and elder are also summarized in Fig. 2. Daily mortality number for all-cause was 9–51 with an average of 33, 5–36 for male with an average of 21, 1–25 for female with an average of 12, 3–32 for young with an average of 18, 4–29 for elder with an average of 15.

Spatial differences

Figure 3 presents spatial distributions of PM10 and PM2.5 in Shenzhen during 2013. PM10 and PM2.5 concentrations were higher in western parts of Shenzhen than in eastern parts. According to National Ambient Air Quality Standard, annual average PM2.5 concentrations at five monitoring stations named Nan’ao, Liyuan, Kuiyong, Meisha and Yantian were within Class 2 limits (annual average PM2.5 concentration < 35 μg/m3), but exceeded these limits at the other six stations.

Fig. 3
figure 3

Spatial distributions of PM10 and PM2.5 in Shenzhen

To represent spatial differences of PM in Shenzhen more directly, we analysed monthly and hourly concentrations of PM10 and PM2.5 at HQC and NA sites, in the downtown and tourist areas of Shenzhen, respectively. These two parts of the city serve distinct and quite different urban functions, and differences in air quality might indicate that they are affected by different pollutant emission sources.

Monthly differences

Figure 4 shows average monthly PM10 and PM2.5 concentrations at HQC and NA from January through November. The results show that PM10 and PM2.5 had similar hourly trends at HQC, which had higher concentrations in January and October and lower concentrations in July. At the NA monitoring station, both PM10 and PM2.5 had higher concentrations in January and October, and lower concentrations in May for PM10 and June for PM2.5. Concentrations of PM10 at HQC were higher than at NA during January–May and September–November, and lower than at NA during May–September. Concentrations of PM2.5 at HQC were higher than at NA during January–June and September–November, and lower than at NA during June–September.

Fig. 4
figure 4

Monthly/hourly PM10 and PM2.5 concentrations at HQC and NA

Hourly differences

Hourly average PM10 and PM2.5 concentrations at HQC and NA in 2013 are also presented in Fig. 4. Changes in hourly concentrations of PM10 and PM2.5 had similar patterns at NA, which had maxima at 20:00 and secondary maxima at 9:00. Minima were from 3:00 to 4:00. Changes in hourly concentration of PM10 and PM2.5 did not show patterns common with HQC, which may be related to changes of pollution source over a day in the downtown area. Hourly concentrations of PM2.5 at NA were lower than at HQC throughout the day, but the hourly PM10 concentration did not show this pattern.

Associations with all-cause mortality

Table 3 presents ER percentages(ERs) (95 % confidence interval [CI]) of daily all-cause mortality number with every 10-μg/m3 increase in PM10 or PM2.5 concentration.

Table 3 Excess risk (ER) percentage for daily all-cause mortality number with every 10-μg/m3 increase in PM concentration

To identify possible time delay of PM10 or PM2.5 pollution exposure and daily mortality number, we analyzed lag effects of air pollutants. ER in the all-cause mortality number with a 10 μg/m3 increase of pollutants for single-day measures, 1–3 days prior to mortality (L0–L3), and moving averages from day 0 and day 1 to day 3 prior to the mortality are also listed in Table 3. When running the models, lag effects of more than 3 days for PM10 and PM2.5 were also considered. However, as little to no relationship was found, the results of that analysis were not included. Gender and age differences were also considered. Unlike cities in northern China, temperature differences in Shenzhen were not significant. Therefore, we did not run seasonal models.

The results showed that the greatest ERs of PM10 and PM2.5 were in 2-day cumulative measures (L02) for the all-cause mortality group, 2-day lag (L2) for females and the young (0–65 years), and L02 for males and the elder (> 65 years). The greatest ERs in the mortality number with a 10-μg/m3 increase of PM10 were 0.61, 0.33, 0.85, 1.55 and 0.21 % for the all-cause mortality group, females, males, elder and young, respectively. The greatest ERs in the mortality number with a 10-μg/m3 increase of PM2.5 were 0.69, 0.76, 0.74, 1.97 and 0.28 % for the same respective groups. ERs of males with increases in PM10 or PM2.5 concentration were greater than those of females, and ERs of the elder were greater than the young with concentration increases of PM10 or PM2.5.

Discussion

This study focused on spatiotemporal patterns and possible associations of PM10 and PM2.5 with all-cause mortality during the first year (2013) of National Ambient Air Quality Standard implementation in Shenzhen, a relatively clean city compared to other cities in China. The objectives were to provide 2013 PM monitoring information of Shenzhen to the general public, to discover possible associations between PM and mortality in a comparatively clean city, and to provide scientific results to researchers in other areas. We also intend to encourage health services and public health policymakers in Shenzhen to consider ideas for real-time public health alerts for air quality, so that vulnerable groups and others affected by air pollution can be appropriately advised. The present study was unique in the following aspects: 1) we analysed daily patterns of PM and air quality during the first year (2013) of National Ambient Air Quality Standard implementation in Shenzhen; 2) based on reliable data sources, hourly/monthly patterns of PM in two functional areas of Shenzhen were addressed; 3) spatial patterns of PM were determined; and 4) to our knowledge, the study is the first to investigate associations between PM and all-cause mortality in Shenzhen.

During 2013, annual average PM10 and PM2.5 concentrations were 61.3 and 39.6 μg/m3, respectively; averages were higher than the median values for the two air pollutants. According to National Ambient Air Quality Standard and the Technical Regulation on Ambient Air Quality Index (on trial) (HJ633-2012), PM2.5 was the major air pollutant in Shenzhen, with 104 days as the primary pollutant and 32 days as a “non-attainment” pollutant. The latter indicates 32 days with heavy PM2.5 pollution. Annual average PM10 and PM2.5 concentrations were 108 and 89 μg/m3 in Beijing, and 72 and 53 μg/m3 in Guangzhou [28]. Compared with heavy PM-polluted cities in China (Beijing, Guangzhou, and others), Shenzhen has good air quality [9, 19, 29]. Shenzhen was ranked 7th among 74 first-stage cities, but its annual average PM2.5 concentration exceeded the Class 2 limit of National Ambient Air Quality Standard [28, 30].

The annual average ratio of PM2.5 to PM10 was 62.4 %, which indicates a high percentage of PM2.5 in ambient air pollution of Shenzhen. PM2.5/PM10 maximized in December–February and April–May, with lower values in June–August. Compared with Beijing the PM2.5/PM10 ratio in Shenzhen was higher than Beijing autumn normal days and lower than haze days and winter normal days in Beijing; the average PM2.5/PM10 ratios in Beijing were correspondingly 0.63, 0.32, 0.70, and 0.66 in autumn haze, autumn normal, winter haze and winter normal days, respectively [31]. The PM2.5/PM10 ratio is 0.575 in Taiwan [32]. These findings may be related to meteorological conditions and pollution sources in the city. PM2.5 concentrations can be affected by both local emissions and contributions of meso-scale origin [33]. Further studies on concentrations and ratios of PM at intercity level should be conducted.

Hourly concentrations of PM10 and PM2.5 had similar patterns in the tourist area (NA monitoring station) but did not have any patterns in common in the downtown area. This may be related to changes of pollution source like traffic emission in downtown over a single day. Hourly concentrations of PM2.5 in the tourist area were lower than downtown throughout the day, which may be attributed to more intensive human activities downtown. PM10 and PM2.5 concentrations were higher in western parts of Shenzhen than eastern parts, which may be related to land use, pollution sources, industrial structure, traffic conditions, and other factors [34]. There should be further study of relationships between driving factors (e.g., spatial distribution of pollutant emissions, pollutant emission intensity, and regional industrial structure) and pollutant concentrations.

Time-series studies estimate that a 10 μg/m3 increase in mean 24-hour PM2.5 concentration increases the ERs of daily cardiovascular mortality by ~0.4 to 1.0 % [35]. Consistent with other studies [2, 11, 12, 36, 37], we found a statistically significant association between PM10 or PM2.5 and daily mortality number. There were lag effects in all the study groups, and ERs of PM2.5 were greater than PM10 in all study groups with concentration increases. ERs in the all-cause mortality number with a 10-μg/m3 increase of PM10 and PM2.5 were 0.61 % (95 % CI: 0.50 %–0.72 %) and 0.69 % (95 % CI: 0.55 %–0.83 %), respectively. A study in Beijing for 2005–2009 showed that a 10 μg/m3 increase in PM2.5 was associated with a 0.65 % rise in all-cause mortality, whereas the same increase in PM10 was associated with an increase of 0.15 % [38]. During 2006–2009 in Guangzhou, increments of 10 μg/m3 in PM10 were associated with a ER of 1.26 % for total non-accidental deaths, and 1.79 % for cardiovascular deaths [39]. During 2007 to 2009 in Tianjin, the effect estimates per 10 μg/m3 increase in PM10 concentrations at the moving average of lags 0 and 1 day in high temperature level were 0.62 % for non-accidental mortality [40]. ERs in our study were greater than those reported for Beijing and smaller than those for Guangzhou. Such inter-city variability in ER estimates may have been influenced by a number of factors, such as demographic and socioeconomic variables, culture, air pollution sources, and geographical and weather conditions [39]. Both temperature and particulate air pollution are associated with increased death risk; and extreme high temperature increased the associations of PM10 with daily mortality [41]. ERs in the all-cause mortality number with 10-μg/m3 increase of PM10 and PM2.5 were 0.33 and 0.76 % for females, respectively, 0.85 and 0.77 % for males, 1.55 and 1.97 % for the elder, and 0.21 and 0.28 % for the young. Males were more sensitive to PM10 or PM2.5 concentration changes than females. The elder appeared to be more affected than the young by PM10 or PM2.5 concentration increase. These findings indicate that PM effects on mortality were stronger among the elderly and on male. Because the seasonal difference was not significant in Shenzhen, we did not consider seasonal associations among PM and mortality in this study. The present study has certain limitations. We considered the target population to be relatively homogeneous and did not consider residence or work location of deaths, owing to a lack of data. Pollutant exposure levels were derived from 11 fixed-site monitoring stations. However, because air pollution varies spatially within a city, averages drawn from these stations may not reflect actual exposure levels. Accurate exposure assessment and a homogeneous target population are important factors to consider in future studies estimating mortality risk from air pollution [39]. Further in-depth studies should require air pollutant composition, pollution emission sources, pollutant emission patterns, time-series of human activity, individual exposure to pollutants, social economy, and human health at the city level.

Conclusion

During 2013, annual average PM10 and PM2.5 concentrations were 61.3 and 39.6 μg/m3 in Shenzhen. PM2.5 failed to meet the Class 2 annual limit of National Ambient Air Quality Standard and was the major air pollutant, with 8.8 % of days having heavy PM2.5 pollution. The annual average PM2.5/PM10 ratio was 62.4 %. Hourly PM2.5 concentrations in the tourist area were lower than downtown throughout the day. PM10 and PM2.5 concentrations were higher in western parts than eastern parts. ERs in the all-cause mortality number increased with PM10 and PM2.5. PM2.5 had higher risks than PM10. PM effects on mortality were stronger among male and the elderly. Our findings provide additional information on air quality monitoring and associations between PM and all-cause mortality, and valuable data for scientific research in Shenzhen. It also contributes to the discussion on further developing environmental health policies in urban China.

References

  1. Dominici F, Greenstone M, Sunstein CR. Particulate Matter Matters. Science. 2014;344(6181):257–9.

    Article  PubMed  PubMed Central  Google Scholar 

  2. Shang Y, Sun Z, Cao J, Wang X, Zhong L, Bi X, et al. Systematic review of Chinese studies of short-term exposure to air pollution and daily mortality. Environ Int. 2013;54:100–11.

    Article  CAS  PubMed  Google Scholar 

  3. Bell ML, Zanobetti A, Dominici F. Evidence on vulnerability and susceptibility to health risks associated with short-term exposure to particulate matter: a systematic review and meta-analysis. Am J Epidemiol. 2013;178:kwt090.

    Article  Google Scholar 

  4. Kim K-H, Kabir E, Kabir S. A review on the human health impact of airborne particulate matter. Environ Int. 2015;74:136–43.

    Article  CAS  PubMed  Google Scholar 

  5. Tang D, Li TY, Chow JC, Kulkarni SU, Watson JG, Ho SSH, et al. Air pollution effects on fetal and child development: A cohort comparison in China. Environ Pollut. 2014;185:90–6.

    Article  CAS  PubMed  Google Scholar 

  6. Zhou M, Liu Y, Wang L, Kuang X, Xu X, Kan H. Particulate air pollution and mortality in a cohort of Chinese men. Environ Pollut. 2014;186:1–6.

    Article  CAS  PubMed  Google Scholar 

  7. Dong G-H, Qian ZM, Xaverius PK, Trevathan E, Maalouf S, Parker J, et al. Association between long-term air pollution and increased blood pressure and hypertension in China. Hypertension. 2013;61(3):578–84.

    Article  CAS  PubMed  Google Scholar 

  8. L-w Z, Chen X, Xue X-d, Sun M, Han B, Li C-p, et al. Long-term exposure to high particulate matter pollution and cardiovascular mortality: A 12-year cohort study in four cities in northern China. Environ Int. 2014;62:41–7.

    Article  Google Scholar 

  9. Chen Z, Wang J-N, Ma G-X, Zhang Y-S. China tackles the health effects of air pollution. Lancet. 2013;382(9909):1959–60.

    Article  PubMed  Google Scholar 

  10. Pope CA, Dockery DW. Health effects of fine particulate air pollution: Lines that connect. J Air Waste Manage Assoc. 2006;56(6):709–42.

    Article  CAS  Google Scholar 

  11. Wagner JG, Allen K, Yang H-y, Nan B, Morishita M, Mukherjee B, et al. Cardiovascular Depression in Rats Exposed to Inhaled Particulate Matter and Ozone: Effects of Diet-Induced Metabolic Syndrome. Environ Health Perspect. 2014;122(1):27–33.

    Article  PubMed  Google Scholar 

  12. Goldberg MS, Burnett RT, Stieb DM, Brophy JM, Daskalopoulou SS, Valois M-F, et al. Associations between ambient air pollution and daily mortality among elderly persons in Montreal, Quebec. Sci Total Environ. 2013;463:931–42.

    Article  PubMed  Google Scholar 

  13. Mahiyuddin WRW, Sahani M, Aripin R, Latif MT, Thuan-Quoc T, Wong C-M. Short-term effects of daily air pollution on mortality. Atmos Environ. 2013;65:69–79.

    Article  Google Scholar 

  14. Burnett RT, Pope III CA, Ezzati M, Olives C, Lim SS, Mehta S, et al. An Integrated Risk Function for Estimating the Global Burden of Disease Attributable to Ambient Fine Particulate Matter Exposure. Environ Health Perspect. 2014;122(4):397–403.

    PubMed  PubMed Central  Google Scholar 

  15. WHO Regional Office for Europe. Air quality guidelines, global update 2005: particulate matter, ozone, nitrogen dioxide, and sulfur dioxide. World Health Organization; 2006. http://www.euro.who.int/en/health-topics/environment-and-health/air-quality/publications/pre2009/airquality-guidelines.-global-update-2005.-particulate-matter,-ozone,-nitrogen-dioxide-and-sulfur-dioxide.

  16. Brunekreef B, Holgate ST. Air pollution and health. Lancet. 2002;360(9341):1233–42.

    Article  CAS  PubMed  Google Scholar 

  17. Huang W, Cao J, Tao Y, Dai L, Lu S-E, Hou B, et al. Seasonal Variation of Chemical Species Associated With Short-Term Mortality Effects of PM2.5 in Xi’an, a Central City in China. Am J Epidemiol. 2012;175(6):556–66.

    Article  PubMed  PubMed Central  Google Scholar 

  18. Cao J, Xu H, Xu Q, Chen B, Kan H. Fine Particulate Matter Constituents and Cardiopulmonary Mortality in a Heavily Polluted Chinese City. Environ Health Perspect. 2012;120(3):373–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Ma Y, Chen R, Pan G, Xu X, Song W, Chen B, et al. Fine particulate air pollution and daily mortality in Shenyang, China. Sci Total Environ. 2011;409(13):2473–7.

    Article  CAS  PubMed  Google Scholar 

  20. Dai W, Gao J, Cao G, Ouyang F. Chemical composition and source identification of PM2.5 in the suburb of Shenzhen, China. Atmos Res. 2013;122:391–400.

    Article  CAS  Google Scholar 

  21. Zhang FY, Wang WY, Lv JM, Krafft T, Xu J. Time-series studies on air pollution and daily outpatient visits for allergic rhinitis in Beijing, China. Sci Total Environ. 2011;409(13):2486–92.

    Article  CAS  PubMed  Google Scholar 

  22. Bhaskaran K, Gasparrini A, Hajat S, Smeeth L, Armstrong B. Time series regression studies in environmental epidemiology. Int J Epidemiol. 2013;42(4):1187–95.

    Article  PubMed  PubMed Central  Google Scholar 

  23. Box GE, Jenkins GM, Reinsel GC. Time series analysis: forecasting and control. 4th ed. New York: John Wiley & Sons; 2008.

  24. Wood S. Generalized additive models: an introduction with R. CRC Press; 2006. https://www.crcpress.com/Generalized-Additive-Models-An-Introduction-with-R/Wood/9781584884743.

  25. Yang Y, Cao Y, Li W, Li R, Wang M, Wu Z, et al. Multi-site time series analysis of acute effects of multiple air pollutants on respiratory mortality: A population-based study in Beijing, China. Sci Total Environ. 2015;508:178–87.

    Article  CAS  PubMed  Google Scholar 

  26. Gasparrini A, Armstrong B, Kenward M. Distributed lag non‐linear models. Stat Med. 2010;29(21):2224–34.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Zhang F, Krafft T, Ye B, Zhang F, Zhang J, Luo H, et al. The lag effects and seasonal differences of air pollutants on allergic rhinitis in Beijing. Sci Total Environ. 2013;442:172–6.

    Article  PubMed  Google Scholar 

  28. China National Environmental Monitoring Center. Environmental Quality Report of China in 2013. 2014.

    Google Scholar 

  29. Huang D, Xu J, Zhang S. Valuing the health risks of particulate air pollution in the Pearl River Delta, China. Environ Sci Pol. 2012;15(1):38–47.

    Article  CAS  Google Scholar 

  30. Ministry of Environmental Protection of the People’s Republic of China M. Report on the State of the Environment in China 2013. 2014.

    Google Scholar 

  31. Gao J, Tian H, Cheng K, Lu L, Zheng M, Wang S, et al. The variation of chemical characteristics of PM2.5 and PM10 and formation causes during two haze pollution events in urban Beijing, China. Atmos Environ. 2015;107:1–8.

    Article  CAS  Google Scholar 

  32. Chu H-J, Huang B, Lin C-Y. Modeling the spatio-temporal heterogeneity in the p M10–PM2.5 relationship. Atmos Environ. 2015;102:176–82.

    Article  CAS  Google Scholar 

  33. Marcazzan GM, Vaccaro S, Valli G, Vecchi R. Characterisation of PM10 and PM2.5 particulate matter in the ambient air of Milan (Italy). Atmos Environ. 2001;35(27):4639–50.

    Article  CAS  Google Scholar 

  34. Zhang J, Ouyang Z, Miao H, Wang X. Ambient air quality trends and driving factor analysis in Beijing, 1983-2007. J Environ Sci (China). 2011;23(12):2019–28.

    Article  CAS  Google Scholar 

  35. Brook RD, Rajagopalan S, Pope III CA, Brook JR, Bhatnagar A, Diez-Roux AV, et al. Particulate Matter Air Pollution and Cardiovascular Disease An Update to the Scientific Statement From the American Heart Association. Circulation. 2010;121(21):2331–78.

    Article  CAS  PubMed  Google Scholar 

  36. Atkinson RW, Kang S, Anderson HR, Mills IC, Walton HA. Epidemiological time series studies of PM2.5 and daily mortality and hospital admissions: a systematic review and meta-analysis. Thorax. 2014;69(7):660–5.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Krall JR, Anderson GB, Dominici F, Bell ML, Peng RD. Short-term exposure to particulate matter constituents and mortality in a national study of US urban communities. Environ Health Perspect. 2013;121(10):1148.

    PubMed  PubMed Central  Google Scholar 

  38. Li P, Xin J, Wang Y, Wang S, Shang K, Liu Z, et al. Time-series analysis of mortality effects from airborne particulate matter size fractions in Beijing. Atmos Environ. 2013;81:253–62.

    Article  CAS  Google Scholar 

  39. Yu ITS, Zhang Y, San Tam WW, Yan QH, Xu Y, Xun X, et al. Effect of ambient air pollution on daily mortality rates in Guangzhou, China. Atmos Environ. 2012;46:528–35.

    Article  CAS  Google Scholar 

  40. Li G, Zhou M, Cai Y, Zhang Y, Pan X. Does temperature enhance acute mortality effects of ambient particle pollution in Tianjin City, China. Sci Total Environ. 2011;409(10):1811–7.

    Article  CAS  PubMed  Google Scholar 

  41. Meng X, Zhang Y, Zhao Z, Duan X, Xu X, Kan H. Temperature modifies the acute effect of particulate air pollution on mortality in eight Chinese cities. Sci Total Environ. 2012;435–436:215–21.

    Article  PubMed  Google Scholar 

Download references

Acknowledgments

The authors wish to thank all the staff members at the Shenzhen Centre for Disease Control and Prevention for their strong support of this study. We thank the Shenzhen Environmental Monitoring Center and Meteorological Bureau of Shenzhen Municipality for providing data. The present study was supported by the National Natural Science Foundation of China (NO. 41401101 & NO.41371118).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Fengying Zhang or Wuyi Wang.

Additional information

Competing interests

The authors declare that they have no competing interest.

Authors’ contribution

FYZ led the study, carried out the time-series studies, analyzed the data and wrote the first draft of the manuscript. XJL participated in all-cause mortality date in Shenzhen. LZ and YY conducted PM monitoring/ meteorological data collecting and spatial analysis. LW and JML assisted with statistical analysis and language editing. WYW and TK helped to conceptualize the study, provided intellectual advice, contributed to data interpretation and helped to revise various drafts of the manuscript. All authors read and approved the final manuscript.

Fengying Zhang and Xiaojian Liu are co-first authors.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, F., Liu, X., Zhou, L. et al. Spatiotemporal patterns of particulate matter (PM) and associations between PM and mortality in Shenzhen, China. BMC Public Health 16, 215 (2016). https://doi.org/10.1186/s12889-016-2725-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12889-016-2725-6

Keywords