Abstract

Particulate matter with an aerodynamic diameter of 10 μm or less (PM10) pollution poses a considerable threat to human health, and the first step in quantifying health impacts of human exposure to PM10 pollution is exposure assessment. Population-weighted exposure level (PWEL) estimation is one of the methods that provide a more refined exposure assessment as it includes the spatiotemporal distribution of the population into the pollution concentration estimation. This study assessed the population weighting effects on the estimated PM10 concentrations in Malaysia for years 2000, 2008, and 2013. Estimated PM10 annual mean concentrations with a spatial resolution of 5 kilometres retrieved from satellite data and population count obtained from the Gridded Population of the World version 4 (GPWv4) from the Centre for International Earth Science Information Network (CIESIN) were overlaid to generate the PWEL of PM10 for each state. The calculated PWEL of PM10 concentrations were then classified based on the World Health Organization (WHO) and the national Air Quality Guidelines (AQG) and interim targets (IT) for comparison. Results revealed that the annual mean PM10 concentrations in Malaysia ranged from 31 to 73 µg/m3 but became generally lower, ranging from 20 to 72 µg/m3 after population weighting, suggesting that the PM10 population exposure in Malaysia might have been overestimated. PWEL of PM10 distribution showed that the majority of the population lived in areas that complied with the national AQG, but were vulnerable to exposure level 3 according to the WHO AQG and IT, indicating that the population was nevertheless potentially exposed to significant health effects from long-term exposure to PM10 pollution.

1. Introduction

Air pollution is a major environmental health threat that poses a complex challenge to public health. Despite various efforts taken to reduce its health impacts, it remains a huge problem due to rapid urbanization and growing population along with economic development. In 2012, an estimated 3.7 million premature deaths worldwide was attributable to air pollution, mainly caused by ischaemic heart disease and stroke, and the remaining 20% of which by respiratory disease [1]. According to the World Health Organization (WHO), Southeast Asia has more than 5% increase of particulate matter (PM) levels between 2008 and 2013. More than 50% of cities and populations in Southeast Asia are experiencing increasing trend of annual PM levels and mostly exceeding the WHO Air Quality Guidelines (AQG) [2].

In Malaysia, notable air pollution sources are from land transportations, industrial emissions, and open burning activities. 70–75% of air pollution in the country are from motor vehicle emissions, while stationary sources of air pollution such as power stations and industrial activities contribute to 20–25% [3]. Malaysia has been experiencing serious haze events since 1994 with more frequent trends over the past few years as reported in the chronology of haze episodes in Malaysia [4]. While open burning of solid wastes and forest fires contributed to only 3–5% of the country’s air pollution, it is notable that significant amounts of the air pollution were result of uncontrolled open burning activities from a neighbouring country [3].

Ambient air pollution is a complex mixture of multicomponents. PM and gaseous pollutants such as sulphur dioxide, nitric dioxide, and ozone are among the various components that may affect human health in the short and long term [5, 6]. Among the various components, PM is one of the air pollutants with considerable concern. Studies have found that PM pollution may contribute to adverse pregnancy outcomes [7, 8], type 2 diabetes [9], dementia [1012], and reduction of life expectancy [13]. And among the different size classes of PM, PM2.5 (particulate matter with an aerodynamic diameter less than 2.5 μm) is most consistently associated with adverse health effects [8, 9, 1315]. However, PM2.5 measurement data are presently limited in Malaysia; therefore, analyses in this study were focused on PM10 (particulate matter with an aerodynamic diameter less than 10 μm) with larger data availability as opposed to PM2.5.

Most PM10 exposure assessment studies have used concentrations of pollutants from monitoring stations and assumed that measured concentrations are representative of the mean exposure of the populations [1618]. However, directly using the ambient air concentrations to directly assess population exposure, without taking into account the disproportionate spatial and temporal distribution of the pollutant and the population, might not give the actual human exposure levels [19]. In addition, studies have found that the actual human exposure levels were in fact different from the mean value due to atmospheric pollutant distribution variations in time and space [20, 21].

Therefore, the aim of this study was to detect the effect of population weighting on PM10 concentration estimations. With the spatiotemporal population and pollution distribution variabilities added to the PM10 concentration estimations, the level of the pollutant concentrations exposed to the population would be different than the mean PM10 concentrations. Population weighting could give a better view of human exposure levels leading to better quantitative assessments which will then provide a stronger basis for critical public health evaluation on the extent of damage caused by PM10 pollution.

2. Materials and Methods

2.1. Site Description

Dataset used in this study covered the whole Malaysia for years 2000, 2008, and 2013. Malaysia is a federation consisting of thirteen states and three federal territories [22, 23]. Eleven states and two federal territories (Kuala Lumpur and Putrajaya) are located in the West of Malaysia, which are further classed into 4 regions (Central region, Northern region, Southern region, and East Coast region) in this study, while the remaining two states (Sarawak and Sabah) and one federal territory (Labuan) are classed as East of Malaysia. The three federal territories hereinafter will be referred to as states for simplicity. The 5 regions classified in this study are illustrated in Figure 1.

Putrajaya was afforded Federal Territory status of becoming Malaysia’s third Federal Territory in 2001 [23, 24]. The city was developed for the shift of government seat and relocation of federal government headquarters from Kuala Lumpur to Putrajaya because of overpopulation in the former [25, 26]. Prior to the declaration, Putrajaya was a territory entirely enclaved within Sepang, a district in the southern part of Selangor [25]. Hence, Putrajaya’s PM10 estimate and population count for year 2000 was based on PM10 and population data within its geographical coordinates.

2.2. PM10 Concentration Estimation

The data on aerosol products of Moderate Resolution Imaging Spectroradiometer (MODIS) Terra was obtained from the National Aeronautics and Space Administration’s (NASA) Level-1 and Atmosphere Archive and Distribution System (LAADS) Distributed Active Archive Centre (DAAC) website [27]. The aerosol products for MOD04, MOD07, and MOD021 were downloaded for the three studied years at a spatial resolution of 5 km. Data on Aerosol Optical Depths (AOD), surface temperature (ST), atmospheric stability (KI), and relative humidity (RH) were then extracted from the aerosol products and projected to the World Geodetic System (WGS) 84 coordinates using Environment for Visualizing Images (ENVI) software version 5.1. Using the Artificial Neural Network (ANN) model from a local study [28], mean annual PM10 concentration estimations were then calculated from the projected outputs in ArcGIS software version 10.1 as follows:where H are hidden layers, and

The calculated PM10 concentration estimations were then spatially interpolated using kriging spatial interpolation to fill in the empty pixels, resampled at a resolution of 0.05° × 0.05°, and validated with ground-based PM10 measurements obtained from the Department of Environment Malaysia (DOE). The results obtained in this study were satisfactory (R2 = 0.6225) between MODIS and the corresponding ground PM10 (Figure 2).

2.3. Population Distribution Data Collection

Population count was obtained from the Gridded Population of the World version 4 (GPWv4) from the Centre for International Earth Science Information Network (CIESIN). GPWv4 is the latest update of the GPW dataset and it demonstrates the spatial relationship of human population and the environment through gridded data products. Earlier versions of GPW have been extensively used in global population studies [27, 2931]. Gridded population data with a spatial resolution of 0.05° × 0.05°, matching that of PM10 concentration dataset, were downloaded for the years 2000, 2005, and 2010 (Figure 2). Since the GPWv4 provides gridded human population estimations in 5-year intervals starting from year 2000, the gridded population estimates for years 2008 and 2013 in this study were calculated in ArcGIS software using the exponential population projection formula as follows:where Px is the population estimate in the target year x, Py is the base population in year y, r is the average annual Malaysian population growth rate (2.0%) [32], and t is the number of years between population counts.

2.4. Population-Weighted Exposure Level to PM10 Calculation

Using ArcGIS software, the PM10 concentration layers were overlaid on the gridded population distribution layers and population-weighted exposure level (PWEL) of PM10 for each state were calculated based on the exposure equation as follows:where Pi is the population in grid i and Ci is its mean annual PM10 concentration.

2.5. Population Exposure Levels Categorization

Using the WHO interim target (IT) guideline as the basis [33], five exposure levels were established in this study to categorize the percentage of area and population exposed to different ranges of PWEL of PM10 concentration (Table 1). The IT guideline was proposed by the WHO with the intention to reduce mortality risks from exposures to air pollutants by promoting progressive pollutant emission control. Exposure level 5 is when PWEL is greater than 70 µg/m3 of PM10, suggesting that the population is at more than 15% higher risk for mortality. Lower exposure levels mean more reduction in mortality risks towards the population. Ultimately, exposure level 1 is when the PWEL is less than or equal to 20 µg/m3 of PM10, meaning that the state has achieved the WHO’s AQG for mean annual PM10 concentrations. Using ArcGIS software, population exposure distribution maps were then created for visualization of the different PWEL of PM10 vulnerability levels in each state during the study period.

Three IT values have been set in the New Malaysia Ambient Air Quality Standard, with an aim of improving air quality in stages and achieving the country’s air quality standard by 2020 [34]. With reference to the national standard and IT values, four exposure levels were established in this study to categorize the percentage of area and population exposed to different ranges of PWEL of PM10 concentration (Table 2). Between them, level 1 exposure indicates the percentage of land and population that were occupying in areas within the national PM10 AQG. Population exposure distribution maps for the study period were then created using ArcGIS software to illustrate the different PWEL of PM10 vulnerability levels in each state.

2.6. Population Density and Human Exposure Level to PM10 Correlation Analysis

Statistical analysis was done to measure the strength and direction of association that existed between population density and PWEL of PM10. Inspection of scatter plots between the two variables showed heteroscedasticity. Therefore, Spearman’s correlation test was run on the data.

3. Results

3.1. Malaysian Population Distribution

Table 3 represents the total geographical area, demographic details, and the population densities of each state over the study period. The spatial distribution of population was uneven in each state and population distribution trend was similar for the three years. Population densities were the highest in the Central region, lower in the East Coast region, and the lowest in East of Malaysia. Over the years, Kuala Lumpur remained as the state with the highest population density limited to only 0.07% of the total area for Malaysia, contrary to Sarawak which has the lowest population density despite having the largest area ratio of 37.65%.

3.2. PM10 Pollution in Malaysia

Table 4 shows the PM10 concentrations in Malaysia and its states before and after population weighting. The mean annual pollution concentration in 2013 for the whole country was seen to be increased as compared to 2000. However, lower pollution concentration estimates were observed in 2008 as compared to 2000 and 2008. In general, PM10 pollution concentrations were found higher in densely populated states, namely, Kuala Lumpur, Putrajaya, Selangor, and Penang.

3.3. Population Exposure to Estimated PM10

During the study period, the annual mean PM10 concentration and population-weighted exposure levels of PM10 in the states ranged from 31 to 73 µg/m³, whereas the latter showed concentrations from 20 to 72 µg/m³ (Table 4). With population weighting, PM10 was generally lower than the mean concentration in most states, except in Kuala Lumpur in 2000; Putrajaya, Pahang, and Terengganu in 2008; and Selangor in 2013 where PWEL of PM10 were slightly higher than the mean concentration.

3.4. Population and Area Ratio in Different Exposure Levels

Table 5 and Figure 3 show the percentage of area and population exposed to different ranges of PWEL of PM10 concentration with reference to the WHO IT and AQG. During the key period, 65.40–73.38% of the Malaysian population were living in level 3 exposure areas where PWEL ranged from 31 to 50 µg/m³ of PM10. In year 2000, only Penang which accounted for 0.31% of the land occupied by 5.05% of the total Malaysian population complied with the WHO guideline for annual mean PM10 concentration, whereas in 2008 and 2013, none of the states achieved the WHO recommended mean annual PM10 concentration levels. None of the areas were seen to reach level 5 exposure with PM10 concentration exceeding 70 µg/m³ in 2000 and 2008. However, population exposure to the pollutant worsened in 2013, where the whole country was above exposure level 2, with Putrajaya accounting for 0.01% of the land and 0.25% of the Malaysian population were vulnerable to level 5 exposure.

Table 6 and Figure 4 show the percentage of area and population exposed to different ranges of PWEL of PM10 concentration with reference to the national PM10 IT and AQG. A majority of the population were seen living in areas that complied with the national AQG where PM10 concentrations were less than or equal to 40 µg/m³. However, the ratio of area and population living in the national recommended mean annual PM10 concentration levels decreased in 2013. Less land and less people resided in exposure level 1 and 2 areas, while more land and more people resided in exposure level 3 and 4 areas were observed in 2013 as compared to 2000.

3.5. Correlation between Population Density and Human Exposure

Results showed that there was low positive correlation between population density and PWEL of PM10, r = 0.442, p = 0.002 (Table 7), suggesting that the two were marginally correlated. The data reflects that as population density increases in value, so does the PWEL of PM10 to some extent. However, the amount of increase was not constant over the whole range of values.

4. Discussion

Malaysia is a rapidly developing and newly industrialized country. Industrial activities have polluted the air with smoke from vehicles, factories, and open burning activities; hence, an increasing trend of PM10 was assumed as years go by from 2000 to 2013. However, this study revealed an overall reduction of PM10 concentrations in 2008. Upon further investigation, it was found that the lower PM10 concentrations in between the key period could be partially attributable to local meteorological conditions and confounding factors such as precipitation, humidity, and economy.

Studies have found that influence of climate, particularly rainfall, contributes to lower PM10 concentrations through washing effect [3538]. It should be noted that the mean annual rainfall in Malaysia was higher in 2008 as compared to 2000 and 2013, where the average yearly precipitation recorded for the country was 2393 mm [39], 10535 mm, and 9109 mm [40, 41] for 2000, 2008, and 2013, respectively. Therefore, the higher precipitation level could have reduced the PM10 concentrations in 2008. Additionally, the 24-hour mean RH that was extracted and calculated from the downloaded aerosol products revealed that the average RH was 62%, 67%, and 65% for years 2000, 2008, and 2013, respectively [42]. Due to the fact that particulate matter tends to accumulate and deposit on the ground rather than becoming airborne when ambient relative humidity is high, the higher RH index as a consequent of higher rainfall in 2008 could also explain the lower PM10 concentration levels in that year. These findings are consistent with several studies that have assessed the linkages between meteorological processes and air pollution and have found that humidity increases as rainfall increases which then subsequently decreases the PM10 concentration [37, 38].

Apart from climate, economy may also influence the variation of air pollution concentrations of a country in several ways. It is noteworthy that Malaysia was affected by the Global Financial Crisis in 2008 after a decade of steady economic growth [43]. Censuses have reported that the country had faced declinations in construction, manufacturing, mining, and quarrying businesses in 2008 in line with the deterioration in external demand [44, 45]. Malaysia’s economy started to recover the next year and in 2013 the government announced that the country’s gross national income per capita then had increased by 50% from 2009 [43]. Similarly, previous studies have also remarked reductions of particulate pollution in relation to economic recessions [46], followed by reversion of the pollution levels after the economy rebounds [3, 47, 48].

One of the variables in the calculation of the population exposure is the distribution of population density, i.e., population in grid i (Pi). It is dependent on the dynamics of local populations which could be seen highly influenced by urbanisation. In Malaysia, significant urban agglomeration is observed more towards the central region with an urbanisation level of 100% in Kuala Lumpur and Putrajaya and 91.4% in Selangor. Apart from the three states, another state with high level of urbanisation is Penang in the northern region (90.8%) [49]. When calculating PWEL of PM10 for each state, it is expected that the exposure increases as the population density increases. Other exposure studies have reported higher PM10 concentrations after population weighting, attributing it to high population density [31, 5052]. However, results from this study revealed the opposite. Overall, lower PWEL of PM10 values were observed in this study indicating that most people were probably exposed to lower PM10 concentrations than the actual mean. Despite the contrary, it should be noted that the other studies were done in China and India where the population densities are very much higher as compared to Malaysia. Other than the diverse population and pollution spatial distributions, the massive population density variation between the countries might have also contributed to the differences in postweighting values.

In 2000, Penang unexpectedly recorded the lowest PWEL of PM10 concentration for the whole study duration. With a mean annual PM10 concentration of 42 µg/m³ and referring to the WHO IT and AQG guideline, Penang should have been categorized in exposure level 3. However, after population weighting, the estimated PM10 was lowered to 20 µg/m³, thereby classifying Penang in exposure level 1. Moreover, Penang is one of the country’s most urbanised states and accommodates among the nation’s highest population densities [49]; thus, higher exposure levels were presumed. Upon further investigation, it was discovered that the state’s economic growth slowed down in 2000 due to weakening of manufacturing dynamism [53]. Likewise, previous studies have shown that attenuation of industrial sources of particulate pollution led to pollution reduction [47, 48, 54]. The combination of the lower spatial distribution of PM10 on top of the high population density could explain the reduced PWEL of PM10 concentration estimations. Therefore, population density and urbanisation level of a region are not the only factors that may reflect PWEL of PM10 values, but rather a combination of demographic and local contributing factors.

Though PWEL of PM10 concentrations in Malaysia were overall lower as compared to the actual mean, results from this study have shown that the population were nevertheless exposed to worsening levels of PM10 concentration. This is regrettably expected that the more areas are being developed, the more people will migrate and reside in these regions susceptible to unhealthy levels of air quality, making them vulnerable to chronic exposure to the coarse pollutant. These findings are consistent with previous population exposure assessment studies which also noted that more people are living in areas prone to higher PM10 concentrations [50, 55, 56].

Industrial emission is one of the major sources of air pollution in Malaysia [57], which is an unfavourable outcome of the intense process of urbanisation and industrial development since this country aimed for rapid economic growth to obtain industrial country status by 2020. Due to this phenomenon, Malaysia has to contend with increasing air pollution and has set its own interim targets with regard to the national capability in air quality management. The cut-off value for PM10 pollution in the national AQG and IT guideline is higher than most developed countries because the country has yet to complete its industrialization programs.

To further address this issue, the DOE under The Ministry of Energy, Science, Technology, Environment and Climate Change Malaysia (MESTECC) has outlined national and international preventive measures. One of the efforts include enforcing laws that prohibit the public from engaging in open burning [58]. Malaysia has also signed the ASEAN Agreement on Transboundary Haze Pollution in 2002 as an endeavour to tackle the issue [59]. Ratification of the treaty is hoped to start a new dimension in addressing forest and land deforestation issues as well as transboundary haze concerns.

To improve air quality status for the betterment of people’s lives, the DOE have also outlined five strategies under the Pelan Tindakan Udara Bersih (Clean Air Action Plan) developed in 2011 [60]. The five strategies are (1) to reduce motor vehicle emissions, (2) reduce air pollution due to land and forest fires, (3) reduce industrial emissions, (4) develop self-efficacy and ability, and (5) strengthen community engagement and awareness. While it may take decades for the outcomes of these strategies to be visible, the efforts must nevertheless be continued. A great example is the Clean Air Act, the law that defines the United States Environmental Protection Agency’s (EPA) responsibilities for protecting and improving its air quality, where actions which have been taken to implement the statutes have reduced the key air pollutant emissions since 1990 by 50% [61]. Therefore, it is evident that persistent and consistent action can definitely lead to dramatic improvements of the air quality.

5. Limitation

This study has used satellite remote sensing observations for PM10 concentration estimations over Malaysia in view of limited number of air pollution monitoring stations coverage for the whole country. Thick cloud covers and satellite coverage gaps are known limitations when using satellite data and due to mid-latitude jet stream instability, extensive cloudiness is a common phenomenon in Southeast Asia [62]. This occurrence reduces the number of cloud-free images available for AOD retrievals [63] which could affect the estimated PM10 concentrations and subsequently influence the PWEL of PM10 values.

Moreover, this study only provides an overview on the level of Malaysian population exposure to coarse particulate matter. The average time spent outdoors by the population and source-orientated assessments should be accounted for more accurate inhaled concentrations and human exposure levels. Such investigations would enable more precise and comprehensive evaluations for better understanding of the true level of human exposure to air pollution.

6. Conclusion

The mean annual population exposure to PM10 in Malaysia for 2000, 2008, and 2013 achieved the national AQG but not the WHO standards, indicating that the majority of the Malaysian population lived in the probable risk zones which can induce significant health effects associated with long-term PM10 exposures. Geographic, demographic, and local confounding factors contributed to the dynamic spatiotemporal variations of both PM10 concentration and population distribution, which then affected the PWEL of PM10 estimation concentrations. Therefore, including population weighting into PM10 concentration estimations could improve human exposure assessment as it takes these variables into account. The real situation would hence be better reflected and policy makers will be able to properly evaluate and establish appropriate control measures to achieve the lowest possible risk of population exposure to air pollution.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors have no conflicts of interest with the materials presented in this paper.

Acknowledgments

This study was financially supported by the Ministry of Health Malaysia Research Grant (MRG) under the National Medical Research Register code NMRR-16-2742-33757. The authors would like to thank the Director General of Health Malaysia for his permission to publish this article and the DOE for providing the PM10 concentrations data. The authors also acknowledge NASA for making MODIS data available to users. The authors also thank Dr. Rafiza Shaharudin for helping with the paper editing.