|
|
||||||||
a Centre for Forest, Landscape and Planning, Royal Veterinary and Agricultural University, DK-1871 Frederiksberg C, Denmark
b Centre for Forest, Landscape and Planning, Danish Forest and Landscape Research Institute, Hørsholm Kongevej 11 DK-2970 Hørsholm, Denmark
* Corresponding author (hv{at}kvl.dk)
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: SOM, soil organic matter
| INTRODUCTION |
|---|
|
|
|---|
Biomass in living plants is accumulated by the autotrophic synthesis of organic compounds. This process incorporates both C and N in the living biomass making the two elements intimately associated in both living and dead organic material. Because of this relationship the C/N ratio is used to characterize living and dead organic matter. The total C accumulated in biomass is generally limited by the availability of N. Forest trees are often adapted to N limited conditions by mechanisms that exert tight controls on the N cycles (Melillo, 1981; Aber et al., 1989). In these N limited ecosystems, increased supplies of plant available N from deposition may result in increased biomass production, and consequently additional fixation of C (Mäkipää et al., 1999). Holland et al. (1997) estimated that N deposition alone might allow for the sequestration of an additional 1.4 x 103 to 2 x 103 Mg C yr-1 worldwide. Therefore, N input was seldom considered harmful in terms of biomass production. However, since the early 1980s, the focus has been on the negative impacts of deposition of atmospheric N in forest ecosystems (Gundersen and Rasmussen, 1995), including acidification, induction of nutrient deficiencies, leaching of essential nutrients and harmful elements, requiring attention and countermeasures (Skeffington, 1990).
Soil organic matter may diminish because of changes in climate, land use, and land management (Harrison et al., 1995; Murillo, 1994; Schlesinger, 1995; Lal et al., 1995). The size of the annual input of organic material and decomposition rates are key factors regulating SOM pool sizes. Because these factors depend on forest productivity and soil fertility, soil properties are the primary factor determining the C content of forest soils within climatic regions (same moisture and temperature regimes as proposed by Soil Taxonomy) (Tate, 1992; Eswaran et al., 1995).
The role of SOM as a C and N reservoir, the sensitivity of SOM to environmental changes, and the use of C/N ratios as an indicator of ecosystem stability have necessitated precise estimations on the soil C and N pools worldwide.
Our objective in this study was to estimate the size and distribution of organic C and N pools in well-drained Danish forest soils. We gathered all data available; soil profiles characteristics, and chemical and physical variables that are related to C and N storage in soils, we then analyzed and discussed these data to establish relationships with soil forming processes.
We examined hypotheses concerning the relationship between accumulation of C and N and soil order; and soil order differences in C and N concentration and C/N ratio changes with soil depth. Further, we examine whether the forest floor C and N storage is influenced by soil order.
| MATERIAL AND METHODS |
|---|
|
|
|---|
Data Acquisition
Data for this study were acquired from 140 forest soil profiles from pedological surveys performed between 1970 and 1998 (Dalsgaard et al., 1981; Dalsgaard et al., 1991; Raulund-Rasmussen, 1993; Breuning Madsen and Olsson, 1995; Raulund-Rasmussen and Vejre, 1995; Vesterdal et al., 1995; Vejre and Hoppe 1998; Vesterdal and Raulund-Rasmussen 1998). Soil profile data included in this study had: detailed profile description, samples collected in pedogenic horizons, total C or N data, and exact location mapped. Precise profile description and horizon depth allowed classification to soil order according to Soil Taxonomy (Soil Survey Staff, 1992). All sites were forested and classified as moderately well drained or better, to avoid hydromorphic and organic soils. See Figure 1
for site locations. Most sites were plantation forests, except 22 sites located in semi-natural or managed beech (Fagus Silvatica L.) forest; all had been forested for the last 100 yr.
|
Soil Analysis and Classification
A total of 106 profiles had data on C concentrations in the O horizon and mineral soil to a depth of 100 cm. Nitrogen concentrations (determined by the Kjeldahl method) were available for 83 profiles allowing the calculation of total N pool to a 100-cm depth. Carbon was analyzed by dry combustion (Matejovic, 1993). Clay (%) was analyzed by the hydrometer method (Gee and Bauder, 1986), pH was measured in 0.01 M CaCl2 (soil/extract 1:2.5), and available P extracted with 0.1 M H2SO4 for 2 h (Murphy and Riley, 1965). When analyzing organic C, carbonate C may constitute a problem, as combustion methods do not discriminate between inorganic and organic C. This problem has either resulted in the exclusion of carbonate containing horizons from analysis, or an overestimation of C. However, only five Alfisols and one Inceptisol included in this study contained carbonates (Calcium Carbonate in glacial till developed on limestone), and the very low C concentrations measured do not indicate carbonate C interference.
Soil profiles were classified to taxonomic order level according to Keys to Soil Taxonomy, using mainly the fourth edition (Soil Survey Staff 1992) but also earlier editions. Thirty-four profiles were classified as Spodosols, 31 Entisols, 23 Inceptisols, 16 Alfisols, and two Mollisols, which were merged with the Alfisols (18 total) based on similar pedogenesis and parent material. Conifers were growing on 28% of the Alfisols and more than 75% of the sites in other soil orders. Fagus sylvatica and Quercus robur were growing on sites with broadleaved species. Soil sampling procedures varied with site. Some samples were taken in the center of horizons; others (approximately 80%) were sampled from each horizon boundary to boundary. Detailed field sampling procedures and profile descriptions are obtainable in the original publications.
The soil horizons were given designations from FAO (FAO Soil resources, Management and Conservation Service, 1990) or Soil Survey Manual (Soil Survey Division Staff, 1993). To homogenize the data, horizon designations were standardized as follows: O horizons consisted of organic material on top of the mineral soil, A horizons are humus-rich mineral soil material in the top profile, E horizons are eluvial horizons, B horizons are illuvial horizons (accumulation of humus, clay, or sesquioxides), and C horizons denote the parent material. Standardization was based on profile descriptions and analytical data of each horizon.
Calculations: Soil Bulk Density, Carbon, and Nitrogen
Data on soil bulk density were available for 53 of the soil profiles. Carbon and N pools were calculated by horizon using concentration, horizon depth, and bulk density, and volumes were summed to a depth of 100 cm to express soil C and N contents in kilogram per square meter. Bulk densities (Db) for remaining profiles were calculated from models specifically developed for this paper by the authors. For this purpose, data on C concentration, bulk density, master horizon, and classification were supplied by L. Tau-Strand, Norway; M. Olsson, Sweden; C.J. Westman, Finland; and K. Raulund-Rasmussen, Denmark from a larger dataset of mineral forest soil horizons in four Nordic countries (N = 844). Plotting bulk density vs. C concentration suggested a curvilinear relationship, and hence the independent variable (C) was transformed to stabilize variances. The square-root function gave the highest degree of explanation in the following ANOVA, in which soil order and master horizon (A, E, B, and C), were used as class variables, and C concentration as a covariate (Eq. [1]).
![]() | [1] |
Soil order (P = 0.02), master horizon (P < 0.0001), and C (P < 0.0001) were significant, but some of the groups (by soil order and master horizon) were too small to yield good parameter estimates. The data were consequently pooled in two groups to secure a fair number of observations in each group: (i) Alfisols (loamy soils, N = 45 soil horizons), and (ii) Spodosols, Inceptisols, Entisols, and not classified soils (sandy soil texture, N = 799 soil horizons); thereby distinguishing loamy sandy soils and loamy soils from sandy soils. The model was reduced to Eq. [2] and analyzed for each dataset.
![]() | [2] |
Significant differences between master horizons were tested in Tukey-adjusted multiple comparisons using proc GLM (SAS Institute, 1997). The A, E, and B horizons in the larger dataset did not differ significantly from each other, and were consequently grouped together for parameter estimation (Table 1).
|
![]() | [3] |
The parameter estimates are presented in Table 1. The models were validated by an independent data set of 104 bulk density determinations in Danish forest soils (Vesterdal and Jørgensen, unpublished data, 1999). In this data set, the residuals (observed - predicted bulk densities) had a normal distribution with a mean of 0.03 g cm-3 and a standard deviation 0.18 g cm-3.
Out of a total of 106 O horizons used in calculation of O horizon: mineral soil C ratios (Table 2), nine had mass determinations based on areal sampling using 25 by 25 cm frames, and 83 were assigned a bulk density based on a study by Vesterdal and Raulund-Rasmussen (1998) of forest floor bulk density in Danish forest stands. These estimates were based on 49 determinations of seven tree species in monocultures including spruce, beech, and oak. Values used are: 0.12 g cm-3 for spruce, 0.06 g cm-3 for beech, and 0.04 g cm-3 for O horizons in oak stands. Fourteen O horizons in broadleaved stands (seven Alfisols, five Inceptisols, and two Entisols) had not been sampled at all mostly because of rapid decomposition. These profiles were given fixed C and N pools in the O horizon of 0.1 kg C m-2 and 0.005 kg N as a LOD (limit of detection), thus allowing us to compare the relative distribution of C in the organic layer versus mineral soil. Soils with no or a very limited organic layer could thus be compared with soils with a large O horizon.
|
|
= 5%) was determined using the Tukey option in Proc GLM (SAS Institute, 1997). Natural log-transformation was used to achieve variance homogeneity, since the residual plot showed increasing residuals with increasing predicted values (Weisberg, 1985). Estimates and 95% confidence limits were transformed back from log-scale to linear scale. Overall population means from the ANOVA (Table 2) were adjusted by MSE2/2 to correct the bias introduced when mean values of log-distributed variables are back-transformed to original scale (Parkin and Robinson, 1994). Estimates of C and N pools and C/N ratio presented in Tables 3 through 6 are back-transformed least squares means, and may be interpreted as median values on original scale. The correlation of C/N ratios with soil depth in master horizons was analyzed in a Spearman rank correlation analysis using the CORR procedure, (SAS Institute, 1990).
|
|
|
|
The distribution of total C and N between the O horizons and the mineral soil is expressed as the ratios CO/C0100 and NO/N0100 of a soil profile, where "O" designates the O horizon, and "0100" designates mineral soil in the entire volume 0 to 100 cm (Table 5). A ratio of 1.0 indicates an equal distribution between the O horizon and the mineral soil. A high ratio indicates the significance of the O horizon. A ratio of zero indicates a soil with no organic horizon. Clay (%) as a weighted average content (by depth of horizon) of the 50- to 100-cm section in mineral soil, P content (g m-2) of the 50- to 100-cm section, and the highest pH within the 100-cm soil depth (generally the C horizon) were used as indicators of soil fertility. Carbon, N, and C/N ratios in the O-horizon, the 0 to 100 cm of mineral soil, and the O horizon + mineral soil (0100 cm) were correlated with the fertility indicators using a Spearman rank correlation analysis in the CORR procedure (SAS Institute, 1990), as variables did not have normal distributions.
| RESULTS |
|---|
|
|
|---|
The total soil C content (O-horizon + 0100 cm) for the Spodosols was significantly greater (P < 0.0001) than in all other soil orders (Table 3); Alfisols had the least. In contrast, the highest N contents occurred in the Alfisols, and the lowest in the Spodosols, other orders were not significantly different. The Spodosol C/N ratio was the greatest; the Alfisol C/N ratio was the least.
The concentrations of C and N decreased downward in the profile from the O horizons (Table 4). The highest C and N concentrations were in the O horizons, whereas the lowest concentrations were found in the C horizons. Among the mineral soil master horizons, A horizons contained the largest concentrations of C and N (Table 4). The E and B master horizons (excluding Bh horizons) differed significantly in C concentration, while the N concentrations were not significantly different. The highest C/N ratios were found in the O and E horizons, followed by A horizons and C horizons. The lowest C/N ratios were found in the B horizons.
The concentrations of C in the B master horizons varied considerably; being highest in the horizons containing illuvial humus complexes; Bh and Bhs, and lowest in the Bt horizons (Table 4). The N concentrations in the five types of B horizons followed a different pattern. Again the highest concentrations were in the Bh and Bhs horizons, whereas the Bs, Bw, and Bt were all similar.
The distribution of C and N between the O horizons and the mineral soil (CO-org/C0100 and NO-org/N0100) showed that Alfisols contained significantly lower amounts of C and N in the O horizon (Fig 2). The greatest C and N content in the O horizon were in the Spodosols (approximately 0.3), although only for total N was significantly different from other soil orders.
The vertical distribution of C and N content with depth and master horizon is shown in Fig. 2. Horizontal lines representing the mean upper limit of the master horizon separate the horizons. Spodosols contained the most C in the O horizons and in the B horizon compared with other soil orders (P < 0.0001). Nitrogen contents in A horizons were significantly higher in Alfisols than in Spodosols and Entisols. Alfisols also had lower C/N ratios than Spodosols in both A and B horizons; Inceptisols and Entisols were intermediate, but significantly different from both Alfisols and Spodosols. Carbon/N ratios generally declined with the depth of the master horizon (Spearman rank correlation r - 0.26, P = 0.0001). Analyzed by soil order, the correlations were -0.06NS for Spodosols, -0.54*** for Alfisols, -0.39*** for Inceptisols, and -0.25** for Entisols (Fig. 2). Even though the overall trend was not significant, Spodosols had lower C/N ratios in A versus C horizons (P = 0.05). The negative correlation was because of higher C/N ratios in O horizon than B horizons in Alfisols (P = 0.003), and in Entisols (P = 0.001).
The soil clay content (%) ranged from 0 to 31% of the mineral soil. Clay was negatively correlated with O horizon C and N contents (P < 0.001), C/N in the O horizons (P < 0.001), and C/N ratio in the mineral soil (P < 0.001) (Table 6); and was positively correlated with N in the mineral soil (P < 0.001). Maximum mineral soil pH within the 100-cm soil depth was negatively correlated with C and N contents in the O horizons (P < 0.05), and the C/N ratio of mineral soils (P < 0.05), but positively with N content in the mineral soil (P < 0.01). Total P content of the 50- to 100-cm depth was negatively correlated with O horizon C (P < 0.001), N (P < 0.01), and C/N ratio (P < 0.01). Phosphorus content was negatively correlated with mineral soil C (P < 0.01) and was positively correlated with mineral soil N (P < 0.05) (Table 6).
| DISCUSSION |
|---|
|
|
|---|
Estimates of total N through the soil profile exist in the literature, but estimates of soil N including N in both O horizon and mineral soil to depth are rare. Carter et al. (1998) found 0.36 to 1.05 kg N m-2 in Canadian agricultural soils. Compared with these values, the Danish forest soils had intermediate N contents. Zinke and Stangenberger (2000) reported median values of 0.61 kg N m-2 to a 1-m depth in dense conifer forest, and 0.27 kg N m-2 in low density forest on the lower western slopes of the Sierra Nevada, CA. It was not stated explicitly whether those values included the O-horizon. Assuming exclusion of the O-horizon these values bracket our estimate of 0.5 kg N m-2 for mineral soil to a 1-m depth (Table 2). Batjes (1996) calculated soil C and N storage in 1 m of mineral soil using worldwide databases, and found rather high soil N contents, for example, 1.39 kg N m-2 for Podzols (Spodosols), 1.03 kg N m-2 for Luvisols (Alfisols), and 0.52 kg N m-2 in Arenosols (Entisols). Other soil orders also had higher mean values for soil N than our data. However, comparison is not straight forward for three reasons: (i) the means reported in Batjes (1996) are probably skewed to the right (Zinke and Stangenberger, 2000), and (ii) the distribution of land uses was not reported, but agricultural soils are probably a prevalent land use in the data, and (iii) the global data reflect the distribution of global climate zones, whereas our estimates pertain to a cool, humid temperate climate. Eriksson and Rosén (1994) estimated total soil N to a 95-cm soil depth (including O-horizon) to be 0.5 to 0.65 kg N m-2 in a tree species trial with conifers in southwestern Sweden. The climate is similar to the Danish climate, and the contents of soil N appear quite similar.
The only deviation from the general pattern of gradually decreasing C and N concentrations with depth occurred in soils with a spodic horizon (Bh, Bhs). The pattern in C and N concentration was O > A > E < Bh, Bhs > C (Table 4). The C concentrations of spodic Bh and Bhs horizons (Table 4) indicate that the podzolization process is responsible for the accumulation of C, and expressed by contents (kg m-2) this is also illustrated in Fig. 2. Accordingly, our analysis showed that podzolized soils contain more C than the other soil orders. Entisol profiles, containing large quantities of C, usually had a horizon with spodic properties, which, because of insufficient thickness, did not meet the formal requirements of a spodic horizon (Soil Survey Staff, 1992). Normally, soils that do not meet the requirements for inclusion in a major soil order end up in the Inceptisol order. However, the widespread Danish coarse-sandy textured soils do not comply with the requirements for a cambic horizon; and they lack the dark color required for an umbric epipedon, since the upper epipedon is often dominated by albic material. Therefore, many podzolized sandy soils of Denmark are classified as Entisols, unless they possess a spodic placic horizon in the Bh or Bhs horizon. Westman et al. (1994) also attributed the podzolization process per se to the immobilization of significant amounts of C (4.011.9 kg m-2) in a study of forest soils in Finland. The cooler climatic conditions and soil ages of >6000 yr in Finland are probably the reasons for higher Spodosol C pools of this study. Hence, in Spodosols and Spodic Udipsamments, the quantity of organic C is probably predictable from information on soil order and suborder. Other studies on Spodosol C pools found C stores of 10.4 kg C m-2 (Stone et al., 1993), 17.4 to 29.8 kg C m-2 (Eswaran et al., 1993), and 9.0 to 11.8 kg C m-2 (Martel and Deshenes, 1976). Gregorich et al. (1995) also found more C in podzolic soils (Spodosols) than in luvisolic soils (Alfisols). The differences are because of variations in the soil forming factors, not least being soil age. Extended periods of time allows for creation of deep spodic horizons. Climate and hence primary production play a major role too.
In contrast to C stores, N stores in this study were highest in Alfisols. Alfisols are generally base rich and have characteristically high rates of decomposition because of high litter quality, that is, high nutrient concentration content (Vesterdal, 1999), and a more favorable environment for decomposers, for example, more stable soil moisture content because of a higher water storage capacity. The Spodosols possess the lowest N stores because of the lower quality of litter (Vesterdal, 1999), and in combination with the huge quantities of C, it explains the higher C/N ratios of this soil order.
The ratios CO-hor/C0100 and NO-hor/N0100 demonstrate the importance of the O horizons as forest ecosystem C and N pools (Table 5). In Spodosols and Entisols, the ratios for C were 0.31 and 0.28, respectively, and for N the ratios were 0.30 and 0.17. The variation in the C and N contents of O horizons was because of variation in thickness and mineral matter content of the O horizon. The O horizon contained from <1 to 150% (data not shown) of total soil C and N present in depth 0 to 100 cm. Grigal and Ohmann (1992) found 55% of total ecosystem C (including the living biomass) in the mineral soil, and 9% in the forest floor, that is, a CO-hor/C0100 ratio of 0.16, which is similar to the Inceptisols in our study. McFee and Stone (1965) found 43% of the N content and 38% of the organic matter content of the total system, including the standing crop, in the O horizon of forest soils. These results show that O horizons should be included in any estimation on soil C and N pools, though this compartment is often underestimated or ignored in pedological studies (Eswaran et al., 1995).
The narrowing C/N ratio downwards in the profile from O to A and B horizons (excluding Bh horizons) may be because of the decomposition stage and age of humus in horizons (Lal et al., 1995). However, E and C horizon C/N ratios were not significantly different from O horizon C/N ratios. Recently shed litter has a high C/N ratio, this ratio decreases during decomposition. As the age of organic matter increases with soil depth (Lal et al., 1995), the C/N ratio consequently decreases. The decreasing C/N ratios with depth and age of organic matter may be attributed to microbial immobilization of N and loss of C through respiration (Aber and Melillo, 1980; Staaf and Berg, 1982).
The inverse relationship between clay content and C content, N content, and C/N ratio in the O horizons found in this study is a result of the O horizons associated with Spodosols. Interestingly, the total amount of C (O horizon + 0100 cm) is also negatively related to clay content. This suggests that fertile Alfisols with humus rich A horizons are surpassed by the sandy nutrient-poor Spodosols in terms of C sequestration. Spodosols have two C sinks; the O horizon and the Bh/Bhs horizon. These results are in contradiction to the findings of Grigal and Ohmann (1992), who concluded that C content was positively related to clay content, probably by its influence on soil moisture and fertility. Several other studies suggested that soil C and soil clay is positively related to clay content (Oades, 1988; Burke et al., 1989), or not related (Carter et al., 1998).
The same sort of relationship is demonstrated by using the extractable P content (g m-2) in depth 50 to 100 cm, as an indicator of site fertility. The C and N contents, the C/N ratio of the O horizons, and the total C content of the soil profile are negatively related to P content. On the other hand, the fertile soils contain the highest amounts of N as demonstrated by the weakly positive relationship between mineral soil N content and P content. The chain of mechanisms behind this may be that fertile soils have high productivity, but also high decomposition rates, resulting in low C contents. However, the fertile sites have relatively high N contents, possibly because humus with low C/N ratios is protected in organo-mineral complexes as suggested by Christensen (2000).
A concern in many soil C and N studies that rely on data from previously conducted studies is that data may contain errors because of differences in the sampling and analytical methods of the original studies (Eswaran et al., 1995). Estimations of C and N pools in this study were potentially compromised by measurement errors in bulk density, soil volume calculations, C and N determinations, and stoniness of soils. A lack of bulk density data is a common problem for estimation of global soil C pools (Carter et al., 1998), and our predictions of missing bulk densities may add to the error. Murillo (1994) estimated the compounded error of soil profile organic C estimates to be 30%, which is a realistic error estimate for our data. The sources of error might introduce a bias, when comparing soil C and N content among studies in the literature where different methodologies have been applied. However, none of these sources of error introduce a bias among soil orders, and with a high number of profiles the error is diminished, thus justifying the use of pedological soil databases for ecological purposes. However, tree species and soil order are confounded. This affects the O-horizon estimates, but it nevertheless reflects the distribution of tree species in Denmark.
When interpreting data from natural systems, the multidimensional features of these systems should be taken into account. Grigal and Ohmann (1992) used a multiple regression technique; 65% of the variation in total forest C was attributed to forest type, stand age, available water capacity, actual evapotranspiration (AET), and soil clay content. Our results suggest that horizon designation, horizon thickness, and soil order are predictors of total soil C and N content; all based on pedology. Percentage of clay, although correlated to C content, is not a predictor of soil development, or the degree of podzolization on sandy soils. The P content of the subsoil (the 50- to 100-cm stratum) seemed to be a suitable indicator, but such data are not generally available in soil surveys. Soil pH did not prove an adequate variable for prediction of soil C and N content. Our study clearly suggests that soil order may predict the levels of C and N. At least a distinction in C and N content between the two major Danish soil orders, Alfisols and Spodosols, is possible under the given climatic conditions. For soils possessing spodic properties, sound estimates on C and N content should also be possible.
If soil C sequestration becomes a primary goal of forest management, afforestation in Denmark should take place on Spodosols or soils with spodic properties. Podzolization mechanisms favor accumulation of large O horizons and the formation of humus-rich spodic horizons. This was also proposed by Cannell and Milne (1995) who concluded that the equilibrium storage of C is greatest when litter decompose slowly and a large fraction is transferred to recalcitrant soil organic matter pools. Our study suggests that C storage in Danish forest soils is most likely because of impeded decomposition in nutrient poor soils, and not driven by high productivity at nutrient-rich soils.
Received for publication June 26, 2001.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
Y. L. Zinn, R. Lal, J. M. Bigham, and D. V. S. Resck Edaphic Controls on Soil Organic Carbon Retention in the Brazilian Cerrado: Texture and Mineralogy Soil Sci. Soc. Am. J., June 8, 2007; 71(4): 1204 - 1214. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Natural Resources and Life Sciences Education |
Vadose Zone Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||