# Tillage Effects on Spatiotemporal Variability of Particulate Organic Matter

Authors: Juhwan Lee; Emilio A. Laca; Chris van Kessel; Dennis E. Rolston; Jan W. Hopmans; Johan Six
Journal: Applied and Environmental Soil Science
 (2009)

---

## Abstract

This study was performed to evaluate effects of no-till (NT) and standard tillage (ST) on POM in two 15-ha neighboring fields from 2003 to 2004. We also evaluated the effects of minimum tillage (MT) on POM after both NT and ST fields were converted to MT in the summer of 2005. We quantified C and N stocks of three size fractions (53–250, 250–1000, and 1000–2000μm) of POM (0–0.15 m depth). The POM-C 53–250 μm and 250–1000 μm fractions decreased by 25% and 36% after six months under ST, whereas relatively little change occurred under NT, suggesting significant tillage effects over the period 2003-2004. Only small changes in POM content then occurred under MT on both fields. Changes in POM-N were similar to POM-C changes upon tillage conversions. This suggests that reduced tillage did not lead to soil C increase compared to ST but may help maintain the level of soil C for a typical California farming system. Short-term, field level variability of POM was primarily affected by tillage and was further influenced by clay content, bulk density, and scale of observation.


---

## Body

## 1. Introduction

There is growing evidence that the amount of soil organic C and N storage tends to increase following the conversion from ST to reduced tillage, such as NT and MT [1, 2]. Therefore, agricultural soils under reduced tillage practices may have the potential to mitigate greenhouse gas (GHG) emissions [3, 4]. However, the effect of reduced tillage on soil C and N changes is still not conclusive, as the relationship between soil organic matter (SOM) and its controlling processes is highly variable across all spatial and temporal scales [4–6]. This is likely because a number of biotic and abiotic variables controlling the loss and sequestration of SOM frequently interact over space and time [7].Many field studies have compared tillage systems using controlled experiments where treatments are randomly replicated at relatively small scales [8]. Patterns and related processes analyzed at these small-scale studies, however, are often difficult to be validated at larger scales due mostly to scale dependency [9]. Recently, on the other hand, the scientific community is recognizing the importance of conducting research at scales of agricultural management [10]. Unfortunately, it is often found to be physically impossible or economically undesirable to replicate large-scale field trials [11, 12], especially in farmer’s fields. Therefore, large-scale variation itself is being embraced because it can be used to test and refine interactions among controlling factors (e.g., soil properties, microclimate, etc.) and agricultural management practices [7, 13]. Understanding these interactions at a larger scale is of primary importance for a comprehensive evaluation of potential benefits and drawbacks of reduced tillage.For unreplicated field studies, one needs to explore alternative statistical methods to compare large-scale systems that are characterized by high variation. A Before-After Control-Impact (BACI) approach can be used to detect treatment-induced changes in field-scale processes [14, 15]. To account for problems of spatial and temporal variation, the field should be paired with a similar field not exposed to the treatment of interest. It is not reasonable to expect that the two paired fields will have a similar degree of changes in processes. Thus, results from BACI are sometimes ambiguous because more than one interpretation is possible on paired observations [16]. However, if obvious differences in spatial and temporal patterns develop after dividing a field into a continuing treatment versus a contrasting treatment then one can reasonably postulate that they are caused by the treatments. More detailed discussions of BACI are presented elsewhere [14, 16–19].At the field scale, it is difficult to detect short-term changes in SOM in response to management changes because of large background levels and high spatial variability of SOM [20, 21]. Studies suggest that measures of POM can provide an indication of longer-term changes of SOM due to tillage effects [22, 23]. The POM fraction can contribute up to 40% of the total SOM in near-surface horizons, primarily derived from root materials in various stages of decomposition [22, 24]. Cambardella and Elliott [22] showed that POM has a faster C turnover rate than whole soil organic C or mineral-associated C. Within the POM fraction, the coarse (250–2000 μm) POM fraction is the most sensitive to changes in residue input and decomposition rate by tillage [1]. The fine (53–250 μm) POM is, on the other hand, relatively less sensitive to changes in management and is a more stabilized POM fraction [23].The effects of reduced tillage on SOM and POM levels have been intensively investigated in the Great Plains and Corn Belt Regions [3, 25, 26]. In contrast, few such studies have been conducted in California [27]. Differences in the magnitude of tillage-induced changes in SOM between sites and regions are still poorly understood for Californian agroecosystems. In particular, no field-scale studies evaluating tillage effects on SOM changes have been performed in California. Therefore, we hypothesized that in a typical California farming system, the change in the rate of soil C cycling, induced by reduced tillage, leads to a predictable increase in soil C sequestration across the field. The objectives of this study were (i) to determine horizontal and temporal changes in POM contents under ST, NT, and MT practices in a Mediterranean agroecosystem and (ii) to determine effects of selected soil physical properties on POM contents across fields under different tillage systems.


## 2. Materials and Methods

### 2.1. Site Description and Sampling Scheme

In 2003, an experimental site was established on a 30-ha agricultural field (38°36’ N, 121°50’ W) located in the Sacramento Valley of California. The site has been land planed for irrigation with a slope of about 1%. Irrigation was primarily by furrow irrigation, but sprinkler irrigation was used for seed germination if necessary (Figure 1). The mean annual temperature and precipitation were 16.1°C and 564 mm, respectively. Two major and one minor soil series occur on the site: Myers clay (fine, montmorillonitic, thermic Entic Chromoxerert), Hillgate loam (fine, montmorillonitic, thermic Typic Palexeralf), and San Ysidro loam (fine, montmorillonitic, thermic Typic Palexeralf) [28].Precipitation, irrigation, and daily mean air temperature at the field site, September 2003–December 2005.
(a)(b)The site was managed under NT management in 2001 and 2002. Following sudangrass (Sorghum vulgare) in 2001 and maize (Zea  mays) in 2002, the field was seeded to winter wheat (Triticum  aestivum). After the harvest of the winter wheat in June 2003, the site was split into two fields in October 2003, with the north half of the site converted to ST management practices and the south half remaining under NT. The ST treatments consisted of one pass each of deep ripping to 45 cm, stubble disking, disking to 15 cm, grading, and listing the seed beds (in 0.75 m rows). The terminology used to describe these tillage systems in California is that presented by Mitchell et al. [29]. Both fields remained fallow until maize was planted on April 12 and 13, 2004. At planting, urea-ammonium nitrate (32–0–0 NPK) was band applied (0.10 m depth) at a rate of 55 kg N ha-1. An additional 21 kg N ha−1 was broadcast applied as 8–24–6 NPK. In addition, 168 kg N ha-1 of N-fertilizer was side dressed at 15-cm depth on May 24 and 25, 2004. Goats grazed the field of wheat and maize residue for 1–4 weeks after harvest and before tillage in 2003 and after harvest in 2004. The whole field was fenced off and divided into four identical areas. Goats spent a similar amount of time in each area. Grazing resulted in a decrease of 53±2% (mean ± standard error) and 63±3% in the amount of aboveground residue C for wheat (2003) and maize (2004), respectively, (data not shown). The difference between the fields was on average 4% in 2003 and 1% in 2004. Manure returns of C were almost negligible (less than 3% of total C input from the crop). As goats were not observed to range widely during short time intervals, we assumed that grazing effects on C input levels were similar for the two fields, and homogeneous within each field.In May 2005, maize stubble was chopped in both fields, with three bed disc (to 15 cm) passes on the north side of the field and two passes on the south. Both fields also had one mulcher pass and pre-emergent herbicide incorporated. Male sunflowers (Helianthus annuus.) were planted on May 16, 2005, and females were planted on May 23. Fertilizer (UAN-32) was applied at a rate of 90 kg N ha-1 in a side-dress application on June 17-18.


### 2.2. Soil and Plant Sampling

Soil and plant samples were taken from the fields at 72 global position system referenced locations on 64 m regular grids (Figure2). In August 2003, the 1 m soil profile (0.05 m in diameter) was sampled at each sampling location two months prior to tillage and split into 0–0.15, 0.15–0.3, 0.3–0.5, 0.5–0.75, and 0.75–1 m depth increments to analyze spatial variations in soil C and N and their controlling factors before the imposition of the ST treatments. Soil cores were taken to a depth of 0–0.15 m at the same locations in April and October 2004 (6 and 12 months under ST versus NT), and in June 2005 (one month after the conversion of both ST and NT to MT). Wheat shoot samples were collected in the fields at the same 72 locations in June 2003 at harvest to estimate residue C input in soil (see below).Figure 2
Layout of sampling points. The sampling points on the 64 m grid are indicated with solid circles for the north and open circles for the south field.


### 2.3. Soil and Plant Analyses and Calculations

Field-moist soil samples were sieved to pass<8 mm, and subsamples were dried at 105°C overnight to determine gravimetric water content and bulk density (BD). The remaining sample was air-dried, mechanically crushed, and sieved to pass through a 2 mm mesh. Plant residues (>2 mm) were removed from the samples, and rock fragments were separated and weighed for BD and water content corrections. Soil particle size distribution was determined by the laser diffraction method using a Beckman-Coulter LS-230 with a 750 nm laser beam [30].Soils were analyzed for C and N stocks of total SOM and POM (i.e., 53–250, 250–1000, and 1000–2000μm size fractions). To separate POM and sand fractions, a 2-mm-sieved, air-dried subsample (30 g) was dispersed by shaking in 0.5% sodium hexametaphosphate solution (100 mL) for 18 hours [22]. The dispersed sample was passed through a series of three sieves (1000, 250, and 53 μm), and all fractions were thoroughly rinsed and transferred to preweighed pans. All POM fractions were oven-dried overnight at 65°C, weighed, ground and stored at room temperature.Plant samples were oven-dried at60°C for 48 hours, separated into residue and grain fractions, and then ground (<0.5 mm). The ground plant samples were further ground using a ball grinder for chemical analysis. Organic C and N content, and δ13C in the bulk soil, POM fractions, and plant samples were measured with a PDZ Europa ANCA CN analyzer connected to a PDZ Europa 20-20 Isotope Ratio Mass Spectrometer (SerCon, Cheshire, UK).The stocks of soil organic C and N were expressed on an equivalent soil mass basis [31], using the following equations:
(1)Msoil=BD·depthMelement=conc·Msoil·1000-conc·(Msoil-Msoilequiv)·1000,
where Melement is the equivalent C or N mass (g m-2), conc is concentration of organic C or N (g kg-1 soil), BD is bulk density (Mg m-3), depth is the depth of horizon (m), the factor 1000 is a factor for unit conversions, Msoil is soil mass to the depth of horizon (Mg m-2), and Msoilequiv is equivalent soil mass (Mg m-2), which is the mass of the least dense sample among the soil samples at each sampling time.To estimate the amount of new C input derived from fresh wheat residue into POM-C during the first eight months of the study (from August 2003 to April 2004), we used theδ13C values for POM fractions and wheat residues. A simple mixing model [32] was used to calculate soil POM-Cnew:
(2)P=(δt1-δt0δw-δt0)soilPOM-Cnew=eq·conct1×P,
where P is the proportion of C derived from wheat residues in the POM fractions between time 0 and 1, δt0 and δt1 are the δ13C values of the POM fraction at time 0 (August 2003) and 1 (April 2004), δw is the average δ13C value of wheat residue (-27.7%), and eq·conct1 is the concentration of POM-C (g C m-2) calculated on an equivalent soil mass basis for time 1.


### 2.4. Statistical Analysis

Variables were tested for non-normality using the D’Agostino-Pearson omnibus test in SAS [33, 34]. Outliers with an absolute value of more than 3 standard deviations from the mean were excluded from the descriptive statistics and geostatistical and multivariate regression analyses [35]. On average, there was less than one outlier for each variable for each sampling time.By accounting for spatial correlation, image maps of POM fractions were produced at each sampling time in Surfer 8 (Golden Software, Golden, CO). Briefly, simple linear regression was first performed on the 64-m regular grid with x and y coordinates as independent variables to identify if there were significant (P<.05) spatial trends in POM fractions. When a significant trend with direction existed, trend values were calculated at all sampling points using the trend model and then subtracted from the measured values. Semivariogram models were defined to determine the spatial structure of the residuals from the trend, and the residuals were then block-kriged and added back to the trends using their respective semivariogram models in GS+ (Gamma Design Software, St. Planinwell, MI). The best model was selected based on visual fit, reduced sums of squares, and r2 of the regression and through cross validation [36]. When there was no significant trend, semivariograms and estimation were made using the measured POM data. Due to a relatively small number (N=36) of samples in each field, anisotropic semivariograms were not considered for each of the data sets [36]. For data sets showing a very weak spatial correlation, point kriging with a linear model was used for interpolation.Repeated measurements over time were used to compare tillage-induced changes in POM-C and -N by the interaction offield (i.e., between the north and south fields) and sampling time (i.e., before versus after potential differences in tillage effects) using a mixed model in SAS [15, 34]. All the POM-C and-N fractions were log-transformed for BACI analysis. Significant main effects of the factors “field” and “sampling time” may indicate the presence of inherent differences and temporal fluctuations across and within the fields that persisted throughout the study. We assumed that these effects of “field” and “sampling time” were not directly controlled by the treatment. Due to difficulty in establishing the control field, we excluded the period of both fields being under MT for the BACI analyses. All BACI model residuals by field or sampling time were checked for independence by visually inspecting histograms and empirical semivariograms, and correlation analysis. A P-value of .02 was used to avoid a type I error.Standardized partial regression coefficients were used to assess the impact of soil texture, BD, and water content on POM. Regression coefficients were obtained by developing a mixed model for each response variable (after log transformation to meet normality) in R using the nlme package [37]. A final model was obtained by backward elimination, starting with a full model that included a full factorial of field, time, sand, clay and BD, minus all interactions of order equal to or greater than four. Location was the random effect. All effects and interactions with a P-value greater than .05 were removed, except when there was a significant interaction involving the effect. Random effects and residuals were checked for normality by visual inspection of quantile plots. Homogeneity of variance was checked by plotting residuals against predicted values. Spatial independence was checked by calculating the empirical semivariograms for residuals and random effects.Due to the unreplicated nature of the tillage treatments, the inference about the effect of ST versus NT on POM was limited to the two paired fields used in this study. When necessary, descriptive statistics were used to compare results between the fields.


## 2.1. Site Description and Sampling Scheme

In 2003, an experimental site was established on a 30-ha agricultural field (38°36’ N, 121°50’ W) located in the Sacramento Valley of California. The site has been land planed for irrigation with a slope of about 1%. Irrigation was primarily by furrow irrigation, but sprinkler irrigation was used for seed germination if necessary (Figure 1). The mean annual temperature and precipitation were 16.1°C and 564 mm, respectively. Two major and one minor soil series occur on the site: Myers clay (fine, montmorillonitic, thermic Entic Chromoxerert), Hillgate loam (fine, montmorillonitic, thermic Typic Palexeralf), and San Ysidro loam (fine, montmorillonitic, thermic Typic Palexeralf) [28].Precipitation, irrigation, and daily mean air temperature at the field site, September 2003–December 2005.
(a)(b)The site was managed under NT management in 2001 and 2002. Following sudangrass (Sorghum vulgare) in 2001 and maize (Zea  mays) in 2002, the field was seeded to winter wheat (Triticum  aestivum). After the harvest of the winter wheat in June 2003, the site was split into two fields in October 2003, with the north half of the site converted to ST management practices and the south half remaining under NT. The ST treatments consisted of one pass each of deep ripping to 45 cm, stubble disking, disking to 15 cm, grading, and listing the seed beds (in 0.75 m rows). The terminology used to describe these tillage systems in California is that presented by Mitchell et al. [29]. Both fields remained fallow until maize was planted on April 12 and 13, 2004. At planting, urea-ammonium nitrate (32–0–0 NPK) was band applied (0.10 m depth) at a rate of 55 kg N ha-1. An additional 21 kg N ha−1 was broadcast applied as 8–24–6 NPK. In addition, 168 kg N ha-1 of N-fertilizer was side dressed at 15-cm depth on May 24 and 25, 2004. Goats grazed the field of wheat and maize residue for 1–4 weeks after harvest and before tillage in 2003 and after harvest in 2004. The whole field was fenced off and divided into four identical areas. Goats spent a similar amount of time in each area. Grazing resulted in a decrease of 53±2% (mean ± standard error) and 63±3% in the amount of aboveground residue C for wheat (2003) and maize (2004), respectively, (data not shown). The difference between the fields was on average 4% in 2003 and 1% in 2004. Manure returns of C were almost negligible (less than 3% of total C input from the crop). As goats were not observed to range widely during short time intervals, we assumed that grazing effects on C input levels were similar for the two fields, and homogeneous within each field.In May 2005, maize stubble was chopped in both fields, with three bed disc (to 15 cm) passes on the north side of the field and two passes on the south. Both fields also had one mulcher pass and pre-emergent herbicide incorporated. Male sunflowers (Helianthus annuus.) were planted on May 16, 2005, and females were planted on May 23. Fertilizer (UAN-32) was applied at a rate of 90 kg N ha-1 in a side-dress application on June 17-18.


## 2.2. Soil and Plant Sampling

Soil and plant samples were taken from the fields at 72 global position system referenced locations on 64 m regular grids (Figure2). In August 2003, the 1 m soil profile (0.05 m in diameter) was sampled at each sampling location two months prior to tillage and split into 0–0.15, 0.15–0.3, 0.3–0.5, 0.5–0.75, and 0.75–1 m depth increments to analyze spatial variations in soil C and N and their controlling factors before the imposition of the ST treatments. Soil cores were taken to a depth of 0–0.15 m at the same locations in April and October 2004 (6 and 12 months under ST versus NT), and in June 2005 (one month after the conversion of both ST and NT to MT). Wheat shoot samples were collected in the fields at the same 72 locations in June 2003 at harvest to estimate residue C input in soil (see below).Figure 2
Layout of sampling points. The sampling points on the 64 m grid are indicated with solid circles for the north and open circles for the south field.


## 2.3. Soil and Plant Analyses and Calculations

Field-moist soil samples were sieved to pass<8 mm, and subsamples were dried at 105°C overnight to determine gravimetric water content and bulk density (BD). The remaining sample was air-dried, mechanically crushed, and sieved to pass through a 2 mm mesh. Plant residues (>2 mm) were removed from the samples, and rock fragments were separated and weighed for BD and water content corrections. Soil particle size distribution was determined by the laser diffraction method using a Beckman-Coulter LS-230 with a 750 nm laser beam [30].Soils were analyzed for C and N stocks of total SOM and POM (i.e., 53–250, 250–1000, and 1000–2000μm size fractions). To separate POM and sand fractions, a 2-mm-sieved, air-dried subsample (30 g) was dispersed by shaking in 0.5% sodium hexametaphosphate solution (100 mL) for 18 hours [22]. The dispersed sample was passed through a series of three sieves (1000, 250, and 53 μm), and all fractions were thoroughly rinsed and transferred to preweighed pans. All POM fractions were oven-dried overnight at 65°C, weighed, ground and stored at room temperature.Plant samples were oven-dried at60°C for 48 hours, separated into residue and grain fractions, and then ground (<0.5 mm). The ground plant samples were further ground using a ball grinder for chemical analysis. Organic C and N content, and δ13C in the bulk soil, POM fractions, and plant samples were measured with a PDZ Europa ANCA CN analyzer connected to a PDZ Europa 20-20 Isotope Ratio Mass Spectrometer (SerCon, Cheshire, UK).The stocks of soil organic C and N were expressed on an equivalent soil mass basis [31], using the following equations:
(1)Msoil=BD·depthMelement=conc·Msoil·1000-conc·(Msoil-Msoilequiv)·1000,
where Melement is the equivalent C or N mass (g m-2), conc is concentration of organic C or N (g kg-1 soil), BD is bulk density (Mg m-3), depth is the depth of horizon (m), the factor 1000 is a factor for unit conversions, Msoil is soil mass to the depth of horizon (Mg m-2), and Msoilequiv is equivalent soil mass (Mg m-2), which is the mass of the least dense sample among the soil samples at each sampling time.To estimate the amount of new C input derived from fresh wheat residue into POM-C during the first eight months of the study (from August 2003 to April 2004), we used theδ13C values for POM fractions and wheat residues. A simple mixing model [32] was used to calculate soil POM-Cnew:
(2)P=(δt1-δt0δw-δt0)soilPOM-Cnew=eq·conct1×P,
where P is the proportion of C derived from wheat residues in the POM fractions between time 0 and 1, δt0 and δt1 are the δ13C values of the POM fraction at time 0 (August 2003) and 1 (April 2004), δw is the average δ13C value of wheat residue (-27.7%), and eq·conct1 is the concentration of POM-C (g C m-2) calculated on an equivalent soil mass basis for time 1.


## 2.4. Statistical Analysis

Variables were tested for non-normality using the D’Agostino-Pearson omnibus test in SAS [33, 34]. Outliers with an absolute value of more than 3 standard deviations from the mean were excluded from the descriptive statistics and geostatistical and multivariate regression analyses [35]. On average, there was less than one outlier for each variable for each sampling time.By accounting for spatial correlation, image maps of POM fractions were produced at each sampling time in Surfer 8 (Golden Software, Golden, CO). Briefly, simple linear regression was first performed on the 64-m regular grid with x and y coordinates as independent variables to identify if there were significant (P<.05) spatial trends in POM fractions. When a significant trend with direction existed, trend values were calculated at all sampling points using the trend model and then subtracted from the measured values. Semivariogram models were defined to determine the spatial structure of the residuals from the trend, and the residuals were then block-kriged and added back to the trends using their respective semivariogram models in GS+ (Gamma Design Software, St. Planinwell, MI). The best model was selected based on visual fit, reduced sums of squares, and r2 of the regression and through cross validation [36]. When there was no significant trend, semivariograms and estimation were made using the measured POM data. Due to a relatively small number (N=36) of samples in each field, anisotropic semivariograms were not considered for each of the data sets [36]. For data sets showing a very weak spatial correlation, point kriging with a linear model was used for interpolation.Repeated measurements over time were used to compare tillage-induced changes in POM-C and -N by the interaction offield (i.e., between the north and south fields) and sampling time (i.e., before versus after potential differences in tillage effects) using a mixed model in SAS [15, 34]. All the POM-C and-N fractions were log-transformed for BACI analysis. Significant main effects of the factors “field” and “sampling time” may indicate the presence of inherent differences and temporal fluctuations across and within the fields that persisted throughout the study. We assumed that these effects of “field” and “sampling time” were not directly controlled by the treatment. Due to difficulty in establishing the control field, we excluded the period of both fields being under MT for the BACI analyses. All BACI model residuals by field or sampling time were checked for independence by visually inspecting histograms and empirical semivariograms, and correlation analysis. A P-value of .02 was used to avoid a type I error.Standardized partial regression coefficients were used to assess the impact of soil texture, BD, and water content on POM. Regression coefficients were obtained by developing a mixed model for each response variable (after log transformation to meet normality) in R using the nlme package [37]. A final model was obtained by backward elimination, starting with a full model that included a full factorial of field, time, sand, clay and BD, minus all interactions of order equal to or greater than four. Location was the random effect. All effects and interactions with a P-value greater than .05 were removed, except when there was a significant interaction involving the effect. Random effects and residuals were checked for normality by visual inspection of quantile plots. Homogeneity of variance was checked by plotting residuals against predicted values. Spatial independence was checked by calculating the empirical semivariograms for residuals and random effects.Due to the unreplicated nature of the tillage treatments, the inference about the effect of ST versus NT on POM was limited to the two paired fields used in this study. When necessary, descriptive statistics were used to compare results between the fields.


## 3. Results and Discussion

### 3.1. Field Level Variability of Soil Properties

The amount of additional C and N that can be sequestered in soil by implementing reduced tillage closely depends on both spatial and temporal variability of soil processes [7]. Thus, by including the spatial and temporal variability of soil C and N stocks as well as soil properties, a more realistic measure of additional C and N changes can be made at the field scale. There was considerable spatial variability among soil C and N and related physical soil properties across the two fields in August 2003 before the tillage conversion (Table 1). The soil heterogeneity of the fields is partially due to the spatial distribution of three soil series across the fields, but also within soil types variability was observed. In the surface zone (0–0.15 m), soil-water content was the most variable among the measured soil properties, with coefficient of variation (CV) values of 27 and 33%, followed by sand content (CV = 20 and 23%). Total C, N, clay content, and BD had less variability (CV = 9 and 18%). Silt content (CV = 7 and 12%) was the least variable and contributed approximately 50% to the total particle mass. Total C and N contents were spatially (P<.05) correlated with clay content and decreasing sand content. Total C and N contents were higher in the north field (1287 g C m-2 and 115.6 g N m-2) than the south field (1100 g C m-2 and 105.8 g N m-2). Within each of the fields, soil-profile C distribution was relatively uniform across the depth of tillage (e.g., 0.15–0.45 m in this study) (Figure 3). On average, the north and south fields contained approximately 1190 and 1000 g C m-2 in the upper 0.5-m depth, respectively. This is a common phenomenon when agricultural soils are subjected to frequent soil disturbance and mixing by intensive tillage. The north field had higher clay (19% at the north field versus 15% at the south field) and silt content (54% versus 48%) and lower sand content (27% versus 37%) than the south field. Water content was higher in the north (0.10 m3m-3) than the south field (0.08 m3m-3). Both fields had a BD of approximately 1.4–1.5 Mg m-3. The measured soil variables were normally distributed at the field scale, except for water content in the south field (data not shown).Table 1
Statistical summary of the surface horizon (0–0.15 m) characteristics in August 2003.
FieldTotal C (g C m-2)Total N (g N m-2)Sand (%)Silt (%)Clay (%)BD (Mg m-3)Water content (m3m-3)NorthMean1287115.627.354.318.61.400.10Minimum99492.820.546.213.20.810.05Maximum1524145.342.759.622.51.880.16CV (%)†910237141827SouthMean1100105.837.048.214.81.470.08Minimum69872.723.337.111.10.980.04Maximum1448128.748.959.620.71.860.14CV (%)17132012161533CV†: coefficient of variation.Figure 3
Total soil carbon content across depths prior to the tillage operations. Bars indicate the standard deviation of the mean.The means for each of the physical properties differed by 9–26% between the fields, except for BD (Table1). Despite the mean differences in the soil properties, the majority of ranges overlapped in both fields, suggesting that the spatial and temporal variability of factors influencing soil organic C and N dynamics may not be different between the fields. Specifically, 61–94% of the measured soil property range for the north field and 53–94% for the south field were similar to each other.


### 3.2. Short-Term Changes in Soil Bulk Density and Water under Tillage Systems

Two months after the tillage conversion in 2003, the mean BD values at the 0.15 m depth were higher in the south field under NT (1.24 Mg m-3) than in the north field under ST (1.00 Mg m-3) (Figure 4). The BD then increased in ST to almost the same value as in the NT. These results are in accordance with the often observed large temporal changes of BD that could result in insignificant tillage effects on BD [31]. Soil BD is often more influenced by soil properties, climatic conditions and their seasonal variations rather than tillage practices [25, 38–41]. Overall, the response of surface horizon BD to tillage treatments greatly depends on time of determination and tillage intensity [42].Seasonal changes in bulk density and volumetric water content. Arrows indicate the time of the tillage operations. The north field was converted from no-till to standard tillage in August 2003, while the south field remained under no-till. Both fields were then converted to minimum tillage in May 2005. Bars indicate the standard deviation of the mean.
(a)(b)Soil-water conditions are closely associated with changes in soil biochemical properties (e.g., mineral N) over space and time [43, 44]. Lee et al. [43] showed that soil-water conditions, which were regulated by tillage and irrigation, could affect the bioavailability and type of soil C and N sources for microbial activity in irrigated agricultural systems of California, hence continuously altering SOM mineralization [45, 46]. Soil-water content prior to the tillage conversion was greater in the north field than the south field (Figure 4). Following tillage in the ST field, soil-water content increased for both fields during the rainy season (November 2003~March 2004) (Figure 4). Over the period of this study, there was no effect of tillage on soil water content. Consequently, the water content differences between the fields under ST and NT were not large enough and consistent in direction to suggest that conversion to reduced tillage can lead to short-term gains in water conservation as has been observed in the Great Plains and the Midwest [47, 48].


### 3.3. Spatial Variation in Particulate Organic Matter

Prior to the tillage conversion, soil POM exhibited large spatial variability within and across the fields (Table2 and Figure 5) and total POM (53–2000 μm) accounted for approximately 21% and 14% of total C and N stocks, respectively, in the top 15 cm of the soils (Tables 1 and 2). The combined C and N stocks in the POM 53–250 μm and 250–1000 μm fractions accounted for 94% and 83% of total POM-C and -N across the fields, respectively. Despite high spatial variability of POM, C and N stocks of all size classes of POM decreased across the north field following the tillage conversion, while there was only little change in the south field with time under continuous NT (Table 2 and Figure 5). Others have also reported that POM levels generally decreased with increasing tillage intensity [49, 50]. This suggests that the spatial distribution of POM stocks in the fields was in part affected by tillage at the field scale. With a subsequent conversion to MT in both fields, little or no change in spatial distribution of POM-C stocks occurred. Image maps of the POM-N fractions exhibited a similar spatial pattern with respect to the tillage treatments (data not shown). Regardless of the tillage treatments, the coefficient of variation of POM-C and -N increased with increasing size fraction because plant residue input was dominantly incorporated into coarse POM compared to fine POM [51].Table 2
Descriptive statistics for soil particulate organic matter (POM)-C and -N fractions, and C to N ratio in the 0–15 cm depth.
POM-C (g C m-2)POM-N (g N m-2)C to N ratioFractionFieldTimeAug. 03Apr. 04Oct. 04June 05Aug. 03Apr. 04Oct. 04June 05Aug. 03Apr. 04Oct. 04June 05>1000  μmNorthMean14.811.320.29.50.740.370.600.4821.232.835.421.2Minimum4.03.13.23.10.130.120.200.209.315.311.09.0Maximum41.241.755.817.02.271.231.571.1933.274.095.357.4CV†556965446268455023407142SouthMean15.921.225.423.80.761.021.191.3519.521.120.318.4Minimum3.03.83.97.50.120.170.210.577.511.19.37.3Maximum49.154.556.946.12.752.182.153.0533.834.931.032.5CV725853496852414432272833250–1000μmNorthMean85.054.754.853.14.632.262.682.9018.824.820.719.0Minimum47.427.925.430.82.621.041.521.7615.315.615.315.5Maximum155.291.680.681.17.845.494.005.1323.132.236.522.9CV33302520323722269132311SouthMean82.089.883.298.63.974.094.255.3021.522.119.718.4Minimum45.336.340.547.31.761.551.973.0214.016.114.111.9Maximum164.5186.6131.7179.17.147.856.698.5130.331.027.123.9CV32392734363524282116161453–250μmNorthMean157.8118.4117.7105.910.556.647.367.9415.118.115.813.3Minimum120.887.076.572.77.214.305.225.6512.513.713.111.9Maximum200.1150.1155.1136.015.0610.949.2210.5219.422.720.414.2CV13141415172114131110104SouthMean137.8144.6139.6149.19.878.778.5910.9014.416.816.213.6Minimum81.087.692.295.46.115.175.367.3611.014.312.510.9Maximum227.0211.4235.1231.316.7312.9011.5615.1622.521.020.715.6CV25202321312017171911128CV†: coefficient of variation.Figure 5
Distribution of soil particulate organic matter C (POM-C) for standard tillage (ST), no-till (NT), and minimum tillage (MT) in the 0–0.15 m depth from August 2003 to June 2005. MT+ indicates one more bed disc (to 0.15 m) pass on the north field.The POM C/N ratios were on average slightly higher in the north field than those in the south field, mostly due to different amount of recently incorporated wheat residues (Table2). In 2003, the C/N ratio of the wheat residue was 86±7 (mean ± standard error) (data not shown). Across the fields and sampling times, the C/N ratios ranged 13–18, 18–25, and 18–35 for the 53–250, 250–1000, and 1000–2000 μm fractions, respectively. This tendency of increasing POM C/N ratios with size suggests that fine POM was more decomposed than coarse POM [24, 52, 53]. There was no measureable enrichment of C relative to N in POM under NT when compared to ST, which is consistent with the results of Cambardella and Elliott [22].


### 3.4. Temporal Variation in Particulate Organic Matter

After the tillage conversion in 2003, C and N stocks of all the POM fractions were significantly higher in NT than ST at the two subsequent sampling times (Tables2, 3 and 4). We found a consistent interactive effect of field by sampling time on both POM-C and -N. Contrast analyses suggest that the differences (P<.01) between tillage systems were significant for all POM fractions in April and October 2004 and the overall 12 months period after the tillage conversion. The major differences in POM fractions between tillage practices appeared between August 2003 and April 2004. In the north field, the average C and N stocks of the POM 53–250, 250–1000, and 1000–2000 μm fractions decreased by approximately 25, 36, and 23% for C and 37, 51, and 50% for N (Table 2). In contrast, both C (≈5%) and N stocks (≈10%) of the POM 53–250 μm and 250–1000 μm fractions slightly increased in the south field, although the magnitude of change was significantly lower than for the north field under ST. Therefore, the contrasting effects of ST and NT on the 250–1000 μm fractions were more evident than in the 53–250 μm fractions. The differences in POM persisted without much change thereafter, confirming preferential loss of POM at the time of intensive tillage management [50, 54]. Overall, the temporal differences in POM-C between the fields induced by tillage were on average 7.8, 34.7, and 44.1 g m-2 for the 1000–2000, 250–1000, and 53–250 μm fractions. The differences in POM-N were 0.66, 2.35, and 2.36 g m-2 for the 1000–2000, 250–1000, and 53–250 μm fractions. In addition, significant field or time effects indicate that the tillage effects are scale-dependent (Table 3).Table 3
Effect of standard tillage and no-till on particulate organic matter (POM)-C (BACI† analysis). All the POM-C fractions were log transformed.
FractionSourceDFSum of squaresMean squareF valuePr>F>1000 μmField11.0421.04212.940.001Time21.2740.6378.800.000Field× time20.6780.3394.690.011Location (field)705.6390.0811.110.294Residual14010.1350.072Contrast‡T0-T119.300.003T0-T211.670.199Overall means16.280.013250–1000μmField10.8320.83237.86<0.0001Time20.3510.17510.32<0.0001Field× time20.5840.29217.19<0.0001Location (field)701.5390.0221.290.100Residual1402.3790.017ContrastT0-T1127.71<0.0001T0-T2123.70<0.0001Overall means134.22<0.000153–250μmField10.0480.0484.030.0487Time20.1450.07315.97<0.0001Field× time20.2600.13028.63<0.0001Location (field)700.8310.0122.61<0.0001Residual1400.6360.005ContrastT0-T1148.79<0.0001T0-T2136.15<0.0001Overall means156.31<0.0001BACI†: Before-After and Control-Impact model.Difference‡ in the interaction of field and time between August 2003 (T0), April (T1), and October 2004 (T2).Table 4
Effect of standard tillage and no-till on particulate organic matter (POM)-N (BACI† analysis). All the POM-N fractions were log transformed.
FractionSourceDFSum of squaresMean squareF valuePr>F>1000 μmField13.9783.97839.33<0.0001Time21.2840.64212.56<0.0001Field× time21.5100.75514.77<0.0001Location (field)707.0870.1011.980.000Residual1407.1060.051Contrast‡T0-T1129.17<0.0001T0-T2111.080.001Overall means125.41<0.0001250–1000μmField10.9200.92040.65<0.0001Time20.8340.41725.22<0.0001Field× time21.1440.57234.61<0.0001Location (field)701.5840.0231.370.060Residual1402.3140.017ContrastT0-T1159.75<0.0001T0-T2142.63<0.0001Overall means167.77<0.000153–250μmField10.1150.1159.430.003Time20.6040.30247.87<0.0001Field× time20.2450.12319.44<0.0001Location (field)700.8540.0121.930.0005Residual1400.8820.006ContrastT0-T1138.16<0.0001T0-T2114.640.000Overall means133.35<0.0001BACI†: Before-After and Control-Impact model.Difference‡ in the interaction of field and time between August 2003 (T0), April (T1), and October 2004 (T2).The results show that the mean differences over time in the amount of each POM fraction were significantly lower in the north than in the south field, due to the level of tillage intensity and time since tillage. Coarse POM was more sensitive to tillage than fine POM, but the largest C and N amounts were lost in the fine POM fraction as soon as tillage operations started. This suggests that fine POM is an important pool of soil C and N sensitive to soil disturbance. No-till and MT appeared to maintain the level of soil POM relative to the level under ST, suggesting a capacity of reduced tillage in relative to ST to delay the loss of labile organic fractions under the Mediterranean climate conditions. We also found no major short-term (<3 years) declines in C and N of the non-POM fraction (equivalent to mineral-associated SOM <53 μm) in the north field when compared with the south field (Tables 2 and 4), leading to a greater proportion of non-POM in bulk SOM under ST than NT [49, 50, 53].


### 3.5. New C Input to Soil in Standard and No-Tillage

Similar to POM-C, the concentration of new C input derived from wheat residue into each POM fraction had substantial spatial variability across the fields, with a range of CVs from 53 to 79% (Table5). There were minor differences in the concentration of new POM-C derived from wheat residue between the fields, at least in the short term However, differences in δ13C in all the POM fractions over time (Figure 6) and the proportion of wheat-derived C/POM-C clearly showed faster turnover of new POM-C after seven months under ST than under NT. Even though more wheat residue C was probably incorporated in the north field at the time of ST, our results indicate that simultaneously less residue-derived C accumulated as POM under ST than under NT. In contrast, NT had relatively high fine POM contents, primarily due to the lack of soil disturbance and increased residue cover [5].Table 5
Wheat-derived C and proportion of wheat-derived C/total fraction C in POM fractions under standard tillage and no-till in the first eight months of study.
Standard tillageNo-tillFraction>1000 μm250–1000μm53–250μm>1000 μm250–1000μm53–250μmSoil POM-Cnew (g C m-2)Mean3.38.014.02.611.310.1Maximum8.527.628.87.029.428.8Minimum0.10.50.30.20.30.2CV†657353796075Skewness‡0.401.430.220.970.670.95Soil POM-Cnew/POM-C (%)Mean32.915.212.417.514.97.3Maximum66.632.223.649.233.119.2Minimum1.71.90.20.50.60.2CV566052805271Skewness-0.100.40-0.110.840.270.65CV†: coefficient of variation (%).Skewness‡ indicates the degree of asymmetry of a distribution.Short-term effect of tillage onδ13C particulate organic matter (POM). Arrows indicate the time of tillage operations. The north field was converted from no-till to standard tillage in August 2003, while the south field remained under no-till. Both fields were then converted to minimum tillage in May 2005. Bars indicate the standard deviation of the mean.
(a)(b)(c)Mean differences in new C input between ST and NT accounted for only 0–8.5% of the tillage-induced differences in POM-C (Tables2 and 5). Therefore, the decomposition of old POM-C by tillage primarily led to total POM changes. Independent of the tillage treatments, the relative proportion of new C incorporated into each POM fraction decreased in the following order: 1000–2000 μm > 250–1000 μm > 53–250 μm. Therefore, NT also reduced the decomposition rate of pre-existing labile C and N in smaller POM fractions.


### 3.6. Field Level Soil Controls on Particulate Organic Matter

#### 3.6.1. Soil Texture

Sand was a significant predictor of the coarse POM fractions, whereas clay was considered significant for controlling the fine POM fraction (Table6). Overall, sand content was positively related to the 1000–2000 μm POM fraction (Table 7). We speculate that this is because soils with high sand content could store more crop residues as coarse POM than soils with low sand content. Clay may play a role in sequestering partly decomposed POM through physical and chemical protection of soil C and N [1, 49, 55], probably causing different responses of the fine POM fraction to varying tillage intensity and texture. However, the fine POM fraction in the soil surface significantly decreased with clay content across the fields (Table 7). This could be explained by relatively low clay contents (<20%) of the fields that have a tendency of fast degradation of POM after tillage is introduced.Table 6
Significance of effects of soil physical properties on particulate organic matter (POM) C and N fractions, accounting for effects of field and time. All POM fraction values were log-transformed.
Effect>1000 μm POM250–1000μm POM53–250μm POMPOM-CP-valueField0.0001.00010.0001Time0.0001.00010.0001Sand0.0416.51010.0801Clay0.1504.35040.0010Bulk density0.0010.00010.0001Field× time0.0001.00010.0001Time× clay0.0043.01550.7712Time× bulk density0.0196.00160.0399Std. Dev. locationNA†.09690.0935Std. Dev. residual0.5762.27170.1499POM-NP-valueField0.0001.00010.0001Time0.0001.00010.0001Sand0.0001.0150.8818Clay0.1786.09120.0004Bulk density0.1767.00110.0001Field× time0.0001.00010.0001Field× clay–‡–0.2774Time× clay––0.7443Field×time×clay––0.0017Std. Dev. location0.1858.09320.0814Std. Dev. residual0.5061.27490.1631NA†: not applicable. The likelihood ratio test indicated that the random effect for location was not necessary and thus it was removed from the model according to Pinheiro and Bates (2000).The‡ effect was removed from the model because it was not significant at the level of 2%.Table 7
Standardized partial regression coefficients relating soil physical properties to particulate organic matter (POM)-C and -N fractions. Numbers in the columns labeled 98%CI are the half-width of the confidence intervals.
FractionTimeSand98%CIClay98%CIBulk density98%CIPOM-C>1000μmAugust 030.030.339-0.100.221April 04-0.230.355-0.480.504October 040.370.343-0.530.497June 050.130.354-0.690.473All0.270.237250–1000μmAugust 03-0.240.321-0.190.196April 04-0.450.334-0.870.447October 040.080.324-0.360.441June 05-0.160.334-0.760.419All-0.010.24753–250μmAugust 03-0.480.351-0.240.188April 04-0.620.362-0.660.431October 04-0.390.353-0.200.425June 05-0.550.362-0.640.402All-0.470.297POM-N>1000μmAll0.220.246-0.160.244-0.120.504250–1000μmAll0.020.23-0.210.228-0.250.27253–250μmAll-0.210.261-0.300.263-0.160.150For the coarse POM-C fractions, studies suggest that the degree of POM dependency on soil texture seems to be regulated by the level to which soil texture varies spatially between the fields with ST, NT, and MT [25, 56]. In addition, the relationships of sand and clay contents with other biotic and abiotic factors were largely affected by time and tillage [43]. As a result, this would significantly confound the dynamics of coarse POM. Our results are consistent with Wander et al. [40], showing that the relationship between soil texture and POM changed inconsistently over time. However, it was not evident what led to this shift in their relationship. In contrast, the relationship between the 53–250 μm POM-C fraction and clay appeared to be temporally consistent at the field scale. Therefore, the effect of clay on the dynamics of fine POM would be predominantly confounded by tillage and presumably independent of spatiotemporal variability in the other factors. Similarly, the effect of sand was positive or zero for the 250–2000 μm POM-N fractions and the effect of clay was consistently negative for the 53–250 μm POM-N fraction.


#### 3.6.2. Soil Bulk Density

The BD was a significant predictor for all POM fractions except for the POM-N> 1000 μm (Table 6). Regardless of the tillage treatments, the BD had a consistently negative effect on POM-C in all fractions at the field scale. Although the effect of BD on POM was greater but more variable than the effects of sand and clay in the majority of the cases, no significant differences were evident (Table 7). Tillage-induced changes in BD could be in part related to POM dynamics. For example, da Silva et al. [39] and Wander et al. [40] showed that there was a trend for decreasing bulk SOM with increasing BD under cultivation. Ball et al. [57] showed that maximum BD was negatively related to water content, total porosity, and clay content. We speculate that the rate of residue decomposition and new POM formation would be limited under lower soil water availability due to a high BD, leading to decreased POM stocks. However, we were not able to find the evidence to support this speculation. The effects of BD on the POM-N fractions were similar to those of the POM-C, but slopes tended to be more constant over time.


#### 3.6.3. Soil Water

The effect of soil water on POM was assessed by adding this factor to the reduced models obtained by backward elimination. Soil water was not significant for any of the response variables. This may be because water content is a too dynamic variable in comparison to POM. Although SOM mineralization can be greatly enhanced with increasing water content [45, 46], the spatial variation in soil conditions (e.g., soil texture) is likely responsible for variable and complex changes in soil chemical and biological properties [44]. Spatial patterns of such soil properties interacting with water were likely highly confounded by tillage.


## 3.1. Field Level Variability of Soil Properties

The amount of additional C and N that can be sequestered in soil by implementing reduced tillage closely depends on both spatial and temporal variability of soil processes [7]. Thus, by including the spatial and temporal variability of soil C and N stocks as well as soil properties, a more realistic measure of additional C and N changes can be made at the field scale. There was considerable spatial variability among soil C and N and related physical soil properties across the two fields in August 2003 before the tillage conversion (Table 1). The soil heterogeneity of the fields is partially due to the spatial distribution of three soil series across the fields, but also within soil types variability was observed. In the surface zone (0–0.15 m), soil-water content was the most variable among the measured soil properties, with coefficient of variation (CV) values of 27 and 33%, followed by sand content (CV = 20 and 23%). Total C, N, clay content, and BD had less variability (CV = 9 and 18%). Silt content (CV = 7 and 12%) was the least variable and contributed approximately 50% to the total particle mass. Total C and N contents were spatially (P<.05) correlated with clay content and decreasing sand content. Total C and N contents were higher in the north field (1287 g C m-2 and 115.6 g N m-2) than the south field (1100 g C m-2 and 105.8 g N m-2). Within each of the fields, soil-profile C distribution was relatively uniform across the depth of tillage (e.g., 0.15–0.45 m in this study) (Figure 3). On average, the north and south fields contained approximately 1190 and 1000 g C m-2 in the upper 0.5-m depth, respectively. This is a common phenomenon when agricultural soils are subjected to frequent soil disturbance and mixing by intensive tillage. The north field had higher clay (19% at the north field versus 15% at the south field) and silt content (54% versus 48%) and lower sand content (27% versus 37%) than the south field. Water content was higher in the north (0.10 m3m-3) than the south field (0.08 m3m-3). Both fields had a BD of approximately 1.4–1.5 Mg m-3. The measured soil variables were normally distributed at the field scale, except for water content in the south field (data not shown).Table 1
Statistical summary of the surface horizon (0–0.15 m) characteristics in August 2003.
FieldTotal C (g C m-2)Total N (g N m-2)Sand (%)Silt (%)Clay (%)BD (Mg m-3)Water content (m3m-3)NorthMean1287115.627.354.318.61.400.10Minimum99492.820.546.213.20.810.05Maximum1524145.342.759.622.51.880.16CV (%)†910237141827SouthMean1100105.837.048.214.81.470.08Minimum69872.723.337.111.10.980.04Maximum1448128.748.959.620.71.860.14CV (%)17132012161533CV†: coefficient of variation.Figure 3
Total soil carbon content across depths prior to the tillage operations. Bars indicate the standard deviation of the mean.The means for each of the physical properties differed by 9–26% between the fields, except for BD (Table1). Despite the mean differences in the soil properties, the majority of ranges overlapped in both fields, suggesting that the spatial and temporal variability of factors influencing soil organic C and N dynamics may not be different between the fields. Specifically, 61–94% of the measured soil property range for the north field and 53–94% for the south field were similar to each other.


## 3.2. Short-Term Changes in Soil Bulk Density and Water under Tillage Systems

Two months after the tillage conversion in 2003, the mean BD values at the 0.15 m depth were higher in the south field under NT (1.24 Mg m-3) than in the north field under ST (1.00 Mg m-3) (Figure 4). The BD then increased in ST to almost the same value as in the NT. These results are in accordance with the often observed large temporal changes of BD that could result in insignificant tillage effects on BD [31]. Soil BD is often more influenced by soil properties, climatic conditions and their seasonal variations rather than tillage practices [25, 38–41]. Overall, the response of surface horizon BD to tillage treatments greatly depends on time of determination and tillage intensity [42].Seasonal changes in bulk density and volumetric water content. Arrows indicate the time of the tillage operations. The north field was converted from no-till to standard tillage in August 2003, while the south field remained under no-till. Both fields were then converted to minimum tillage in May 2005. Bars indicate the standard deviation of the mean.
(a)(b)Soil-water conditions are closely associated with changes in soil biochemical properties (e.g., mineral N) over space and time [43, 44]. Lee et al. [43] showed that soil-water conditions, which were regulated by tillage and irrigation, could affect the bioavailability and type of soil C and N sources for microbial activity in irrigated agricultural systems of California, hence continuously altering SOM mineralization [45, 46]. Soil-water content prior to the tillage conversion was greater in the north field than the south field (Figure 4). Following tillage in the ST field, soil-water content increased for both fields during the rainy season (November 2003~March 2004) (Figure 4). Over the period of this study, there was no effect of tillage on soil water content. Consequently, the water content differences between the fields under ST and NT were not large enough and consistent in direction to suggest that conversion to reduced tillage can lead to short-term gains in water conservation as has been observed in the Great Plains and the Midwest [47, 48].


## 3.3. Spatial Variation in Particulate Organic Matter

Prior to the tillage conversion, soil POM exhibited large spatial variability within and across the fields (Table2 and Figure 5) and total POM (53–2000 μm) accounted for approximately 21% and 14% of total C and N stocks, respectively, in the top 15 cm of the soils (Tables 1 and 2). The combined C and N stocks in the POM 53–250 μm and 250–1000 μm fractions accounted for 94% and 83% of total POM-C and -N across the fields, respectively. Despite high spatial variability of POM, C and N stocks of all size classes of POM decreased across the north field following the tillage conversion, while there was only little change in the south field with time under continuous NT (Table 2 and Figure 5). Others have also reported that POM levels generally decreased with increasing tillage intensity [49, 50]. This suggests that the spatial distribution of POM stocks in the fields was in part affected by tillage at the field scale. With a subsequent conversion to MT in both fields, little or no change in spatial distribution of POM-C stocks occurred. Image maps of the POM-N fractions exhibited a similar spatial pattern with respect to the tillage treatments (data not shown). Regardless of the tillage treatments, the coefficient of variation of POM-C and -N increased with increasing size fraction because plant residue input was dominantly incorporated into coarse POM compared to fine POM [51].Table 2
Descriptive statistics for soil particulate organic matter (POM)-C and -N fractions, and C to N ratio in the 0–15 cm depth.
POM-C (g C m-2)POM-N (g N m-2)C to N ratioFractionFieldTimeAug. 03Apr. 04Oct. 04June 05Aug. 03Apr. 04Oct. 04June 05Aug. 03Apr. 04Oct. 04June 05>1000  μmNorthMean14.811.320.29.50.740.370.600.4821.232.835.421.2Minimum4.03.13.23.10.130.120.200.209.315.311.09.0Maximum41.241.755.817.02.271.231.571.1933.274.095.357.4CV†556965446268455023407142SouthMean15.921.225.423.80.761.021.191.3519.521.120.318.4Minimum3.03.83.97.50.120.170.210.577.511.19.37.3Maximum49.154.556.946.12.752.182.153.0533.834.931.032.5CV725853496852414432272833250–1000μmNorthMean85.054.754.853.14.632.262.682.9018.824.820.719.0Minimum47.427.925.430.82.621.041.521.7615.315.615.315.5Maximum155.291.680.681.17.845.494.005.1323.132.236.522.9CV33302520323722269132311SouthMean82.089.883.298.63.974.094.255.3021.522.119.718.4Minimum45.336.340.547.31.761.551.973.0214.016.114.111.9Maximum164.5186.6131.7179.17.147.856.698.5130.331.027.123.9CV32392734363524282116161453–250μmNorthMean157.8118.4117.7105.910.556.647.367.9415.118.115.813.3Minimum120.887.076.572.77.214.305.225.6512.513.713.111.9Maximum200.1150.1155.1136.015.0610.949.2210.5219.422.720.414.2CV13141415172114131110104SouthMean137.8144.6139.6149.19.878.778.5910.9014.416.816.213.6Minimum81.087.692.295.46.115.175.367.3611.014.312.510.9Maximum227.0211.4235.1231.316.7312.9011.5615.1622.521.020.715.6CV25202321312017171911128CV†: coefficient of variation.Figure 5
Distribution of soil particulate organic matter C (POM-C) for standard tillage (ST), no-till (NT), and minimum tillage (MT) in the 0–0.15 m depth from August 2003 to June 2005. MT+ indicates one more bed disc (to 0.15 m) pass on the north field.The POM C/N ratios were on average slightly higher in the north field than those in the south field, mostly due to different amount of recently incorporated wheat residues (Table2). In 2003, the C/N ratio of the wheat residue was 86±7 (mean ± standard error) (data not shown). Across the fields and sampling times, the C/N ratios ranged 13–18, 18–25, and 18–35 for the 53–250, 250–1000, and 1000–2000 μm fractions, respectively. This tendency of increasing POM C/N ratios with size suggests that fine POM was more decomposed than coarse POM [24, 52, 53]. There was no measureable enrichment of C relative to N in POM under NT when compared to ST, which is consistent with the results of Cambardella and Elliott [22].


## 3.4. Temporal Variation in Particulate Organic Matter

After the tillage conversion in 2003, C and N stocks of all the POM fractions were significantly higher in NT than ST at the two subsequent sampling times (Tables2, 3 and 4). We found a consistent interactive effect of field by sampling time on both POM-C and -N. Contrast analyses suggest that the differences (P<.01) between tillage systems were significant for all POM fractions in April and October 2004 and the overall 12 months period after the tillage conversion. The major differences in POM fractions between tillage practices appeared between August 2003 and April 2004. In the north field, the average C and N stocks of the POM 53–250, 250–1000, and 1000–2000 μm fractions decreased by approximately 25, 36, and 23% for C and 37, 51, and 50% for N (Table 2). In contrast, both C (≈5%) and N stocks (≈10%) of the POM 53–250 μm and 250–1000 μm fractions slightly increased in the south field, although the magnitude of change was significantly lower than for the north field under ST. Therefore, the contrasting effects of ST and NT on the 250–1000 μm fractions were more evident than in the 53–250 μm fractions. The differences in POM persisted without much change thereafter, confirming preferential loss of POM at the time of intensive tillage management [50, 54]. Overall, the temporal differences in POM-C between the fields induced by tillage were on average 7.8, 34.7, and 44.1 g m-2 for the 1000–2000, 250–1000, and 53–250 μm fractions. The differences in POM-N were 0.66, 2.35, and 2.36 g m-2 for the 1000–2000, 250–1000, and 53–250 μm fractions. In addition, significant field or time effects indicate that the tillage effects are scale-dependent (Table 3).Table 3
Effect of standard tillage and no-till on particulate organic matter (POM)-C (BACI† analysis). All the POM-C fractions were log transformed.
FractionSourceDFSum of squaresMean squareF valuePr>F>1000 μmField11.0421.04212.940.001Time21.2740.6378.800.000Field× time20.6780.3394.690.011Location (field)705.6390.0811.110.294Residual14010.1350.072Contrast‡T0-T119.300.003T0-T211.670.199Overall means16.280.013250–1000μmField10.8320.83237.86<0.0001Time20.3510.17510.32<0.0001Field× time20.5840.29217.19<0.0001Location (field)701.5390.0221.290.100Residual1402.3790.017ContrastT0-T1127.71<0.0001T0-T2123.70<0.0001Overall means134.22<0.000153–250μmField10.0480.0484.030.0487Time20.1450.07315.97<0.0001Field× time20.2600.13028.63<0.0001Location (field)700.8310.0122.61<0.0001Residual1400.6360.005ContrastT0-T1148.79<0.0001T0-T2136.15<0.0001Overall means156.31<0.0001BACI†: Before-After and Control-Impact model.Difference‡ in the interaction of field and time between August 2003 (T0), April (T1), and October 2004 (T2).Table 4
Effect of standard tillage and no-till on particulate organic matter (POM)-N (BACI† analysis). All the POM-N fractions were log transformed.
FractionSourceDFSum of squaresMean squareF valuePr>F>1000 μmField13.9783.97839.33<0.0001Time21.2840.64212.56<0.0001Field× time21.5100.75514.77<0.0001Location (field)707.0870.1011.980.000Residual1407.1060.051Contrast‡T0-T1129.17<0.0001T0-T2111.080.001Overall means125.41<0.0001250–1000μmField10.9200.92040.65<0.0001Time20.8340.41725.22<0.0001Field× time21.1440.57234.61<0.0001Location (field)701.5840.0231.370.060Residual1402.3140.017ContrastT0-T1159.75<0.0001T0-T2142.63<0.0001Overall means167.77<0.000153–250μmField10.1150.1159.430.003Time20.6040.30247.87<0.0001Field× time20.2450.12319.44<0.0001Location (field)700.8540.0121.930.0005Residual1400.8820.006ContrastT0-T1138.16<0.0001T0-T2114.640.000Overall means133.35<0.0001BACI†: Before-After and Control-Impact model.Difference‡ in the interaction of field and time between August 2003 (T0), April (T1), and October 2004 (T2).The results show that the mean differences over time in the amount of each POM fraction were significantly lower in the north than in the south field, due to the level of tillage intensity and time since tillage. Coarse POM was more sensitive to tillage than fine POM, but the largest C and N amounts were lost in the fine POM fraction as soon as tillage operations started. This suggests that fine POM is an important pool of soil C and N sensitive to soil disturbance. No-till and MT appeared to maintain the level of soil POM relative to the level under ST, suggesting a capacity of reduced tillage in relative to ST to delay the loss of labile organic fractions under the Mediterranean climate conditions. We also found no major short-term (<3 years) declines in C and N of the non-POM fraction (equivalent to mineral-associated SOM <53 μm) in the north field when compared with the south field (Tables 2 and 4), leading to a greater proportion of non-POM in bulk SOM under ST than NT [49, 50, 53].


## 3.5. New C Input to Soil in Standard and No-Tillage

Similar to POM-C, the concentration of new C input derived from wheat residue into each POM fraction had substantial spatial variability across the fields, with a range of CVs from 53 to 79% (Table5). There were minor differences in the concentration of new POM-C derived from wheat residue between the fields, at least in the short term However, differences in δ13C in all the POM fractions over time (Figure 6) and the proportion of wheat-derived C/POM-C clearly showed faster turnover of new POM-C after seven months under ST than under NT. Even though more wheat residue C was probably incorporated in the north field at the time of ST, our results indicate that simultaneously less residue-derived C accumulated as POM under ST than under NT. In contrast, NT had relatively high fine POM contents, primarily due to the lack of soil disturbance and increased residue cover [5].Table 5
Wheat-derived C and proportion of wheat-derived C/total fraction C in POM fractions under standard tillage and no-till in the first eight months of study.
Standard tillageNo-tillFraction>1000 μm250–1000μm53–250μm>1000 μm250–1000μm53–250μmSoil POM-Cnew (g C m-2)Mean3.38.014.02.611.310.1Maximum8.527.628.87.029.428.8Minimum0.10.50.30.20.30.2CV†657353796075Skewness‡0.401.430.220.970.670.95Soil POM-Cnew/POM-C (%)Mean32.915.212.417.514.97.3Maximum66.632.223.649.233.119.2Minimum1.71.90.20.50.60.2CV566052805271Skewness-0.100.40-0.110.840.270.65CV†: coefficient of variation (%).Skewness‡ indicates the degree of asymmetry of a distribution.Short-term effect of tillage onδ13C particulate organic matter (POM). Arrows indicate the time of tillage operations. The north field was converted from no-till to standard tillage in August 2003, while the south field remained under no-till. Both fields were then converted to minimum tillage in May 2005. Bars indicate the standard deviation of the mean.
(a)(b)(c)Mean differences in new C input between ST and NT accounted for only 0–8.5% of the tillage-induced differences in POM-C (Tables2 and 5). Therefore, the decomposition of old POM-C by tillage primarily led to total POM changes. Independent of the tillage treatments, the relative proportion of new C incorporated into each POM fraction decreased in the following order: 1000–2000 μm > 250–1000 μm > 53–250 μm. Therefore, NT also reduced the decomposition rate of pre-existing labile C and N in smaller POM fractions.


## 3.6. Field Level Soil Controls on Particulate Organic Matter

### 3.6.1. Soil Texture

Sand was a significant predictor of the coarse POM fractions, whereas clay was considered significant for controlling the fine POM fraction (Table6). Overall, sand content was positively related to the 1000–2000 μm POM fraction (Table 7). We speculate that this is because soils with high sand content could store more crop residues as coarse POM than soils with low sand content. Clay may play a role in sequestering partly decomposed POM through physical and chemical protection of soil C and N [1, 49, 55], probably causing different responses of the fine POM fraction to varying tillage intensity and texture. However, the fine POM fraction in the soil surface significantly decreased with clay content across the fields (Table 7). This could be explained by relatively low clay contents (<20%) of the fields that have a tendency of fast degradation of POM after tillage is introduced.Table 6
Significance of effects of soil physical properties on particulate organic matter (POM) C and N fractions, accounting for effects of field and time. All POM fraction values were log-transformed.
Effect>1000 μm POM250–1000μm POM53–250μm POMPOM-CP-valueField0.0001.00010.0001Time0.0001.00010.0001Sand0.0416.51010.0801Clay0.1504.35040.0010Bulk density0.0010.00010.0001Field× time0.0001.00010.0001Time× clay0.0043.01550.7712Time× bulk density0.0196.00160.0399Std. Dev. locationNA†.09690.0935Std. Dev. residual0.5762.27170.1499POM-NP-valueField0.0001.00010.0001Time0.0001.00010.0001Sand0.0001.0150.8818Clay0.1786.09120.0004Bulk density0.1767.00110.0001Field× time0.0001.00010.0001Field× clay–‡–0.2774Time× clay––0.7443Field×time×clay––0.0017Std. Dev. location0.1858.09320.0814Std. Dev. residual0.5061.27490.1631NA†: not applicable. The likelihood ratio test indicated that the random effect for location was not necessary and thus it was removed from the model according to Pinheiro and Bates (2000).The‡ effect was removed from the model because it was not significant at the level of 2%.Table 7
Standardized partial regression coefficients relating soil physical properties to particulate organic matter (POM)-C and -N fractions. Numbers in the columns labeled 98%CI are the half-width of the confidence intervals.
FractionTimeSand98%CIClay98%CIBulk density98%CIPOM-C>1000μmAugust 030.030.339-0.100.221April 04-0.230.355-0.480.504October 040.370.343-0.530.497June 050.130.354-0.690.473All0.270.237250–1000μmAugust 03-0.240.321-0.190.196April 04-0.450.334-0.870.447October 040.080.324-0.360.441June 05-0.160.334-0.760.419All-0.010.24753–250μmAugust 03-0.480.351-0.240.188April 04-0.620.362-0.660.431October 04-0.390.353-0.200.425June 05-0.550.362-0.640.402All-0.470.297POM-N>1000μmAll0.220.246-0.160.244-0.120.504250–1000μmAll0.020.23-0.210.228-0.250.27253–250μmAll-0.210.261-0.300.263-0.160.150For the coarse POM-C fractions, studies suggest that the degree of POM dependency on soil texture seems to be regulated by the level to which soil texture varies spatially between the fields with ST, NT, and MT [25, 56]. In addition, the relationships of sand and clay contents with other biotic and abiotic factors were largely affected by time and tillage [43]. As a result, this would significantly confound the dynamics of coarse POM. Our results are consistent with Wander et al. [40], showing that the relationship between soil texture and POM changed inconsistently over time. However, it was not evident what led to this shift in their relationship. In contrast, the relationship between the 53–250 μm POM-C fraction and clay appeared to be temporally consistent at the field scale. Therefore, the effect of clay on the dynamics of fine POM would be predominantly confounded by tillage and presumably independent of spatiotemporal variability in the other factors. Similarly, the effect of sand was positive or zero for the 250–2000 μm POM-N fractions and the effect of clay was consistently negative for the 53–250 μm POM-N fraction.


### 3.6.2. Soil Bulk Density

The BD was a significant predictor for all POM fractions except for the POM-N> 1000 μm (Table 6). Regardless of the tillage treatments, the BD had a consistently negative effect on POM-C in all fractions at the field scale. Although the effect of BD on POM was greater but more variable than the effects of sand and clay in the majority of the cases, no significant differences were evident (Table 7). Tillage-induced changes in BD could be in part related to POM dynamics. For example, da Silva et al. [39] and Wander et al. [40] showed that there was a trend for decreasing bulk SOM with increasing BD under cultivation. Ball et al. [57] showed that maximum BD was negatively related to water content, total porosity, and clay content. We speculate that the rate of residue decomposition and new POM formation would be limited under lower soil water availability due to a high BD, leading to decreased POM stocks. However, we were not able to find the evidence to support this speculation. The effects of BD on the POM-N fractions were similar to those of the POM-C, but slopes tended to be more constant over time.


### 3.6.3. Soil Water

The effect of soil water on POM was assessed by adding this factor to the reduced models obtained by backward elimination. Soil water was not significant for any of the response variables. This may be because water content is a too dynamic variable in comparison to POM. Although SOM mineralization can be greatly enhanced with increasing water content [45, 46], the spatial variation in soil conditions (e.g., soil texture) is likely responsible for variable and complex changes in soil chemical and biological properties [44]. Spatial patterns of such soil properties interacting with water were likely highly confounded by tillage.


## 3.6.1. Soil Texture

Sand was a significant predictor of the coarse POM fractions, whereas clay was considered significant for controlling the fine POM fraction (Table6). Overall, sand content was positively related to the 1000–2000 μm POM fraction (Table 7). We speculate that this is because soils with high sand content could store more crop residues as coarse POM than soils with low sand content. Clay may play a role in sequestering partly decomposed POM through physical and chemical protection of soil C and N [1, 49, 55], probably causing different responses of the fine POM fraction to varying tillage intensity and texture. However, the fine POM fraction in the soil surface significantly decreased with clay content across the fields (Table 7). This could be explained by relatively low clay contents (<20%) of the fields that have a tendency of fast degradation of POM after tillage is introduced.Table 6
Significance of effects of soil physical properties on particulate organic matter (POM) C and N fractions, accounting for effects of field and time. All POM fraction values were log-transformed.
Effect>1000 μm POM250–1000μm POM53–250μm POMPOM-CP-valueField0.0001.00010.0001Time0.0001.00010.0001Sand0.0416.51010.0801Clay0.1504.35040.0010Bulk density0.0010.00010.0001Field× time0.0001.00010.0001Time× clay0.0043.01550.7712Time× bulk density0.0196.00160.0399Std. Dev. locationNA†.09690.0935Std. Dev. residual0.5762.27170.1499POM-NP-valueField0.0001.00010.0001Time0.0001.00010.0001Sand0.0001.0150.8818Clay0.1786.09120.0004Bulk density0.1767.00110.0001Field× time0.0001.00010.0001Field× clay–‡–0.2774Time× clay––0.7443Field×time×clay––0.0017Std. Dev. location0.1858.09320.0814Std. Dev. residual0.5061.27490.1631NA†: not applicable. The likelihood ratio test indicated that the random effect for location was not necessary and thus it was removed from the model according to Pinheiro and Bates (2000).The‡ effect was removed from the model because it was not significant at the level of 2%.Table 7
Standardized partial regression coefficients relating soil physical properties to particulate organic matter (POM)-C and -N fractions. Numbers in the columns labeled 98%CI are the half-width of the confidence intervals.
FractionTimeSand98%CIClay98%CIBulk density98%CIPOM-C>1000μmAugust 030.030.339-0.100.221April 04-0.230.355-0.480.504October 040.370.343-0.530.497June 050.130.354-0.690.473All0.270.237250–1000μmAugust 03-0.240.321-0.190.196April 04-0.450.334-0.870.447October 040.080.324-0.360.441June 05-0.160.334-0.760.419All-0.010.24753–250μmAugust 03-0.480.351-0.240.188April 04-0.620.362-0.660.431October 04-0.390.353-0.200.425June 05-0.550.362-0.640.402All-0.470.297POM-N>1000μmAll0.220.246-0.160.244-0.120.504250–1000μmAll0.020.23-0.210.228-0.250.27253–250μmAll-0.210.261-0.300.263-0.160.150For the coarse POM-C fractions, studies suggest that the degree of POM dependency on soil texture seems to be regulated by the level to which soil texture varies spatially between the fields with ST, NT, and MT [25, 56]. In addition, the relationships of sand and clay contents with other biotic and abiotic factors were largely affected by time and tillage [43]. As a result, this would significantly confound the dynamics of coarse POM. Our results are consistent with Wander et al. [40], showing that the relationship between soil texture and POM changed inconsistently over time. However, it was not evident what led to this shift in their relationship. In contrast, the relationship between the 53–250 μm POM-C fraction and clay appeared to be temporally consistent at the field scale. Therefore, the effect of clay on the dynamics of fine POM would be predominantly confounded by tillage and presumably independent of spatiotemporal variability in the other factors. Similarly, the effect of sand was positive or zero for the 250–2000 μm POM-N fractions and the effect of clay was consistently negative for the 53–250 μm POM-N fraction.


## 3.6.2. Soil Bulk Density

The BD was a significant predictor for all POM fractions except for the POM-N> 1000 μm (Table 6). Regardless of the tillage treatments, the BD had a consistently negative effect on POM-C in all fractions at the field scale. Although the effect of BD on POM was greater but more variable than the effects of sand and clay in the majority of the cases, no significant differences were evident (Table 7). Tillage-induced changes in BD could be in part related to POM dynamics. For example, da Silva et al. [39] and Wander et al. [40] showed that there was a trend for decreasing bulk SOM with increasing BD under cultivation. Ball et al. [57] showed that maximum BD was negatively related to water content, total porosity, and clay content. We speculate that the rate of residue decomposition and new POM formation would be limited under lower soil water availability due to a high BD, leading to decreased POM stocks. However, we were not able to find the evidence to support this speculation. The effects of BD on the POM-N fractions were similar to those of the POM-C, but slopes tended to be more constant over time.


## 3.6.3. Soil Water

The effect of soil water on POM was assessed by adding this factor to the reduced models obtained by backward elimination. Soil water was not significant for any of the response variables. This may be because water content is a too dynamic variable in comparison to POM. Although SOM mineralization can be greatly enhanced with increasing water content [45, 46], the spatial variation in soil conditions (e.g., soil texture) is likely responsible for variable and complex changes in soil chemical and biological properties [44]. Spatial patterns of such soil properties interacting with water were likely highly confounded by tillage.


## 4. Conclusion

Field level changes in soil POM was confounded by large spatial variability within and across the ST, NT, and MT fields, due in part to the influence of the natural field heterogeneity. Nevertheless, ST significantly accelerated changes in soil surface POM-C and -N in the initial period of the tillage conversion, but little seasonal changes occurred thereafter. Although coarse POM was more sensitive to tillage than fine POM, tillage-induced differences between the fields were greater for fine POM than for coarse POM. Accordingly, the relative contribution of total POM-C and -N to total SOM decreased in ST compared to NT. This suggests that fine POM is an important component of soil C and N cycling, at least in the short term. The faster turnover of new C in coarse POM under ST than NT indicates that reduced intensity tillage can improve SOM status by stabilizing pre-existing labile C and N sources in coarse POM.The short-term (<3 years) changes in POM-C, -N, and new C input affected by tillage were small relative to the spatiotemporal variability of non-POM or bulk SOM within each field. We found that the POM-C and -N themselves co-varied in space as a function of tillage at the field scale. Sand was a significant predictor of the short-term dynamics of coarse POM, whereas clay was considered significant for controlling the fine POM fraction. The BD was a significant predictor for all POM fractions except for the POM-N  > 1000 μm. However, the effect of BD on POM was usually temporally variable and independent of tillage. Although tillage affected BD as much as POM and its stocks were significantly related to BD, little is known about the mechanisms relating BD and soil POM-C and -N pools. Moreover, the control of clay or BD on coarse POM tended to vary greatly with time, probably due to the confounding effects of tillage and temporally variable factors, such as soil water. However, the relationship between fine POM and clay or BD was irrespective of field or sampling time, suggesting dominant tillage effects that are not significantly influenced by the spatiotemporal variation of other biophysical factors. In conclusion, in a typical California farming system, reduced tillage may not necessarily lead to soil C increase compared to ST but maintain the level of soil C. In addition, clay content and BD were identified as key factors influencing the spatiotemporal variability of POM.



---

*Source: 219379-2009-11-08.xml*