Assessing recovery from acidification of European surface waters in the year 2010: evaluation of projections made with the MAGIC model in 1995

Sizeable net export of base cations from a Carpathian flysch catchment indicates their geogenic origin while the 26Mg/24Mg, 44Ca/40Ca and 87Sr/86Sr isotope ratios in runoff are indistinguishable from atmospheric input

Nutrient imbalances may negatively affect the health status of forests exposed to multiple stress factors, including drought and bark beetle calamities. We studied the origin of base cations in runoff from a small Carpathian catchment underlain by base-poor flysch turbidites using magnesium (Mg), calcium (Ca) and strontium (Sr) isotope composition of 10 ecosystem compartments. Our objective was to constrain conclusions drawn from long-term hydrochemical monitoring of inputs and outputs. Annual export of Mg, Ca and Sr exceeds 5-to-15 times their atmospheric input. Mass budgets per se thus indicate sizeable net leaching of Mg, Ca and Sr from bedrock sandstones and claystones. Surprisingly, δ26Mg, δ44Ca and 87Sr/86Sr isotope ratios of runoff were practically identical to those of atmospheric deposition and soil water but significantly different from bedrock isotope ratios. We did not find any carbonates in the studied area as a hypothetical, easily dissolvable source of base cations whose isotope composition might corroborate the predominance of geogenic base cations in the runoff. Marine carbonates typically have lower δ26 Mg and 87Sr/86Sr ratios, and silicate sediments often have higher δ26Mg and 87Sr/86Sr ratios than runoff at the study site. Mixing of these two sources, if confirmed, could reconcile the flux and isotope data.

Soil solution data from Bohemian headwater catchments record atmospheric metal deposition and legacy pollution

Soil solution chemistry depends largely on mineralogy and organic matter properties of soil horizons with which they interact. Differing lithologies within a given catchment area can influence variability in soil cation exchange capacities and affect solute transport. Zero-tension and tension lysimeters were used to evaluate the fast transport of solutes in the topsoil vs. slow diffusional matrix flow at the subsoil of three contrasting lithology catchments in a mid-elevation mountain forest. Our aim was to test the feasibility of lysimeters' hydrochemical data as a gauge for legacy subsoil pollution. Due to contrasting lithologies, atmospheric legacy pollution prevailing at the soil-regolith interface is differently yet consistently reflected by beryllium, lead, and chromium soil solution concentrations of the three catchments. Geochemical (dis)equilibrium between the soil and soil matrix water governed the hydrochemistry of the soil solutions at the time of collection, potentially contributing to decreased dissolved concentrations with increased depths at sites with higher soil pH. A complementary isotopic δ18O runoff generation model constrained potential seasonal responses and pointed to sufficiently long water-regolith interactions as to permit important seasonal contributions of groundwater enriched in chemical species to the topsoil levels. Our study also reflects subsoil equilibration with atmospheric solutes deposited at the topsoil and thus provides guidance for evaluating legacy pollution in soil profiles derived from contrasting lithology.

Forest growth responds more to air pollution than soil acidification

The forests of central Europe have undergone remarkable transitions in the past 40 years as air quality has improved dramatically. Retrospective analysis of Norway spruce (Picea abies) tree rings in the Czech Republic shows that air pollution (e.g. SO2 concentrations, high acidic deposition to the forest canopy) plays a dominant role in driving forest health. Extensive soil acidification occurred in the highly polluted "Black Triangle" in Central Europe, and upper mineral soils are still acidified. In contrast, acidic atmospheric deposition declined by 80% and atmospheric SO2 concentration by 90% between the late 1980s and 2010s. In this study we oserved that annual tree ring width (TRW) declined in the 1970s and subsequently recovered in the 1990s, tracking SO2 concentrations closely. Furthermore, recovery of TRW was similar in unlimed and limed stands. Despite large increases in soil base saturation, as well as soil pH, as a result of repeated liming starting in 1981, TRW growth was similar in limed and unlimed plots. TRW recovery was interrupted in 1996 when highly acidic rime (originating from more pronounced decline of alkaline dust than SO2 from local power plants) injured the spruce canopy, but recovered soon to the pre-episode growth. Across the long-term site history, changes in soil chemistry (pH, base saturation, Bc/Al soil solution ratio) cannot explain observed changes in TRW at the two study sites where we tracked soil chemistry. Instead, statistically significant recovery in TRW is linked to the trajectory of annual SO2 concentrations or sulfur deposition at all three stands.

Mercury methylation and methylmercury demethylation in boreal lake sediment with legacy sulphate pollution

Sulphate and dissolved organic matter (DOM) in freshwater systems may regulate the formation of methylmercury (MeHg), a potent neurotoxin that biomagnifies in aquatic ecosystems. While many boreal lakes continue to recover from decades of elevated atmospheric sulphate deposition, little research has examined whether historically high sulphate concentrations can result in persistently elevated MeHg production and accumulation in aquatic systems. This study used sediment from a historically sulphate-impacted lake and an adjacent reference lake in northwestern Ontario, Canada to investigate the legacy effects of sulphate pollution, as well as the effects of newly added sulphate, natural organic matter (NOM) of varying sulphur content and a sulphate reducing bacteria (SRB) inhibitor on enhancing or inhibiting the Hg methylation and demethylation activity (K_meth and K_demeth) in the sediment. We found that K_meth and MeHg concentrations in sulphate-impacted lake sediment were significantly greater than in reference lake sediment. Further adding sulphate or NOM with different sulphur content to sediment of both lakes did not significantly change K_meth. The addition of a SRB inhibitor resulted in lower K_meth only in sulphate-impacted sediment, but methylation was not entirely depressed. Methylmercury demethylation potentials in sediment were consistent across lakes and experimental treatments, except for some impacts related to SRB inhibitor additions in the reference lake sediment. Overall, a broader community of microbes beyond SRB may be methylating Hg and demethylating MeHg in this system. This study reveals that legacies of sulphate pollution in boreal lakes may persist for decades in stimulating elevated Hg methylation in sediment.

Diverse effects of accelerating climate change on chemical recovery of alpine lakes from acidic deposition in soil-rich versus scree-rich catchments

The current recovery of mountain lakes from atmospheric acidification is increasingly affected (both accelerated and/or delayed) by climate change. We evaluated long-term trends in the ionic composition of 30 lakes situated in the alpine zone of the Tatra Mountains, and compared the rates of their recovery with model (MAGIC) simulations done 20 years ago for the 2003-2020 period. The observed recovery was faster than the model forecast, due to greater reductions in acidic deposition than projected. Trends in water composition were further modified by climate change. Rising temperatures increased the length of the growing season and retention of inorganic N and SO₄^2- more in soil-rich compared with soil-poor catchments. In contrast, elevated precipitation and an increase in rainfall intensity reduced water residence time in soils, and consequently reduced N retention, especially in soil-poor catchments. It is likely that increases in rainfall intensity and annual number of days without snow, along with air temperatures fluctuating around the freezing point elevated the physical erosion of rocks, especially in high-elevation, steep, and scree-rich areas where rocks are not thermally insulated and stabilized by soils. Weathering of exposed accessory calcite in the eroded granodiorite bedrock was a source of Ca²⁺ and HCO₃^-, while S-bearing minerals likely contributed to lake water SO₄^2- and partly mitigated its deposition-related decrease in scree-rich catchments. The extent of climate effects on changes in the water composition of alpine lakes recovering from acidic deposition thus depended on elevation and cover of soil and scree in catchments. Our results highlight the need for incorporating dominant climate-related process into existing process-based models to increase their reliability in predicting the future development of lake water composition.

Winter arsenic pollution in 10 forest ecosystems in the mountainous border regions of the Czech Republic

Arsenic (As) concentrations and deposition fluxes were measured in snow and rime at 10 mountain-top sites near the borders between the Czech Republic and Austria, Germany, Poland, and Slovakia during three consecutive winter seasons (2009-2011). Our study was performed at a time following several decades of sharply decreasing regional atmospheric pollution and following the 2006 implementation of stricter air quality standards across Europe. Our objective was to compare vertical and horizontal depositions of soluble and insoluble As forms throughout the Czech Republic and define a recent Central European As pollution gradient. Arsenic soluble in weak nitric acid contributed 83 to 85% to the total As deposition, with the remaining 17-15% bound to stable particulate forms. The highest As deposition rates were recorded in the eastern Czech Republic near the borders with Poland and Slovakia. Complementary hydrochemical monitoring in four mountain-slope catchments situated near selected main study sites revealed a further decrease in open-area As deposition by the end of 2018 in the east of the country. In contrast, spruce canopy throughfall flux did not change significantly between 2009-2011 and 2016-2018. The site-specific relative roles of coal-burning-derived and ore-smelting-derived atmospheric As are discussed.

Controls on δ26Mg variability in three Central European headwater catchments characterized by contrasting bedrock chemistry and contrasting inputs of atmospheric pollutants

Magnesium isotope ratios (²⁶Mg/²⁴Mg) can provide insights into the origin of Mg pools and fluxes in catchments where Mg sources have distinct isotope compositions, and the direction and magnitude of Mg isotope fractionations are known. Variability in Mg isotope compositions was investigated in three small, spruce-forested catchments in the Czech Republic (Central Europe) situated along an industrial pollution gradient. The following combinations of catchment characteristics were selected for the study: low-Mg bedrock + low Mg deposition (site LYS, underlain by leucogranite); high-Mg bedrock + low Mg deposition (site PLB, underlain by serpentinite), and low-Mg bedrock + high Mg deposition (site UDL, underlain by orthogneiss). UDL, affected by spruce die-back due to acid rain, was the only investigated site where dolomite was applied to mitigate forest decline. The δ²⁶Mg values of 10 catchment compartments were determined on pooled subsamples. At LYS, a wide range of δ²⁶Mg values was observed across the compartments, from -3.38 ‰ (bedrock) to -2.88 ‰ (soil), -1.48% (open-area precipitation), -1.34 ‰ (throughfall), -1.19 ‰ (soil water), -0.99 ‰ (xylem), -0.95 ‰ (needles), -0.82 ‰ (bark), -0.76 ‰ (fine roots), and -0.76 ‰ (runoff). The δ²⁶Mg values at UDL spanned 1.32 ‰ and were thus less variable, compared to LYS. Magnesium at PLB was isotopically relatively homogeneous. The δ²⁶Mg systematics was consistent with geogenic control of runoff Mg at PLB. Mainly atmospheric/biological control of runoff Mg was indicated at UDL, and possibly also at LYS. Our sites did not exhibit the combination of low-δ²⁶Mg runoff and high-δ²⁶Mg weathering products (secondary clay minerals) reported from several previously studied sites. Six years after the end of liming at UDL, Mg derived from dolomite was isotopically undetectable in runoff.

Dynamic Modeling and Target Loads of Sulfur and Nitrogen for Surface Waters in Finland, Norway, Sweden, and the United Kingdom

The target load concept is an extension of the critical load concept of air pollution inputs to ecosystems. The advantage of target loads over critical loads is that one can define the deposition and the point in time (target year) when the critical (chemical) limit is no longer violated. This information on the timing of recovery requires dynamic modeling. Using a well-documented dynamic model, target loads for acidic deposition were determined for 848 surface waters across Finland, Norway, Sweden, and the United Kingdom for the target year 2050. In the majority of sites ( n = 675), the critical ANC-limit was predicted to be achieved by 2050; however, for 127 sites, target loads were determined. In addition, 46 sites were infeasible, i.e., even a reduction of anthropogenic deposition to zero would not achieve the limit by 2050. The average maximum target load for sulfur was 38% lower than the respective critical load across the study lakes ( n = 127). Target loads on a large regional scale can inform effects-based emission reduction policies; the current assessment suggests that reductions beyond the Gothenburg Protocol are required to ensure surface water recovery from acidification by 2050.

Long-term changes (1990-2015) in the atmospheric deposition and runoff water chemistry of sulphate, inorganic nitrogen and acidity for forested catchments in Europe in relation to changes in emissions and hydrometeorological conditions

The international Long-Term Ecological Research Network (ILTER) encompasses hundreds of long-term research/monitoring sites located in a wide array of ecosystems that can help us understand environmental change across the globe. We evaluated long-term trends (1990-2015) for bulk deposition, throughfall and runoff water chemistry and fluxes, and climatic variables in 25 forested catchments in Europe belonging to the UNECE International Cooperative Programme on Integrated Monitoring of Air Pollution Effects on Ecosystems (ICP IM). Many of the IM sites form part of the monitoring infrastructures of this larger ILTER network. Trends were evaluated for monthly concentrations of non-marine (anthropogenic fraction, denoted as x) sulphate (xSO₄) and base cations x(Ca+Mg), hydrogen ion (H⁺), inorganic N (NO₃ and NH₄) and ANC (Acid Neutralising Capacity) and their respective fluxes into and out of the catchments and for monthly precipitation, runoff and air temperature. A significant decrease of xSO₄ deposition resulted in decreases in concentrations and fluxes of xSO₄ in runoff, being significant at 90% and 60% of the sites, respectively. Bulk deposition of NO₃ and NH₄ decreased significantly at 60-80% (concentrations) and 40-60% (fluxes) of the sites. Concentrations and fluxes of NO₃ in runoff decreased at 73% and 63% of the sites, respectively, and NO₃ concentrations decreased significantly at 50% of the sites. Thus, the LTER/ICP IM network confirms the positive effects of the emission reductions in Europe. Air temperature increased significantly at 61% of the sites, while trends for precipitation and runoff were rarely significant. The site-specific variation of xSO₄ concentrations in runoff was most strongly explained by deposition. Climatic variables and deposition explained the variation of inorganic N concentrations in runoff at single sites poorly, and as yet there are no clear signs of a consistent deposition-driven or climate-driven increase in inorganic N exports in the catchments.

Morphological and immunohistochemical reactions of the larval epidermis in the Italian newt (Lissotriton italicus) after exposure to low pH

Mounting evidence suggests that amphibians are globally and currently the most threatened group of vertebrates and different causes might be responsible for this phenomenon. Acidification of water bodies is a global environmental issue that has been proposed as a possible cause for amphibian populations decline. Indeed, it has been widely demonstrated that low pH may exert harmful effects on amphibians, either directly or by increasing the adverse effects of other stressors. Surprisingly only few studies documented the response of amphibian integument to acidic pH conditions and no data are available on the effects of a non-lethal level of pH onto the amphibian larval epidermis. The present study showed that acidic pH (4.5) condition has severe effects on the epidermis of the Italian newt (Lissotriton italicus, formerly Triturus italicus) inducing both morphological and functional alterations. The increase of mucus is the first evident effect of acid injury, followed by the flattening of the epithelium and the appearance of a keratinized shedding layer. The immunolabeling of cytokeratins substantially changes acquiring an adult-like pattern. Also aquaporin 3 and iNOS expression modify their distribution according to a change of the histological features of the epidermis. These results clearly indicate that a short-term exposure to a sub-lethal pH disrupts the epidermis morphology and function in L. italicus larvae. Since the skin exerts a prominent role in both respiration and osmoregulation, the described alterations may adversely affect the overall ionic balance, with a long chain of cascading effects significantly decreasing newts survival probabilities when environmental pH lowering occurs.

Acid rain footprint three decades after peak deposition: Long-term recovery from pollutant sulphate in the Uhlirska catchment (Czech Republic)

The granitic Uhlirska headwater catchment with a size of 1.78km² is located in the Jizera Mountains in the northern Czech Republic and received among the highest inputs of anthropogenic acid depositions in Europe. An analysis of sulphate (SO₄^2-) distribution in deposition, soil water, stream water and groundwater compartments allowed to establish a SO₄^2- mass-balance (deposition input minus surface water export) and helped to evaluate which changes occurred since the last evaluation of the catchment in 1997. The determined SO₄^2- concentrations decreased in the following order: wetland groundwater>groundwater from 20m below ground level (bgl)>groundwater from 30m bgl>stream water>groundwater from10m bgl>hillslope soil water>wetland soil water>bulk deposition with median values of 0.24, 0.21, 0.17, 0.15, 0.11, 0.07, 0.03 and 0.01mmolL^-1, respectively. Our results show that average deposition reductions of 62% did not result in equal changes of the sulphate mass-balance, which changed by only 47%. This difference occurs because sulphate originates from internal sources such as the groundwater and soil water. The Uhlirska catchment is subject to delayed recovery from anthropogenic acid depositions and remains a net source of stored sulphur even after three decades of declining inputs. The wetland groundwater and soil water provide environmental memories of legacy pollutant sulphate. Elevated stream water sulphate concentrations after the unusually dry summer 2015 imply importance of weather and climate patterns for future recovery from acidification.

Long-Term Trends in Acid Neutralizing Capacity under Increasing Acidic Deposition: A Special Example of Eutrophic Taihu Lake, China

While North America and Europe have been recovering from acidification, China is experiencing impacts of acid deposition. The Taihu region is a seriously impacted area by acid rain in China, with the average rainfall pH < 5. However, the acid neutralizing capacity (ANC) and pH of Taihu Lake have significantly increased over the past 60 years (p < 0.05). Analyses showed that watershed neutralization by carbonates and in-lake alkalinization by algae activities were the two major reactions responsible for the increase. In the Taihu basin, the dominant carbonate bedrocks are the major source of base cations (particularly Ca²⁺ and Mg²⁺) and act as the acidification buffer. In addition, our field measurements across the lake showed that the pH values were significantly higher in algal bloom waters than in areas without blooms. This observation was further supported by our statistical analysis showing that the Taihu ANC and pH were significantly correlated with the chlorophyll increase (p < 0.05; 1985-2015). However, our regression analysis indicated that the base cations in the watershed would be depleted by the early 2040s if the acid deposition continues at the current rate. Our results suggest that interactions between human accelerated weathering, watershed geochemistry, and in-lake algae activities significantly impact the water chemistry of the lake. We urgently recommend an "integrated and balanced" recovery plan for the lake ecosystem.

Cu-Zn isotope constraints on the provenance of air pollution in Central Europe: Using soluble and insoluble particles in snow and rime

Copper (Cu) and zinc (Zn) isotope ratios can be used to fingerprint sources and dispersion pathways of pollutants in the environment. Little is known, however, about the potential of δ⁶⁵Cu and δ⁶⁶Zn values in liquid and solid forms of atmospheric deposition to distinguish between geogenic, industrial, local and remote sources of these potentially toxic base metals. Here we present Cu-Zn deposition fluxes at 10 mountain-top sites in the Czech Republic, a region affected by extremely high industrial emission rates 25 years ago. Additionally, we monitored isotope composition of Cu and Zn in vertical and horizontal atmospheric deposition at two sites. We compared δ⁶⁵Cu and δ⁶⁶Zn values in snow and rime, extracted by diluted HNO₃ and concentrated HF. Cu and Zn isotope signatures of industrial pollution sources were also determined. Cu and Zn deposition fluxes at all study sites were minute. The mean δ⁶⁵Cu value of atmospheric deposition (-0.07‰) was higher than the mean δ⁶⁵Cu value of pollution sources (-1.17‰). The variability in δ⁶⁵Cu values of atmospheric deposition was lower, compared to the pollution sources. The mean δ⁶⁶Zn value of atmospheric deposition (-0.09‰) was slightly higher than the mean δ⁶⁶Zn value of pollution sources (-0.23‰). The variability in δ⁶⁶Zn values of atmospheric deposition was indistinguishable from that of pollution sources. The largest isotope differences (0.35‰) were observed between the insoluble and soluble fractions of atmospheric deposition. These differences may result from different sources of Cu/Zn for each fraction. The difference in isotope composition of soluble and insoluble particles appears to be a promising tool for pollution provenance studies in Central Europe.

Pollution control enhanced spruce growth in the "Black Triangle" near the Czech-Polish border

Norway spruce (Picea abies (L.) Karst.) stands in certain areas of Central Europe have experienced substantial dieback since the 1970s. Understanding the reasons for this decline and reexamining the response of forests to acid deposition reduction remains challenging because of a lack of long and well-replicated tree-ring width chronologies. Here, spruce from a subalpine area heavily affected by acid deposition (from both sulfur and nitrogen compounds) is evaluated. Tree-ring width measurements from 98 trees between 1000 and 1350m above sea level (a.s.l.) reflected significant May-July temperature signals. Since the 1970s, acid deposition has reduced the growth-climate relationship. Efficient pollution control together with a warmer but not drier climate most likely caused the increased growth of spruce stands in this region, the so-called "Black Triangle," in the 1990s.