key: cord-0857651-lo4hqokv
authors: Abbasi, Sajjad; Keshavarzi, Behnam
title: Source identification of total petroleum hydrocarbons and polycyclic aromatic hydrocarbons in PM10 and street dust of a hot spot for petrochemical production: Asaluyeh County, Iran
date: 2019-02-28
journal: Sustainable Cities and Society
DOI: 10.1016/j.scs.2018.11.015
sha: 942557ee1d2f2f47dea421b787f9aee2da8ba9e7
doc_id: 857651
cord_uid: lo4hqokv

Abstract Total petroleum hydrocarbons (TPH) and polycyclic aromatic hydrocarbons (PAHs) are important pollutants that affect public health in urban areas, especially in developing and oil-rich countries such as Iran. This assesses the relationship between TPH and PAHs in street dust and suspended dust, and investigates toxicity level in the urban environment of the most important petrochemical center in Iran. For this purpose, 21 and 48 street dust samples were collected for TPH analysis and PAH analysis, respectively, in Asaluyeh County. Moreover, seven air dust samples were taken for PAH analysis. TPH concentrations ranged between 240 and 4400 mg kg−1, with a mean of 1371.43 mg kg−1. The maximum ∑PAH concentration (6016.3 mg kg−1) was detected in a petrochemical complex while the minimum ∑PAHs content (16.93 mg kg−1) was measured in an urban area. The mean concentrations of total PAHs in street dust particles were 491.35 mg kg−1 in summer and 304.04 mg kg−1 in winter. The results indicated that PAH concentration in summer was higher. PAH sources were identified using both PAHs ratios and robust statistical methods such as Generalized Estimating Equations (GEE), backward GEE, logistic regression, principal components analysis (PCA) in conjunction with multiple linear regression (MLR) and positive matrix factorization (PMF). The results showed that PAH species generally originate from pyrogenic sources and about 0.08% of TPH was typically PAHs. However, petrogenic sources of PAHs in the industrial areas were 11.2 times more abundant than in urban areas. Also, backward GEE model demonstrated that TPH is more influenced by HMW PAHs, particularly indene. Estimated incremental lifetime cancer risk (ILCR) was higher than 10-4, showing that Asaluyeh inhabitants (especially children and indoor workers) are probably exposed to cancer risk, particularly through dermal contact and dust ingestion.

Urban design and geographic location are two important factors in amplifying effects of environmental contaminants. Over the last few decades, several studies were carried out to link urban quality with urban design. For instance, Nikolopulou et al. (2004) investigated microclimate in multiple open urban spaces. Based on air ventilation assessment, the development of urban planning guidelines was considered by Ng (2009) . Oke (1988) also investigated optimum street canyon characteristics according to diverse climate-related criteria in urban regions and demonstrated a number of useful relationships between the geometry and the microclimate of urban street canyons. Previous studies are potentially helpful in establishing guidelines governing street dimensions to be used by urban designers. On the other hand, rapid urbanization, population growth and consequent accumulation of contaminants pose major challenges for people and planners. For more explanation, in 2003, Hong Kong was hit by severe acute respiratory syndrome, which led to numerous fatalities. Thus, a number of focused studies were conducted with a recommendation being utilization of air ventilation systems (Ng, 2009) .

Asaluyeh County is the most important industrial center in Iran. The county is limited to the Persian Gulf in the south and a mountain range in the north (Abbasi, Keshavarzi, Moore, & Mahmoudi, 2018; Abbasi, Keshavarzi, Moore, Turner et al., 2019) . The northern barrier prevents the southwest-northeast wind current to disperse the pollution load and hence resulting in increasing the air pollution potential. Thus, it is important to select appropriate methods to investigate contaminants and their effects on human health. Accelerated urban and industrial development and a huge petrochemical complex are already reported to be responsible for the high level of pollutants in Asaluyeh port (Alahverdi & Savabieasfahani, 2012) .

Generally, atmospheric pollution is a primary concern in urban areas with particulate matter (PM) leading to the most deaths globally among environmental threats (WHO, 2006) . PM 10 is capable of reaching conductive tracheobronchial airways through inhalation (Harrison & Yin, 2000) . On the other hand, street dust serves as both source and sink for atmospheric PM. Street dust is generated from soil, cement, sea salt (Zhao, Feng, Zhu, & Jianhui, 2006) , evaporation of vehicle fuel (Duffy, Nelson, Ye, & Weeks, 1999) , direct emission from vehicle exhaust (Maricq, 2007) , pieces of rubbers Abbasi, Keshavarzi, Moore, Turner et al., 2019) , abrasion of the street surface and the resuspension of accumulated dust on streets (Thorpe & Harrison, 2008) . Street dust thus serves as an indicator of urban life quality and reflects contamination of different environments like soil, water and air (Wong, Li, & Thornton, 2006) .

TPH is a term used to describe a wide family of chemical compounds that originally derive from crude oil (Todd, Chessin, & Colman, 1999) . When TPH is adsorbed by street dust, rainfall can wash TPH and transfer it to surface water. However, there are risks for shallow aquifers, which may be used for private wells for drinking water purposes (Amadi, Abbey, & Nma, 1996; Albers, 1995; Todd et al., 1999) . Also, different TPH compounds may be separated from the original mixture (depending on the chemical properties), evaporated, dissolved into the groundwater or attached to soil particles. Hence, human can be exposed to these compounds through respiration, ingestion and dermal pathways (Todd et al., 1999) . TPH provides an estimate of the concentration of higher molecular weight hydrocarbons, while PAHs provide an estimate of more toxic oil components (Wade et al., 2011) .

PAHs are a great concern as some carcinogenic species (Chrysikou, Gemenetzis, & Samara, 2009 ) are always present in the soil, air, and sediments (Giacalone, Gianguzza, Mannino, Orecchio, & Piazzese, 2004) . PAHs are introduced into the environment by both anthropogenic sources and/or natural processes (Junker et al., 2000) . Anthropogenic sources are divided into petrogenic and pyrogenic ones. Pyrogenic PAHs often originate from incomplete combustion of fossil fuels, biomass and municipal wastes (Yunker et al., 2002) . Petrogenic PAHs are spills of petroleum-derived products and the diagenesis of organic matter in anoxic environment (Lima, Farrington, & Reddy, 2005) .

In general, 2-and 3-ring PAH compounds are more volatile ) (Cui et al., 2016; Guo et al., 2009; Liu, Liu, Chen, Wang, & Cao, 2013) and mostly occur in the gas phase of the atmosphere and can deposit on the surface of particles via wet deposition, while less volatile i.e., 5-6 ring PAHs ) (Cui et al., 2016; Guo et al., 2009; Liu et al., 2013) mostly tend to deposit on the surfaces of particles via wet and dry deposition (Mannino & Orecchio, 2008) . Among the hazardous contaminants, PAHs are serious because of their mutagenicity and carcinogenicity (National Research Council (NRC), 1983) . Street dust contaminated with PAHs indicates a higher health risk to children (Keshavarzi, Sajjad Abbasi, Moore, Delshab, & Soltani, 2017; Krugly et al., 2014) .

As already mentioned, PAH species have different characteristics and several factors, including temperature and humidity, affect their concentration. In order to identify sources, robust statistical methods such as Generalized Estimating Equations (GEE), backward GEE, logistic regression and Positive Matrix Factorization (PMF) models must be employed. Each method offers particular advantages that depend on the nature of data and the desired research focus (Hu, Goldberg, Hedeker, Flay, & Pentz, 1998) . The reasons for selecting these statistical methods in the current study are given in Section 2.3.1.

The specific objectives of this study are to: (1) measure the concentrations of PAHs and TPH in street dust of urban and industrial areas in the winter and summer; (2) measure the concentrations of PAHs in the PM 10 fraction of street dust samples in urban and industrial areas;

(3) identify various sources of PAHs by GEE, backward GEE, logistic regression, PMF, diagnostic ratios and principal components analysis (PCA) in conjunction with multiple linear regression (MLR); (4) determine how TPH concentration is affected by different members of PAHs; and (5) evaluate human exposure risk to PAHs in the Asaluyeh street dust considering different scenarios.

Iran is located in the mid-latitudinal belt of arid and semiarid regions. Asaluyeh County is located in Bushehr province, south Iran ( Fig. 1) . Average annual temperature and rainfall in Asaluyeh are 25.9°C and 35 mm, respectively. Asaluyeh is the biggest gas field in the world shared by Iran and Qatar. It contains 51,000 km 3 gas in place and 56 billion barrels (8.9 × 10 9 m 3 ) condensate in place on both sides. The Iranian side accounts for 10% of the world's and 60% of Iran's total gas reserves. Many petrochemical plants, gas powerhouses and semi-heavy industries are active in Asaluyeh. Related industries include electrical and electronic industries, instrumentation industries and chemical industries. A population of about 73,958 inhabitants live in Asaluyeh.

Before sampling, amber glass containers were first washed with soapy water, and then with distilled water, n-hexane, and finally dried at 180°C for 3 h in an oven. Furthermore, for more certainty, screw aluminium foil caps were used to provide a gas-tight seal. Then, 48 and 21 street dust samples, each representing an area of 5-10 m 2 , were collected from the surface of different streets for PAHs and TPH analyses, respectively ( Fig. 1, Table S1 ). Regarding PAH samples, 43 and 5 street dust samples were collected in summer and winter, respectively. Approximately 100 g of dust samples were collected in each site and passed through a 2-mm sieve and transferred to an ice chest, immediately. The gathered samples were rapidly transported to the laboratory of the Isfahan University of Technology.

Two sampling stations were considered for (re-)suspended dust: (i) the first four days in an industrial area within the PSEEZ (R1), and (ii) the second three days in an urban area (R2) ( Fig. 1 and Table S1 ). In each station two samplers were used, one for PAHs analysis and the other for PM concentration. Samples of particulate matter (PM) were collected using an ECHO PM ambient filter sampler (TECORA, Italy) during almost ten successive days (seven samples) at both sites on PTFE filter papers (2 μm pore size and 46.2 mm diameter; Tisch Scientific, USA). To avoid local disturbance from, for example, automobiles, the sampler was set up at a height of 3-4 m, and in accordance with United States Environmental Protection Agency (USEPA) reference methods for PM 10 , the sampling flow rate was set at 16.67 L min −1 (EPA, 2006) . After 24 h, the filters were transferred to plastic petri dishes already washed with deionized water and covered with Al sheet and finally stored at 4 ᵒC until being transported to the laboratory for PAHs extraction.

For the determination of PM 10 PAHs, filter extracts were analyzed according to the USEPA method TO 13 A (USEPA, 1999) . PM10 filters were extracted in Soxhlet with 125 mL of dichloromethane (CH2Cl2) for 18 h (USEPA, 1999). Then, for eliminating possible interferences, the extracts were separated/pre-concentrated through the cleanup procedure using a silica gel column (activated with 5% of distilled H2O) (Teixeira, Mattiuzi, Agudelo-Castañeda, de Oliveira Garcia, & Wiegand, 2013) . The cleanup was accomplished for each sample with four solvent fractions with different polarities: first fraction: 20 mL of hexane (aliphatic); second fraction: 10 mL of hexane, 10 mL of dichloromethane; third fraction: 15 mL of hexane, 5 mL of dichloromethane; and fourth fraction: 20 mL of dichloromethane. Rotary evaporation followed by a mild stream of pure nitrogen gas was conducted for concentrating volumes for fractions 2 and 3. Finally, 1 mL of dichloromethane was added for subsequent analysis.

For PAH analysis of street dust, approximately 5 g of each dried sample was spiked with surrogate standards (Pyrene-D10, lot: 10,510 semi-volatile internal standards), and an ultrasonic bath (KUDOS, SK3210LHC model) was used during the extraction. Pyr-D10 was used as internal standard for the calibration of all species of PAHs. The samples were then extracted using a 30-mL mixture of organic solvents (n-hexane and dichloromethane [DCM] in a 1:1 v/v) for 30 min at room temperature. Amorphous sodium sulfate is added to the solution and then 2 mL of dried extracts were taken for the next step. A silica gel column was used to clean up the extract similar to PM10 samples. The samples were analyzed, for PAHs, following the EPA 3550B (extraction), EPA 3630C (clean up), and EPA 8310 (determination) procedures, and EPA 418 for TPH.

The samples were analysed for PAHs using a Hewlett-Packard (HP) 1090 high-performance liquid chromatography (HPLC) system equipped with a HP-1046 fluorescence detector. Considering the recommendation by the International Agency for Research on Cancer . The extracts were concentrated to 1 ml by a rotary vacuum evaporator. 20 μl of each extract was used for PAHs analysis. The mobile phase was acetonitrile/water in gradient mode at a flow rate of 1 ml/min and the temperature was set at 35°C. Generally, average recoveries for PAHs and TPH in dust were approximately 88-95% for all 16 measured PAHs and TPH. The detection limit (DL) of PAHs for dust ranged from 0.05 to 0.4 μg/kg and precision was 4-8%. Reagent blanks, analytical duplicates/replicates, and analysis of the standard reference material (Dr. Ehrenstorfer GmbH Alkanes-Mix 10, and Sigma-Aldrich Co. LLC EPA 525 PAH Mix A and EPA 525 PAH Mix B) were processed. 

EPA's Positive Matrix Factorization (PMF) Model is a mathematical receptor model that provides scientific support for the development and review of environmental contaminations such as PAHs. The current study employed EPA PMF v5.0.14.

Generally, the descriptive statistics, PCA-MLR, GEE, backward GEE and logistic regression were done using SPSS version 19.0 software. PCA-MLR is a statistical procedure that uses an orthogonal transformation to convert a set of observations of possibly correlated variables (entities of which take on various numerical values) into a set of values of linearly uncorrelated variables called principal components. Also, MLR was used to quantify source contributions to samples based on pollutants. The primary function of PCA was a reduction in the number of variables while retaining as much as possible of the original information; therefore, variables with similar characteristics were grouped into factors (Fang et al., 2004) . When PCA was performed with Varimax normalized rotation, each principal component score contained information on all PAHs combined into a single number, while the loadings indicated the relative contribution for PAHs species made to that score (Kong et al., 2012; Yongming, Peixuan, Junji, & Posmentier, 2006) .

Regarding the GEE and backward GEE, the variables in this research were the response (dependent) variable. Also, the variables (PAH species) , and their interaction effects were continuous predictor (independent) variables. For example, if we have five variable, such as A, B, C, D and E, the general equation of GEE is as in Eq. (1).

where i s are model parameters (coefficients) and t is the random component of the model, which follows a random distribution with a mean value of 0. As can be seen in Eq.

(1), this model contains the main effects and all 2-way interaction effects of the predictors A, B, C, D and E. However, wherever the effects were not significant (Sig > 0.05), backward GEE and GEE were used to remove ineffective parameters. The final model had maximum accuracy.

Logistic regression analysis was carried out to compare petrogenic and pyrogenic sources of PAHs in industrial and urban areas. All statistical analyses were performed using the SPSS statistical package software (version 19). A probability of Sig < 0.05 was considered statistically significant. Generally, based on Sig, one can find if the relation between parameters is significant or not. When the level of Sig is lower than 0.05 it indicated significant correlation and as such (Sig < 0.5), beta (B) is used to determine the relationship intensity.

Each statistical method has its own advantages and disadvantages and in the current study several methods were considered to get more accurate results. For this purpose PCA-MLR, PMF and GEE were performed. Generally, PCA's key advantages are (i) its low noise sensitivity, (ii) the decreased requirements for capacity and memory, and (iii) increased efficiency given the processes taking place in a smaller dimensions (Karamizadeh et al., 2013) . Other advantages of PCA include: (i) Lack of redundancy of data given the orthogonal components; (ii) Reduced complexity in images' grouping with the use of PCA; (iii) Smaller database representation since only the trainee images are stored in the form of their projections on a reduced basis; (iv) Reduction of noise since the maximum variation basis is chosen and so the small variations in the background are ignored automatically (Phillips et al., 2005; Srinivasulu Asadi et al., 2010) ; (v) Results cannot be weighted to account for uncertainties in the measured data; (vi) PCA models cannot properly handle missing data or values below the detection limit (both of which commonly occur in environmental measurements); and (vii) MLR analysis of the factor scores was used to quantify source contributions to samples based on pollutants (Larsen & Baker, 2003; Ashayeri et al., 2018) . PCA's key disadvantages of PCA are: (i) The covariance matrix is difficult to be evaluated in an accurate manner (Phillips et al., 2005) ; (ii) Even the simplest invariance could not be captured by the PCA unless the training data explicitly provides this information (Li et al., 2008) ; and (iii) PCA is not as powerful as other statistical methods in small datasets (Karamizadeh et al., 2013) .

It is worth mentioning that weaknesses and limitations of PMF method according to EPA reports are important and include: 1) PMF models require large datasets on measured concentrations (preferable > 100 samples), 2) Analysis is limited by the accuracy, precision, and range of species measured at the receptor (e.g. ambient monitoring) sites, 3) A determination must be made of how many 'factors' to retain, 4) Emission sources have to be deduced by interpreting these factors, 5) Information is needed on source profiles or existing profiles in order to verify the representativeness of the calculated source profiles and uncertainties in the estimated source contributions, 6) The method relies on many parameters and initial conditions and model input; results are sensitive to the pre-set parameters.

A specific advantage of GEE is its ability to robustly estimate variances of the regression coefficient for data exhibiting high correlation between repeated measurements (Abbasi, Keshavarzi, Moore, Mahmoudi, 2018; Ballinger, 2004; Ghisletta & Spini, 2004; Hu et al., 1998) . Other advantages include: 1) The GEE algorithm has been incorporated into many major statistical software packages, including SAS, STATA, R, and S-PLUS, 2) the model is robust to the misspecification of correlation structure, because the parameter estimate remain consistent, 3) when the working correlation structure is correctly specified, the parameter estimates from GEE are efficient (Liang & Zeger, 1986; Khajeh-Kazemi et al., 2011) . GEE's limitations include: 1) in the case misspecification of the working correlation matrix the parameter estimates from GEE can be inefficient (Liang & Zeger, 1986) , 2) the consistency of the estimated parameter depends only on the correct specification of the mean model; but there is no universally accepted test for goodness-of-fit for mean model in GEE that extends beyond binary dependent variable (Barnhart & Williamson, 1998) , 3) GEE parameter estimates are sensitive to outliers or contaminated data and in that cases GEE fails to give consistent estimators and more seriously will lead to incorrect conclusion (Preisser & Qaqish, 1996 Qu & Song, 2004) , 4) In GEE due to lack of an objective function it is complicated to use model selection criterion (Pan, 2001) , 5) if the correlation parameters are not consistently estimated, GEE fails to produce consistent estimators (Crowder, 1986 (Crowder, , 1995 .

The incremental lifetime cancer risk (ILCR) resulting from exposures to PAHs in Asaluyeh's street dust was calculated using the USEPA standard model (residential/ Commercial and industrial scenarios) (USEPA, 1991) . This model is used in the current study to assess exposure risk of adults and children to PAHs in street dust. The following assumptions underlie the model applied in this study: (a) Human beings are exposed to urban surface dust through three main pathways: ingestion, inhalation, and dermal contact with dust particles; (b) Intake rates and particle emission can be approximated by those developed for soil particles; (c) Some exposure parameters of people in the observed areas are similar to those of reference populations; (d) The total carcinogenic risk could be computed by summing the individual risks calculated for the three exposure ways; (e) The cancer risk is assessed based on exposure under a specific type of land use pattern over the entire lifetime. The doses received through the ingestion (ing), inhalation (inh) and dermal contact (der) of PAHs were calculated following the guidance (USEPA, 1989).

where CS is the sum of the converted PAH levels (μg kg −1 ) based on the toxic equivalents of BaP using the toxic equivalency factor (TEF) listed in Table 1 , AT is the average lifespan (day), BW is body weight (kg), CSF is carcinogenic slope factor (mg kg −1 day −1 ) −1 , ED is the exposure duration (year), EF is the exposure frequency (day year -1 ), PEF is the particle emission factor (m 3 kg −1 ) and ABS is the dermal adsorption fraction. IR Inhalation is the inhalation rate (m 3 day −1 ), IR Ingestion is the dust intake rate (mg day −1 ), AF is the dermal adherence factor (mg cm −2 ) and SA is the dermal surface exposure (cm 2 day −1 ). CSF Ingestion , CSF Dermal , and CSF Inhalation of BaP were addressed as 7.3, 25, and 3.85 (mg kg -1 day −1 ), respectively, as determined by the cancercausing potential of BaP (USEPA, 1994) . Other parameters referred in the model for children (1-6 years old) and adults (7-31 years old) (in the residential/commercial and industrial areas) are based on the risk assessment guidance of USEPA and related publications (Table S2) .

As shown in Table 1 , large variability in TPH concentrations of sampling sites was observed. Urban areas typically exhibited higher TPH concentrations than industrial sites. The maximum, minimum and mean TPH are 4400, 240 and 1371.43 mg/kg, respectively. The high concentration of TPH is probably related to petroleum products, oil carrier trucks and automobiles. The TPH concentrations in refineries and petrochemical units are lower than those of the urban areas and far from the authors' expectation (Fig. S1 ). This can be the result of many factors such as washing the connecting roadways. Such roadways in refineries and petrochemical units (industrial areas) are regularly washed, and hence pollutants adsorbed to street dust are removed (Kennedy, 2003; Zhang, Dou et al., 2008) . Consequently, the concentration of pollutants in such places is lower than in regions where street dust remains on the street surface for a long time. According to Adeniyi and Afolabi (2002) , street dust pollution can be identified when TPH concentration is higher than 10-100 mg/kg, indicating that all collected samples in Asaluyeh are polluted. Hence, measuring TPH subunits such as PAHs seemed reasonable (Fig. S1 ).

Descriptive statistics of PAH concentrations in summer are presented in Table 1 . The highest ∑PAHs (6016.3 μg/kg) was found in P3 (a (Nisbet & LaGoy, 1992) .

sampling site in a petrochemical company), while minimum ∑PAHs (16.93 μg/kg) occurred in P42 (a sampling site in an urban area). During summer and winter, the mean concentration of total PAHs in street dust particles was 491.35 and 304.04 μg/kg, respectively. In general, ∑PAHs concentrations at petrochemical stations are higher compared to urban areas. Contrary to many previous studies that show a higher concentration of PAHs in winter, in the current study, total PAHs are higher in the summer (Fig. 2a) . As already mentioned during winter, meteorological conditions like shorter daylight periods, low temperature, low volatilization, lower atmospheric mixing height, reduced vertical dispersion due to thermal inversion, increased use of heaters, enhanced sorption to particles at lower temperature and low photochemical reactivity of PAHs would increase PAHs content during winter (Baalbaki et al., 2018; Gope, Masto, George, & Balachandran, 2018; Liu, Cheng, Ou et al., 2007; Periera-Netto, Krauss, Cunha, & Rego, 2006) . Nonetheless, this is not the case everywhere and depends on many other factors. For example, in the current study, heating devices are not frequently used and the weather is not that cold in winter. It was also found that PAHs concentration during monsoon season is lower than summer because the heavy washout and runoff can wash street dust and their adsorbed PAHs (Gope et al., 2018) . Therefore, rainfall and runoff can be important for washing street dust in the winter, rather than PAHs themselves. According to Nowak, Crane, Stevens, Hoehn, and Walton (2008) , precipitation helps wash away some surface contaminants, but PAHs are not readily washable. For example, Škrbić, Đurišić-Mladenović, Živančev, and Tadić (2019) reported similar results to those of the current study. The higher concentrations of PAHs in street dust sampled in summer might be attributed to the additional pollution sources in summer such as the emission from numerous motorcycles that are not driven during the cold months (Škrbić et al., 2019) . It is worth mentioning that PAHs emission from motorcycles was about 10 to 16 times greater than those from vehicles equipped with catalyst and low emission automobiles (Pham, Kameda, Toriba, & Hayakawa, 2013) . In addition, natural gas and diesel fuel are used in Asaluyeh for generating electricity. Considering hot and humid weather in the summer, a sharp . Both petrochemical and other industrial areas indicated high average concentrations of low molecular weight (LMW) PAHs (2 and 3 rings), while only urban areas demonstrated remarkable mean content of high molecular weight (HMW) PAHs (4 rings and above). The main pollution source of PAHs is believed to be motor vehicles in urban areas of the current study. Generally, the LMW PAH species are produced by crude oil products at low to moderate temperatures, such as wood and coal combustion (Wang et al., 2011) , but the HMW compounds are generated, for instance, when fuel is burned in engines at high temperatures (Mastral & Callen, 2000) . Molecular weight of 4 rings PAHs lies between those of LMW and HMW PAHs and their characteristics are comparable to the properties of HMW and LMW PAHs (Table 1) . Therefore, they may occur both in industrial regions (with high potential for generating LMW PAHs) and in urban areas (with high potential for producing HMWPAHs).

On average, the percentage of PAHs with 2, 3, 4, 5 and 6 rings are 1. 98, 30.90, 29.81, 33 .27 and 4.04%, respectively (Table 1 and Fig. 3) . The abundance of PAHs based on ring numbers in the street dust of Asaluyeh follows the order of 5 rings > 3 rings > 4 rings > 6 rings > 2 rings. The 5 ring compounds indicated that industrial and traffic pollutants are more abundant (Wilcke, Amelung, Martius, Garcia, & Zech, 2000). Obviously, PAH compounds with 2 and 3 rings have less molecular weight and hence are more volatile than PAHs with 4 to 6 rings Liu, Cheng, Ou et al., 2007) . Therefore, the greater concentration of PAH compounds with 2 and 3 rings at the mentioned stations indicates their recent release from pollution sources in the environment. High concentration of the 5 and 6 ring species probably reflects incomplete fuel combustion in industry and petroleum carrier trucks. These compounds have a strong tendency to deposit and bond to soil or dust particles (Khalili, Scheff, & Holsen, 1995) . Also, the mean of combustion PAHs (the sum of Flu, Pyr, BaA, Chr, B(b + k)F, and BaP (Rogge, Hildemann, Mazurek, Cass, & Simoneit, 1993; Rajput & Lakhani, 2009 )) to total PAHs in the street dust samples is 0.64 μg kg −1 and constitutes a significant proportion of the total PAHs concentration. According to International Agency for Research on Cancer (IARC), carcinogenic PAHs include Nap, BaA, Chr, BeP, BbF, BkF, BaP, DbahA, and Ind. The mean ratio of carcinogenic PAHs to total PAHs (CANPAHs/PAHs) is lower than the mean ratio of the non-carcinogenic PAHs to total PAHs (NCANPAHs/PAHs) indicating higher abundance of the non-carcinogenetic PAHs (Table 1) . Furthermore, considering Fig. 2c , the concentration of carcinogenic PAHs is higher than the value of non-carcinogenic PAHs, especially in regions with high traffic load. The carcinogenic PAHs concentration is also higher in summer (Fig. 2c) . Fig. 2b depicts lower LMW PAHs concentration in winter because as already pointed out LMW PAHs are easily washed by surface runoff in this season. In Previous studies the bioavailable fractions of organic pollutants is reported as the sum of water-soluble and acidsoluble fractions which that negatively correlate with the log K ow of organic pollutants. Therefore, HMW PAHs are more stable than LMW PAHs and rainfall washes LMW more efficiently (Wang, Zhu, & Zhang, 2015) (Fig. 2b) . Table S3 reveals that the concentrations of PAHs in Asaluyeh street dust are more varied and higher than Isfahan (Iran), Beijing (China), Cairo (Egypt) and Ulsan (South Korea), and lower than Tehran and Bushehr (Iran), Manchester (UK), Guangzhou (China), New Delhi (India) and Istanbul (Turkey). Generally speaking, comparison of contamination in different regions and sampling sites is almost impossible. The reason is that many factors (including climatically and geological) can greatly affect the results. Hence, one must always bear in mind that origin of pollution may also be geogenic rather than anthropogenic. Furthermore, the intensity of precipitation varies in different locations and thus may wash street dust and pollutants differently. Finally, different sampling and analytical methods must be taken into account when a comparison is made. It is worth noting that similar studies were conducted in Asaluyeh, Bushehr and Bandar Abbas cities (south of Iran) with almost the same climate, weather, and sampling methods. The main reason for differing concentrations of PAHs in Bandar Abbas are anthropogenic activity such as traffic load and industrial activity. More specifically, the industrial activity (petroleum units) and traffic load in Asaluyeh are higher than Bushehr and Bandar Abbas resulting in higher PAHs concentration in Asaluyeh. Thus, traffic load and industrial activity are two important factors that affect contamination of street dust. As mentioned, Asaluyeh is limited to the Persian Gulf in the south and mountains to the north. Thus, the sweeping effect of wind blowing from sea toward mountains (the northern barrier) highlights air and street dust pollution due to lack of dispersion. This sweeping effect of wind can be severely reduced in areas surrounded by mountains. Short and long-term exposure to PM 10 in outdoor and indoor environments can result in respiratory and cardiovascular diseases placing an economic burden on society due to medical expenses or even death of experienced workers (Martins & da Graça, 2018) .

The PM fraction was obtained by subtracting the pre-sampling weight of filter from the post-sampling weight using an analytical microbalance (Model LIBROR AEL-40SM; SHIMADZU Co., Kyoto, Japan) with a 1 μg sensitivity. Since PAHs were not investigated in Asaluyeh where the weather is warm and humid, for more certainty in PAHs results, sampling campaign was conducted in winter over approximately a ten-day period (27/1/2018 to 6/2/2018). The results revealed that the PM concentration varies in different days depending on weather condition. For example on 28/1/2018, high wind speed increased resuspended dusts. However, PM concentration and temperature generally had a positive correlation. Furthermore, it is evident that increasing temperature decreases soil moisture and consequently increases resuspended dust (Figs. 4 and S2) . The maximum (6.25 μg/m³), minimum (0.72 μg/m³) and average (2.38 μg/m³) concentrations of PM10 did not exceed the USEPA (100 μg/m³) and WHO (50 μg/m³) air quality standards. Also, according to PM10 concentration, calculated air quality indices (AQIs) ranged from 0.66 to 5.78 with an average of 2.2 indicating "good" category. Fig. 4 shows variations of the PAH concentrations in the (re-)suspended dust samples collected daily in industrial and urban areas. The highest ∑PAHs is 4.73 ng/m 3 in N4 sampling site (a petrochemical company), while the minimum ∑PAHs was 2.84 ng/m 3 in U1 station (an urban area). The average air dust PAH content in the study area is 3.75 ng/m 3 with non-carcinogenetic PAHs being more abundant than carcinogenetic ones. The results indicate higher total PAHs and LMW PAHs concentrations in the industrial areas (petrochemical units) compared with urban areas. Among combustion-derived PAHs, HMW-PAHs are abundantly generated at high temperature (automobile motors), while the petroleum-derived residues have relatively high LMW- PAHs content (Mai et al., 2003) . Figs. 3 and 4 reveal that the decreasing trend of LMW PAHs concentrations is: air dust > winter street dust > summer street dust. This clearly shows that two major factors affecting LMW PAHs concentration include (i) temperature (ii) volatility. It is no surprise that the filter dust has relatively high levels of LMW. LMW (2-3 rings) PAHs are known to be relatively more volatile resulting in their less abundance in settled dust compared to air particles. Also, given the lower temperature in winter, higher LMW PAHs concentration was observed. PAH concentration and temperature have a negative relationship since PAHs (especially LMW PAHs) are temperature sensitive and with increasing temperature, PAHs volatilize. It is interesting to know that PM and PAHs concentrations have a negative correlation owing to dilution of the PAHs concentration by increased suspended dust concentration (Figs. 4 and S2) .

The skewness, standard deviation and Sig value (< 0.05) of Shapiro-Wilk test confirmed non-normality of PAHs (Table 1) . Thus, considering the non-normal distribution of PAH compounds, nonparametric statistical tests were chosen to conduct statistical analyses.

Factor analysis is an appropriate statistical method for determining the relationship between components. In this study, PCA as an extraction method, Varimax rotation for better representation of data and eigenvalues greater than 1 were used to extract factors. By using this statistical test, the number of variables can be reduced and presented in several main components, each containing several variables. The data were also analyzed by KMO test for evaluating the number of samples and Bartlett test for measuring similarity of data prior to performing factor analysis. The KMO and Bartlett values were more than 0.6 and lower than 0.3, respectively, while the total variance of the components was more than 70% of the total data variance. Hence, suitability of factor analysis was approved for reducing the number of variables and identifying hidden structures between the data. Regarding the results of PCA (Table 2) , the first component represented 67.6% of the total variance and included LMW PAHs (naphthalene, acenaphthene, fluorine, phenanthrene and anthracene) that mostly originate from oil leaks that evaporate and degrade rather easily Liu, Cheng, Ou et al., 2007) Liu, Cheng, Ou et al., 2007) . The presence of pyrene and chrysene in this group indicates that their origin is petroleum burning. Burning fossil fuels commonly leads to the formation of a large amount of pyrene, benzo[a]pyrene, benzo[b]fluoranthene and benzo[k]fluoranthene, all included in the second component (Hwang, Kusch, Barral, & Huffaker, 2003; Rajput & Lakhani, 2009). Benzo[g,h,i] perylene and indene have both similar molecular weight and tendency to be absorbed into dust particles. Generally, the MLR results indicated that the first and second components include 74.05% and 25.95% of the total PAHs concentration in Asaluyeh, respectively (Table 2) . A fact that reflects LMW PAHs being usually generated in petrochemical units in the industrial areas probably make the main proportion of PAHs in the study area.

The results of GEE analysis confirmed PCA-MLR analysis (Table 3) . According to the Sig and B of PAHs in GEE analysis, no strong positive correlation exists between LMW and HMW PAH subgroups reflecting their diff ;erent origins and physicochemical characteristics. LMW compounds have high atmospheric mobility and may have been transported from remote sites via atmospheric transportation, while Bold value signifies: P value < 0.05.

GEE matrix for PAHs concentration. HMW compounds are more likely associated with airborne particulates and only moved short distances (Augusto et al. 2010; Chung et al., 2007; Keshavarzi et al. 2017) . Also, strong positive correlations exist between members of each group. For example, a strong positive correlation exists between LMW members. According to Ravindra et al. 2006 ). According to Table 3 , none of the analyzed PAH species is related to TPH. By implementing the backward GEE model, it was indicated that the Nap, Chr, BbF and Ind have the greatest effect on TPH (Table 4) . For more accuracy, the backward GEE model was also implemented on Nap, Chr, BbF and Ind group against TPH ( Table 5 ). The results revealed Ind effect on the TPH levels as Ind is a six-ring HMW PAH with the molecular weight of 276. Thus, Ind is more stable in the environment and usually makes a fraction of TPH concentration; however, it does not necessarily mean that the rest of PAH species are ineffective.

To improve source identification of the PAHs, PMF model was run. The Asaluyeh street dust and air dust particles are potentially influenced by three sources of PAHs including, (i) pyrogenic, (ii) petrogenic, and (iii) geogenic sources (soils, dust, minerals and etc.). The PMF model was run 20 times for the three factors in which PAHs with a signal-to-noise (S/N) ratio between 0.2 and 2, and missing or values below an MDL greater than 50% were considered as weak variables. Fig. 5 indicates the source profiles of PMF model in the street dust and air dust, respectively. In terms of street dust, Benzo(k)fluoranthene, Benzo(a)pyrene, Dibenz(a,h)anthracene, Benzo(ghi)perylene and Indene are grouped in the first factor with high factor loading values which feasibly originated from urban sources and include HMW PAHs. HMW compounds have similar chemical and physical properties and originate from vehicle motors. In the second factor which can be considered as representing industrial sources (refinery and petrochemical units), most LMW compounds constitute the main PAHs. On the other hand, the third factor demonstrates characteristics of 4 ring PAHs which are similar to HMW and LMW PAHs and are classified as moderate molecular weight (MMW) PAHs. They probably originated from mixed pyrogenic and petrogenic sources.

Regarding air dust, three geogenic, pyrogenic and petrogenic factors were identified. The first factor includes PM (geogenic). The second and third factors include petrogenic and pyrogenic sources, respectively. HMW PAHs are mostly pyrogenic and originate from fossil fuel combustion (often in urban areas) and LMW PAHs are mostly petrogenic and originate from oil leaks (often petrochemical units in industrial areas). Also, the temperature is effective on all three factors, which is due to the effect of temperature on (re-)suspended dust generation and sensitiveness of PAHs to temperature Due to the fact that the transport and fates of PAHs of street dust and (re-)suspended dust in urban ecosystems such as Asaluyeh County should be identified, an estimation of the partitioning of PAHs between street dust and (re-)suspended dust may give beneficial insights. Generally, three scenarios may occur that include (i) direct PAHs attachment to air dust, (ii) street dust PAHs evaporate and then attach to air dust, (iii) resuspension of street dust as air dust. In the first two scenarios, LMW PAHs concentration is expected to be higher than HMW PAHs concentration. Hence, Fig. 6 reveals that the LMW PAHs concentration in the air dust is higher than street dust. The atmospheric LMW PAHs (2 and 3 rings) mainly occur in the vapour phase, while HMW PAHs (5 rings or more) are often limited to particles. Regarding atmospheric temperature MMW PAHs (4 rings) partition between the particulate and vapour phases (Srogi, 2007) . Therefore, according to Fig. 6 , PAH species display different transmission values in street dust Bold value signifies: P value < 0.05. 

Diagnostic ratios of PAH species are commonly used to distinguish pyrogenic and petrogenic sources of PAH compounds. These ratios can explain potential PAHs emission sources. Generally, LMW/HMW, Phe/ Ant, Ant/(Ant + Phe), Flu/(Flu + Pyr), BaA/(BaA + Chr) and Ind/ (Ind + BghiP) ratios are used for this purpose (Table S4) . LMW/HMW ratio < 1 indicates pyrogenic and > 1 shows a petrogenic source (Zhang, Dou et al., 2008) . Phe/Ant ratio > 15 and < 10 reveals petrogenic and pyrogenic sources, respectively (Takada, Onda, Harada, & Ogura, 1991) . Furthermore, Ant/(Ant + Phe) > 0.1 and < 0.1 are indicative of pyrogenic and petrogenic sources (Pies et al., 2008) . Similarly Flu/(Flu + Pyr) < 0.4 ratio indicate a petroleum source, while ratios between 0.4 and 0.6 indicate gasoline emissions, while ratios between 0.6 and 0.7 indicate diesel emissions sources (Yunkera et al., 2002) . If BaA/(BaA + Chr) is lower than < 0.2, it reveals petrogenic PAHs, and ranges from 0.22 to 0.55, 0.38 to 0.64 and 0.35 to 0.50 show gasoline emissions, diesel emissions and petrogenic origin, respectively (Simcik, Eisenreich, & Lioy, 1999; Sicre et al., 1987; Wang, Tian, Yang, Liu, & Yi-Fan, 2009; Yunkera et al., 2002) . Also, ratio of Ind/(Ind + BghiP) < 0.2; 0.2 to 0.5; and > 0.5 indicate petroleum and petrogenic; petroleum combustion; and coal, grass, and wood sources, respectively (Zhang, Luo, Wong, Zhao, & Zhang, 2006) . Table S4 presents the calculated ratios in this study. As it can be seen, PAHs in Asaluyeh street dust, generally originate from pyrogenic sources, but a closer look at each sampling site (Fig. S3) , reveals that PAHs sources (in both summer and winter) are mostly petrogenic in the industrial area and pyrogenic in urban areas, and probably originate from diesel and fossil fuel combustion and also automobiles, gasoline (Table S4 and Fig. S3 ). About (re-)suspended dust, all samples illustrated diesel emissions (except N2 sample that showed a mixed source). In the industrial areas, in addition to petrochemical units, the movement of vehicles is effective and the air pollution can move a long distance (Fig. S3) . For more certainty, logistic regression analysis was used in the current study. A benefit of logistic regression analysis compared with other statistical methods (such as t-test) is that it takes into account the number of polluted stations, while other statistical methods are based on the concentration of pollutants. In the current study, logistic regression analysis was used to compare pyrogenic and petrogenic sources of PAHs between industrial and urban areas (Table 6 ). It can be seen that there is a significant difference between industrial and urban areas in terms of PAHs origin (Sig = 0.04). Accordingly, Exp (B) indicated that petrogenic PAHs in industrial areas are 11.2 times more abundant than in urban areas. The reason is that in industrial areas, petrochemical and refining processes and consequently petrogenic sources for PAHs are more common. Also, combustion of fossil fuels by motor vehicles in urban areas produces pyrogenic PAHs.

Human health risk assessment was carried out assuming that adults and children are exposed to PAHs (based on residential/commercial and industrial scenarios) (Table S2) . Table 1 shows calculated toxic equivalency concentration (TEQ) for Asaluyeh street dust along with toxic equivalency factors (TEF) of PAH species. TEQ is used in the incremental lifetime cancer risk (ILCR) formula. Generally, the acceptable value for ILCR is equal to or lower than 10 −6 while ILCR between 10 -4 and 10 -6 indicates potential risk, and above 10 -4 , represents high health risk potential (US EPA, 2008) .

The results revealed that TEQ in summer is higher than winter (Fig. 2d) obviously as already mentioned PAHs concentration in summer is higher. Based on residential/commercial and industrial scenarios for dermal contact and ingestion exposure pathways, the average estimated ILCRs for both adults and children is more than 10 −4 and for inhalation pathway, lower than 10 -6 (Table 7) . Therefore, dermal and ingestion pathways are the main routes of exposure to PAHs and pose a high potential health risk, while exposure through inhalation pathway is acceptable ( Table 7) . The results also indicated that mean cancer risk is higher than 10 −4 and carcinogenic risk is high and worrying. Also, the results revealed that estimated ILCR for residential/ commercial areas is higher than in industrial areas and that in the industrial areas dermal and inhalation pathways for indoor workers are more hazardous. Finally, according to different groups and scenarios, children are exposed to higher health risk (Table 7) .

TPH and PAHs distribution in Asaluyeh County revealed that the reason for pollutants accumulation in the street dust is not confined to human activity. The streets enriched with dust are a dynamic environment and can change easily. For instance, human activity, precipitation, wind blow, street wash out and management plans all affect street dust. On the other hand, pollutants are affected by temperature, sunlight, human activity, technology level, remediation operation and several other factors. Therefore, all above-mentioned factors should be considered in environmental monitoring. However, in the current study in order to obtain more precise results, strong statistical methods such as GEE, backward GEE, logistic regression, PCA-MLR and PMF were used. Logistic regression demonstrated petrogenic sources for PAHs in the industrial and pyrogenic sources in the urban areas. Since TPH is used as an organic pollutants indicator including PAHs, backward GEE model was used and demonstrated that TPH is more influenced by HMW PAHs, especially Indene. Also, decreasing trend of LMW PAHs concentrations is: air dust > winter street dust > summer street dust. This clearly shows that two major factors affecting LMWPAHs concentration include (i) temperature (ii) volatility. Eventually, this study will help to design effective control and mitigation strategies to improve ambient air quality by providing information about source pollutions and create engineering mechanisms for air ventilation. Finally, the proposed statistical methods provide a valuable and reliable means of environmental evaluation. 

Fractionation, source identification and risk assessment of potentially toxic elements in street dust of the most important center for petrochemical products

Distribution and potential health impacts of microplastics and microrubbers in air and street dusts from Asaluyeh County

Investigation of microrubbers, microplastics and heavy metals in street dust: A study in Bushehr City

Determination of total petroleum hydrocarbons and heavy metals in soils within the vicinity of facilities handling refined petroleum products in Lagos Metropolis

Metal pollution in seaweed and related sediment of the Persian Gulf

Petroleum and individual polycyclic aromatic hydrocarbons

Chronic effects of oil spill on soil properties and microflora of a rainforest ecosystem in Nigeria

Presence of polycyclic aromatic hydrocarbons in sediments and surface water from Shadegan wetland-Iran: A focus on source apportionment, human and ecological risk assessment and sediment-water exchange

Lichens as an integrating tool for monitoring PAH atmospheric deposition: A comparison with soil, air and pine needles

Comparison of atmospheric polycyclic aromatic hydrocarbon levels in three urban areas in Lebanon

Goodness-of-fit tests for GEE modeling with binary responses

A statistical PCA method for face recognition

Wintertime size distribution of polycyclic aromatic hydrocarbons (PAHs), polychlorinated biphenyls (PCBs) and organochlorine pesticides (OCPs) in the urban environment: Street-vs rooftop-level measurements

Pollutants in Hong Kong soils: Polycyclic aromatic hydrocarbons

On consistency and inconsistency of estimating equations

On the use of a working correlation matrix in using generalized linear models for repeated measures

Sediment-water exchange, spatial variations, and ecological risk assessment of polycyclic aromatic hydrocarbons (PAHs) in the Songhua River

A comparative study of face recognition with principal component analysis and cross-correlation technique

Characteristics, toxicity, and source apportionment of polycylic aromatic hydrocarbons (PAHs) in road dust of Ulsan

Speciated hydrocarbon profiles and calculated reactivities of exhaust and evaporative emissions from 82 in-use lightduty Australian vehicles

National ambient air quality standards for particulate matter; Final rule. Environmental Protection Agency Part II 40 CFR Part 50

Comparison of population-averaged and subject-specific approaches for analyzing repeated binary outcomes Am

Polycyclic aromatic hydrocarbons study in Taichung

Using generalized estimating equations for longitudinal data analysis

Polycyclic aromatic hydrocarbons in sediments of marine coastal lagoons in Messina, Italy: extraction and GC/MS analysis, distribution and sources

Exposure and cancer risk assessment of polycyclic aromatic hydrocarbons (PAHs) in the street dust of Asansol city

Distribution, partitioning and sources of polycyclic aromatic hydrocarbons in Daliao River water system in dry season

Particulate matter in the atmosphere: Which particle properties are important for its effect on health

Spindle orientation in Saccharomyces cerevisiae depends on the transport of microtubule ends along polarized actin cables

Airborne particle number profiles, particle mass distributions and particle-bound PAH concentrations within the city environment of basel: An assessment as part of the BRISKA project

An overview of principal component analysis

The effects of road transport on freshwater and marine ecosystems

Polycyclic aromatic hydrocarbons in street dust of Bushehr City, Iran: Status, source, and human health risk assessment

Comparison of generalized estimating equations and quadratic inference functions in superior versus inferior Ahmed glaucoma valve implantation

PAH source fingerprints for coke ovens, diesel and, gasoline engines, highway tunnels, and wood combustion emissions

Distribution and sources of polycyclic aromatic hydrocarbons in size-diff ;erentiated re-suspended dust on building surfaces in an Oilfeld City

Characterization of particulate and vapor phase polycyclic aromatic hydrocarbons in indoor and outdoor air of primary schools

Source apportionment of polycyclic aromatic hydrocarbons in the urban atmosphere: A comparison of three methods

Longitudinal data analysis using generalized linear models

Pollution characteristics and ecological risk of polycyclic aromatic hydrocarbons (PAHs) in surface sediments of the southern part of the Haihe River system in China

Characterization, identification of road dust PAHs in central Shanghai areas

Characterization, identification of road dust PAHs in Central Shanghai Areas

Combustion-derived polycyclic aromatic hydrocarbons in the environment-a review

Distribution of polycyclic aromatic hydrocarbons in the Coastal Region off ; Macao, China: Assessment of input sources and transport pathways using compositional analysis

Polycyclic aromatic hydrocarbons (PAHs) in indoor dust matter of Palermo (Italy) area: Extraction, GC-MS analysis, distribution and sources

Chemical characterization of particulate emissions from diesel engines: A review

Impact of PM2.5 in indoor urban environments: A review

A review on polycyclic aromatic hydrocarbon (PAH) emissions from energy generation

Policies and technical guidelines for urban planning of high-density cities -air ventilation assessment (AVA) of Hong Kong

Key action 4, city of tomorrow and cultural heritage, from the program energy, environment and sustainable development, within the fifth framework program of the EU

Toxic equivalency factors (TEFs) for polycyclic aromatic hydrocarbons (PAHs)

A ground-based method of assessing urban forest structure and ecosystem services

Polycyclic aromatic hydrocarbons: Evaluation of sources and effects. Committee on pyrene and selected analogues, board on toxicology and environmental health hazard, commission on life sciences

Street design and urban canopy layer climate

An introduction to generalized estimating equations and an application to assess selectivity effects in a longitudinal study on very old individuals

Overview of the face recognition grand challenge. Computer vision and pattern recognition947-954 CVPR

Polycyclic aromatic hydrocarbons levels in street dust in the central area of Niteroi city

Polycyclic aromatic hydrocarbons and nitropolycyclic aromatic hydrocarbons in particulates emitted by motorcycles

Characterization and source identification of polycyclic aromatic hydrocarbons (PAHs) in River Bank Soils

Deletion diagnosis for generalized estimating equations

Robust regression for clustered data with application to binary responses

Assessing robustness of generalized estimating equations and quadratic inference functions

Measurements of polycyclic aromatic hydrocarbons at an industrial site in India

Seasonal and site-specific variation in vapour and aerosol phase PAHs over flanders (Belgium) and their relation with anthropogenic activities

Sources of fine organic aerosol. 3. Road dust, tire debris, and organometallic brake lining dust: Roads as sources and sinks

Aliphatic and aromatic hydrocarbons in different sized aerosols over the Mediterranean Sea: Occurrence and origin

Source apportionment and source/ sink relationships of PAHs in the Coastal Atmosphere of Chicago and Lake Michigan

Seasonal occurrence and cancer risk assessment of polycyclic aromatic hydrocarbons in street dust from the Novi Sad city, Serbia. The Science of the Total Environment

Monitoring of environmental exposure to polycyclic aromatic hydrocarbons: A review

Distribution and sources of polycyclic aromatic hydrocarbons (PAHs) in street dust from the Tokyo Metropolitan area. The Science of the Total Environment

Polycyclic aromatic hydrocarbons study in atmospheric fine and coarse particles using diagnostic ratios and receptor model in urban/industrial region

Sources and properties of non-exhaust particulate matter from road traffic: A review. The Science of the Total Environment

Toxicological profile for total petroleum hydrocarbons (TPH)

Risk assessment guidance for superfund. human health evaluation manual

Atlas asbestos mine superfund site. U.S. EPA releases exposure and risk assessment for clear creek management area

Benzo[a]pyrene (BaP) CASRN 50-32-8

Risk assessment guidance for superfund, volume 1, human health evaluation manual (Part B, development of risk-based preliminary remediation goals

Determination of polycyclic aromatic hydrocarbons (PAHs) in ambient air using gas chromatography/mass spectrometry (GC/MS). compendium method TO-13ACompendium of methods for the determination of toxic organic compounds in ambient air center for environmental research information

Analyses of water samples from the Deepwater Horizon oil spill: Documentation of the subsurface plume. Monitoring and Modeling the deepwater horizon oil spill: A

A new speciation scheme of soil polycyclic aromatic hydrocarbons for risk assessment

Application of positive matrix factorization to identify potential sources of PAHs in soil of Dalian

Polycyclic aromatic hydrocarbons (PAHs) in urban surface dust of Guangzhou, China: Status, sources and human health risk assessment. The Science of the Total Environment

Health risk of particulate matter from long-range transboundary air pollution. Joint WHO/Convention Task Force on the health aspects

Biological sources of polycyclic aromatic hydrocarbons (PAHs) in the Amazonian Rain Forest

Urban environmental geochemistry of trace metals

Multivariate analysis of heavy metal contamination in urban dusts of Xi'an, Central China. The Science of the Total Environment

PAHs in the Fraser River Basin: A critical appraisal of PAH ratios as indicators of PAH source and composition

PAHs in the Fraser River Basin: A critical appraisal of PAH ratios as indicators of PAH source and composition

Distributions and concentrations of PAHs in Hong Kong soils

Source diagnostics of polycyclic aromatic hydrocarbons in urban road runoff, dust, rain and canopy throughfall

Emission of polycyclic aromatic hydrocarbons from indoor straw burning and emission inventory updating in China

Characterizations of resuspended dust in six cities of North China

Analysis of correlation structures using generalized estimating equation approach for longitudinal binary data

Rapid assessment of exposure to particulate emissions from surface contamination sites

Preliminary assessment of heavy metal contamination in water and wastewater from Asaluyeh Port (Persian Gulf)

A comparative analysis of several vehicle emission models for road freight transportation

Geochemistry and risk assessment of street dust in Luanda, Angola: A tropical urban environment

Contamination level, source identification and risk assessment of potentially toxic elements (PTEs) and polycyclic aromatic hydrocarbons (PAHs) in street dust of an important

Chemical speciation, human health risk assessment and pollution level of selected heavy metals in urban street dust of Shiraz

Distribution, source identification and health risk assessment of soil heavy metals in urban areas of Isfahan province

Prepared for EPA office of health and environmental assessment

Heavy metals and polycyclic aromatic hydrocarbons: Pollution and ecological risk assessment in street dust of Tehran

Ecological and human health hazards of heavy metals and polycyclic aromatic hydrocarbons (PAHs) in road dust of Isfahan metropolis, Iran. The Science of the Total Environment

Development of an improved dry and wet deposition collector and the atmospheric deposition of PAHs onto Ulsan Bay

Supplemental guidance for developing soil screening levels for superfund sites

Supplemental guidance for developing soil screening levels for superfund sites

Exposure factors handbook

Human exposure to soil contamination: A qualitative and quantitative

Characteristics of PAHs adsorbed on street dust and the correlation with specific surface area and TOC

Evaluating bioavailability of organic pollutants in soils by sequential ultrasonic extraction procedure

Metals associated with stormwater-relevant brake and tire samples

This research was financially supported by Bushehr Environmental Protection Office. The authors also wish to express their gratitude to the Research Committee and Medical Geology Center of Shiraz University for logistic and technical assistance. It was an honor to work with Professor Farid Moore who takes his time for editing our manuscript.

Supplementary material related to this article can be found, in the online version, at doi:https://doi.org/10. 1016/j.scs.2018.11.015.