Journal of Environmental Quality 30:1747-1756 (2001)
© 2001 American Society of Agronomy, Crop Science Society of America, and Soil Science Society of America
TECHNICAL REPORT
Surface Water Quality
Aluminum Output Fluxes from Forest Ecosystems in Europe
A Regional Assessment
N. B. Dise*,a,
E. Matznerb,
M. Armbrusterb and
J. MacDonalda
a Dep. of Earth Sciences, The Open Univ., Milton Keynes, MK7 6AA, UK
b Dep. of Soil Ecology, Bayreuth Institute for Terrestrial Ecosystem Research (BITÖK), Univ. of Bayreuth, D-95440 Bayreuth, Germany
* Corresponding author (n.b.dise{at}open.ac.uk)
Received for publication May 4, 2000.
 |
ABSTRACT
|
|---|
Data from 89 forested catchments and plots across Europe were used to define empirical relationships between aluminum leaching and input fluxes of major ions, output fluxes of major ions, ecosystem parameters such as soil pH, and combinations of these. Forests that release dissolved Al to seepage or surface waters are located primarily in areas receiving the highest loading of acid rain, and the output flux of Al shows the highest correlations to the throughfall flux of inorganic nitrogen, the output fluxes of NO-3, H+, and SO2-4, and the mineral soil pH. If the speciation of Al is taken to be Al3+ (an overestimate), Al is released in a nearly 1:1 molar charge ratio with the sum of NO-3 and SO2-4 in runoff or seepage water over a wide range of base-poor bedrock types and acid deposition across Europe. The empirical data point to a threshold range of N deposition of 80 to 150 mmolc N m-2 yr-1 and a (less clearly defined) range of S deposition of 100 to 200 mmolc SO2-4 m-2 yr-1 above which Al released from forests exceeds 100 mmolc Al m-2 yr-1. Within this threshold range, the sites that release little or no dissolved Al are those that continue to assimilate input N and/or have high soil pH (>4.5).
 |
INTRODUCTION
|
|---|
ALUMINUM IN SOLUTION is one of the most important ions determining the functioning of terrestrial and aquatic ecosystems. The toxic effects of ionic aluminum on aquatic and terrestrial biota are well documented (e.g., Baker and Schofield, 1982; Verbost et al., 1995; Cronan and Grigal, 1995) and have received increased attention since the recognition of the effects of acidic deposition on surface waters in the 1970s.
In surface water, ionic Al can be toxic to fish, invertebrates, and a wide variety of other freshwater organisms. In forests, Ulrich (1995) and Matzner and Murach (1995) concluded that high levels of aluminum ions reduce the fine root biomass of trees in relation to the aboveground biomass. In addition, the roots in soil with high ionic Al are often concentrated in the upper organic layer (where aluminum is complexed), which may be only a few centimeters thick. Such a shallow root distribution can lead to drought stress and greatly increased susceptibility to wind throw. Because of these effects, aluminum was hypothesized to be a major agent of forest damage in acid depositionaffected forests in central Europe (Ulrich, 1995). Other potential effects of Al in the soil include interfering with ion transport through root hairs and reducing diversity of soil biota.
Elevated soil solution Al is an effect of prolonged soil acidification; the Al concentration of rainfall and snowmelt is nearly always extremely low. Instead, Al is released from the soil if the pH of the soil solution declines below about 4.5. The so-called "Al-buffer-range" (Ulrich, 1981) is characterized by the release of Al ions into the soil solution as an effect of H+ buffering.
The Al chemistry of soils and soil solutions and its regulation by solid phase properties have been addressed in several papers (e.g., Dahlgren et al., 1989; Mulder and Stein, 1994; Berggren and Mulder, 1995). Aluminum is most toxic in its reactive ionic forms: Al
+2, Al(OH)2+, and Al3+. In general, the higher the charge, the more toxic the species. Monomeric aluminum, Al3+, frequently dominates in low-pH soils containing little organic matter. Although upper soil horizons may have high concentrations of total aluminum, this is often complexed with organic matter so the (toxic) reactive aluminum concentration is low.
Once in the soil solution, aluminum ions are transported with seepage to deeper layers or to ground waters and surface waters. In contact with solutions of pH > 5, Al3+ ions can act as cation acids by the precipitation of aluminum hydroxides and the related formation of H+ (Dahlgren et al., 1989). Thus, beside the direct toxic effects to biota, the transport of aluminum ions with seepage to deeper soil layers and aquatic systems and their acidification is a major detrimental effect of Al.
The modeling and prediction of Al output from acid soils has received attention in numerous papers. So far, modeling has followed a deterministic approach using equilibrium reactions of H+ and soil solids with a given pHpAl dependence. For mineral soils gibbsite dissolution often is used. More recently it was shown that the release of organically bound Al may regulate the soil solution concentration, especially in very acid soils (Berggren and Mulder, 1995; Matzner et al., 1998).
In general, prediction of the H+ and Al concentrations in soil solutions using equilibria approaches is very difficult because of the uncertainty related to the thermodynamic constants used (Schecher and Driscoll, 1987). Furthermore, disequilibria of soil solutions with respect to defined Al-containing solids has often been found (Matzner and Prenzel, 1992; Franken et al., 1995) as well as long-term changes in the Al chemistry (Mulder and Stein, 1994). Nevertheless, deterministic models predicting Al output are used today on both local and regional scales (e.g., MAGIC, SAFE, SMART), including the calculation of critical loads for acidic deposition (deVries et al., 1994; Warfinge and Sverdrup, 1992).
Here we follow a different approach. By using long-term data on fluxes of Al and other mineral ions from a large number of case studies, we define empirical relationships between ecosystem parameters and the input of ions from atmospheric deposition on the one side, and Al output with seepage and runoff on the other. Such relationships have already been developed for NO-3 output, where it has been shown that N input, the C to N ratio of the forest floor, and the mineral soil pH were largely able to explain the regional variation of NO-3 output (Dise et al., 1998a).
Considering the clear relationship between soil acidity and soil solution Al, and the demonstrated toxic affects of Al on freshwater biota, the challenge is to assess to what extent Al leaching is already initiated across European forests, and to what extent this is related to the pattern of acid (N and S) deposition and ecosystem changes in soil processes. In this paper we analyze a database of 89 forested plots and catchments across Europe to establish the regional pattern of Al output from these forests and relate output patterns to other ecosystem parameters such as the level of acid deposition, the soil chemistry, and the input and output fluxes of other ions. The empirical relations between Al output and rates of atmospheric deposition might also be of use for the evaluation of critical loads. Our hypothesis is that sites that release Al to seepage or runoff waters are those that receive the highest loading of acidifying ions (SO2-4, NO-3, NH+4), release the largest levels of these ions (SO2-4, NO-3), and/or show suppressed soil pH and other symptoms of soil acidification.
 |
METHODS
|
|---|
A subset of the Indicators of Forest Ecosystem Functioning (IFEF) database (Dise et al., 1998b) of European forested plots and catchments is used for this study. The database was compiled from published reports and contacting individual research groups across Europe (Table 1). The database has average annual fluxes of major mineral elements with throughfall, seepage, and runoff, as well as various ecosystem characteristics (e.g., chemistry of soil and vegetation) and site information (e.g., management history, average precipitation). Not all information is available for all sites.
For this study, 53 plots and 36 catchments were chosen under the criteria of at least one year's output flux of Al and throughfall Cl- fluxes within ±50% of runoffseepage Cl- fluxes (Table 2). The latter is a check that the water budget is not greatly out of balance, and assumes that Cl- is a conservative ion. Input and output fluxes of ions were calculated by the authors of each study using models appropriate to the individual sites. For plots, output fluxes of ions were generally considered to be those from the deepest lysimeter depths measured; in most plots this is the C horizon. The depths used for calculating output fluxes averaged 86 cm and ranged from 17 to 140 cm, with the shallowest soils in Scandinavian sites.
View this table:
[in this window]
[in a new window]
|
Table 2. Summary of data used in the Indicators of Forest Ecosystem Functioning (IFEF) Al database. Fluxes are annual mean values. (min.) = Minumum number of sites.
|
|
Sites are distributed across northern and central Europe, and 69 of the sites are coniferous. Most of the measurements are averages of several years within the period 19841995, although some very long-term data are also used (for example, Solling, Germany, representing up to 15 yr of measurement). Differing measurement periods may introduce a slight error when comparing regional trends in Al output fluxes, but would not play a role when comparing relationships between Al output fluxes and input, output, or ecosystem parameters that were measured at the same time at each site. Sites known to be influenced by calcareous bedrock or soils are excluded from the analysis, as we wished to consider only those forests at potential risk for Al release. Thus, the pH values of the soil B horizons are in the range of 3.5 to 5.8, with the majority of sites at pH < 5.0 (Table 2; also Fig. 4). However, we were not able to exclude all potential sources of strong buffering, and several Czech catchments in our database release relatively high amounts of calcium in streamwater, suggesting either some catchmentground water sources of Ca, or a significant calcium input in dust.
Output fluxes of Al in seepage water or runoff (henceforth abbreviated as Alout) are examined in relation to latitude, longitude, soil characteristics such as pH and organic horizon depth, and the input and output fluxes of water, nitrogen (both NH+4 and NO-3), sulfur (SO2-4), and other major ions (Table 2). Simple statistics (correlations, simple and multiple regressions) are used to identify significant relationships between Alout and (i) output fluxes of other ions, (ii) input fluxes of other ions, (iii) ecosystem parameters, and (iv) combinations of these. The ion fluxes are expressed as mmolc (mmol charge) m-2 yr-1 to test hypotheses about the composition of runoff or seepage water in forested ecosystems, that is, whether SO2-4, NO-3, Cl-, or none of these anions consistently accompany Al3+ in solution.
In order to test these hypotheses, we need to assume a charge on the Al ions. We do not have information on Al species from any of the sites, and the Alout in our database is total dissolved Al. In this study we assign the Al ions a charge of 3, bearing in mind that this will be an overestimate of the actual charge for many sites, and may be a major overestimate in sites with highly organic soils. In general, the more acid sites will be expected to have a high proportion of ionic Al as Al3+ (Cronan, 1994). Fluxes are thus calculated as 3 x (Altot) and expressed either as that or as mmolc Al m-2 yr-1 (again with the recognition that this is an upper estimate).
We consider output fluxes of >100 mmolc Al m-2 yr-1 (9 kg Al ha-1 yr-1) as the threshold value for initiation of the risk of soil acidification and potential ecosystem degradation as described in the introduction of this paper. The exact value is arbitrary, but the approximate value of 7 to 15 kg Al ha-1 yr-1 in streamwater or leachate differentiated in the ALBIOS study (Aluminum in the Biosphere; Cronan, 1994) between sites in which some effects of ecosystem acidification (e.g., base saturation <1015%) were apparent and those that were relatively unaffected by acidification. We also in this paper define sites that leach <35 mmolc NO3 m-2 yr-1 (<5 kg NO3N ha-1 yr-1) as N retaining, regardless of the N input.
We use regression analyses to identify significant relationships between Alout and other input, output, or catchment variables, and to identify variables that show a similar relationship to Alout (e.g., slopes of regression lines statistically identical). This is useful for testing hypotheses about the conditions under which forests begin to release ionic Al. In addition to this, the regression models using input fluxes or ecosystem variables can potentially be used as predictors of Alout flux; but in this case the assumptions of regression analyses (normality of the data; random scatter of residuals) must be more strictly adhered to. We do not attempt any predictions of Al fluxes in this paper. Potential nonlinear relationships are examined by logarithmic or exponential transformation of variables. In some cases we consider subsets of the data (plots only, conifer sites only, similar bedrock) to determine if they behave differently with respect to Alout from the sites considered collectively.
Our database contains both plot and catchment sites. It may be expected that the chemical and physical fates of some ions of biogeochemical interest would be different between plots and catchments. In particular, flowpaths in catchments tend to be longer than in plots, allowing more time for biological removal of nitrogen and buffering reactions on acids and Al. Furthermore, Al complexation in peaty layers and processes in riparian zones such as SO2-4 reduction and denitrification increase solution alkalinity and reduce Alout fluxes in catchments as compared with plot studies.
However, in our study, we were not able to draw any inferences about the behavior of plots versus catchments per se. The reason is that plots and catchments tend to be segregated by location and receive different amounts of acid deposition, both of which confound the detection of trends based simply on their morphology. Catchments were dominant in areas with higher relief, especially Scandinavia and southern Europe, and these received on average one-third the inorganic N (see Fig. 3a) and one-half the SO2-4 deposition (Fig. 3b) of the plots, which dominated in central Europe. (An exception to this is six catchments located in the Czech Republic, which received intermediate to high N and S deposition, similar to that of many of the plots. Some of these sites also released relatively high amounts of Ca in streamwater, as described earlier.)
Because of these correlations, in cases where we consider plots and catchments separately we cannot say whether they behave differently because they receive different amounts of acid deposition, are exposed to different climates, have different hydrology, or are affected by other factors (or combinations of these factors).
 |
RESULTS
|
|---|
Overall Distribution of Aluminum Output
Forests that release the largest amounts of Al to seepage or runoff are predominantly located across northwestern Germany (up to 590 mmolc Al m-2 yr-1) and the Netherlands (210870 mmolc Al m-2 yr-1), with lower Al output in sites located in eastern Germany (up to 480 mmolc Al m-2 yr-1), southwestern Sweden (up to 340 mmolc Al m-2 yr-1), western Denmark (up to 270 mmolc Al m-2 yr-1), and southern Norway (up to 170 mmolc Al m-2 yr-1) (Fig. 1). The lowest levels of Al output (<130 mmolc Al m-2 yr-1) occur in sites located in the rest of Scandinavia, Wales, Ireland, the Czech Republic, and southern Germany.

View larger version (104K):
[in this window]
[in a new window]
|
Fig. 1. 3 x Altot flux in seepage water or runoff survey sites; data are the average of at least one year of regular measurement. Note: sites are located at the center of the column.
|
|
Output Flux Variables Related to Aluminum Output Flux
The output flux variable related most strongly to Alout in the data is the output of NO-3 (NO-3 out; Fig. 2a, Table 3; r2 = 0.72). The r2 value is slightly higher when the output of total inorganic N is used (
out: r2 = 0.73, N = 81, p < 0.0001), but the two regression lines are not significantly different because nearly all of the output flux of inorganic N is as NO-3. The relationship between Alout and NO-3 out is much stronger than that of Alout and any other output ion (Table 3), or for Alout and water fluxes alone (not significant). Fitting a logarithmic or exponential relationship does not improve the proportion of variability in the data explained.

View larger version (16K):
[in this window]
[in a new window]
|
Fig. 2. (a) 3 x Altot flux in seepage water or runoff versus NO-3 flux in seepage water or runoff. (b) 3 x Altot flux in seepage water or runoff versus SO2-4 flux in seepage water or runoff. (c) 3 x Altot flux in seepage water or runoff versus the combined flux in seepage water or runoff.
|
|
View this table:
[in this window]
[in a new window]
|
Table 3. Single output flux variables significantly correlated to Al output flux. Fluxes expressed as mmolc m-2 yr-1.
|
|
From the regression equation using all of the data, mean fluxes of Al and NO-3 are related in a maximum of a ca. 1.5 Al to 1 NO-3 equivalent ratio (Table 3, Fig. 2a). With the exception of two sites, forest soils that release >100 mmolc m-2 yr-1 NO-3 release as much or more Al (Fig. 2a). Of these two exceptions, one site is a Czech catchment with high Ca efflux (155 mmolc Ca2+ m-2 yr-1), indicating a local strong source of Ca in input dust or through weathering, and the other is a first generation pine plantation on former heathland, which may have high NO-3 release due to disturbance but low levels of ionic Al due to organic complexation. Overall, however, there is a strong pattern in which sites that leach significant amounts of NO-3 also release significant amounts of Al. However, the reverse is not true: even at NO-3 out of <100 mmolc m-2 yr-1, substantial amounts of Al may be released from forests into leachate or runoff.
Well after NO-3, the output ions that correlated next best to the output flux of Al are SO2-4 out and H+out (r2 = 0.42 with N = 84 and r2 = 0.46 with N = 44, respectively, Table 3). Transforming the data does not improve the fit. Although many sites release Al in a roughly molar equivalent proportion to SO2-4, there is considerable scatter from sites that release relatively high amounts of SO2-4 (>100 mmolc SO2-4 m-2 yr-1), but retain Al almost completely (Fig. 2b). These sites are sulfate-saturated but either retain nitrogen, leach high amounts of calcium, or both. All three of the largest outliers (leaching >400 mmolc SO2-4 m-2 yr-1 and retaining nearly all Al) leach high levels of Ca (370800 mmolc Ca2+ m-2 yr-1), suggesting a local strong source of Ca buffering. One of these sites is also strongly N retaining (Nout = 31 mmolc N m-2 yr-1). Although several other output ions (notably NH+4 out and K+out) show significant relationships to Alout, none other than NO-3 out, SO2-4 out, and H+out show a higher regression coefficient than 0.2 (Table 3).
The best-fit two-variable regression of output ions to Alout combines NO-3 out and H+out (R2 = 0.79, N = 43, p < 0.0001). However, with nearly the same R2 value, and almost twice the number of sites, the regression combining NO-3 out and SO2-4 out is arguably a stronger relationship:
Considering only the plots (which overall receive higher acid deposition), and those sites not suspected of a local source of calcium, the slope of the relationship between Alout and the combined fluxes of NO-3 out and SO2-4 out is slightly less than 1 (Fig. 2c):
Thus, equivalents of Al are released in a maximum ratio of 0.9:1.0 with equivalents of SO2-4 and NO-3 in the sites that are most influenced by acid deposition. This ratio is a maximum estimate because we overestimate the Al3+ flux by not accounting for other Al species of lower charge. However, for most of these fairly to highly acid sites we expect the ionic aluminum to be dominated by Al3+ (Cronan, 1994).
The highest proportion of variability in output Al flux for all sites together (plots + catchments) can be explained by combining the output fluxes of NO-3, SO2-4, and H+:
 |
Throughfall Fluxes Related to Aluminum Output Flux
The best regression fit for Alout using input fluxes is with the total throughfall flux of inorganic N
, hereafter referred to as Nin (Fig. 3a, Table 4). Slightly poorer, but still highly significant relationships occur between Alout and the input of the two inorganic nitrogen ions (NO-3, NH+4) considered separately (Table 4). The relationship between Alout and throughfall input of nitrogen shows a broadly positive trend with wide scatter (Fig. 3a). Fitting a log relationship does not improve the proportion of variability in Alout explained. Alout shows poorer relationships to SO2-4 in (Fig. 3b) or (Nin + SO2-4 in) (Fig. 3c) because there are a number of sites that receive high levels of SO2-4 (and high combined levels of N + SO2-4) but release low levels of Al (<100 mmolc m-2 yr-1). These include the six Czech catchments discussed in the introduction of this paper.
View this table:
[in this window]
[in a new window]
|
Table 4. Single input flux variables significantly correlated to Al output flux. Fluxes expressed as mmolc m-2 yr-1.
|
|
The relationship shown in Fig. 3a shows that the great majority of forested plots or catchments receiving less than about 80 mmolc N m-2 yr-1 in throughfall (ca. 11 kg N ha-1 yr-1) leach less than 100 mmolc Al m-2 yr-1 (9 kg Al ha-1 yr-1) in seepage or runoff water. Conversely, the risk of leaching >100 mmolc Al m-2 yr-1 begins at throughfall inputs of inorganic N above about 80 to 150 mmolc N m-2 yr-1 (1121 kg N ha-1 yr-1). Many of the sites that receive N deposition at the high end of this range but do not leach Al also show high N retention (as seen by comparing Fig. 3a with Fig. 2a)an extreme example of this is a forested plot in southern Sweden receiving an estimated 230 mmolc m-2 yr-1 N but leaching only 13 mmolc m-2 yr-1 Al (Fig. 3a). With a high organic horizon C to N ratio (33) and high mineral soil pH (5.4) this site leaches only 21 mmolc m-2 yr-1 N. Above 250 mmolc N m-2 yr-1 (35 kg N ha-1 yr-1) all sites released more than >100 mmolc Al m-2 yr-1.
The most significant two-variable regression of input ions on Al output combined NO-3 flux in throughfall and NH+4 flux in throughfall:
The only other combination of two input ions that was more significant than the ions alone was with K+in and SO2-4 in (R2 = 0.39, N = 71, p < 0.004).
Ecosystem Characteristics Related to Aluminum Output Flux
Of the ecosystem characteristics we considered (Table 2), the strongest correlation to Alout was the B horizon soil (H+) (r2 = 0.37; Table 5) followed by the C horizon (H+) (r2 = 0.32) and the A horizon (H+) (r2 = 0.17). Another way of considering the relationship is to relate the log of Alout to the soil pH (Table 5). The regional data show a broadly negative relationship between subsoil pH and Alout (Fig. 4). Sites with relatively high-pH soils (pH > 4.5) in general release less Al than sites with relatively low-pH soils. Sites that release the highest levels of Al (Alout > 300 mmolc m-2 yr-1) all have soil B horizon pH < 4.5. Considering only plots improves this relationship significantly (Table 5, Fig. 4), in large part because of the removal of several catchments that appear to be influenced by either a ground water Ca2+ source or highly organic soils.
View this table:
[in this window]
[in a new window]
|
Table 5. Ecosystem variables significantly correlated to Alout flux (concentrations in µmolc L-1). Fluxes expressed as mmolc m-2 yr-1.
|
|
Although soil pH is clearly related to Al leaching, for the potential of predicting Al fluxes from input or site characteristics, adding soil pH does not significantly improve the regression between input fluxes of inorganic N and Alout (Table 4), which is still the best potential predictor equation. This is because Nin and soil pH are highly correlated.
Related to soil pH, there is some suggestion that bedrock type modifies the dominant effect of acid input on Alout, although data are limited. For any given input of acid precursor ions (Nin + SO2-4 in), the flux of Al from plots dominated by loamy sands and marls is generally (but not always) lower that from poorer sandstones and metamorphic rocks, which is in turn generally (but not always) lower than that from more resistant igneous rocks (Fig. 5a). (We considered plots only in this analysis as we were more confident that the bedrock type described from plots related directly to the seepage chemistry described).

View larger version (20K):
[in this window]
[in a new window]
|
Fig. 5. (a) 3 x Altot flux in seepage water versus the combined (NO-3 + NH+4 + SO2-4) flux in throughfall; data divided into major bedrock classes. Plots only. (b) 3 x Altot flux in seepage water versus the combined (NO-3 + SO2-4) flux in seepage water; data divided into major bedrock classes. Plots only.
|
|
This relationship is, however, not due to higher buffering capacity from the sites with more weatherable bedrock (base cation fluxes are not significantly different for the different bedrock types). Rather it is due to higher N (and to a lesser extent, S) retention at the sites with more weatherable soils. This can be seen in Fig. 5b, in which the shift from acid ion precursor input (Nin + SO2-4 in) (Fig. 5a) to acid ion output (NO-3 out + SO2-4 out) (Fig. 5b) collapses most of the points close to the 1:1 line. Points representing the more weatherable bedrocks move farther to the left, that is, show a greater difference between acid ion input and acid ion output (higher retention). The net result is that the Alout fluxes of nearly all sites follow the output of NO3 and SO4, regardless of bedrock type (Fig. 5b).
 |
DISCUSSION
|
|---|
Analyzing regional trends from data compilations is inevitably subject to high uncertainty. In particular, the data we use come from different studies and time periods, often using different methods for different purposes. The quality of data and local effects such as amounts of dust deposition, site history, ground water contribution to runoff, or deposition history cannot be assessed in detail. Despite this, when a large number of case studies are considered together, patterns emerge that may be due to ecosystem responses to a limited number of dominant processes (e.g., Dise et al., 1998a; Gundersen et al., 1998; deVries et al., 1999). Local or site-specific effects then become part of the unexplained variation in the data.
In the database the overall amounts of Al released to deeper soil layers or to the aquatic environment can be high: in 40% of the sites it exceeds 100 mmolc m-2 yr-1 and in highly polluted areas it can exceed 500 mmolc m-2 yr-1 (Fig. 1 and 3). In addition, there is a clear pattern of enhanced Al output across areas, such as central Germany and the Netherlands (Fig. 1), which have been identified as receiving the highest loads of acid deposition (EMEP, 1994). The areas releasing intermediate amounts of Al (western Denmark, southwestern Sweden, southern Norway) receive intermediate amounts of acid deposition, and those releasing little or no Al receive the lowest levels of acid deposition. This is supported by the regressions, which show that a large part of the variability in the output flux of Al across European forests is explained by the input fluxes of the acid-precursor ions NO-3, NH+4, and, to a lesser extent, SO2-4 (Table 4, Fig. 3). Thus, the first result from this study is the demonstration that, across Europe, Al is released preferentially from forests receiving acid (especially N) deposition.
Furthermore, the correlations demonstrated between Alout and the mobile anions NO-3 out, SO2-4 out, and H+out (Fig. 2, Table 3), and Alout and soil pH (Fig. 4, Table 5), indicate that Al is released preferentially from sites that not only receive acid deposition, but are acidified. The high correlations between the output fluxes of these acid ions is not due to autocorrelation with water fluxes, since there is no relationship between output fluxes of Al and discharge, or Alout and conservative ions such as Cl- (Table 3). Indeed, some of the sites with the highest discharge (e.g., in Scandinavia) have the lowest Al fluxes. Overall, the strong linear relationship between Alout and the combined fluxes of SO2-4 out and NO-3 out (Fig. 2c), across a variety of bedrock types (Fig. 5b), suggest that buffer mechanisms other than Al release are only of minor importance in the base-poor soils included in our data, especially in sites that are influenced by acid deposition.
Relationships between high NO-3 and SO2-4 fluxes and high Al fluxes in seepage and runoff have been shown in individual catchment or plot studies over the course of a season (especially during spring runoff or high-flow acid pulses) (e.g., McAvoy, 1989) and for concentrations rather than fluxes (deVries et al., 1999); we now demonstrate this as a regional pattern across European forests. Overall, these findings confirm the basic concepts of the effects of acid deposition on soil chemistry and, especially, on the release of Al (Reuss, 1991; Ulrich and Sumner, 1991).
The relation of Alout to the input of different chemical ions in throughfall (Table 4, Fig. 3) showed that of all the major ions, the input of inorganic N
showed the highest relationship to the output flux of ionic aluminum. The threshold value of N deposition between sites that overwhelmingly retain Al (as defined as leaching <100 mmolc m-2 yr-1) and sites that overwhelmingly leach Al is 80 to 150 mmolc N m-2 yr-1 (1121 kg N ha-1 yr-1) (Fig. 3a). Within that range, sites that show higher N retention, higher S retention, and/or higher soil pH have lower losses of Al.
There is no other input ion, combination of input ions, or site characteristic explored that so clearly differentiates between sites that are and are not at current risk of acidification (as defined in this paper by release of Al > 100 mmolc m-2 yr-1). In particular, total inputs of SO2-4, or of the acid precursor ions
(Fig. 3b,c) showed much poorer relationships to Alout. Why nitrogen input in particular appears to be such an important factor for Al leaching is not clear. It may be that N saturation and resulting acidification is more closely related to overall forest "health" than S saturation, which is widespread across forest ecosystems in Europe. Determining the processes underlying these relationships should be a priority of future research in this area.
Mean regional critical loads on a European scale have been calculated from deterministic models as 48 mmolc m-2 yr-1 for N, 77 mmolc m-2 yr-1 for S, and 80 mmolc m-2 yr-1 for acidity (deVries et al., 1994). The empirical data point to a threshold level of N deposition of 80 to 150 mmolc N m-2 yr-1 and a (less clearly defined) threshold level of S deposition of 100 to 200 mmolc SO2-4 m-2 yr-1, above which Al release from forests exceeds 100 mmolc Al m-2 yr-1. Considering that critical loads are intended to protect ecosystems over the long term and therefore have a built-in protection aspect, the data from our study match the average given by deVries et al. (1994) quite well. However, our results suggest that N + S retention, rather than base cation buffering, is the most important factor for distinguishing whether forests on base-poor soils subject to acid deposition are acidified and release dissolved aluminum to seepage or runoff.
 |
CONCLUSIONS
|
|---|
The analyses of the database of 89 forested plots and catchments across Europe suggest the following:
(i) Forests that leach the highest amount of dissolved ionic aluminum are located in central Europe, primarily Germany, the Netherlands, and the western Czech republic. A scarcity of data from other central European countries, such as Poland, Russia, and the Baltic states, as well as France and southern Europe, does not allow us currently to make inferences about these regions.
(ii) Among input or site variables, the output flux of aluminum is most closely correlated to the input flux in throughfall of inorganic nitrogen (Fig. 3a, Table 4), and then to the mineral soil pH (Fig. 4, Table 5).
(iii) If a forest leaches NO-3 at levels >100 mmolc m-2 yr-1 it is also highly likely to leach Al at levels >100 mmolc m-2 yr-1 (Fig. 2a, Table 3). The same does not necessarily hold for sites leaching SO2-4 (Fig. 2b).
(iv) Over a broad range of sites and base-poor bedrocks, equivalents of Al are released in a ca. 0.9:1.0 ratio with equivalents of NO-3 + SO2-4 (Fig. 2c).
(v) There is some suggestion that soil parent material mitigates Al release (Fig. 5a); however, this is because the richer sites have higher nitrogen and sulfur retention (Fig. 5b) rather than significantly higher buffering capacity due to weathering of base cations.
 |
ACKNOWLEDGMENTS
|
|---|
Financial support for this work was provided by the German Ministry for Education, Science, Research and Technology (BMBF, Grant BEO-51-0339476) to N.B. Dise while working as a guest scientist at BITÖK, University of Bayreuth, and by the Open University Research Development Fund to J. MacDonald. We also thank the many researchers across Europe who have contributed to the IFEF database.
 |
NOTES
|
|---|
M. Armbruster, current address: Technical Institute of Dresden, Postfach 1117, D-01735 Tharandt, Germany. J. MacDonald, current address: Scottish Environmental Protection Agency (SEPA), Stirling, FK94TR, UK.
 |
REFERENCES
|
|---|
- Baker, J.P., and C.L. Schofield. 1982. Aluminum toxicity to fish in acidic waters. Water Air Soil Pollut. 18:289309.[ISI]
- Berggren, D., and J. Mulder. 1995. The role of organic matter in controlling Al solubility in acidic mineral soil horizons. Geochim. Cosmochim. Acta 59:41674180.
- Bredemeier, M., K. Blank, J.-Y. Xu, A. Tietema, A.W. Boxman, B. Emmet, F. Moldan, P. Gundersen, P. Schleppi, and R.F. Wright. 1998. Inputoutput balances at the NITREX sites. For. Ecol. Manage. 101:5764.
- Cronan, C.S. 1994. Aluminum biogeochemistry in the ALBIOS forest ecosystems: The role of acidic deposition in aluminium cycling. p. 5181. In D.L. Godbold and A. Hütterman (ed.) Effects of acid rain on forest processes. WileyLiss, New York.
- Cronan, C.S., and D.F. Grigal. 1995. Use of Ca/Al ratios as indicators of stress in forest ecosystems. J. Environ. Qual. 24:209226.[Abstract/Free Full Text]
- Dahlgren, R.A., C.T. Driscoll, and D.C. McAvoy. 1989. Aluminium precipitation and dissolution rates in Spodosol B horizons in the northeastern USA. Soil Sci. Soc. Am. J. 53:10451052.[Abstract/Free Full Text]
- DeVries, W., G.J. Reinds, H.D. Deelstra, J.M. Klap, and E.M. Vel. 1999. Soil solution chemistry. Chapter 8. In Intensive monitoring of forest ecosystems in Europe. Forest Intensive Monitoring Coordinating Inst. (FIMCI) EC-UNECE, Brussels, Belgium.
- DeVries, W., G.J. Reinds, and M. Posch. 1994. Assessment of critical loads and their exceedance on European forests using a one layer steady state model. Water Air Soil Pollut. 72:357394.
- Dise, N.B., E. Matzner, and M. Forsius. 1998a. Evaluation of organic horizon C:N ratio as an indicator of nitrate output in conifer forests across Europe. Environ. Pollut. 102:453456.
- Dise, N.B., E. Matzner, and P. Gundersen. 1998b. Synthesis of nitrogen pools and fluxes from European forest ecosystems. Water Air Soil Pollut. 105:143154.
- EMEP. 1994. Transboundary acidifying pollution in Europe: Calculated fields and budgets 19851993. EMEP/MSC-W Rep. 1/94 [Online]. Available at http://www.emep.int/ladm.html (verified 21 Feb. 2001).
- Forsius, M., S. Kleemola, M. Starr, and T. Ruoho-Airola. 1995. Ion mass budgets for small forested catchments in Finland. Water Air Soil Pollut. 79:1938.
- Fottova, D., and I. Skorepova. 1998. Changes in mass elements fluxes and their importance for critical loads: Geomon network, Czech Republic. Water Air Soil Pollut. 105:365376.
- Franken, G., M. Pijpers, and E. Matzner. 1995. Al chemistry of the soil solution in an acid forest soil as influenced by percolation rate and soil structure. Eur. J. Soil Sci. 46:613619.
- Gundersen, P. 1995. Nitrogen deposition and leaching in European forestsPreliminary results from a data compilation. Water Air Soil Pollut. 85:11791184.
- Gundersen, P., B.A. Emmett, O.J. Kjonaas, C.J. Koopmans, and A. Tietema. 1998. Impact of nitrogen deposition on nitrogen cycling in forests: A synthesis of NITREX data. For. Ecol. Manage. 101: 3755.
- Hauhs, M., K. Rost-Siebert, G. Raben, T. Paces, and B. Vigerust. 1989. Summary of European data. p. 5.15.37. In J.L. Malanchuk and J. Nilsson. (ed.) The role of nitrogen in the acidification of soils and surface waters. Miljörapport 89 (10). Council of Nordic Ministers.
- Johnson, D.W. 1995. Temporal patterns in beech forest soil solutions: Field and model results compared. Soil Sci. Soc. Am. J. 59:17321740.[Abstract/Free Full Text]
- Kleemola, S., and M. Forsius (ed.) 1996. International cooperative programme on integrated monitoring of air pollutant effects on ecosystems. 5th Annu. Rep. The Finnish Environ. 27. Finnish Environ. Inst., Helsinki, Finland.
- Kleemola, S., and M. Forsius (ed.) 1997. International cooperative programme on integrated monitoring of air pollutant effects on ecosystems. 6th Annu. Rep. The Finnish Environ. 116. Finnish Environ. Inst., Helsinki, Finland.
- Kleemola, S., and M. Forsius (ed.) 1999. International cooperative programme on integrated monitoring of air pollutant effects on ecosystems. 8th Annu. Rep. The Finnish Environment 325. Finnish Environ. Inst., Helsinki, Finland.
- Kölling, C. 1999. Air pollution and its effects on Bavarian forests. Results of the element budget calculations at the Bavarian Forest Climate Station 19911998 (In German). Berichte aus der Bayerischen Landesanstalt für Wald und Forstwirtschaft 22:146.
- Kreutzer, K., C. Beier, M. Bredemeier, K. Blank, T. Cummins, E.P. Farrell, N. Lammersdorf, L. Rasmussen, A. Rothe, P.H.B. de Visser, W. Weis, T. Weis, and Y.-J. Xu. 1998. Atmospheric deposition and soil acidification in five coniferous forest ecosystems: A comparison of the control plots of the EXMAN sites. For. Ecol. Manage. 101:125142.
- Matzner, E., and D. Murach. 1995. Soil changes induced by air pollutant deposition and their implication for forests in Central Europe. Water Air Soil Pollut. 85:6373.
- Matzner, E., M. Pijpers, W. Holland, and B. Manderscheid. 1998. Al in soil solutions of forest soils: Influence of water flow and soil Al-pools. Soil Sci. Soc. Am. J. 62:445454.[Abstract/Free Full Text]
- Matzner, E., and J. Prenzel. 1992. Acid deposition in the German Solling area: Effects on soil solution chemistry and Al-mobilization. Water Air Soil Pollut. 61:221234.
- McAvoy, D.C. 1989. Episodic response of aluminium chemistry in an acid sensitive Massachusetts catchment. Water Resour. Res. 25:233240.
- Mulder, J., and A. Stein. 1994. The solubility of Al in acidic forest soils: Long-term changes due to acid deposition. Geochim. Cosmochim. Acta 58:8594.
- Paces, T. 1985. Sources of acidification in central Europe estimated from elemental budgets in small basins. Nature 315:3134.
- Reuss, J.O. 1991. The transfer of acidity from soils to surface waters. p. 203218. In B. Ulrich and M.E. Sumner (ed.) Soil acidity. SpringerVerlag, Berlin, Germany.
- Schecher, W.D., and C.T. Driscoll. 1987. An evaluation of uncertainty associated with Al equilibrium calculations. Water Resour. Res. 23: 525534.
- Ulrich, B. 1981. Ecological grouping of soils according to their chemical characteristics. (In German.) Z. Pflanzenernähr. Bodenkd. 144 Band Heft3 S. 289305.
- Ulrich, B. 1995. The history and possible causes of forest decline in central Europe, with particular attention to the German situation. Environ. Rev. 3:262276.
- Ulrich, B., and M.E. Sumner. 1991. Soil acidity. SpringerVerlag, Berlin, Germany.
- Verbost, P.M., M.H.G. Berntssen, F. Kroglund, E. Lydersen, H.E. Witters, B.O. Rosseland, B. Salbu, and S.E. Wendelaar Bonga. 1995. The toxic mixing zone of neutral and acidic river water: Acute aluminum toxicity in brown trout (Salmo trutta L.). Water Air Soil Pollut. 85:341346.
- Warfinge, P., and H. Sverdrup. 1992. Calculating critical loads of acid deposition with PROFILE: A steady-state soil chemistry model. Water Air Soil Pollut. 63:119143.
This article has been cited by other articles:

|
 |

|
 |
 
K. A. Aronsson and N. G. A. Ekelund
Biological Effects of Wood Ash Application to Forest and Aquatic Ecosystems
J. Environ. Qual.,
September 1, 2004;
33(5):
1595 - 1605.
[Abstract]
[Full Text]
[PDF]
|
 |
|