Next Article in Journal
Assessing the Impact of Cumulus Parameterization Schemes on Simulated Summer Wind Speed over Mainland China
Next Article in Special Issue
The Contents of Potentially Toxic Elements and Emission Characteristics of PM2.5 in Soil Fugitive Dust around Six Cities of the Yunnan-Guizhou Plateau in China
Previous Article in Journal
Insight on Poleward Moisture and Energy Transport into the Arctic from ERA5
Previous Article in Special Issue
Genotoxicity of PM2.5 and PM1.0 Particulates on Human Peripheral Blood Lymphocytes in Manila, Philippines
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Population Health Risks Assessment from Air Pollution Exposure in an Industrialized Residential Area in Greece

by
Michail Pachoulis
1,
Thomas Maggos
1,
Panagiotis Panagopoulos
1,
Maria Dasopoulou
1,
Dimitra Balla
1,
Asimina Stamatelopoulou
1,
Manousos Ioannis Manousakas
2,3,
Konstantinos Eleftheriadis
2 and
Dikaia Saraga
1,*
1
Atmospheric Chemistry & Innovative Technologies Laboratory, NCSR Demokritos, INRASTES, Aghia Paraskevi, 15310 Athens, Greece
2
Environmental Radioactivity Laboratory, NCSR Demokritos, INRASTES, Aghia Paraskevi, 15310 Athens, Greece
3
Laboratory of Atmospheric Chemistry, Paul Scherrer Institute (PSI), 5232 Villigen, Switzerland
*
Author to whom correspondence should be addressed.
Atmosphere 2022, 13(4), 615; https://doi.org/10.3390/atmos13040615
Submission received: 2 March 2022 / Revised: 2 April 2022 / Accepted: 6 April 2022 / Published: 11 April 2022
(This article belongs to the Special Issue Health Effects of Urban Atmospheric Aerosols)

Abstract

:
Industrial activities nearby residential areas lead to poor local air quality. Therefore, short-term exposure to an aggravated environment and the subsequent health effects should be the subject of further research. The purpose of this study is to estimate the health risks resulting from such exposure in population groups living in an industrialized area. The risk estimation was performed using different approaches suggested in relative literature. Monitoring of the air quality in an industrial zone of Attica was carried out including 24-h measurements of PM2.5 and analysis of their chemical composition for Polycyclic Aromatic Hydrocarbons and heavy metals (Pb, Cd, As, Ni, Hg, Cu, Zn). Samples of Volatile Organic Compounds were also collected. Health effects on different population subgroups were estimated for the targeted pollutants through different mathematical approaches provided by the literature, taking into consideration different parameters (e.g., age, gender, exposure duration). Inhalation rate and body weight were important parameters to estimate the exposure dose of people, and they can vary greatly depending on the age, gender, and daily activity of the person under consideration. The results indicated that the risk for potential carcinogenic and non-carcinogenic effects varies depending on the applied methodology. In any case, the acceptable limits for cancer risk provided by the OEHHA, EPA, and WHO were not exceeded.

1. Introduction

Air quality deterioration has become a serious matter of concern due to increased anthropogenic and natural emissions [1,2], leading to increasing cases of acute air pollution episodes and exceedances of the air quality standards on a global scale [3,4]. Public health has been seriously affected by air pollution during the last decades, and the problem is expected to intensify in the future. According to the guidelines of the World Health Organization (WHO), seven million people die early from air pollution every year worldwide, while nine out of 10 people breathe air containing high levels of pollutants [5,6]. Nowadays, there is much more evidence to support these claims, as well as better our understanding of the way the air pollutants affect public health even in lower concentrations, as it is estimated that the number of deaths and the loss of healthy years of life due to air pollution has not been reduced [7].
Suspended Particulate Matter (PM) is considered to be one of the most serious pollutants and has been classified as a ‘Group 1’ contaminant (according to the International Agency for Research on Cancer, Group 1 includes substances that have sufficient evidence of carcinogenicity in humans, while Group 2 includes substances for which the degree of evidence of carcinogenicity in humans is almost sufficient) [8,9].
Consequently, an increasing number of studies have focused on providing data about air quality deterioration and the resulting health effects on the population. This relation has been extensively documented [10,11] in studies examining indoor and outdoor exposure to hazardous substances in children and adults (e.g., PAHs & heavy metals) [12,13,14,15,16]. The main target of those studies was to assess the exposure to ambient PM associated with increased mortality and morbidity, bio-accessibility through inhalation exposure and lifetime lung cancer risk [17,18], disability-adjusted life years (DALYs), cardiopulmonary diseases, and risk factors including a wide variety of health effects for people of all ages [4,19,20]. However, there is still uncertainty about the mechanisms through which air pollutants influence human health; for example, it is not clear if health implications are the result of synergistic or individual effects of pollutants, which makes the assessment a complicated process.
The current study aimed to assess the risk of exposure of a population of an industrialized residential area (Elefsina, Greece) to atmospheric pollutants originating mainly from industrial activities. For this scope, estimation of the risk of carcinogenic and non-carcinogenic effects in the human body was carried out, based on a number of different mathematical approaches reported in the literature. A range of impact assessments was achieved and a comparison of their strengths and weaknesses was attempted with the scope of contributing to a better understanding of the population health risk caused by air pollution. The effect of considering various parameters (i.e., age, gender, potential cancer factor, daily exposure) on the risk estimation outcome, was also investigated. To the best of our knowledge, this is the first time that a comparative study on the population health risks assessment from air pollution exposure in an industrialized residential area has been conducted.

2. Materials and Methods

2.1. Location Description

An integrated sampling campaign was carried out at three sites in the wider area of Elefsina (Figure 1 and Figure 2; Table 1) during the winter and summer period of 2019 (February–October). Elefsina is a municipality of the West Attica Prefecture, located at a distance of about 20 km from the center of Athens, Greece. It covers an area of 20 sq km and has a population of 24,910 inhabitants (2011 census). The area is characterized by intense industrial activity (i.e., oil refinery). The sampling sites selection followed EU Directive 2008/50/EC while taking into consideration the characteristics of the surroundings of each area. All sites were selected to be located within the residential area of Elefsina. More specifically, E1 and E2 were located close to the industrial zone (refiners) while the E3 site was near the port (Figure 2).
All sampling and monitoring equipment was installed in schoolyards. For all three locations, meteorological parameters (wind speed and direction) were continuously recorded using a portable anemometer. During the sampling period, the prevailing wind direction was W-NW with wind speed ranging from 0.70 to 4.07 m/s while temperature ranged from 5 to 21 °C (February–May) and 20 to 31 °C (June–October). The relative humidity ranged between 42 and 88%.

2.2. Sampling and Chemical Analysis

The 24-h (starting at 8 a.m.) PM2.5 samples (n = 180) were collected on 47 mm Tissue Quartz 25000QAO PALL membrane filters with the use of low volume (2.3 m3/h) samplers (Derenda, Leckel). The determination of the PM mass was conducted according to EN 12,341:2014. After PM mass concentration determination, PM samples were analyzed for 22 PAHs (naphthalene, acenaphthylene, acenaphthene, fluorene, phenanthrene, anthracene, fluoranthene, pyrene, benzo(a)anthracene, chrysene, benzo(b)fluoranthene, benzo(k)fluoranthene, benzo(a)pyrene, indeno(123-cd)pyrene, dibenzo(ah)anthracene, benzo(ghi)perylene, benzo(b)fluorene, benzo(e)pyrene, perylene, 1,2-dimethylnaphthalene, 1-methylnaphthalene, 1-methylphenanthrene, 2-methylnaphthalene) and seven heavy metals (Pb, Hg, Ni, Cr, As, Cd & Be). The determination of heavy metals was based on EN14902:2005 European Standard using an ED-XRF (Energy Dispersive X-Ray Fluorescence) Epsilon 5, Panalytical, analyzer for the analysis of Pb, Hg, Ni, Cr while for As, Cd, Be an Atomic Absorption Spectrometer (AAS) Varian 220, GTA 110 was used. The sampling and the analysis of the PAHs were performed according to EN15549:2008 using a gas chromatograph (GC Agilent 5975C) coupled to a mass spectrometry detector (Agilent 7890A MS. The values of LOD, LOQ, as well as the expanded uncertainty at 95% confidence level and k = 2 that has been estimated experimentally (Uexp) for each PAH, are summarized in the Supplementary Material, Table S1. The analytical procedures are described elsewhere [16,21].
In parallel, 330 samples of VOCs were collected during the summer and winter periods of 2019 (February–October), covering the factor of seasonal fluctuations in their concentration levels (Table 2). Samples were analyzed using a gas chromatography system (Agilent GC6890) coupled with an FID detector and a Thermal Desorption System (Gerstel TDSA) according to EN 16017:2001 (Table S2). Samples were collected in preconditioned glass tubes filled with Tenax TA (Gerstel) and were analyzed within 1 day after sampling. LOD for VOCs (benzene) analysis was 0.02 ng/l while the expanded uncertainty (% Uexp, 95%, k = 2) of the analytical procedure was 11.7%.

2.3. Health Risk Assessment (Methodological Approach)

The process to determine if a substance, chemical or not, constitutes a potential risk to human health is a complex procedure. The dose-response assessment refers to the characterization of the correlation between exposure to an agent and the incidence of an adverse health effect in the exposed population. In the present study, a literature review on methods used to evaluate these potential health risks through different mathematical approaches was conducted. Based on this review, different hypothetical scenarios of human exposure were selected to be applied in the case of Elefsina’s population. Specifically, the following representative cases were selected
The case of a male and a female living and working at the place of exposure, i.e., 24 h/day.
The case of a male and a female living at the place of exposure (their place of residence) for 14 h/day, and working outside of that region,
The cases of a child in the age groups of 0–2 and 2–16 years old, living at the place of exposure, i.e., 24 h/day.
The exposure period was selected to be 350 days/year in all cases. Furthermore, specific values of the parameters IR and BW (Inhalation Rate & Body Weight) were selected according to literature [22,23], for the residents of Elefsina. Particularly:
  • For males: IR = 16.4 (m3/day) and BW = 76 (kg),
  • For females: IR = 12.6 (m3/day) and BW = 63 (kg),
  • For children of age 2–16 years old: IR = 10.8 (m3/day) and BW = 32.5 (kg),
  • For children of age 0–2 years old: IR = 4.9 (m3/day) and BW = 10.3 (kg).
  • For the calculation of the examined parameters, the average concentration of the considered pollutants given in Table 3 was used.

2.3.1. Carcinogenic Risk Assessment Methodological Approaches (Group A)

The probability of developing cancer over a lifetime in the people who live in the Elefsina area was estimated based on methods reported in the literature. In quantitative carcinogenic risk assessment, the dose-response correlation is expressed in terms of a potency slope, which is used to calculate the probability of carcinogenic risk associated with an estimated exposure (95th percent upper confidence limit of the slope of the dose response curve). The cancer risk (CR) value below 10−6–10−4 is considered acceptable [13,14,24,25,26,27], whereas 10−6 is considered the most tolerable risk [28]. In this study, CR characterization was conducted by applying the methodologies described in the following paragraphs (A.I–A.IV). As a final step, the cancer risk of each pollutant was summed up for the calculation of the Total Cancer Risk and then the overall Risk of Cancer was converted to “chances per million”
(Total Cancer Risk) (1 × 106) = Total Cancer Risk in chances per million
Approaches used for the estimation of potential cancer risk

Methodology A.I

In the study of R.M. Maertens (2008) [14] a risk assessment was conducted to evaluate the excess lifetime cancer risk associated with exposure to seven PAHs (Table S3), using the method proposed by [12] Collins et al., 1998. To be more specific, the following equation was used to estimate the lifetime cancer risk for these seven PAHs Equation (1):
Lifetime   Cancer   Risk = i = 1 n C i × PEF i × IR × EF × SF × AF BW × 1000
where C: concentration (µg/g) of each PAH, which are categorized as probable human carcinogens (B2) according to U.S. EPA classifications [14], PEFi: the potency of each PAH in correlation to benzo[a]pyrene, according to Collins et al., 1998 [12] (Potency equivalency factor), SF: Slope factor for carcinogenic (mg kg−1 day−1)−1, indicating the probability of cancer occurring per unit of inhaled PAHs over a lifetime. Maertens (2008) examined the cancer risk in association with ingestion instead of inhalation, which differentiates the values of SF, IR and BW, so were selected according to Handbook “Guidance Manual for Preparation of Health Risk Assessment” from the OEHHA, 2015 [24]. The PEF value remains the same, while IR (m3/day) was the same as mentioned previously (M: 16.4, W: 12.6, Child2–16: 10.8 & Child0–2: 4.9), as well as for BW (Kg), respectively, (M: 76, W: 63, Child2–16: 32.5 & Child0–2: 10.3), EF: exposure factor, the average proportion of daily exposure of people through inhalation rate (h/day), AF: adjustment factor for exposure dose modification in early life stages, as young children are more vulnerable to the effects of chemical toxins according to EPA, 2003 & 2005 [29,30].
Specifically:
  • AF0–2: 10-fold adjustment, for children between 0–2 years old
  • AF2–16: 3-fold adjustment, for children between 3–16 years old
  • AF16–80: no adjustment needed for children >16 years old and adults

Methodology A.II

A previous study by R.M. Maertens (2004) [13], also evaluating the excess lifetime cancer risk caused by exposure to the specific PAHs, applied the same methodology without including the adjustment factor (AF). The evaluation of cancer risk was estimated in two steps, as described by Equations (2) and (3):
Lifetime   average   daily   exposure   dose = C   ×   IR   ×   EF BW  
Lifetime   Cancer   Risk = Lifetime   average   daily   exposure   dose   mg / kg / day × Slope   Factor   mg / kg / day 1 × Potency   Equivalency   Factor
The remaining parameters were selected similarly to the methodology A.I.

Methodology A.III

In the study of M.M. Jackson (2005) [31], the contribution of road traffic to the air pollution levels in Dar-es-Salaam city of Tasmania was examined; a risk assessment was conducted for the people who were in the nearby area. In the present study, the potential cancer risk for the residents in the region Elefsina was estimated for a variety of pollutants, including some heavy metals (Pb & Ni). In addition, according to U.S. EPA risk guidelines from 1986 [32], Pb has been classified, with a group B2 classification, as a potentially carcinogenic substance. The risk for cancer (R) was calculated from the following equations Equations (4) and (5):
CDI = CA × IR × ED × EF × L BW × AT × 365
R =   CDI × P
where P is the cancer potency for the pollutants (mg·kg−1·day−1)−1, that were selected according to US-EPA Regional Screening Level (RSL) Summary Supplementary 2019 [33], except for Ni and Pb. The p-value for Pb was set equal to 4.2 × 10−2, according to the study of M. Jackson et al., 2005 [31] complying with OEHHA, 2015 (Appendices) [34]. Regarding Ni, it was set equal to 9.1 × 10−1 [34]. The CDI (chronic daily intake) dose (mg·kg−1·day−1) was calculated from Equation (4), with CA: the concentration of pollutant in the air (mg/m3), IR the inhalation rate (m3/h), EF the exposure frequency (days/year), ED the exposure duration per day (h/day), L the length of the exposure (years), BW the average body weight (kg) and AT the averaging time of 80 years length for carcinogens (ICRP, 2002) [22].

Methodology A.IV

Gao, (2019) [5] and Hong et al., 2020 [17], evaluated the potential cancer risk during inhalation for nine PAHs detected in PM2.5. In Gao’s 2019 study, the bioaccessibility of these PAHs was investigated, by employing a physiological extraction test with simulated lung fluids [Gamble’s solution and artificial lysosomal fluid (ALF)]. With the term bioaccessibility (%), we refer to the proportion of the pollutants which are important contributors to the effects on the human body. In the present study, this practical application was not feasible. The estimated risk for respiratory cancer (CR) was evaluated through the following Equation (6):
Cancer   Risk = i = 1 n C PAHi × RPF i × UR BaP
where CPAHi: the concentration of each PAH detected in PM2.5, RPF: the relative potency factor of each PAH can be determined in correlation to benzo[a]pyrene, according to the U.S. EPA (2010) Integrated Risk Information System [35] (Table 4), which can be found in the supplementary data of Gao et al. 2019 [5], URBaP (Unit Risk) is the probability of the maximum theoretical limit of the number of people with cancer in the respiratory system caused by inhalation at an equivalent concentration of 1 μg/m3 BaP over a lifetime of 70 years. Two different URBaP values were selected for the inhalation cancer risk assessment according to the California Environmental Protection Agency (CalEPA) (1.10 × 10−6 per ng/m3 based on the data for respiratory tract tumors from inhalation exposure in hamsters) and the World Health Organization (8.70 × 10−5 per ng/m3 based on an epidemiological study on coke-oven workers in Pennsylvania) (OEHHA, 1993, 2005; WHO, 2000, 2010) [36,37,38,39].
It is clear from the description of the different methodologies that there are key differences in the approaches for the evaluation of the potential cancer risk. The first three approaches include common parameters (PEF, SF & P). However, the A.I and A.II differ in the adjustment factor (AF), which modifies the estimated cancer risk for children as a function of age. In the A.III method, the parameters of Daily exposure (hours/day), exposure duration (days per year), and the length of the exposure (years) of how many years the exposure took place were added. On the contrary, the A.IV methodology presents significant differences from the others as exposure time, AF, IR and BW parameters are not included, while UR parameter is included.

2.3.2. Risk Assessment for Carcinogenesis and non-Cancer Effects Approaches (Group B)

Characterization of the risk of developing cancer and non-cancer effects over a lifetime was estimated for people in Elefsina, using the methods found in our literature review.
Cancer risk (CR) with a value below 10−6–10−4 is considered acceptable, as already mentioned. The non-cancer risk was estimated according to U.S.EPA (2004, 2009) [40,41], is expressed as the hazard quotient (HQ) and the hazard index (HI). HQ is calculated for each pollutant and HI corresponds to the sum of the individual HQs [41]. The acceptable values for HI and HQ are considered to be lower than unity (HI & HQ < 1), meaning that the greater the HQ and HI values are, the higher the probability of developing non-cancer effects in humans (U.S. EPA, 1989) [42]. In the final step of risk assessment, modeled concentrations and exposure data, which are determined through exposure evaluation and are combined with potency factors and Reference Exposure Levels (REL’s) are developed through a dose-response curve estimated assuming continuous lifetime exposure to substances [40,41,42,43].
The following methodologies were applied in order to estimate the potential cancer risk and non-cancer effects.

Methodology B.I

Using Air Toxics Hot Spots Program Guidance Manual for Preparation of Health Risk Assessment (Guidance Manual) from the Office of Environmental Health Hazard Assessment [24,34] (OEHHA, 2015; OEHHA Appendices, 2015), the evaluation of potential inhalation cancer risk and non-cancer acute hazard quotient (HQ) and hazard index (HI) was carried out. The risk of cancer was calculated from the following Equation (8) and the Dose through inhalation (mg/kg/d) from Equation (7):
Dose air = Cair BR BW A EF 1 × 10 6
Cancer   Risk = InhalationDose mg kg day Cancer   Potence kg day mg ASF FAH ED years AT years
where Cair: concentration of pollutant (μg/m3), BR/BW: daily breathing rate adjusted with body weight (L/kg BW-day), A: inhalation absorption factor (unitless, A = 1 according to OEHHA, 2015) [24], EF: exposure frequency (unit less, days/365 days), and 1 × 10−6: conversion micrograms to milligrams and liters to cubic meters, Cancer Potency was the cancer potency of a pollutant (mg·kg−1·day−1)−1, and was selected according to the OEHHA, 2015 [24] (Table S5), the ASF and FAH were variables that were used only when estimating residential cancer risk, so they were not applied, ED: the years of exposure duration and AT: averaging time period over which exposure duration (always 70 years) is averaged, (BR and BW for adults have used the averages values of men and women combined: 14.5 and 69, respectively).
Risk assessment for non-cancer Acute Hazard Indices was estimated. Firstly, the HQAcute for each pollutant was determined (9) with the 1-h maximum corresponding concentration for each pollutant and the acute reference exposure level (REL) for each pollutant correlated with the target organ system(s), according to OEHHA (2015) [24] (Table S6). Then, HQAcute for each pollutant was summed up for the calculation of the cumulative HI.
Acute   Hazard   Quotient = Maximum 1   h   Concentration Acute   REL

Methodology B.II

In the study by Tianjie Shao et al., 2018 [15], an assessment method of human exposure risk proposed by the U.S. EPA (Supplemental Guidance for Developing Soil Screening Levels for Superfund Sites EPA, 2002; L. Ferreira-Baptista et al., 2005) [44,45], was implemented, to estimate carcinogenic and non—carcinogenic exposure risk. The calculation of the potential cancer risk (Risk) was estimated by using the following Equations (10) and (11):
Risk = LADD × SF
LADD inh = C × EF PEF × AT × InhR child × ED child BW child + InhR adult × ED adult BW adult
where SF: the slope cancer factor for each pollutant, LADDinh: refers to the daily average exposure for life through the inhalation, C: concentration of pollutant (mg/kg), EF: exposure frequency (days/year), PEF: Particulate emission factor 1.36 × 10−9 (m3/kg). “This factor represents an estimate of the relationship between soil contaminant concentrations and the concentration of these contaminants in the air as a consequence of particle suspension” [44], AT: mean exposure time (70 years × 365 days), ED: exposure time (years), InhR: Inhalation rate (m3·day−1) and BW: weight per citizen (kg). The non-carcinogenic risk was evaluated from Equations (12)–(14):
HI = HQ i
HQ = ADD RfD
ADD inh = InhR × EF × ED PEF × BW × AT
where HQ is a non-carcinogenic risk factor that characterizes the non-carcinogenic risk of a single contaminant, ADDinh is a non-carcinogenic risk factor of a single contaminant from a pathway inhalation (mg·kg−1·day−1), RfD is the reference dose for the pathway inhalation (mg·kg−1·day−1), AT: mean exposure time (ED years × 365 days) and the remaining parameters were the same as those in Equation (11). In this study, the RfD was selected according to EPA, 2002 [44] and non-cancer effects were estimated only for Ni because there were no reference values for other pollutants.

Methodology B.III

In the study of Chalvatzaki et al., (2019) [26], health risk indices caused by inhalation were estimated. Specifically, the cancer risk (CR), non–cancer effects (HQ), and mortality cases (RR) from PM10 and PM2.5 were calculated. The probability of developing cancer risk (CR) was calculated by using Equation (15):
CR = CDI × CSF
where CSF is the cancer slope factor (mg/kg day)−1 and CDI is the chronic daily intake (mg/kg/day) which were calculated from the following Equations (16) and (17), respectively:
CDI = C a × IR × ET × EF × ED BW × AT
CSF = IUR × BW IR d × 10 3
Ca is the contaminant concentration (mg/m3), IR is the inhalation rate (m3/h), ET is the exposure time (h/day), EF is the exposure frequency (days/year), ED is the exposure duration (years), BW is the body weight (kg) and AT is the averaging time (70 years × 365 days/year), IUR is the inhalation unit risk for each pollutant according to U.S.EPA Regional Screening Level (RSL) Summary Table (2019) [33], except for the case of Pb, where the value IURPb: 1.2 × 10−5 proposed by U.S.EPA Regional Screening Level (RSL) Summary Table (2017) [46] was used, and IRd is the daily inhalation rate (m3/day). Finally, the number of new cancer cases (I) per lifetime (80 years) in Elefsina was estimated using Equation (18):
I = N × CR
CR is the cancer risk probability estimated in this work and N is the number of people in the target city (Elefsina), which according to the latest population census is 29,000 residents. The non-cancer risk assessment was performed using Equations (19) and (20):
HQ = CDI / RFD
RFD = RFC × IR d BW
where RFD is the reference dose (mg/kg/day), CDI was calculated as in the case of cancer risk Equation (16) with the only difference being the AT calculation (ED years × 365 days/year), RFC is the reference concentration according to Summary Table (2019) [33], except for the case of Pb, where the value RFCPb: 2 × 10−4 by the U.S. EPA Regional Screening Level Summary Table (2017) [46] was used (Table S7).
The methodology for the estimation of mortality cases was performed according to Ostro et al., 2004 [19]. Relative Risk (RR) for all-cause mortality, for all ages from PM10, was estimated by (21) & (22):
RR = exp β X X 0
AF = RR 1 RR
where X is the annual mean concentration of PM10 (µg/m3), and X0 is the baseline concentration of PM10 (µg/m3). For the annual mean concentrations (X), the average concentration of PM2.5 from E1, E2 and E3 was used (X = 14.4 µg/m3) while for the baseline concentration, the average concentration of PM2.5 from E3 (X0 = 3.2 µg/m3) was taken into account. Due to the lack of PM10 measurements, values of the PM10/PM2.5 ratio from the literature were used (for the region of Europe is 0.73; Ostro et al., 2004; Cardaba et al., 2014) [19,47]. Thus, the concentrations for PM10 were Χ = 19.3 and X0 = 10 μg/m3, β is the coefficient of the risk function (0.0008; 95% confidence interval (CI):0.0006–0.0010), AF is the attributable fraction (Equation (22)) which was used to estimate the proportion of deaths from a disease (e.g., lung cancer), and which could have been avoided if PM concentrations were reduced to background concentration, (i.e., concentrations that would exist without any human activities) [19,26]. Finally, the number of attributable deaths (AI) was calculated using Equation (23):
AI = AF × I
where I is the total number of deaths in the target population; for the region of Elefsina, this was equal to 222 for the year 2019.

Methodology B.IV

In the study of Megido et al., (2017) [27], a different approach for the assessment of human exposure risk (cancer and non—cancer effects) was proposed (U.S. EPA, 2003, 2009) [35,41]. The determination of cancer risk through inhalation (CRinh) was conducted from the following Equations (24) and (25):
CR inh = IUR × EC inh
where IUR is the Inhalation Unit Risk (µg/m3)−1 [33], while ECinh is the exposure concentration through inhalation (µg/m3), which was estimated by using the following Equation (25):
EC inh = C × ET × EF × ED AT n
C is the concentration of pollutant (μg/m3), ET is the exposure time (h/day), EF is the exposure frequency through the year (days/year), ED is the exposure duration (years) and ATn is the averaging time (70 years × 365 days/year).
The non-cancer risk assessment was evaluated through the hazard quotient (HQinh):
HQ inh = EC inh RfC i × 1000 × μ g mg
RfCi is the inhalation reference concentration (mg/m3) for each pollutant [33], and ECinh is the exposure concentration through inhalation as estimated from Equation (25) with the only difference of changing the calculation of parameter AT (ED years × 365 days/year).
There are key differences in the above-mentioned approaches (group B) for the risk assessment for carcinogenesis and non-cancer effects. The B.I methodology differs (a) on the values used for the parameters BR, BW, ED, and CR (which were set according to OEHHA) and (b) on the estimation of the cancer risk, as no differentiation in the gender, as well as adjustment of the daily exposure (hours/day), occurs. Additionally, in methodologies B.II, B.III, and B.IV, the parameters for the cancer potency factor and non-cancer factor are different. Method B.I is the only one in which the estimation of the Acute Hazard Quotient for each Target Organ System(s) is feasible. Contrary to B.I, the BII method includes alternative types (LADDinh, PEF, RfD), but lacks the distinction of age and gender. The B.III method is very similar to the A.III methodology for cancer risk assessment, while the only parameter missing in both is the AF parameter. In addition, with this methodology (B.III), the number of new cases of cancer for humans (I) per lifetime (80 years) can be estimated. Contrary to B.III, the parameters AF, daily exposure (IR), body weight (BW), and people’s gender are missing from the B.IV method. However, the outcomes from the two methods are similar (Table 5 and Table 6), possibly due to erroneous validation of these parameters in the B.III method.

2.3.3. Risk Assessment Approach Based on Methodologies Combination

In the present study, a combination of methodologies, which were presented in the previous paragraphs, was attempted in order to apply a different approach for estimating the risk of developing cancer over a lifetime for the residents of the Elefsina area (Table S4).

Methodology C.I

From a combination of the methodologies described in the studies of Xu et al., 2018 and Farris et al., 2014 [25,48], the potential cancer effect caused by exposure to air pollution was estimated. In particular, in Xu et al., 2018, the daily exposure dose (Dw) by air pollutants was calculated according to the following Equation (27):
D w = i = 1 C i × j = 1 T ij × IR j
C: the concentration of pollutant (in μg/m3), Tij is the exposure time (in h/h), IR is the inhalation rate adjusted with the bodyweight of people (m3/day-kg) [23]. Secondly, the potential cancer effects (Risk) from the study of Farris et al., 2014 [25] were based on the third edition of the book “Encyclopedia of Toxicology (2014)”, describing the cancer potency factor parameter (CPF) for quantifying the risk of chemical factors which are evaluated as carcinogenic (Equation (28)). ‘Dose’ is the exposure dose-adjusted with BW (mg/kg-day) and CPF is the cancer potency factor for each pollutant according to OEHHA, 2019 [43] (Table S3).
Risk = Dose × CPF

Methodology C.II

In this case, the potential cancer effect was estimated based on the combination of the methodologies of Maertens et al., 2004 and Farris et al., 2014, [13,24]. The difference from the previous methodology is the estimation of the Exposure Dose (or DW) of people to air pollutants, which was calculated through Equation (29), instead of Equation (27), Maertens, 2004. Subsequently, the probability of developing cancer (Risk) was calculated similarly to the previous methodology (C.I) [25,48].
Lifetime   average   daily   exposure   dose = C   ×   IR   ×   EF BW
Risk = Dose × CPF

3. Results and Discussion

The potential cancer and non cancer risk assessments for each methodology are presented in Table 5 and Table 6, and the estimated risks for each pollutant are presented in the Appendix A (Table A1, Table A2, Table A3, Table A4, Table A5, Table A6, Table A7, Table A8, Table A9, Table A10, Table A11, Table A12, Table A13 and Table A14).

3.1. Cancer Risk Assessment

The cancer risk (through inhalation) was estimated, using the methodologies described in the previous paragraphs, for each pollutant and then summed to estimate the total lifetime cancer risk, for the residents of Elefsina (Table 5). The results were different according to the methodology and the formulas applied in each case; the examined pollutants and the parameters are taken into account in each case (e.g., the age, the parameters of Exposure Duration, Cancer Potency).
Risk estimation using approaches A.I and A.II [13,14] resulted in low potential risks values (Table A1 and Table A2). The results for adults (men versus women) were similar differing to those for children (0–2years vs. 2–16years group) which were different due to the presence of the adjustment factor (AF). As a result, the calculated risk included the AF parameter, which was higher for every PAH for both age groups (Table 5 and Table 6). It was also observed that the risk of developing cancer increases as age decreases due to IR and BW factors. The higher value was calculated for benzo(a)pyrene (1.09 × 10−6) and dibenz[a,h]anthracene (1.95 × 10−6) (due to high PEF and C) and the lower one for chrysene (9.09 × 10−12).
On the contrary, the estimation performed using the methodology A.III, [31] showed conflicting results. The estimated cancer risk was higher for adults, compared to children as the length of the exposure time for adults (L) was larger, a parameter that is absent from approaches A.I and A.II. Furthermore, in the A.III method, the parameter of exposure frequency (EF) was included and the risk was not estimated only for PAHs. As a result, the total lifetime cancer risk (18.6 cases per 106 residents) and the individual risks for each pollutant were higher (Table A3). The higher potential risk was estimated for benzene (1.27 × 10−5) and the lower one for chrysene (6.88 × 10−12), similarly to AI and AII methods.
According to approach A.IV, the total lifetime cancer risk was estimated for two different UR cases (1.1 × 10−6, OEHHA [37] and 8.7 × 10−5, WHO [39]), while there were no differences for estimations related to age or gender of the residents (Table A4). Moreover, the parameters for daily exposure duration, exposure frequency, and the years of exposure (duration) are missing from this methodology. It is worthy to be mentioned that the RPFi parameter had a significant influence on the estimation of risk. Specifically, DBahA had the lowest concentration (0.02 ng/m3) among PAHs, but the corresponding estimated cancer risk was the second highest, because DBahA had a high value of the parameter RPFi compared to the other PAHs. Additionally, in the A.IV approach, the parameter (RPFi) for the cancer potency is different for each PAH [35]. In contrast, in the A.III method, it is the PEFi parameter that is used for the cancer risk estimation, leading to different results.
As a result, the estimated risks were higher with the A.IV methodology, with BkF (2.01 × 10−5) and DBahA (1.74 × 10−5) pollutants being the highest. Furthermore, in A.IV higher cancer assessment was estimated with the UR according to WHO (2010) for all PAHs, except for the case where UR was used according to OEHHA (2005).
According to the B.I approach, the total inhalation cancer risk as a result of an 80-year exposure to a variety of pollutants was estimated to be 33 chances per million people (Table 6). The risk increased with increasing age, while cancer risk for adults was higher in contrast to children (Table A5). The same assessment was conducted based on the methodologies B.II, B.III, and B.IV and the increase in the estimated risk was proportional to the increase in human age [15,27,28]. In the B.I methodology, the highest cancer risk of all pollutants was estimated for benzene (2.24 × 10−5, for Adults) and the lower one for chrysene (4.81 × 10−11, for Children0–2 years), (Table A5). With the application of methodology B.II, the lower total lifetime cancer risk was estimated (0.05 chances per million, Table A7). A significant difference between the two methodologies was the absence of age and gender parameters, resulting in no differentiation in the potential estimated risks between women, men, and children. Furthermore, the B.I method lacks the parameter of the daily exposure dose. The highest cancer risk was estimated for Ni (4.47 × 10−8) and the lowest for chrysene (9.78 × 10−12). In general, a parameter that significantly differentiates the outcomes of the several approaches is the Cancer Potency of each pollutant.
According to the B.III methodology, the Total Cancer Risk was estimated to be 11.7 cases per 106 humans and the risk assessment from long-term exposure are presented in Table A10. For men and women, the cancer risk was the same (9.3 cases per million) while for children0–16 was 2.3 cases per million (Table 5 and Table A8). The potential risk for men and women calculated for all pollutants gave the same results. It is important to mention that in the B.III methodology, the parameters IRd and BW were included in the equation for the Chronic daily intake dose Equation (16) and in the equation for the Cancer slope factor Equation (17). As a result, the fraction of BW and IR did not differentiate the estimated risk for women and men (Table A8). Similar results were observed for the assessment of cancer and non-cancer effects, which means that the parameters BW and IR were not properly expressed. For example, as shown in Table 5 and Table 6 or Table A8 and Table A9, the estimated values for men and women are the same, even though the parameter of gender is included. It is well known that the IR clearly constitutes a crucial parameter for the calculation of the daily intake dose, and consequently for the risk estimation.
Furthermore, the duration of the daily exposure was proven to be a critical parameter for potential cancer risk calculation (Table A8). Moreover, the number of new cancer cases (I) per lifetime (80 years) in Elefsina was found to be low (0.28 cases for adults and 0.07 for children), due to the limited number of residents (29,900). However, methodology B.IV presented significant differences (IR,BW) from B.III, the results of the cancer risks assessment were found to be the same as those of B.III methodology. In contrast to B.III, the parameters of IR and BW were absent from the B.IV methodology (Table A11). Nevertheless, the results among these two methodologies were similar, verifying that in the B.III methodology the parameters IR and BW do not lead to differentiation in the estimated risks, because they were used for the calculation of both CDI and CSF did not differentiate the assessment. In addition, in B.IV the parameter IUR was included, instead of the CSF parameter from the B.III methodology, in which the CSF Equation (17) was calculated from the IR, BW and IUR. For both methodologies B.III and B.IV the highest cancer risk was found for benzene (8.32 × 10−6 for adults) and the lowest for chrysene (1.24 × 10−11 for children0–16 years).
The combination of the methods (C.I & C.II) leads to the highest total cancer risk, i.e., 147.3 cases per million people (Table 5). The potential cancer risk was observed to increase as the age-range of people decreases, due to IR and BW factors. The daily exposure duration was proved to be an important parameter with 58% fewer cases (18.11 cases per million) of cancer estimated for men having 14 h of exposure per day (Table A14) compared to those men who were exposed for 24 h per day (31.04 cases per million). The comparison of the results from the cancer risk assessment as presented in the following Figure 3 and the difference of these estimated risks between benzene and other pollutants is noteworthy.
Specifically, if the potential cancer risk related to benzene is removed from the total cancer risk of all pollutants, then the cancer chances decrease significantly (children0–2 yrs: 68.44 cases per million with Benzene, 2.31 cases without Benzene, 0.25 without Benzene, Ni and Pb) and if these three pollutants are removed from the total assessment, the estimated cancer risks are also significantly lower. This result is independent of the range of age of the people (kids & adults), which could be also concluded from Table 5 or Table A13. It should be noted that the parameter of the length of the exposure duration (years) of humans was not included in the C.I methodology.

3.2. Non-Cancer Effects

The non-cancer risk (through inhalation) for the residents of Elefsina was estimated and compared, using the methodologies, for each pollutant (HQ) and then summed for the total lifetime non-cancer effects (HI), (Table 6 and Table 7). There were no exceedances οf the acceptable limits of HI and HQ, for all the methodologies used (HI & HQ < 1).
The results indicate several differences depending on the applied methodology (Table 8). Firstly, according to methodology B.I [24] the risk for non-cancer acute health impacts of each pollutant was estimated and expressed as a Hazard Quotient (individual substances) or a Hazard Index (multiple substances) for the target organ system(s). As shown in Table 6, the non-cancer risk was low, due to low concentrations of the pollutants. As a result, the HI values for every target organ system were low (Table A6), especially for the Respiratory, Nervous Systems, and Eyes (2.28 × 10−4), while the higher risk was calculated for Immune, Hematologic Systems, and Reproductive/Development (1.72 × 10−2). In this methodology, the non-cancer risk did not differentiate between age groups and gender, while the parameters of exposure duration were absent. On the other hand, by using method B.II, it was possible to estimate the total non-cancer risk for children and adults, based on their exposure to Ni. The risk was estimated to be 3.59 × 10−6 for children and 2.33 × 10−6 for adults (Table 6). According to the methodology B.III, the non-cancer effects were assessed for several pollutants, for different age ranges and gender of people, as well as for various exposure periods (Daily exposure, Exposure days/year & length of years exposure duration). The Index for the total non-cancer risk (HI) for the daily exposure time of 24 h per day was estimated to be 1.05, while the estimated HQ separately for children, women, and men was the same, 0.5269 (Table A9). As previously mentioned, IR and BW parameters did not affect the estimated risk value, as they are included in the equation for the calculation of reference dose (RFD) for each pollutant. In addition, for the case of children, the results coincide with those of adults, despite the fact that the duration of the exposure is practically different. This can be attributed to the parameter AT (ED years × 365 days/year), since for the calculation of the non-cancer effects, in AT parameter was included the ED (31), while in the calculation of cancer assessment the AT parameter was not included (80 years × 365 days/year).
Equation (16), Non-cancer: AT= ED years × 365 days/year, Cancer: AT= 80 years × 365 days/year.
The daily exposure duration proved to be a critical parameter for assessing the non-cancer impacts. The higher potential risk among the pollutants (HQ) was calculated for nickel (0.21) for children and the lower one for toluene (3.99 × 10−4) for adults. In the assessment for the non-cancer risk according to the methodology by B.IV [27] (Table A12), while there were significant differences in the parameters included in both methodologies, the results of the assessment were the same as the results by B.III [26]. In this methodology (B.IV), the parameters IR, BW, and RFD were not included in the equation. A situation that confirms the lack of influence of these parameters in B.III Equations (17) and (20), whether included or not. As a result, the highest and the lowest non-cancer risk for B.III & B.IV was estimated for Nickel (0.21) for children and for Toluene (3.99 × 10−4) for adults (women & men).

4. Conclusions

This study aimed to estimate the health risks caused by exposure to air pollution for people living in an urban, industrialized area of Greece. Health effects on different population subgroups were estimated for targeted pollutants (PM2.5, PAHs, heavy metals and VOCs) through different mathematical approaches provided by the literature and taking into consideration different parameters, such as age, gender and daily average exposure time. The main conclusions of the study are summarized in the following:
  • The estimation of the risk for potential cancer and non-cancer effects varies depending on the applied methodology; however, the acceptable limits for cancer risk provided by OEHHA, EPA and WHO were not exceeded in any case except for one approach, as well as only one case (C.I & C.II) which were found to be in excess of the acceptable limits of cancer risk (10−6–10−4) for the total lifetime.
  • Several methodologies lack critical parameters (e.g., IR, BW, AF) which are important for the estimation of the exposure dose of people, depending on their age, gender, and the daily activity of people. In six approaches in which there was the possibility to estimate the exposure specifically for women, the results were lower compared to those for men (with the exception of methods B.III & B.IV). The main reason was the existence of IR and BW parameters.
  • The duration of the exposure to air pollution (i.e., hours per day, days per year and number of years) is a major factor influencing the estimated risks and was differentiated as a function of the time the people were exposed to the pollutants, regardless of people gender, age or other parameters.
  • Age adjustment factor (AF) proves to be a critical parameter for cancer risk assessment as, during the early stages of life, people are more vulnerable to the effects of chemical toxins.
  • Benzene comprises a critical factor as its estimated cancer risk was significant, mainly due to its high concentration and its high cancer slope factor, as well. Additionally, pollutants with low concentrations (DBahA, BaP) have been associated with high potential risks, due to their high value of parameters for the corresponding risk (CSF,P or RFD, HQ).
  • In the majority of the methodologies used, the younger the age of people who were studied, the lower the estimated risk for cancer. In cases A.I, A.II, A.IV, C.I & C.II where the duration of exposure was not taken under consideration, the risk for children was equal to or higher than that of adults.

Limitations—Future Study

The assessment was performed according to literature references, as certain parameters used in risk calculation formulas were not available for all examined pollutants (i.e., cancer slope factor, reference dose). This could comprise a limitation of the present study but at the same time a challenge for future research. Moreover, it is important that there are no widely accepted values for all these parameters and as a result, different values of the parameters are selected in every study and consequently, lead to different assessments. Additionally, the collection of air pollution data of higher temporal (i.e., daily and seasonal variation) and spatial resolution (network of sampling sites) would contribute to a more robust risk estimation, as the influence of environmental and human activity factor variability would be taken into account. For this purpose, the scientific community needs to focus on more comprehensive mathematical models and approaches, which will take into consideration all the parameters of the relation between exposure to air pollution and population health.
The limited number of sampling sites, seasonal influences, changes in population vulnerabilities, and influences of lifestyle characteristics.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/atmos13040615/s1; Table S1. Values of the limit of detection (LOD), the limit of quantitation (LOQ) and the expanded uncertainty at 95% confidence level and k = 2 that has been estimated experimentally (Uexp) for each detected PAH; Table S2. Description of gaseous pollutant recording methods; Table S3. PAHs, PEF & SF (OEHHA, 2015 & 2019); Table S4. Studies were used for risk assessment of cancer and non cancer effects; Table S5. Cancer potency factor, OEHHA (2015 & 2019); Table S6. Concentrations, Acute RELs & Target Organ System(s) (OEHHA Appendices, 2015); Table S7. Parameters for non – cancer risks.

Author Contributions

Conceptualization, T.M. and D.S.; Methodology, T.M., M.P. and A.S.; Investigation, M.P. and P.P.; Analysis, M.D., D.B., M.I.M. and K.E.; Writing—Original Draft Preparation, M.P.; Writing—Review & Editing, D.S. and T.M. All authors have read and agreed to the published version of the manuscript.

Funding

Not applicable.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the principals of the schools for providing the sampling sites.

Conflicts of Interest

The authors declare sole responsibility for the research results.

Glossary

PMParticulate Matter
IARCInternational Agency for Research on Cancer
PAHsPolycyclic Aromatic Hydrocarbons
DALYsDisability Adjusted Life Years
VOCsVolatile Organic Compounds
SVOCsSemivolatile Organic Compounds
IRInhalation Rate
BWBody weight
CRCancer Risk
OEHHAOffice of Environmental Health Hazard Assessment
PEFPotency Equivalency Factor
SFSlope Factor
EFExposure Factor
AFAdjustment Factor
EDExposure Duration
ATAveraging Time
RPFRelative Potency Factor
ASFAge Sensitivity Factor
FAHFraction At Home
RfDReference Dose
RRRelative Risk
CDIChronic Daily Intake
IURInhalation Unit Risk
ECExposure Concentration

Appendix A. Tables with Results of Assessment

Table A1. Cancer risk assessment A.I method (Maertens, 2008).
Table A1. Cancer risk assessment A.I method (Maertens, 2008).
Lifetime Cancer Risk:
PAHsC (ng/m3)Children 0–2 YearsChildren 2–16 YearsWomen 16–80 YearsMen 16–80 Years
14 h/Day24 h/Day14 h/Day24 h/Day
BaP0.0591.09 × 10−62.29 × 10−72.68 × 10−84.60 × 10−82.90 × 10−84.97 × 10−8
BaA0.0539.83 × 10−92.06 × 10−92.41 × 10−104.13 × 10−102.60 × 10−104.46 × 10−10
BbF0.2895.36 × 10−81.12 × 10−81.31 × 10−92.25 × 10−91.42 × 10−92.43 × 10−9
BkF0.0791.47 × 10−83.07 × 10−93.59 × 10−106.16 × 10−103.88 × 10−106.65 × 10−10
CHRY0.1082.00 × 10−104.20 × 10−114.91 × 10−128.42 × 10−125.30 × 10−129.09 × 10−12
DBahA0.021.95 × 10−64.09 × 10−74.78 × 10−88.20 × 10−85.16 × 10−88.85 × 10−8
I123cdP0.1352.50 × 10−85.25 × 10−96.14 × 10−101.05 × 10−96.63 × 10−101.14 × 10−9
Σ3.15 × 10−66.60 × 10−77.72 × 10−81.32 × 10−78.33 × 10−81.43 × 10−7
Cases/1063.150.660.080.130.080.14
Table A2. Cancer risk assessment A.II method (Maertens, 2004).
Table A2. Cancer risk assessment A.II method (Maertens, 2004).
Lifetime Cancer Risk
PAHsC (ng/m3)Children 0–2 YearsChildren 2–16 YearsWomen 16–80 YearsMen 16–80 Years
14 h/Day24 h/Day14 h/Day24 h/Day
BaP0.0591.09 × 10−77.65 × 10−82.68 × 10−84.60 × 10−82.90 × 10−84.97 × 10−8
BaA0.0539.83 × 10−106.87 × 10−102.41 × 10−104.13 × 10−102.60 × 10−104.46 × 10−10
BbF0.2895.36 × 10−93.75 × 10−91.31 × 10−92.25 × 10−91.42 × 10−92.43 × 10−9
BkF0.0791.47 × 10−91.02 × 10−93.59 × 10−106.16 × 10−103.88 × 10−106.65 × 10−10
CHRY0.1082.00 × 10−111.40 × 10−114.91 × 10−128.42 × 10−125.30 × 10−129.09 × 10−12
DBahA0.021.95 × 10−71.36 × 10−74.78 × 10−88.20 × 10−85.16 × 10−88.85 × 10−8
I123cdP0.1352.50 × 10−91.75 × 10−96.14 × 10−101.05 × 10−96.63 × 10−101.14 × 10−9
Σ3.15 × 10−72.20 × 10−77.72 × 10−81.32 × 10−78.33 × 10−81.43 × 10−7
Cases/1060.310.220.080.130.080.14
Table A3. Cancer risk assessment A.III method (Jackson, 2005).
Table A3. Cancer risk assessment A.III method (Jackson, 2005).
Lifetime Cancer Risk
PollutantsC (ng/m3)Children 0–16 YearsWomen 16–80 YearsMen 16–80 Years
14 h/Day24 h/Day14 h/Day24 h/Day
BaP0.0593.76 × 10−95.28 × 10−99.05 × 10−95.70 × 10−99.77 × 10−9
BaA0.0533.38 × 10−104.74 × 10−108.13 × 10−105.12 × 10−108.77 × 10−10
BbF0.2891.84 × 10−92.59 × 10−94.43 × 10−92.79 × 10−94.78 × 10−9
BkF0.0795.03 × 10−117.07 × 10−111.21 × 10−107.63 × 10−111.31 × 10−10
CHRY0.1086.88 × 10−129.67 × 10−121.66 × 10−111.04 × 10−111.79 × 10−11
DBahA0.021.27 × 10−91.79 × 10−93.07 × 10−91.93 × 10−93.31 × 10−9
I123cdP0.1358.60 × 10−101.21 × 10−92.07 × 10−91.30 × 10−92.23 × 10−9
1Methyl0.5771.07 × 10−91.50 × 10−92.57 × 10−91.62 × 10−92.77 × 10−9
Ni4.382.54 × 10−73.57 × 10−76.12 × 10−73.85 × 10−76.60 × 10−7
Pb8.032.15 × 10−83.02 × 10−85.17 × 10−83.26 × 10−85.58 × 10−8
Benzene13904.87 × 10−66.84 × 10−61.17 × 10−57.38 × 10−61.27 × 10−5
Σ5.83 × 10−61.10 × 10−51.89 × 10−51.15 × 10−51.98 × 10−5
Cases/1065.27.212.47.813.4
Table A4. Cancer risk assessment A.IV method (Gao, 2019).
Table A4. Cancer risk assessment A.IV method (Gao, 2019).
PollutantsC (ng/m3)Lifetime Cancer Risk
UR (OEHHA. 2005)UR (WHO. 2010)
BaP0.0596.49 × 10−85.13 × 10−6
BaA0.0531.17 × 10−89.22 × 10−7
BbF0.2892.54 × 10−72.01 × 10−5
BkF0.0792.61 × 10−92.06 × 10−7
CHRY0.1081.19 × 10−89.40 × 10−7
DBahA0.022.20 × 10−71.74 × 10−5
I123cdP0.1351.04 × 10−88.22 × 10−7
BghiP0.2222.20 × 10−91.74 × 10−7
FLA0.1361.20 × 10−89.47 × 10−7
Σ5.90 × 10−74.67 × 10−5
Cases/1060.5946.66
Table A5. Cancer risk assessment B.I method (OEHHA, 2015).
Table A5. Cancer risk assessment B.I method (OEHHA, 2015).
Lifetime Cancer Risk
PollutantsC (ng/m3)Children 0–2 YearsChildren 2–16 YearsAdults 16–80 YearsTotal Cancer Risk 80 Years
BaP0.0592.63 × 10−91.28 × 10−83.71 × 10−85.26 × 10−8
BaA0.0532.36 × 10−101.15 × 10−93.34 × 10−94.73 × 10−9
BbF0.2891.29 × 10−96.29 × 10−91.82 × 10−82.58 × 10−8
BkF0.0793.52 × 10−101.72 × 10−94.97 × 10−97.04 × 10−9
CHRY0.1084.81 × 10−112.35 × 10−106.80 × 10−109.63 × 10−10
DBahA0.029.36 × 10−104.58 × 10−91.32 × 10−81.87 × 10−8
I123cdP0.1356.01 × 10−102.94 × 10−98.50 × 10−91.20 × 10−8
Nickel4.384.55 × 10−82.23 × 10−76.43 × 10−79.11 × 10−7
Lead8.033.85 × 10−91.88 × 10−85.44 × 10−87.71 × 10−8
Benzene13901.59 × 10−67.76 × 10−62.24 × 10−53.18 × 10−5
Σ1.64 × 10−68.03 × 10−62.32 × 10−53.29 × 10−5
Cases/1061.648.023.233
Table A6. HQ for each Target organ System(s) & Acute HI. B.I method (OEHHA, 2015).
Table A6. HQ for each Target organ System(s) & Acute HI. B.I method (OEHHA, 2015).
SubstanceC (μg/m3)Total HQAcuteImmune SystemReproductive/ DevelopmentHematologic SystemRespiratory SystemEyesNervous System
Benzene1.395.1 × 10−21.72 × 10−21.72 × 10−21.72 × 10−2
Toluene3.569.62 × 10−52.41 × 10−52.41 × 10−52.41 × 10−52.41 × 10−5
Xylenes 13.081.40 × 10−44.67 × 10−54.67 × 10−54.67 × 10−5
Xylenes 22.761.25 × 10−44.18 × 10−54.18 × 10−54.18 × 10−5
Xylenes 37.533.42 × 10−41.14 × 10−41.14 × 10−41.14 × 10−4
Total Acute Hazard Index 1.72 × 10−21.72 × 10−21.72 × 10−22.27 × 10−42.27 × 10−42.28 × 10−4
Table A7. Risk assessment B.II method (Tianjie, 2018).
Table A7. Risk assessment B.II method (Tianjie, 2018).
PollutantsC (mg/kg)Lifetime Cancer Risk Non Carcinogenic Risk
ADD0–16 YearsADD16–80 YearsHQ0–16 YearsHQ16–80 Years
Ni315.74.47 × 10−8 7.40 × 10−84.80 × 10−83.59 × 10−62.33 × 10−6
BaP4.464.14 × 10−11NO RfDNO RfD
BaA3.914.81 × 10−10
BbF21.592.66 × 10−9
BkF5.937.30 × 10−11
CHRY7.959.78 × 10−12
DBahA1.491.83 × 10−9
I123cdP10.201.26 × 10−9
Σ5.11 × 10−8
Cases/1060.05
Table A8. Cancer risk assessment B.III method (Chalvatzaki, 2019).
Table A8. Cancer risk assessment B.III method (Chalvatzaki, 2019).
Lifetime Cancer Risk
PollutantsC (ng/m3)Children 0–16 YearsWomen 16–80 YearsMen 16–80 Years
14 h/Day24 h/Day14 h/Day24 h/Day
PAHBaP0.0596.79 × 10−91.58 × 10−82.72 × 10−81.58 × 10−82.72 × 10−8
BaA0.0536.10 × 10−101.42 × 10−92.44 × 10−91.42 × 10−92.44 × 10−9
BbF0.2893.33 × 10−97.76 × 10−91.33 × 10−87.76 × 10−91.33 × 10−8
BkF0.0799.09 × 10−112.12 × 10−103.64 × 10−102.12 × 10−103.64 × 10−10
CHRY0.1081.24 × 10−112.90 × 10−114.97 × 10−112.90 × 10−114.97 × 10−11
DBahA0.022.30 × 10−95.37 × 10−99.21 × 10−95.37 × 10−99.21 × 10−9
I123cdP0.1351.55 × 10−93.62 × 10−96.21 × 10−93.62 × 10−96.21 × 10−9
Heavy metalsNi4.382.18 × 10−75.10 × 10−78.74 × 10−75.10 × 10−78.74 × 10−7
Pb8.031.85 × 10−84.31 × 10−87.39 × 10−84.31 × 10−87.39 × 10−8
VOCBenzene13902.08 × 10−64.85 × 10−68.32 × 10−64.85 × 10−68.32 × 10−6
Σ2.33 × 10−65.44 × 10−69.32 × 10−65.44 × 10−69.32 × 10−6
Cases/1062.35.49.35.49.3
I (Elefsina) *0.070.160.280.160.28
* Cases cancer risk for Elefsina (29900 residents).
Table A9. Non-cancer risk assessment B.III method (Chalvatzaki, 2019).
Table A9. Non-cancer risk assessment B.III method (Chalvatzaki, 2019).
Hazard Quotient
PollutantsC (ng/m3)Children 0–16 YearsWomen 16–80 YearsMen 16–80 Years
14 h/Day24 h/Day14 h/Day24 h/Day
BaP0.0590.02830.01650.02830.01650.0283
Ni4.380.21000.12250.21000.12250.2100
Pb8.030.03850.02250.03850.02250.0385
Benzene13900.04440.02590.04440.02590.0444
Toluene35636.83 × 10−43.99 × 10−46.83 × 10−43.99 × 10−46.83 × 10−4
Trimethylben23003.68 × 10−22.14 × 10−23.68 × 10−22.14 × 10−23.68 × 10−2
m-Xylene30772.95 × 10−21.72 × 10−22.95 × 10−21.72 × 10−22.95 × 10−2
o-Xylene27602.65 × 10−21.54 × 10−22.65 × 10−21.54 × 10−22.65 × 10−2
p-Xylene75337.22 × 10−24.21 × 10−27.22 × 10−24.21 × 10−27.22 × 10−2
Cyclohexane71871.15 × 10−36.70 × 10−31.15 × 10−36.70 × 10−31.15 × 10−3
Nonane7733.71 × 10−22.16 × 10−23.71 × 10−22.16 × 10−23.71 × 10−2
Hazard Index0.52690.30740.52690.30740.5269
Table A10. Risk assessment from long-term exposure to PM, B.III method (Chalvatzaki, 2019).
Table A10. Risk assessment from long-term exposure to PM, B.III method (Chalvatzaki, 2019).
ParametersAll Cause Mortality (PM10)
(β = 0.0008)
Cardiopulmonary Mortality (PM2.5)
(β = 0.15515)
Lung Cancer Mortality (PM2.5)
(β = 0.23218)
ER (95% CI)0.0080.2230.352
AF (95% CI)0.0080.1830.260
Deaths (95% CI)5.76 × 10−51.36 × 10−31.93 × 10−3
Table A11. Cancer risk assessment B.IV method (Megido, 2017).
Table A11. Cancer risk assessment B.IV method (Megido, 2017).
Lifetime Cancer Risk:
PollutantsC (ng/m3)Children 0–16 YearsAdults 16–80 Years
14 h/Day24 h/Day
PAHBaP0.0596.79 × 10−91.58 × 10−82.72 × 10−8
BaA0.0536.10 × 10−101.42 × 10−92.44 × 10−9
BbF0.2893.33 × 10−97.76 × 10−91.33 × 10−8
BkF0.0799.09 × 10−112.12 × 10−103.64 × 10−10
CHRY0.1081.24 × 10−112.90 × 10−114.97 × 10−11
DBahA0.022.30 × 10−95.37 × 10−99.21 × 10−9
I123cdP0.1351.55 × 10−93.62 × 10−96.21 × 10−9
Heavy metalsNi4.382.18 × 10−75.10 × 10−78.74 × 10−7
Pb8.031.85 × 10−84.31 × 10−87.39 × 10−8
VOCBenzene13902.08 × 10−64.85 × 10−68.32 × 10−6
Σ2.33 × 10−65.44 × 10−69.32 × 10−6
Cases/1062.35.49.3
Table A12. Non-cancer risk assessment B.IV method (Megido, 2017).
Table A12. Non-cancer risk assessment B.IV method (Megido, 2017).
Hazard Quotient:
PollutantsC (ng/m3)Children 0–16 YearsAdults 16–80 Years
14h/Day24h/Day
PAHBaP0.0590.02830.01650.0283
Heavy metalsNi4.380.21000.12250.2100
Pb8.030.03850.02250.0385
VOCBenzene13900.04440.02590.0444
Toluene35636.83 × 10−43.99 × 10−46.83 × 10−4
Trimethylben23003.68 × 10−22.14 × 10−23.68 × 10−2
m-Xylene30772.95 × 10−21.72 × 10−22.95 × 10−2
o-Xylene27602.65 × 10−21.54 × 10−22.65 × 10−2
p-Xylene75337.22 × 10−24.21 × 10−27.22 × 10−2
Cyclohexane71871.15 × 10−36.70 × 10−31.15 × 10−3
Nonane7733.71 × 10−22.16 × 10−23.71 × 10−2
Hazard Index0.52690.30740.5269
Table A13. Cancer risk assessment C.I method (Xu, 2018 / Farris, 2014).
Table A13. Cancer risk assessment C.I method (Xu, 2018 / Farris, 2014).
Lifetime Cancer Risk
C (ng/m3)Children 0–2 YearsChildren 2–16 YearsWomen 16–80 YearsMen 16–80 Years
14 h/Day24 h/Day14 h/Day24 h/Day
BaP0.0591.09 × 10−77.65 × 10−82.68 × 10−84.60 × 10−82.90 × 10−84.97 × 10−8
BaA0.0539.83 × 10−96.87 × 10−92.41 × 10−94.13 × 10−92.60 × 10−94.46 × 10−9
BbF0.2895.36 × 10−83.75 × 10−81.31 × 10−82.25 × 10−81.42 × 10−82.43 × 10−8
BkF0.0791.47 × 10−81.02 × 10−83.59 × 10−96.16 × 10−93.88 × 10−96.65 × 10−9
CHRY0.1082.00 × 10−91.40 × 10−94.91 × 10−108.42 × 10−105.30 × 10−109.09 × 10−10
DBahA0.023.90 × 10−82.72 × 10−89.57 × 10−91.64 × 10−81.03 × 10−81.77 × 10−8
I123cdP0.1352.50 × 10−81.75 × 10−86.14 × 10−91.05 × 10−86.63 × 10−91.14 × 10−8
Ni4.381.90 × 10−61.32 × 10−64.65 × 10−77.97 × 10−75.02 × 10−78.60 × 10−7
Pb8.031.60 × 10−71.12 × 10−73.93 × 10−86.75 × 10−84.25 × 10−87.28 × 10−8
Benzene13906.61 × 10−54.62 × 10−51.62 × 10−52.78 × 10−51.75 × 10−53.00 × 10−5
Σ6.84 × 10−54.78 × 10−51.68 × 10−52.88 × 10−51.81 × 10−53.10 × 10−5
Cases/10668.4447.8016.7828.7718.1131.04
Without Benzene2.311.610.570.970.611.05
Without Benz, Ni, Pb0.250.180.060.110.070.12
Table A14. Cancer risk assessment C.II method (Maertens, 2004/ Farris, 2014).
Table A14. Cancer risk assessment C.II method (Maertens, 2004/ Farris, 2014).
Lifetime Cancer Risk:
C (ng/m3)Children 0–2 YearsChildren 2–16 YearsWomen 16–80 YearsMen 16–80 Years
14 h/Day24 h/Day14 h/Day24 h/Day
BaP0.0591.09 × 10−77.65 × 10−82.68 × 10−84.60 × 10−82.90 × 10−84.97 × 10−8
BaA0.0539.83 × 10−96.87 × 10−92.41 × 10−94.13 × 10−92.60 × 10−94.46 × 10−9
BbF0.2895.36 × 10−83.75 × 10−81.31 × 10−82.25 × 10−81.42 × 10−82.43 × 10−8
BkF0.0791.47 × 10−81.02 × 10−83.59 × 10−96.16 × 10−93.88 × 10−96.65 × 10−9
CHRY0.1082.00 × 10−91.40 × 10−94.91 × 10−108.42 × 10−105.30 × 10−109.09 × 10−10
DBahA0.023.90 × 10−82.72 × 10−89.57 × 10−91.64 × 10−81.03 × 10−81.77 × 10−8
I123cdP0.1352.50 × 10−81.75 × 10−86.14 × 10−91.05 × 10−86.63 × 10−91.14 × 10−8
Ni4.381.90 × 10−61.32 × 10−64.65 × 10−77.97 × 10−75.02 × 10−78.60 × 10−7
Pb8.031.60 × 10−71.12 × 10−73.93 × 10−86.75 × 10−84.25 × 10−87.28 × 10−8
Benzene13906.61 × 10−54.62 × 10−51.62 × 10−52.78 × 10−51.75 × 10−53.00 × 10−5
Σ6.84 × 10−54.78 × 10−51.68 × 10−52.88 × 10−51.81 × 10−53.10 × 10−5
Cases/10668.4447.8016.7828.7718.1131.04
Without Benzene2.311.610.570.970.611.05
Without Benz, Ni, Pb0.250.180.060.110.070.12

References

  1. Olstrup, H.; Johansson, C.; Forsberg, B.; Åström, C.; Orru, H. Seasonal Variations in the Daily Mortality Associated with Exposure to Particles, Nitrogen Dioxide, and Ozone in Stockholm, Sweden, from 2000 to 2016. Atmosphere 2021, 12, 1481. [Google Scholar] [CrossRef]
  2. Emmanouil, C.; Drositi, E.; Vasilatou, V.; Diapouli, E.; Krikonis, K.; Eleftheriadis, K.; Kungolos, A. Study on particulate matter air pollution, source origin, and human health risk based of PM10 metal content in Volos City, Greece. Toxicol. Environ. Chem. 2016, 99, 691–709. [Google Scholar] [CrossRef]
  3. Scorer, R.S. The meteorological scene. In Air Pollution Meteorology; Woodhead Publishing: Sawston, UK, 2002; pp. 1–20. ISBN 9781898563938. [Google Scholar] [CrossRef]
  4. WHO (World Health Organization). WHO Ambient Air Pollution: A Global Assessment of Exposure and Burden of Disease; WHO: Geneva, Switzerland, 2016; ISBN 9789241511353. Available online: https://www.who.int/phe/publications/air-pollution-global-assessment/en/ (accessed on 14 February 2022).
  5. Gao, P.; Guo, H.; Wang, S.; Guo, L.; Xinge, Y.; Yaoe, C.; Jiae, L.; Fana, Q.; Hang, J. In Vitro investigations of high molecular weight polycyclic aromatic hydrocarbons in winter airborne particles using simulated lung fluids. Atmos. Environ. 2019, 201, 293–300. [Google Scholar] [CrossRef]
  6. WHO (World Health Organization). 7 Million Premature Deaths Annually Linked to Air Pollution; WHO Press Release: Geneva, Switzerland, 2014; Available online: https://www.who.int/mediacentre/news/releases/2014/air-pollution/en/#:~:text=25%20March%202014%20%7C%20Geneva%20%2D%20In,result%20of%20air%20pollution%20exposure (accessed on 14 February 2022).
  7. World Health Organization. WHO Global Air Quality Guidelines: Particulate Matter (PM2.5 and PM10), Ozone, Nitrogen Dioxide, Sulfur Dioxide and Carbon Monoxide; World Health Organization: Geneva, Switzerland, 2021; Available online: https://apps.who.int/iris/handle/10665/345329. (accessed on 1 March 2022).
  8. IARC. Outdoor Air Pollution. IARC Monographs on the Evaluation of Carcinogenic Risks to Humans. International Agency for Research on Cancer, France; World Health Organization: Geneva, Switzerland, 2016; Volume 109, Available online: http://monographs.iarc.fr/ENG/Monographs/vol109/mono109.pdf (accessed on 14 February 2022).
  9. Sun, Y.; Hu, X.; Wu, J.; Lian, H.; Chen, Y. Fractionation and health risks of atmospheric particle-bound As and heavy metals in summer and winter. Sci. Total Environ. 2014, 493, 487–494. [Google Scholar] [CrossRef] [PubMed]
  10. Lim, S.S.; Vos, T.; Flaxman, A.D.; Danaei, G.; Shibuya, K.; Adair-Rohani, H.; AlMazroa, M.A.; Amann, M.; Anderson, H.R.; Andrews, K.G.; et al. A comparative risk assessment of burden of disease and injury attributable to 67 risk factors and risk factor clusters in 21 regions, 1990–2010: A systematic analysis for the Global Burden of Disease Study 2010. Lancet 2012, 380, 2224–2260. [Google Scholar] [CrossRef] [Green Version]
  11. Hänninen, O.; Knol, A. European Perspectives on Environmental Burden of Disease Estimates for Nine Stressors in Six European Countries (2011); National Institute for Health and Welfare (THL): Helsinki, Finland, 2011; ISBN 978-952-245-413-3. [Google Scholar]
  12. Collins, J.F.; Brown, J.P.; Alexeeff, G.V.; Salmon, A.G. Potency equivalency factors for some polycyclic aromatic hydrocarbons and polycyclic aromatic hydrocarbon derivatives. Regul. Toxicol. Pharmacol. 1998, 28, 45–54. [Google Scholar] [CrossRef] [PubMed]
  13. Maertens, R.M.; Bailey, J.; White, P.A. The mutagenic hazards of settled house dust: A review. Mutat. Res. 2004, 567, 401–425. [Google Scholar] [CrossRef]
  14. Maertens, R.M.; Yang, X.; Zhu, J.; Gagne, R.W.; Douglas, G.R.; White, P.A. Mutagenic and carcinogenic hazards of settled house dust I: Polycyclic aromatic hydrocarbon content and excess lifetime cancer risk from preschool exposure. Environ. Sci. Technol. 2008, 42, 1747–1753. [Google Scholar] [CrossRef] [Green Version]
  15. Shao, T.; Pan, L.; Chen, Z.; Wang, R.; Li, W.; Qin, Q.; He, Y. Content of Heavy Metal in the Dust of Leisure Squares and Its Health Risk Assessment—A Case Study of Yanta District in Xi’an. Int. J. Environ. Res. Public Health 2018, 15, 394. [Google Scholar] [CrossRef] [Green Version]
  16. Stamatelopoulou, A.; Dasopoulou, M.; Bairachtari, K.; Karavoltsos, S.; Sakellari, A.; Maggos, T. Contamination and Potential Risk Assessment of Polycyclic Aromatic Hydrocarbons (PAHs) and Heavy Metals in House Settled Dust Collected from Residences of Young Children. Appl. Sci. 2021, 11, 1479. [Google Scholar] [CrossRef]
  17. Hong, W.J.; Jia, H.; Yang, M.; Li, Y.F. Distribution, seasonal trends, and lung cancer risk of atmospheric polycyclic aromatic hydrocarbons in North China: A three-year case study in Dalian city. Ecotoxicol. Environ. Saf. 2020, 196, 110526. [Google Scholar] [CrossRef] [PubMed]
  18. Gao, P.; Guo, H.; Zhang, Z.; Ou, C.; Hang, J.; Fan, Q.; He, C.; Wu, B.; Feng, Y.; Xing, B. Bioaccessibility and exposure assessment of trace metals from urban airborne particulate matter (PM10 and PM2.5) in simulated digestive fluid. Environ. Pollut. 2018, 242, 1669–1677. [Google Scholar] [CrossRef] [PubMed]
  19. Ostro, B.; Prüss-üstün, A.; Campbell-lendrum, D.; Corvalán, C.; Woodward, A. Outdoor Air Pollution: Assessing the Environmental Burden of Disease at National and Local Levels; Environmental Burden of Disease Series, No. 5; World Health Organization Protection Human Environmental: Geneva, Switzerland, 2004; Available online: https://www.who.int/quantifying_ehimpacts/publications/ebd5/en/ (accessed on 14 February 2022).
  20. Burnett, R.T.; Pope, A., III; Ezzati, M.; Olives, C.; Lim, S.S.; Mehta, S.; Shin, H.H.; Singh, G.; Hubbell, B.; Brauer, M.; et al. An Integrated Risk Function for Estimating the Global Burden of Disease Attributable to Ambient Fine Particulate Matter Exposure. Environ. Health Perspect. 2014, 122, 4. [Google Scholar] [CrossRef] [PubMed]
  21. Manousakas, M.; Diapouli, E.; Papaefthymiou, H.; Kantarelou, V.; Zarkadas, C.; Kalogridis, A.-C.; Karydas, A.G.; Eleftheriadis, K. XRF characterization and source apportionment of PM10 samples collected in a coastal city. X-ray Spectrom. 2018, 47, 190–200. [Google Scholar] [CrossRef]
  22. ICRP. International Commission on Radiological Protection. In Basic Anatomical and Physiological Data for Use in Radiological Protection: Reference Values; Valentin, J., Ed.; ICRP: Stockholm, Sweden, 2002; Volume 32, ISBN 008 0442668. Available online: http://www.icrp.org/publication.asp?id=icrp%20publication%2089 (accessed on 14 February 2022).
  23. U.S. Environmental Protection Agency. Exposure Factors Handbook; National Center for Environmental Assessment: Washington, DC, USA, 2011.
  24. OEHHA. Office of Environmental Health Hazard Assessment, Air Toxics Hot Spots Program. In Guidance Manual for Preparation of Health Risk Assessment; California Environmental Protection Agency U.S.: Sacramento, CA, USA, 2015; 231. Available online: http://oehha.ca.gov/air/hot_spots/2015/2015GuidanceManual.pdf (accessed on 14 February 2022).
  25. Farris, F.F.; Ray, S.D. Cancer Potency Factor. In Encyclopedia of Toxicology, 3rd ed.; Elsevier: Amsterdam, The Netherlands, 2014; Volume 1, pp. 642–644. [Google Scholar] [CrossRef]
  26. Chalvatzaki, E.; Chatoutsidou, S.E.; Lehtomäki, H.; Almeida, S.M.; Eleftheriadis, K.; Hänninen, O.; Lazaridis, M. Characterization of human health risks from particulate air pollution in selected European cities. Atmosphere 2019, 10, 1–16. [Google Scholar] [CrossRef] [Green Version]
  27. Megido, L.; Suárez-Peña, B.; Negral, L.; Castrillón, L.; Fernández-Nava, Y. Suburban air quality: Human health hazard assessment of potentially toxic elements in PM10. Chemosphere 2017, 177, 284–291. [Google Scholar] [CrossRef]
  28. Huang, M.; Wang, W.; Chan, C.Y.; Cheung, K.C.; Man, Y.B.; Wang, X.; Wong, M.H. Contamination and risk assessment (based on bioaccessibility via ingestion and inhalation) of metal(loid)s in outdoor and indoor particles from urban centers of Guangzhou, China. Sci. Total Environ. 2014, 479–480, 117–124. [Google Scholar] [CrossRef]
  29. U.S. EPA; Cogliano, J.; Flowers, L.; Valcovic, L.; Barton, H.; Tracey Woodruff, T.; Choksi, N. Supplemental Guidance for Assessing Cancer Susceptibility from Early-Life Exposure to Carcinogens; EPA/630/R-03/003; U.S. Environmental Protection Agency: Washington, DC, USA, 2003; p. 20460. Available online: www.epa.gov/ncea/raf/cancer2003.htm (accessed on 14 February 2022).
  30. U.S. EPA; Barton, H.; Cogliano, J.; Firestone, M.P.; Flowers, L.; Valcovic, L.; Setzer, R.W.; Woodruff, T. Supplemental Guidance for Assessing Susceptibility from Early-Life Exposure to Carcinogens; EPA/630/R-03/003F; U.S. Environmental Protection Agency: Washington, DC, USA, 2005; p. 20460.
  31. Jackson, M.M. Roadside concentration of gaseous and particulate matter pollutants and risk assessment in dar-Es-Salaam, Tanzania. Environ. Monit. Assess. 2005, 104, 385–407. [Google Scholar] [CrossRef]
  32. U.S. EPA. Guidelines for Carcinogen Risk Assessment; EPA/630/R-00/004; Risk Assessment Forum U.S. Environmental Protection Agency: Washington, DC, USA, 1986. Available online: https://cfpub.epa.gov/ncea/risk/recordisplay.cfm?deid=54933 (accessed on 14 February 2022).
  33. U.S. EPA. Regional Screening Level (RSL) Summary Table. 2019. Available online: https://www.epa.gov/risk/regional-screening-levels-rsls-generic-tables (accessed on 14 February 2022).
  34. OEHHA Appendices. Office of Environmental Health Hazard Assessment, Air Toxics Hot Spots Program; Guidance Manual for Preparation of Health Risk Assessment, Appendix A-I, 245. 2015. Available online: https://oehha.ca.gov/air/crnr/notice-adoption-air-toxics-hot-spots-program-guidance-manual-preparation-health-risk-0 (accessed on 14 February 2022).
  35. U.S. EPA. Development of a Relative Potency Factor (RPF) Approach for Polycyclic Aromatic Hydrocarbon (PAH) Mixtures; External Review Draft; Integrated Risk Information System; U.S. EPA: Washington, DC, USA, 2010.
  36. OEHHA. Benzo[a]pyrene as a Toxic Air Contaminant. In Part B. Health Effects of Benzo[a]pyrene; California Environmental Protection Agency, Office of Environmental Health Hazard Assessment, Air Toxicology and Epidemiology Section: Berkeley, CA, USA, 1993. [Google Scholar]
  37. OEHHA. Air Toxics Hot Spots Program Risk Assessment Guidelines; California Environmental Protection Agency, Office of Environmental Health Hazard Assessment, Air Toxicology and Epidemiology Section: Oakland, CA, USA, 2005.
  38. World Health Organization. Air Quality Guidelines Copenhagen, Regional Office for Europe Second edition. Air Qual. Guidel. 2000, 22, 1–8. [Google Scholar] [CrossRef]
  39. World Health Organization. WHO Guidelines for Indoor Air Quality: Selected Pollutants; Regional Office for Europe: København, Denmark, 2010; Available online: https://apps.who.int/iris/handle/10665/260127 (accessed on 14 February 2022).
  40. U.S. EPA. Risk Assessment Guidance for Superfund: Volume I—Human Health Evaluation Manual (Part E, Supplemental Guidance for Dermal Risk Assessment) (Final); EPA/540/R/99/005; Office of Superfund Remediation and Technology Innovation, U.S. Environmental Protection Agency: Washington, DC, USA, 2004. Available online: https://www.epa.gov/sites/production/files/2015-09/documents/part_e_final_revision_10-03-07.pdf (accessed on 14 February 2022).
  41. U.S. EPA. Risk Assessment Guidance for Superfund Volume I: Human Health Evaluation Manual (Part F, Supplemental Guidance for Inhalation Risk Assessment); Office of Superfund Remediation and Technology Innovation Environmental Protection Agency: Washington, DC, USA, 2009; Volume I, pp. 1–68. Available online: https://www.epa.gov/sites/production/files/2015-09/documents/partf_200901_final.pdf (accessed on 14 February 2022).
  42. U.S. EPA. Risk Assessment Guidance for Superfund: Volume I—Human Health Evaluation Manual (Part A) (Interim Final); EPA/540/1e89/002; Office of Emergency and Remedial Response, U.S. Environmental Protection Agency: Washington, DC, USA, 1989.
  43. OEHHA Appendices. Office of Environmental Health Hazard Assessment Appendix A: Hot Spots Unit Risk and Cancer Potency Values. 2019. Available online: https://oehha.ca.gov/air/crnr/technical-support-document-cancer-potency-factors-2009 (accessed on 14 February 2022).
  44. U.S. EPA. Supplemental Guidance for Developing Soil Screening Levels for Superfund Sites; Office of Emergency and Remedial Response, U.S. Environmental Protection Agency: Washington, DC, USA, 2002; Volume 106, p. 20460.
  45. Ferreira-Baptista, L.; De Miguel, E. Geochemistry and risk assessment of street dust in Luanda, Angola: A tropical urban environment. Atmos. Environ. 2005, 39, 4501–4512. [Google Scholar] [CrossRef] [Green Version]
  46. U.S. EPA. Regional Screening Level (RSL) Summary Table; 2017. Available online: https://www.epa.gov/risk/regional-screening-levels-rsls-generic-tables (accessed on 1 March 2022).
  47. Cárdaba, A.M.; Moreno, M.M.; Medina, A.A.; Capitán, M.A.; Vaquer, F.C.; Gómez, A.A. Health impact assessment of air pollution in Valladolid, Spain. BMJ Open 2014, 4, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Xu, T.; Hou, J.; Cheng, J.; Zhang, R.; Yin, W.; Huang, C.; Zhu, X.; Chen, W.; Yuan, J. Estimated individual inhaled dose of fine particles and indicators of lung function: A pilot study among Chinese young adults. Environ. Pollut. 2018, 235, 505–513. [Google Scholar] [CrossRef] [PubMed]
Figure 1. The sampling site and the industrial area in Eleusina nearby Athens (source: Google Maps).
Figure 1. The sampling site and the industrial area in Eleusina nearby Athens (source: Google Maps).
Atmosphere 13 00615 g001
Figure 2. Spatial representation of the sampling sites in the municipality of Elefsina (1, 2, 3, Source: Google Maps).
Figure 2. Spatial representation of the sampling sites in the municipality of Elefsina (1, 2, 3, Source: Google Maps).
Atmosphere 13 00615 g002
Figure 3. Percentage of cancer risk assessment (C.I method).
Figure 3. Percentage of cancer risk assessment (C.I method).
Atmosphere 13 00615 g003
Table 1. Description of Sampling and Coordinating Points.
Table 1. Description of Sampling and Coordinating Points.
Sampling Point CodeSampling Point FeaturesRegionCoordinates
XY
E1Primary School MandraElefsina (Mandra)38°3′7″ B23°31′35″ A
E2Primary School I ElefsinaElefsina38°3′10″ B23°31′51″ A
E3Primary School II ElefsinaElefsina38°2′26″ B23°32′4″ A
Table 2. Number of samples per sampling site.
Table 2. Number of samples per sampling site.
Sampling SitePM2.5PAHVOCsHeavy Metals
E1934213143
E2442211522
E343278425
Total1809133090
Table 3. Average concentration of measured pollutants in Elefsina area.
Table 3. Average concentration of measured pollutants in Elefsina area.
Chemical SubstanceConcentration (ng/m3)
PAHsBenzo(a)pyrene0.059
Benz(a)anthracene0.053
Benzo(b)fluoranthene0.289
Benzo(k)fluoranthene0.079
Chrysene0.108
Dibenzo(a,h)anthracene0.023
Indeno(1,2,3-c,d)pyrene0.135
1-methylnaphthalene0.577
Heavy metalsNickel4.380
Lead8.031
VOCBenzene1390
Table 4. Relative potency factor, (EPA, 2010) [35].
Table 4. Relative potency factor, (EPA, 2010) [35].
PAHRelative Potency Factor
benzo[a]pyrene1.000
benz[a]anthracene 0.200
benzo[b]fluoranthrene0.800
benzo[k]fluoranthrene0.030
chrysene0.100
dibenz[a,h]anthracene10.00
Indeno[1,2,3-c,d]pyrene0.070
benzo[g,h,i]perylene0.009
fluoranthene0.080
Table 5. Potential cancer risk assessment for Elefsina area.
Table 5. Potential cancer risk assessment for Elefsina area.
Lifetime Cancer Risk Cases/106
MethodologyChildren 0–2 YearsChildren 2–16 YearsWomen 16–80 YearsMen 16–80 YearsΣ0–80 yrs
(24 h/day)
14 h/Day24 h/Day14 h/Day24 h/Day
A.I3.150.660.080.130.080.143.95 1
A.II0.310.220.080.130.080.140.68 1
A.III5.27.212.47.813.418.6 1
A.IV0.59 (UR OEHHA, 2005)0.59
46.66 3
46.66 (UR WHO, 2010)
B.I1.648.023.233.84 2
B.II0.050.05 3
B.III2.35.49.35.49.311.7 1
B.IV2.35.49.311.7 2
C.I68.4447.8016.7828.7718.1131.04147.3 1
C.II68.4447.8016.7828.7718.1131.04147.3 1
1 Calculation of total inhalation lifetime cancer risk were used the estimated cancer risk from all pollutants for ages intervals children0–2, children2–16 (or children0–16) and Men16–80 (exposure duration 24 h/day). 2 According to these methodologies, there were not differences to estimations for the gender of adults. 3 These methodologies did not include any separation for the age of the people.
Table 6. Non-cancer risk assessment for Elefsina area.
Table 6. Non-cancer risk assessment for Elefsina area.
Hazard Quotient
MethodologyChildren 0–2 YearsChildren 2–16 YearsWomen 16–80 YearsMen 16–80 YearsHazard Index
(24 h/day)
14 h/Day24 h/Day14 h/Day24 h/Day
B.IHQ for each Target organ System(s) & Acute HI (OEHHA, 2015)
B.II3.59 × 10−62.33 × 10−65.92 × 10−6 *
B.III0.52690.30740.52690.30740.52691.05 **
B.IV0.52690.30740.52691.05 **
* In this methodology there was no separation of gender and duration of daily exposure people were exposed. ** For the calculation of total inhalation lifetime non-cancer risk, the estimated non-cancer risk from all pollutants for ages intervals Children0–2, Children2–16 (or Children0–16) and Men16–80 (exposure duration 24 h/day) was used.
Table 7. HQ for each Target organ System(s) & Acute HI (OEHHA, 2015).
Table 7. HQ for each Target organ System(s) & Acute HI (OEHHA, 2015).
SubstanceC (μg/m3)Total HQAcuteImmune SystemReproductive/DevelopmentHematologic SystemRespiratory SystemEyesNervous System
Benzene1.395.1 × 10−21.72 × 10−21.72 × 10−21.72 × 10−2
Toluene3.569.62 × 10−52.41 × 10−52.41 × 10−52.41 × 10−52.41 × 10−5
Xylenes 13.081.40 × 10−44.67 × 10−54.67 × 10−54.67 × 10−5
Xylenes 22.761.25 × 10−44.18 × 10−54.18 × 10−54.18 × 10−5
Xylenes 37.533.42 × 10−41.14 × 10−41.14 × 10−41.14 × 10−4
Total Acute Hazard Index1.72 × 10−21.72 × 10−21.72 × 10−22.27 × 10−42.27 × 10−42.28 × 10−4
Table 8. Table of methodologies including parameters.
Table 8. Table of methodologies including parameters.
MethodologyGenderAgesAFIR & BWExposure
(h/Day)
Exposure
(Day/Years)
Exposure
(Years)
A.IXX
A.IIXXX
B.IXX
B.IIXXX
A.IIIX
B.IIIX
B.IVXXX
A.IVXXXXXXX
C.IXXX
C.IIXXX
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pachoulis, M.; Maggos, T.; Panagopoulos, P.; Dasopoulou, M.; Balla, D.; Stamatelopoulou, A.; Manousakas, M.I.; Eleftheriadis, K.; Saraga, D. Population Health Risks Assessment from Air Pollution Exposure in an Industrialized Residential Area in Greece. Atmosphere 2022, 13, 615. https://doi.org/10.3390/atmos13040615

AMA Style

Pachoulis M, Maggos T, Panagopoulos P, Dasopoulou M, Balla D, Stamatelopoulou A, Manousakas MI, Eleftheriadis K, Saraga D. Population Health Risks Assessment from Air Pollution Exposure in an Industrialized Residential Area in Greece. Atmosphere. 2022; 13(4):615. https://doi.org/10.3390/atmos13040615

Chicago/Turabian Style

Pachoulis, Michail, Thomas Maggos, Panagiotis Panagopoulos, Maria Dasopoulou, Dimitra Balla, Asimina Stamatelopoulou, Manousos Ioannis Manousakas, Konstantinos Eleftheriadis, and Dikaia Saraga. 2022. "Population Health Risks Assessment from Air Pollution Exposure in an Industrialized Residential Area in Greece" Atmosphere 13, no. 4: 615. https://doi.org/10.3390/atmos13040615

APA Style

Pachoulis, M., Maggos, T., Panagopoulos, P., Dasopoulou, M., Balla, D., Stamatelopoulou, A., Manousakas, M. I., Eleftheriadis, K., & Saraga, D. (2022). Population Health Risks Assessment from Air Pollution Exposure in an Industrialized Residential Area in Greece. Atmosphere, 13(4), 615. https://doi.org/10.3390/atmos13040615

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop