The MAGIC simulation of surface water acidification at, and first year results from, the Bear Brook Watershed Manipulation, Maine, USA

The catchments of East and West Bear Brooks, Maine, USA, have been hydrologically and chemically monitored for 3.5 years. Stream chemistries and hydrographs are similar. These clear water streams are low in ANC (0-70 microeq litre(-1)), with variations caused by changing concentrations of base cations, SO4, NO3 and Cl. The latter range between 90-120, 0-40 and 65-75 microeq litre(-1), respectively. The West Bear catchment is being treated with six applications per year of dry (NH4)2SO4 at 1800 eq ha(-1) year(-1). After one year of treatment, the response of the stream chemistry and the response modelled by MAGIC are similar. Retentions of NH4 and SO4 are nearly 100% and greater than 80%, respectively. The additional flux of SO4 is compensated principally by an increased Ca concentration. Episodes of high discharge in the treated catchment are now characterized by lower ANC and pH, and higher Al than prior to the manipulation. Concentrations of NO3 have increased about 10 microeq litre(-1) during the dormant season, presumably due to additional nitrification of N from NH4. Discharge-chemistry relationships indicate that changes in stream chemistry, except for NO3, are dominated by ion exchange reactions in the upper part of the soil profile.     

Simulating the long-term chemistry of an upland UK catchment: major solutes and acidification

CHUM-AM was used to investigate changes in soil and water chemical variables in four moorland sub-catchments in Cumbria UK, to which non-marine S deposition has declined by 65% since the 1970s. The principal processes represented in the model comprise N and S uptake and release, water movements, the binding of cations by soil organic matter, chemical interactions in solution, and chemical weathering. CHUM-AM reproduced reasonably well the current soil pH and pools of N and S, and changes in streamwater chemistry over the period 1970-2000, notably decreases in the concentrations of alkaline earth cations and sulphate, and increases in pH. The model also predicts streamwater pH-Al relationships in agreement with observations. Predictive calculations suggest that constant atmospheric deposition of N at present rates will lead to N saturation and re-acidification, whereas a 50% reduction in N would stabilise soil and streamwater pH at about the present levels.     

Modeling recovery of Swedish ecosystems from acidification

Dynamic models complement existing time series of observations and static critical load calculations by simulating past and future development of chemistry in forest and lake ecosystems. They are used for dynamic assessment of the acidification and to produce target load functions, that describe what combinations of nitrogen and sulfur emission reductions are needed to achieve a chemical or biological criterion in a given target year. The Swedish approach has been to apply the dynamic acidification models MAGIC, to 133 lakes unaffected by agriculture and SAFE, to 645 productive forest sites. While the long-term goal is to protect 95% of the area, implementation of the Gothenburg protocol will protect approximately 75% of forest soils in the long term. After 2030, recovery will be very slow and involve only a limited geographical area. If there had been no emission reductions after 1980, 87% of the forest area would have unwanted soil status in the long term. In 1990, approximately 17% of all Swedish lakes unaffected by agriculture received an acidifying deposition above critical load. This fraction will decrease to 10% in 2010 after implementation of the Gothenburg protocol. The acidified lakes of Sweden will recover faster than the soils. According to the MAGIC model the median pre-industrial ANC of 107 microeq L(-1) in acid sensitive lakes decreased to about 60 microeq L(-1) at the peak of the acidification (1975-1990) and increases to 80 microeq L(-1) by 2010. Further increases were small, only 2 microeq L(-1) between 2010 and 2040. Protecting 95% of the lakes will require further emission reductions below the Gothenburg protocol levels. More than 7000 lakes are limed regularly in Sweden and it is unlikely that this practice can be discontinued in the near future without adverse effects on lake chemistry and biology.     

Heterogeneous response of central European streams to decreased acidic atmospheric deposition

We investigated the relations between mountain streamwater chemistry and atmospheric pollution in an arca of 1611 km2 of the Czech Republic by comparing concentrations of SO4. NO3, Cl, Ca and the pH at 432 localities at the time of high (1984-1986) and relatively low (1996-2000) acidic atmospheric deposition. Medians of Cl. SO4, and NO3 decreased by 17, 96 and 60 microeq l(-1), or by 23, 17 and 49%, respectively, during 12 +/- 2 years. The decreased Cl corresponds to decreased emission of industrial Cl (as HCl). The decreased stream SO4 was proportionally lower than the 71% decrease in S-emissions due to leaching of previously accumulated SO4 from soils and non-zero background concentrations. Decreases of NO3 up to 60% in streamwaters with pH < or = 6 was greater than the decrease of N emission in central Europe, about 35%. Extensive regional decrease of NO3 is surprising and is probably described for the first time. The difference in NO3 concentrations between the two periods was probably enhanced by (a) an increase of mineralisation of forest floor in the mid-1980s and (b) by higher uptake of N in the late-1990s. The median pH of the 432 streams did not change but the pH of the sub-population with pH < 6 in the mid-1980s recovered substantially. The pH of circum-neutral streams (pH > 6.5) decreased even as acidic atmospheric deposition decreased.     

Response of the Plastic Lake catchment, Ontario, to reduced sulphur deposition

As a consequence of decreases in the emission rate of sulphur in eastern North America in the late 1970s and early 1980s, sulphate deposition in central Ontario declined by about 40%, but has remained constant for about six years. Plastic Lake, a small, dilute lake on the Precambrian shield that the authors have studied since 1979, acidified between the start of the study and about 1986, but since then has not changed. The authors also monitored the chemistry of streamwater draining the Plastic Lake catchment. Water quality of runoff from an upland site improved rapidly (pH and alkalinity increased, SO4(2-) and Al decreased), but two factors offset these improvements. A small wetland area downstream reversed most of these changes, resulting in a constant output of strong acid from the catchment. In addition, in extremely dry years (1983, 1987, 1989) there were very high concentrations of SO4(2-) in the streamwater, suggesting substantial re-oxidation of reduced S in the catchment.     

Acid-sensitive waters of the English Lake District: a steady-state model of streamwater chemistry in the upper Duddon catchment

Data on deposition and streamwater chemistry, obtained for the upper catchment of the River Duddon in the 1970s and 1980s, are reviewed. These data, together with soil chemical data, are used to deduce key processes in the deposition-catchment interaction, the analysis being based on current concepts of acidification. The processes are incorporated into a steady-state model that allows streamwater compositions to be calculated. The large baseflow pH range (5-7) of Duddon streams is accounted for in the model by a range of base cation weathering rates. Other processes invoked are evapotranspiration, the uptake of nitrogen by plants, dissolution of Al(OH)(3) in the mineral soil, precipitation of Al(OH)(3) in the baserock zone and in streamwater, Al(3+) hydrolysis, and reactions of the carbonate system. Both cation exchange and sulphate adsorption are ignored, because they are assumed to influence rates-of-change between steady states, but not steady-state water compositions per se. The model can be used to estimate variations in streamwater composition with flow. Model calculations suggest that a 50% decrease in depositional acidifying components (sulphur oxides and NH(4)(+)) would result in increases of up to 1 pH unit in streamwaters with present-day baseflow pH values of 5 or less. It appears that water quality in the upper Duddon is currently more sensitive to inputs of NH(4)(+) than of H(2)SO(4). To improve the reliability of model predictions, more information is required on (a) the pH dependence of base cation weathering, (b) transformations involving nitrogen, (c) aluminium chemistry and (d) partial pressures of CO(2) in soil and baserock.     

Surface water acidification in the South Pennines I. Current status and spatial variability

The South Pennines, an area of acid-sensitive geology at the centre of a major industrial region, have undergone perhaps the most severe historic exposure to sulphur and nitrogen deposition in the UK. This study addresses a lack of existing research on the region by presenting the findings of a survey of 62 surface waters sampled during a 1-week period in April 1998. Results indicate that acidification in the region is acute; 27 of the sampled surface waters had a negative acid neutralising capacity (ANC) and 28 had a pH below 5.0. Minimum recorded pH values were below 4.0. Non-marine sulphate levels were extremely high (median 222 microeq/1), and widespread high nitrate concentrations (median 41 microeq/1) suggest that soils in the region as a whole may be at an advanced stage of nitrogen saturation. A consistent relationship was identified between site acidity and the balance between the major weathering-derived cations, calcium and magnesium, and sulphate. This could in turn be linked to catchment soil type and land use, with the most acidic conditions occurring in peat-dominated catchments, where weathering is minimal and the influence of atmospheric deposition most pronounced. Percentage of peat in each catchment was the single best predictor of surface water acidity. Nitrate concentrations, although not a dominant control on acidity, varied significantly according to land use. Elevated concentrations were observed in catchments containing forestry, due to enhanced deposition inputs, and in catchments containing improved land, linked to fertiliser use. Ammonium concentrations, although low at most catchments, were a significant component of the inorganic nitrogen total in a number of surface waters draining waterlogged peat catchments.     

Buffering processes in a boreal dissolved organic carbon-rich stream during experimental acidification

The role of organic acids on surface water acidity as well as their buffering during anthropogenic acidification and subsequent recovery was studied in a field experiment on a total organic carbon (TOC)-rich stream draining the Svartberget catchment in northern Sweden. H(2)SO(4) was added to the stream to increase SO(4)(2-) concentration by 90 microeq l(-1) for 30 h. About 60% of the added H(+) was buffered by protonation of organic acids, another 20% was buffered by base cations released from the surface of the stream channel and only ca. 20% of the added acid remained unbuffered. TOC concentrations (27 mg l(-1)), and site density of carboxylic groups--8.6 microeq (mg TOC)(-1)--remained stable during the experiment. Two models of organic acid dissociation (a triprotic model and a monoprotic pH-dependent pKa model) were fitted to the experimental results. These models explained the observed variations in organic anions, but the model parameters were quite different from those reported by studies in Northern America and Central Europe. This experiment had substantially more buffering effect of TOC between pH 4.4 and 5.3, which is an environmentally important pH range.     

Mechanisms of episodic acidification in low-order streams in Maine, USA

Five factors contribute to episodic depressions in pH and ANC during hydrologic events in low-order streams in Maine: (1) increases of up to 50 microeq litre(-1) NO3; (2) increases of up to 75 microeq litre(-1) organic acidity; (3) increases of as much as 0.3 in the anion fraction of SO4; (4) as much as 100 microeq litre(-1) acidity generated by the salt-effect in soils; and (5) typically < or = 40% dilution by increased discharge. In conjunction with increased discharge, factors 1, 2 or 4 appear necessary to depress pH to less than 5.0. The chemistry of individual precipitation events is irrelevant to the generation of acidic episodes, except those caused by high loading of neutral salts in coastal regions. Increases in discharge, but not necessarily in dilution of solutes, in combination with the chronically high SO4 from atmospheric deposition, provide the antecedent chemical conditions for episodic acidification. Differences in antecedent moisture conditions determine the processes that control output of either ANC or acidifying agents to aquatic systems.     

Afforestation, seasalt episodes and acidification--a paired catchment study in western Norway

As acid deposition has declined during the past 15-20 years in western Norway, afforestation and episodic seasalt deposition have become factors of increasing importance in explaining the mobilization of toxic aluminum (Aln+) to rivers and lakes. We conducted a paired catchment at four sites in western Norway across a gradient in acid deposition to evaluate the importance of afforestation and seasalt episodes. Streamwater was sampled intensively before, during and after seasalt episodes over a three-year period. A seasalt episode in January 2003 caused considerable impact on the streamwater chemistry. pH dropped and concentrations of Aln+ increased due to cation exchange of Na+ ions for H+ and Aln+ in the soil. The response was larger in streams draining the catchments which receive high acid deposition and in those afforested with spruce as compared with adjacent catchments in native birch. The results indicate that acid pulses induced by episodic inputs of seasalts are exacerbated by land use change from native birch to planted spruce.     

Long-term trends and spatial variability in nitrate leaching from alpine catchment-lake ecosystems in the Tatra Mountains (Slovakia-Poland)

Relationships between catchment characteristics of 31 alpine lakes and observed trends in lake water concentrations of nitrate were evaluated in the Tatra Mountains. Nitrate concentrations increased from background levels <4 microeql(-1) in the 1930s to maxima (up to 55 microeql(-1)) in the 1980s, after which they declined to 4-44 microeql(-1) by the late 1990s. In-lake nitrate concentrations correlated negatively with parameters characterising catchment-weighted mean pools (CWM; kgm(-2)) of soil, i.e. with percent land cover with meadow and soil depth, and positively with grade of terrain, annual precipitation, and the highest elevation in the catchment. The CWM pool of soil and annual precipitation explained together 65% of the current spatial variability in nitrate concentrations. Denitrification and direct N deposition on surface area explained 14% of the variability. Increased atmospheric N deposition and declining net N retention in soils were responsible for long-term changes in nitrate concentrations. Long-term decline in %N retention in soils decreased along with the estimated decline in C:N ratios (from 21 to 18 on average during the last 70 years). An empirical model linking nitrate concentrations in different types of alpine Tatra Mountain lakes to four independent variables (CWM soil pool, annual precipitation, increased N deposition, and average trend in soil C:N ratios) explained 80% of the observed spatial and temporal nitrate variability over the period 1937-2000.     

Contemporary (1979-1988) and inferred historical status of headwater lakes in North Central Ontario, Canada

Fifty-six headwater Canadian Shield lakes were repetitively sampled from 1979 to 88 to determine their response to changes in acidic deposition of the period. Annual wet sulphate loadings varied between 38 and 83 meq m(-2), with highest deposition in the late 1970s followed by somewhat lower but variable deposition in the 1980s. Median pH of the lakes increased 0.42 pH units from 1979 to 1985 and decreased by 0.15 units between 1985 and 1988. Short water renewal times (x=1.1 y) promoted rapid equilibration. Since lake were so responsive to changes in SO4(2-) inputs, they were at or near steady state at all times. Comparison of predicted original pH and ANC with 1979 data indicate a median decline of 0.45 pH units and a loss of 34 microeq litre(-1). ANC. Four of 9 lakes were found to be historically fishless, based on the continued presence of Chaoborus americanus in sediment cores. The remaining five lakes historically had fish populations, but fish were not collected in 1979 when pH ranged betwen 4.6 and 5.3. By 1987, fish species were found in five of these lakes where pH had increased on average by 0.9 pH units. Our data indicate that water quality improvements could allow for the reinvasion or resumption of recruitment for a significant number of Ontario lakes.     

Modelling future soil chemistry at a highly polluted forest site at Istebna in Southern Poland using the "SAFE" model

The multi-layer dynamic model SAFE was applied to the forested catchment Istebna (Southern Poland), to study recovery from acidification. Environmental pollution in the area has been historically high. The model uses data from an intensive monitoring plot established in 1999 in a spruce stand, which was planted in 1880. Observations showed that the soil was depleted of base cations. The measured base saturation in 1999 was between 5 and 8% in the different soil layers. Model predictions assuming full implementation of the UNECE 1999 Gothenburg Protocol and present day base cation deposition show that the base saturation will slowly increase to 20% by 2100. Despite large emission reductions, Istebna still suffers from the very high loads of acidifying input during the past decades. Soil recovery depends on future emissions especially on base cation deposition. The recovery will be even slower if the base cation deposition decreases further.     

Soil solution nitrogen and cations influenced by (NH4)2SO4 deposition in a coniferous forest

The effects of chronically enhanced (NH(4))(2)SO(4) deposition on ion concentrations in soil solution and ionic fluxes were investigated in a Picea abies plot at Grizedale forest, NW England. Soil cores closed at the base and containing a ceramic suction cup sampler were 'roofed' and watered every 2 weeks with bulk throughfall collected in the field. Treatments consisted of the inclusion of living roots from mature trees in the lysimeters and increasing (NH(4))(2)SO(4) deposition (NS treatment) to ambient + 75 kg N ha(-1) a(-1). Rainfall, throughfall and soil solutions were collected every 2 weeks during 18 months, and analysed for major cations and anions. NO(3)(-) fluxes significantly increased following NS treatment, and were balanced by increased Al(3+) losses. Increased SO(4)(2-) concentrations played a minor role in controlling soil solution cation concentrations. The soil exchange complex was dominated by Al and, during the experimental period, cores of all treatments 'switched' from Ca(2+) to Al(3+) leaching, leading to mean [Formula: see text] molar ratios in soil solution of NS treated cores of 0.24. The experiment confirmed that the most sensitive soils to acidification (through deposition or changing environmental conditions) are those with low base saturation, and with a pH in the lower Ca, or Al buffer ranges.     

Quantifying the role of forest soil and bedrock in the acid neutralization of surface water in steep hillslopes

The role of soil and bedrock in acid neutralizing processes has been difficult to quantify because of hydrological and biogeochemical uncertainties. To quantify those roles, hydrochemical observations were conducted at two hydrologically well-defined, steep granitic hillslopes in the Tanakami Mountains of Japan. These paired hillslopes are similar except for their soils; Fudoji is leached of base cations (base saturation <6%), while Rachidani is covered with fresh soil (base saturation >30%), because the erosion rate is 100-1000 times greater. The results showed that (1) soil solution pH at the soil-bedrock interface at Fudoji (4.3) was significantly lower than that of Rachidani (5.5), (2) the hillslope discharge pH in both hillslopes was similar (6.7-6.8), and (3) at Fudoji, 60% of the base cations leaching from the hillslope were derived from bedrock, whereas only 20% were derived from bedrock in Rachidani. Further, previously published results showed that the stream pH could not be predicted from the acid deposition rate and soil base saturation status. These results demonstrate that bedrock plays an especially important role when the overlying soil has been leached of base cations. These results indicate that while the status of soil acidification is a first-order control on vulnerability to surface water acidification, in some cases such as at Fudoji, subsurface interaction with the bedrock determines the sensitivity of surface water to acidic deposition.     

Comparative impacts of forest harvest and acid precipitation on soil and streamwater acidity

Annual H+ budgets were calculated for a catchment with a 65-year-old spruce-fir forest, and for an adjacent catchment at three years after a whole-tree harvest. In the harvested catchment sinks for H+ due to weathering and mineralization reactions were 85% greater, and sources of H+ due to biomass uptake were 60% greater than in the undisturbed catchment. There was an overall net consumption of 550 eq H+ ha(-1) year(-1) at year 3 after harvest compared to a net generation of 520 eq H+ ha(-1) year(-1) for the 65-year-old forest. Soil pH increased by 0.2-0.4 units soon after harvest, but there was little change in pH of streamwater from the harvested watershed. The whole-tree harvest resulted in a total production of about 30,000 eq H+ ha(-1) due to biomass removal. In contrast, wet and dry deposition at rates measured in this study could add more than 50,000 eq H+ ha(-1) in the 65-year period before the next harvest. Reducing the intensity of harvest may lessen long-term impacts of these sources of H+ on acidification of soils and streams.     

Assessment of Uncertainty in Long-Term Mass Balances for Acidification Assessments: A MAGIC Model Exercise

Long-term (1860–2010) catchment mass balance calculations rely on models and assumptions which are sources of uncertainty in acidification assessments. In this article, we report on an application of MAGIC to model acidification at the four Swedish IM forested catchments that have been subject to differing degrees of acidification stress. Uncertainties in the modeled mass balances were mainly associated with the deposition scenario and assumptions about sulfate adsorption and soil mass. Estimated base cation (BC) release rates (weathering) varied in a relatively narrow range of 47–62 or 42–47 meq m⁻² year⁻¹, depending on assumptions made about soil cation exchange capacity and base saturation. By varying aluminum solubility or introducing a dynamic weathering feedback that allowed BC release to increase at more acidic pHs, a systematic effect on predicted changes in acid neutralizing capacity (ΔANC ca. 10–41 μeq l⁻¹) and pH (ca. ΔpH = 0.1–0.6) at all sites was observed. More robust projections of future changes in pH and ANC are dependent on reducing uncertainties in BC release rates, the timing, and extent of natural acidification through BC uptake by plants, temporal changes in soil element pools, and fluxes of Al between compartments.     

Reversal of acidification in tributaries of the River Duddon (English Lake District) between 1970 and 1998

Long-term changes in stream water chemistry in the upper Duddon catchment (southwest Lake District, UK) were investigated. Ten streams were sampled and analysed weekly during 1998, and the results compared with data for the early 1970s and 1986. The waters exhibited a range of pH, average values for 1998 being 5.04-7.04. For all the streams, the average pH in 1998 was greater than that during 1971-73. Statistical analysis was carried out, using the 1970s data to estimate the magnitude of inter-annual variation, and taking discharge into account on the basis of antecedent rainfall. The results showed that for two of the streams the pH increase was significant at the 2.5% level, while for a further three it was significant at the 10% level. Comparison of the 1998 concentrations of nitrate and non-marine sulphate with data obtained for five streams in 1973-74 showed that average nitrate concentration had increased from 11 to 20 microeq dm(-3) while that of non-marine sulphate had decreased from 94 to 50 microeq dm(-3). For four of the streams, comparisons were also made between the 1998 data and those for 1986. In three cases, pH in 1998 was generally higher, and Al generally lower, than the values for 1986, but in the fourth case little difference was evident. The present results support observations for five nearby standing waters, strengthening the evidence for a general reversal of acidification in the southwest part of Lake District, due to a decline in the deposition of pollutant sulphur.     

Will reduced sulphur emissions under the Second Sulphur Protocol lead to recovery of acid sensitive sites in UK?

A conceptual model of the combined effects of acid deposition and land-use, Model of Acidification of Groundwater In Catchments (MAGIC), was applied to 21 upland sites in the UK Acid Waters Monitoring Network (AWMN) to assess the likely future recovery in response to the latest international agreements controlling anthropogenic sulphur emissions throughout Europe. Future estimates of sulphur deposition were generated by the Hull Acid Rain Model (HARM), based on the agreed reductions outlined in the Second Sulphur Protocol. The results indicate only a limited degree of recovery in surface-water chemistry at all sites over the next 50 years; moreover, a continuing decline in soil base status is predicted to occur at 70% of sites, resulting in longer term reacidification of surface-water at 38% of sites. However, compared with a 'business as usual' scenario the recovery is pronounced, although acidified sites will require further reductions in acidic deposition if recovery to pre-industrial chemical conditions are to be achieved. Furthermore, land-use scenarios at afforested sites suggest that replanting of felled forest will lead to a further increase in acidification. This strengthens the argument that plantation forestry should be avoided in areas considered geologically sensitive to acidic deposition.     

Terrestrial ecosystem recovery--modelling the effects of reduced acidic inputs and increased inputs of sea-salts induced by global change

The reduced emissions of acidifying sulfur and nitrogen in Europe since the late 1970s will be further reduced when the Gothenburg protocol is fully implemented by 2010. Here we address the consequences for the recovery of acidified terrestrial ecosystems using the acidification model MAGIC applied to 3 large-scale "clean rain" experiments, the so-called roof experiments at Risdalsheia, Norway; G?rdsj?n, Sweden, and Klosterhede, Denmark. Implementation of the Gothenburg protocol will initiate recovery of the soils at all 3 sites by rebuilding base saturation. The rate of recovery is small and base saturation increases less than 5% over the next 30 years. A climate-induced increase in storm severity will increase the sea-salt input to the ecosystems. This will provide additional base cations to the soils and more than double the rate of the recovery, but also lead to strong acid pulses following high sea-salt inputs as the deposited base cations exchange with the acidity stored in the soil. Future recovery of soils and runoff at acidified catchments will thus depend on the amount and rate of reduction of acid deposition, and in the case of systems near the coast, the frequency and intensity of sea-salt episodes as well.