Tools for determining critical levels of atmospheric ammonia under the influence of multiple disturbances

Lichens as Bioindicators of Global Change Drivers

In recent decades, the scientific community has put the spotlight on the severe impacts that environmental stressors are producing on ecosystem functioning worldwide [...].

NH3 concentrations below the current critical level affect the epiphytic macrolichen communities - Evidence from a Northern European City

Acidophytic, oligotrophic lichens on tree trunks are widely considered to be the most sensitive biota to elevated concentrations of atmospheric ammonia (NH₃). We studied the relationships between measured NH₃ concentrations and the composition of macrolichen communities on the acidic bark of Pinus sylvestris and Quercus robur and the base-rich bark of Acer platanoides and Ulmus glabra at ten roadside and ten non-roadside sites in Helsinki, Finland. NH₃ and nitrogen dioxide (NO₂) concentrations were higher at the roadside than non-roadside sites indicating traffic as the main source of NH₃ and nitrogen oxides (NO_x). The diversity of oligotrophs on Quercus was lower at the roadside than non-roadside sites, while that of eutrophs was higher. The abundance and presence of oligotrophic acidophytes (e.g., Hypogymnia physodes) decreased with increasing NH₃ concentration (2-year means = 0.15-1.03 μg m^-3) especially on Q. robur, while those of eutrophic/nitrophilous species (e.g., Melanohalea exasperatula, Physcia tenella) increased. The abundance of some nitrophytes seemed to depend only on bark pH, i.e., their abundances were highest on Ulmus, which had the highest average bark pH. Overall, the results of lichen bioindicator studies may depend on tree species (bark pH) and lichen species used in calculating indices describing the air quality impact. Nevertheless, Quercus is recommended to be used to study the impact of NH₃ alone and in combination with NO_x on lichen communities, because the responses of both oligotrophic acidophytes and eutrophic species can already be observed at NH₃ concentrations below the current critical level.

Synthesis of lichen response to gaseous nitrogen: ammonia versus nitrogen dioxide

The dominant chemical form of nitrogen pollution in the atmosphere in the U.S. is shifting from oxidized nitrogen, primarily from combustion of fossil fuels, to reduced nitrogen from agricultural animal waste and fertilizer applications. Does it matter to lichens? In this synthesis, we characterize U.S. air concentrations of the most ubiquitous gaseous forms of reduced and oxidized nitrogen, NO2 and NH3, respectively, and their direct effects on lichens. In the U.S., the 3-year average (2017-2019) of the annual mean for each monitoring site ranges up to 56.4 μg NO2 m-3 (~30 ppb) and 6 μg NH3 m-3 (~9 ppb). The spatial coverage of current routine monitoring of NO2 and NH3 likely does not accurately represent exposures of NO2 to ecosystems in rural areas or capture spikes of NH3 concentrations proximal to intensive agriculture, which are documented to exceed 700 μg NH3 m-3 (~1000 ppb) for short durations. Both NO2 and NH3 can act as nutrients to lichens, but as exposures rise, both can cause physiological stress and mortality that then change community composition and diversity. There is a growing body of evidence that lichen community composition is altered at current levels of exposure in the U.S. with estimated no effect or lowest effect concentrations from <1-3 μg m-3 NO2 and <1 μg m-3 NH3. Better spatial characterization of both NO2 and NH3 concentrations, especially near intensive agriculture, would help to characterize the extent of the impacts across the U.S. These findings are discussed in the context of U.S. air pollution policy.

Resilience of Epiphytic Lichens to Combined Effects of Increasing Nitrogen and Solar Radiation

Lichens are classified into different functional groups depending on their ecological and physiological response to a given environmental stressor. However, knowledge on lichen response to the synergistic effect of multiple environmental factors is extremely scarce, although vital to get a comprehensive understanding of the effects of global change. We exposed six lichen species belonging to different functional groups to the combined effects of two nitrogen (N) doses and direct sunlight involving both high temperatures and ultraviolet (UV) radiation for 58 days. Irrespective of their functional group, all species showed a homogenous response to N with cumulative, detrimental effects and an inability to recover following sunlight, UV exposure. Moreover, solar radiation made a tolerant species more prone to N pollution’s effects. Our results draw attention to the combined effects of global change and other environmental drivers on canopy defoliation and tree death, with consequences for the protection of ecosystems.

How much matching there is in functional, phylogenetic and taxonomic optima of epiphytic macrolichen communities along a European climatic gradient?

Adopting an integrative approach that explicitly includes the different facets of biodiversity is crucial to assess the response of biological communities to changing environments. The identification of the optimal climatic conditions where communities maximize their functional, phylogenetic and taxonomic diversity is useful to compare whether the optima of the different facets of biodiversity match. Using a wide climatic gradient across Europe, we quantified the functional, phylogenetic and taxonomic diversity of epiphytic macrolichen communities, which are valuable early-warning ecological indicators. We ordinated 22 environmental variables and simultaneously illustrated non-parametric regressions of the diversity metrics against the climatic space using the 'hilltop plot' method to detect the climatic conditions in which the different diversity facets peaked and to compare the match between them. Functional diversity predicted at least part of the peaks of phylogenetic and taxonomic diversity, but phylogenetic and taxonomic hotspots did not overlap. Epiphytic macrolichen communities maximized their functional and phylogenetic diversity in the southernmost forests, with the Mediterranean region appearing as a biodiversity hotspot. Regarding the studied traits, photobiont type and growth form showed clearly defined optima while the quantitative physiological traits and families' optima did not show this pattern in response to climate. The different facets of biodiversity were not surrogates of each other highlighting the need for an integrative approach to assess the effect of environmental changes on communities and to establish conservation priorities. As functional traits mediated the response of lichen communities to climate, preserving high functional diversity might indirectly preserve high phylogenetic and taxonomic diversity. Relevant ecological indicators useful to develop rapid assessment methods to evaluate the effects of climatic changes include the photobiont type and growth form. The lack of relation between quantitative traits and climate call for further research to unveil their role as ecological indicators of small-scale variables or as effect traits.

Modeling the provision of air-quality regulation ecosystem service provided by urban green spaces using lichens as ecological indicators

The UN Sustainable Development Goals states that urban air pollution must be tackled to create more inclusive, safe, resilient and sustainable cities. Urban green infrastructures can mitigate air pollution, but a crucial step to use this knowledge into urban management is to quantify how much air-quality regulation can green spaces provide and to understand how the provision of this ecosystem service is affected by other environmental factors. Considering the insufficient number of air quality monitoring stations in cities to monitor the wide range of natural and anthropic sources of pollution with high spatial resolution, ecological indicators of air quality are an alternative cost-effective tool. The aim of this work was to model the supply of air-quality regulation based on urban green spaces characteristics and other environmental factors. For that, we sampled lichen diversity in the centroids of 42 urban green spaces in Lisbon, Portugal. Species richness was the best biodiversity metric responding to air pollution, considering its simplicity and its significative response to the air pollutants concentration data measured in the existent air quality monitoring stations. Using that metric, we then created a model to estimate the supply of air quality regulation provided by green spaces in all green spaces of Lisbon based on the response to the following environmental drivers: the urban green spaces size and its vegetation density. We also used the unexplained variance of this model to map the background air pollution. Overall, we suggest that management should target the smallest urban green spaces by increasing green space size or tree density. The use of ecological indicators, very flexible in space, allow the understanding and the modeling of the provision of air-quality regulation by urban green spaces, and how urban green spaces can be managed to improve air quality and thus improve human well-being and cities resilience.

Do Different Teams Produce Different Results in Long-Term Lichen Biomonitoring?

Lichen biomonitoring programs focus on temporal variations in epiphytic lichen communities in relation to the effects of atmospheric pollution. As repeated surveys are planned at medium to long term intervals, the alternation of different operators is often possible. This involves the need to consider the effect of non-sampling errors (e.g., observer errors). Here we relate the trends of lichen communities in repeated surveys with the contribution of different teams of specialists involved in sampling. For this reason, lichen diversity data collected in Italy within several ongoing biomonitoring programs have been considered. The variations of components of gamma diversity between the surveys have been related to the composition of the teams of operators. As a major result, the composition of the teams significantly affected data comparability: Similarity (S), Species Replacement (R), and Richness Difference (D) showed significant differences between “same” and “partially” versus “different” teams, with characteristics trends over time. The results suggest a more careful interpretation of temporal variations in biomonitoring studies.

δ15N of lichens reflects the isotopic signature of ammonia source

Although it is generally accepted that δ¹⁵N in lichen reflects predominating N isotope sources in the environment, confirmation of the direct correlation between lichen δ¹⁵N and atmospheric δ¹⁵N is still missing, especially under field conditions with most confounding factors controlled. To fill this gap and investigate the response of lichens with different tolerance to atmospheric N deposition, thalli of the sensitive Evernia prunastri and the tolerant Xanthoria parietina were exposed for ten weeks to different forms and doses of N in a field manipulation experiment where confounding factors were minimized. During this period, several parameters, namely total N, δ¹⁵N and chlorophyll a fluorescence, were measured. Under the experimental conditions, δ¹⁵N in lichens quantitatively responded to the δ¹⁵N of released gaseous ammonia (NH₃). Although a high correlation between the isotopic signatures in lichen tissue and supplied N was found both in tolerant and sensitive species, chlorophyll a fluorescence indicated that the sensitive species very soon lost its photosynthetic functionality with increasing N availability. The most damaging response to the different N chemical forms was observed with dry deposition of NH₃, although wet deposition of ammonium ions had a significant observable physiological impact. Conversely, there was no significant effect of nitrate ions on chlorophyll a fluorescence, implying differential sensitivity to dry deposition versus wet deposition and to ammonium versus nitrate in wet deposition. Evernia prunastri was most sensitive to NH₃, then NH₄⁺, with lowest sensitivity to NO₃^-. Moreover, these results confirm that lichen δ¹⁵N can be used to indicate the δ¹⁵N of atmospheric ammonia, providing a suitable tool for the interpretation of the spatial distribution of NH₃ sources in relation to their δ¹⁵N signal.

Mapping ammonia risk on sensitive habitats in Ireland

The aim of this study was to provide a simple, cost-effective, risk-based map of terrestrial areas in Ireland where environmental quality may be at risk from atmospheric ammonia. This risk-based approach identifies Natura 2000 sites in Ireland at risk from agricultural atmospheric ammonia, collating best available data using Geographical Information Systems (GIS). In mapping ammonia risk on sensitive habitats (MARSH), the method identifies sources of ammonia, classifying them on a scale of risk from 0 to 5. These sources are subsequently summed based on a weighting determined by their contribution to national emissions divided by their potentially impacted area. A Pearson's correlation coefficient of 0.72 allows for concentrations from United Kingdom's FRAME modelling to be applied to the MARSH model, which are corrected based on recent monitoring. Applying Designation Weighted Indicators (DWI), the MARSH model predicts that 80.7, 34.3 and 5.9% of Natura 2000 sites in Ireland may exceed ambient concentrations of 1, 2, and 3 μg/m³, respectively. A Nitroindex map of Ireland based on available lichen records was also developed and is presented as part of this study. This Nitroindex was used to identify areas where impacts have already been recorded, thus informing the classification of sites "at-risk". The combination of both the MARSH and Nitroindex models ascertains which Natura 2000 sites are most at risk, thereby providing valuable data to relevant authorities. The MARSH model acts as a first step towards screening and assessing Natura 2000 sites most at risk from atmospheric ammonia, providing a tool to demonstrate compliance with the National Emissions Ceilings Directive.

Modelling local uncertainty in relations between birth weight and air quality within an urban area: combining geographically weighted regression with geostatistical simulation

In this study, we combine known methods to present a new approach to assess local distributions of estimated parameters measuring associations between air quality and birth weight in the urban area of Sines (Portugal). To model exposure and capture short-distance variations in air quality, we use a Regression Kriging estimator combining air quality point data with land use auxiliary data. To assess uncertainty of exposure, the Kriging estimator is incorporated in a sequential Gaussian simulation algorithm (sGs) providing a set of simulated exposure maps with similar spatial structural dependence and statistical properties of observed data. Following the completion of the simulation runs, we fit a geographically weighted generalized linear model (GWGLM) for each mother's place of residence, using observed health data and simulated exposure data, and repeat this procedure for each simulated map. Once the fit of GWGLM with all exposure maps is finished, we take the distribution of local estimated parameters measuring associations between exposure and birth weight, thus providing a measure of uncertainty in the local estimates. Results reveal that the distribution of local parameters did not vary substantially. Combining both methods (GWGLM and sGs), however, we are able to incorporate local uncertainty on the estimated associations providing an additional tool for analysis of the impacts of place in health.

Mapping Portuguese Natura 2000 sites in risk of biodiversity change caused by atmospheric nitrogen pollution

In this paper, we assess and map the risk that atmospheric nitrogen (atN) pollution poses to biodiversity in Natura 2000 sites in mainland Portugal. We first review the ecological impacts of atN pollution on terrestrial ecosystems, focusing on the biodiversity of Natura 2000 sites. These nature protection sites, especially those located within the Mediterranean Basin, are under-characterized regarding the risk posed by atN pollution. We focus on ammonia (NH3) because this N form is mostly associated with agriculture, which co-occurs at or in the immediate vicinity of most areas of conservation interest in Portugal. We produce a risk map integrating NH3 emissions and the susceptibility of Natura 2000 sites to atN pollution, ranking habitat sensitivity to atN pollution using expert knowledge from a panel of Portuguese ecological and habitat experts. Peats, mires, bogs, and similar acidic and oligotrophic habitats within Natura 2000 sites (most located in the northern mountains) were assessed to have the highest relative risk of biodiversity change due to atN pollution, whereas Natura 2000 sites in the Atlantic and Mediterranean climate zone (coastal, tidal, and scrubland habitats) were deemed the least sensitive. Overall, results allowed us to rank all Natura 2000 sites in mainland Portugal in order of evaluated risk posed by atN pollution. The approach is of great relevance for stakeholders in different countries to help prioritize site protection and to define research priorities. This is especially relevant in countries with a lack of expertise to assess the impacts of nitrogen on biodiversity and can represent an important step up from current knowledge in such countries.

Eco-physiological response of Hypnum cupressiforme Hedw. to increased atmospheric ammonia concentrations in a forest agrosystem

Ammonia (NH₃) emissions are linked to eutrophication, plant toxicity and ecosystem shifts from N to P limitation. Bryophytes are key components of terrestrial ecosystems, yet highly sensitive to N deposition. Hence, physiological responses of mosses may be indicative of NH₃-related impacts, and thus useful to foresee future ecosystem damages and establish atmospheric Critical Levels (CLEs). In this work, samples of Hypnum cupressiforme Hedw. were seasonally collected along a well-defined NH₃ concentration gradient in an oak woodland during a one-year period. We performed a comprehensive evaluation of tissue chemistry, stoichiometry, metabolic enzymes, antioxidant response, membrane damages, photosynthetic pigments, soluble protein content and N and C isotopic fractionation. Our results showed that all the physiological parameters studied (except P, K, Ca and C) responded to the NH₃ gradient in predictable ways, although the magnitude and significance of the response were dependent on the sampling season, especially for enzymatic activities and pigments content. Nutritional imbalances, membrane damages and disturbance of cellular C and N metabolism were found as a consequence to NH₃ exposure, being more affected the mosses more exposed to the barn atmosphere. These findings suggested significant implications of intensive farming for the correct functioning of oak woodlands and highlighted the importance of seasonal dynamics in the study of key physiological processes related to photosynthesis, mosses nutrition and responses to oxidative stress. Finally, tissue N showed the greatest potential for the identification of NH₃-related ecological end points (estimated CLE=3.5μgm^-3), whereas highly scattered physiological responses, although highly sensitive, were not suitable to that end.

Deriving nitrogen critical levels and loads based on the responses of acidophytic lichen communities on boreal urban Pinus sylvestris trunks

The deposition of reactive nitrogen (N) compounds currently predominates over sulphur (S) deposition in most of the cities in Europe and North America. Acidophytic lichens growing on tree trunks are known to be sensitive to both N and S deposition. Given that tree species and climatic factors affect the composition of epiphytic lichen communities and modify lichen responses to air pollution, this study focused on the impact of urban air pollution on acidophytes growing on boreal conifer trunks. The study was performed in the Helsinki metropolitan area, southern Finland, where annual mean nitrogen dioxide (NO₂) concentrations range from 4-5μgm^-3 to >50μgm^-3. In addition, background forest sites in southern and northern Finland were included. The results demonstrated elevated N contents (≥0.7%) in Hypogymnia physodes and Platismatia glauca at all the sites where the species occurred. In the Helsinki metropolitan area, a higher frequency of green algae+Scoliociosporum chlorococcum and reduced numerical frequencies of other indicator lichen species (e.g. Pseudevernia furfuracea, Bryoria spp., Usnea spp.) were associated with elevated atmospheric concentrations of NO₂ and particulate matter containing N, as well as elevated concentrations of inorganic N in bark. The N isotope values (δ¹⁵N) of lichens supported the uptake of oxidized N mainly originating from road traffic. Sulphur dioxide (SO₂) also negatively affected the most sensitive species, despite the current low levels (1-4μgm^-3yr^-1). Critical levels of 5μgNO₂m^-3yr^-1 and 0.5μgNH₃m^-3yr^-1, and a critical load of 2-3kgNha^-1yr^-1 are proposed for protecting the diversity of boreal acidophytes. This study calls for measurements of the throughfall of various N fractions in urban forest ecosystems along precipitation and temperature gradients to verify the proposed critical levels and loads.

The influence of growth form and substrate on lichen ecophysiological responses along an aridity gradient

In this paper, we investigated whether growth form and substrate in lichens influence their physiological responses along an aridity gradient. Thalli of the foliose lichen Parmotrema perlatum and the fruticose lichen Ramalina canariensis were transplanted in selected rural/forested sites of Southern Portugal characterized by a different aridity index. Physiological parameters including photosynthetic performances, assimilation pigments, ergosterol content and sample viability were measured prior to exposure (winter) and after 6-month exposure (summer). Photosynthetic performances were also investigated in common native foliose and fruticose epiphytic lichens and in fruticose terricolous species. Both transplanted and native lichens showed a decrease in photosynthetic performances in summer and lower performances in sites classified as drier and higher performances in humid forested sites. No relevant differences occurred between epiphytic foliose and fruticose growth forms. However, terricolous fruticose samples showed a significant difference in humid and drier sites and between winter and summer, probably due to microclimatic conditions similarly to other biological crusts.

Using nitrogen concentration and isotopic composition in lichens to spatially assess the relative contribution of atmospheric nitrogen sources in complex landscapes

Reactive nitrogen (Nr) is an important driver of global change, causing alterations in ecosystem biodiversity and functionality. Environmental assessments require monitoring the emission and deposition of both the amount and types of Nr. This is especially important in heterogeneous landscapes, as different land-cover types emit particular forms of Nr to the atmosphere, which can impact ecosystems distinctively. Such assessments require high spatial resolution maps that also integrate temporal variations, and can only be feasibly achieved by using ecological indicators. Our aim was to rank land-cover types according to the amount and form of emitted atmospheric Nr in a complex landscape with multiple sources of N. To do so, we measured and mapped nitrogen concentration and isotopic composition in lichen thalli, which we then related to land-cover data. Results suggested that, at the landscape scale, intensive agriculture and urban areas were the most important sources of Nr to the atmosphere. Additionally, the ocean greatly influences Nr in land, by providing air with low Nr concentration and a unique isotopic composition. These results have important consequences for managing air pollution at the regional level, as they provide critical information for modeling Nr emission and deposition across regional as well as continental scales.

Ecological impacts of atmospheric pollution and interactions with climate change in terrestrial ecosystems of the Mediterranean Basin: Current research and future directions

Mediterranean Basin ecosystems, their unique biodiversity, and the key services they provide are currently at risk due to air pollution and climate change, yet only a limited number of isolated and geographically-restricted studies have addressed this topic, often with contrasting results. Particularities of air pollution in this region include high O₃ levels due to high air temperatures and solar radiation, the stability of air masses, and dominance of dry over wet nitrogen deposition. Moreover, the unique abiotic and biotic factors (e.g., climate, vegetation type, relevance of Saharan dust inputs) modulating the response of Mediterranean ecosystems at various spatiotemporal scales make it difficult to understand, and thus predict, the consequences of human activities that cause air pollution in the Mediterranean Basin. Therefore, there is an urgent need to implement coordinated research and experimental platforms along with wider environmental monitoring networks in the region. In particular, a robust deposition monitoring network in conjunction with modelling estimates is crucial, possibly including a set of common biomonitors (ideally cryptogams, an important component of the Mediterranean vegetation), to help refine pollutant deposition maps. Additionally, increased attention must be paid to functional diversity measures in future air pollution and climate change studies to establish the necessary link between biodiversity and the provision of ecosystem services in Mediterranean ecosystems. Through a coordinated effort, the Mediterranean scientific community can fill the above-mentioned gaps and reach a greater understanding of the mechanisms underlying the combined effects of air pollution and climate change in the Mediterranean Basin.

Geostatistical uncertainty of assessing air quality using high-spatial-resolution lichen data: A health study in the urban area of Sines, Portugal

In most studies correlating health outcomes with air pollution, personal exposure assignments are based on measurements collected at air-quality monitoring stations not coinciding with health data locations. In such cases, interpolators are needed to predict air quality in unsampled locations and to assign personal exposures. Moreover, a measure of the spatial uncertainty of exposures should be incorporated, especially in urban areas where concentrations vary at short distances due to changes in land use and pollution intensity. These studies are limited by the lack of literature comparing exposure uncertainty derived from distinct spatial interpolators. Here, we addressed these issues with two interpolation methods: regression Kriging (RK) and ordinary Kriging (OK). These methods were used to generate air-quality simulations with a geostatistical algorithm. For each method, the geostatistical uncertainty was drawn from generalized linear model (GLM) analysis. We analyzed the association between air quality and birth weight. Personal health data (n=227) and exposure data were collected in Sines (Portugal) during 2007-2010. Because air-quality monitoring stations in the city do not offer high-spatial-resolution measurements (n=1), we used lichen data as an ecological indicator of air quality (n=83). We found no significant difference in the fit of GLMs with any of the geostatistical methods. With RK, however, the models tended to fit better more often and worse less often. Moreover, the geostatistical uncertainty results showed a marginally higher mean and precision with RK. Combined with lichen data and land-use data of high spatial resolution, RK is a more effective geostatistical method for relating health outcomes with air quality in urban areas. This is particularly important in small cities, which generally do not have expensive air-quality monitoring stations with high spatial resolution. Further, alternative ways of linking human activities with their environment are needed to improve human well-being.

Evaluating green infrastructure in urban environments using a multi-taxa and functional diversity approach

Forested areas within cities host a large number of species, responsible for many ecosystem services in urban areas. The biodiversity in these areas is influenced by human disturbances such as atmospheric pollution and urban heat island effect. To ameliorate the effects of these factors, an increase in urban green areas is often considered sufficient. However, this approach assumes that all types of green cover have the same importance for species. Our aim was to show that not all forested green areas are equal in importance for species, but that based on a multi-taxa and functional diversity approach it is possible to value green infrastructure in urban environments. After evaluating the diversity of lichens, butterflies and other-arthropods, birds and mammals in 31 Mediterranean urban forests in south-west Europe (Almada, Portugal), bird and lichen functional groups responsive to urbanization were found. A community shift (tolerant species replacing sensitive ones) along the urbanization gradient was found, and this must be considered when using these groups as indicators of the effect of urbanization. Bird and lichen functional groups were then analyzed together with the characteristics of the forests and their surroundings. Our results showed that, contrary to previous assumptions, vegetation density and more importantly the amount of urban areas around the forest (matrix), are more important for biodiversity than forest quantity alone. This indicated that not all types of forested green areas have the same importance for biodiversity. An index of forest functional diversity was then calculated for all sampled forests of the area. This could help decision-makers to improve the management of urban green infrastructures with the goal of increasing functionality and ultimately ecosystem services in urban areas.