Abstract: High yields of rice (Oryza sativa L.) and wheat (Triticum aestivum L.) achieved in India following the green revolution have tended to decline after 5 to 8 years of continuous production in the same field. We hypothesized that degradation of soil physical quality was responsible for the yield decline because nutrient, pesticide and irrigation management were relatively well optimized. A soil physical quality assessment model was thus developed using 4 years of data from tillage and crop residue experiments conducted on a silty clay loam (Aquic hapludoll) at the Crop Research Centre, Pantnagar, India. Tillage for rice increased Bulk Density (BD) and Puddling Index (PI) and decreased saturated hydraulic conductivity (Ks), porosity (f) and Infiltration Rate (IR) over the 4 years. Changes in soil properties were significantly higher in the plots puddled by 4 Passes of Rotavator (PR) and in Conventional Puddling (CP) plots than in reduced puddling (ReP) and Direct Seeding Without Puddling (DSWP) plots. The crop residue management and wheat tillage treatments did not significantly influence soil physical properties in the 4 years time but the negative effect of puddling on soil properties was more in the Residue Removed (RR) plots than in the Residue Incorporated (RI) plots. Rice yield was highest in PR, which was at par with that in CP and ReP plots and wheat yield was highest in DSWP, which was at par with that in ReP, irrespective of tillage levels for wheat. Decrease in rice and wheat yields over the 4 years was not significant. Regression of grain yield on soil physical properties showed significant effects of BD and Ks on rice yield and that of IR and f on wheat yield. These relationships were used to develop a Soil Physical Quality Index (SPQI) model. Values of SPQI for treatments with grain yields that were statistically equal to maximum yield ranged from 0.70 to 0.80 for rice soils and 0.75 to 1.0 for wheat soils. Overall, tillage and residue treatments did not significantly degrade soil physical quality during the first 4 years, but projections based on the SPQI suggest they could in 5 to 9 years depending on the practices used for soil puddling and residue management. In DSWP, SPQI for rice was significantly lower from the first year itself (which indicates need for tillage) but was in the optimum range for wheat even under zero tillage. The proposed SPQI model has been verified using the published data on rice-wheat system under different soil conditions across the country.
Introduction
Mission to produce more food, feed and fiber encouraged rapid adoption of new technologies in agriculture without in-depth evaluation of their long-term impact on soil productivity and environmental quality (Doran et al., 1996). As a result of this, high yields of rice (Oryza sativa L.) and wheat (Triticum aestivum L.) observed in India in the initial years have tended to decline in about 5 to 8 years of continuous rice-wheat system (Tripathi, 1992a; Hobbs, 2001), questioning the sustainability of soil quality condition with the existing package of management practices. Whereas nutrient, pesticide and irrigation schedules have been relatively optimized for a cropping system, tillage and residue management, which are highly variable in different soils and environmental conditions, need thorough optimization (Karlen et al., 1997; Sharma et al., 2004; Tripathi et al., 2005). Karlen et al. (1997) observed that no-tillage or crop residue application improved soil quality compared with the practice of intensive tillage or crop residue removal. No-till system improved earthworm population, organic carbon sequestration and water stable aggregates but increased bulk density compared to the conventional tillage or reduced tillage (Arshad et al., 1999; Tebrágge and Dáring, 1999). However, Lal (1999) did not find any significant difference in bulk density between no-tillage and minimum tillage or plough tillage in 25 years of no-till system with continuous maize (Zea mays L.).
A higher amount of C and N observed in the 0 to 0.3 m layer of the rangeland soils subjected to 11 years of grazing than in the soil prohibited for grazing showed the evidence of improved soil quality due to greater opportunities for plant residue incorporation and less loss through oxidation (Manley et al., 1995). Liebig et al. (2004) also reported the improved soil quality and agricultural sustainability by residue incorporation with reduced tillage management. Tillage affects organic C, erosion and many physical, chemical and biological properties and processes used as indicators of soil quality (Karlen, 2004).
Larsen and Pierce (1991) defined soil quality as the capacity to function within the ecosystem boundaries and to interact positively with the environment external to that ecosystem. In agricultural context, the soil quality is aimed at maximizing production without adversely affecting the environment. Potential physical, chemical and biological indicators with minimum data set have been suggested to evaluate soil quality (Doran et al., 1996; Karlen et al., 2001). Efforts have been made to identify critical soil functions, select indicators for those functions and develop appropriate scoring functions to interpret the indicators for various soil resources (Karlen et al., 1994). Wander and Bollero (1999) showed that particulate organic matter was a promising indicator of soil quality. Masciandaro and Ceccanti (1999) observed that humic pool, which is often preserved even where most of the labile soil organic matter (>86%) is lost by tillage, was a potential indicator of soil quality. Similarly methane oxidation rate, which is enhanced with structural damage, could also be a potential soil quality indicator (Ball et al., 1999). Lapen et al. (2004) observed that air-filled porosity based soil quality indicator, such as least limiting water range, was best for maize throughout the growing season. On slopy lands, erosion risk was the reliable indicator of soil quality for sustainable management (Anken et al., 2004). Soil quality test kits (Liebig et al., 1996; Sarrantonio et al., 1996), farmer-based scorecards (Romig et al., 1996), soil management assessment framework (Karlen et al., 1997), soil resource management programs (Walter et al., 1997) and visual soil assessment procedures (Shepherd, 2000) have been developed and used to determine soil quality. Using standardized scoring functions approach, Hussain et al. (1999) reported that soil quality index is sensitive to local practices and responsive to inherent soil conditions and therefore, useful for making management recommendations. However, a more systematic approach is needed to objectively select the indicators for determining soil quality (Wang et al., 2003). There is no single correct indicator or index value because of inherent differences among soils, crops and climates (Andrew et al., 2002; Karlen et al., 2001; Sojka et al., 2003). Characterization of soil quality in terms of economic yield and changes in soil properties caused by tillage and other management practices is promising (Acton and Gregorich, 1995; Doran et al., 1996) but a limited information concerning tillage induced changes on soil properties are monitored (Gantzer and Blake, 1978; Hill, 1990; Hussain et al., 1999; Tripathi et al., 2003). In this study a method has been described to evaluate soil quality for rice and wheat from soil physical properties and grain yield as affected by tillage and residue management. The chemical and biological quality parameters were not examined due to their dependence on a number of interacting parameters and lack of expertise to separate them.
Materials and Methods
Experiments were conducted to generate data on tillage and crop residue effects on soil physical properties and yield of rice and wheat during 1999 to 2003 at the Crop Research Centre, Pantnagar (28°26N, 77°29E, 243.8 m above msl), India, in a silty clay loam of subgroup Aquic hapludoll (Deshpande et al., 1971). The area lies as east-west belt adjoining the Himalayan mountains and is characterized by non-saline natural water tables and sub-humid tropical conditions (Tripathi, 1997). The 4-year average rainfall during the rice and wheat seasons was 1715 and 638 mm, respectively. The water table fluctuated between 0.05 to 0.4 m during the rice season and 0.5 to 0.97 m depths during the wheat season.
Experimental Details
Rice experiment was laid out in split plot design in 12x10 m plots in 4
replications. Main plots consisted of 4 tillage treatments viz., PR (puddling
by 4 passes of rotavator), CP (conventional puddling by 4 passes of tractor
drawn local puddler with cage wheel), ReP (reduced puddling by 2 passes of rotavator)
and DSWP (direct seeding without puddling). The subplot treatments were Residue
Incorporation (RI) and Residue Removal (RR) in 6x10 m plot sizes. Tillage treatments
for wheat were Zero Tillage (ZT) and Conventional Tillage (CT) with RI and RR
superimposed over the plots of rice tillage treatments in split-split plot design
with rice tillage as main plots and wheat tillage as subplots and residue management
as sub-subplots. In ZT, wheat was planted by seed drill without any land preparatory
tillage. The CT for wheat included ploughing by one-way disc plough followed
by 3 harrowing, leveling, planking and planting by seed drill. The plot size
for wheat became 6x5 m. The rice crop was fertilized with 150 kg N ha-1,
26 kg P ha-1, 33 kg K ha-1 and 9 kg Zn ha-1
and the wheat crop was fertilized with 120 kg N ha-1, 26 kg P ha-1
and 33 kg K ha-1. Half of the nitrogen and full P, K and Zn were
applied as basal dose and the remaining dose of nitrogen was applied in two
splits. The rice crop was irrigated with 75 mm of water 3 days after the disappearance
of ponded water from the surface of the PR plot (Tripathi et al., 1986).
The wheat crop was irrigated based on climate-crop-soil data (Tripathi, 1992b).
The physical properties measured were Puddling Index (PI), Bulk Density (BD), saturated hydraulic conductivity (Ks), Infiltration Rate (IR) and porosity (f). The bulk density was measured by core method and the saturated hydraulic conductivity by constant head permeameter method at 30 Days After Sowing or Transplanting (DAS/DAT) and at harvesting of rice and wheat using soil cores collected from 0 to 0.05 m tilled layer. The infiltration rate was measured at 30 DAS/DAT and harvesting of the crops using double ring infiltrometer. The porosity was determined from BD and particle density data. The puddling index was determined as the volume ratio of soil settled in 24 h in a 100 mL measuring cylinder to that within 5 min after puddling. Harvesting the central 2x2 m area, threshing and air-drying the grain to constant moisture content was grain yield. The 4-year data were analyzed in a split-split plot design taking year as main factor and evaluating significance (p=0.05) by F-test (Gomez and Gomez, 1984).
Soil Physical Quality Index Model
A model to calculate Soil Physical Quality Index (SPQI) was developed from
the data on soil physical properties and yields of rice and wheat collected
during 1999 to 2003. The model utilizes multiple linear regressions of crop
yields on soil physical properties to test the relative significance of soil
properties on crop yield as
(1) |
where Y is grain yield of the crop; X1, X2 , Xn are different soil properties; and a, b1, b2, ., bn are constants. Those physical properties whose coefficients (b1, b2, bn) in Eq. (1) were found significant by t-test were selected for calculating the SPQI. The selected physical properties and yields were then normalized as X'j by the relation
(2) |
Maximum value of soil physical properties used and experimental maximum yields of rice and wheat are given in Table 1. The X'j was 1 for the values of Xio exceeding Xim. The normalization was done to extend the applicability of the model under diverse soils and climates. Since maximum yield is variety, soil and area specific, use of location specific maximum values of the parameters would be more appropriate. The normalized values of the selected physical properties were individually subjected to regression with normalized yields to determine their coefficient of determination, R2, whose proportionate variation was obtained and expressed as Ai.
(3) |
where n is number of properties Xi found significant by t-test. The soil physical quality index, SPQI, was then expressed as sum of the products of Ai and X'j as:
(4) |
Table 1: | Maximum value of soil physical properties and rice-wheat yields used in the calculation of SPQI |
Thereafter, SPQI was calculated for the normalized values of the physical properties and plotted against the corresponding normalized yields to determine SPQI at the observed maximum yield. A second-degree polynomial fitted the data. The value of SPQI at the observed maximum yield of the crop and the range of SPQI corresponding to yields statistically similar to the maximum yield represented the optimum range of SPQI. The SPQI of different treatments was then plotted against years to determine soil physical quality projections to follow the temporal trends under different management systems/treatments.
Results and Discussion
Dynamics of Soil Physical Properties
Soil physical properties are responsive to tillage and crop residue management
in rice-wheat system (Fig. 1 to 5). Tillage for rice (PR,
CP and ReP) caused a significant increase in PI and BD and decrease in Ks, f
and IR but tillage for wheat (ZT and CT) and residue management showed no significant
effect. Changes in soil physical properties were significantly higher in PR
and CP plots than in ReP and DSWP but changes in a treatment with time, due
to tillage for rice or wheat and residue management, were not significant in
the four years. However, increase in PI from 0.71 to 0.86, 0.69 to 0.81 and
0.57 to 0.71 in PR, CP and ReP, respectively, in 4 years, amply indicates increasing
breakdown of aggregates with time depending upon tillage intensity. Residue
Incorporation (RI) resisted the aggregate breakdown such that PI and BD in RR
plots were always higher than in RI (Fig. 1 and 2).
Increased water stable aggregates by residue incorporation and reduced tillage
have also been reported by Bajpai et al. (1995), Tebrágge and
Dáring (1999) and Kharub et al. (2004). Puddling effects on soil
physical properties persisted even in the wheat season and the effects were
more in RR than in RI plots and in ZT than in CT plots. Overall changes in physical
properties were rapid in the first 3 years.
Fig. 1: | Puddling index (PI) over the years in 2 (ReP-RI), 3 (CP-RI), 4 (PR-RI), 6 (ReP-RR), 7 (CP-RR), 8 (PR-RR) plots of rice |
Fig. 2: | Bulk density (BD) of the tilled layer over the years in different tillage and residue treatments of rice (A) and in zero tillage (ZT) and conventional tillage (CT) conditions of wheat after rice in 1 (DSWP-RI), 2 (ReP-RI), 3 (CP-RI), 4 (PR-RI), 5 (DSWP-RR), 6 (ReP-RR), 7 (CP-RR), 8 (PR-RR) plots |
Dynamics of Grain Yield
Rice yield in PR plot was highest and statistically at par with that in
CP and ReP (Fig. 6).
Fig. 3: | Hydraulic conductivity (Ks) of the tilled layer over the years in different tillage and residue treatments of rice (A) and in zero tillage (ZT) and conventional tillage (CT) conditions of wheat after rice in 1 (DSWP-RI), 2 (ReP-RI), 3 (CP-RI), 4 (PR-RI), 5 (DSWP-RR), 6 (ReP-RR), 7 (CP-RR), 8 (PR-RR) plots |
Fig. 4: | Porosity (f) of the tilled layer over the years in different tillage and residue treatments of rice (A) and in zero tillage (ZT) and conventional tillage (CT) conditions of wheat after rice in 1 (DSWP-RI), 2 (ReP-RI), 3 (CP-RI), 4 (PR-RI), 5 (DSWP-RR), 6 (ReP-RR), 7 (CP-RR), 8 (PR-RR) plots |
But wheat yield was lowest in the PR plot of rice and highest in the DSWP plot in which it was at par with that in ReP. In DSWP, rice yield was significantly lower than in other treatments from the first year itself. Sharma et al. (2005) also observed that tillage necessary to maximize rice yield deteriorates soil physical condition for wheat. Total average yield of rice and wheat was highest from RI plots of ReP-CT treatment combination.
Fig. 5: | Infiltration rate (IR) over the years in different tillage and residue treatments of rice (A) and in zero tillage (ZT) and conventional tillage (CT) conditions of wheat after rice in 1 (DSWP-RI), 2 (ReP-RI), 3 (CP-RI), 4 (PR-RI), 5 (DSWP-RR), 6 (ReP-RR), 7 (CP-RR), 8 (PR-RR) plots |
Both rice and wheat yields followed a declining trend over the years. The decrease in yield was more in RR than in RI plots. Also, wheat yield decline in CT was greatly slowed down in the fourth year indicating its superiority over continuous ZT.
Fig. 6: | Grain yield over the years in different tillage and residue treatments of rice (A) and in zero tillage (ZT) and conventional tillage (CT) conditions of wheat after rice in 1 (DSWP-RI), 2 (ReP-RI), 3 (CP-RI), 4 (PR-RI), 5 (DSWP-RR), 6 (ReP-RR), 7 (CP-RR), 8 (PR-RR) plots |
Tillage and crop residue effects on soil physical properties and yields suggest that wheat crop liked improvement or preservation of inherent condition of soil aggregates but rice liked puddling that caused breakdown of aggregates into soft mud and soil compaction. Crop residues improved rice-wheat yields and resisted deterioration of soil aggregates. Relative importance of soil physical properties on grain yield as influenced by tillage and residue management was quantified by multiple linear regressions (Eq. 1) as:
(5) |
(Calculated t-value for BD = 6.38, Ks = 2.23, IR = 0.40, f = 0.16 and PI = 0.29 and Table value = 2.03; number of observations = 96) and
(6) |
(Calculated t-value for BD = 0.58, Ks = 0.47, IR = 3.78 and f = 4.54; and Table value = 2.03; number of observations = 128).
Analysis showed that combined effects of the properties considered accounted for at least 60% variations in rice yield and 62% variations in wheat yield. The t-test indicated a significant effect of BD and Ks on rice yield (Eq. 5) and that of IR and f on wheat yield (Eq. 6).
Soil Physical Quality Index Model
Equation (5) and (6) formed the basis
of the Soil Physical Quality Index (SPQI) model. As described in the methodology,
proportionate variation of R2, obtained from regression of the normalized
significant physical properties (Eq. 2) on normalized yields
(Fig. 7), was determined (Eq. 3) and expressed
as in Eq. (4) to obtain where BD', Ks', IR' and f ' are their
normalized values calculated using Eq. (2).
Fig. 7: | Relationships between normalized values of (A) bulk density and rice yield, (B) hydraulic conductivity and rice yield, (c) porosity and wheat yield and (D) infiltration rate and wheat yield |
Fig. 8: | Normalized yield of (A) rice and (B) wheat as affected by Soil Physical Quality Index (SPQI) |
(7) |
(8) |
Thereafter, SPQI was calculated using Eq. (7 and 8) and plotted against the corresponding normalized yields (Fig. 8) to obtain the optimum range of SPQI as described in the methodology.
As shown in Fig. 8, optimum range of SPQI for optimum yields was from 0.70 to 0.80 for rice soils and from 0.75 to 1.0 for wheat soils. A lower SPQI for rice soils corresponds to higher level of puddling and yield, whereas a lower SPQI for wheat soils corresponds to higher yield and minimum tillage in rice-wheat system. The values of SPQI and rice yield indicate that soil physical quality remained optimum during the 4 years of rice-wheat system in all the puddled plots (Fig. 9). Extrapolation of the curves (equations) in Fig. 9, indicated that soil physical quality for rice would not deteriorate in 8, 6 and 5 years in ReP-RI, CP-RI and PR-RI plots and in 6, 5 and 5 years in ReP-RR, CP-RR and PR-RR plots, respectively. The SPQI for rice in DSWP plot remained in the suboptimal range from the first year itself. But for wheat, it was optimum in all the 4 years in DSWP and ReP plots of rice both under ZT and CT conditions and may remain optimum until 8 and 6 years in DSWP-RI and ReP-RI plots and until 7 and 5 years in DSWP-RR and ReP-RR plots under ZT conditions and until 9 and 6 years in DSWP-RI and ReP-RI plots and until 7 and 5 years in DSWP-RR and ReP-RR plots under CT conditions, respectively (Fig. 9).
Fig. 9: | Changes in Soil Physical Quality Index (SPQI) over the years in different tillage and residue treatments of rice (A) and in zero tillage (ZT) and conventional tillage (CT) conditions of wheat after rice in Y1 (DSWP-RI), Y2 (ReP-RI), Y3 (CP-RI), Y4 (PR-RI), Y5 (DSWP-RR), Y6 (ReP-RR), Y7 (CP-RR), Y8 (PR-RR) plots |
Results thus show that soil physical quality deterioration was (i) least (not significant) with optimum yields of rice in ReP-RI plot for 8 years and that of wheat in the same plot under CT-RI condition for 6 years, (ii) faster in highly puddled PR-RR and CP-RR plots of rice and in ZT-RR plot of wheat and (iii) resisted by residue incorporation.
Karlen et al. (1997) and Liebig et al. (2004) also reported improvements in soil quality and sustainability of cropping systems by residue incorporation and reduced tillage.
Verification of the Soil Physical Quality Index Model
Equations (7 and 8) were used to calculate
SPQI using the data from the experiments on rice-wheat system conducted at New
Delhi (Aggarwal et al., 1999), Pantnagar (Kumar and Tripathi, 1990; Bharadwaj
et al., 1994; Bajpai and Tripathi, 2000; Sharma et al., 2004;
Lal, 2005) and Kharagpur (Kar, 2002). As shown in Table 2,
the SPQI for optimum grain yield was in the range of 0.74 to 0.80 for rice soils
and 0.77 to 0.98 for wheat soils, which is well within the proposed optimum
range of SPQI calculated from the model for rice and wheat soils. Results thus
indicate that SPQI in the range of 0.70 to 0.80 for rice soils and 0.75 to 1.0
for wheat soils could be used as indicators of the optimum range of SPQI for
optimum productivity and quality of rice and wheat soils.
Table 2: | Soil physical quality index calculated from the published data of different locations in India |
*Optimum rang of SPQI for rice soil = 0.70 to 0.80 and for wheat soil = 0.75 to 1.0, **Significantly high yield |
Conclusions
A model to calculate SPQI from grain yield and soil physical properties of the tilled layer as affected by tillage and residue management has been developed. The model uses multiple linear regressions to identify significant soil physical properties affecting grain yield and to develop expressions for determining SPQI for soils under rice-wheat cropping system. The parameters have been normalized to enable soil physical quality assessment under diverse soil conditions. The model has been verified on variety of soils in different regions of India under rice-wheat system. Optimum SPQI ranged from 0.70 to 0.80 for rice soils and from 0.75 to 1.0 for wheat soils. The temporal trends in SPQI indicated that soil physical quality for rice-wheat system did not significantly degrade by tillage and residue treatments in the first 4 years but may degrade in 6 to 8 years even in the most optimum treatment combination of ReP-CT under RI in this silty clay loam. The SPQI model can thus be used to determine soil physical quality to keep track of sustainability of the cropping system and changes in soil physical properties resulting from management.
Acknowledgments
Authors are grateful to Indian Council of Agricultural Research and NATP for funding this ad hoc research project.