Air quality assessment in different urban areas from Rio Grande do Sul state , Brazil , using lichen transplants

Based on the need to monitor air pollution and on the importance of doing it with biological organisms, the present work had as main objective to assess air quality in urban areas through lichen transplants. Besides, we also tested whether there were differences among the studied sites, aiming to point out the main pollution characteristics of the areas. The monitoring was carried out in seven cities, which varied from mainly rural to mainly industrial, in Rio Grande do Sul, southern Brazil. The foliose lichen Parmotrema tinctorum was used as the biomonitor species. Physiological tests and contents of sulphur and heavy metals were measured in the lichen thalli. Other pollutants assessed through air samplers were available for some cities, and modeled concentrations of fine particulate matter and nitrogen oxides for all sites. The effect of air pollution negatively reflected on the physiology parameters of the lichen and the concentration of all heavy metals increased after exposure, especially in the industrial sites. We generated a profile of air pollutants and pointed out that Charqueadas and Esteio had the worst air quality. Since they have high demographic densities, the results are of great public health concern and could be guidelines for future studies.


INTRODUCTION
Air pollution from industrial biomass and fossil fuel burning can spread far beyond the limits of the country that generated it (Akimoto 2003).Globally spread, these emissions are changing the world climate.According to the latest report of the IPCC (Intergovernmental Panel on Climate Change) (Pachauri et al. 2014) recent anthropogenic emissions of greenhouse gases are the highest in history and together with other anthropogenic drivers (such as heavy metals or organic compounds), are extremely likely to have been the dominant cause of the observed warming.In addition to climate change, environmental pollution is also one of the greatest concerns in public health (Cansaran-Duman et al. 2010).Recently the World Health Organization (WHO) called attention to the great impact of air pollution on children's health and mortality worldwide (WHO 2017).
It is very important to adopt local actions to control and monitor air quality, not only at the regional scale of population health, but also considering the global impact of air pollution.The use of biological monitors such as trees (Domingos et al. 2015, Alves et al. 2001), grasses (Crittenden andRead 1979, Sandrin et al. 2013), mosses (Bignal et al. 2008, Gerdol et al. 2014) or lichens (Carreras and Pignata 2002, Loppi and Frati 2006, Riddell et al. 2012) have proven to be efficient.Here we will focus on lichens, which are long considered good tools to monitor air pollution and air quality (Hawksworth et al. 2005).
The use of lichens for monitoring atmospheric pollution is known since independent observations made in England, Munich and Paris in the 1800s documented that these organisms were disappearing from urban areas (Nash III 2008).So far, many studies had shown and proved that lichens are indeed good biomonitors and bioindicators of air pollution of different kinds and origins (Conti and Cecchetti 2001, Käffer et al. 2012, Pinho et al. 2011).
One of the possible ways to monitor air pollution with lichens is the active monitoring method.It consists on the exposure of a welldefined biomonitor species under relatively controlled conditions (Szczepaniak and Biziuk 2003).This strategy allows an immediate and accessible evaluation, which is convenient where there is no other way to measure the atmospheric pollution (Aras et al. 2011), as is the case with most developing countries.
Because of the constant need to monitor air pollution and the importance of using biological organisms, our work had as main objective to assess air quality in urban areas with industrial and rural influences.We looked at changes in lichen biology: 1) pollution absorbed by the biomonitor species; 2) physiological damages observed on the thalli; and used other biomonitoring tools: 3) automatic/ semi-automatic air samplers of pollutants and 4) pollution data models for the areas.We also compared cities to characterize the main pollutants of each site.
These cities have great importance for the state economy hosting many industries (Esteio, Charqueadas and Montenegro), a Petrochemical Complex (Triunfo), and an important agricultural area of family farming (Santo Antônio da Patrulha, Maquiné and Caraá).A description of the main industrial activities, the demographic density and the rural practices of the monitored cities are detailed in Koch et al. (2016).They also have a similar climate, elevation, and some of them (Charqueadas, Esteio and Santo Antônio da Patrulha) have air-monitoring systems with automatic and semi-automatic air samplers.The elevation of sampling areas ranges from 10 to 75 m above sea level and the climate is subtropical humid, Cfa type according to the updated Köppen-Geiger classification (Peel et al. 2007).

SAMPLE PREPARATION AND LICHEN EXPOSURE TREATMENTS
The foliose epiphytic lichen Parmotrema tinctorum (Despr.ex Nyl.) Hale was chosen as the biomonitor species.This species has a known tolerance to air pollution, having been used in other monitoring studies in southern Brazil (Käffer et al. 2012) and is abundant throughout the country (Spielmann and Marcelli 2009).The lichen P. tinctorum is naturally found in all the studied cities, but we chose to transplant the thalli from a less polluted site to standardize the exposure and due to its low abundance in some sites.Samples of P. tinctorum were collected mainly on trunks of planted Eucalyptus spp. in Tapes (Figure 1) to be later transplanted into potentially polluted sites.The municipality of Tapes is located at least 80 km distant from large urban centers and thus far from important sources of air pollution, being agriculture and cattle raising the main economic sources of the region.Before being exposed, the thalli of P. tinctorum were left in the laboratory for acclimation during ten days (Käffer et al. 2012).
Lichens were placed in 55 x 35 cm synthetic nettings in standardized structures ca.120 cm from the ground and shaded to reduce sunlight exposure (Figure 2).They were placed vertically, which was the same position they were found on the trees.When available, lichen transplants were set close to existing automatic/semi-automatic air samplers.Before exposing the samples of P. tinctorum, we performed chemical control analyses to test for metal and sulphur contents.Prior physiological damages were also assessed in some samples of these thalli.Transplants were exposed from March to October 2013 and from April to November 2014.Sub-samples of these materials were taken after two, five and seven months of exposure, in all seven cities monitored, in both years.We gathered data for 14 months and six sub-samples in each city.Maquiné and Caraá were the control sites, based on the characteristics of low urbanization and the absence of large industries in these municipalities.Maquiné was only monitored in 2014, after Caraá, supposed to have lower air pollution, showed some important pollutants concentration in 2013.

PHYSIOLOGICAL ANALYSES
Each of the six lichen sub-samples were tested for physiological responses to atmospheric pollution, and analyzed prior to exposure, for control.All pieces of bark were removed and the lichen thalli were cleaned with a brush.Part of the thalli was macerated and the percentages of live, dead and plasmolysed algal cells were counted under the microscope using the neutral-red test at 5% (Le Blanc 1971, Calvelo andLiberatore 2004).With the same macerate samples we analyzed Chlorophyll a and b according to the method proposed by Boonpragob (2002), which consisted of test tubes with 10 mL of 96% ethanol and 5 cm 2 of lichens.Absorbance in the wavelengths of 649 and 665 nm was measured with spectrophotometer Digimed DME-21.Chlorophyll contents and the percentage of live, dead and plasmolysed algal cells are methods for analyzing physiological damages on the lichen photobiont.Organic carbon content, which measures physiological damage on the lichen mycobiont, was analyzed through wet combustion using the Walkley-Black method with external heat (APHA 2005), on the Laboratory of Soil Analyses at Universidade Federal do Rio Grande do Sul.

CHEMICAL ANALYSES
We tested total content of the following pollutants in the lichen thalli: sulphur (S), copper (Cu), zinc (Zn), iron (Fe), manganese (Mn), chromium (Cr), nickel (Ni), lead (Pb), vanadium (V) and aluminum (Al).Analyses were made through inductively coupled plasma optical emission spectrometry Roessler), but in some cases data were not available for all the monitoring period and did not include all parameters (PM 10 , SO 2, O 3 , CO and NO 2 ).Sample data were compared with the standards of air quality in Brazil (CONAMA 1990) and with the guidelines of World Health Organization (WHO 1999(WHO , 2006)).
We also assessed modeled atmospheric seven cities) differed according to air quality, we performed a PCA (Principal Component Analysis) using a matrix of the sites described by metals, sulphur and modeled pollutants.The result scores were plotted in a biplot and variables scores were scaled by the square root of eigenvalues.This analysis was made with the "rda" function of CRAN software R using package Vegan (R Core Team 2015).

LICHEN PHYSIOLOGICAL RESPONSES
Overall, monitored sites had only slight differences regarding physiological damages on the thalli of Parmotrema tinctorum.The organic carbon contents (OC) after exposure did not show any clear differences among cities, and neither did the chlorophyll contents (Table I) I).The percentage of dead algal cells in 2013 had a similar, but opposite, pattern to the live cells.However, the same did not happen in 2014, when the highest values in seven months were found in Esteio (ES) and the lowest in TR, the latter showing also an expressive higher value of plasmolysed algal cells in 2014 (59.8%) (Table I).Most air pollutants analyzed with the lichen biomonitor, the air samplers and the modeled data had different concentrations in industrial and rural cities (Table II, Figures 3 and 4).
Regarding the metals absorbed by the lichen P. tinctorum, differences were found especially in the concentration of copper, zinc, iron, chromium and nickel, with higher values in CH and MO in general.There was no clear difference among sites for the other heavy metals and sulphur contents.However, when comparing these values with the ones before exposure, in the control sample, we observed that almost all pollutants were more concentrated after the exposure period (Figure 3).
High contents of pollutants were even found in the rural areas (CA, MA and SA): suplhur, zinc, chromium, nickel, lead, vanadium and aluminum had similar or higher concentrations than those found in some industrial sites (Figure 3).It is important to note that chromium had a considerable peak in CA and SA in 2013 (more than 160 mg/ kg), which did not repeat in 2014.Moreover, sulphur peaked at the beginning of 2014, reaching all sampled cities.At the industrial sites (MO, TR, CH and ES) all pollutants increased after exposure when compared to the control (Figure 3).
Data from the air samplers were not available for all monitored cities during the seven months when the lichens were exposed (2013 and 2014).However, based on the available data we found differences among some cities (Table II).The modeled contents of NO x and PM 2.5 (Figure 4) also varied among cities, being higher in ES and CH, both industrial cities.
CH had higher values of PM 10 than SA, similar to the modeled PM 2.5 , where even higher values were found for ES (Figure 4).Particulate matter in CH also exceeded the recommended values for health security in both periods evaluated (WHO 2006)      (SA).NO x also appeared in higher values in ES and CH (Figure 4).We observed greater values of ozone in CH in both years when compared to Esteio (ES), while carbon monoxide did not show any difference.

MAIN POLLUTION COMPOSITION IN EACH CITY
Cities were grouped based on pollution characteristics through a principal component analysis (Figure 4).Axis 1 explained 43.4% of the observed variance and axis 2, 28.8%.Caraá (CA) was more related to higher contents of chromium, while Charqueadas (CH), is closely related to many different metals.A greater relation of Mn, Cu and modeled NO X and PM 2.5 was found in Esteio (ES).Zn was more important in Triunfo (TR) and Montenegro (MO).

DISCUSSION
The present study provides a description of the atmosphere quality at sites with industrial and rural influences.Further, it consists of a pioneer study on air pollution biomonitoring with lichens in most of the sampled cities, showing relevant concentrations of some air pollutants with negative impact on public health and environmental quality.

PHYSIOLOGICAL RESPONSES
The effect of exposure to pollutants on Parmotrema tinctorum was visible in most of the thalli.The overall pattern was an increase in the physiology parameters, with the exception of the percentages of live algal cells.Distinctively, the sites of Montenegro and Maquiné had the highest *Air quality standards (annual means, except when mentioned differently) according to the Brazilian legislation (CONAMA 1990) and the World Health Organization guidelines (WHO 1999(WHO , 2006)).**Only one month available.Abbreviation: m/d = missing data.
Values shown are the mean of the three exposure periods (two, five and seven months after exposure) and the standard deviations of the same data.
percentages of live algal cells after the exposure period, while Triunfo had the lowest values.It is also noteworthy that Triunfo had higher values of dead and plasmolysed algal cells.While some works have noted that the percentage of dead algal cells is related to air pollution (Calvelo and Liberatore 2004, Käff er et al. 2012), here we show that live algal cells can also be a useful parameter.Many authors have demonstrated that lichens are highly sensitive to contents of particulate matter, PAHs (Polycyclic Aromatic Carbons), NO 2 and SO 2 , showing physiological damages when exposed to these pollutants (van Herk et al. 2003, Riddell et al. 2012, Käffer et al. 2012).We did not evaluate organic compounds, however their presence is likely refl ecting in loss of lichen vitality in the industrial cities, namely PAHs, dioxins, furans and others, which tend to bioaccumulate as a result of atmospheric deposition (Augusto et al. 2013).Although we did not measure organic pollution, we hypothesize that important damage on lichen thalli in Triunfo could be due to the eff ect of these pollutants from the petrochemical complex nearby the study area, since most metal concentrations were not higher than in other industrial areas.

AIR POLLUTION MONITORING
Almost all the pollutants absorbed by P. tinctorum increased after exposure, a pattern also found in other studies (Käffer et al. 2012).Lichens are known to accumulate and retain many heavy metals in quantities that exceed their physiological requirements because they are able to keep metals Zinc contents were higher in sites with industrial activities such as Esteio and Montenegro.Esteio has a high demographic density (2918 inhabitants/km 2 , one of the highest in the state of Rio Grande do Sul), which results in an intense daily traffic.According to Minganti et al. (2003) the main potential sources of zinc are indeed industrial emissions and traffic flow.Sulphur peaks were registered at the beginning of 2014 in all sampled sites.These peaks are unlikely to be related to local pollution events, but instead to some global emissions as maybe a volcanic eruption (Thomas and Prata 2011) or another important source.
The observed peak of chromium in the monitored sites of Caraá and Santo Antônio da Patrulha recorded in 2013 could be either related to leather waste burning in any of the small shoe factories in these cities or could have been released (accidentally or not) by some metalworking industry, since it only occurred in one month during the monitored period.In the urban environment of Porto Alegre, the capital and the largest and most populated city of Rio Grande do Sul, Käffer et al. (2012) found 19.0 mg/kg as the mean highest concentration of chromium in P. tinctorum, while in our study that amount was more than three times higher (Caraá = 73.0mg/kg).Charqueadas also had higher contents of chromium (averaging 17.4 mg/kg), but constant throughout the monitoring period.In Europe, a review of chromium averages in mosses between 1990 and 2000 showed that no country had higher concentrations than 8.5 mg/kg during this period and the trend expected for the next years was a reduction in chromium emissions (Harmens et al. 2007).
Despite only having data for PM 10 , SO 2 , NO 2 , O 3 and CO at some sites, this information was valuable to understand the differences among them.It is important to note that all these pollutants were below the levels required by the National Patterns of Air Quality (CONAMA 1990).However, if we consider the guidelines of the World Health Organization for particulate matter (PM 10 ), Charqueadas exceeded the suggested limit for health security, which is an annual mean of 20 μg/ m 3 (WHO 2006) and Santo Antônio da Patrulha is also very close to reaching the limit.Charqueadas has a greater demographic density and more intense industrial activity, while Santo Antônio da Patrulha has not many industries, but lately has had a large increment in motor vehicles (82% increase in seven years) (Silva et al. 2015).
Particulate matter may include a wide range of chemical substances, such as elemental and organic carbon compounds; oxides of silicon, aluminum and iron; trace-metals; sulfates; nitrates and ammonia (Bergamaschi et al. 2007).The great danger of high values of PM 10 is due to its association to increased rates of bronchitis, reduced lung functions (WHO 2000), and also an increase of cancer risk, since many compounds of the particulate matter are also mutagenic (Coronas et al. 2008(Coronas et al. , 2009)).
In a recent study, Coronas et al. (2008) analyzed the mutagenic activity of airborne particulate material (PM 10 ) in Triunfo through the Salmonella/ microsome assay, and found positive genotoxic activity, which they related to the presence of PAHs (Polycyclic Aromatic Hydrocarbons).The concentration of PM 10 found by this study in 2008 was even higher than was registered here (up to 100 μg/m 3 ).In Esteio, Coronas et al. (2009) found mutagenic activity in all samples of particulate matter (PM 10 ), using the same method.Similar results were found by Pereira et al. (2013) for Montenegro and Santo Antônio da Patrulha.In the latter city, Silva et al. (2015) recorded peaks of 64 μg/m 3 of PM 2.5 , much higher than the limit defined by WHO (25 μg/m 3 for this parameter), as well as high values of PAHs.
Almost all other greenhouse gases differed among sites, with the exception of Carbon Monoxide.This gas is highly related to vehicular emissions (Mayer 1999) and has been controlled by the Brazilian government directly at car industries, through several actions.Since 2003, when this program started, two thirds of emissions have been reduced in the country (IBAMA 2014).As we have observed in this study, it is indeed true that emissions of CO are controlled, since their concentrations did not change between sites and are much lower than the WHO recommendations (WHO 1999).
As for the other gases, high contents were found in Charqueadas (SO 2 and O 3 ) and Esteio (NO 2 ), which is comparable to the results found through modeled data of NO x .The higher concentration of SO 2 in Charqueadas might be due to the coalfueled power plant, a well-known source of SO 2 , heavy metals and other elements (Garty 2000).The higher concentration of ozone could be a result of pollutants coming from elsewhere, since this greenhouse gas arises from the oxidation of nitrogen oxides, carbon monoxides and/or volatile organic compounds (Jenkin and Clemitshaw 2000).Nitrogen dioxide (NO 2 ), which had the highest values in Esteio, derives mainly from the burning of fossil fuels, and so it is closely related to motor vehicle traffic (Mayer 1999).Esteio has indeed one of the highest demographic densities in our study (2918 inhabitants/km 2 ), almost eighteen times the density of Charqueadas (163 inhabitants/km 2 ) (IBGE 2014).

MAIN POLLUTION COMPOSITION IN EACH CITY
We created an air pollution profile for the seven monitored cities using data for metals, sulphur, modeled PM 2.5 and NO x concentrations.It showed the most polluted sites to be Charqueadas and Esteio, with high contents of many pollutants as already discussed.The former has a thermoelectric power plant; the latter, important petrochemical industries; in both cities there are large steel manufacturing industries.These, together with other kinds of industries and also traffic emissions, are probably the main air pollution sources in these cities, which had also high concentrations of gases (SO 2 , NO 2 and O 3 ) and PM 10 .
The cities of Triunfo and Montenegro, mainly industrial, also accounted for the observed air pollution.Triunfo holds a large Petrochemical Complex, while Montenegro is situated at the main direction of winds coming from this complex.However, organic compounds that were not evaluated would probably better separate these sites from the rural ones (SA, MA and CA).Caraá (CA) seemed to be highly influenced by the great peak of chromium in 2013, while Santo Antônio da Patrulha (SA) and Maquiné (MA) were characterized by lower concentrations of the analyzed pollutants.

CONCLUSIONS
As lichens are one of the most valuable biomonitors of atmospheric pollution and monitoring this pollution is extremely important for public and environmental health, this study contributes to the knowledge of the air quality in southern Brazil.
Based on the profile of air pollution we generated for the sites we could point out that Charqueadas and Esteio had the worst air quality in the region.Since both cities have high demographic densities, these results are of great public health concern.Moreover, it is important to highlight that air pollution is not limited to dense populated and industrial cities.We showed here that even apparently unpolluted sites, such as small rural cities with family farming as the main economic activity, could be affected by external (or internal) sources of pollutants.
Although it is still hard to determine the exact sources of elements reaching the lichen thalli, the knowledge on the distribution patterns of these pollutants is indeed of high value for society and for decision making such as where to direct attention to health or where the air quality is already NATÁLIA MOSSMANN KOCH et al. too saturated.Thus, this study could be a tool for environmental health protection and a guideline for future studies.

Figure 2
Figure 2 -a.Thalli of Parmotrema tinctorum exposed in the standardized structures.b.Detail of P. tinctorum in the synthetic nettings after three months of exposure.
concentrations of NO X and PM 2.5 , from the Chemistry Coupled Aerosol-Tracer Transport Modeling System (CCATT-BRAMS, Freitas et al. 2009, Longo et al. 2013), for each study site.Data were provided by the Group of Modeling of the Atmosphere and its Interfaces from Universidade Federal de Pelotas, state of Rio Grande do Sul (http://ccatt.ufpel.edu.br).Because the model had been recently implemented on the state, data were only available for 2014 on.STATISTICAL ANALYSES Aiming to test whether monitored sites (i.e.
. Regarding the percentage of algal live cells, Triunfo (TR) was the site with the lowest values at the end of the exposure period in 2013 and 2014 (29.8% and 10.4%, respectively), while Montenegro (MO) showed the highest values after seven months of exposure in 2013 (46.1%).The highest values for Maquiné (MA) and Charqueadas (CH) were in 2014 (both around 46%).This parameter also showed a decrease in almost all sites after exposure, with the exception of MA and MO in 2014 (Table = Organic Carbon; LC = Live Cells; DC = Dead Cells; PC = Plasmolysed Cells; R = Rural; I = Industrial; TA = Tapes; b.e.= before exposure; MA = Maquiné; CA = Caraá; SA = Santo Antônio da Patrulha; MO = Montenegro; CH = Charqueadas; ES = Esteio; TR = Triunfo; nd = not detected; T0 = contents prior to exposure; T1 = after two months; T2 = after five months; T3 = after seven months

Figure 3 -
Figure 3 -Contents of sulphur and metals absorbed by the lichen Parmotrema tinctorum during seven months of exposure in 2013 and 2014.The red column represents the basal values, prior to exposure in the cities.T0 = contents prior to exposure; T1 = after two months; T2 = after fi ve months; T3 = after seven months.

Figure 3 (
Figure 3 (continuation)-Contents of sulphur and metals absorbed by the lichen Parmotrema tinctorum during seven months of exposure in 2013 and 2014.The red column represents the basal values, prior to exposure in the cities.T0 = contents prior to exposure; T1 = after two months; T2 = after fi ve months; T3 = after seven months.

TABLE I Lichen physiological parameters measured to evaluate the damages on the photobiont (Chl a, Chl b and percentage of live, plasmolysed and dead cells) and on the mycobiont (percentage of organic carbon) after exposure. Thalli of Parmotrema tinctorum were exposed during seven months in 2013 and 2014.
(Table II).Measured contents of SO 2 , showed differences among CH, ES and SA, with Charqueadas (CH) showing higher concentrations than Esteio (ES) and Santo Antônio da Patrulha *Abbreviations: OC