Next Article in Journal
Spatial Differences in the Effect of Communities’ Built Environment on Residents’ Health: A Case Study in Wuhan, China
Next Article in Special Issue
The Influence of Environmental Polycyclic Aromatic Hydrocarbons (PAHs) Exposure on DNA Damage among School Children in Urban Traffic Area, Malaysia
Previous Article in Journal
A Global Perspective of Racial–Ethnic Inequities in Dental Caries: Protocol of Systematic Review
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Distribution of Air Pollution, Hotspots and Sources in an Urban-Industrial Area in the Lisbon Metropolitan Area, Portugal—A Biomonitoring Approach

1
Centro de Ciências e Tecnologias Nucleares, Instituto Superior Técnico, Universidade de Lisboa, Estrada Nacional 10, 2695-066 Bobadela, Portugal
2
ESTSetúbal/IPS and CINEA, IPS Campus, Polytechnic Institute of Setúbal, 2914-508 Setúbal, Portugal
3
CESAM—Centre for Environmental and Marine Studies, Department of Environment and Planning, University of Aveiro, 3810-193 Aveiro, Portugal
4
Laboratory for Heritage Science, Institute for Nuclear Research, H-4026 Debrecen, Hungary
*
Author to whom correspondence should be addressed.
Int. J. Environ. Res. Public Health 2022, 19(3), 1364; https://doi.org/10.3390/ijerph19031364
Submission received: 23 November 2021 / Revised: 19 January 2022 / Accepted: 21 January 2022 / Published: 26 January 2022
(This article belongs to the Special Issue 2nd Edition of Integrated Human Exposure to Air Pollution)

Abstract

:
This study aimed to understand the influence of industries (including steelworks, lime factories, and industry of metal waste management and treatment) on the air quality of the urban-industrial area of Seixal (Portugal), where the local population has often expressed concerns regarding the air quality. The adopted strategy was based on biomonitoring of air pollution using transplanted lichens distributed over a grid to cover the study area. Moreover, the study was conducted during the first period of national lockdown due to COVID-19, whereas local industries kept their normal working schedule. Using a set of different statistical analysis approaches (such as enrichment and contamination factors, Spearman correlations, and evaluation of spatial patterns) to the chemical content of the exposed transplanted lichens, it was possible to assess hotspots of air pollution and to identify five sources affecting the local air quality: (i) a soil source of natural origin (based on Al, Si, and Ti), (ii) a soil source of natural and anthropogenic origins (based on Fe and Mg), (iii) a source from the local industrial activity, namely steelworks (based on Co, Cr, Mn, Pb, and Zn); (iv) a source from the road traffic (based on Cr, Cu, and Zn), and (v) a source of biomass burning (based on Br and K). The impact of the industries located in the study area on the local air quality was identified (namely, the steelworks), confirming the concerns of the local population. This valuable information is essential to improve future planning and optimize the assessment of particulate matter levels by reference methods, which will allow a quantitative analysis of the issue, based on national and European legislation, and to define the quantitative contribution of pollution sources and to design target mitigation measures to improve local air quality.

Graphical Abstract

1. Introduction

Industrial emissions have an important contribution to particulate matter (PM) in urban-industrial areas [1,2]. To solve environmental problems in this type of areas it is crucial to understand and identify prevailing emission sources to promote targeted and successful mitigation measures.
It is known that particulate matter may have adverse health effects, since PM may contain potentially toxic elements (PTEs) [3,4,5]. Thus, the determination of the chemical composition of PM aims not only to determine the sources of the PM sampled at receptor sites but also to identify its potential health hazards.
In the last years, the episodic deposition of airborne particulates onto homes and properties was a common complaint from the residents of Aldeia de Paio Pires-Seixal (municipality of Seixal of the Lisbon metropolitan area), an urban-industrial area in Portugal [6]. To answer the concerns raised by the population, the local Council of Seixal promoted a set of actions, such as the assessment of the chemical composition of the settled dust to determine its sources and potential health hazards [7]. In this previous study, the settled dust was collected in January 2019, characterized by micro-PIXE, and the influence of the steel industries was identified due to the content of Fe, Cr, and Mn, along with a minor traffic influence.
PM emissions from steelworks are a complex mixture of stationary and diffuse emissions, associated with the process and with operations such as stocking and transportation of raw materials and slags [1,8,9]. Atmospheric PM from iron and steel industries have high concentrations of As, Cd, Cr, Cu, Fe, Mn, Ni, Se, V, and Zn [10,11,12]. Moreover, steelworks are often close to other industries and traffic zones, making it difficult to distinguish between the contributions of the processes.
Traditionally, air pollution studies are performed through instrumental techniques that are limited to a small number of sampling stations [13]. Biomonitoring offers several advantages over standard sampling methods since biomonitors can be used in vast areas with many monitoring points while requiring little maintenance and are low cost [14]. Biomonitors can be native organisms in the ecosystem or they can be collected in an unpolluted site and transplanted to the area of interest [15,16].
Lichens have the ability to accumulate elements, in correlation with atmospheric levels. In fact, since lichens have no roots, they are dependent on the atmosphere for the uptake of water and mineral substances [17]. Lichens are extensively used in biomonitoring studies [13,18,19,20,21,22] and for spatial mapping of air contaminants [14,23,24], allowing the identification of pollution hotspots.
Continuing the set of actions promoted by the local Council of Seixal to understand the population exposure and the associated health risks of the settled dust events, the present study was conducted for the Council of Seixal and it aimed to assess the spatial distribution of air pollution, to identify its hotspots and its sources, by the use of the biomonitoring technique of air pollution using transplanted lichens in the Aldeia de Paio Pires-Seixal. The elemental characterization of the transplanted lichens was performed to assess the potential pollution sources, through the analysis of their chemical tracers. The Micro-X-ray Fluorescence technique was applied since it has been previously used to quantify metals in lichens and has several advantages: it is nondestructive, presents no need for sample digestion, and accuracy and reproducibility are at least equal to spectroscopic methods [25].

2. Materials and Methods

2.1. Study Area

This study was carried out in the parish of União das Freguesias do Seixal, Arrentela e Aldeia de Paio Pires (UFSAAPP), in the municipality of Seixal (Portugal), which is located in the peninsula of Setúbal and it is a part of the Lisbon Metropolitan Area (Portugal), next to the Nature Reserve of the Tagus Estuary (Figure 1). The municipality of Seixal is one of the most densely populated municipalities in Portugal, with 167,294 inhabitants in 95.5 Km2 [26] (PORDATA, 2021).
Besides including a significant network of highways and national roads (A2-IP7, A33, and EN10), the study area (UFSAAPP) comprises small and medium-sized industries and an industrial park where the following facilities are located (as shown in Figure 1): industry A–that is a steelwork that manufactures galvanized sheet metal (with Cr passivation) and cold rolled sheet, with an installed capacity of 800,000 tons per year [27]; industry B–that is steelwork with an installed capacity of 215 tons per hour for the production of steel, where the main processes consist of an Electric Arc Furnace for steelmaking and hot rolling [28]; industry C–that it is a lime factory, which manufactures lime by calcination of limestone in a coke kiln [29]; and industry D that is a company focused on metal waste (such as iron and aluminum) management and treatment (recovered iron is a raw material for industry B).

2.2. Transplantation and Sampling

Samples of the lichen Flavoparmelia caperata (L.) Hale were collected from olive trees at about 1.5 m above the soil, in Montargil (39°03′24″ N, 8°10′36″ W) (as shown in Figure 1), on 22 January 2020. Montargil is a rural area considered clean from an air pollution point of view [31]. Lichens were collected using powderless gloves and stored temporarily in paper envelopes to be transported to the laboratory. In the laboratory, four lichen samples were separated randomly, as reference base levels. Lichens were placed in nylon mesh bags and exposed in the study area, fixed to an appropriate substrate using nylon string, at about 1.8 m above the soil (Figure 2). The exposure period of the samples was approximately four and a half months (total of 137 days), from 1 February to 17 June 2020.
A georeferenced grid of 4.55 km × 6.82 km, with 77 cells of 650 m × 620 m, was drawn for the lichens exposure, between the coordinates −9.11, 38.65 and −9.05, 38.59 (upper left and lower right corner of the grid), corresponding to the extremes of UFSAAPP (Figure 3).
Lichens samples were distributed over the study area based on two strategies: (1) participation of the local population by exposing the transplanted lichens in outdoor of their dwellings, and (2) to have good coverage of the study area, lichen samples were exposed in selected sites to cover the established grid cells that were not covered by the strategy 1. For the population engagement, a public presentation session was held on 1 February 2020, to explain the study and to identify volunteers. A total of 88 lichen samples were exposed but only 63 could be retrieved at the end of the exposure period and subsequently analyzed (retrieved lichen samples are shown in Figure 3, left).

2.3. Meteorological Data

Meteorological data for the exposure period were obtained from the Barreiro-Lavradio weather station next to the study area (38°40′28″ N, 9°2′51″ W, it is the nearest weather station, and it is located at 5.5 km from the study area) and it was supplied by the Portuguese Institute of Sea and Atmosphere (Instituto Português do Mar e da Atmosfera-IPMA).
During the exposure period, the average temperature in the study area was 16.9 ± 3.2 °C (ranging from 9.7 to 25.1 °C) and the average relative humidity was 80 ± 8% (ranging from 63 to 98%), with a total of 36 rainy days registered.
The wind rose for the exposure period is presented in Figure 3, where a predominance of winds from North (N), Northwest (NW), and West (W) was registered.

2.4. Chemical Analysis

At the laboratory, lichens were cleaned and cleared of extraneous material (such as dust, bark remaining, leaf debris, fungus contamination, and degraded material), and rinsed for 5 s with demineralized water, freeze-dried, and ground into powder in a ball mill with PTFE capsules under liquid nitrogen, for homogenization [32]. Pellets with an average thickness of 12 mm, were then prepared for elemental characterization by micro-X-ray Fluorescence [33]. For each lichen sample, a total of three pellet replicates were done. Figure S1 provides an overview of the procedure for the preparation of pellets.
Micro-X-ray Fluorescence (micro-XRF) analysis was conducted in the Laboratory for Heritage Science (Debrecen, Hungary), ATOMKI, Debrecen, Hungary using a Bruker M4 Tornado micro-XRF equipment (Bruker Corporation, Billerica, MA, USA). All measurements were made in a vacuum (~20 mbar). Rh excitation source set to 50 kV with 300 µA current was applied. The beam was focused down to 25 µm using polycapillary optics. To set the optimal measurement conditions, different filter combinations were applied on selected samples: (1) no filter, (2) 12.5 µm Al, (3) 100 µm Al + 25 µm Ti [34]. Test measurements to detect elements like Cd, Sn, Sb were carried out with W (50 kV, 700 µA) excitation source. Since Cd, Sn, Sb were under detection limit even with the W tube, an Rh excitation source was applied for the measurements. To achieve the result for the widest possible elemental range (Na-U), measurements were performed without a filter.
X-ray spectra were collected simultaneously with 2 detectors (XFlash® SDD with 30 mm2 active surfaces, Be window). Spectra and maps of ~6 mm × 7 mm areas were acquired to get the best estimate of the average composition of the samples and to eliminate the possible fluctuations in concentrations due to sampling inhomogeneity. For each lichen sample, only one of the three replicate pellets was measured. Regarding the quality control of the measurements, randomly selected pellets were re-measured 2 more times on different days with the same settings. In the case of 5 random pellets, both sides were measured, and in the case of 6 samples, all of their replicate prepared pellets were analyzed (three replicates per sample). In all cases, the resulting concentration data were within 2 sigmas.
For the quantification, the fundamental parameter (FA) method [35] was applied using the MQuant built-in software [36]. The composition of cellulose (C6H10O5)n was set as an “unknown” matrix.
The elemental composition of exposed and unexposed biomonitors was assessed for a total of 20 elements (Al, As, Br, Ca, Co, Cr, Cu, Fe, K, Mg, Mn, Pb, Rb, S, Se, Si, Sr, Ti, Zn, and Zr).

2.5. Statistical Analysis

Statistical analysis was performed using STATISTICA software version 13. The variables in the data set exhibited a non-normal distribution, and so the Mann-Whitney U nonparametric statistics was applied [37], at a significance level of 0.05, for independent groups, to suggest whether samples come from the same population or not. One Sample Wilcoxon nonparametric test was also applied, at a significance level of 0.05, to indicate if there is a significant difference between the median of a sample group and an hypothesized value [37]. Spearman correlations [37] were used to understand the associations between parameters.
The enrichment factor (EF) is used for identifying the crustal and non-crustal origin of elements [38], and it is applied to particulate matter collected in filters [14,39], plants/lichens [32,40], and sediments [41]. The enrichment factor of the element was determined for each sample, using the reference values of soil composition (Table S1) [42] and taking Si as the crustal reference element [43], according to Equation (1):
E F X = X S i L i c h e n   X S i S o i l
where X S i L i c h e n   is the ratio between the concentrations of element X and Si in the lichen and X S i S o i l is the ratio between the reference concentrations of the element X and Si in Soil. To account for the local variation in the soil composition in the EF calculations, elements with EF values between 1 and 10 are considered to be of a crustal origin [32], while elements with EF values greater than 10 were considered to be enriched by other sources rather than the crust [44].
To compare the mean elemental concentrations assessed in the exposed lichens with the unexposed lichens, the calculation of the Contamination Factor (CF) was conducted, where CF for a specific element is the ratio between the exposed and the unexposed lichen [45]. CF is also known as Exposure to Control ratio, which has been divided in five classes [46], namely (i) 0–0.25: severe loss, (ii) 0.25–0.75: loss, (iii) 0.75–1.25: normal, (iv) 1.25–1.75: accumulation, and (v) >1.75: severe accumulation.
Geostatistical modeling maps were built using ArcGIS 10.1 software (ESRI, 2020, Redlands, CA, USA), for studying the geospatial distribution of the elemental concentrations in the study area. The mapping was performed using the IDW (Inverse Distance Weighted) tool, which has been applied in similar studies [45]. IDW calculates the cell value for the unmeasured location, averaging the sampled data around each processing cell; the closer a measured point is to the center of the prediction cell, the more weight it will have [47]. In the overall equation for IDW, v0 represents the estimated value at point 0 (Equation (2)), vi the value in the known point I, di is the distance to point I, S is the number of known points applied in the estimate, and k is assumed to equal to 2. A set of 12 neighboring samples was chosen to represent the spatial variation of the elements under study.
v 0 = i = 1 S v i 1 d i k i = 1 S 1 d i k
The range between the maximum and minimum concentrations obtained for each element was subdivided into five bands, corresponding to 0–20, 20–40, 40–60, 60–80, and 80–100 percentiles.

3. Results and Discussion

3.1. Elemental Characterisation

Table 1 presents the mean elemental concentrations in unexposed lichens and in the transplanted lichens after the exposure period. The major elements assessed in the exposed lichens were (by decreasing order): Ca (22.6%) > Fe (1.18%) > K (0.79%) > Si (0.62%) > Al (0.26%) > S (0.16%) > Mg (0.12%) > Ti (0.12%), with the remaining elements contributing with less than 0.10% to the total mass of the lichen samples (namely, Zn > Sr > Mn > Br > Cu > Pb > Zr > Se > Rb > Co > Cr > As).
Compared with studies reported in the literature that used the same lichen species in industrial environments, it is possible to conclude that the concentrations found in our study area were always higher than in Sines (Portugal), an industrial area, with heavy chemical industry, a coal power plant and a deep water harbor [48,49].
It is important to highlight that the exposure period of the transplanted lichens in the present study was between February and June 2020, a period when the COVID-19 pandemic situation reached Europe. In fact, due to the global COVID-19 pandemic situation, on 19 March 2020, Portugal decreed the state of national emergency, which included mandatory confinement, movement restriction of citizens, and closure of non-essential business. As a consequence of the nationwide lockdown that occurred from 19 March to 31 May 2020, concentrations in mainland Portugal presented a mean reduction of 41% for NO2 and 18% for PM10, considering the 20 air quality monitoring stations analyzed [50], following a tendency observed all over the world [51].
In the present study area, a drastic decrease of concentrations of PM10 and NO2 was registered in April 2020 (−40.3% and −44.0%, respectively), compared to the previous six years (a significant difference, with p-value = 0.000) [52]. In May 2020, a lower but still significant decrease of PM10 (−16.7%, p-value = 0.000) and NO2 (−17.4%, p-value = 0.048) was also registered. These reductions were accentuated again in June 2020 for both pollutants: PM10 (−38.0%, p-value = 0.000) and NO2 (−29.0%, p-value = 0.000). Despite the reduction of emissions from anthropogenic activities (such as traffic) due to the national lockdown (which reflects in pollutants, such as PM10 and NO2), the industries (namely, steelworks) that exist in the study area never stopped operating during this period. Therefore, the present biomonitoring study has the potential to assess and evaluate the contribution of these specific industries since they were the main sources of air pollution in the area during the exposure period of the transplanted lichens.

3.2. Enrichment and Contamination Factors

Figure 4 presents the mean EFs assessed for each element in the exposed lichens.
The elements Al and Mg were found to have a predominantly crustal origin (with EFs always below 10), while the remaining elements showed the contribution of non-crustal emissions to their levels. The elements Fe, Ti, and K, typically associated with a soil source [43], presented EFs slightly above 10 (with means of 10.5 ± 2.3, 11.7 ± 1.6, and 14.0 ± 2.9, respectively), indicating also anthropogenic contributions to their levels. A mix of different anthropogenic sources may have contributed to lichens enrichments, given the urban-industrial location of the study area.
Figure 5 provides the contamination factors regarding all the studied elements and it is observed a severe accumulation for Cr (CF = 7.11 ± 0.25) and an accumulation for Zn (CF = 1.75 ± 0.43), Mn (CF = 1.48 ± 0.50) and Pb (CF = 1.29 ± 0.67). However, a significant difference between exposed and unexposed elemental contents of lichens was only found for Cr (p-value = 0.034) and Zn (p-value = 0.022), but not for Mn (p-value = 0.064) neither Pb (p-value = 0.190). However, a significant difference between the medians of elemental contents of exposed and unexposed lichens was found for all the elements, except for Cu and Sr, as ascertained by the One Sample Wilcoxon test (at 0.05 significance level). Arsenic showed a loss (CF = 0.51 ± 2.02), with significant difference between exposed and unexposed levels (p-value = 0.031), as well as K (CF = 0.89 ± 0.18, p-value = 0.033). There is evidence that metals (both mineral and PTEs) may be lost due to leaching under certain meteorological conditions [48].
The CF ratios found in the present study are within the same intervals of loss/accumulation when compared to CF ratios for lichens exposed in the urban outdoor environment [31], except for Cr and Zn, which presented higher CF ratios in the present study.

3.3. Identification of Emission Sources

3.3.1. Spearman Correlations

Table 2 presents the positive significant Spearman correlations between the elements assessed in the exposed lichens, where values in bold represent the strongest correlations. The possible emission sources which contribute to the atmospheric PM in the study area can be qualitatively identified from the correlation matrix.
Given the significant positive correlations between Al-Si (r = 0.97), Al-Fe (r = 0.87), Al-Ti (r = 0.84), Si-Fe (r = 0.87), Si-Ti (r = 0.84), Ti-Fe (r = 0.84), it is concluded that these elements come from a same source, namely soil, since Al, Si, and Ti are typical soil elements [53,54]. In fact, Al, Fe, Si, and trace elements, such as Rb, are associated with feldspars, quartz, micas and Ti is associated with the titanite silicate [43].
In South European regions, such as Portugal, it is known that atmospheric PM can have a contribution to dust transport episodes from the interior of the Iberian Peninsula and the Saharan desert [55]. It was already found that Fe could reach high concentrations during these episodes of long-transportation of dust [56]. In fact, during the exposure period of the lichens, some Saharan dust events were identified (e.g., from 18 to 21 March 2020) [50].
An important amount of Al, Si, Mg, and Ti was also found in sinter plant emissions (fugitive and emissions coming from the cooling area) [1], consisting of internally mixed aluminosilicates/metallic particles, suggesting also an association between these elements and sinter plant emissions [12,53]. However, the steelworks installed in the study area do not include sintering in their processes, since scrap is used as raw material and not iron ore, and, therefore, this specific source is excluded in the study area.
Strong correlations were found between Fe-Mn (r = 0.77), Fe-Cr (r = 0.81), Cr-Mn (r = 0.92), Cr-Zn (r = 0.73), and Mn-Zn (r = 0.81) in the exposed lichens, which may indicate the existence of a common emission source, assigned to the steelworks in the study area, considering that previous studies also found strong correlations between Fe-Mn-Zn near iron and steel industry [8,57]. Typically, Fe, Cr, and Mn are considered tracers of the iron and steel industry [43]. Therefore, in the present study, despite a natural origin (soil), Fe also has the contribution of this industrial activity. This influence of different sources on Fe was found in other studies [58].
In fact, especially when scrap is employed as a raw material in the steel industry, levels of Cr and Mn can be higher [1]. Zn is also considered a tracer of the steel industry [43,59], being especially related to the use of galvanized scrap [60].
Finally, the association of Pb with the steel industry has been reported by several authors [1,61]. In the present study, significant but weak correlations between Pb-Fe (r = 0.38), Pb-Cr (r = 0.37), Pb-Mn (r = 0.41), and Pb-Zn (r = 0.48) were found, which may suggest some association of Pb emissions with the steelworks, but also the existence of other sources.
A high correlation was found for Fe and Co (r = 0.86), which suggests that Co probably is originated from the steelworks, as already found in previous studies [8].
The Electric Arc Furnace (EAF) steelmaking plants, as the one installed in the study area, typically present fugitive emissions from raw materials and slag handling and storage piles, in addition to the emission from the furnace and rolling [8].
According to the Best Available Techniques (BAT) for Iron and Steel Production [62], emissions from EAF steelmaking plants include the following metals, after abatement: Hg, Pb, Cr, Ni, Zn, Cd, and Cu (where Cr and Ni are, naturally, higher in the production of stainless steel). Typical emissions have concentrations of 0.5–50 mg dust/Nm3. Emissions from the secondary metallurgy (ladle metallurgy, ingot and continuous casting, oxygen blow unit) include the following metals, after abatement: Pb, Co, Ni, Se, Te, Sb, Cr, Cu, Mn, V, Sn [62].
The primary constituents of the slag produced in the EAF process are Ca, Fe, Al, Mg, Mn, Si, and Cr, where calcium and magnesium oxides are added as fluxing agents [63]. In the present study, no correlation was found between Ca-Fe, Ca-Al, Ca-Mg, Ca-Mn, Ca-Si, or Ca-Cr, and, therefore, EAF slag is not identified as a probable source.
A moderate correlation has been found between Cu-Cr (r = 0.54) and Cu-Zn (r = 0.50), which may suggest multiple emission sources. These elements are usually associated with emissions from road traffic, namely through tire (Cr, Zn) and brake wear (Zn, Cu, Ba) and motor oil (Zn) [38,43,64]. Moderate correlations between Cu, Cd, and Zn were also found regarding car traffic [65].
USEPA-SPECIATE and SPECIEUROPE databases were used to identify emission sources, selecting source profiles from the steel industry (general and specific processes) and traffic (break and tire wearing, highway vehicles) [54]. Additionally, a comparison was also done with the elemental ratios determined in the settled dust collected in the study area in 2019 [7]. Figure 6A shows a visible high correlation between Fe and Mn concentrations assessed in the exposed lichens (r2 = 0.80) and a notable association to the steelmaking processes according to the databases of source profiles. It is also verified some similarity between the Fe/Mn ratio in the exposed lichens and the settled dust collected previously in the study area (11.5 and 7.0, respectively).
Figure 6B shows a very high correlation between Cr and Mn measured in lichens (r2 = 0.97) and a notable similarity between the Cr/Mn ratio in the exposed lichens and the settled dust (0.30 and 0.36, respectively), supporting that both elements are associated with the steelworks emissions.
Furthermore, Figure 6C,D evidence the association of Cr, Zn, and Cu to traffic emissions. Concurrently it is verified that Cr/Zn ratios in the exposed lichens and in the settled dust (0.1 and 3.6, respectively) are substantially different, suggesting that Zn emissions in the study area are substantially affected by traffic and not only by the steelworks.

3.3.2. Concentrations versus Distance to Steelworks

The mean elemental concentrations were calculated for the lichens exposed in a circle with center in industry B, which is the biggest steelwork, and with 1 km radius (“(0–1)”), and in consecutive concentric crowns of 1 or 2 km radius, namely, “(1–2)”, “(2–3)”, “(3–4)” and “(4–6)”, where the edges of each interval are the distance (in km) from the established center in industry B. Figure 7 shows the mean concentrations of selected elements regarding the distance from industry B. It should be mentioned that there is a distance of around 2 km from both steelworks in the study area (namely, industries A and B), and thus both facilities are included (and are expected to directly influence the samples) in the 2 km radius (namely the samples located in “(0–1)” and “(1–2)”).
As expected, elements associated with the steel industry (Fe, Cr, Mn, Zn, and Pb) registered the highest concentrations in the lichens exposed near the steelworks, within a 1 or 2 km radius from industry B, and with concentrations decreasing with distance, as shown by Figure 7. Mean Zn levels showed a peak at 4 km from industry B, and this may reflect the traffic influence since it is in the distance range of the nearby highway.
Comparing the mean concentrations of exposed lichens in the 1 km radius circle, to the control sample exposed at a 6580 m distance from industry B, levels decreased more sharply for Cr (6.1 times), than for Zn (2.5), Pb (1.9), Mn (1.7), and Fe (1.5).
The same analysis regarding the variability of EFs with the distance to industry B was also conducted, as shown by Figure 8. The same trend was observed, namely, the EFs of Mn and Fe decrease on a regular basis with the distance; in the case of Zn and Cr, there is an increase in the EF at a distance of 3–4 km from the industry B, consistent with the traffic highway influence.
The ratio between the EFs for samples exposed at a distance less than and more than 1 km from the steelworks’ area is shown in Figure 9. This ratio is higher for Cr (1.77), Zn (1.50), and Pb (1.45), i.e., close to steelworks there was an increase in EF for these elements, as well as for Mn and Fe. Nevertheless, only the EFs of the elements Cr (p = 0.015), Mn (p = 0.035), and Fe (p = 0.031) presented significant differences, between the samples located at less than and more than 1 km from the steelworks.

3.3.3. Spatial Distribution Patterns

As a final step in the identification of emission sources, the spatial distribution of elemental concentrations in the exposed lichens was mapped (Figure 10, Figure 11 and Figure 12). The similarities of the spatial distributions obtained between elements and the correlations found above, support the grouping of elements by types of emission source: natural origin from the soil (Al, Si, and Ti—Figure 10, top); a mixture of soil with natural and anthropogenic origins (Fe and Mg—Figure 10, bottom); anthropogenic resulting from the industrial activity, namely steelworks (Co, Cr, Mn, Pb, and Zn—Figure 11); anthropogenic from road traffic (Cu and Zn—Figure 12, top) and biomass burning (Br and K—Figure 12, bottom).
There is a close resemblance of the spatial distributions of soil origin elements Al, Si, and Ti, as expected from the high correlation coefficients obtained for these elements. The relatively high concentrations of these elements nearby the steelworks is certainly due to the resuspension of particles of geological origin, by the wind and heavy duty traffic [66].
In the studied area, Fe and Mg concentrations may result from the natural emission of soil and the activity of the steelworks, which may provide these elements from fugitive emissions. This hypothesis is supported by the similarity of the spatial distribution of these elements compared to Co, Cr, Mn, Cr, and Zn, associated with the steelworks. Regarding Pb, an additional hot spot can be observed in the north of the study area, but no industry could be identified. Additionally, the spatial distribution of Cu evidences high concentrations near roads due to the influence of traffic [43].
As the wind blew predominantly from the north (N), northwest (NW), and west (W) (Figure 3) during the exposure period, the highest concentrations of the elements related to the steelworks would be expected downwind of this emission source, that is in the south-southeast sector of the sampling grid, as is indeed apparent from the spatial distribution in Figure 11.
The S enrichment in the industrial area may be associated with metal smelting [43] and, in fact, good correlations were obtained above for S-Fe (r = 0.62), S-Cr (r = 0.62), and S-Mn (r = 0.62); but may also derive from local emissions of SOx related to coal/coke burning [43], since there is a coke kiln at a lime company located in the steelworks park-industry C. The significant correlation S-Cu (r = 0.52) also suggests some association with traffic emissions [12], which may be due to heavy-duty traffic in the steelworks area (Figure 12).
Elements Br and K present a spatial distribution with some similarities (and a significant correlation, r = 0.54), indicating a common anthropogenic source, such as biomass burning [43]. This is in accordance with the analysis of the national inventory of emissions for the municipality of Seixal, which identified residential combustion as one of the main emission sources in the region, concerning both PM10 and PM2.5, besides industry and road transport [67].
Figure S2 provides the spatial distribution in the study area of the chemical elements that were not associated with an identified source, such as Ca. Lime and calcium carbide are mineral additions used in the production process of industry B [28], and industry C is a factory of lime products and derivatives, with both industries located in the steelworks park. Nevertheless, high levels of Ca were not obtained in the steelworks area, but in the southeast of the grid (Figure S2). Since there is no cement industry or quarry in the study area, Ca enrichment is probably due to the physiological characteristics of the lichens [32], as granules with calcium oxalate may occur in lichens due to environmental factors [68].
The spatial distribution of the elements discussed previously is in accordance with the Air Quality Map of Seixal, according to which the highest values of particulate matter (PM10 and PM2.5) were registered in the industrial areas and main roads of the municipality [67]. The authors concluded that the air quality index regarding PM10 is “average” generically throughout the municipality, except in the industrial area where it is “weak” or “bad”; for PM2.5, the air quality index varies from generically “weak” to “very bad” in the industrial area.

4. Conclusions

The present study aimed to perform a biomonitoring study of air pollution using lichens in an urban-industrial area where the local population has expressed several concerns regarding air quality in recent years. The methodology applied in this study allowed to involve the local population and to obtain a spatial distribution of the chemical elements absorbed by the lichens during a period where the activities in the area were mainly the industries, due to the COVID-19 national lockdown.
Using a set of different analysis approaches to the chemical content of the exposed transplanted lichens in the study area, it was possible to identify hotspots of air pollution in the area (taking advantage of the possibility to assess the spatial variability due to the use of a grid of biomonitors) and to define the potential sources affecting the local air quality.
A total of five different sources were identified as contributing to the local air quality: (i) a soil source of natural origin (based on Al, Si, and Ti), (ii) a soil source of natural and anthropogenic origins (based on Fe and Mg), (iii) a source from the local industrial activity, namely steelworks (based on Co, Cr, Mn, Pb, and Zn); (iv) a source from the road traffic (based on Cr, Cu, and Zn), and (v) a source of biomass burning (based on Br and K).
The impact of the industries located in the study area on the local air quality was identified (namely, the steelworks). To identify the quantitative contribution of such sources and to design targeted mitigation measures, future efforts should be conducted to assess particulate matter levels by reference methods, which will allow a more comprehensive and quantitative analysis of the issue, taking into account the national and European legislation.
Despite biomonitoring being a technique known for decades, this study showed again the potentialities and advantages of such an approach (used as a complementary strategy to reference methodologies), providing spatial information of air pollution. This information is crucial to define future strategies, based on knowledge: understand exposure levels of the population to air pollution, their potential health impacts, and to define mitigation measures to improve air quality.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijerph19031364/s1, Figure S1: Procedure for preparation of pellets: (A) used ball mill RETSCH, (B) homogeneised powder of lichen sample, (C) used Pelletiser SPECAC, and (D) pellet of a lichen sample. Figure S2. Spatial distribution of elements without a specific identified source. Table S1. Reference values of soil composition (in mg·Kg−1) defined by Mason and Moore [42], used for the calculation of the Enrichment factors.

Author Contributions

L.A., C.A.G. and A.R.J. were responsible for the management of volunteers; L.A., C.A.G., A.R.J., N.C. and S.M.A. were responsible for the methodology of biomonitoring using transplanted lichens from a clean area (regarding exposure and collection of samples); L.A., C.A.G., A.R.J. and I.D. performed the experimental work at the laboratory regarding the treatment of samples and their preparation for further analysis; Z.K. performed the chemical analysis of samples; L.A. and C.A.G. performed the data analysis, including statistical, and wrote the original draft. N.C. and S.M.A. contributed to the review and editing of the original draft. S.M.A. was also responsible for the coordination and management of the project. All authors have read and agreed to the published version of the manuscript.

Funding

N. Canha acknowledges the funding by national funds through FCT-Fundação para a Ciência e Tecnologia, I.P. (Portugal) for his IST-ID contract (IST-ID/098/2018). The FCT support is also acknowledged by C2TN/IST authors (UIDB/04349/2020+UIDP/04349/2020) and by CESAM author (UIDB/50017/2020+UIDP/50017/2020).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated and analysed during the current study are not publicly available but are available from the corresponding author on reasonable request.

Acknowledgments

Acknowledgements are due to Câmara Municipal do Seixal, for funding and all the support given in the sampling process. Instituto Português do Mar e da Atmosfera (IPMA) is also gratefully acknowledged for providing the meteorological data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hleis, D.; Fernández-Olmo, I.; Ledoux, F.; Kfoury, A.; Courcot, L.; Desmonts, T.; Courcot, D. Chemical profile identification of fugitive and confined particle emissions from an integrated iron and steelmaking plant. J. Hazard. Mater. 2013, 250–251, 246–255. [Google Scholar] [CrossRef] [PubMed]
  2. Karagulian, F.; Belis, C.A.; Dora, C.F.C.; Prüss-Ustün, A.M.; Bonjour, S.; Adair-Rohani, H.; Amann, M. Contributions to cities’ ambient particulate matter (PM): A systematic review of local source contributions at global level. Atmos. Environ. 2015, 120, 475–483. [Google Scholar] [CrossRef]
  3. Harrison, R.M.; Yin, J. Particulate matter in the atmosphere: Which particle properties are important for its effects on health? Sci. Total Environ. 2000, 249, 85–101. [Google Scholar] [CrossRef]
  4. Loomis, D.; Grosse, Y.; Lauby-Secretan, B.; El Ghissassi, F.; Bouvard, V.; Benbrahim-Tallaa, L.; Guha, N.; Baan, R.; Mattock, H.; Straif, K. The carcinogenicity of outdoor air pollution. Lancet Oncol. 2013, 14, 1262–1263. [Google Scholar] [CrossRef]
  5. Pope, C.A., III; Dockery, D.W. Health Effects of Fine Particulate Air Pollution: Lines that Connect. J. Air Waste Manag. Assoc. 2012, 2247, 709–742. [Google Scholar] [CrossRef] [PubMed]
  6. Chaíça, I. Poeira negra que cobre Paio Pires não é nociva para a saúde. Origem Permanece Desconhecida. Público, 11 May 2019. Available online: https://www.publico.pt/2019/05/11/local/noticia/poeira-negra-cobre-paio-pires-nao-inalavel-sao-precisos-estudos-saber-onde-vem-1872369(accessed on 9 January 2022).
  7. Justino, A.R.R.; Canha, N.; Gamelas, C.; Coutinho, J.T.T.; Kertesz, Z.; Almeida, S.M.M. Contribution of micro-PIXE to the characterization of settled dust events in an urban area affected by industrial activities. J. Radioanal. Nucl. Chem. 2019, 322, 1953–1964. [Google Scholar] [CrossRef]
  8. Mohiuddin, K.; Strezov, V.; Nelson, P.F.; Stelcer, E. Characterisation of trace metals in atmospheric particles in the vicinity of iron and steelmaking industries in Australia. Atmos. Environ. 2014, 83, 72–79. [Google Scholar] [CrossRef]
  9. Silva, A.V.; Oliveira, C.M.; Canha, N.; Miranda, A.I.; Almeida, S.M. Long-Term Assessment of Air Quality and Identification of Aerosol Sources at Setúbal, Portugal. Int. J. Environ. Res. Public Health 2020, 17, 5447. [Google Scholar] [CrossRef]
  10. Machemer, S.D. Characterization of Airborne and Bulk Particulate from Iron and Steel Manufacturing Facilities. Environ. Sci. Technol. 2004, 38, 381–389. [Google Scholar] [CrossRef]
  11. Querol, X.; Viana, M.; Alastuey, A.; Amato, F.; Moreno, T.; Castillo, S.; Pey, J.; de la Rosa, J.; Sánchez de la Campa, A.; Artíñano, B.; et al. Source origin of trace elements in PM from regional background, urban and industrial sites of Spain. Atmos. Environ. 2007, 41, 7219–7231. [Google Scholar] [CrossRef]
  12. Taiwo, A.M.; Beddows, D.C.S.; Calzolai, G.; Harrison, R.M.; Lucarelli, F.; Nava, S.; Shi, Z.; Valli, G.; Vecchi, R. Receptor modelling of airborne particulate matter in the vicinity of a major steelworks site. Sci. Total Environ. 2014, 490, 488–500. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Almeida, S.M.; Lage, J.; Freitas, M.D.C.; Pedro, A.I.; Ribeiro, T.; Silva, A.V.; Canha, N.; Almeida-Silva, M.; Sitoe, T.; Dionisio, I.; et al. Integration of biomonitoring and instrumental techniques to assess the air quality in an industrial area located in the coastal of central Asturias, Spain. J. Toxicol. Environ. Health Part A Curr. Issues 2012, 75, 1392–1403. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Lage, J.; Almeida, S.M.; Reis, M.A.; Chaves, P.C.; Ribeiro, T.; Garcia, S.; Faria, J.P.; Fernández, B.G.; Wolterbeek, H.T. Levels and Spatial Distribution of Airborne Chemical Elements in a Heavy Industrial Area Located in the North of Spain. J. Toxicol. Environ. Health Part A Curr. Issues 2014, 77, 856–866. [Google Scholar] [CrossRef] [PubMed]
  15. Fränzle, O. Chapter 2 Bioindicators and environmental stress assessment. Trace Met. Other Contam. Environ. 2003, 6, 41–84. [Google Scholar] [CrossRef]
  16. Markert, B.A.; Breure, A.M.; Zechmeister, H.G. Chapter 1 Definitions, strategies and principles for bioindication/biomonitoring of the environment. Trace Met. Other Contam. Environ. 2003, 6, 3–39. [Google Scholar] [CrossRef]
  17. Gür, F.; Yaprak, G. Biomonitoring of metals in the vicinity of Soma coal-fired power plant in western Anatolia, Turkey using the epiphytic lichen, Xanthoria parietina. J. Environ. Sci. Health Part A Toxic/Hazard. Subst. Environ. Eng. 2011, 46, 1503–1511. [Google Scholar] [CrossRef]
  18. Bermejo-Orduna, R.; McBride, J.R.; Shiraishi, K.; Elustondo, D.; Lasheras, E.; Santamaría, J.M. Biomonitoring of traffic-related nitrogen pollution using Letharia vulpina (L.) Hue in the Sierra Nevada, California. Sci. Total Environ. 2014, 490, 205–212. [Google Scholar] [CrossRef]
  19. Canha, N.; Freitas, M.D.C.; Almeida, S.M. Contribution of short irradiation instrumental neutron activation analysis to assess air pollution at indoor and outdoor environments using transplanted lichens. J. Radioanal. Nucl. Chem. 2019, 320, 129–137. [Google Scholar] [CrossRef]
  20. Conti, M.E.; Cecchetti, G. Biological monitoring: Lichens as bioindicators of air pollution assessment—A review. Environ. Pollut. 2001, 114, 471–492. [Google Scholar] [CrossRef]
  21. Garty, J. Biomonitoring atmospheric heavy metals with lichens: Theory and application. CRC. Crit. Rev. Plant Sci. 2001, 20, 309–371. [Google Scholar] [CrossRef]
  22. Leonardo, L.; Damatto, S.R.; Gios, B.R.; Mazzilli, B.P. Lichen specie Canoparmelia texana as bioindicator of environmental impact from the phosphate fertilizer industry of São Paulo, Brazil. J. Radioanal. Nucl. Chem. 2014, 299, 1935–1941. [Google Scholar] [CrossRef]
  23. Bozkurt, Z. Determination of airborne trace elements in an urban area using lichens as biomonitor. Environ. Monit. Assess. 2017, 189, 573. [Google Scholar] [CrossRef] [PubMed]
  24. Cruz, A.M.J.; Freitas, M.D.C.; Canha, N.; Verburg, T.G.; Almeida, S.M.; Wolterbeek, H.T. Spatial mapping of the city of Lisbon using biomonitors. Int. J. Environ. Health 2012, 6, 1. [Google Scholar] [CrossRef]
  25. Bontempi, E.; Bertuzzi, R.; Ferretti, E.; Zucca, M.; Apostoli, P.; Tenini, S.; Depero, L.E. Micro X-ray fluorescence as a potential technique to monitor in-situ air pollution. Microchim. Acta 2008, 161, 301–305. [Google Scholar] [CrossRef]
  26. PORDATA Base de Dados Portugal Contemporâneo. Available online: http://www.pordata.pt/Municipios/Ambiente+de+Consulta/Tabela (accessed on 22 September 2021).
  27. APA—Agência Portuguesa do Ambiente. Licença Ambiental Lusosider Aços Planos, S.A.; APA: Amadora, Portugal, 2008. [Google Scholar]
  28. APA—Agência Portuguesa do Ambiente. Licença Ambiental LA no 658_1.1_2017—SN Seixal_Siderurgia Nacional, S.A.; APA: Amadora, Portugal, 2017. [Google Scholar]
  29. APA—Agência Portuguesa do Ambiente. Licença Ambiental Microlime—Produtos de Cal e Derivados, S.A.; APA: Amadora, Portugal, 2011. [Google Scholar]
  30. Google Earth Pro 2021, version 7.3.4.8248. Available online: https://www.google.com/intl/pt-PT/earth/versions/#earth-pro (accessed on 22 September 2021).
  31. Canha, N.; Almeida, S.M.; Freitas, M.C.; Wolterbeek, H.T. Indoor and outdoor biomonitoring using lichens at urban and rural primary schools. J. Toxicol. Environ. Health Part A Curr. Issues 2014, 77, 900–915. [Google Scholar] [CrossRef]
  32. Canha, N.; Almeida-Silva, M.; Freitas, M.C.; Almeida, S.M.; Wolterbeek, H.T. Lichens as biomonitors at indoor environments of primary schools. J. Radioanal. Nucl. Chem. 2012, 291, 123–128. [Google Scholar] [CrossRef] [Green Version]
  33. Nečemer, M.; Kump, P.; Ščančar, J.; Jaćimović, R.; Simčič, J.; Pelicon, P.; Budnar, M.; Jeran, Z.; Pongrac, P.; Regvar, M.; et al. Application of X-ray fluorescence analytical techniques in phytoremediation and plant biology studies. Spectrochim. Acta Part B At. Spectrosc. 2008, 63, 1240–1247. [Google Scholar] [CrossRef]
  34. Pessanha, S.; Samouco, A.; Adão, R.; Carvalho, M.L.; Santos, J.P.; Amaro, P. Detection limits evaluation of a portable energy dispersive X-ray fluorescence setup using different filter combinations. X-ray Spectrom. 2017, 46, 102–106. [Google Scholar] [CrossRef]
  35. Sitko, R.; Zawisz, B. Quantification in X-Ray Fluorescence Spectrometry. In X-ray Spectroscopy; InTech: Nappanee, IN, USA, 2012; pp. 137–162. [Google Scholar] [CrossRef] [Green Version]
  36. Carvalho, P.M.S.; Pessanha, S.; Machado, J.; Silva, A.L.; Veloso, J.; Casal, D.; Pais, D.; Santos, J.P. Energy dispersive X-ray fluorescence quantitative analysis of biological samples with the external standard method. Spectrochim. Acta Part B At. Spectrosc. 2020, 174, 105991. [Google Scholar] [CrossRef]
  37. Miller, J.N.; Miller, J.C.; Miller, R.D. Statistics and Chemometrics for Analytical Chemistry, 7th ed.; Pearson: London, UK, 2018; ISBN 9781292186719. [Google Scholar]
  38. Belis, C.A.; Karagulian, F.; Larsen, B.R.; Hopke, P.K. Critical review and meta-analysis of ambient particulate matter source apportionment using receptor models in Europe. Atmos. Environ. 2013, 69, 94–108. [Google Scholar] [CrossRef]
  39. Mohd Tahir, N.; Poh, S.C.; Suratman, S.; Ariffin, M.M.; Shazali, N.A.M.; Yunus, K. Determination of trace metals in airborne particulate matter of Kuala Terengganu, Malaysia. Bull. Environ. Contam. Toxicol. 2009, 83, 199–203. [Google Scholar] [CrossRef] [PubMed]
  40. Canha, N.; Freitas, M.C.; Anawar, H.M.; Dionísio, I.; Dung, H.M.; Pinto-Gomes, C.; Bettencourt, A. Characterization and phytoremediation of abandoned contaminated mining area in Portugal by INAA. J. Radioanal. Nucl. Chem. 2010, 286, 577–582. [Google Scholar] [CrossRef] [Green Version]
  41. Pekey, H. The distribution and sources of heavy metals in Izmit Bay surface sediments affected by a polluted stream. Mar. Pollut. Bull. 2006, 52, 1197–1208. [Google Scholar] [CrossRef] [PubMed]
  42. Mason, B.; Moore, C.B. Principles of Geochemistry; Wiley: New York, NY, USA, 1982. [Google Scholar]
  43. Calvo, A.I.; Alves, C.; Castro, A.; Pont, V.; Vicente, A.M.; Fraile, R. Research on aerosol sources and chemical composition: Past, current and emerging issues. Atmos. Res. 2013, 120–121, 1–28. [Google Scholar] [CrossRef]
  44. Almeida-Silva, M.; Canha, N.; Vogado, F.; Baptista, P.C.; Faria, A.V.; Faria, T.; Coutinho, J.T.; Alves, C.; Almeida, S.M. Assessment of particulate matter levels and sources in a street canyon at Loures, Portugal—A case study of the REMEDIO project. Atmos. Pollut. Res. 2020, 11, 1857–1869. [Google Scholar] [CrossRef]
  45. Fernández, J.A.; Carballeira, A. A comparison of indigenous mosses and topsoils for use in monitoring atmospheric heavy metal deposition in Galicia (northwest Spain). Environ. Pollut. 2001, 114, 431–441. [Google Scholar] [CrossRef]
  46. Frati, L.; Brunialti, G.; Loppi, S. Problems related to lichen transplants to monitor trace element deposition in repeated surveys: A case study from central Italy. J. Atmos. Chem. 2005, 52, 221–230. [Google Scholar] [CrossRef]
  47. Watson, D.F.; Philip, G.M. A Refinement of Inverse Distance Weighted Interpolation. Geoprocessing 1985, 2, 315–327. [Google Scholar]
  48. Godinho, R.M.; Verburg, T.G.; Freitas, M.C.; Wolterbeek, H.T. Accumulation of trace elements in the peripheral and central parts of two species of epiphytic lichens transplanted to a polluted site in Portugal. Environ. Pollut. 2009, 157, 102–109. [Google Scholar] [CrossRef]
  49. Pacheco, A.M.G.; Freitas, M.C.; Baptista, M.S.; Vasconcelos, M.T.S.D.; Cabral, J.P. Elemental levels in tree-bark and epiphytic-lichen transplants at a mixed environment in mainland Portugal, and comparisons with an in situ lichen. Environ. Pollut. 2008, 151, 326–333. [Google Scholar] [CrossRef]
  50. Gama, C.; Relvas, H.; Lopes, M.; Monteiro, A. The impact of COVID-19 on air quality levels in Portugal: A way to assess traffic contribution. Environ. Res. 2021, 193. [Google Scholar] [CrossRef] [PubMed]
  51. Liu, F.; Wang, M.; Zheng, M. Effects of COVID-19 lockdown on global air quality and health. Sci. Total Environ. 2021, 755, 142533. [Google Scholar] [CrossRef] [PubMed]
  52. Gamelas, C.; Abecasis, L.; Canha, N.; Almeida, S.M. The Impact of COVID-19 Confinement Measures on the Air Quality in an Urban-Industrial Area of Portugal. Atmosphere 2021, 12, 1097. [Google Scholar] [CrossRef]
  53. Almeida, S.M.M.; Lage, J.; Fernández, B.; Garcia, S.; Reis, M.A.A.; Chaves, P.C.C. Chemical characterization of atmospheric particles and source apportionment in the vicinity of a steelmaking industry. Sci. Total Environ. 2015, 521–522, 411–420. [Google Scholar] [CrossRef] [PubMed]
  54. Lage, J.; Wolterbeek, H.T.; Reis, M.A.; Chaves, P.C.; Garcia, S.; Almeida, S.M. Source apportionment by positive matrix factorization on elemental concentration obtained in PM10 and biomonitors collected in the vicinities of a steelworks. J. Radioanal. Nucl. Chem. 2016, 309, 397–404. [Google Scholar] [CrossRef]
  55. Rodríguez, S.; Querol, X.; Alastuey, A.; Kallos, G.; Kakaliagou, O. Saharan dust contributions to PM10 and TSP levels in Southern and Eastern Spain. Atmos. Environ. 2001, 35, 2433–2447. [Google Scholar] [CrossRef]
  56. Almeida, S.M.; Freitas, M.C.; Pio, C.A. Neutron activation analysis for identification of African mineral dust transport. J. Radioanal. Nucl. Chem. 2008, 276, 161–165. [Google Scholar] [CrossRef]
  57. Dai, Q.L.; Bi, X.H.; Wu, J.H.; Zhang, Y.F.; Wang, J.; Xu, H.; Yao, L.; Jiao, L.; Feng, Y.C. Characterization and source identification of heavy metals in ambient PM10 and PM2.5 in an integrated Iron and Steel industry zone compared with a background site. Aerosol Air Qual. Res. 2015, 15, 875–887. [Google Scholar] [CrossRef] [Green Version]
  58. Mazzei, F.; Lucarelli, F.; Marenco, F.; Nava, S.; Prati, P.; Valli, G.; Vecchi, R. Elemental composition and source apportionment of particulate matter near a steel plant in Genoa (Italy). Nucl. Instrum. Meth. B 2006, 249, 548–551. [Google Scholar] [CrossRef]
  59. Sammut, M.L.; Noack, Y.; Rose, J. Zinc speciation in steel plant atmospheric emissions: A multi-technical approach. J. Geochem. Explor. 2006, 88, 239–242. [Google Scholar] [CrossRef]
  60. Timonen, K.L.; Wiikinkoski, T.; Ruuskanen, A.R. Source contributions to PM2.5 particles in the urban air of a town situated close to a steel works. Atmos. Environ. 2003, 37, 1013–1022. [Google Scholar] [CrossRef]
  61. Tsai, J.H.; Lin, K.H.; Chen, C.Y.; Ding, J.Y.; Choa, C.G.; Chiang, H.L. Chemical constituents in particulate emissions from an integrated iron and steel facility. J. Hazard. Mater. 2007, 147, 111–119. [Google Scholar] [CrossRef]
  62. Remus, R.; Aguado Monsonet, M.; Roudier, S. DSL Best Available Techniques (BAT) Reference Document:for:Iron and Steel Production:Industrial Emissions Directive 2010/75/EU:(Integrated Pollution Prevention and Control); Joint Research Centre: Luxembourg, 2012. [Google Scholar]
  63. Proctor, D.M.; Fehling, K.A.; Shay, E.C.; Wittenborn, J.L.; Green, J.J.; Avent, C.; Bigham, R.D.; Connolly, M.; Lee, B.; Shepker, T.O.; et al. Physical and chemical characteristics of blast furnace, basic oxygen furnace, and electric arc furnace steel industry slags. Environ. Sci. Technol. 2000, 34, 1576–1582. [Google Scholar] [CrossRef]
  64. Almeida-Silva, M.; Canha, N.; Freitas, M.C.; Dung, H.M.; Dionísio, I. Air pollution at an urban traffic tunnel in Lisbon, Portugal-an INAA study. Appl. Radiat. Isot. 2011, 69, 1586–1591. [Google Scholar] [CrossRef] [PubMed]
  65. Zgłobicki, W.; Telecka, M.; Skupiński, S.; Pasierbińska, A.; Kozieł, M. Assessment of heavy metal contamination levels of street dust in the city of Lublin, E Poland. Environ. Earth Sci. 2018, 77, 774. [Google Scholar] [CrossRef] [Green Version]
  66. Hernández-Mena, L.; Murillo-Tovar, M.; Ramírez-Muñíz, M.; Colunga-Urbina, E.; De La Garza-Rodríguez, I.; Saldarriaga-Noreña, H. Enrichment factor and profiles of elemental composition of PM2.5 in the city of Guadalajara, Mexico. Bull. Environ. Contam. Toxicol. 2011, 87, 545–549. [Google Scholar] [CrossRef]
  67. Universidade de Aveiro—Departamento de Ambiente e Ordenamento, Carta de Qualidade do Ar para o Município do Seixal; University of Aveiro: Aveiro, Portugal, 2019.
  68. Osyczka, P.; Boroń, P.; Lenart-Boroń, A.; Rola, K. Modifications in the structure of the lichen Cladonia thallus in the aftermath of habitat contamination and implications for its heavy-metal accumulation capacity. Environ. Sci. Pollut. Res. 2018, 25, 1950–1961. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location of the study area: (left) Framework of the study area (black rectangle) and the clean area (black dot, Montargil) in Portugal mainland; (center) Seixal municipality location; (right) location of industries A, B, C and D within the study area (geographical data obtained from Google Earth Pro [30]).
Figure 1. Location of the study area: (left) Framework of the study area (black rectangle) and the clean area (black dot, Montargil) in Portugal mainland; (center) Seixal municipality location; (right) location of industries A, B, C and D within the study area (geographical data obtained from Google Earth Pro [30]).
Ijerph 19 01364 g001
Figure 2. Exposed transplanted Flavoparmelia caperata lichens in nylon bags (right) that were placed in the study area (left).
Figure 2. Exposed transplanted Flavoparmelia caperata lichens in nylon bags (right) that were placed in the study area (left).
Ijerph 19 01364 g002
Figure 3. Grid of the spatial distribution of the transplanted and retrieved lichens after the exposure period (green dots) in the study area (left), with identification of the industrial area described before (black rectangle) (geographical data obtained from Google Earth Pro [30]), and (right) wind rose of the prevailing winds monitored in the weather station described in Section 2.3 (data obtained from the Portuguese Institute of Sea and Atmosphere, IPMA), during the exposure period.
Figure 3. Grid of the spatial distribution of the transplanted and retrieved lichens after the exposure period (green dots) in the study area (left), with identification of the industrial area described before (black rectangle) (geographical data obtained from Google Earth Pro [30]), and (right) wind rose of the prevailing winds monitored in the weather station described in Section 2.3 (data obtained from the Portuguese Institute of Sea and Atmosphere, IPMA), during the exposure period.
Ijerph 19 01364 g003
Figure 4. Mean EFs (and standard deviation) of the elements assessed in the exposed lichens. The red line represents the threshold of 10, the accepted minimum for the enrichment from a non-crustal source. Regarding the box plot, the square represents the mean, upper and lower times sign (×) shows the maximum and minimum values, the box provides the 25 percentile, the median and the 75% percentile, and the lower and upper minus sign (−) shows the 5 and 95 percentiles, respectively.
Figure 4. Mean EFs (and standard deviation) of the elements assessed in the exposed lichens. The red line represents the threshold of 10, the accepted minimum for the enrichment from a non-crustal source. Regarding the box plot, the square represents the mean, upper and lower times sign (×) shows the maximum and minimum values, the box provides the 25 percentile, the median and the 75% percentile, and the lower and upper minus sign (−) shows the 5 and 95 percentiles, respectively.
Ijerph 19 01364 g004
Figure 5. Average CFs of the elements present in the exposed samples. *–Significant statistical differences between exposed and unexposed elemental contents of lichens by a Mann-Whitney U test, at 0.05 significance level. +—Not significant statistical difference between exposed and unexposed elemental contents of lichens by the One Sample Wilcoxon test, at 0.05 significance level. Regarding the box plot, the square represents the mean, upper and lower times sign (×) shows the maximum and minimum values, the box provides the 25 percentile, the median and the 75% percentile, and the lower and upper minus sign (−) shows the 5 and 95 percentiles, respectively.
Figure 5. Average CFs of the elements present in the exposed samples. *–Significant statistical differences between exposed and unexposed elemental contents of lichens by a Mann-Whitney U test, at 0.05 significance level. +—Not significant statistical difference between exposed and unexposed elemental contents of lichens by the One Sample Wilcoxon test, at 0.05 significance level. Regarding the box plot, the square represents the mean, upper and lower times sign (×) shows the maximum and minimum values, the box provides the 25 percentile, the median and the 75% percentile, and the lower and upper minus sign (−) shows the 5 and 95 percentiles, respectively.
Ijerph 19 01364 g005
Figure 6. Relation between elements in the exposed lichens, in the settled dust and according to USEPA-SPECIATE and SPECIEUROPE profiles for iron and steel industry and traffic: (A) Mn vs. Fe, (B) Mn vs. Cr, (C) Cu vs. Zn, and (D) Zn vs. Cr. The red line presents the linear regression regarding the data of the exposed lichens.
Figure 6. Relation between elements in the exposed lichens, in the settled dust and according to USEPA-SPECIATE and SPECIEUROPE profiles for iron and steel industry and traffic: (A) Mn vs. Fe, (B) Mn vs. Cr, (C) Cu vs. Zn, and (D) Zn vs. Cr. The red line presents the linear regression regarding the data of the exposed lichens.
Ijerph 19 01364 g006
Figure 7. Mean elemental concentration versus distance to the steelworks.
Figure 7. Mean elemental concentration versus distance to the steelworks.
Ijerph 19 01364 g007
Figure 8. Mean Enrichment Factors (EFs) versus distance to the steelworks. The blue line represents the threshold of 10, the accepted minimum for the enrichment from a non-crustal source.
Figure 8. Mean Enrichment Factors (EFs) versus distance to the steelworks. The blue line represents the threshold of 10, the accepted minimum for the enrichment from a non-crustal source.
Ijerph 19 01364 g008
Figure 9. Ratio between EF in samples located at less than and more than 1 km from steelworks. * Significant statistical difference between the EF in the two groups by a Mann-Whitney U test, at 0.05 significance level.
Figure 9. Ratio between EF in samples located at less than and more than 1 km from steelworks. * Significant statistical difference between the EF in the two groups by a Mann-Whitney U test, at 0.05 significance level.
Ijerph 19 01364 g009
Figure 10. Spatial distribution of soil-related elements in the study area: (top) natural origin and (bottom) a mixture of natural and anthropogenic origins.
Figure 10. Spatial distribution of soil-related elements in the study area: (top) natural origin and (bottom) a mixture of natural and anthropogenic origins.
Ijerph 19 01364 g010
Figure 11. Spatial distribution of elements related with anthropogenic origin, namely industry.
Figure 11. Spatial distribution of elements related with anthropogenic origin, namely industry.
Ijerph 19 01364 g011
Figure 12. Spatial distribution of elements related to an anthropogenic origin, namely from road traffic (top) and biomass burning (bottom).
Figure 12. Spatial distribution of elements related to an anthropogenic origin, namely from road traffic (top) and biomass burning (bottom).
Ijerph 19 01364 g012
Table 1. Mean elemental mass fractions (in mg·Kg−1, dry weight) in unexposed and exposed lichens samples and comparison with values from other studies.
Table 1. Mean elemental mass fractions (in mg·Kg−1, dry weight) in unexposed and exposed lichens samples and comparison with values from other studies.
Present Study(Godinho et al., 2009) [48](Pacheco et al., 2008) [49]
Unexposed LichensExposed LichensIndustrial Area, Sines 1Industrial Area, Sines 1
ElementMean ± SDMinMaxMean ± SDMinMaxMean ± SDMean ± SD
Al2310 ± 190214025502590 ± 450177037901040 ± 1463170
As58.5 ± 24.836.088.030.0 ± 28.411250.33 ± 0.111.33
Br192 ± 17170211180 ± 1614322114.6 ± 1.813.7
Ca208,000 ± 23,000180,000232,000226,000 ± 15,000187,000257,0008890 ± 16904390
Co47.0 ± 10.532.056.054.9 ± 16.2241300.68 ± 0.161.33
Cr7.33 ± 0.587.008.0052.1 ± 90.836076.4 ± 0.47.1
Cu140 ± 16122158163 ± 90102654n.d.n.d.
Fe9860 ± 1480843011,70011,800 ± 3800650030,600651 ± 1373190
K8910 ± 690820097107900 ± 1120441011,7503877 ± 7374170
Mg1130 ± 90102012201230 ± 2708302670992 ± 791320
Mn274 ± 33246314406 ± 297225230324.0 ± 0.528.2
Pb101 ± 1582117130 ± 11181967n.d.n.d.
Rb69.8 ± 17.048.086.069.7 ± 15.1451424.4 ± 1.323.4
S1390 ± 120128015601580 ± 41010703950n.d.n.d.
Se77.5 ± 6.058.086.074.7 ± 5.75886n.d.0.49
Si5460 ± 370501057906240 ± 1280381011,220n.d.n.d.
Sr553 ± 36508586565 ± 74415737n.d.n.d.
Ti1010 ± 10091011401150 ± 2306901840102 ± 1274
Zn347 ± 16333365608 ± 4542802270140 ± 34134
Zr83.8 ± 17.159.098.078.5 ± 19.743156n.d.n.d.
1 after 4-month exposure; SD—standard deviation; n.d.—not determined.
Table 2. Spearman correlations between elements in the exposed lichens, significant at the 0.05 level. Values in bold represent strong correlations (correlation coefficient above 0.70), while weak and medium correlations refer to correlation coefficients between 0.30–0.49 and 0.50–0.70, respectively.
Table 2. Spearman correlations between elements in the exposed lichens, significant at the 0.05 level. Values in bold represent strong correlations (correlation coefficient above 0.70), while weak and medium correlations refer to correlation coefficients between 0.30–0.49 and 0.50–0.70, respectively.
ElementsAlAsBrCaCoCrCuFeKMgMnPbRbSSeSiTiZnZr
Al 0.33 0.710.62 0.870.290.500.56 0.660.56 0.970.840.420.67
As 0.52
Br 0.430.320.540.420.25 0.460.490.350.390.37 0.38
Ca
Co 0.640.270.860.280.410.590.340.620.49 0.720.640.440.60
Cr 0.540.81 0.390.920.370.320.62 0.620.640.730.42
Cu 0.39 0.350.570.32 0.52 0.320.370.500.34
Fe 0.510.770.380.670.62 0.870.840.600.65
K 0.44 0.610.36 0.350.33 0.34
Mg 0.44 0.460.48 0.480.500.290.57
Mn 0.410.330.62 0.580.620.810.43
Pb 0.27 0.48
Rb 0.48 0.660.730.270.58
S 0.610.610.580.41
Si 0.840.420.68
Ti 0.550.66
Zn 0.29
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Abecasis, L.; Gamelas, C.A.; Justino, A.R.; Dionísio, I.; Canha, N.; Kertesz, Z.; Almeida, S.M. Spatial Distribution of Air Pollution, Hotspots and Sources in an Urban-Industrial Area in the Lisbon Metropolitan Area, Portugal—A Biomonitoring Approach. Int. J. Environ. Res. Public Health 2022, 19, 1364. https://doi.org/10.3390/ijerph19031364

AMA Style

Abecasis L, Gamelas CA, Justino AR, Dionísio I, Canha N, Kertesz Z, Almeida SM. Spatial Distribution of Air Pollution, Hotspots and Sources in an Urban-Industrial Area in the Lisbon Metropolitan Area, Portugal—A Biomonitoring Approach. International Journal of Environmental Research and Public Health. 2022; 19(3):1364. https://doi.org/10.3390/ijerph19031364

Chicago/Turabian Style

Abecasis, Leonor, Carla A. Gamelas, Ana Rita Justino, Isabel Dionísio, Nuno Canha, Zsofia Kertesz, and Susana Marta Almeida. 2022. "Spatial Distribution of Air Pollution, Hotspots and Sources in an Urban-Industrial Area in the Lisbon Metropolitan Area, Portugal—A Biomonitoring Approach" International Journal of Environmental Research and Public Health 19, no. 3: 1364. https://doi.org/10.3390/ijerph19031364

APA Style

Abecasis, L., Gamelas, C. A., Justino, A. R., Dionísio, I., Canha, N., Kertesz, Z., & Almeida, S. M. (2022). Spatial Distribution of Air Pollution, Hotspots and Sources in an Urban-Industrial Area in the Lisbon Metropolitan Area, Portugal—A Biomonitoring Approach. International Journal of Environmental Research and Public Health, 19(3), 1364. https://doi.org/10.3390/ijerph19031364

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