Special Session on Better Air Quality in Asia (II)

Gita Pati Humairoh1, Arie Dipareza Syafei This email address is being protected from spambots. You need JavaScript enabled to view it.1, Muhayatun Santoso2, Rachmat Boedisantoso1, Abdu Fadli Assomadi1, Joni Hermana1

1 Departement of Environment Engineering, Institut Teknologi Sepuluh Nopember, Sukolilo, Surabaya 60111, Indonesia
2 Center for Applied Nuclear Science and Technology, National Energy Agency of Indonesia

Received: November 14, 2019
Revised: May 12, 2020
Accepted: May 14, 2020

 Copyright The Author(s). This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are cited.

Download Citation: ||https://doi.org/10.4209/aaqr.2019.11.0590  

  • Download: PDF

Cite this article:

Humairoh, G.P., Syafei, A.D., Santoso, M., Boedisantoso, R., Assomadi, A.F. and Hermana, J. (2020). Identification of Trace Element in Ambient Air Case Study: Industrial Estate in Waru, Sidoarjo, East Java. Aerosol Air Qual. Res. 20: 1910–1921. https://doi.org/10.4209/aaqr.2019.11.0590


  • Atmospheric heavy metals come from a mix of anthropogenic and natural activities.
  • Smelting activities contribute significantly towards Pb contamination.
  • Pollutants from vehicles make up a total composition of trace elements in ambient.
  • Particulate matter in near ocean area is contributed by sea salt presence.


Metal is one of the air pollutants found in air particulates. The presence of heavy metals in air can be due to human activities or natural factors. Heavy metals can affect human health, causing respiratory disease and even death. The purpose of this study was to determine daily particulate matter (PM) concentrations in ambient air at the Waru Industrial Estate, analyze the results, and then characterize and estimate the locations of pollutant sources. PM was collected for 24-hour periods with gent stacked filter units. Filters were analyzed via X-ray fluorescence (XRF) to find concentrations of metal particles. The measurement data were analyzed via principal component analysis (PCA) and the conditional probability function (CPF) method in order to identify and estimate the industrial pollutant sources that contribute to these metal particles being in the ambient air. Results arrange PM2.5 concentrations from 2.65 to 32.68 µg m–3, with an average daily concentration of 17.67 ± 7.29 µg m–3, whereas PM10 concentrations ranged from 14.69 to 72.27 µg m–3, with an average daily concentration of 40.70 ± 13.78 µg m–3. The elements identified with XRF were Na, Mg, Al, Si, S, K, Ca, Ti, Cr, Mn, Fe, Cu, Zn, Pb, and Cl. The PCA results explain that there are four main components (factors) which then become a potential source of pollutants, namely the first is a marker of industrial activity. the second indicates the activity of the metal smelting industry. third shows the contribution of sea salt. Fourth shows the biomass combustion emissions. The results of the CPF method show that the first factor originates from north to northeast with a probability of 0.5. The second factor comes from west to southwest with a probability of 0.4. The third factor comes from northeast to east with a probability of 0.4–0.45.

Keywords: Particulate matter; Metal; Industry; PCA; CPF.


Pollution caused by metal particles in Indonesia has worsened due to increased economic activity, rapid population growth, urbanization, and industrialization. Metal is one of the air pollutants found in air particulates. The presence of metals in the air can be due to human activities or natural factors such as volcanoes, ocean aerosols, chemical reactions in the atmosphere, and so forth (Santoso and Lestiani, 2014). Some studies have shown that ±20 metal particles are present in the air, including Al, Ag, As, Br, Ca, Co, Cr, Cs, Fe, Hg, I, K, La, Mg, Mn, Na, Sb, Sc, Sm, Ti, Si, Ni, S, V, Sr, and Zn (Santoso and Lestiani, 2014). The metal particles are PM2.5, which is particulate matter containing fine particles with diameters of less than 2.5 µm (WHO, 2007). are emitted by either anthropogenic or natural sources (Deshmukh et al., 2010).

Fine particles that contain metals can affect human health if they are inhaled and penetrate the lungs. Particles can cause acute respiratory infections (ARIs), symptoms of anemia, barriers to growth, a weak immune system, symptoms of autism, lung cancer, and even premature death (Mukhtar et al., 2013). Schaumann et al. (2004) reported that fine particles containing metals could trigger the oxidant radical generation of brochoalveolar lavage cells. Cytokine concentrations in brochoalveolar lavage fluid were elevated after the introduction of PM2.5 into the Hettstedt area (Schaumann et al., 2004). Hassanvand et al. (2015) revealed that Si, Fe, Zn, Al, and Pb were the most abundant metal materials measured in a retirement home and school dormitory in Tehran, Iran. A multi-city study conducted in Europe (29 cities) and the United States (20 cities) reported a short-term mortality effect for PM10 of 0.62%, with a 0.46 percentage point increase for every 10 µg m–3 increase in PM10 concentration (24 hours, on average) (Samet et al., 2009). A Harvard study of 24 cities in the US indicated that long-term exposure to PM resulted in an increase in the occurrence of respiratory disease in children (Dockery et al., 1993).

Over the past several decades, multivariate analysis has been widely used to identify sources of particles in ambient air, one of which is the application of principal component analysis (PCA). PCA analysis was carried out in identifying sources of PM2.5 emissions which produced five factors of sources of PM2.5 emissions in the City of Anglet located in the southwest of France. The five factors are factor 1: soil dust, factor 2: combustion, factor 3: industry, factor 4: vehicle, factor 5: sea (Chavent et al., 2009). To estimate the location of local pollutant sources can be done by combining elemental concentration data with meteorological data of wind direction and speed, better known as the CPF (Conditional Probability Function) method. The CPF method was used by Pekney et al. (2006) in Pittsburgh, the results of his study showed that the source of PM2.5 pollutants is likely to come from the southeast with a probability of 0.70, in the southeast there is a coal-fired power plant, which allows the combustion of rock these embers can contribute to PM2.5.

PM2.5 concentration in Surabaya for 2012 to 2014 was 15.05 µg m–3, so PM2.5 had exceeded the quality standard yearly, based on PP 41 1999, U.S. EPA and WHO (Ahmad and Santoso, 2016). The results of astudy on PM containing heavy metals in the Serpong industrial area, which was characterized by fine particulates (PM2.5), identified five factors in the spread of Pb, namely, the smelting industry (9.61%), soil dust (17%), a mixture of metal smelting and sea salt industries (13.02%), transportation (44.36%), and burning of biomass (22.58%). PM quality studies in Surabaya found 12 metal elements in PM2.5 and PM2.5-10 including Na, Al, Si, S, K, Ca, Ti, Mn, Fe, Cu, Zn, and Pb (Santoso et al., 2016).

Increasing industrial development activities have the consequence of increasing waste generated by the industry, including air pollutants that can change ambient air quality (Mukono, 2011). Industrial activities are complex and involve various processes. The use of fuel, incineration processes or the combustion of raw materials with high temperatures are generally found in industrial activities (Zannaria et al, 2009). Emissions from industry are believed to be one of the main contributors to pollutant sources to air particulates (Santoso et al., 2013).

One of the main industrial zones in the south of Surabaya is in the District of Waru. Waru is also known as the buffer industry center of Surabaya, and many important industries were previously based in this sub-district. For example printing, jewelry making, nail factories, milk factories, cosmetics companies, biscuit companies to insulation factories and steel companies. The number of industries in the Waru area is feared to make the concentration of metal pollutants become large and have a negative impact on public health. A study conducted by Mukhtar et al, 2013 concerning the content of heavy metals in ambient air in several cities in Indonesia, that Surabaya identified elements of Pb, Na, Al, S, K, Fe, Mg, Mn, Si, Ca and Zn. The highest Pb concentration is found in Surabaya which is 2664 ng m–3, Pb can be generated from industrial emissions and motor vehicles. Besides Pb, in Surabaya also the highest S, Zn and Mn elements were identified with an average value of 400 ng m–3, 344.78 ng m–3 and 14.07 ng m–3. Based on the results of the above studies, a study on PM pollutant quality which focuses on particulate metal elements in one of the industrial zones in Waru, Sidoarjo, East Java is needed. This study will include the identification of PM elements as well as their characterizations and estimated sources


Sampling Method

The sampling was conducted in March–May 2019 around the Waru Industrial Estate, Sidoarjo, East Java, Indonesia. Waru is one of the main industrial zones in the south of Surabaya, precisely in the Waru District, Sidoarjo Regency, East Java Province, Indonesia. This district is bordered by Surabaya City, in this district there is the Purabaya Terminal, the largest bus terminal in Indonesia. The north side of the district is the Waru Roundabout, which is the main gate of Surabaya City from the southwest (Mojokerto / Madiun / Kediri / Solo / Yogyakarta) and from the south (Malang / Banyuwangi). This sub-district is a border area between South Sidoarjo and Surabaya and is an area that is experiencing rapid development. In addition to its strategic location, with a variety of potentials such as in the industrial sector from the metal industry and shoes/sandals, trade, as well as regional small and medium businesses as well as adequate human resource support, the Waru District is one of the strategic areas for development economy. Fig. 1 shows the location of the sampling site around the Waru Industrial Estate (7°20′43.64′′S, 112°45′21.55′′E). The gent stack had a place data height of 5.5 meters from the ground within 1–2 km of industries that pollute the air (SNI.19-7119.6-2005). The sample analysis was carried out at the Center for Applied Science and Technology: National Nuclear Energy Agency (PSTNT-BATAN) in Bandung.

Fig. 1. Location of the sampling site in the Waru Industrial Estate.
Fig. 1.
 Location of the sampling site in the Waru Industrial Estate.

The air particulate test sample was obtained by using the gent stack, which was able to capture PM2.5 and PM2.5-10. Samples were taken every 24 hours over a period of 35 days from March 14, 2019, to May 18, 2019. Overall, 35 samples were collected at the site.

The Gent SFU Sampler consists of two main parts namely the black receptacle and the pump system which regulates with the timer. Gent SFU Sampler is a dichotomous sampler sampler, which is a type of sampler that can collect air particles with a size of less than 2.5 µm and particulates with a size of 2.5–10 µm. Gent stack had an impactor-stage inlet-cut-off point system on the head of the outer part of its container that only allowed dust below 10 µm in diameter to enter. These dust particles first encountered coarse filter size of 2.5–10 µm, and those particles able to pass through the first filter then accumulated on a fine filter of less than 2.5 µm in size. The flow rate used in the gent stack was 18 L min–1 (Hopke et al., 1997). Wind speed and direction were assessed using a Kestrel 5500 at the same elevation as the gent stack. The average wind speed ranged from 0 to 6 m s–1 (with a mean of 1 m s–1), the average temperature ranged from 23.2°C to 37.7°C, and the relative humidity ranged from 48.3% to 100% (mean value: 81.47%).

All filters were conditioned for 24 hours in a clean room equipped with a dehumidifier and an air conditioner to control air humidity in the range of 40–50% and temperatures in the range of 18–25°C, respectively. Filters were then weighed using a Mettler Toledo micro balance with a precision of up to six digits because the samples weighed very little.

Each filter set was then arranged in a filter cassette with the coarse filter placed on the outside or at the top. The fine filter was placed on the inside of the filter cassette. After the cassette filters were filled, they were connected to the Gent SFU Sampler for sampling. Samples were collected using a Gent SFU capable of collecting particulate matter in the PM2.5-10 and PM2.5 size fractions (Hopke et al., 1997). Sampling was performed at a flow rate of 15–18 L min–1 for 24 hours over 3 month periods. There were a total of 35 pairs of PM2.5-10 and PM2.5.

Sample Analysis

The PM2.5 concentrations were determined via the gravimetric method, namely, by balancing the samples on the fine filters with a fine, smooth flange weight. For PM10, the concentrations were obtained by using both the fine and rough filters. X-ray fluorescence (XRF) spectroscopy was used to determine mass concentrations and identifications of the metal particle elements in the samples. The XRF analysis method is based on the emission of X-rays by atoms which are excited by high-energy photons such as X-rays and γ-rays (Santoso et al., 2010). XRF can detect various elements in air particles and produce a data set of 25–30 elements in the identification of sources of air pollutants (Santoso et al., 2014). The principle of XRF spectroscopic analysis is to utilize the emission of radiation emitted from secondary characteristics X-rays that characterize each element of a material. The secondary characteristic X-rays are produced by the process of excitation or ionization of the innermost skin electrons due to the emission of primary X-ray radiation emission at certain energies towards a target called the photoelectric effect. Secondary X-ray radiation is then detected by an X-ray detector into electrical pulses. The pulses are reinforced by the initial amplifier and the final amplifier. Electric pulses that have been strengthened by the final amplifier are used as input for the Analog to Digital Converter (ADC) to be converted into digital numbers and then processed by the computer into an area of the spectrum of analysis results (Rixson et al., 2015)

The XRF used for metal analysis contained in the PM was the XRF spectrometer Epsilon5 PANalytical. The XRF Epsilon5 is a type of XRF Energy Dispersive System (EDXRF) that can separate the characteristic radiation energy originating directly from the sample.

Statistical Analysis

Determination of pollutant sources was carried out with the help of principal component analysis (PCA) and conditional probability function (CPF). The purpose of PCA is to obtain the maximum number of components that can explain the variations in the data (Thurston and Spengler, 1989). The source profile is then estimated by finding a profile with the same characteristics as in the literature (Reff et al., 2007). The CPF determine pollutant sources based on plots of wind direction and speed.

PCA is part of a multivariate technique that can reduce dimensions in data groups by identifying new data from several independent variables that explain the covariance of the original data (Lau et al., 2012). Principal component analysis (PCA) and cluster analysis (CA) were carried out using the SPSS statistical software package version 13.0 for Windows. The PCA model is a statistic that identifies data patterns and shows differences and similarities in data. The main purpose of a PCA application is to obtain a number of component values that can explain the maximum data inheritance (Thurston and Spengler, 1989). The PCA model does not change, reduce, or add data but only replaces variables that are correlated with data that are not correlated through linear combinations of original variables. The first variable resulting from dimension reduction should explain almost all data variability. The result of PCA is the eigenvalues and eigenvector of the matrix (Henry and Hidy, 1979). The correlation matrix always has at least one non-zero eigenvalue, and all these non-zero values are positive. The matrix of the diagonal transformation is derived from the eigenvector as a column/row depending on the original vector of the variables. The principal component (PC) has the greatest association with the eigenvalue. This value is a linear combination of variables that produces the maximum value of the total data. These characteristics were seen through the element markers in each type of source. In general, pollutants can be characterized singly or in combination (Seinfeld and Pandis, 2006). After the elements in PM are identified, the XRF measurement results are then analyzed using PCA. PCA was carried out using SPSS 17.0 to identify sources that contribute to air pollution around the Waru Industrial Estate, Sidoarjo. The XRF measurement data is the standard for generating all the variables of identical variation and calculating the associated eigenvalues. PCs consisting of different sources were analyzed by applying orthogonal rotation varimax with an identical number of eigenvalues. The number of PCs is determined by extracting them with a PC eigenvalue greater than 1, the Kaiser criterion. It contains the original variable value (source of emissions) on each PC counted. Loading values greater than 0.5 are considered to be strong correlations between PCs and emission sources.

PCA is widely used to reduce data dimension in metal distribution studies, it extracts a small number of latent factors, referred to principal components to analyze relationships among observed variables (Yeung et al., 2003; Han et al., 2006; Meza-Figueroa et al., 2007). The concentrations of heavy metals investigated in this study vary by different orders of magnitude, and PCA assists the identification of pollutant sources. Despite that PCA suffers in terms of method in which its eigenvector analysis is based on unweighted least-square fit to the data, which may not provide best estimators for parameters of interest (Hopke and Jaffe, 2020), its usage remains important to provide first step of a more complex process, provided with other tool such as PMF, in quantifying the sources (Chavent et al., 2009).

The CPF determine pollutant sources based on plots of wind direction and speed. To identify possible sources of heavy metals based on the wind speed and direction, a CPF receptor model is used. The conditional probability function (CPF) seen in Eq. (1) calculates the probability that a source is located within a particular wind direction sector, ∆θ:

where n∆θ is the number of times that the wind passed through the direction sector ∆θ, and m∆θ is the number of times that the source contribution peaked while the wind passed through sector ∆θ (Ashbaugh et al., 1985). To use the CPF method, the 24-hour averaged source contribution data were applied to all 30-minute wind direction averages measured at the site on each date. All time periods with wind speeds less than 1 m s1 were removed from the data set. To calculate m∆θ, the highest 25% of the source contribution concentrations were used. CPF is useful in determining the direction of a source from a receptor site (Begum et al., 2010).

Estimation of the location of local pollutant sources is done by combining the polluted source mass concentration data with meteorological data on wind direction and speed hour period. The second analysis of the data uses the CPF method with the help of Microsoft Excel software. The CPF method is to divide the number of events in a particular wind direction with the largest contribution of factor mass concentration with wind direction and velocity above 1 m s–1 to the total factor mass concentration in the direction and wind speed greater than 1 m s–1.


Particulate Mass Concentrations

A total of 35 samples of PM2.5 and PM10 aerosols were collected from March 2019 to May 2019 at the Waru Industrial Estate sampling site. The concentrations of PM2.5 and PM10 ranged from 2.65 to 32.68 µg m–3 and from 14.69 to 72.27 µg m–3, respectively. Respective daily mean concentrations of 18.26 ± 6.40 µg m–3 and 43.59 ± 14.53 µg m–3 were recorded for PM2.5 and PM10.

When compared with the daily quality standard (24 hours) of 25 µg m–3 published by the WHO, there were three days for which the PM2.5 concentrations exceeded the quality standard, namely, on April 14 and 15,2019 and May 6, 2019, with concentrations of 28.44, 32.68, and 25.05 µg m–3, respectively. There were decreases in the concentration on several days, namely, on March 17, 20, 23 and April 3, 7, 16, 18, 28, 2019, to 6.00, 4.81, and 13.46 µg m–3 and 11.22, 13.91, 2.65, 7.78 and 11.40 µg m–3, respectively. For PM10, when its concentrations are compared with the daily standard (24 hours) of 50 µg m–3 published by the WHO, there were 12 days for which the concentrations exceeded the quality standard, namely, on March 14, 24, 28 2019, April 4, 14, 15, 2019, and May 4, 6, 8, 9, 11, 13, 2019, with PM10 concentrations of 66.53, 53.05, and 55.03, 56.06, 52.41, and 68.71 µg m–3; and 56.60, 65.65, 72.277, 53.74, 53.18, and 55.76 µg m–3, respectively. There were decreases in the concentration on several days, namely, on March 17, 20 and April16, 18, 19, 2019, of 21.74 and 30.83 µg m–3, and 14.69, 17.25, and 31.48 µg m–3, respectively. Based on the calculations of the PM10 concentrations during the study period around the Waru Industrial Estate, Sidoarjo, the PM10 concentration had an average concentration value of 43.59 µg m–3, with the lowest concentration value of 14.69 µg m–3 and the highest concentration value of 72.27 µg m–3. When compared with the Government Regulation No. 41/1999 concerning air pollution control quality standard for PM10 of 150 µg m–3 over a 24-hours period, the concentration of PM10 in the vicinity of the Waru Industrial Area, Sidoarjo, is still below the quality standard. This result is consistent with the research conducted by Ahmad and Santoso (2016), which showed that the average PM10 concentration in Surabaya during the period from October 2012 to February 2014 was 30.41 µg m–3, well below 1999’s 150 µg m–3, the U.S. EPA’s 150 µg m–3, and the WHO’s 50 µg m–3 (Ahmad and Santoso, 2016). This finding is also supported by Pitaloka and Adriyani (2016), who showed that the PM10 concentrations in the Ready Mix Concrete area of Waru, Sidoarjo amounting to 0.032 mg Nm–3 did not exceed the threshold value set by the Ministry of Man power and Transmigration’s Regulation No. PER.13/MEN/X/2011 of 0.26 mg Nm–3 over 24 hours.

The mean of the values of the PM2.5 and PM10 concentrations were compared with six cities (Table 1). The mean PM2.5 concentrations in Waru, Sidoarjo, Bandung, Jakarta, Palangkaraya, Serpong, and Yogyakarta were 18. 26, 18.35, 16.50, 7.74, 16.50, 7.74, 16.68, and 8.78 µg m–3, respectively, while the mean PM10 concentrations were 43.59, 35.81, 41.58, 17.58, 32.05, and 20.13 µg m–3, respectively. The PM2.5 concentrations for Waru, Sidoarjo, Bandung, Jakarta, and Serpong are higher compared with those for Palangkaraya and Yogyakarta, indicating that Waru, Sidoarjo, Bandung, Jakarta, and Serpong are more polluted than Palangkaraya and Yogyakarta. The PM10 concentrations in Waru, Sidoarjo, Jakarta, Bandung, and Serpong were also higher than they were for Yogyakarta and Palangkaraya. It should be noted that the amount of industrial activity in Yogyakarta and Palangkaraya was less than that in Waru, Sidoarjo, Jakarta, Bandung, and Serpong. The rapid growth of industry, the population, and motor vehicle usage in Waru, Sidoarjo, Bandung, Jakarta, and Serpong has affected their environmental quality considerably. The imbalance between the number of motor vehicles and available roads has created traffic problems in almost every section of these cities and resulted in air quality deterioration (Santoso et al., 2013).

Table 1. The mean of the PM2.5 and PM10 concentrations.

The decreasing and increasing patterns seen in the PM2.5 and PM10 concentrations around the Waru Industrial Estate occurred on almost the same dates, as can be seen in Fig. 2. The pattern could have been influenced by meteorological parameters at the time of the measurements. When the PM concentration is high, the wind speed ranges from 0.5 to 2.1 m s–1, and the dominant wind direction comes from the west with a temperature of 27°C–30°C and a humidity of 70%–85%, whereas when the PM concentration is low, the wind speed ranges between 2.1–3.6 m s–1 and the dominant wind direction comes from the east with a temperature of 27°C–29°C and a humidity of 70%–80%. These results are consistent with research conducted by Proias et al. (2010) that concluded that wind direction and speed correlate significantly with the average daily PM concentration. During the measurement period, it rained for several days, such as on April16, 2019, when it rained at a temperature of 27°C, and the PM concentration experienced a small decrease. Rain greatly affects the concentration of PM because it causes air borne particulates to float in the air bound to rain drops, which then descend toward the surface of the earth (Holst et al., 2008).

Fig. 2. Mass concentrations of PM2.5 and PM10.
Fig. 2.
 Mass concentrations of PM2.5 and PM10.

The correlation test for the PM2.5 and PM10 concentrations showed significantly high correlations between their concentrations. Detemination coefficient (R2) gives the meaning of the percentage effect of mass concentration PM2.5 on PM10. The relationship pattern between the two is the root of the value R2. In Fig. 3, the PM10 value of 67.50% is influenced by changes in the mass concentration variable PM2.5 while 32.50% is influenced by other variables, as shown in Fig. 3.

Fig 3. Correlations for PM2.5 concentrations with PM10 concentrations
Fig 3.
 Correlations for PM2.5 concentrations with PM10 concentrations

Multi-elemental Analysis of Particulates

he elements identified by XRF spectroscopy for the PM2.5 fine particle fractions and PM2.5-10 coarse particle fractions were Na, Mg, Al, Si, S, K, Ca, Ti, Cr, Mn, Fe, Cu, Zn, Pb, and Cl. The box and whisker plots in Fig.4 illustrate the distributions for each element. Multi-element analysis reveals that the dominant substance in PM2.5 was S and that the dominant metal in PM2.5-10 was Si, as shown in Table 2.

Fig. 4. (a). Box and whisker plot for PM2.5. (b) Box and whisker plot for PM2.5-10.
Fig. 4.
 (a). Box and whisker plot for PM2.5. (b) Box and whisker plot for PM2.5-10.

Table 2. Mass concentrations of elements in PM2.5 and PM2.5-10.

Based on the XRF test results, there are 15 elements in PM2.5 and PM2.5-10. Of these 15 elements, Pb, Mn, and Cr already have quality standards for ambient air in several countries. The National Ambient Quality Standards (NAAQS) of the USEPA state that the Pb concentration within a measurement period of three months should not exceed 0.15 µg m–3. The average Pb concentrations measured in PM2.5 and PM2.5-10 for the sampling period were 0.1655 µg m–3 and 0.0862 µg m–3, respectively. Hence, the Pb measured in PM2.5 exceeded the established quality standards NAAQS, and the measurement results in one sample for PM2.5 and PM2.5-10 indicated Pb concentrations of 0.6471 µg m–3 and 0.3112 µg m–3, respectively, also exceeded the NAAQS quality standard. If the Pb concentrations are compared instead with the Indonesian standard during the 24-hours measurement time of 2 µg m–3, then one of them exceeded this standard. The WHO’s 24-hours standard for Mn concentrations is 0.018 µg m–3, while the Mn concentrations measured by PM2.5 and PM2.5-10 were 0.0094 µg m–3 and 0.023 µg m–3, respectively, meaning that the quantity of Mn in PM2.5-10 exceeded the established quality standards. The 24-hours standard for Cr concentrations of the Ontario Ministry of the Environment Ambient Air Quality Criteria/Texas Commission on Environment Quality (OAQC/TCEQ) is 1.5 µg m–3, while the Cr concentrations measured by PM2.5 and PM2.5-10 are 0.0016 µg m–3 and 0.0038 µg m–3, respectively, indicating that the quantity of Cr has not exceeded the established quality standards.

Based on the XRF plots used in the identification of the metal elements in PM2.5 (Fig. 4), the seven elements with the highest concentrations (in order from greatest to smallest concentrations) are S > Na > K > Zn > Si > Fe > Pb. S has the highest concentration at 2080.2 ng m–3, with a mean value of 1090 ng m–3. This element is one of the markers of vehicle activity conducted with low quality fuel. Emissions from the transportation sector that contain high quantities of S are transformed in the atmosphere into SO2. This finding is backed by de Bruin et al. (2006) and the International Council on Clean Transportation (ICCT). A report by Miller and Jin. (2018) states that Indonesia has a high sulfur average, especially because it uses EURO II vehicles that run on diesel fuel with high S content (> 500 ppm).

Based on the XRF plots for the identification of metal elements in PM2.5-10 (Fig. 4), the seven elements with the highest concentrations are Si > Fe > Ca > Na > Al > Cl > S. The Si element has the highest concentration of 1305 ng m–3, with a mean value of 2767.8 ng m–3. In general, the Si element is derived from soil dust particles but can also come from industrial sources that release Si, depending on area conditions. This research was carried out in the Waru Industrial Zone, which has approximately 235 companies with a variety of industrial activities, including the manufacture of paper, glass, dye, and cosmetics. Some of these industrial activities can contribute to the Si found in PM2.5-10 around the Waru Industrial Area, Sidoarjo. As a comparison, the concentration of Si in PM around the Industrial Zone of the city of Raipur, Central India is very high, ranging between 6.6–102 µg m–3 (Petel et al., 2016).

Principal Component Analysis (PCA)

After the elements in PM are identified, the XRF measurement data are then analyzed using PCA to identify sources that contribute to air pollution around the Waru Industrial Area, Sidoarjo. PCA has been used widely in the study of particulate matter composition in the air (Hopke, 1999). PCA is widely used to reduce data dimension in metal distribution studies, it extracts a small number of latent factors, referred to principal components to analyze relationships among observed variables (Yeung et al., 2003; Han et al., 2006; Meza-Figueroa et al., 2007). The PCA results can explain how much information is explained by four main components (factors) based on the values of the eigenvalues (root characteristics) for each major component.

Initial eigenvalue values greater than one identifies factors one to four, which then become the potential sources of pollutants. Markers are selected for the factors based on loading values > 0.5, and the loading values for the elements meeting this criterion are given in bold type in Table 3.

Table 3. Factor analysis for PCA.

The first factor consists of the elements Fe, Ca, Al, and Si, Ti, Mn, Cr, Mg, and Cu. Elements Fe, Ca, Al, and Si are markers of industrial activities. Paper, glass, cosmetics, and electronics industries which can contribute to the source of these pollutants, especially the Si concentrations, are located around the sampling locations. The elements Si, Fe, Ca, and Al are markers of industrial activities around the Industrial Zone of the city of Raipur, Central India (Patel et al.,2016). Si is usually used for various purposes in various sizes, depending on the application, in the tire, ceramics, rubber, glass, textiles, paper, cosmetics, electronics, toothpaste, and other industries (NIOSH, 2002). The Ti and Mn elements in the first factor are markers derived from soil dust (de Bruin et al., 2006). The soil dust can come from activities around the sampling location. A field about 60 meters from the sampling location is being used for the construction of housing. The first factor also includes Cr, Mg, and Cu, which indicate soil as well as the vehicle sources (de Bruin et al., 2006). Around the Waru Industrial Estate, vehicles are used to transport goods, and vehicles from outside also pass through the streets within the Industrial Estate because they lead to the airport and out of town. After work hours, the streets are congested with frequent traffic jams. These vehicle activities can contribute to this pollutant source. It can be concluded that the first factor is a combination of industrial activities, soil dust, and vehicle activities. The correlation matrix in Table 3 shows the strength of the relationships between elements, and the values in bold indicate strong relationships between variables (correlation > 0.5). It can be seen from the correlation matrix that the elements Al, Si, Ca, Ti, and Fe are strong, so it is possible that these elements come from one source, most likely soil dust.

The second factor consists of the elements S and Pb, where Pb can come from activities in the metal smelting industry. There is a metal smelting industry within 2–6 km of the sampling location, so the Pb in the PM mostly likely came from this industry. Pb in PM can also come from iron and steel industry activities (Dai et al., 2015). S in this second factor indicates that there is a contribution from diesel vehicles (Chueinta et al., 2000; Begum et al., 2014). Around this sampling location, many diesel vehicles are used as means of transporting industrial products/services that will be distributed to the Surabaya area or other cities outside East Java. Indonesia still has a high Scontent in diesel fuel (5,000 ppm) (Santoso et al., 2008).

The third factor consists of Na and Cl, which indicates that there is a contribution of sea salt to air pollution around the Waru Industrial Area. The contribution of sea salt comes from sea water evaporation (Mukhtar et al., 2013). Na evaporation usually coincides with Cl (Santoso et al., 2008).

The fourth factor consists of K, which indicates biomass combustion emissions (Santoso et al., 2008). K is present in biomass smoke and is the result of an incomplete, open combustion process, such as open air burning of garbage, forest fires, and cigarette smoking

Conditional Probability Function (CPF)

The following are the results of using the CPF method, with the help of Microsoft Excel, for estimating PM pollutant sources based on wind direction and velocity data. The results are in the form of a radar plot graph that uses 16 cardinal directions as a prediction of the direction of the pollutant source. The plot is based on the element that has the highest concentration and contributes to each factor.

The first factor consists of the elements Fe, Ca, Al, Si, Ti, Mn, Cu, Cr, Zn, and Mg from industrial activities and soil dust. Of these, Si, Al, and Fe have the highest concentrations. Si, Fe, and Al (Fig. 5) are estimated to originate from the north to northeast with a probability of 0.5. An industrial estate lies north to northeast from the sampling site. The industrial processes there can release Si, Fe, and Alparticles into the air, so it is entirely possible that the estate could be the source.

The second factor consists of Pb and S, which come from the smelting industry and diesel vehicles. S has a higher concentration than Pb. The results of overlaying CPF radar plot graphs for the elements Pb and S can be seen in Fig. 5, respectively. The results indicate that the Pb pollutant sources (non-ferrous metal industry activities) are from the west to southwest with a probability of 0.4. A steel smelting plantis 2–6 km from the sampling location in that direction, and in the southwest, there is also a metal village 2.3 km from the sampling location. Surabaya is flanked by industrial regions, such as Gresik and Sidoarjo, in that direction, which house Rungkut Industri Raya and Brebek Industry as well as steel smelting (Ahmad and Santoso, 2016).

In addition to contributions from the non-ferrous metal industry, pollutant emission source factors in the second factor can also come from diesel vehicles, because high concentrations of S are present. The result of overlaying the CPF radar plot graph for the S element can be seen in Fig. 5. An industrial estate lies in the north to northeast direction for the sampling location. Many vehicles pass through the industrial estate because its roads lead to the airport and out of town. These S concentrations are likely the result of high S concentrations in motor vehicle fuels used by the many vehicles passing through the area. Thus, the locally formed S condenses on to the surfaces of the existing particles (McMurry and Friedlander, 1979).

Fig. 5. CPF plot for the estimated location of the source of pollutants.
Fig. 5.
 CPF plot for the estimated location of the source of pollutants.

The third factor consists of the elements Na and Cl, which are contributions from sea salt. The CPF radar-plot-overlay results for Na and Cl can be seen in Fig. 5. The sea salt arrived from the northeast to the east with a probability of 0.4–0.45. Sea salt is made possible because the area borders the Madura Strait to the east, this conclusion is supported by Ahmad and Santoso (2016).


PM2.5 and PM10 concentrations around the Waru Industrial Area ranged from 2.65 to 32.68 µg m–3 and from 14.69 to 72.27 µg m–3, respectively. Mass concentrations for PM2.5 and PM10 at the time of measurement in the form of time series did not exceed the Indonesian quality standard but did exceed the WHO quality standard. XRF spectroscopy identified 15 metal elements in fine particles (PM2.5) and coarse particles (PM2.5-10), namely, Na, Mg, Al, Si, S, K, Ca, Ti, Cr, Mn, Fe, Cu, Zn, Pb, and Cl. The multi-element concentrations were analyzed using PCA to produce four components of PM pollutant sources around the Waru Industrial Area, Sidoarjo, East Java. Amixture of industrial activities, soil dust, and vehicle activity contributed to the first factor; metal smelting and diesel vehicle emissions contributed to the second factor; sea salt contributed to the third factor; biomass combustion emissions contributed to the fourth factor. Estimations of the locations of pollutant sources using the CPF method show that the PM pollutants come from industrial activities, vehicle activities, metal smelting activities, and sea salt emissions.


The author would like to thank the Department of Environmental Engineering, Faculty of Civil, Environmental and Geo Engineering, Institute of Technology Sepuluh Nopember (ITS) for facility and academic support of this research. Appreciation and gratitude are also extended to the Directorate of Research and Community Service, Directorate General of Research and Development, Ministry of Research, Technology and Higher Education of the Republic of Indonesia for providing project funds. This manuscript is part of Air Quality Monitoring Project under the contract: 992/PKS/ITS/2019 based on 6/E/KPT/2019. A thank you is also for the Center of Nuclear Technology for Materials and Radiometry, BATAN, Bandung, which provided help in supervising the research, and in particular sample preparation and data analysis.


  1. Ahmad, E.F. and Santoso, M. (2016). Analisis karaterisasi konsentrasi dan komposisi partikulat udara (Studi case: Surabaya) J. Kim. Valensi 2: 97–103. (in Indonesian) http://dx.doi.org/10.15408/jkv.v0i0.3602

  2. Ashbaugh, L., Malm, W. and Sadeh, W. (1985). A residence time probability analysis of sulfur concentrations at grand Canyon National Park. Atmos. Environ. 19: 1263–1270. https://doi.org/10.1016/0004-6981(85)90256-2

  3. Begum, B.A., Biswas, S.K., Markwitz, A. and Hopke, P.K. (2010). Identification of sources of fine and coarse particulate matter in Dhaka, Bangladesh. Aerosol Air Qual. Res. 10: 345–353. https://doi.org/10.4209/aaqr.2009.12.0082

  4. Begum, B.A., Kim, E., Biswas, S.K. and Hopke, P.K. (2014). Investigation of sources of atmospheric aerosol at urban and semi-urban areas in Bangladesh. Atmos. Environ. 38: 3025–3038. https://doi.org/10.1016/j.atmosenv.2004.02.042

  5. Chavent, M., Guégan, H., Kuentz, V., Patouille, B. and Saracco, J. (2009). PCA and PMF based methodology for air pollution sources identification and apportionment. Environmetrics 20: 928–942. https://doi.org/10.1002/env.963

  6. Chueinta, W., Hopke, P.K. and Paatero, P. (2000). Investigation of sources of atmospheric aerosol at urban and suburban residential areas in Thailand by positive matrix factorization. Atmos. Environ. 34: 3319–3329. https://doi.org/10.1016/S1352-2310(99)00433-1

  7. Dai, Q.L., Bi, X.H., Wu, J.H., Zhang, Y.F., Wang, J., Xu, H., Yao, L., Jiao, L. and Feng, Y.C. (2015). Characterization and source identification of heavy metals in ambient PM10 and PM2.5 in an integrated iron and steel industry zone compared with a background site. Aerosol Air Qual. Res. 15: 875–887. https://doi.org/10.4209/aaqr.2014.09.0226

  8. de Bruin, Y.B., Koistinen, K., Yli-Tuomi, T., Kephalopoulos, S. and Jantunen, M. (2006). A review of source apportionment techniques and marker substances, available for identification of personal exposure indoor and outdoor sources of chemicals. European Commission Directorate General Joint Research Center, Italy.

  9. Deshmukh, D.K., Deb, M.K., Tsai, Y.I. and Mkoma, S.L. (2010). Atmospheric ionic species in PM2.5 and PM1 aerosols in the ambient air of eastern central India. J. Atmos. Chem. 66: 81–100. https://doi.org/10.1007/s10874-011-9194-1

  10. Dockery, D.W., Pope Ш, C.A., Xu, X., Spengler, J.D., Ware, J.H., Fay, M.E., Ferris, B.G. and Speizer, F.E. (1993). An association between air pollution and mortality in six US cities. N. Engl. J. Med. 329: 1753–1759. https://doi.org/10.1056/NEJM199312093292401

  11. Han, Y., Du, P., Cao, J. and Posmentier, E.S. (2006). Multivariate analysis of heavy metal contamination in urban dusts of Xi’an, Central China. Sci. Total Environ. 355: 176–186. https://doi.org/10.1016/j.scitotenv.2005.02.026

  12. Hassanvand, M.S., Naddafi, K., Faridi, S., Nabizadeh, R., Sowlat, M.H., Momeniha, F., Gholampour, A., Arhami, M., Kashani, H., Zare, A., Niazi, S., Rastkari, N., Nazmara, S., Ghani, M. and Yunesian, M. (2015). Characterization of PAHs and metals in indoor/outdoor PM10/PM2.5/PM1 in a retirement home and a school dormitory. Sci. Total Environ. 527–528: 100–110. https://doi.org/10.1016/j.scitotenv.2015.05.001

  13. Henry, R.C. and Hidy, G.M. (1979). Multivariate analysis of particulate sulfate and other air quality variables by principal components. Part I. Annual data form Los Angeles and New York. Atmos. Environ. 13: 1581–1596.

  14. Holst, J., Mayer, H. and Holst, T. (2008). Effect of meteorological exchange conditions on PM10 concentration. Meteorol. Z. 17: 273–282. https://doi.org/10.1127/0941-2948/2008/0283

  15. Hopke, P.K., Xie, Y., Raunemaa, T., Biegalski, S., Landsberger, S., Maenhaut, W., Artaxo, P. and Cohen, D. (1997). Characterization of the gent stacked filter unit PM10 sampler. Aerosol Sci. Technol. 27: 726–735. https://doi.org/10.1080/02786829708965507
  16. Hopke, P.K. (1999). An introduction to source receptor modeling. In Elemental analysis of airborne particles, Landsberger, S. and Creatchman, M. (Eds.), Gordon and Breach Science Publishers, Amsterdam, pp. 273–315. 

  17. Hopke, P.K. and Jaffe, D.A. (2020). Letter to the editor: Ending the use of obsolete data analysis methods. Aerosol Air Qual. Res. 20: 688–689. https://doi.org/10.4209/aaqr.2020.01.0001

  18. Lau, A., Yuan, Z., Turner, J. and Yadav, V. (2012). Integrated data analysis and characterization of particulate matter in Hong Kong. Environmenttal Protection Departement HKSAR.

  19. McMurry, P. and Friedlander, S.K. (1979). New particle formation in the presence of an aerosol. Atmos. Environ. 13: 1635–1651. https://doi.org/10.1016/0004-6981(79)90322-6

  20. Meza-Figueroa, D., De la O-Villanueva, M. and De la Parra, M.L. (2007). Heavy metal distribution in dust from elementary schools in Hermosillo, Sonora, México. Atmos. Environ. 41: 276–288. https://doi.org/10.1016/j.atmosenv.2006.08.034

  21. Miller, J. and Jin, L. (2018). Global Progress Toward Soot-free Diesel Vehicles in 2018. ICCT (Internasional Council on Clean Transportation) report. https://theicct.org/publications/global-progress-toward-soot-free-diesel-vehicles-2018

  22. Mukhtar, R., Esron, H., Hari, W., Muhayatun, S. and Syukria, K. (2013). Chemical components PM2.5 and PM10 in ambient air in serpong-tangerang. Ecolab. 7: 1–48.

  23. Mukono, H. (2011). Health aspects of air pollution. Pusat Penerbitan dan Percetakan Unair (AUP), Surabaya

  24. National Institute for Occupational Safety and Health (NIOSH) (2002). Hazard review. Health Effects of Occupational Exposure to Respirable Crystalline Silica. DHHS Publication No. 2002-129. National Institute for Occupational Safety and Health, Cincinnati.

  25. Patel, K.S., Gupta, S., Ramteke, S., Rajhans, K.P., Nava, S. and Lucarelli, F. (2016). Silica particulate pollution in central India. J. Environ. Prot. 7: 170–175. https://doi.org/10.4236/jep.2016.72015

  26. Pekney, N.J, Davidson, C.I., Zhou, L. and Hopke, P.K. (2006). Aplication of PSCF and CPF to PMF- modeled sources of PM2.5 in Pittsburgh. Aerosol Sci. Technol. 40: 952–961. https://doi.org/10.1080/02786820500543324

  27. Pitaloka, A.P. and Adriyani, R. (2016). Paparan PM10 dan keluhan kesehatan mata pekerja bagian produksi PT. Varia Usaha Beton, Sidoarjo. J. Ilm. keperawatan 2: 62–71. (in Indonesian) https://journal.stikespemkabjombang.ac.id/index.php/jikep/article/view/27

  28. Proias, T., Nastos, P.T., Larissi, I.K. and Paliatsos, A.G. (2010). PM10 concentrations related to meteorology in Volos, Greece. AIP Conf. Proc. 1203: 1091–1096. https://doi.org/10.1063/1.3322316

  29. Reff, A., Elberly, S.I., and Bhave, P.V. (2007). Receptor modeling of ambient particulate matter data using positive matrix factorization: Review of existing methods. J. Air Waste Manage. Assoc. 57: 146–154. https://doi.org/10.1080/10473289.2007.10465319

  30. Rixson, L. Riani, E. and Santoso, M. (2015). Characterization of long term exposure of particulate matter at Puspiptek Serpong-South Tangerang. J. Ilmiah AI & R 11: 51–64. (in Indonesian)

  31. Samet, J.M., Zeger, S., Dominici, F., Curriero, F., Coursac, I., Dockery, D.W., Schwartz, J. and Zanobetti, A. (2009). The National morbidity, mortality, and air pollution study. Part II: Morbidity and mortality from air pollution in the United States. Res. Rep. Health Eff. Inst. 94: 5–70.

  32. Santoso, M., Hopke, P.K., Hidayat, A. and Diah, D.L. (2008). Sources identification of the atmospheric aerosol at urban and suburban sites in Indonesia by positive matrix factorization. Sci. Total Environ. 397: 229–237. https://doi.org/10.1016/j.scitotenv.2008.01.057

  33. Santoso, M., Lestiani, D.D. and Hopke, P.K. (2013). Atmosfir black carbon in PM2.5 in Indonesian cities. J. Air Waste Manage. Assoc. 63: 1022–1025. https://doi.org/10.1080/10962247.2013.804465

  34. Santoso, M. and Lestiani, D.D. (2014). Application of ED XRF in supporting national program of air quality improvement in Indonesia. XRF Newsletter 26: 8–13.

  35. Santoso, M., Marselina, M., Lestiani, D.D. and Mukhtar, R. (2016). Characteristics of ambient and ambient air particulates around non-formal industrial estates. J. Sains Teknol. Nuklir Indones. 17: 49–58. (in Indonesian)

  36. Santoso, M.S., Taftazani, A., Darsono, M.R., Lestiani, D.D., Damastuti, E. and Kurniawati, S. (2010). The role of nuclear analysis techniques in health and the environment. Jurnal Seminar Nasional Keselamatan Kesehatan dan Lingkungan VI, Jakarta.

  37. Schaumann, F., Borm, P.J.A., Herbrich, A., Knoch, J., Pitz, M., Schins, R.P.F., Luettig, B., Hohlfeld, J.M., Heinrich, J. and Krug, N. (2004). Metal-rich ambient particles (Particulate Matter2.5) cause airway inflammation in healthy subjects. Am. J. Respir. Crit. Care. Med. 170: 898–903. https://doi.org/10.1164/rccm.200403-423OC

  38. Seinfeld, J.H., and Pandis, S.N. (2006). Atmospheric chemistry and physics from air pollution to climate change. John Wiley and Sons Inc., Canada.

  39. Thurston, G.D. and Spengler, J.D. (1989). A quantitative assessment of source contributions to inhalable particulate matter pollution in metropolitan Boston. Atmos. Environ. 19: 9–25. https://doi.org/10.1016/0004-6981(85)90132-5

  40. World Health Organization (WHO) (2017). Health risks of heavy metals from long-range transboundary air pollution. WHO Regional Office for Europe, Copenhagen. https://apps.who.int/iris/handle/10665/107872

  41. Yeung, Z.L.L., Kwok, R.C.W. and Yu, K.N. (2003). Determination of multi-element profiles of street dust using energy dispersive X-ray fluorescence (EDXRF). Appl. Radiat. Isot. 58: 339–346. https://doi.org/10.1016/S0969-8043(02)00351-2

  42. Zannaria, N.D., Rosmini, D. and Santoso, M. (2009). Chemical characteristics of exposed particulate exposure. J. Sains Teknol. Nuklir Indones. IX: 37–50. (in Indonesian)

Aerosol Air Qual. Res. 20 :1910 -1921 . https://doi.org/10.4209/aaqr.2019.11.0590  

Share this article with your colleagues 


Subscribe to our Newsletter 

Aerosol and Air Quality Research has published over 2,000 peer-reviewed articles. Enter your email address to receive latest updates and research articles to your inbox every second week.

79st percentile
Powered by
   SCImago Journal & Country Rank

2023 Impact Factor: 2.5
5-Year Impact Factor: 2.8

Aerosol and Air Quality Research partners with Publons

CLOCKSS system has permission to ingest, preserve, and serve this Archival Unit
CLOCKSS system has permission to ingest, preserve, and serve this Archival Unit

Aerosol and Air Quality Research (AAQR) is an independently-run non-profit journal that promotes submissions of high-quality research and strives to be one of the leading aerosol and air quality open-access journals in the world. We use cookies on this website to personalize content to improve your user experience and analyze our traffic. By using this site you agree to its use of cookies.