Chemical composition of rainwater under two events of aerosol transport: A Saharan dust outbreak and wildfires

A one-year campaign of joint sampling of aerosols and precipitation, carried out in León, Spain, allowed to study the impact of two special events that affected the air quality in the north of the country, on rainfall in the city: a period with wildfires and a Saharan dust intrusion. The wildfires that occurred in northern Portugal and northwestern Spain in August 2016 affected the chemistry of rainfall on 15 August 2016, causing an increase in concentrations of NH₄⁺, Na⁺, Cl^-, K⁺, Mg²⁺, Ca²⁺, SO₄^2- and NO₃^- and in the concentrations of organic acids, which was reflected in the levels of soluble and insoluble organic carbon. This led to acidification of rainwater (pH = 4.8). The second precipitation event was registered between 11 and 14 February 2017, during which the rainwater was collected in four daily fractions (P₁, P₂, P₃ and P₄). The rain sample of 12 February (P₂) coincided with a Saharan dust intrusion that reached northern Iberia that day. The chemical composition of P₂ showed an increase in the Ca²⁺ (>800%), Mg²⁺ (71%), Cl^- (62%), and SO₄^2- (33%) concentrations, with respect to P₁. The input of crustal elements to the atmosphere helped to neutralize the P₂ rainwater, causing pH values higher than 6.5. Once the dust intrusion left the north of the Peninsula, the composition of rainwater P₃ and P₄ revealed a mixture of marine contribution with local anthropogenic emissions, as well as a decrease in ion concentrations and conductivity, and an increase in pH values.

Stream chemistry response to changing nitrogen and sulfur deposition in two mountain areas in the Iberian Peninsula

Sulfur (S) and nitrogen (N) increasing anthropogenic emissions in the last century has arisen wide concern on the ecological effects of S and N deposition. In this paper, we use bulk deposition and stream water measurements in the central Pyrenees (PYR-C and PYR-AT sites) and Montseny (MSY-TM0) covering different time lengths in the period 1983-2017 to investigate how these mountain environments respond to ongoing changes of regional emissions to the atmosphere. PYR-C, in spite of its position far away from urban and industrial areas, presented higher SO₄-S, NO₃-N and NH₄-N bulk deposition than the Montseny site closer to Barcelona and the inclusion of dry deposition only reversed this pattern for NO₃-N. This indicates that distance to pollution sources does not protect these mountain sites from a considerable impact of pollution. Time-trends in SO₄^2- and NO₃^- concentrations in bulk deposition were similar between sites: SO₄^2- monotonically decreased, while NO₃^- increased until the mid-2000 s and decreased thereafter. In the period 1983 to 2017, SO₂ emissions in Europe (EU-28) decreased by 95%, while in the SO₄^2- concentrations in bulk deposition declined by 35-50% in Pyrenees and Montseny respectively and SO₄^2- concentrations in the streams by 25-35%, respectively. Other sources of SO₄^2- (e.g. episodic African dust) may explain the different reduction rate between anthropogenic emissions and bulk deposition. Net S budget was positive for MSY-TM0 (indicating flushing from the catchment) and negative for the PYR-C site (indicating retention), while it was close to zero for the other Pyrenean site, but in the PYR-C site net retention showed a significant increasing trend tending to lower retention in recent years. Bulk N deposition in the Pyrenees was lower but stream concentrations and export was higher than at Montseny, this leading to less N retention in the Pyrenean sites. However, the MSY-TM0 site showed a trend towards less N retention in recent years. This was driven by higher exports during the wet months, which would correspond to a first stage of N saturation according Stoddard's classification.

Denitrification rates in lake sediments of mountains affected by high atmospheric nitrogen deposition

During the last decades, atmospheric nitrogen loading in mountain ranges of the Northern Hemisphere has increased substantially, resulting in high nitrate concentrations in many lakes. Yet, how increased nitrogen has affected denitrification, a key process for nitrogen removal, is poorly understood. We measured actual and potential (nitrate and carbon amended) denitrification rates in sediments of several lake types and habitats in the Pyrenees during the ice-free season. Actual denitrification rates ranged from 0 to 9 μmol N₂O m⁻² h⁻¹ (mean, 1.5 ± 1.6 SD), whereas potential rates were about 10-times higher. The highest actual rates occurred in warmer sediments with more nitrate available in the overlying water. Consequently, littoral habitats showed, on average, 3-fold higher rates than the deep zone. The highest denitrification potentials were found in more productive lakes located at relatively low altitude and small catchments, with warmer sediments, high relative abundance of denitrification nitrite reductase genes, and sulphate-rich waters. We conclude that increased nitrogen deposition has resulted in elevated denitrification rates, but not sufficiently to compensate for the atmospheric nitrogen loading in most of the highly oligotrophic lakes. However, there is potential for high rates, especially in the more productive lakes and landscape features largely govern this.

Towards a better understanding of long-term wood-chemistry variations in old-growth forests: A case study on ancient Pinus uncinata trees from the Pyrenees

Dendrochemical studies in old forests are still underdeveloped. Old trees growing in remote high-elevation areas far from direct human influence constitute a promising biological proxy for the long-term reconstructions of environmental changes using tree-rings. Furthermore, centennial-long chronologies of multi-elemental chemistry at inter- and intra-annual resolution are scarce. Here, we use a novel non-destructive method by applying Micro X-ray fluorescence (μXRF) to wood samples of old Pinus uncinata trees from two Pyrenean high-elevation forests growing on acidic and basic soils. To disentangle ontogenetic (changes in tree age and diameter) from environmental influences (e.g., climate warming) we compared element patterns in sapwood (SW) and heartwood (HW) during the pre-industrial (1700-1849) and industrial (1850-2008) periods. We quantified tree-ring growth, wood density and relative element concentrations at annual (TRW, tree-ring) to seasonal resolution (EW, earlywood; LW, latewood) and related them to climate variables (temperature and precipitation) and volcanic eruptions in the 18th and 19th centuries. We detected differences for most studied elements between SW and HW along the stem and also between EW and LW within rings. Long-term positive and negative trends were observed for Ca and K, respectively. Cl, P and S showed positive trends during the industrial period. However, differences between sites were also notable. Higher values of Mg, Al, Si and the Ca/Mn ratio were observed at the site with acidic soil. Growing-season temperatures were positively related to growth, maximum wood density and to the concentration of most elements. Peaks in S, Fe, Cl, Zn and Ca were linked to major volcanic eruptions (e.g., Tambora in 1815). Our results reveal the potential of long-term wood-chemistry studies based on the μXRF non-destructive technique to reconstruct environmental changes.

Nitrogen-cycling genes in epilithic biofilms of oligotrophic high-altitude lakes (central Pyrenees, Spain)

Microbial biofilms in oligotrophic environments are the most reactive component of the ecosystem. In high-altitude lakes, exposed bedrock, boulders, gravel, and sand in contact with highly oxygenated water and where a very thin epilithic biofilm develops usually dominate the littoral zone. Traditionally, these surfaces have been considered unsuitable for denitrification, but recent investigations have shown higher biological diversity than expected, including diverse anaerobic microorganisms. In this study, we explored the presence of microbial N-cycling nirS and nirK (denitrification through the conversion of NO2(-) to NO), nifH (N2 fixation), anammox (anaerobic ammonium oxidation), and amoA (aerobic ammonia oxidation, both bacterial and archaeal) genes in epilithic biofilms of a set of high-altitude oligotrophic lakes in the Pyrenees. The concentrations of denitrifying genes determined by quantitative PCR were two orders of magnitude higher than those of ammonia-oxidizing genes. Both types of genes were significantly correlated, suggesting a potential tight coupling nitrification-denitrification in these biofilms that deserves further confirmation. The nifH gene was detected after nested PCR, and no signal was detected for the anammox-specific genes used. The taxonomic composition of denitrifying and nitrogen-fixing genes was further explored by cloning and sequencing. Interestingly, both microbial functional groups were richer and more genetically diverse than expected. The nirK gene, mostly related to Alphaproteobacteria (Bradyrhizobiaceae), dominated the denitrifying gene pool as expected for oxygen-exposed habitats, whereas Deltaproteobacteria (Geobacter like) and Cyanobacteria were the most abundant among nitrogen fixers. Overall, these results suggest an epilithic community more metabolically diverse than previously thought and with the potential to carry out an active role in the biogeochemical nitrogen cycling of high-altitude ecosystems. Measurements of activity rates should be however carried out to substantiate and further explore these findings.

High bacterial diversity in epilithic biofilms of oligotrophic mountain lakes

Benthic microbial biofilms attached to rocks (epilithic) are major sites of carbon cycling and can dominate ecosystem primary production in oligotrophic lakes. We studied the bacterial community composition of littoral epilithic biofilms in five connected oligotrophic high mountain lakes located at different altitudes by genetic fingerprinting and clone libraries of the 16S rRNA gene. Different intra-lake samples were analyzed, and consistent changes in community structure (chlorophyll a and organic matter contents, and bacterial community composition) were observed along the altitudinal gradient, particularly related with the location of the lake above or below the treeline. Epilithic biofilm genetic fingerprints were both more diverse among lakes than within lakes and significantly different between montane (below the tree line) and alpine lakes (above the tree line). The genetic richness in the epilithic biofilm was much higher than in the plankton of the same lacustrine area studied in previous works, with significantly idiosyncratic phylogenetic composition (specifically distinct from lake plankton or mountain soils). Data suggest the coexistence of aerobic, anaerobic, phototrophic, and chemotrophic microorganisms in the biofilm, Bacteroidetes and Cyanobacteria being the most important bacterial taxa, followed by Alpha-, Beta-, Gamma-, and Deltaproteobacteria, Chlorobi, Planctomycetes, and Verrucomicrobia. The degree of novelty was especially high for epilithic Bacteroidetes, and up to 50 % of the sequences formed monophyletic clusters distantly related to any previously reported sequence. More than 35 % of the total sequences matched at <95 % identity to any previously reported 16S rRNA gene, indicating that alpine epilithic biofilms are unexplored habitats that contain a substantial degree of novelty within a short geographical distance. Further research is needed to determine whether these communities are involved in more biogeochemical pathways than previously thought.

Altitudinal distributions of BDE-209 and other polybromodiphenyl ethers in high mountain lakes

The present study shows the occurrence of 2,2',3,3',4,4',5,5',6,6'-decabromodiphenyl ether (BDE-209) in microbial biofilms of Pyrenean and Tatra high mountain lakes despite its low vapor pressure and high hydrophobicity. Aerosol air transport is therefore a feasible mechanism for BDE-209 accumulation in sites up to 2688 m above sea level. This compound and other PBDEs exhibit altitudinally-dependent distribution involving higher concentrations with increasing mountain lake elevation. However, the apparently very high enthalpies of the concentration gradients observed, including BDE-209, suggest that bacterial anaerobic debromination also plays a significant role in the resulting altitudinal distributions. This microbial mechanism explains the relative abundances of PBDEs and their within lake differences between rocky and sediment microbial biofilms, thereby showing that the altitudinal pattern observed is not purely due to water temperature control on bacterial activity but also to changes in the availability of anaerobic microenvironments which increase with increasing lake productivity at lower altitudes.

A method for upscaling soil parameters for use in a dynamic modelling assessment of water quality in the Pyrenees

Dynamic modelling of hydrochemistry is a valuable tool to study and predict the recovery of surface waters from acidification, and to assess the effects of confounding factors (such as delayed soil response and changing climate) that cause hysteresis during reversal from acidification. The availability of soil data is often a limitation for the regional application of dynamic models. Here we present a method to upscale site-specific soil properties to a regional scale in order to circumvent that problem. The method proposed for upscaling relied on multiple regression models between soil properties and a suite of environmental variables used as predictors. Soil measurements were made during a field survey in 13 catchments in the Pyrenees (NW Spain). The environmental variables were derived from mapped or remotely sensed topographic, lithological, land-cover, and climatic information. Regression models were then used to model soil parameters, which were supplied as input for the biogeochemical model MAGIC (Model for Acidification of Groundwater In Catchments) in order to reconstruct the history of acidification in Pyrenean lakes and forecast the recovery under a scenario of reduced acid deposition. The resulting simulations were then compared with model runs using field measurements as input parameters. These comparisons showed that regional averages for the key water and soil chemistry variables were suitably reproduced when using the modelled parameters. Simulations of water chemistry at the catchment scale also showed good results, whereas simulated soil parameters reflected uncertainty in the initial modelled estimates.

The chemistry of precipitation, throughfall and stemflow in two holm oak (Quercus ilex L.) forests under a contrasted pollution environment in NE Spain

Atmospheric deposition was studied through measurements of bulk deposition, throughfall and stemflow at two sites of contrasted exposure to pollution in the Montseny mountains (Northeastern Spain). To explore the contribution of local or distant sources at both sites, concentration data and precipitation amounts (log transformed) for both bulk deposition and net throughfall were fitted by linear regression. These models indicated the more important contribution of washout scavenging processes and dry deposition at the pollution exposed site. This is relevant in the context of Mediterranean holm oak forests: up to now, most of the studies have been conducted in protected remote sites and were little representative of the fluxes in forests close to industrial activity, traffic, agriculture and residential areas.