Simulating the long-term chemistry of an upland UK catchment: heavy metals

CHUM-AM was used to investigate the behaviours of atmospherically-deposited heavy metals (Ni, Cu, Zn, Cd and Pb) in three moorland sub-catchments in Cumbria UK. The principal processes controlling cationic metals are competitive partitioning to soil organic matter, chemical interactions in solution, and chemical weathering. Metal deposition histories were generated by combining measured data for the last 30 years with local lake sediment records. For Ni, Cu, Zn and Cd, default parameters for the interactions with organic matter provided reasonable agreement between simulated and observed present-day soil metal pools and average streamwater concentrations. However, for Pb, the soil binding affinity in the model had to be increased to match the observations. Simulations suggest that weakly-sorbing metals (Ni, Zn, Cd) will respond on timescales of decades to centuries to changes in metal inputs or acidification status. More strongly-sorbing metals (Cu, Pb) will respond over centuries to millennia.     

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.     

Recovery from acidification in central Europe--observed and predicted changes of soil and streamwater chemistry in the Lysina catchment,Czech Republic

The geochemical model MAGIC was applied to estimate streamwater and soil chemistry between 1851 and 2030 at the Lysina catchment, an acid-sensitive granitic catchment covered by planted Norway spruce monoculture in the western Czech Republic. The total deposition of sulfur to the catchment was 164 meq m(-2) in 1991, but had declined to 52 meq m(-2) by 2000. Although SO2 emissions in the region declined by 90% compared to the 1980s, acidification recovery was small within the period 1990-2000. Stream pH increased only slightly (from 3.92 to 4.07), although SO4 concentration declined sharply from 568 microeq l(-1) (1990) to 232 microeq l(-1) (2000). Organic acids played an important role in streamwater buffering. According to the MAGIC prediction using deposition measured in 1999-2000, streamwater pH will increase to 4.3 and soil base saturation will increase to 6.2% by 2030 (from 5.7% in 2002). Pre-industrial pH was estimated to be 5.5 and soil base saturation 24.7%. The loss of base cations (Ca, Mg, Na, K) was caused predominantly by atmospheric acidity, but intensive forestry was responsible for approximately one third of the net base cation loss via accumulation in harvested biomass. Severely damaged sites, under continued pressure from forestry, will not return to a good environmental status in the near future (if ever) when the acid deposition input is only partially reduced.     

Long-term changes of acidifying deposition in Finland (1973-2000)

The long-term changes of acidifying deposition in Finland during the period 1973-2000 were studied using bulk deposition data from 19 stations belonging to the national monitoring network. The regional-scale approach (southern, central and northern Finland) was used for trend assessment with respect to implementation of European sulphur (S) emission reduction amendments involving deposition changes prior to (1973-1985) and after (1986-2000) the agreements (S protocols in 1985 and 1994). There were no marked changes in sulphate deposition between the 1970s and 1980s and consistent trends in 1973-1985 were not observed. Deposition of nitrogen (N) compounds, particularly NO3-N, were increasing between the 1970s and 1980s. Deposition of base cations exhibited a slight decline throughout the 1970s and 1980s. Decrease of calcium and magnesium deposition without corresponding decrease in sulphate resulted in increased acidifying potential (AP) of deposition. Due to successful implementation of S (and N) emission reduction measures, sulphate deposition has decreased substantially (30% in northern and up to 60% in southern Finland) since the late 1980s. N deposition also decreased, but less than S deposition. Base cation deposition has also declined substantially, but this decline appeared to be leveling off during the 1990s, accounting for the decrease of AP in deposition. The observed deposition pattern is in agreement with the on-going biochemical recovery of acidified small Finnish lakes taking place since the early 1990s.     

Recovery of Soil Water,Groundwater, and Streamwater From Acidification at the Swedish Integrated Monitoring Catchments

Recovery from anthropogenic acidification in streams and lakes is well documented across the northern hemisphere. In this study, we use 1996–2009 data from the four Swedish Integrated Monitoring catchments to evaluate how the declining sulfur deposition has affected sulfate, pH, acid neutralizing capacity, ionic strength, aluminum, and dissolved organic carbon in soil water, groundwater and runoff. Differences in recovery rates between catchments, between recharge and discharge areas and between soil water and groundwater are assessed. At the IM sites, atmospheric deposition is the main human impact. The chemical trends were weakly correlated to the sulfur deposition decline. Other factors, such as marine influence and catchment features, seem to be as important. Except for pH and DOC, soil water and groundwater showed similar trends. Discharge areas acted as buffers, dampening the trends in streamwater. Further monitoring and modeling of these hydraulically active sites should be encouraged.     

Impacts of afforestation on water quality trends in two catchments in mid-Wales

The impact of afforestation on stream-water chemistry for two catchments in central Wales is investigated. Trends in water chemistry are evaluated with forest age progressing at 8-15 years and 20-27 years. To assess the possible exacerbating effect of afforestation on surface-water acidification two moorland catchments are used as controls. Absolute differences in inter-catchment streamwater chemistry are evident. Differences in acidity between the moorland and the young forest site reflect site preparation and modifications to water pathways. Differences in the observed present day chemistry between the young and old forest result from past interactions between deposition changes during forest development. An estimate of the future impact of forest growth was attained from model predictions. The model was calibrated to present day stream-water chemistry and incorporates, cation uptake, evapotranspiration/ concentration effect and increased scavenging of occult and dry deposition. Predictions suggest the pH will decline to 5.3 when the young forest progresses to 20 years of age, which is still much higher than the pH recorded at the older forest site at the age of 20 years.     

Dynamic modelling of atmospherically-deposited Ni, Cu, Zn, Cd and Pb in Pennine catchments (northern England)

Simulation modelling with CHUM-AM was carried out to investigate the accumulation and release of atmospherically-deposited heavy metals (Ni, Cu, Zn, Cd and Pb) in six moorland catchments, five with organic-rich soils, one with calcareous brown earths, in the Pennine chain of northern England. The model considers two soil layers and a third layer of weathering mineral matter, and operates on a yearly timestep, driven by deposition scenarios covering the period 1400-2010. The principal processes controlling heavy metals are competitive solid-solution partitioning of solutes, chemical interactions in solution, and chemical weathering. Agreement between observed and simulated soil metal pools and surface water concentrations for recent years was generally satisfactory, the results confirming that most contemporary soil metal is from atmospheric pollution. Metals in catchments with organic-rich soils show some mobility, especially under more acid conditions, but the calcareous mineral soils have retained nearly all anthropogenic metal inputs. Complexation by dissolved organic matter and co-transport accounts for up to 80% of the Cu in surface waters.     

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.     

Soil acidification and foliar nutrient status of Ontario's deciduous forest in 1986 and 2005

To assess the impacts of the decline in sulphur (S) deposition over the past 20 years in Ontario, soil chemistry and sugar maple (Acer saccharum Marsh) foliar chemistry were measured at 17 sites in south and central Ontario in 2005 and compared with archived samples collected in 1986. Foliar S concentrations were lower in 2005, reflecting the decline in S deposition whereas foliar N remained unchanged, reflecting the lack of change in N deposition in Ontario. Mineral soil pH, exchangeable base cations were lower in 2005 whereas total S, N and cation exchange capacity (CEC) were unchanged. Foliar concentrations of Ca were positively related to soil Ca levels in the A-horizon and were lower in 2005. Despite evidence of increasing soil acidity and losses of calcium, foliar base cation concentrations at most sites were adequate for sugar maple and forest health in terms of canopy appearance (Decline Index) improved.     

Simulating Dissolved Organic Carbon Dynamics at the Swedish Integrated Monitoring Sites with the Integrated Catchments Model for Carbon,INCA-C

Surface water concentrations of dissolved organic carbon ([DOC]) are changing throughout the northern hemisphere due to changes in climate, land use and acid deposition. However, the relative importance of these drivers is unclear. Here, we use the Integrated Catchments model for Carbon (INCA-C) to simulate long-term (1996–2008) streamwater [DOC] at the four Swedish integrated monitoring (IM) sites. These are unmanaged headwater catchments with old-growth forests and no major changes in land use. Daily, seasonal and long-term variations in streamwater [DOC] driven by runoff, seasonal temperature and atmospheric sulfate (SO₄ ²⁻) deposition were observed at all sites. Using INCA-C, it was possible to reproduce observed patterns of variability in streamwater [DOC] at the four IM sites. Runoff was found to be the main short-term control on [DOC]. Seasonal patterns in [DOC] were controlled primarily by soil temperature. Measured SO₄ ²⁻ deposition explained some of the long-term [DOC] variability at all sites.     

Decreased acid deposition and the chemical recovery of Killarney,Ontario, lakes

Lakes in Killarney Park near Sudbury, Ontario, Canada, have shown dramatic water quality changes including general increases in pH and alkalinity, and decreases in SO4(2-), base cations and metals. While some lakes have recovered to pH > 6.0, many are still highly acidic despite decades of improvement. Very high historical S deposition related to emissions from the Sudbury metal smelters dominated the acidification process in this region. However, since the implementation of substantial S emission controls (90%) at the smelters, the Sudbury emissions are no longer the major source of S deposition in the Sudbury area. Wet deposition of SO4(2-) and SO4(2-) concentrations in lakewaters at Killarney now approach values in the Dorset, Ontario, area, about 200 km from Sudbury. This suggests that the S deposition to the Killarney area is now primarily from long-range transport, not from local sources. Studies of Killarney lakes are revealing the complex nature of the chemical recovery process. As lake acidity decreases, other changes including decreased Ca2+ concentrations, increased transparency, and altered thermal regimes may potentially affect some of these ecosystems. It is clear that continuing assessments of the recovery of Killarney lakes, within a multiple-stressor framework, are needed.     

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.     

Empirical and simulated critical loads for nitrogen deposition in California mixed conifer forests

Empirical critical loads (CL) for N deposition were determined from changes in epiphytic lichen communities, elevated NO(3)(-) leaching in streamwater, and reduced fine root biomass in ponderosa pine (Pinus ponderosa Dougl. ex Laws.) at sites with varying N deposition. The CL for lichen community impacts of 3.1kg ha(-1) year(-1) is expected to protect all components of the forest ecosystem from the adverse effects of N deposition. Much of the western Sierra Nevada is above the lichen-based CL, showing significant changes in lichen indicator groups. The empirical N deposition threshold and that simulated by the DayCent model for enhanced NO(3)(-)leaching were 17kg N ha(-1) year(-1). DayCent estimated that elevated NO(3)(-) leaching in the San Bernardino Mountains began in the late 1950s. Critical values for litter C:N (34.1), ponderosa pine foliar N (1.1%), and N concentrations (1.0%) in the lichen Letharia vulpina ((L.) Hue) are indicative of CL exceedance.     

Modelling long-term cation supply in acidified forest stands

A dynamic soil chemistry model was used to explain the observed decrease in soil base saturation between 1949 and 1984 at three stands in southern Sweden. The results show that acid deposition has caused soil acidification. The model, SAFE (Soil Acidification in Forest Ecosystems), includes the fundamental physical processes such as leaching and accumulation, and chemical processes such as cation exchange, mineral weathering, nutrient uptake and solute equilibrium reactions. The sources and sinks of base cations in the soil system were quantified, showing that weathering, deposition of base cations and depletion of exchangeable base cations supply cations to the soil solution in similar amounts in the upper 1 m during the acidification phase. This demonstrates that budget studies alone cannot be used to distinguish between long-term capacity to resist acidification, represented by weathering, from short-term buffering caused by cation exchange.     

Critical loads and steady-state chemistry for streams in the state of Maryland

The critical loads to streams, steady-state stream chemistry and catchment chemical weathering rate in 73 catchments has been determined in the state of Maryland, USA. It was calculated with the PROFILE model from chemical limits for biological indicators, soil mineralogy, soil texture, annual average temperature, average soil moisture, net long-term uptake of base cations and nitrogen to the vegetation, annual precipitation and runoff and deposition of sulphur and nitrogen precursors of acid deposition. The results show a full range of critical loads from very low values in the sensitive catchments of western Maryland and the Coastal Plain on the Chesapeake Bay, to insensitive catchments in the Fredrick Valley and Ridge and the Piedmont plain. The critical loads will be used as an input to an integrated regional assessment of the quantitative sensitivity of streams to acid rain, and the assessment of regional stream alkalinity response to different abatement strategies. The mapping of steady-state stream chemistry indicates that streams in Maryland are still acidfying under the present deposition load. Land-use seems to play an important role in maintaining neutral pH in many of the streams of Maryland.     

Base cation depletion, eutrophication and acidification of species-rich grasslands in response to long-term simulated nitrogen deposition

Pollutant nitrogen deposition effects on soil and foliar element concentrations were investigated in acidic and limestone grasslands, located in one of the most nitrogen and acid rain polluted regions of the UK, using plots treated for 8-10years with 35-140kg Nha(-2)y(-1) as NH(4)NO(3). Historic data suggests both grasslands have acidified over the past 50years. Nitrogen deposition treatments caused the grassland soils to lose 23-35% of their total available bases (Ca, Mg, K, and Na) and they became acidified by 0.2-0.4 pH units. Aluminium, iron and manganese were mobilised and taken up by limestone grassland forbs and were translocated down the acid grassland soil. Mineral nitrogen availability increased in both grasslands and many species showed foliar N enrichment. This study provides the first definitive evidence that nitrogen deposition depletes base cations from grassland soils. The resulting acidification, metal mobilisation and eutrophication are implicated in driving floristic changes.     

The Solling roof revisited--slow recovery from acidification observed and modeled despite a decade of "clean-rain" treatment

Soil chemistry under the Solling clean-rain roof was simulated using the dynamic multi-layer soil chemistry model SAFE, including sulfate adsorption. Soil was sampled in order to parameterize the pH and sulfate concentration dependent sulfate adsorption isotherm used in SAFE. Modeled soil solution chemistry was compared to the 14 year long time-series of monthly measurements of soil solution data at 10 and 100cm depth. The deposition of N and S under the roof has been reduced by 68% and 53%, respectively, compared to the surrounding area. Despite this the soil solution concentrations of sulfate are still high (a median of 420mumol(c)/L at 100cm depth between 2000 and 2002) and the soil base saturation low (approximately 3% in the whole profile in 1998). Sulfate adsorption is an important process in Solling. The soil capacity to adsorb sulfate is large, the modeled adsorbed pool in 2003 down to 100cm was 1030kg S/ha, and the measured sulfate concentration is high, due to release of adsorbed sulfate. The addition of sulfate adsorption improved the modeled sulfate dynamics although the model still slightly underestimated the sulfate concentration at 100cm. Model predictions show no recovery, based on the criteria of Bc/Al ratio above 1 in the rooting zone, before the year 2050, independent of future deposition cuts.     

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.     

Sampling and physico-chemical analysis of precipitation: a review

Wet deposition is one of two processes governing the transfer of beneficial and toxic chemicals from the atmosphere on to surfaces. Since the early 1970s, numerous investigators have sampled and analyzed precipitation for their chemical constituents, in the context of "acidic rain" and related atmospheric processes. Since then, significant advances have been made in our understanding of how to sample rain, cloud and fog water to preserve their physico-chemical integrity prior to analyses. Since the 1970s large-scale precipitation sampling networks have been in operation to broadly address regional and multi-regional issues. However, in examining the results from such efforts at a site-specific level, concerns have been raised about the accuracy and precision of the information gathered. There is mounting evidence to demonstrate the instability of precipitation samples (e.g. with N species) that have been subjected to prolonged ambient or field conditions. At the present time precipitation sampling procedures allow unrefrigerated or refrigerated collection of wet deposition from individual events, sequential fractions within events, in situ continuous chemical analyses in the field and even sampling of single or individual rain, cloud and fog droplets. Similarly analytical procedures of precipitation composition have advanced from time-consuming methods to rapid and simultaneous analyses of major anions and cations, from bulk samples to single droplets. For example, analytical techniques have evolved from colorimetry to ion chromatography to capillary electrophoresis. Overall, these advances allow a better understanding of heterogeneous reactions and atmospheric pollutant scavenging processes by precipitation. In addition, from an environmental perspective, these advances allow better quantification of semi-labile (e.g. NH4+, frequently its deposition values are underestimated) or labile species [e.g. S (IV)] in precipitation and measurements of toxic chemicals such as Hg and PCBs (polychlorinated biphenyls). Similarly, methods now exist for source-receptor studies, using for example, the characterization of reduced elemental states and/or the use of stable isotopes in precipitation as tracers. Future studies on the relationship between atmospheric deposition and environmental impacts must exploit these advances. This review provides a comprehensive and comparative treatment of the state of the art sampling methods of precipitation and its physico-chemical analysis.     

Nitrogen deposition in California forests: a review

Atmospheric concentrations and deposition of the major nitrogenous (N) compounds and their biological effects in California forests are reviewed. Climatic characteristics of California are summarized in light of their effects on pollutant accumulation and transport. Over large areas of the state dry deposition is of greater magnitude than wet deposition due to the arid climate. However, fog deposition can also be significant in areas where seasonal fogs and N pollution sources coincide. The dominance of dry deposition is magnified in airsheds with frequent temperature inversions such as occur in the Los Angeles Air Basin. Most of the deposition in such areas occurs in summer as a result of surface deposition of nitric acid vapor (HNO3) as well as particulate nitrate (NO3-) and ammonium (NH4+). Internal uptake of gaseous N pollutants such as nitrogen dioxide (NO2), nitric oxide (NO), HNO3, peroxyacetyl nitrate (PAN), ammonia (NH3), and others provides additional N to forests. However, summer drought and subsequent lower stomatal conductance of plants tend to limit plant utilization of gaseous N. Nitrogen deposition is much greater than S deposition in California. In locations close to photochemical smog source areas, concentrations of oxidized forms of N (NO2, HNO3, PAN) dominate, while in areas near agricultural activities the importance of reduced N forms (NH3, NH4+) significantly increases. Little data from California forests are available for most of the gaseous N pollutants. Total inorganic N deposition in the most highly-exposed forests in the Los Angeles Air Basin may be as high as 25-45 kg ha(-1) year(-1). Nitrogen deposition in these highly-exposed areas has led to N saturation of chaparral and mixed conifer stands. In N saturated forests high concentrations of NO3- are found in streamwater, soil solution, and in foliage. Nitric oxide emissions from soil and foliar N:P ratios are also high in N saturated sites. Further research is needed to determine the ecological effects of chronic N deposition, and to develop appropriate management options for protecting water quality and managing plant nutrient resources in ecosystems which no longer retain excess N.