Published in Soil Sci. Soc. Am. J. 68:895-903 (2004).
© 2004 Soil Science Society of America
677 S. Segoe Rd., Madison, WI 53711 USA
DIVISION S-5PEDOLOGY
Spatial Distribution of Soil Carbon in Southern New England Hardwood Forest Landscapes
Aletta A. Davisa,
Mark H. Stolt*,b and
Jana E. Comptonc
a U.S. Forest Service, Center for Forested Wetlands Research, Charleston, SC
b Dep. of Natural Resources Science, University of Rhode Island, Kingston, RI
c USEPA, NHEERL, Western Ecology Division, Corvallis, OR
* Corresponding author (mstolt{at}uri.edu).
 |
ABSTRACT
|
|---|
Understanding soil organic C (SOC) spatial variability is critical when developing C budgets, explaining the cause and effects of climate change, and for basic ecosystem characterization. We investigated delineations of four soil series to elucidate the factors that affect the size, distribution, and variability of SOC pools from horizon to landscape scales. These soils, classified as Udipsamments, Dystrudepts, Endoaquepts, and Haplosaprists, were sampled along random transects to a depth of 1 m. In very poorly and poorly drained soils 75 and 45% of total SOC was found below 30 cm, respectively. In contrast, only 30% of the total SOC could be accounted for below 30 cm in the well and excessively drained soils. Soils formed in outwash and young alluvium sequestered a greater portion of total SOC within the subsoil, while soils formed in loess held approximately 70% of the SOC within O and A horizons. Total SOC contents among the four soil types differed significantly (p < 0.001), with the wetter soils having greater accumulations of C. Soil C pools ranged from 110 Mg C ha1 in the excessively drained Psamments (double the mean national value) to 586 Mg C ha1 in the very poorly drained Saprists (3060% lower than the mean national value). The two-fold differences between our data and the national averages support the need for regional assessments of soil C pools. Based on the coefficient of variation (CV) values, there appears to be nearly as much variability in the SOC pool within a delineation (CVs ranged 9 to 30%) as among delineations (CVs ranged from 15 to 31%) for the same soil type. Since significant differences were found for total SOC among delineations of the same soil type, we concluded that sampling from a significant number of delineations of the same series will provide a more accurate representation of SOC for scaling to the landscape or region than sampling at multiple locations within a single representative delineation.
Abbreviations: CV, coefficient of variation SOC, soil organic C
 |
INTRODUCTION
|
|---|
ASSESSMENTS OF THE DISTRIBUTION of C within and among soil types are critical to developing an understanding of the cause and effect relationships between climate or land use change, and release of CO2 to the atmosphere (Schimel et al., 1994). Carbon distribution data are available from many sources at a variety of scales. Pedon data collected as a part of the USDA-NRCS Progressive Soil Survey have been used to estimate SOC pools in a number of studies (Franzmeier et al., 1985; Davidson and Lefebvre, 1993; Eswaran et al., 1993; Kern, 1994; Homann et al., 1995, 1998; Grossman et al., 1998). Similar conclusions and recommendations were drawn from each of these works and included: (i) taxonomic grouping is the best approach to estimate soil C pools; (ii) lower taxonomic categories (more detail) in Soil Taxonomy are more reliable predictors of C pools than are higher categories; and (iii) there is a strong need for estimates of C variability within and among pedons of the same soil type.
Although the national soil survey database is quite comprehensive, additional information would increase the usefulness of the database for predicting regional or national C pools. Of particular use would be data on the contribution of O horizons and subsoils of forest soils to the C pool (Eswaran et al., 1993; Stone et al., 1993; Grossman et al., 1998; Lal et al., 1998). In addition, a measurement of spatial variability in soil properties among taxonomically similar soils that are important in the calculation of C pools (i.e., percentage of C, rock fragment content, and bulk density), and an understanding of the confounding effects of climate, land use, species age and composition, and parent material on soil C pools would aid in the interpretation of the C estimates (Eswaran et al., 1993; Kern, 1994; Lathrop et al., 1995; Bouwman and Leemans, 1995; Batjes, 1996).
In this study, we examined the factors that influence soil C pools at the landscape scale. Landscape unit boundaries were defined based on soil survey delineations of four soil series typically found in southern New England. The objectives were to assess the size, variability, and distribution of the soil C pool in hardwood forests of this region. We addressed the following questions: (1) How does soil C content vary among soils with similar forest classification (hardwood forest) but varying soil properties? (2) Is there more variability in total soil C content within a delineation than among mapping units of the same series? (3) Is variability systematically related to soil characteristics such as drainage class or parent material?
 |
MATERIALS AND METHODS
|
|---|
Study Area
This research was focused on forested areas within the Pawcatuck River watershed in southwestern Rhode Island (Fig. 1)
. This 64000 ha watershed is dominated by glacial landforms composed of till and outwash materials of late Pleistocene age that are representative of southern New England (Schafer and Hartshorn, 1965; Quinn, 1971). Outwash covers most of the valley landscape in which the present day streams have shifted or meandered across depositing alluvial materials. In addition, portions of outwash landscapes are covered with a silt mantle (Flint, 1930; Schafer, 1981; Rector, 1981; Lawson, 1995), considered to be loess (Wright and Sautter, 1988; Stolt, 1998).
Potential research sites were identified by overlaying geographic information systems (GIS) land use/land cover and soils coverages (August et al., 1995) across the Pawcatuck River watershed to identify mapping units of soils formed in outwash, aeolian, alluvial, and organic parent materials that were dominated by hardwood forests. Soil GIS coverages were developed from the 1:15840 soil survey of Rhode Island (Rector, 1981). Over 70 delineations (mapping units), identified as Windsor, Enfield, Raypol, and Carlisle (see Table 1 for family classifications) and mapped as consociations, were examined in a field reconnaissance survey. These four soil types were chosen because they are commonly mapped in southern New England and provide a diverse range of drainage classes and parent materials (Table 1).
Three delineations representative of the Enfield and Carlisle series, and two units of the Raypol and Windsor series were chosen for detailed study. These sites were chosen because the soils within each delineation met the range of characteristics established for each soil series, the landscapes were dominated by hardwood forests, and the age of the trees was relatively uniform within the delineation. The largest trees (based on diameter at breast height) within 25 m of the soil transects (see below) of each of the dominant species were cored to determine the age of the stand. A minimum of 10 increment cores was collected from trees at each site and stand age was determined by averaging the four oldest trees of those selected (Table 1). We avoided unusually large trees because they might represent a period before the majority of the stand developed.
Sampling Design
Soils were collected along two randomly chosen transects within each of the 10 delineations (Young et al., 1991). Sampling locations were spaced at 25-m intervals for a total of five sampling points (0, 25, 50, 75, and 100 m) per transect (10 per delineation). Samples at each location were collected by horizon to a depth of 1 m (except in a few cases). Organic horizons of the mineral series were collected from a 15 by 15 cm area intact and separated into fibric (Oi), hemic (Oe), and sapric (Oa) horizons when possible. Mineral samples were collected using a 4.7-cm diam. split-core sampler. Fibric and hemic horizons of the organic series (Carlisle) were collected from a 15 by 15 cm area intact and separated into the two horizon materials. Sapric materials were collected by horizon using a Macauley peat sampler.
All samples were analyzed for bulk density and organic matter content. Bulk density calculations were corrected for coarse fragment (>2 mm) mass and volume (Blake and Hartge, 1986). Organic matter contents were determined for all of the soil samples via loss on ignition (Nelson and Sommers, 1996). Samples from four of the mapping units (one for each series) were analyzed for total C using an automated C and N analyzer (Carlo-Erba, Milan, Italy) (Nelson and Sommers, 1996). These acid soils are free of carbonates. Organic C data from these four mapping units (550 samples) were used to establish horizon-specific and soil-series specific C/OM regression equations, which were used to determine C concentrations in samples from the remaining six units (Davis, 2001). To calculate total SOC (Mg C ha1), C data were pooled by Oi and Oe (Oie), Oa, A, B, and C horizons in the mineral soils. In the Carlisle soils, data were pooled by Oie horizons, and upper (050 cm) and lower (50100 cm) Oa horizons. When calculating means of the pooled data we weighted the subhorizon data that collectively forms the pooled master horizon by the thickness of each subhorizon. In certain pedons particularly loose or dense materials, or an abundance of coarse fragments, prevented us from sampling the C horizons. We estimated the C horizon portion of the total SOC stored in these pedons by using the average bulk density and percentage of C values from the C horizons of the pedons within the mapping unit that could be sampled.
Descriptive statistics such as means and coefficients of variation (standard deviation/mean) and comparative (ANOVA, and Tukey's multiple comparison) analyses were performed in SPSS version 10 (SPSS Inc., 1999). A nested ANOVA (with mapping unit effects nested within soil type) was used to compare total SOC contents in the upper 1 m of the soils. Four one-way ANOVAs were run to determine if there were significant differences in total SOC content among delineations of the same series. The Bonferonni adjustment was used to correct p-values for the increased probability of Type I errors associated with comparing the results of multiple comparisons.
 |
RESULTS AND DISCUSSION
|
|---|
Site, Pedon, and Horizon Characteristics
Acer rubrum and mixed Quercus species were the dominant vegetation at the study sites (Table 1). Silt loam, deposited as loess or alluvium, was the predominant texture of the Raypol and Enfield sola. In these soils, sandy or sandy skeletal outwash occurred at the base of the solum such that these soils had strongly contrasting family particle-size classes (Table 1). Windsor soils developed in outwash and were dominated by loamy sand and sand textures.
The average solum thickness of the mineral soils ranged from 66 to 88 cm (Table 2). The combined Oi and Oe horizons (Oie) ranged in thickness from 2 to 11 cm and had similar bulk density values (0.070.12 g cm3) and SOC concentrations (400500 g kg1) regardless of the soil or site (Table 2). Windsor soils lacked Oa horizons and less than half of the Enfield pedons (13 of 29) contained Oa horizons. The Raypol soils had the thickest Oa horizons of the mineral soils (mean: 9 cm). Carlisle soils were dominated by thick (>1 m) accumulations of sapric (Oa) materials. The SOC concentrations of the Carlisle Oa horizons were significantly higher, and the bulk density values significantly lower, than the Oa horizons of the mineral soils (Table 3). The mineral soil Oa horizons are thin by nature and there is a greater likelihood that denser mineral material will be mixed in with the organic materials as a result of pedoturbation. Therefore, these horizons tend to have higher bulk density and the lower SOC concentrations. Overall, the wetter soils had greater SOC storage in O horizons, and this difference is a function of the C contents in Oa horizons. By comparison, C contents in O horizons in Oregon (10 Mg C ha1, Homann et al., 1995), the Great Lakes Region (17 Mg C ha1, Grigal and Ohmann, 1992), and on average for temperate forests (21 Mg C ha1, Vogt et al., 1995) are lower than the average C contents in the O horizons of the mineral soils we studied (2142 Mg C ha1).
View this table:
[in this window]
[in a new window]
|
Table 2. Weighted means (Mean) and coefficients of variation (CV) for thickness, bulk density, soil organic C (SOC) content, and total SOC storage for the master horizons, sola, and upper 1 m of the 10 study sites. Means were weighted based on the thickness of the subhorizons that together form the master horizon.
|
|
View this table:
[in this window]
[in a new window]
|
Table 3. Average weighted means (Mean) and coefficients of variation (CV) for thickness, bulk density, SOC content, and total soil organic C (SOC) storage for the master horizons, sola, and upper 1 m of the four soil types. Data represent an average of the data from each site within a soil type.
|
|
Soil organic C concentrations for A and B horizons ranged from 25 to 63 and 3 to 14 g kg1, respectively (Table 2). Carbon concentrations in the A horizons increased as drainage class became wetter such that the Enfield (well) and Raypol (poorly drained) A horizons had significantly higher SOC levels than the excessively drained Windsor (Table 3). Raypol soils had significantly higher SOC concentrations in the B horizons, suggesting that drainage class affects SOC content throughout the solum. The higher SOC concentrations in the Raypol B horizons resulted in significantly lower bulk density values compared with B horizons in the Windsor soils. The thickness of the A and B horizons or the solum the three mineral soils were not significantly different, however, the total SOC held in these horizons and the solum increased as the drainage class became wetter.
Total C storage is a function of horizon thickness, bulk density, and SOC concentration. Of these parameters, the one that showed the most variability depended upon the soil, horizon, and whether the comparisons were made within sites (Table 2) or among sites of the same series (Table 3). For example, CV for horizon thickness for the Oie horizons ranged from 3 to 25% within individual sites (Table 2) and between 19 and 79% when means were calculated based on all of the data for each series (Table 3). By comparison, CVs for SOC concentration of the Oie horizons were all <20% and for means calculated within a site the CVs were all <10%. Much greater variability (CVs of 2757%) was reported for SOC concentrations in the top 20 cm of forests soils in Oregon (Homann et al., 1995) and Wisconsin (Kabrick et al., 1997), respectively, suggesting that a number of factors such as soil type, forest type, management history, stand age, or climate are also important in explaining SOC variability. Average CVs for the parameters used to calculate total SOC for a horizon or profile were highest (4155%) for the poorly drained Raypol soils (Table 3). This variability is a reflection of the two types of parent materials of these soils (alluvium and loess). Loess tends to have minimal variability associated with the depositional process (Mausbach et al., 1980). In contrast, alluvial materials are deposited as layers in response to repeated flooding which can result in considerable variability with respect to depth-related properties.
Distribution of Soil Organic Carbon within Profiles
We examined the average proportion of total SOC stored within each master horizon among delineations and soil parent materials (Fig. 2 and 3)
. The O horizons in mineral soils accounted for 20 to 40% of the total soil C pool. Based on the ages of the forests, the O horizons of these mineral soils are probably only 60 to 100 yr old. We did not track the land-use history of these sites. However, based on soil properties we can suggest some likely uses. All of the mineral soils contained Ap horizons indicating that these soils were used for agriculture purposes in the past. Enfield and Windsor soils have minimal amounts of coarse fragments, have good drainage throughout the solum, and are presently used for crops and turfgrass in southern New England. Raypol soils were also cultivated at some point, however, because of periodic saturation in the upper solum and possible flooding, these sites were likely used more often for pasture than cultivation. Therefore, our C storage data suggest that within the first 60 to 100 yr of forest regrowth, following a land use change from agriculture, the development of O horizons may lead to a substantial increase in SOC.

View larger version (46K):
[in this window]
[in a new window]
|
Fig. 2. Proportion of the total soil organic C (SOC) pool stored in upper 1 m by master horizon of the three mineral soil sites. Windsor soils have formed in outwash and are excessively drained; Enfield soils have formed in loess over outwash and are well drained; and Raypol soils are poorly drained having formed in loess (Rt2) and alluvium (Rt138).
|
|

View larger version (35K):
[in this window]
[in a new window]
|
Fig. 3. Proportion of the total soil organic C (SOC) pool stored in the upper 1 m by master horizon based on soil parent materials.
|
|
The A horizons generally contained a similar proportion of the total SOC as the O horizons (Fig. 2 and 3). Only at the drier Windsor Bg site does the proportion of SOC contained in the A horizon appear considerably less than in the O horizons. Subsurface horizons (B and C horizons) contained between 26 and 55% of the SOC within the mineral soils (Fig. 2 and 3). Nearly 40% of the total SOC stored in the Windsor soils was stored in the B horizons. Torn et al. (1997) found similar amounts of SOC (3370% of the total pool) in subsurface horizons. These studies emphasize the importance of including SOC in the subsoil horizons to the total SOC.
Drainage class influences the distribution of SOC throughout the profile, with wetter soils storing a greater portion of SOC deeper in the profile (Zdruli et al., 1995; Rapalee et al., 1998). The effect of drainage class on the distribution of C with depth can be clearly observed when comparing the excessively well-drained Windsor and well-drained Enfield soils to the poorly drained Raypol and very poorly drained Carlisle soils (Fig. 4) . Over 75% of the C stored in the upper meter of the Carlisle soil is below 30 cm. In comparison, Raypol soils store about 45% of the total SOC pool below this depth, and in the better drained Windsor and Enfield soils approximately 30% of the total SOC is present below 30 cm. The similarity in the distributions of total SOC with depth for the Windsor and Enfield soils (Fig. 4) contrasts to how SOC is proportioned among horizons (Fig. 2). Windsor soils (Psamments) store the largest proportion of the SOC pool in the B horizons (Fig. 2), whereas in the Enfield soils (coarse-silty over sandy) the A or O horizons contain the largest proportion of the SOC pool. These data suggest that parent material affects the distribution of SOC among horizons and that the effects of drainage class on the proportions of the SOC pool with depth are more important on the wetter end of the spectrum.

View larger version (24K):
[in this window]
[in a new window]
|
Fig. 4. Cumulative proportion of the total soil organic C (SOC) content with depth for upper 1 m of the four soil types studied. Data were averaged at 5-cm intervals for all of the pedons for each soil type.
|
|
Total Soil Organic Carbon Content Within and Among Soil Types
Total SOC content, defined on the upper 1 m of the soil including the organic horizons, ranged from 104 to 223 Mg C ha1 in the mineral soils and between 494 and 645 Mg C ha1 in the organic Carlisle soils (Tables 2 and 3). Analysis of variance showed a significant difference in total SOC among both soil types and sites (Table 4). The size of the total SOC pool followed a drainage class trend with Carlisle soils containing significantly greater (p < 0.001) amounts of total SOC than the three mineral soils and the Raypol soils also having significantly higher total SOC than the better drained Enfield and Windsor soils (Table 3). The SOC pools reported for forest soils in New Hampshire (160 Mg C ha1; Huntington et al., 1988) fall within the range of SOC in the mineral soils of this study. Total SOC in the Enfield series (136 Mg C ha1) was similar to SOC content in well-drained till-derived soils (115130 Mg C ha1) in Rhode Island (Hooker and Compton, 2003). Soil C content in the Carlisle series was approximately 20% lower than estimated pools of SOC in the upper meter of organic soils in the north central USA (740 Mg C ha1; Franzmeier et al., 1985) and the estimated global mean for wetlands (720 Mg C ha1; United States Department of Energy, 1999).
View this table:
[in this window]
[in a new window]
|
Table 4. Summary of nested analysis of variance with site nested within soil type for total soil organic C (SOC) (Mg ha1) within the upper meter.
|
|
There was no significant difference at the 0.05 probability level in SOC content between the two Windsor mapping units (Table 5). Total SOC stored varied significantly within mapping units of the Enfield (p = 0.003), Raypol (p = 0.001), and Carlisle (p = 0.004) soils, suggesting that soil and site conditions, such as basal area or geomorphic setting affects the amount of SOC in these soils. In Enfield sites, the differences in SOC content among sites could be attributed to lower SOC concentrations (3% lower) in the Oie horizon and the absence of an Oa horizon in the GS2 site, and lower SOC concentrations (ranging from 0.55% lower) throughout the Oie and Oa horizons of the Rt95 site (Table 2). Thus, seemingly minor differences in O horizon characteristics can have a significant effect on SOC pools at the landscape scale. Each of the forest stands occupying the Enfield landscapes are approximately 80 to 85 yr old (Table 1), thus stand age cannot explain differences in SOC storage in the Enfield soils. Basal area for the GS1 site, which has significantly higher SOC content than the other two Enfield sites, is nearly 60% greater (10 m2 ha1) than the GS2 site. A greater basal area may equate to higher yearly C additions to the soil through litterfall and root turnover, and thus a larger C pool in the O and A horizons.
View this table:
[in this window]
[in a new window]
|
Table 5. Summary of one-way analysis of variance comparing total SOC (Mg ha1) in the upper meter means within soil types.
|
|
The large difference in SOC content between the two Raypol units was attributed to a thinner O horizon and lower SOC concentrations throughout the solum of the alluvial Raypol site (Rt138; Table 2). The Rt138 site is located on a floodplain geomorphic setting and by far has the lowest SOC content within the Oie horizons of the mineral soils. Soils on floodplains are not as stable as the upland areas, and flooding events may remove some of the lighter organic fractions while depositing organic matterpoor silty mineral material, resulting in a lower SOC pool. The Rt2 site in the present land use appears to be nearly twice as old as the Rt138 site. The additional 50 yr of continued forested conditions may account for some of the additional C in the soil pool.
The reasons for the 150 Mg C ha1 difference among the Carlisle sites are not apparent. Trees at the Rt110 site are the youngest of the three sites and the Rt110 site has the largest SOC pool (Table 2), suggesting that stand age is not a factor. Likewise, even though the Carlisle mapping unit with the lowest SOC pool also had the smallest basal area (18.6 m2 ha1), this trend did not follow for the other two sites, suggesting that differences in basal area could not be used to account for variability in SOC content. Observations of change in the peat morphology with depth suggest that over time a number of different plant communities may have occupied these landscapes. For example, at the University of Rhode Island site reddish brown (2.5YR 3/3) Oa horizons occur interdispersed within the more typical black sapric materials. Reddish brown sapric materials in southern New England wetlands are indicative of vegetation dominated by white cedar (Chamaecyparis thyoides). These differences in plant communities likely reflect changes in hydrologic conditions, which may affect C concentrations and total SOC contents.
Measures of variability for SOC (Mg ha1) content in the solum and upper meter of these soils were similar or lower than for individual horizons (Tables 2 and 3). Variation in SOC content in the upper 1 m within individual mapping units ranged from 9 to 30% (Table 2). Average measures of variability for the four soil types were similar to the individual sites and ranged from 15 to 31% (Table 3). Windsor soils had the lowest average CV for SOC content in the upper 1 m (15%), followed by the Carlisle (20%), with the Enfield and Raypol having similar amounts of variability (29 and 31%, respectively). Huntington et al. (1988) reported slightly higher CVs (38%) for SOC content within mixed hardwood forest soils in New Hampshire. Lower within-site variability in SOC content (CVs of 016%) was observed in forest soils throughout the Northwest (Homann et al., 1995).
Comparison with National and Global Mean Values
A number of studies at regional, national, and global scales have used soil taxonomic groupings and soil survey data to estimate area-wide SOC pools (Table 6). Total SOC content within the Windsor units (110 Mg C ha1) was approximately double the global and national values predicted for soils with the same classification (soil order, suborder, and great group), but were similar to estimated SOC contents in Regosols (Entisols) in Canada. Soil C content in the Enfield series was slightly lower than mean values estimated for Udepts (Ochrepts) in Maine (Davidson and Lefebvre, 1993) and Dystrudepts in the Pacific Northwest (Kern et al., 1998). Both the Enfield and Raypol total SOC estimates were similar to global values for their corresponding World Taxonomy suborders, Dystric and Gleyic Cambisols, respectively. Raypol soils had higher SOC contents than would have been predicted based on their USDA taxonomic classification. Soil C content in the Carlisle series (586 Mg C ha1) was 30 to 60% lower than national and global values for Histosols. Similar disagreements between national averages and regional estimates were reported for Spodosols in Maine (Davidson and Lefebvre, 1993) and Florida (Stone et al., 1993). These findings suggest national and global aggregations of soils data, based on soil taxonomic classification, may not accurately reflect regional or local SOC pools.
View this table:
[in this window]
[in a new window]
|
Table 6. Total soil organic C (SOC) (Mg ha1) in the upper 1 m of the Windsor, Enfield, Raypol, and Carlisle soils compared with total SOC values from the literature for soils of similar taxonomic classification. Locations in parentheses under the Criteria heading indicate whether data are reported as local, regional, national, or global averages.
|
|
Implications for Scaling-up Carbon Estimates
At the landscape scale we found that C pools varied from 100 to 650 Mg ha1. Considering that the vegetation was hardwood forest at each site, and that there were significant differences in SOC pools among the various soil types, the use of coarse vegetation classes to group landscapes is not adequate for estimating SOC pools. Therefore, any comprehensive approach to inventory SOC pools will likely involve national soil survey data. How these soils data are to be grouped remains in question. Our comparisons with global and national averages suggest that grouping by taxonomic order is not an acceptable approach. For example, the Histosols that we examined stored 30 to 60% less SOC than the global or national averages. Similar discrepancies were observed for the mineral soils. One of the core problems with using the order approach is exemplified by the 50 Mg ha1 (or 2736%) difference in SOC between the poorly drained (187 Mg ha1) and well-drained (136 Mg ha1) Inceptisols. To reduce the effect of this coarse-scaling problem, Kern (1994) suggested that SOC data be averaged by great group because at this level specific information regarding soil properties believed to affect C storage such as drainage, climate, and coarse textures are incorporated in the classification. Region-wide, the great group approach may be useful in grouping soils, but even at this level we still found substantial differences between national averages and our SOC pools for the same great groups.
Soils examined in this study were chosen based on series level criteria since this level of classification provides the most detailed information about a soil. Most map units within second-order (county-wide) soil surveys are consociations. With the digitizing of most of the soil surveys in process or completed, and the additional SSURGO certification of many surveys in the works, estimating C pools using a GIS with series level attribute pedon data appears to be a possibility. The spatial variability data collected in this study provide some insight into how the soil survey data can be best used to estimate SOC pools within a region. Three of the four soil series we examined were found to have significantly different pools of C among delineations of the same series. These differences suggest that data from a single typical pedon may not provide the true mean SOC pool stored in these soils. We sampled 10 pedons within each of two or three delineations for each series. Based on the CV values, there appears to be nearly as much variability in the SOC pool within (CV's ranged 930%) as among delineations (CV's ranged from 1531%) for the same land cover and soil type. These spatial variability data suggest that, instead of sampling at multiple locations within a single delineation to obtain the average SOC for a map unit, a more useful approach may be to sample from a significant number of delineations of the same series. This map unit SOC pool could then be applied to each of the map unit delineations in the soil survey area, recognizing that land use differences may result in additional variability to the overall C content estimate.
Received for publication March 4, 2003.
 |
REFERENCES
|
|---|
- August, P.V., A.J. McCann, and C.L. LaBash. 1995. Geographic Information System in Rhode Island. Natural Resources Facts, Factsheet 951. Univ. Rhode Island. Kingston.
- Batjes, N. 1996. Total carbon and nitrogen in the soils of the world. Eur. J. Soil Sci. 47:151163.
- Blake, G.R., and K. Hartge. 1986. Bulk density. p. 363375. In A. Klute (ed.) Methods of soil analysis. Part 1. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Bouwman, A.E., and R. Leemans. 1995. The role of forest soils in the global carbon cycle. p. 503526. In W.W. McFee and J.M. Kelly (ed.) Carbon forms and functions in forest soils. SSSA, Madison, WI.
- Davidson, E.A., and P.A. Lefebvre. 1993. Estimating regional carbon stocks and spatially covarying edaphic factors using soil maps at three scales. Biogeochemistry 22:107131.
- Davis, A.A. 2001. Carbon storage, variability, and flux in forest soils of southern New England. Master Thesis. University of Rhode Island, Kingston.
- Eswaran, H., E. Van Den Berg, and P. Reich. 1993. Organic carbon in soils of the world. Soil Sci. Soc. Am. J. 57:192194.
- Flint, R.F. 1930. The glacial geology of Connecticut: Connecticut geological and natural history survey bulletin 47.Connecticut Geological and Natural History Survey, Hartford, CT.
- Franzmeier, D.P., G.D. Lemme, and R.J. Miles. 1985. Organic carbon in soils of North Central United States. Soil Sci. Soc. Am. J. 49:702708.[Abstract/Free Full Text]
- Grigal, D.F., and L.F. Ohmann. 1992. Carbon storage in upland forests of the Lake states. Soil Sci. Soc. Am. J. 56:935943.[Abstract/Free Full Text]
- Grossman, R.B., D.S. Harms, M.S. Kuzila, S.A. Glaum, S.L. Hartung, and J.R. Fortner. 1998. Organic carbon in deep alluvium in Southeast Nebraska and Northeast Kansas. p. 4556. In R. Lal et al. (ed.) Soil processes and the carbon cycle. Adv. Soil Sci. CRC Press, Boston, MA.
- Homann, P.S., P. Sollins, H.N. Chappell, and A.G. Stangenberger. 1995. Soil organic carbon in a mountainous, forested region: Relation to soil characteristics. Soil Sci. Soc. Am. J. 59:14681475.[Abstract/Free Full Text]
- Homann, P.S., P. Sollins, M. Fiorella, T. Thorson, and J.S. Kern. 1998. Regional soil organic carbon storage estimates for western Oregon by multiple approaches. Soil Sci. Soc. Am. J. 62:789796.[Abstract/Free Full Text]
- Hooker, T.D., and J.E. Compton. 2003. Forest ecosystem C and N accumulation during the first century after agricultural abandonment. Ecological Applications 13:299313.
- Huntington, T.G., D.F. Ryan, and S.P. Hamburg. 1988. Estimating soil nitrogen and carbon pools in a Northern hardwood forest ecosystem. Soil Sci. Soc. Am. J. 52:11621167.[Abstract/Free Full Text]
- Kabrick, J.M., M.K. Clayton, and K. McSweeney. 1997. Spatial patterns of carbon and texture on drumlins in Northeastern Wisconsin. Soil Sci. Soc. Am. J. 61:541548.[Abstract/Free Full Text]
- Kern, J.S. 1994. Spatial patterns of soil organic carbon in the contiguous United States. Soil Sci. Soc. Am. J. 58:439455.[Abstract/Free Full Text]
- Kern, J.S., D.P. Turner, and R.F. Dodson. 1998. Spatial patterns in soil organic carbon pool size in the Northwestern United States. p. 2944. In R. Lal et al. (ed.) Soil processes and the carbon cycle. Adv. Soil Sci. CRC Press, Boston, MA.
- Lal, R., J. Kimble, and R. Follett. 1998. Knowledge gaps and researchable priorities. p. 595604. In R. Lal et al. (ed.) Soil processes and the carbon cycle. Adv. Soil Sci. CRC Press, Boston, MA.
- Lathrop, R.G., Jr., J.D. Aber, and J.A. Bognar. 1995. Spatial variability of digital soil maps and its impact on regional ecosystem modeling. Ecol. Modell. 82:110.
- Lawson, D. 1995. Quaternary geology and deglacial history of south central Rhode Island. M.S. thesis. Dep of Geology, Univ. of Rhode Island, Kingston.
- Mausbach, M.J., B.R. Brasher, R.D. Yeck, and W.D. Nettleton. 1980. Variability of measured properties in morphologically matched pedons. Soil Sci. Soc. Am. J. 44:358363.[Abstract/Free Full Text]
- Nelson, D.W., and L.E. Sommers. 1996. Total carbon, organic carbon, and organic matter. p. 9611010. In D.L. Sparks et al. (ed.) Methods of soil analysis. Part 3. Book Series 5. SSSA, Madison, WI.
- Quinn, A.W. 1971. Bedrock geology of Rhode Island. USGS Bull. no. 1295.U.S. Geological Survey, Reston, VA.
- Rapalee, G., S.E. Trumbore, E.A. Davidson, J.W. Harden, and H. Veldhuis. 1998. Soil carbon stocks and their rates of accumulation and loss in a boreal forest landscape. Global Biogeochem. Cycles 12:687701.
- Rector, D.D. 1981. Soil survey state of Rhode Island. USDA-SCS, U.S. Gov. Print. Office, Washington DC.
- Schafer, J.P., and J.H. Hartshorn. 1965. The Quaternary of New England. p. 113125 In H.E. Wright, and D.G. Frey (ed.) The Quaternary of the United States. Princeton University Press, Princeton, NJ.
- Schafer, J.P. 1981. Glacial geology in southern Rhode Island. p. 229235 In J.C. Boothroyd and O.D. Hermes (ed.) Guidebook to geologic field studies in Rhode Island and adjacent areas. New England Intercollegiate Geologic Conference, University of Rhode Island, Kingston.
- Schimel, D.S., I.G. Enting, M. Heimann, T.M. Wigley, D. Raynaud, D. Alves, and U. Siegenthaler. 1994. CO2 and the carbon cycle. p. 3571. In J.T. Houghton et al. (ed.). Climate change. Cambridge Univ. Press, New York.
- SPSS, Inc.1999. SPSS base 10.0 applications guide. SPSS, Inc., Chicago, IL.
- Stone, E.L., W.G. Harris, R.B. Brown, and R.J. Kuehl. 1993. Carbon storage in Florida Spodosols. Soil Sci. Soc. Am. J. 57:179182.[Abstract/Free Full Text]
- Stolt, M.H. 1998. Aeolian soils in Rhode Island. Glacial geology in southern Rhode Island. p. 112. In J.C. Boothroyd and O.D. Hermes (ed.) Guidebook to geologic field studies in Rhode Island and adjacent areas. New England Intercollegiate Geologic Conference, University of Rhode Island, Kingston.
- Torn, M.S., S.E. Trumbore, O.A. Chadwick, P.M. Vitousek, and D.M. Hendricks. 1997. Mineral control of soil organic carbon storage and turnover. Nature 389:170173.
- United States Department of Energy. 1999. Carbon sequestration state of the science: Carbon sequestration research and development. Office of Science, Office of Fossil Fuel Energy, U.S. Department of Energy, National Technical Information Service, Springfield, VA.
- Vogt, K.A., D.J. Vogt, S. Brown, J.P. Tilley, R.L. Edmonds, W.L. Silver, and T.G. Siccama. 1995. Dynamics of forest floor and soil organic matter accumulation in boreal, temperate, and tropical forests. p. 159177. In A.F. Bouwman (ed.) Soils and the greenhouse effect. John Wiley & Sons, Chichester, NY.
- Wright, W.R., and E.H. Sautter. 1988. Soils of Rhode Island Landscapes. Univ. of Rhode Island Agric. Exp. Stn. Bull. 429, University of Rhode Island, Kingston.
- Young, F.J., J.M. Maatta, and R.D. Hammer. 1991. Confidence intervals for soil properties within map units. p. 213229. In M.J. Mausbach and L.P. Wilding (ed.) Spatial variabilities of soils and landforms. SSSA Spec. Pub. 28. SSSA, Madison, WI.
- Zdruli, P., H. Eswaran, and J. Kimble. 1995. Organic carbon content and rates of sequestration in soils of Albania. Soil Sci. Soc. Am. J. 59:16841687.[Abstract/Free Full Text]
This article has been cited by other articles:

|
 |

|
 |
 
S. A. Wills, C. L. Burras, and J. A. Sandor
Prediction of Soil Organic Carbon Content Using Field and Laboratory Measurements of Soil Color
Soil Sci. Soc. Am. J.,
March 12, 2007;
71(2):
380 - 388.
[Abstract]
[Full Text]
[PDF]
|
 |
|