|
|
||||||||
a Dep. of Natural Resource Ecology & Management, Iowa State Univ., Ames IA 50011
b USDA-ARS, National Soil Tilth Lab., Ames IA 50011
c Dep. of Agronomy, Iowa State Univ., Ames IA 50011
* Corresponding author (arussell{at}iastate.edu)
| ABSTRACT |
|---|
|
|
|---|
Abbreviations:
b, bulk density CC, continuous corn CCOA, corncornoatalfalfa CE, Carlo-Erba COAA, cornoatalfalfaalfalfa CS, cornsoybean MAP, mean annual precipitation PMC, potential mineralization of carbon POC, particulate organic carbon SIC, soil inorganic carbon SOC, soil organic carbon SOM, soil organic matter TN, total nitrogen WB, Walkley-Black
| INTRODUCTION |
|---|
|
|
|---|
Within the current conceptual framework, management practices that optimize cropping systems and N fertilization are believed to offer the greatest potential for increasing SOC storage in agricultural soils (Lal et al., 1999; Lal, 2002). For example, analysis of long-term experiments indicated that increasing crop rotation complexity increased SOC sequestration by 20 g C m2 yr1, on average (West and Post, 2002). In long-term experiments in Canada, SOC sequestration rates were 50 to 75 g C m2 yr1 in well-fertilized soils with optimum cropping systems (Dumanski et al., 1998). In contrast, long-term experiments in the northern Great Plains (ND) have shown that fertilizer N increased crop residue returns, but generally did not increase SOC sequestration (Halvorson et al., 2002). Timing and intensity of tillage also must be taken into account in the design of best management practices for maximizing SOC sequestration (Studdert and Echeverría, 2000).
Changes in SOC due to management practices are difficult to quantify because changes occur slowly, are relatively small compared with the total SOC pool size, and vary both spatially and temporally (Paustian et al., 1997; Russell et al., 2004). Soil C fractions, such as POC and potentially mineralizable carbon (PMC), are expected to be more responsive to management than total SOC (Powlson and Jenkinson, 1981; Cambardella and Elliott, 1992). Thus, these fractions may serve as harbingers of future changes in total SOC that are presently undetectable. Whereas PMC provides an index of labile SOC, POC is a measure of SOC that turns over on intermediate time scales (Cambardella and Elliott, 1992). Resistant SOC generally has longer turnover times, and thus provides information about the long-term potential for SOC sequestration (Paul et al., 2001). Because these soil C fractions are characterized by their turnover times, their pool sizes provide insight into the effects of management practices that could not be gleaned from studies of total SOC alone.
Model simulations indicate that reduced tillage and greater inputs of organic matter can return soils to precultivation levels of SOC (Paustian et al., 1998). Modeled potential C sequestration, however, may not be readily attainable (Ingram and Fernandes, 2001) due to various limiting factors (e.g., climate, time, and less-than-optimum management). Realistic estimates of SOC sequestration potential are best derived from careful analysis of long-term, replicated experiments in sites without confounding histories. Currently, there is a need for more data for some of the most productive and intensively managed soils of the world, Mollisols in the midwestern Corn Belt. Well-managed long-term experimental sites are rare, but thanks to the diligence of our predecessors we had access to two sites with experiments designed for testing long-term effects of N fertilization and cropping systems under conventional management. The objectives of this study were to quantify the long-term effects of N fertilization and cropping system on SOC stocks and fractions, to determine the rate of change in SOC during the last 12 yr, and to advance understanding of the mechanisms by which these management systems influence soil C sequestration.
| MATERIALS AND METHODS |
|---|
|
|
|---|
At both sites, all cropping systems were tile drained and rain fed. Primary tillage in the fall following corn and alfalfa was chisel plowing at Nashua and moldboard plowing at Kanawha. At both sites, the soils were always disked once in the spring, hence spring disking was the only tillage used for soybean residues. Nitrogen fertilizer treatments were applied in spring before disking, using granulated urea at rates of 0, 90, 180, and 270 kg N ha1. At Kanawha, N rates had been increased in 1971, and again in 1984 to keep pace with contemporary practices. The N rates were 0, 34, 68, and 136 kg N ha1 from 1954 until 1970. In 1971, the 34-kg N rate was changed to 202 kg N ha1. A final change of N rates was made in 1984, when the 68-kg rate was changed to 90 kg, the 136-kg rate to 180 kg, and the 202-kg rate was changed to 270 kg. Since 1984, the N rates evaluated have been 0, 90, 180, and 270 kg N ha1, the same as at Nashua. These N-rate treatments will be referred to as 0-N, 90-N, 180-N, and 270-N. An undisturbed prairie (Hayden Prairie), located 50 km north of Nashua, was also sampled to provide a reference point. The predominant soil series were similar to those at Nashua. The dominant plant species included Poa pratensis L., Kentucky bluegrass; Andropogon gerardii Vitman, big bluestem; Heliopsis helianthoides (L.) Sweet, smooth oxeye; Rosa arkansana Porter, prairie rose; Calamagrostis canadensis Michx. (Beauv.), Bluejoint reedgrass; and Thalictrum dasycarpum Fisch. & Avé-Lall, purple meadow rue (Christensen, 1996). The Kanawha site is 120 km west of Nashua. We had no reference for the Kanawha site, because native prairie on similar soils no longer exists near Kanawha. Mean annual precipitation (MAP) was 806 mm at Kanawha, ranging from 533 mm (1989) to 1221 mm (1993) during a 53-yr period (data from Clarion station) (19512003) (Iowa Environmental Mesonet, 2004). At Nashua, MAP was 847 mm, with a range of 421 mm (1989) to 1304 mm (1999) (Charles City station). Hayden Prairie MAP was 845 mm and ranged from 492 mm (1989) to 1161 mm (1999) (Osage station).
Sampling Protocol
For an intensive study at the Nashua site in 2001, soil samples (015 cm) were collected monthly (AprilNovember) from plots of the 0- and 180-N treatments of the CC, CS, and CCOA systems. For an extensive study in 2002, soil profile samples (05, 515, 1530, 3050, 5075, and 75100 cm) were collected from all N treatments and cropping systems at both sites during a single postharvest sampling (in October) before fall tillage, and also from the native prairie (in late August). Soil cores were taken using a 3.2-cm-diam. probe in 2001, a 4.1-cm-diam. probe in 2002 for the agricultural plots, and a 6.0-cm core for the prairie. Six soil cores were taken from each plot, at points randomly selected within a design stratified by position: in the cornrow, the midpoint between cornrows, and halfway between those two positions. The six cores, collected at Nashua and Kanawha from each plot and soil depth, were bulked into a composite sample. The six cores that were taken from the prairie were analyzed separately (because replicates were not applicable). To avoid edge effects, no samples were collected within 1 m of the borders of the contiguous plots. All samples were air-dried, roots were removed, rock masses and volumes were determined, and soil was passed through a 2-mm sieve. Subsamples were dried at 105°C to determine conversion factors to a 105°C dry-weight basis.
For the intensive study, we measured SOC, total nitrogen (TN), PMC, POC, and resistant C (non-acid-hydrolyzable), for samples collected from April to November for one depth increment (015 cm). For the extensive study, we measured bulk density (
b), SOC, TN, and soil inorganic carbon (SIC) for each depth increment, while PMC was measured only for the 5- to 15-cm depth.
Laboratory Evaluations
Total soil C and TN were determined by dry combustion using a Carlo-Erba NA1500 NSC elemental analyzer (Haake Buchler Instruments, Paterson, NJ) (CE method). Total SOC concentration was calculated as the difference between total C and SIC, which was measured by the modified pressure-calcimeter method (Sherrod et al., 2002). Storage of SOC was calculated as the product of
b, soil thickness, and SOC concentration. Bulk density was determined for each depth increment by the soil core method, using the samples collected in October 2002 (Blake and Hartge, 1986). Measurements were corrected to exclude rocks. Sampling of unequal soil masses across treatments and across time can introduce bias into calculations of SOC stocks (Ellert et al., 2001). To minimize this bias, we used mean
b across all plots within sites and depths in our calculations of SOC stocks, and sampled extensively to characterize each site when
b would be least variable among treatments (postharvest, but before tillage). Plots at both sites were uniformly level, so differences in losses due to erosion were assumed to be equal. In the agricultural plots, the layer most sensitive to treatment effects, 0 to 15 cm, was situated well within the homogenized plow layer, such that samples at fixed depths within this layer were comparable among treatments and across time. Low
b in the prairie resulted in a soil sample mass of only 6.57 kg, compared with an average of 10.12 kg per fixed volume (0- to 100-cm depth) in the Nashua agricultural plots. To compare soils on an equivalent-mass basis with the prairie soils, SOC stocks in the agricultural soils were summed to a depth of 67.7 cm.
Historical SOC data were not available from the beginning of either experiment; however, published data from a 1990 soil sampling for selected plots were available (Robinson, 1993). To assess changes in SOC across time, we compared our SOC concentrations for 0 to 15 cm measured in 2002 with concentrations reported by Robinson (1993) for the same depth and treatments in 1990: CC (0-N and 180-N), CS (180-N), CCOA (180-N), and COAA (180-N, Kanawha only). At Nashua, measurements in 1990 and 2002 correspond to Years 11 and 23 of the experiment. At Kanawha, the 0-N treatment had been in effect since 1954, so these measurements correspond to Years 36 (by 1990) and 48 (by 2002) without N fertilization. In the 1990 study, SOC was determined by the Walkley-Black (WB) method (Robinson et al., 1996), whereas we used dry combustion, the CE method. Thus, to make comparisons between sample dates, it was necessary to determine site-specific conversion factors between the two methods. We analyzed SOC by both WB and CE in two cropping systems (CS and CCOA), two N levels (0- and 180-N), and two depths (05 and 515 cm) for all blocks at Nashua (24 samples total) and Kanawha (16 samples). We assumed a constant factor of 1.724 for conversion of soil organic matter (SOM) determined by WB to SOC (Nelson and Sommers, 1996). The relationships between SOC (g kg1) estimated by the two methods were: Nashua, CE = 1.05 x WB 6.19 (R2 = 0.989); Kanawha, CE = 0.892 x WB 5.56 (R2 = 0.984). Historical data for WB conversion factors were not available, so the untested assumption was that there was no change.
Potentially mineralizable C was assayed by measures of CO2C released during 28-d laboratory incubations at 23°C (Paul et al., 2001, Russell et al., 2004). For the 2001 samples, field-moist soil was sieved (4 mm) and kept refrigerated until incubations began, within 48 h of collection. For each sample time and plot, 20 g of soil were incubated in vials placed in a one-pint (473 mL) Ball jar. Approximately 1.5 mL of distilled water was added to the bottom of each Ball jar to keep the atmosphere hydrated. Duplicates were run on all samples. Soil moisture was low at all sampling times from May to November, so distilled water was added to adjust samples to 50% water-filled pore space before incubations in those months. Subsamples were dried at 105°C to determine moist-to-dry-weight conversion factors. Rate of CO2C released was measured periodically (before CO2 concentrations reached 4%) by flushing the jar headspace through an infrared gas analyzer for 2001 samples (IRGA, Rosemount Analytical, Inc., Orrville, OH) (Paul et al., 2001). In 2002, the same principles but slightly different methods were applied to accommodate the extensiveness of the field sampling and to conform to protocols of a multisite project. The depths sampled were 0 to 15 cm in 2001 and 5 to 15 cm in 2002. Because these soils were mold-board and chisel plowed, the top 15 cm was well homogenized; preliminary SOC data indicated that the two sampled layers did not differ significantly. In the 2002 sampling, air-dried 2-mm-sieved soils were used in the assays of 5-g samples with three laboratory replicates using an IRGA (Model 3300, Automated Custom Systems, Orange, CA). Soil samples were hydrated to 50% water-filled pore space before incubation. The PMC was high in the prairie soils, so to avoid surpassing the measurement capacity of the IRGA, sample masses were reduced to 3 g, and six replicates were analyzed. Incubation-tube atmosphere was kept hydrated by piping the air supply through a water bath.
Particulate organic C was quantified according the method of Cambardella and Elliott (1992). A single 10-g sample per plot was first dispersed by shaking 18 h on a reciprocal shaker in 30 mL of 5 g L1 sodium metaphosphate. The sample was rinsed thoroughly with distilled water and passed through a 53-µm sieve. The (silt + clay) material that passed through the sieves was captured in a receiving pan, dried at 70°C, ground with mortar and pestle, and analyzed for SOC and TN. The amount retained on the sieve (POC plus sand) was dried at 50°C and weighed. The POC was calculated as the difference between SOC and silt + clay SOC, corrected for the amount recovered. Particulate OC was analyzed in the samples taken from April to July.
The passive or resistant soil C fraction was determined by refluxing 1 g soil in 5 mL 6 M HCl at 115°C for 18 h in a temperature-controlled digestion block (Paul et al., 2001). Refluxed samples were washed in 1000 mL deionized water, dried at 55°C, and ground in a mortar and pestle for analyses of C and N. Values were corrected for non-C content lost during hydrolysis. Although material in this fraction is relatively old, it may also contain newer SOC from modern plant lignin that is not hydrolyzed (Paul et al., 2001).
Statistical Analyses and Calculations
Analyses of variance for testing cropping system and N fertilization effects on measurements were conducted for a split-plot randomized block design using PROC GLM of SAS (SAS Institute, 1990) with all effects treated as fixed (Littell et al., 1991). Treatment means were compared using Tukey's studentized (HSD) multiple-comparison test. Planned comparisons within a cropping system were unfertilized (0-N) vs. fertilized treatments (90-N, 180-N, and 270-N). At Kanawha, planned comparisons within an N treatment consisted of contrasts between cropping systems with and without alfalfa in the rotation (CC and CS vs. CCOA and COAA). For the response variables measured at monthly intervals in 2001 (POC, PMC, and resistant C), ANOVA was first conducted using a repeated measures design to test for differences among the monthly sampling times (Littell et al., 1991). Relationships between PMC and SOC were assessed using correlation analysis (SAS Institute, 1990).
Gross and net C sequestration were calculated as described by Schlesinger (2000), using a factor of 1.4 to represent the full accounting of emissions of CO2 associated with all aspects involved in N fertilizer production, transportation, and application (Cole et al., 1993; Izaurralde et al., 1998). Nitrogen fertilization effects on SOC were calculated as the difference between SOC in fertilized and unfertilized plots sampled in 2002 because no data were available from the beginning of the experiment.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
b (Fig. 1). In the 0- to 15-cm depth interval, the prairie's SOC (66.1 ± 1.2 g kg1) was 2.8 times higher than the highest SOC for the cropping systems at Nashua, CCOA, whether fertilized or not (22.7 ± 0.3 g C kg1 in 0-N and 23.6 ± 0.7 g C kg1 in 270-N). Throughout the entire 100-cm profile,
b was substantially lower in the prairie, and C/N tended to be higher.
|
|
High variability in SOC, relative to differences among treatments, contributed to the lack of statistical significance. At Kanawha, other factors contributed to the lack of significant effect of N fertilization: only two replicates for each treatment, high variability between the two blocks, and a high buffering capacity of the soil due to the high soil clay content. Omay et al. (1997) found no significant increase in SOC with N fertilization at two long-term sites in Kansas. Neither N fertilization nor cropping system influenced concentrations of SIC at Kanawha (Fig. 1), the only site in which soils contained carbonates.
Effects of Cropping System on Soil Organic Carbon Stocks
At both sites, cropping systems influenced SOC stocks, but the effects depended on N fertilization and depth (Table 1). At Nashua, in the 0- to 15-cm depth, mean SOC stocks (across all N treatments) increased significantly (P = 0.032) for cropping systems in the following order: CS << CC < CCOA. Soil OC stocks were significantly lower in CS than in CCOA at the 0- and 270-N level (Table 1). At Kanawha, the effect of cropping system was not significant (P = 0.21 for 015 cm). In pooled comparisons, however, cropping systems that contained alfalfa in the rotation (CCOA and COAA) had significantly higher SOC stocks than systems without alfalfa (CC and CS) at all N levels except 180-N (Table 1). Omay et al. (1997) had similar results in that SOC was lower under CS compared with CC at two long-term sites in Kansas, on irrigated soils with lower SOC than at our sites. In the 15- to 100-cm depth, across all N treatments, the effect of cropping system was not significant at either site (P = 0.222 and 0.181, Nashua and Kanawha, respectively). Within the 270-N treatment, however, systems that contained alfalfa had significantly greater SOC stocks than did CS (both sites) and CC (Kanawha only) (Table 1).
The fact that tillage, as well as cropping system, influences SOC stocks is well documented (Lal et al., 1997; Robertson et al., 2000; Hao et al., 2002). In our study, because crop effects are confounded with tillage-frequency effects, we cannot discern whether the increase in SOC stocks in systems containing oat and alfalfa was caused by traits of the crops in the sequence, or the absence of tillage during 2 to 3 yr of the 4-yr sequence.
Change in Soil Organic Carbon and Net Carbon Sequestration
During the 12 yr between the 1990 and 2002 samplings at Nashua, the rates of SOC loss (015 cm) in the CS (180-N) and CC (0-N), and gains in CC and CCOA (180-N) did not differ significantly from zero (Table 2). During the same 12 yr at Kanawha, SOC increased for all systems, but this increase was significant only in the fertilized CC and COAA cropping systems. The two sites had similar trends among cropping systems in that the fertilized CS system had the lowest increase in SOC at Kanawha and the greatest decline at Nashua, even in comparison with unfertilized CC. Similarly, Studdert and Echeverría (2000) found that SOC loss increased with soybean in the rotation.
|
15% of CO2 emissions from the state's fossil fuel emissions in 1997 (USEPA, 2001). These researchers estimated that Iowa's agricultural soils could sequester two to three times more C if the 50% of cropland now managed under conventional, intensive tillage were converted to no-till practices. They also predicted soil C gains of 0.25 Mg ha1 yr1 by changing from intensive- to moderate-till practices. The trends in our Nashua data were consistent with this latter prediction, although ours was not a straightforward test regarding effects of lower tillage frequency because residue inputs also differed (Table 2). Our analysis of changes in soil C stocks was limited to the last 12 yr of the experiments and the nine treatments for which we have historical data. To compare net SOC sequestration (015 cm) among all fertilized treatments during the entire course of the experiment, we used the unfertilized treatments within each site and cropping system as the baseline for comparison of SOC stocks in the fertilized treatments. In general, with such small differences in soil C storage between fertilized and unfertilized treatments, the calculated emissions of CO2 from fertilizer production nearly equaled or exceeded the increase in soil C accumulated in 17 out of 21 of our cropping systems under the three fertilization levels (Table 3). By Schlesinger's (2000) method of calculation, as much as 793% of the increase in soil C stored was released as CO2 during the fertilizer production. The analysis was not meaningful in six of the systems due to a negative difference in SOC stocks (i.e., SOC stocks were higher under the unfertilized treatment). The only system that showed a net gain in SOC at Nashua was the CS under 90-N treatment. At Kanawha, the CC (90-N and 180-N) and COAA (all three N levels) had net gains in SOC. With N fertilizer added only once every 4 yr, the COAA system probably achieved a favorable balance because the CO2 emissions associated with fertilizer production and use were low due to lower N inputs, and did not offset the increase in SOC from the cropping system. Our analyses by Schlesinger's (2000) method differed from recent studies of C emissions in agriculture in that we accounted for CO2 emissions associated only with N fertilizer, and our conversion factor (1.4) was higher than in other studies (West and Marland, 2002). Nevertheless, our estimates of net C flux were only 2% higher than the average reported for corn under conventional tillage (West and Marland 2002), and our estimates of CO2 emissions fall within the range expected due to variability among sites (Marland et al., 2003).
|
|
Soil Organic Carbon Mineralization
Soil C/N ratios provide insight into SOC quality and relative decomposability. Soil C/N ratios were higher under native prairie than the agricultural systems, and remained relatively constant with depth, whereas C/N ratios decreased substantially with depth in agricultural systems (Fig. 1). Among the agricultural soils at Nashua, C/N ratios were higher under CC than CS and CCOA (Table 5). Variability in soil C/N between blocks was high at Kanawha, such that differences were not significant. The high C/N ratios in the prairie soils and the soils under CC at Nashua suggest that the SOC was relatively less decomposed in these systems relative to the CS and CCOA systems. However, C/N ratios tell only part of the story.
|
|
Our data are consistent with the concept that SOC sequestration is influenced by both inputs of organic C and decomposition losses (Tables 4, 6). At Nashua, SOC stocks were lowest in CS (86 Mg C ha1, mean across all N treatments, 0100 cm), the cropping system with the lowest residue inputs (estimated from long-term yield data reported by Mallarino and Pecinovsky, 1999) and relatively high PMC in spring. Soil C stocks were intermediate (102 Mg C ha1) in CC, with intermediate yields and low PMC in spring. Stocks were highest (108 Mg C ha1) in CCOA, the system with the highest yields and relatively high PMC. Trends were similar at Kanawha, with SOC stocks high in CCOA and COAA (179 and 167 Mg C ha1) and low in CS and CC (138 and 133 Mg C ha1). Although yields were generally higher with N fertilization, higher PMC in spring may have offset SOC gains, as evidenced by the general lack of significant increase in SOC stocks with N fertilization.
Soil Organic Carbon Fractions
Both POC and resistant SOC increased in the order CS < CC < CCOA (Table 7). Particulate OC was significantly lower under CS than CC and CCOA (180-N only), and resistant C was significantly higher in CCOA than CS (0-N only). Neither POC nor resistant C was influenced by N fertilization. Liebig et al. (2002) studied several soil properties in a long-term cropping system in Nebraska and observed that N fertilization (180 kg ha1 N) increased POC in only one of four cropping systems examined, the CC system. In our study, resistant C averaged 67% of SOC in the cropped systems, and 91% in the prairie. Particulate OC and resistant SOC concentrations in the prairie were 2.6 and 3.9 times higher, respectively, than the highest values measured for the cropping systems, the unfertilized plots of the CCOA system (Table 7).
|
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| NOTES |
|---|
|
|
|---|
Received for publication July 1, 2004.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
J. L. Pikul Jr., J. M. F. Johnson, T. E. Schumacher, M. Vigil, and W. E. Riedell Change in Surface Soil Carbon under Rotated Corn in Eastern South Dakota Soil Sci. Soc. Am. J., October 30, 2008; 72(6): 1738 - 1744. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. J. Purakayastha, D. R. Huggins, and J. L. Smith Carbon Sequestration in Native Prairie, Perennial Grass, No-Till, and Cultivated Palouse Silt Loam Soil Sci. Soc. Am. J., February 15, 2008; 72(2): 534 - 540. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A. Khan, R. L. Mulvaney, T. R. Ellsworth, and C. W. Boast The Myth of Nitrogen Fertilization for Soil Carbon Sequestration J. Environ. Qual., October 24, 2007; 36(6): 1821 - 1832. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. F. Grant, T. J. Arkebauer, A. Dobermann, K. G. Hubbard, T. T. Schimelfenig, A. E. Suyker, S. B. Verma, and D. T. Walters Net Biome Productivity of Irrigated and Rainfed Maize Soybean Rotations: Modeling vs. Measurements Agron. J., October 15, 2007; 99(6): 1404 - 1423. [Abstract] [Full Text] [PDF] |
||||
![]() |
G. E. Varvel Soil Organic Carbon Changes in Diversified Rotations of the Western Corn Belt Soil Sci. Soc. Am. J., February 2, 2006; 70(2): 426 - 433. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. E. Russell, D. A. Laird, and A. P. Mallarino Nitrogen Fertilization and Cropping System Impacts on Soil Quality in Midwestern Mollisols Soil Sci. Soc. Am. J., January 6, 2006; 70(1): 249 - 255. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||