- Open Access
Spatial variability of excess mortality during prolonged dust events in a high-density city: a time-stratified spatial regression approach
International Journal of Health Geographics volume 16, Article number: 26 (2017)
Dust events have long been recognized to be associated with a higher mortality risk. However, no study has investigated how prolonged dust events affect the spatial variability of mortality across districts in a downwind city.
In this study, we applied a spatial regression approach to estimate the district-level mortality during two extreme dust events in Hong Kong. We compared spatial and non-spatial models to evaluate the ability of each regression to estimate mortality. We also compared prolonged dust events with non-dust events to determine the influences of community factors on mortality across the city.
The density of a built environment (estimated by the sky view factor) had positive association with excess mortality in each district, while socioeconomic deprivation contributed by lower income and lower education induced higher mortality impact in each territory planning unit during a prolonged dust event. Based on the model comparison, spatial error modelling with the 1st order of queen contiguity consistently outperformed other models. The high-risk areas with higher increase in mortality were located in an urban high-density environment with higher socioeconomic deprivation.
Our model design shows the ability to predict spatial variability of mortality risk during an extreme weather event that is not able to be estimated based on traditional time-series analysis or ecological studies. Our spatial protocol can be used for public health surveillance, sustainable planning and disaster preparation when relevant data are available.
Dust events are extreme pollution events that can induce adverse health effects in global cities, as studies have found for major cities in the United States, Australia, Asia and Europe [6, 12, 16, 31, 34, 40, 49]. Previous studies extensively applied temporally stratified models to quantify mortality risk during dust events, with promising results [31, 40, 49]; for example, one study showed that extreme events can lead to a 16% increase in dust mortality in a downwind city . However, there has been no study to investigate the spatial variability of mortality risk during a prolonged dust event. In contrast to a lack of research on this issue, extensive environmental health studies on extreme weather and pollution have pointed out the necessity of predicting spatial variability of mortality and morbidity [2, 20, 23, 27, 29], for the purpose of measuring community vulnerability and public health planning. Estimating community vulnerability is particularly important to a high-density city, as the urban morphology of a high-density city influences air pollutant dispersion , resulting in extreme conditions with severe health risk, in particular at the district level. The prolonged effect of dust combined with urban morphology is expected to induce an additional effect on mortality in a downwind city.
In order to measure community vulnerability, both governmental health guidelines and previous studies have suggested investigating the environmental and socioeconomic factors that can additionally elevate the health risk [11, 17, 19, 24, 38, 41, 42]. Some studies have also proposed combining significant environmental and socioeconomic factors to pinpoint the hotspots of health risk [26, 44, 47, 50]. In order to enhance health planning to minimize adverse health effects of prolonged dust events, this study develops a set of protocols for (1) evaluating potential environmental and socioeconomic factors that can elevate mortality risk during a prolonged dust event, (2) including spatial influences of neighboring communities to adjust for environmental and socioeconomic effects on mortality risk, and (3) locating communities with higher mortality risk for disaster risk management during future dust episodes. The approach developed from this study could be applied to other regions where data on city-specific environmental and socioeconomic factors are available.
Urban and climate settings of Hong Kong
Hong Kong is a typical high-density city located in a subtropical region. There have been ten reported days with dust events in the past decade in Hong Kong, including 2 days in 2006 (Apr 16–17, 2006), 4 days in 2009 (Apr 27–30, 2009), and 4 days in 2010 (Mar 23–26, 2010) [53, 54]. Two of these three dust events (8 of 10 days) were prolonged dust events with ≥3 consecutive dusty days. Significant mortality risk was observed on those days, with 7% increase in all-cause mortality and 7% increase in cardiorespiratory mortality during a dusty day. There was also significant air pollution during those dusty days, with average PM10–2.5 concentrations 147.6% higher than the days without dust. There is also an extreme population pattern in Hong Kong. The population density of Hong Kong is approximately 6500 persons per km2. The significant clustering of the urban population in Hong Kong potentially introduces significant intra-urban differences in mortality risk, due to the built environment and demographic structure [13, 33, 52].
Evaluation of community and environmental characteristics related to mortality risk
Mortality for dusty days and non-dusty days was calculated for each of 287 tertiary planning units (TPU), which is the smallest spatial unit in Hong Kong with mortality and census data available, in order to measure region-specific mortality risks across the entire territory of Hong Kong. The all-cause mortality dataset for dusty days was retrieved from mortality data of the Hong Kong Census and Statistics Department, based on 8 dust days (Apr 27–30, 2009 and Mar 23–26, 2010) associated with prolonged dust events (≥3 consecutive dust days), and by excluding all traffic-related deaths (ICD-10 codes V01–V99) during these prolonged dust episodes. Mortality on non-dusty days for each TPU was used to represent the baseline mortality. Deaths of the same weekday of four control weeks before and four control weeks after each dusty day were used to represent the mortality on non-dusty days, in order to minimize bias of seasonality and weekday/weekend effect; and this was divided by the numbers of control weeks for the purpose of comparison with the total mortality on all dusty days in Hong Kong.
To evaluate the potential community and environmental factors that influence mortality risk during a dust storm in Hong Kong, multivariate linear regression was firstly applied to estimate total mortality during the prolonged dust events. Six variables were included as independent variables to evaluate the environmental and socioeconomic effects on mortality in each TPU: (1) sky view factor (SVF), (2) percentage of vegetation, (3) land surface temperature (LST), (4) percentage of low education, (5) percentage of low income, and (6) percentage of elderly:
where total mortality is the total mortality within 8 days of each TPU, average SVF is the average SVF of each TPU, % vegetation is the percentage of vegetation of each TPU, average LST is the average LST of each TPU, % low education is the percentage of low education population of each TPU, % low income is the percentage of low income population of each TPU and % elderly is the percentage of elderly of each TPU.
Urban geometric characteristics can be depicted by different parameterization indices such as building height, building density, frontal area index (FAI), planar area index (PAI), height/width ratio (H/W), and SVF. The SVF is an indicator representing combinations of building height, building density and topography . The SVF is a ratio to measure the openness of a particular area within an urban setting and in general a terrestrial landscape, which has significant implications for the incoming and outgoing radiation . SVF has been widely used in urban climate research [15, 21, 25, 28, 35, 45, 56], especially to improve spatial models of air pollution prediction [18, 39]. In this study, the SVF was used to locate high-density environments that can potentially trap air pollutants and prohibit air ventilation during prolonged dust events. Average SVF of each TPU was calculated based a raster-based SVF image of Yang et al.  derived from airborne Lidar data (Fig. 1). The SVFs at both rooftop and ground levels of this raster-based image were estimated and the spatial resolution of the airborne LIDAR data is 1 m. The building GIS data of Hong Kong were used to calculate the SVF for vertical facets using the planar area index (PAI) and the frontal area index (FAI).
Vegetation coverage (measured in percentage) can potentially influence or absorb ground-level air pollution in each TPU. In this study, the territory-wide vegetation coverage was estimated using the land use and land cover map derived from the Planning Department of Hong Kong (Fig. 2).
Land surface temperature (LST) images are commonly used to represent spatial variations of surface temperature that can affect health risk [29, 36, 55]. Landsat Thematic Mapper TM 5 on March 25, 2010 was used to estimate LST to demonstrate typical temperature variations during prolonged dust events in Hong Kong. Average LST (Fig. 3) was estimated using an improved urban emissivity model based on the SVF .
Lower education is associated with higher social vulnerability during air pollution events [32, 52], and lower income is related to low socioeconomic status, which may induce adverse health effects on a day with heavy pollution . The elderly are identified as one of the major age groups that are highly vulnerable during days with heavy air pollution [8, 32]. Therefore, the percentages of low education, low income and elderly were extracted from the 2006 census data of Hong Kong, and were used to represent the socioeconomic influence of each TPU. The percentage of low education was calculated based on the percentage of persons who had a primary school education or less (Fig. 4). The percentage of low income was the percentage of persons who were unpaid or had monthly income lower than HKD$10,000 (Fig. 5). The percentage of elderly was the percentage of persons aged ≥65 in each TPU (Fig. 6). Two TPUs with missing data of low education, low income and elderly were excluded from this study.
Finally, the predicted total mortality and the 95% confidence interval (CI) were estimated to represent the additional effect of mortality risk on both dusty days and non-dusty days (baseline) from the spatial variability of each variable. Excess mortality between dusty days and non-dusty days contributed by each spatial factor was also reported in this study (Fig. 7).
Including neighboring effects for mortality risk estimation
To include neighboring effects on mortality risk of each TPU, a spatial error model was applied and was compared with the results of multivariate linear regression. The spatial error model incorporated spatial autocorrelation in a regression error term (Lambda) to adjust spatial dependence in a multivariate linear regression model . In our study, the spatial error model weights the neighboring TPUs to spatially adjust socioeconomic and environmental influences on total mortality. To spatially weight the TPUs, we applied the queen contiguity method. This method weights all spatial neighbors with shared borders and corners , based on a spatial distance with the order of contiguity. This study applied the queen contiguity method to adjust the mortality risk of each TPU based on all surrounding TPUs. To evaluate the appropriate spatial distance for adjustment, the 1st to 3rd orders of contiguity (lag 1–lag 3) were used to estimate total mortality and for comparison with the linear results. All significant environmental and socioeconomic variables were used to construct the spatial error models with the 1st–3rd order of contiguity, and the multivariate linear regression for comparison. The Akaike information criterion (AIC)  was then adopted to compare the models, in which lower AIC indicated better model performance. Total mortality and CIs representing the additional effects of mortality risk from the spatial variability of each variable were also reported, and the differences between models were further evaluated. Finally, predicted mortality change (increase or decrease) in numbers of deaths of each TPU from the appropriate models for both dusty and non-dusty days were used to illustrate the spatial variability of relative mortality across Hong Kong during prolonged dust events.
Contributions of socioeconomic and environmental influences to mortality risk
There were in total of 802 decedents reported on all case days, and 6331 decedents reported on all control days. Based on the multivariate linear models, the SVF, percentage of low education, percentage of low income, and percentage of elderly show varied contributions to local mortality risk of each TPU during prolonged dust events (Table 1). Among all, percentage of low income is the highest risk factor of a community during prolonged dust events. A TPU with 10% more low-income population is found to have 12.5% higher mortality during a dusty day than a non-dusty day. At a TPU with 10% more population who education level was primary school or below, there is also 6.3% higher mortality during prolonged dust events than days without dust. In contrast, SVF has a negative association with mortality risk. A TPU with 10% higher SVF has 5.3% less mortality during a day with a dust storm. This indicates that a TPU with a high-density built environment generally has higher risk during prolonged dust events, while lower-density environments with higher SVF have less mortality risk.
In addition, there is no observation of an increase of total mortality in the TPUs with a higher percentage of elderly. Similar results have been found in other Hong Kong studies; for example, Chan et al.  estimated community vulnerability with census data and found that those aged ≥75 had significantly lower mortality risk than those aged <75, especially for the male population. Spatial differences in temperature and vegetation did not have significant contributions to mortality risk during prolonged dust events.
Including spatial influences for mortality risk estimation
Compared to the non-spatial model for dusty days using all variables with an AIC of 1692.3, the non-spatial model using only significant environmental and socioeconomic variables (SVF, % lower education, % lower income and % elderly) has a lower AIC of 1689.6, indicating a better model for prediction. By using all significant environmental and socioeconomic variables, it is observed that inclusion of neighboring effects as spatial influential factors has enhanced mortality risk estimation. By comparing all models for predicting mortality during dusty days with and without incorporating spatial autocorrelation (Table 2), the model considering the 1st order of queen contiguity (lag 1) has the best performance. The AIC of this lag-1 model is 1679.97, and it is the lowest among the others. There is a positive value of the regression error term of the lag-1 model (Lambda: 0.3), indicating that including neighboring effects of surrounding TPUs results in less spatial error for mortality risk prediction. It is also important to note that only including the 1st- or 2nd-order queen contiguity in modelling can enhance mortality risk estimation for prolonged dust events in Hong Kong. The 3rd-order queen contiguity does not improve the modelling, based on Lambda reported with the spatial error model. We also repeated the analyses for mortality predictions of non-dusty days (Table 3). We found similar results for model comparison, with the lag-1 model the best for predicting mortality during days without prolonged dust events (AIC: 1582.04).
By using the lag-1 model to include the neighboring effects (Table 4), areas with 10% higher SVF will have 5.3% less mortality risk than TPUs during a prolonged dust event, while a TPU with 10% more low-education population will have 6.7% higher mortality during prolonged dust events compared to non-dusty days, with all these reaching statistical significance.
Based on the comparison of spatial and non-spatial models, we applied a spatial error model incorporating the 1st order of queen contiguity to predict total mortality on dusty days and non-dusty days, and a predicted change of total mortality as relative risk of each TPU is reported in this study (Fig. 8). The mortality risk map indicates that rural areas with low-density environments have a potential decrease in mortality during prolonged dust events compared to non-dusty days. In contrast, the TPUs with high-density environments and high socioeconomic deprivation, such as TPUs in Tuen Mun, Sham Shui Po, Wong Tai Sin and Kwun Tong, generally have a higher increase in mortality during prolonged dust events compared to non-dusty days. These TPUs are predicted to have 0.1–0.5 more deaths in a period of 8 dust days than the control periods, controlling for SVF, percentage of low education, percentage of low income, percentage of elderly, and spatial autocorrelation.
This study applied a spatial regression approach to estimate spatial variability of mortality risk across a high-density city during prolonged dust events. Based on this approach, the influence of the built environment is highlighted by the negative association between SVF and mortality increase. This result indicates that high-density urban areas may trap air pollutants during days with dust storms, resulting in poorer air quality and severely increasing health risk; while areas with more openness allow better air ventilation and dispersion, therefore less health risk attributable to air pollution can be found in these areas during dusty days. Influence from socioeconomic deprivation is also determined by the positive association between total mortality, percentage of low education, and percentage of low income. In contrast, percentage of elderly of a TPU does not have a positive association with mortality risk. This might be due to the presence of more health facilities in communities with higher percentages of elderly, which reduces the mortality risk of such neighborhoods, while elderly in TPUs with lower percentages of older population may not benefit from such facilities, therefore increasing their risk. In conclusion, these findings are innovative, because previous studies only temporally stratified the dust mortality [12, 14, 31, 40, 49], without understanding the intra-urban difference in mortality risk during dust events.
In the context of spatial health planning, health risk mapping can characterize vulnerability of specific populations in a specific region , for the purpose of supporting health authorities, policymakers, and city officials to determine future health protocols in different communities . This mapping technique has been widely used along with governmental actions to develop public health surveillance. For example, the City of Toronto in Canada initiated a heat vulnerability mapping project for minimizing summer risks , and Vancouver Coastal and Fraser Health Authorities gave impetus to the development of the Vancouver Area Neighborhood Deprivation Index (VANDIX) for general health risk estimation . Mapping mortality risk adjusted for environmental and socioeconomic factors can help target a single disaster episode for comprehensive health planning. This is necessary because the general health vulnerability index can be somewhat useful, but may not be able to fully describe the spatial variability of a particular health risk . One example is that VANDIX is related to heat mortality in the Vancouver area, but it is necessary to adjust it to pinpoint heat risks with accuracy . Therefore, previous health studies indicate the need to calibrate spatial vulnerability assessments with health outcome data [5, 11, 26, 47, 50], while mortality data will be the most appropriate dataset for demonstrating disaster episodes. Therefore, mapping mortality risk during prolonged dust events is essential, since such spatial assessment can be used for local government action, as well as serving as a regional protocol for developing similar health indices in other cities.
Furthermore, previous research mostly developed health indices based on a simple spatial overlay technique ; there were issues related to spatial autocorrelation that these studies did not consider, and that might create potential biases in results. Our study adopting spatial error modelling is test-proven, with promising results showing that spatial autocorrelation can significantly improve the accuracy of predicting mortality during dust episodes. Its findings are similar to those of other spatial epidemiologic literature showing that analysis with spatial autocorrelation can help predict spatial variability of health risks [37, 51].
One limitation of this study is that the prolonged dust events were isolated episodes in Hong Kong. Based on only two dust events in Hong Kong, it was not able to employ time-series analysis for a more comprehensive spatio-temporal assessment. Application of an alternative method such as a time-stratified approach for estimating the standard mortality ratio of each small neighborhood is also problematic, since comparing rare death cases on non-dusty days in each neighborhood may create extreme estimation, resulting in statistical bias. To avoid the statistical bias, previous health studies generally applied spatial delineation techniques to stratify socioeconomic or environmental data by groups [27, 48, 52]. This method can capture spatial differences between groups, but is still insufficient to estimate the individual risk of each district. Our approach is applicable for the present case study, because we applied spatial regression to predict mortality on both dusty and non-dusty days for comparison. With the support of spatial regressions to compare total mortality between dusty and non-dusty days, the results of this study can be used to demonstrate the additional mortality effect in each district due to spatial variability of environmental and socioeconomic factors during the isolated but fatal dust events.
For future study, inclusion of spatial data on air pollution exposure may enhance mortality risk mapping. However, existing pollution mapping methods such as land use regression are limited by the spatiotemporal coverage of the data, which may not be able to demonstrate extreme cases such as a prolonged dust event. In addition, there is an accuracy issue in using such mapping methods in a high-density city, because a complex urban built environment influences air ventilation, and as a result produces bias in pollution mapping . Misuse of air pollution maps for a spatial study can induce a significant ecological fallacy, especially since community vulnerability is already influenced by the adverse effect due to aggregate-level data, instead of the association with individual-level response . In order to tackle this issue, some studies have started to use moderate-resolution satellite images for mapping Aerosol Optical Depth (AOD) or Aerosol Optical Thickness (AOT) in order to demonstrate spatiotemporal variation in air pollution. However, there are varied associations between AOD/AOT and particulate matters (the main components of dust), depending on the size of particulate matters and spatial locations. While fine particulate matter is a common pollutant contributing to health risk in typical non-dust scenarios , there are also studies finding that PM10 or PM10–2.5 concentration may severely increase the mortality risk during a dusty day [16, 31, 40]. Therefore, further investigation is needed on how to use AOD/AOT to map spatiotemporal variations of both fine and coarse particulate matters for determining how air pollution exposures actually influence mortality risk during a dusty day. A future study should, therefore, combine an existing city-based mapping method with satellite images to improve the spatio-temporal modelling of air pollution exposure, at the same time increasing the spatial quality, because moderate or low spatial resolution of satellite images itself may increase potential spatial biases of modelling.
This study applied a spatial regression approach to estimate the spatial variability of mortality during prolonged dust events. The results indicated that spatial difference in built environment (SVF) and socioeconomic status (low education and low income) will increase the mortality in a community during dust events. This study also demonstrates there is a need to include spatial autocorrelation in modelling in order to improve the accuracy of prediction. Finally, the mortality risk map can be used to locate at-risk communities and vulnerable populations for developing health protocols in Hong Kong.
frontal area index (FAI)
planar area index (PAI)
land surface temperature
Akaike information criterion
sky view factor
tertiary planning unit
Vancouver Area Neighborhood Deprivation Index
aerosol optical depth
aerosol optical thickness
Aminipouri M, Knudby A, Ho HC. Using multiple disparate data sources to map heat vulnerability: Vancouver case study. Can Geogr. 2016;60(3):356–68.
Alexeeff SE, Schwartz J, Kloog I, Chudnovsky A, Koutrakis P, Coull BA. Consequences of kriging and land use regression for PM2.5 predictions in epidemiologic analyses: insights into spatial variability using high-resolution satellite data. J Eposure Sci Environ Epidemiol. 2015;25(2):138–44.
Anselin L, Rey S. Properties of tests for spatial dependence in linear regression models. Geogr Anal. 1991;23(2):112–31.
Anselin L, Sridharan S, Gholston S. Using exploratory spatial data analysis to leverage social indicator databases: the discovery of interesting patterns. Soc Indic Res. 2007;82(2):287–309.
Bao J, Li X, Yu C. The construction and validation of the heat vulnerability index, a review. Int J Environ Res Public Health. 2015;12(7):7220–34.
Barnett AG, Fraser JF, Munck L. The effects of the 2009 dust storm on emergency admissions to a hospital in Brisbane, Australia. Int J Biometeorol. 2012;56(4):719–26.
Bell N, Hayes MV. The Vancouver Area Neighbourhood Deprivation Index (Vandix): a census-based tool for assessing small-area variations in health status. Can J Public Health. 2012;103:S28–32.
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(6):kwt090.
Brown, M. J., Grimmond, S., & Ratti, C. (2001, December). Comparison of methodologies for computing sky view factor in urban environments. In International Society of Environmental Hydraulics Conference, Tempe, AZ.
Burnham KP, Anderson DR. Multimodel inference understanding AIC and BIC in model selection. Sociol Methods Res. 2004;33(2):261–304.
Carrier M, Apparicio P, Kestens Y, Séguin AM, Pham H, Crouse D, Siemiatycki J. Application of a Global Environmental Equity Index in Montreal: diagnostic and further implications. Ann Am Assoc Geogr. 2016;106(6):1268–85.
Chan CC, Ng HC. A case-crossover analysis of Asian dust storms and mortality in the downwind areas using 14-year data in Taipei. Sci Total Environ. 2011;410:47–52.
Chan EYY, Goggins WB, Kim JJ, Griffiths SM. A study of intracity variation of temperature-related mortality and socioeconomic status among the Chinese population in Hong Kong. J Epidemiol Community Health. 2012;66(4):322–7.
Chen YS, Sheen PC, Chen ER, Liu YK, Wu TN, Yang CY. Effects of Asian dust storm events on daily mortality in Taipei, Taiwan. Environ Res. 2004;95(2):151–5.
Chen L, Ng E, An X, Ren C, Lee M, Wang U, He Z. Sky view factor analysis of street canyons and its implications for daytime intra-urban air temperature differentials in high-rise, high-density urban areas of Hong Kong: a GIS-based simulation approach. Int J Climatol. 2012;32(1):121–36.
Crooks JL, Cascio WE, Percy MS, Reyes J, Neas LM, Hilborn ED. The Association between dust storms and daily non-accidental mortality in the United States, 1993–2005. Environ Health Perspect. 2016;124(11):1735–43.
Cutter SL. The vulnerability of science and the science of vulnerability. Ann Assoc Am Geogr. 2003;93(1):1–12.
Eeftens M, Beekhuizen J, Beelen R, Wang M, Vermeulen R, Brunekreef B, Hoek G. Quantifying urban street configuration for improvements in air pollution models. Atmos Environ. 2013;72:1–9.
Ebi KL, Kovats RS, Menne B. An approach for assessing human health vulnerability and public health interventions to adapt to climate change. Environ Health Perspect. 2006;114:1930–4.
Evans J, van Donkelaar A, Martin RV, Burnett R, Rainham DG, Birkett NJ, Grimmond S. Urbanization and global environmental change: local effects of urban warming. Geogr J. 2007;173(1):83–8.
Gal T, Lindberg F, Unger J. Computing continuous sky view factors using 3D urban raster and vector databases: comparison and application to urban climate. Theoret Appl Climatol. 2009;95(1–2):111–23.
Grimmond S. Urbanization and global environmental change: local effects of urban warming. Geogr J. 2007;173(1):83–8.
Hattis D, Ogneva-Himmelberger Y, Ratick S. The spatial variability of heat-related mortality in Massachusetts. Appl Geogr. 2012;33:45–52.
Health Canada. Adapting to extreme heat events: guidelines for assessing health vulnerability. 2011. https://www.canada.ca/en/health-canada/services/environmental-workplace-health/reportspublications/climate-change-health/adapting-extreme-heat-events-guidelines-assessing-healthvulnerability-health-canada-2011.html. Accessed 15 Jan 2017.
Ho HC, Knudby A, Sirovyak P, Xu Y, Hodul M, Henderson SB. Mapping maximum urban air temperature on hot summer days. Remote Sens Environ. 2014;154:38–45.
Ho HC, Knudby A, Huang W. A spatial framework to map heat health risks at multiple scales. Int J Environ Res Public Health. 2015;12(12):16110–23.
Ho HC, Knudby A, Walker BB, Henderson SB. Delineation of spatial variability in the temperature-mortality relationship on extremely hot days in greater Vancouver, Canada. Environ Health Perspect. 2017;125(1):66–75.
Hodul M, Knudby A, Ho HC. Estimation of continuous urban sky view factor from landsat data using shadow detection. Remote Sens. 2016;8(7):568.
Hondula DM, Davis RE, Leisten MJ, Saha MV, Veazey LM, Wegner CR. Fine-scale spatial variability of heat-related mortality in Philadelphia County, USA, from 1983–2008: a case-series analysis. Environ Health. 2012;11(1):1.
Jerrett M, Burnett RT, Beckerman BS, Turner MC, Krewski D, Thurston G, Gapstur SM. Spatial analysis of air pollution and mortality in California. Am J Respir Crit Care Med. 2013;188(5):593–9.
Johnston F, Hanigan I, Henderson S, Morgan G, Bowman D. Extreme air pollution events from bushfires and dust storms and their association with mortality in Sydney, Australia 1994–2007. Environ Res. 2011;111(6):811–6.
Kan H, London SJ, Chen G, Zhang Y, Song G, Zhao N, Chen B. Season, sex, age, and education as modifiers of the effects of outdoor air pollution on daily mortality in Shanghai, China: the Public Health and Air Pollution in Asia (PAPA) Study. Environ Health Perspect. 2008;116(9):1183.
Kandt J, Chang SS, Yip P, Burdett R. The spatial pattern of premature mortality in Hong Kong: How does it relate to public housing? Urban Stud. 2016;0042098015620341.
Kashima S, Yorifuji T, Tsuda T, Eboshida A. Asian dust and daily all-cause or cause-specific mortality in western Japan. Occup Environ Med. 2012;69(12):908–15.
Lai PC, Choi CC, Wong PP, Thach TQ, Wong MS, Cheng W, Wong CM. Spatial analytical methods for deriving a historical map of physiological equivalent temperature of Hong Kong. Build Environ. 2016;99:22–8.
Laaidi K, Zeghnoun A, Dousset B, Bretin P, Vandentorren S, Giraudet E, Beaudeau P. The impact of heat islands on mortality in Paris during the August 2003 heat wave. Environ Health Perspect. 2012;120(2):254.
Lorant V, Thomas I, Deliège D, Tonglet R. Deprivation and mortality: the implications of spatial autocorrelation for health resources allocation. Soc Sci Med. 2001;53(12):1711–9.
Morss RE, Wilhelmi OV, Meehl GA, Dilling L. Improving societal outcomes of extreme weather in a changing climate: an integrated perspective. Annu Rev Environ Resour. 2011;36(1):1.
Mueller MD, Wagner M, Barmpadimos I, Hueglin C. Two-week NO 2 maps for the City of Zurich, Switzerland, derived by statistical modelling utilizing data from a routine passive diffusion sampler network. Atmos Environ. 2015;106:1–10.
Neophytou AM, Yiallouros P, Coull BA, Kleanthous S, Pavlou P, Pashiardis S, Laden F. Particulate matter concentrations during desert dust outbreaks and daily mortality in Nicosia, Cyprus. J Eposure Sci Environ Epidemiol. 2013;23(3):275–80.
Pearce JR, Richardson EA, Mitchell RJ, Shortt NK. Environmental justice and health: a study of multiple environmental deprivation and geographical inequalities in health in New Zealand. Soc Sci Med. 2011;73(3):410–20.
Qi X, Hu W, Mengersen K, Tong S. Socio-environmental drivers and suicide in Australia: Bayesian spatial analysis. BMC Public Health. 2014;14(1):1.
Rey G, Jougla E, Fouillet A, Hémon D. Ecological association between a deprivation index and mortality in France over the period 1997–2001: variations with spatial scale, degree of urbanicity, age, gender and cause of death. BMC Public Health. 2009;9(1):1.
Rinner C, Patychuk D, Bassil K, Nasr S, Gower S, Campbell M. The role of maps in neighborhood-level heat vulnerability assessment for the city of Toronto. Cartogr Geogr Inf Sci. 2010;37(1):31–44.
Scarano M, Sobrino JA. On the relationship between the sky view factor and the land surface temperature derived by Landsat-8 images in Bari, Italy. Int J Remote Sens. 2015;36(19–20):4820–35.
Shi Y, Lau KKL, Ng E. Developing street-level PM2.5 and PM10 land use regression models in high-density Hong Kong with urban morphological factors. Environ. Sci. Technol. 2016;50(15):8178–87.
Shrestha R, Flacke J, Martinez J, Van Maarseveen M. Environmental health related socio-spatial inequalities: identifying “hotspots” of environmental burdens and social vulnerability. Int J Environ Res Public Health. 2016;13(7):691.
Smargiassi A, Goldberg MS, Plante C, Fournier M, Baudouin Y, Kosatsky T. Variation of daily warm season mortality as a function of micro-urban heat islands. J Epidemiol Commun Health. 2009;63(8):659–64.
Tobías A, Pérez L, Díaz J, Linares C, Pey J, Alastruey A, Querol X. Short-term effects of particulate matter on total mortality during Saharan dust outbreaks: a case-crossover analysis in Madrid (Spain). Sci Total Environ. 2011;412:386–9.
Tomlinson CJ, Chapman L, Thornes JE, Baker CJ. Including the urban heat island in spatial heat health risk assessment strategies: a case study for Birmingham, UK. Int J Health Geogr. 2011;10(1):1.
Tsai PJ, Lin ML, Chu CM, Perng CH. Spatial autocorrelation analysis of health care hotspots in Taiwan in 2006. BMC Public Health. 2009;9(1):1.
Wong CM, Ou CQ, Chan KP, Chau YK, Thach TQ, Yang L, Hedley AJ. The effects of air pollution on mortality in socially deprived urban areas in Hong Kong, China. Environ Health Perspect. 2008;116(9):1189.
Wong MS, Nichol JE, Holben B. Desert dust aerosols observed in a humid tropical city: Hong Kong. Int J Remote Sens. 2010;31(4):1043–51.
Wong MS, Xiao F, Nichol JE, Fung J, Kim J, Campbell J, Chan PW. A multi-scale hybrid neural network retrieval model for dust storm detection, a study in Asia. Atmos Res. 2015;158–159:89–106.
Wong MS, Peng F, Zou B, Shi WZ, Wilson GJ. Spatially analyzing the inequity of the Hong Kong urban heat island by socio-demographic characteristics. Int J Environ Res Public Health. 2016;13(3):317.
Yang J, Wong MS, Menenti M, Nichol J. Modeling the effective emissivity of the urban canopy using sky view factor. ISPRS J Photogramm Remote Sens. 2015;105:211–9.
Yang J, Wong MS, Menenti M, Nichol J, Voogt J, Krayenhoff ES, Chan PW. Development of an improved urban emissivity model based on sky view factor for retrieving effective emissivity and surface temperature over urban areas. ISPRS J Photogramm Remote Sens. 2016;122:30–40.
Yuan C, Ng E, Norford LK. Improving air quality in high-density cities by understanding the relationship between air pollutant dispersion and urban morphologies. Build Environ. 2014;71:245–58.
HCH developed the model design. MSW, HCH and TC conducted the data analysis. YL provided input of public health analysis. JY supported the data preparation. WS provided input of spatial analysis. MSW and HCH wrote the manuscript. All authors read and approved the final manuscript.
The authors are grateful for the mortality data from the Hong Kong Census and Statistics Department, airborne LIDAR data from the Hong Kong Civil Engineering and Development Department, land use and land cover data from the Planning Department of Hong Kong, and building GIS data from the Hong Kong Lands Department.
The authors declare that they have no competing interests.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
This research was supported in part by a grant from the Hong Kong Polytechnic University (Grant PolyU 1-ZE24). Dr. M.S. Wong acknowledges the support in part by a grant from General Research Fund (Project ID: 15205515); Grant PolyU 1-ZVAJ from the Faculty of Construction and Environment, the Hong Kong Polytechnic University; and the Grants PolyU 1-ZVFD from the Research Institute for Sustainable Urban Development, the Hong Kong Polytechnic University.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.