Human health impact assessment and temporal distribution of trace elements in Copșa Mică- Romania

The present study aims to analyze the temporal variations of PM10 and to assess the health risk indexes caused by trace elements from particulate matter (PM10) via inhalation, ingestion, and dermal absorption by adults and children in Copșa Mică (Romania) during 2009–2019. The results revealed a high multi-annual mean concentration of PM10 and trace elements. The analyzed air pollutants showed a decreasing trend during the studied years, therefore 44.11%, 43.48%, 36.07%, 16.02%, and 15.80% lower values were observed for As, Cd, Ni, PM10, and Pb, respectively, due to environmental regulations. The daily exceedance percentage of Pb and Cd was very high, representing 21.74% and 11.26%, followed by PM10 and As concentrations with 4.72% and 3.92%. The ratio between the trace element concentration measured in Copșa Mică and the country average was 2.46, 4.01, 2.44 and 10.52 times higher for As, Cd, Ni and Pb. The calculated Hazard Quotient values via inhalation were higher than the safe limit (1), which accounted 1.81, 3.89 and 4.52, for As, Cd and Ni, respectively, indicating that the trace elements might present a non-carcinogenic risk to both adults and children. Furthermore, the concentration of all studied trace elements in Copșa Mică showed cancer risk for adults via inhalation and dermal absorption as well.

Particulate matter (PM 10 ) is an important air pollutant consisting of small particles with an aerodynamic diameter less than or equal to a nominal 10 µm having a significant impact on human health, causing serious disease (respiratory, cardiovascular) and premature death worldwide [1][2][3] . The particulate matters originated from different (industrial, traffic) emissions may contain several toxic trace elements, such as As, Cd, Cr, Cu, Zn, Pb and Ni. These toxic elements may enter the body via three different entry routes: ingestion, inhalation and skin absorption 4 . Many studies have been published on the topic of inhalable particles associated with trace elements that have an increased effect on lung and cardiopulmonary morbidity and mortality 5,6 . The annual acceptable limit for the most important air pollutants was established in the Air Quality Standards including PM 10 (20 µg m -3 ), Pb (0.5 µg m -3 ), As (6 ng m -3 ), Cd (5 ng m -3 ) and Ni (20 ng m -3 ) as well 7 .
The trace elements concentration is strongly associated with the source types 8 , hence the most common source of As, Cd and Pb, are coal combustion 9 and metal smelting industry, which are also responsible for the emission of trace elements 10,11 . The studied town Copșa Mică was best known as one of the most polluted cities in Europe where the carbon black, lead and zinc production industry heavily polluted the ecosystems. Due to this pollution, in Copşa Mică the life expectancy is with 9 years below the Romanian average, where the spread of lung cancer, lead poisoning, bronchitis, rickets, finger-twitching, learning difficulties, impotence, asthma, stunted growth, depression and alcoholism are very common morbidities 12, 13 . Furthermore, the temporal variations of PMs and trace elements highly depend on the emission sources and specific meteorological conditions 14,15 .
The trace elements concentration of particulate matter measured near industrial areas has been studied and reported worldwide, in the USA by Landis et al. 16 and Xia and Gao 5 , in Korea by Kim et al. 17 , in Greece by Manalis et al. 18 , in Italy by Mazzei et al. 19 , Prodi et al. 20 and Amodio et al. 21 , in Turkey by Cetin et al. 22 , in Spain by Querol 23 , in China by Zhang et al. 24 Tian et al. 25 and Du et al. 4 , in Taiwan by Jiun-Horng et al. 26 , in Poland by Pastuszka et al. 27 , in France by Hleins et al. 28 , in Australia by Mohiuddin et al. 29 , in the United Kingdom by Taiwo et al. 30 , in Romania by Dunea et al. 6 and Proorocu et al. 31 . According to the findings reported, the level of trace elements near industrial areas can reach values which may represent serious health risks at local and www.nature.com/scientificreports/ regional levels. Therefore, it is mandatory to decipher the effects on human health, and in many cases the complex interactions are yet to be discovered. The main objective of this research study was to analyze the heavy metals (As, Cd, Ni, Pb) temporal distribution determined from the PM 10 and to assess the human health effect using the Hazard Quotient (HQ) and cancer risk (CR) methods, where three exposure pathways were considered: inhalation, ingestion and dermal absorption.

Study area
Sampling site. Copșa Mică (46°6′45″N 24°13′50″E) is a small industrial town located in Sibiu County, Romania, and was the most polluted city in Europe in 1991 due to heavy emissions from zinc, cadmium and lead refineries and carbon black plants. According to the literature, the uncontrolled emission over the years had a negative effect on air quality, soil, water, plant, animal and humans as well 32 . As a consequence of the uncontrolled air and soil pollution in recent decades has a demonstrable effect, since the trace elements were detected in the food chain, namely in honey 33 .

Materials and methods
In the present study, the temporal variations of daily PM 10 concentration and associated trace elements As, Cd, Ni, Pb were examined from January 2009 to August 2019. The raw data regarding daily pollutants concentrations were obtained from the National Air Quality Monitoring Network (www. calit ateaer. ro). Therefore, during the studied period 3049 PM 10 , 2272 As, 2300 Cd, 2270 Ni and 2272 daily Pb sample was analyzed. Romania's climate is temperate continental transitioning, with four distinct seasons. In order to evaluate the seasonality aspect of pollution, the concentrations were grouped into four groups as follows; spring (March, April and May), summer-warm period (June, July and August), autumn (September, October and November) and winter-cold period (December, January and February).
The industrial monitoring station RO-SB-3(RO0186A) is situated in Copșa Mică with coordinate: latitude 46.11°N and longitude: 24.23°E, and 286 m asl. Trace elements were determined from the PM 10 particulate matter fraction according to the SR EN 14,902 "Ambient air quality from the fraction of PM 10 particulate matter " reference method. As described in this method, the collected particulate matter via aspiration was pre-treated in a microwave oven in a closed vessel, using concentrate nitric acid and hydrogen peroxide (30%). From the obtained solution, the trace element concentrations were measured with inductively coupled plasma mass spectrometry (ICP-MS).

Statistical analysis.
In order to decipher the temporal differences, descriptive statistics, monthly and annual trends were calculated, and the results were presented using box-plot diagrams. The study uses the 25th percentile of data as background levels. The box-plot used the median, the Q1 and Q3 quartiles, and the minimum and maximum data points to convey the level, spread, and symmetry of a distribution of data values 34 .
Using R (R 3.6.2) statistical program, the Spearman correlation analysis was carried to understand the relationship between the monthly mean pollutants (As, Cd, Ni, Pb, PM 10 ) concentration and meteorological parameters (precipitation quantity, temperature, relative humidity, wind speed). The hierarchical cluster analysis method (Centroid Linkage, Correlation Coefficient Distance) was used to classify the PM 10 and the trace elements (As, Cd, Ni, Pb) in into groups or clusters based on their similarities. The results were analyzed with Minitab 17.3.1 statistical software and the outcomes are presented in a dendrogram, while for the Principal Component Analysis the IBM SPSS Statistics 22. was applied.
The health risk posed by heavy metals in PM 10 . According to the United States Environmental Protection Agency (USA-EPA), the human health effect caused by the heavy metals from PM 10 was calculated via ingestion (CDI-chemical daily intake), inhalation (EC-exposure concentration), and dermal absorption (DAD-dermal absorption dose (EPA) for both groups children (ch) and adults (ad) 35 . Regarding the concentration of the pollutant, the multiannual mean value was used throughout the entire studied period (2009-2019).
(1) The trace elements (As, Cd, Ni, Pb), non-carcinogenic and carcinogen effect on human health was separately calculated for the children and adults, based on the different characteristics and exposure time. In order to assess the non-carcinogenic and carcinogenic risks the hazard quotient (HQ) and cancer risk (CR) approaches were used. For interpretation of HQ values: less than 1, no adverse health effects are expected as a result of exposure, while in the case of HQ greater than 1, then adverse health effects are possible. CR represents the increased probability of incident of tumor diseases above the general average due to the impact of carcinogenic compound's effects. The evaluation of CR for carcinogenic chemicals is taken into consideration and represents a risk, when the values varies from 10 −4 to 10 −6 , representing that the cancer development during a human lifetime (70 years) is 1/10,000 or 1/1,000,000, respectively 35 . Values lower than 10 −6 for individual chemicals and pathways show no cancer risks. Generally speaking, a cumulative cancer risk higher than 10 −4 is not accepted, and the maximum tolerable value is 10 −5 .

Results and discussions
Statistical analysis. Through the statistical analysis the daily datasets were used. Within the studied period, the average concentration of PM 10 in Copșa Mică was found to be 24.62 µg m -3 , which is 23.1% higher than the WHO's acceptable limit (20 µg m -3 ) 7 . The average pollutant concentration between 2009 and 2019 in the studied area was 1.65 ng m -3 for As, 2.37 ng m -3 for Cd, 5.5 ng m -3 for Ni and 0.32 µg m -3 for Pb, respectively. During the studied period, the highest trace element concentration in PM 10 was recorded for Pb, meanwhile the lowest one was observed in case of As (Table 1).
The results revealed that the exceedance percentage, based on the WHO Air Quality guideline, of Pb (> 0.5 µg m -3 ) and Cd (> 5 ng m -3 ) was very high. The daily permissible concentration was exceeded with 21.74% and 11.26%, respectively, followed by the PM 10 (> 20 µg m -3 ) and As (6 ng m -3 ) concentration with 4.72% and 3.92%, while the lowest exceedance percentage was found for Ni (> 20 ng m -3 ), representing 0.79%.
Temporal distribution of pollutants. Monthly and annual pollutant distribution. The multiannual monthly variation in PM 10 and trace elements concentrations in Copșa Mică are presented in Fig. 1. According to the results, for all studied elements (PM 10 , As, Cd, Pb) except Ni, the minimum concentration was found in the warm period, meanwhile the maximum concentration in the cold period.
The seasonal variation of the pollutants is strongly related to different factors, mainly to the meteorological condition and emission sources, therefore the higher PM 10 concentrations in the winter period can be explained (4) HQ ing = CDI/RfDo    (Fig. 2.). One of the main reasons of decreasing air pollution level in Romania and Copșa Mică as well is, that in 2007 the country joined the European Union and became a full member state, hence stricter Environmental Protection Regulations were implemented in order to address the industrial pollution issues. During this period a modernization process has also taken place and part of the industry has been closed.
Where: CM and RO represent the multiannual average concentrations in Copșa Mică and Romanian. Compared to the multiannual country average, the Pb concentration in Copșa Mică was 10.52 times higher than the country average; even though all countries worldwide have phased out using leaded petrol by 2012, in Copșa Mică the Pb concentration determined from the PM 10 was very high. The Pb pollution during the last decades caused by Sometra (specialized in lead, zinc and other nonferrous metals) still has a visible effect even nowadays. Since via soil resuspension, very high Pb concentrations were measured from PM 10 samples in the studied area 39 . Moreover, similar to Pb the Cd concentration was 4.01 times higher compared to the country average, followed by the As and Ni with 2.46 and 2.44 ratio, respectively. However, the PM 10 concentration did not differ significantly from the national average, was only 1.054 times higher than the average.
Analyzing the data on an annual breakdown, it can be seen that the highest pollutant concentration was detected in 2015. According to the 2015 annual report, we found that during this period massive rehabilitation www.nature.com/scientificreports/ works were carried out in different sectors including the water-sewerage-, natural gas network, and local streets/ roads repairs as well 13 . Thanks to these activities dust emissions was very high, thus the particulate matter concentration in the air was sharply growing, which was the main cause of high PM 10 levels recorded in 2015, hence the Ni concentration was also increasing. Besides the fact the precipitation amounts were considerable in 2015, they had a torrential character, after which they were accompanied by periods with persistent anticyclonic systems that induced an accentuated static stability, causing frequent periods of thermal inversion. Therefore, the washout effect on trace elements concentrations was not as significant as expected. During thermal inversions, the colder air layers were blocked under the hot air, thus preventing the formation of convection currents (ascending) and blocking the emitted noxious substances, which favors the horizontal distribution and accumulation of pollutants, especially under stable condition (no vertical mixing, air stagnation). In order to decipher the level of pollution in Copșa Mică and the relationship with different locations worldwide, the multiannual trace element concentration measured in the study site was compared with those reported in different regions of the world ( Table 2).
The comparative results revealed that the measured As concentration in Copșa Mică (1.65 ng m -3 ) was higher than what was observed in Appalachia (USA) and Spain. The Cd (2.37 ng m -3 ) concentration was the second highest after India, and the Ni concentration (5.5 ng m -3 ) was higher than what was reported in Spain. Moreover, the Pb concentration was by far the highest in Copșa Mică (320 ng m -3 ).
Correlation analysis. Spearman correlation analysis was carried out between the pollutant concentration (PM 10 , As, Cd, Ni, Pb) and meteorological parameters (Prec, Temp, RH, Wind speed), using the monthly average concentrations to identify the common sources 44 . The correlation coefficients between two studied parameters were considered significant if the P < 0.05 and r ≥ + 0.27, and r ≤ − 0.27 (Fig. 3).  Figure 3. Spearman correlation coefficient matrix. The figure was prepared using R (3.6.2) statistical program. www.nature.com/scientificreports/ As expected, the results show that the most significant correlation was between As-Ni and As-Cd (r = 0.46, r = 0.68). Furthermore, significant positive correlation was found between PM 10 and As, Cd, Pb (r = 0.48, r = 0.46, r = 0.32) and between Pb and Cd (r = 0.34) as well. The Spearman correlation shows the extent to which the magnitude of one variable determines the magnitude of the other variable, as well as the direction and strength of the relationship, hence from the significant positive correlation we can conclude whether the two variables are related or not. The correlation level of Cd-Ni, and Ni-Pb was lower than the significant level, indicating that sources of these elements are different and more diverse. The correlation matrix also indicates that significant negative correlation was between the temperature and the Cd, Pb and PM 10 concentrations (r = − 0.31, r = − 0.46, r = − 0.42), which could be attributed to the thermal inversion. Due to the fact that the days with precipitation and no precipitation were not analyzed separately, the negative correlation between the precipitation and As, Cd and Pb was not significant. The wind speed also showed a significant negative correlation with the As and Pb (r = − 0.47, r = − 0.58).
Cluster and principal component analysis. Using the monthly mean concentration, hierarchical cluster analysis was performed for trace elements (As, Cd, Ni, Pb) and PM 10 to evaluate the potential contributing sources of heavy metals. According to the hierarchical cluster analyses, the variables were distributed in two different clusters (Fig. 4.). The HCA also revealed that As and Cd belong to 1.1 sub-cluster, which means that they are coming from different sources and industrial emissions 45 . The Pb and PM 10 form a separate sub-cluster, 1.2. and 1.3, respectively, meaning a different source, while Ni is in cluster 2, which means that can derive from motor vehicle exhaust 46 .
In order to identify the origin as well as the common sources of heavy metals from the PM 10 samples, Principal Component Analysis was carried out. In the variable statement we include the first two principal components. The adequacy of the Kaiser-Meyer-Olkin (KMO) measure of sampling was 0.67, followed by the execution of the PCA, meanings that the tested samples show medium adequacy.
According to the results, two components were extracted from the component matrix, accounting for 70.66% of the overall variance (Table 3, Fig. 5). Factor 1 contains As, Cd, Pb, PM 10 and represents 48.15% of the total variance, while factor two was represented by Ni.
The most important sources of As and Cd, are related to the coal-burning, diesel fuel, lubricating oil and tire wear 47 while the most important Pb sources are related to soil dust 48 . The origin of Ni represents factor 2 (0.797), which is due to oil burning 49 . According to the annual reports, regarding the state of the environment, the main source of county-level Pb emission is the metallurgical industry-92%-including the lead production industry from Copșa Mică (SOMETRA), the industry of other mineral products (brick), followed by domestic and institutional heating and traffic 50 . On the other hand, the highest share in Ni emission, at county-level, belongs to the industry of other mineral products, such as (brick production)-79%, followed by commercial and institutional heating, domestic heating, road transport, equipment and mobile machinery used in industries, constructions and agriculture as well.
Health risk assessment of toxic metals in PM 10 . Non-cancer risk and cancer risk assessment. For the three studied exposure pathways (ingestion, inhalation, dermal absorption), the calculated HQ and CR values are presented in Table 4. According to the HQ results, the largest amount of trace elements is generally absorbed by the body through inhalation, while via ingestion and dermal absorption, for both adults and children, the HQ www.nature.com/scientificreports/ values were lower than the safe limit. Duo to the fact that the data for inhalation exposure route were not available, the HQ values for Pb were not calculated. Through the inhalation pathway, the HQ values were higher than the safe limit for both adults and children (= 1). The highest non-carcinogenic risk was detected for Ni (4.52), in case of both groups adults and children, while the non-carcinogenic risk value for Cd and As was of 3.89 and 1.81, respectively. However, by taking into consideration the sum of the three trace elements (HI) the intensification parameter was 1.02E + 01, showing non-carcinogen health risk by inhaling a mixture of trace elements. The results revealed that the carcinogen risk via inhalation for adults in the case of As, Cd and Pb was higher than the acceptable limit (1E−06). The total cancer risk for adults and children via inhalation was observed to   www.nature.com/scientificreports/ be 9.34E−05 and 2.34E−05, respectively. There is a significant difference between children and adult cancer risk due to their activity and exposure time differences 51 . Via dermal absorption and inhalation the cancer risk found was higher than the safety limit (1E−06); on the other hand, adults had a higher probability of experiencing carcinogenic risk compared to the children. According to the calculations, in the case of As, adults show cancer risk via dermal absorption (6.55E−03) and inhalation (4.00E−05) while the Cd cancer risk for adults was 1.28E−03 and 2.40E−05, respectively. Via dermal absorption and inhalation, the sum of all the elements under consideration was 9.42E−02 and 9.34E−05 for adults, showing cumulative cancer risk, if the exposure is to a mixture of elements. The reason why the dermal absorption is so high may be related to the big skin contact surface. Furthermore, the cancer risk via ingestion for children and adults was below the minimum acceptable level (1E-06) for all trace elements and cumulative values as well, thus presenting a negligible carcinogen risk.

Conclusions
The present study focused on the temporal distribution of trace elements originated from the PM 10 in Copșa Mică, and to assess the human health effects. During the studied period (2009-2019), the temporal distribution of PM 10 , As, Cd and Pb concentrations variation was significant, the minimum value was recorded in the warm season and the maximum concentration in the cold period. The multiannual mean concentration of PM 10 was 24.62 µg m -3 exceeding the EU Air Quality admissible concentration (20.00 µg m -3 ) by 23.1%. Pb was found as the dominant metal, followed by the Ni, Cd and As. The Pb concentration was higher than the acceptable annual level in over one-fifth of the total studied days. Although more than one decade ago (2007) the Sometra Industrial company closed most of the activities, the sign of past pollution is still present and can be demonstrated even nowadays, since the trace elements in Copșa Mică such as As, Cd, Ni and Pb, are 2.46, − 4.01, 2.44, 10.52 times higher than the country average. Each element analyzed in this study could present non-carcinogenic risk via inhalation while taking together (multiple elements) a significant carcinogenic risk was revealed on adults. Based on the health risk assessment calculation, the highest Hazard Quotient was found via inhalation while exposed to mixture of trace elements, the summarized Hazard Index showed an increased level. Via dermal absorption and inhalation, a potential carcinogenic risk, exceeding the carcinogen acceptable level, was detected indicating an elevated risk of cancer to the inhabitants (adults) in the studied area. In this context, Romania needs to improve its environmental protection measures and procedures to remediate the heavily polluted industrial regions. www.nature.com/scientificreports/ Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.