Print version ISSN 0103-9016
Sci. agric. (Piracicaba, Braz.) vol.63 no.6 Piracicaba Nov./Dec. 2006
SOILS AND PLANT NUTRITION
Critérios estatísticos na seleção de parâmetros para avaliação da resistência do solo à penetração
Tairone Paiva LeãoI; Álvaro Pires da SilvaII, *
IUniversity of Tennessee, Department of Earth and Planetary Sciences, 306 Earth and Planetary Sciences Building, TN - 37996-1410 - USA
IIUSP/ESALQ - Depto. de Ciência do Solo, C.P. 9 - 13418-900 - Piracicaba, SP-Brasil
Measurements of soil penetration resistance (SR) have been frequently used for the evaluation of soil structural quality for plant growth. However, different data analysis approaches have been used, without a previous evaluation of their statistical quality. In this study we tested the hypothesis that the mean is the parameter with best statistical properties to evaluate alterations in soil penetration resistance in response to soil use and management, as compared to other SR statistical parameters. Undisturbed (5 × 5 cm) soil cores were collected from three sampling sites with different degrees of compaction: an undisturbed site under native scierophylous forest (NC); a site under short-duration grazing with post-grazing residue maintained at 2.0 to 2.5 Mg (Total Dry Matter) TDM ha-1 (LR); and a site under short-duration grazing with post-grazing residue maintained at 3.0 to 3.5 Mg TDM ha-1 (HR). The statistical quality of the parameters from undisturbed soil samples of SR profiles: mean (x), median (M), maximum (max), percentage of linear penetrability at 2 MPa (PLP2MPa), and the parameters from Probit analysis intercept (n) and slope (m) was evaluated using the ANOVA and LSD tests. Results from the F ratio, P > F values and LSD tests show that mean, median and maximum were the parameters with better statistical properties as criteria to evaluate alterations in soil penetration resistance in response to soil use and management as compared to other statistical SR parameters, validating the hypothesis of the research.
Key words: LSD, probit analysis, means comparison test, soil physical quality assessment
A quantificação da resistência do solo à penetração (RP) é freqüentemente utilizada na avaliação da qualidade estrutural do solo para o crescimento de plantas. Diferentes abordagens têm sido utilizadas na análise de dados de RP, e na maioria dos casos, sem uma avaliação prévia da qualidade estatística dos parâmetros. Neste trabalho foi testada a hipótese de que a média é o parâmetro com melhores propriedades estatísticas como critério para avaliação de alterações na resistência à penetração do solo em resposta ao uso e manejo do solo, em relação a outros parâmetros estatísticos derivados de conjuntos de dados de RP. Amostras indeformadas de solo (5 × 5 cm) foram coletadas em locais com diferentes graus de compactação: Cerrado nativo (CN); pastejo rotacionado, com nível de resíduo pós-pastejo mantido entre 2,0 e 2,5 Mg (Matéria Seca Total) MST ha-1 (RB); e pastejo rotacionado, com nível de resíduo pós-pastejo mantido entre 3,0 e 3,5 Mg MST ha-1 (RA). A qualidade estatística dos parâmetros obtidos dos perfis de RP: média (x), mediana (M), máximo (max), percentagem linear de penetrabilidade à 2 MPa (PLP2MPa), e os parâmetros da análise Probit: intercepto (n) e inclinação (m) foi avaliada por meio de ANOVA e pelo teste da Diferença Mínima Significativa (DMS). Os resultados da razão F, P > F e do teste DMS indicaram que a média, mediana e máximo são os parâmetros com melhores propriedades estatísticas como critério para avaliar as alterações na resistência à penetração do solo em resposta ao uso e manejo, confirmando a hipótese testada.
Palavras-chave: DMS, análise Probit, teste de comparação de médias, avaliação da qualidade física
An ideal methodology to characterize the resistance of soil matrix (SR) to deformation by a growing root should take into account as many of the unique features of root development as possible. These would include: root size, a lubrication system, an ability to develop both axial and radial pressures, rate of root tip advancement, and mainly the ability to seek out pores of appropriate size and resistance to deformation (Groenevelt et al., 1984; Whiteley & Dexter, 1984; Bengough & MacKenzie, 1994). Statistical parameters such as the mean (i.e. arithmetic average) or the maximum SR value of a soil layer can only give a rough idea of the actual limitation of the soil matrix to root growth.
To minimize the influence of the soil structural heterogeneity on the soil structural quality evaluation, Groenevelt et al. (1984) defined the Percentage Linear Penetrability (PLP) as the percentage of linear trajectory through soil, for which the actual tip resistance is lower than a set critical value. The PLP concept was further developed by incorporating linear regression analysis on Probit (Finney, 1952; Tietjen, 1986) transformed PLP data (Perfect et al., 1990). The regression procedures resulted in two empirical parameters, the slope (m) and the intercept (n) of the lines. It has been suggested that the slope, m, reflects the dispersion in the aggregate penetration resistance, while the intercept, n, reflects the soil macroporosity (Perfect et al., 1990). The mean is the statistical parameter most commonly used to evaluate resistance to penetration in soils. In this paper we tested the hypothesis that the mean is the parameter with better statistical properties as criteria to evaluate alterations in soil penetration resistance in response to soil use and management as compared to the SR statistical parameters: median, maximum, m, n, and PLP at the critical value of 2 MPa (PLP2MPa). The objective of this research was to evaluate comparatively the statistical quality of the parameters from the penetration resistance profile: mean, median, maximum, m, n, and PLP at the critical value of 2 MPa (PLP2MPa).
MATERIAL AND METHODS
The study was conducted with soil cores collected in the Mato Grosso do Sul state, Brazil (20º26' 48'' S, 54º43'19'' W). According to Köppen's classification, the climate is defined as a transition between Cfa (subtropical humid) and Aw (Tropical wet-dry), with a mean annual precipitation of 1500 mm and mean temperature of 22ºC. The soil is classified as Typic Acrudox with 399 g kg-1 clay, 66 g kg-1 silt, and 535 g kg-1 sand.
Three sampling sites were selected for the study:
(i) Native Cerrado (NC): an undisturbed savanah site under native scierophylous forest called "Cerradão";
(ii) Short-duration Grazing with Lower Post-Grazing Residue (LR): the site (0.18 ha) cultivated with Tanzania grass (Panicum maximum cv. Tanzania) was grazed for seven days followed by resting periods of 35 days. This site is part of an intensive short-duration grazing system experiment established in 1999. The experiment was designed to maintain soil base saturation from 45 to 50%; available phosphorous levels in Mehlich-1 (Nelson et al., 1953) from 4 to 8 mg dm-3; and potassium from 60 to 80 mg dm-3, at the 0 to 0.2 m layer. The nitrogen fertilization rate was 100 kg N ha-1 yr-1. The post-graze residue at this site was maintained at 2.0 to 2.5 Mg (Total Dry Matter) TDM ha-1. Average annual stocking rates were 4.12 (Animal Units) AU ha-1 and 2.26 AU ha-1 for wet and dry seasons, respectively.
(iii) Short-duration Grazing with Higher Post-Grazing Residue (HR): Same features as the LR site, except that the post-graze residue at this site was maintained at 3.0 to 3.5 Mg TDM ha-1. Average annual stocking rates were 4.8 AU ha-1 and 2.26 AU ha-1 for wet and dry seasons, respectively.
Twenty-seven undisturbed (5 × 5 cm) soil cores were collected at each site. The cores were subdivided in three groups of nine samples. Each group was dried to three water content ranges, namely pre-field capacity range (bfc), field capacity range (fcr) and post-field capacity range (pfc). The bfc consisted of samples equilibrated at the tension of 4 KPa, fcr at 10 KPa and pfc at 0.5 MPa. After equilibration, the cores were stored in a refrigerator for one month and used in the evaluation of the soil penetration resistance (SR) profiles. The SR profile is the set of penetration resistance values collected along the vertical trajectory of the penetrometer in an undisturbed soil sample (Figure 1).
These profiles were collected using an electronic penetrometer device directly connected to a computer parallel port. The cone penetrometer used has a semi-angle of 30º and 3.86 mm of basal diameter. The penetration rate was 1.0 cm min-1, with values recorded every 0.7 s. Immediately after the SR tests, soil cores were oven dried at 105ºC for 24 hours and volumetric water content (qV) and bulk density (Db) were calculated.
For each SR profile the statistical parameters: mean (), median (M), maximum (max), and the parameters from Probit analysis slope (m), intercept (n) and percentage linear penetrability at the critical value of 2 MPa (PLP2MPa) (Groenevelt et al., 1984; Perfect et al., 1990) were calculated. For each water content range, the set of SR statistical parameters for the sampling sites (NC, HR, LR) was evaluated using the F ratio from ANOVA (Gravetter & Wallnau, 1995) and mean comparisons by Least Significant Difference test (LSD) (Hsu, 1996). All statistical analyses were performed using the SAS/STAT software package (SAS Institute, 1999).
RESULTS AND DISCUSSION
Mean values for volumetric water content (qv), dry bulk density (Db), and for the statistical parameters from the soil penetration resistance (SR) profiles: mean (), median (M), maximum (max), slope (m), intercept (n) and percentage linear penetrability at the critical value of 2 MPa (PLP2MPa) are presented in Table 1. Parameters are grouped by sampling site and water content range, therefore the statistical parameters from the SR profiles are evaluated not only according to their sensitivity to soil use and management, but the influence of soil water content over these relationships was also investigated.
The sites used in this study were selected so that three different and clearly distinct degrees of compaction could be used in SR parameter evaluations. Mean Db values were lower in the NC site, Db = 0.97 g cm-3, and increased with the stocking rate in the grazed sites, Db = 1.30 g cm-3 in HR and Db = 1.42 g cm-3 in LR. For clay Oxisols such the one evaluated here, these Db values are usually associated with loose soil (NC) (Normally undisturbed native vegetation), slightly compacted (HR), and highly compacted sites (LR), respectively. Once stated that the difference in compaction degree between the sites was previously known, the statistical quality of each SR profile parameter was initially evaluated by the magnitude of the F ratio from the ANOVA and by its statistical significance level (P > F) (Gravetter & Wallnau, 1995).
According to the classical statistical theory, the higher the F ratio and the lower the P > F value the better is the statistical quality of the indicator (Gravetter & Wallnau, 1995). However, it is worth noting that not only the statistical quality of the parameter was investigated but also its physical significance in the process was taken into account in this discussion (Webster, 2001). The F ratio, P > F, and r2 from the ANOVA used for comparison of SR profile parameters between sampling sites are presented in Table 2, Table 3, and Table 4, for the water content ranges: bfc, fcr and pfc respectively.
At the three evaluated water content ranges the higher F ratio and lower P > F value were found for the statistical parameters , max, M and PLP2MPa (Tables 2 to 4). The and max are the SR parameters most commonly used for evaluation of soil physical quality, or compaction degree (Corrêa & Reichardt, 1995; Sojka et al., 2001; Villamil et al., 2001). The results presented here only confirm that these parameters are adequate to evaluate the SR in soil quality studies. However, the PLP2MPa may be a more meaningful parameter since it indicates the percentage of the linear trajectory (PLP) through a soil sample (or layer) that would be penetrable with resistance values lower than 2 MPa (Groenevelt et al., 1984). Mean, maximum and median values give only a rough idea of the degree of mechanical resistance of a soil. For mean and maximum, besides the fact that they represent a whole soil layer with a single value, they can be highly affected by extreme observations. Despite not being affected by extreme observations, the median has not been frequently used as an indicator parameter of soil resistance, which is probably related to the fact that it cannot give a true representation of the range of resistances in a soil sample/profile nor give an indication of high SR values in layers/spots which can limit root growth under field conditions. The PLP, on the other hand, estimates the fraction of the linear trajectory through a soil that is readily penetrable by a growing root, given a certain critical penetration resistance value. Therefore, the PLP could give a better picture of the degree of penetrability of a soil layer than any other indexes (Groenevelt et al., 1984).
The critical SR value of 2 MPa for the PLP was chosen since there is experimental evidence that root growth is severelly restricted at this condition (Taylor et al., 1966; Materechera et al., 1991). Thus at the PLP2MPa a given percentage of the soil profile would be penetrable at values lower then 2 MPa. For any water content range, the PLP2MPa was greater at the NC site (ranging from 88.74 to 98.78%) and much lower at the LR site (ranging from 0 to 26.5%). However it is not known which percentage of the soil profile would need to be penetrable under a certain mechanical resistance for adequate root growth to occur. What is known is that plant roots are able to "search" for rooting paths in pores, fissures and zones of lower mechanical resistance through soil (Groenevelt et al., 1984; Whiteley & Dexter, 1984; Dexter, 1986a; 1986b). It has been suggested that the parameters m (slope) and n (intercept) from the Probit regression analysis could represent this soil structural heterogeneity (Perfect et al., 1990). The m parameter would reflect the spread in distribution of aggregate strengths while n would reflect the macroporosity (Perfect et al., 1990). However, their statistical quality as criteria for evaluation of SR differences among different soil management strategies was lower than the other evaluated parameters (Tables 2 to 4).
To further investigate the ability of the parameters as criteria to categorize SR among sampling sites with different compaction degrees, mean comparison tests were performed for the SR parameters with higher F ratio and lower P < F value. Tables 5, 6, 7 and 8 present the LSD means comparison test for PLP2MPa, , max and M soil resistance parameters, respectively, between sampling sites. Theoretically, the more sensitive the parameter to SR alterations by soil management, the more easily the LSD test could identify statistically significant differences among sites. The three parameters from the SR profile, PLP2MPa, , max and M identified the same trend of soil physical quality decrease at any water content range (Tables 5 to 8). However, the interpretation of PLP2MPa results is different, since the sampling sites with best soil physical quality will have higher values (Table 5). Therefore, at any water content range, the PLP2MPa value was lower at LR, but they did not differ statistically between the NC and HR sites for the bfc and fcr water content ranges. At pfc the PLP2MPa did not differ statistically between LR and HR sites. This behavior can be explained by the fact that in compacted sites the SR increases sharply with the decrease in water content while the magnitude of this increase is lower in loose soil (Imhoff et al., 2000).
The same trends were observed between sampling sites for the , max and M. At any water content range the , max and M values were greater at the LR site. As one might expect, the LR site was associated with the worst soil structural properties, a trend which was also identified by other physical parameters, like Db (Table 1), range of critical bulk density values and least limiting water range (Leão et al., 2004). At higher water contents (bfc and fcr) the LSD test for the max and PLP2MPa parameters identified the same statistical trend (Tables 5 and 7). The SR properties were better and not statistically different for the NC and HR sites. On the other hand, even at with higher water content, the LR site soil had low structural quality as evaluated by any SR profile parameter (Tables 5 to 8). At any water content range, M and were the more sensitive parameters to assess changes in soil resistance characteristics between sites (Tables 6 and 8).
With drier soil (pfc) the sensitivity of the LSD test was equal for the , max, and M parameters (Tables 6, 7 and 8). At this water content range (pfc) the LSD test for the , max and M parameters identified the three sites as statistically different, the value of the parameter decreasing from LR to NC site. The SR characteristics in a soil are mainly influenced by bulk density and water content and by their interactions (Sojka et al., 2001). However, as one might observe in surface plots of SR as a function of both bulk density and water content (Imhoff et al., 2000; Sojka et al., 2001) the influence of bulk density over SR characteristics increases as the soil water content decreases. The decrease in SR with the increase in qv is due to the reduction in cohesion and in the angle of internal friction as soil water content increases (Camp & Gill, 1969). The increase in SR with Db is attributed to the effect of compaction over soil matrix, and the increase in interparticle friction as soil particles become closer with soil compaction (Vepraskas, 1984; Sojka et al., 2001). At lower water contents, , max and M were more reliable in identifying these relationships than any other SR profile parameter.
To test if the PLP2MPa could be estimated from easily measured soil physical properties (i.e. water content, bulk density, and resistance parameters) for the specific conditions of the soil under evaluation, a simple stepwise multiple regression exercise was carried out. The PLP2MPa was inputted in the stepwise regression model as a function of bulk density, volumetric water content and the soil resistance parameters: maximum, minimum, mean and median. The PLP2MPa for the soil under evaluation can be estimated with a high degree of accuracy from bulk density and mean soil resistance values (i.e. arithmetic average).
where: PLP2MPa = Percentage of linear penetrability at the soil resistance of 2MPa (%); Db = Soil bulk density (g cm-3); SR = Mean soil resistance (MPa); Db > 0; 0 < PLP2MPa < 100; SR > 0.
The mean, maximum and median are the parameters with better statistical properties as criteria to evaluate alterations in soil penetration resistance in response to soil use and management as compared to the statistical SR parameters: m, n, and PLP at the critical value of 2 MPa (PLP2MPa). The traditional indicator of SR, the mean, is adequate to evaluate the differences in SR among sites with different use and management. However, it would be more effective at lower water contents. A new parameter PLP2MPa is also proposed, but its sensitivity to SR changes should be further investigated.
BENGOUGH, A.G.; MACKENZIE, C.J. Simultaneous measurements of root force and elongation for seedling pea roots. Journal of Experimental Botany, v.45, p.95-102, 1994. [ Links ]
CAMP, C.R.; GILL, W.R. The effect of drying on soil strength parameters. Soil Science Society of America Proceedings, v.33, p.641-644, 1969. [ Links ]
CORREA, J.C.; REICHARDT, K. Efeito do tempo de uso das pastagens sobre as propriedades de um Latossolo Amarelo da Amazônia Central. Pesquisa Agropecuária Brasileira, v.30, p.107-114, 1995. [ Links ]
DEXTER, A.R. Model experiments on the behaviour of roots at the interface between tilled seed-bed and compacted sub-soil. II. Entry of pea and wheat roots into sub-soil cracks. Plant and Soil, v.95, p.135-147, 1986a. [ Links ]
DEXTER, A.R. Model experiments on the behavior of roots at the interface between tilled seed-bed and compacted sub-soil. III. Entry of pea and wheat roots into cylindrical biopores. Plant and Soil, v.95, p.149-161, 1986b. [ Links ]
FINNEY, D.J. Probit analysis: a statistical treatment of the sigmoid response curve. Cambridge: University Press, 1952. 318p. [ Links ]
GRAVETTER, F.J.; WALLNAU, L.B. Essential statistics for the behavioral sciences. 2.ed. Minneapolis: West Publishing Company, 1995. 431p. [ Links ]
GROENEVELT, P.H.; KAY, B.D.; GRANT, C.D. Physical assessment of a soil with respect to rooting potential. Geoderma, v.34, p.101-114, 1984. [ Links ]
HSU, J.C. Multiple comparisons: theory and methods. London: Chapman & Hall, 1996. 277p. [ Links ]
IMHOFF, S.; SILVA, A.P.; TORMENA, C.A. Aplicações da curva de resistência no controle da qualidade física de um solo sob pastagem. Pesquisa Agropecuária Brasileira, v.35, p.1493-1500, 2000. [ Links ]
LEÃO, T.P.; SILVA, A.P.; MACEDO, M.C.M.; IMHOFF, S.; EUCLIDES, V.P.B. Intervalo hídrico ótimo na avaliação de sistemas de pastejo contínuo e rotacionado. Revista Brasileira de Ciencia do Solo, v.28, p.1-8, 2004. [ Links ]
MATERECHERA, S.A.; DEXTER, A.R.; ALSTON, A.M. Penetration of very strong soils by seeding of different plant species. Plant and Soil, v.135, p.31-41, 1991. [ Links ]
NELSON, W.L.; MEHLICH, A.; WINTERS, E. The development evaluation and use of soil tests for phosphorus availability. In: PIERRE, W.H., NORMAN, A.G. (Ed.). Soil fertilizer phosphorus in crop nutrition. New York: Academic Press, 1953. p.153-188. [ Links ]
PERFECT, E.; GROENEVELT, P.H.; KAY, B.D.; GRANT, C.D. Spatial variability of soil penetrometer measurements at the mesoscopic scale. Soil and Tillage Research, v.16, p.257-271, 1990. [ Links ]
SAS INSTITUTE. SAS/STAT user's guide: version 8. Cary, 1999. [ Links ]
SOJKA, R.E.; BUSSCHER, W.J.; LEHRSCH, G.A. In situ strength, bulk density, and water content of a durinodic xeric haplocalcid soil. Soil and Tillage Research, v.166, p.520-529, 2001. [ Links ]
TAYLOR, H.M.; ROBERSON, G.M.; PARKER, J.J. Soil strength-root penetration relations to coarse textured materials. Soil Science, v.102, p.18-22, 1966. [ Links ]
TIETJEN, G.L. A topical dictionary of statistics. New York: Chapman and Hall, 1986. 171p. [ Links ]
VEPRASKAS, M.J. Cone index of loamy sands as influenced by pore size distribution and effective stress. Soil Science Society of America Journal, v.48, p.1220-1225, 1984. [ Links ]
VILLAMIL, M.B.; AMIOTTI, N.M.; PEINEMANN, N. Soil degradation related to overgrazing in the semi-arid southern caldenal of Argentina. Soil Science, v.166, p.441-452, 2001. [ Links ]
WEBSTER, R. Statistics to support soil research and their presentation. European Journal of Soil Science, v.52, p.331-340, 2001. [ Links ]
WHITELEY, G.M.; DEXTER, A.R. The behaviour of roots encountering cracks in soil. II. Development of a predictive model. Plant and Soil, v.77, p.141-149, 1984. [ Links ]
Received May 17, 2006
Accepted October 06, 2006