Efeito de propriedades intrínsecas do solo em avaliações de qualidade do solo

rbcs Revista Brasileira de Ciência do Solo Rev. Bras. Ciênc. Solo 1806-9657 Sociedade Brasileira de Ciência do Solo Viçosa, MG, Brazil A avaliação da qualidade do solo depende do uso de indicadores e índices construídos com as propriedades dele. Entretanto, essas propriedades podem variar entre tipos de uso do solo, bem como interferir nas propriedades do solo usadas para construir os índices. Objetivou-se com este estudo determinar a capacidade de propriedades intrínsecas do solo (teor de argila e óxido de ferro) para explicar as propriedades variáveis (soma de bases, acidez potencial, carbono orgânico, porosidade e densidade) em diferentes usos do solo (floresta nativa, plantio direto e agricultura convencional) em pequenas propriedades familiares no sul do Brasil. Os resultados evidenciaram que as cinco propriedades do solo estudadas podem ser incluídas em avaliações de qualidade e não são influenciadas pelo teor de argila e óxido de ferro. Sugeriu-se que, para solos pouco intemperizados, ricos em filossilicatos do tipo 1:1 e 2:1, se o teor de argila estiver em um intervalo de aproximadamente 200 g kg-1 e o de óxido de ferro for menor que aproximadamente 15 g kg-1, o uso do solo é o maior fator para controlar as propriedades físico-químicas da camada superficial do solo. DIVISÃO 2 - PROCESSOS E PROPRIEDADES DO SOLO COMISSÃO 2.4 - QUÍMICA DO SOLO The effect of intrinsic soil properties on soil quality assessments Efeito de propriedades intrínsecas do solo em avaliações de qualidade do solo Alessandro Samuel-RosaI; Ricardo Simão Diniz DalmolinII; Pablo MiguelIII; Jovani ZalamenaIV; Deborah Pinheiro DickV IDoctoral student, Pos-graduation in Agronomy - Soil Science, Rio de Janeiro Rural Federal University. BR 465, Km 7, Zona Rural. CEP 23890-000 Seropédica (RJ), Brazil. CNPq scholarship. E-mail: [e-mail] IIProfessor, Soil Science Departament, Santa Maria Federal University - UFSM. Av. Roraima, 1000, Camobi. CEP 97105-900 Santa Maria, RS. CNPq grantee. E-mail: [e-mail] IIIProfessor, College of Agronomy Eliseu Maciel, Pelotas Federal University. Campus Universitário, s/n. CEP 96010-900 Pelotas (RS), Brazil. E-mail: [e-mail] IVPos Doctoral student, Pos-graduation in Soil Science, UFSM. FAPERGS scholarship. E-mail: [e-mail] VProfessor, Chemistry Institut, Rio Grande do Sul Federal University. Av. Paulo Gama, 110, Bairro Farroupilha. CEP 90040-060 Porto Alegre (RS), Brazil. CNPq grantee. E-mail: [e-mail] ABSTRACT The assessment of soil quality is based on indicators and indices derived from soil properties. However, intrinsic soil properties may interfere with other soil properties that vary under different land uses and are used to calculate the indices. The aim of this study was to assess the extent to which intrinsic soil properties (clay and iron oxide contents) explain variable soil properties (sum of bases, potential acidity, organic carbon, total porosity, and bulk density) under different land uses (native forest, no-tillage and conventional agriculture) on small family farms in Southern Brazil. The results showed that the five properties evaluated can be included in soil quality assessments and are not influenced by the clay and iron oxide contents. It was concluded that for little weathered 1:1 and 2:1 phyllosilicate rich-soils, if the difference between the maximum and the minimum clay content under the different land uses is less than about 200 g kg-1 and the iron oxide content less than about 15 g kg-1, the physico-chemical soil properties in the surface layer are determined mostly by the land use. Index terms: soil quality, family farms, ANCOVA, spurious correlation. RESUMO A avaliação da qualidade do solo depende do uso de indicadores e índices construídos com as propriedades dele. Entretanto, essas propriedades podem variar entre tipos de uso do solo, bem como interferir nas propriedades do solo usadas para construir os índices. Objetivou-se com este estudo determinar a capacidade de propriedades intrínsecas do solo (teor de argila e óxido de ferro) para explicar as propriedades variáveis (soma de bases, acidez potencial, carbono orgânico, porosidade e densidade) em diferentes usos do solo (floresta nativa, plantio direto e agricultura convencional) em pequenas propriedades familiares no sul do Brasil. Os resultados evidenciaram que as cinco propriedades do solo estudadas podem ser incluídas em avaliações de qualidade e não são influenciadas pelo teor de argila e óxido de ferro. Sugeriu-se que, para solos pouco intemperizados, ricos em filossilicatos do tipo 1:1 e 2:1, se o teor de argila estiver em um intervalo de aproximadamente 200 g kg-1 e o de óxido de ferro for menor que aproximadamente 15 g kg-1, o uso do solo é o maior fator para controlar as propriedades físico-químicas da camada superficial do solo. Termos de indexação: qualidade do solo, agricultura familiar, ANCOVA, correlação espúria. INTRODUCTION The clearing of natural vegetation for agriculture leads to alterations in soil properties and, in several cases, to environmental degradation (Bewket & Stroonijdeb, 2003; Muniz et al., 2011). This phenomenon calls for a concept of soil quality, which is the capacity of a soil to function within the limits of an ecosystem, interacting positively with the environment external to that ecosystem (Larson & Pierce, 1991). A soil with good quality fulfils a number of functions (carbon storage, water reservoir, nutrient cycling, growth medium, among others), while simultaneously protecting other environmental components and sustaining human health (Doran & Parkin, 1994). In this sense, assessing the soil quality is one of the main tasks of soil scientists (Basher, 1997). However, in many regions there is still a lack of studies on soil quality (Reichert et al., 2003). One of these is the hillslope area of central Rio Grande do Sul (RS), Brazil, with mostly small family farms. These areas were first settled in the late XIX century by European immigrants, when native vegetation was replaced by agriculture. Today, farming is based on the integration of crop and livestock. There is little investment in machinery and fertilizers because of the complex topography, thin soils, small property sizes and low income of the farmers. On the other hand, higher-income farms with deeper soils and better road access adopted new and more intensive agricultural practices and activities, requiring investments in machinery and fertilizers. But soil conservation practices are insufficient, characterizing unsustainable agriculture and increasing the potential of land degradation (Samuel-Rosa et al., 2011). The soil quality assessment has to be based on the use of indicators and indices that reflect soil processes, integrate different soil properties and are sensitive to environmental changes (Doran & Parkin, 1994). Many soil properties can be included in such soil quality indices (Larson & Pierce, 1991; Doran & Parkin, 1994; Islam & Weil, 2000; Geissen & Guzman, 2006; Clemens et al., 2010). However, the heterogeneous character of soils must also be taken into consideration: intrinsic soil properties such as particle-size distribution and the type and content of oxides may vary among land uses. This variation can have a pronounced effect on the soil properties used to develop indices of soil quality (Miralles et al., 2009), hampering the evaluation of the true effect of land use. The aim of this study was to assess the extent to which intrinsic soil properties (clay and iron oxide contents) explain variable soil properties (sum of bases, potential acidity, organic carbon, total porosity, and bulk density) under different land uses (native forest, no-tillage and conventional agriculture) on small family farms in Southern Brazil. MATERIAL AND METHODS The data used in this study were measured on small family farms in the hillslope areas of central Rio Grande do Sul (RS), Brazil. The local topography is complex, with a relief ranging from undulated to mountainous. The altitude varies between 240 and 450 m asl. The climate is classified as Cfa (humid subtropical climate without a defined dry season), with a mean annual temperature of 19.2 oC and mean annual pluvial precipitation of 1,708 mm, quite well distributed throughout the year (Maluf, 2000). Geologically, the area consists predominantly of igneous rocks. The natural vegetation (Atlantic semi-deciduous forest) is dominated by Patagonula americana, Parapiptadenia rigida, Myrocarpus frondosus, Cordia trichotoma, Cedrela fissilis, Apuleia leiocarpa and Luehea divaricata. Three land uses were studied: - native forest: nearly all areas under native forest today had been used for agriculture in the past century. European immigrants cleared the original native forest in the late XIX century and practiced conventional agriculture for several decades. Since the 1970s, areas with shallow (<50 cm) and stony soils were gradually abandoned, allowing the regrowth of native vegetation; - no-tillage: >20 years of conventional agriculture - cultivation of potato, Solanum tuberosum - followed by 4 years of no-tillage - soybean, Glycine max, and ryegrass, Lolium multiflorum; - conventional agriculture: >20 years of conventional agriculture - mayze, Zea mays, and ryegrass. The soils are little weathered. Leptosols prevail in native forest and Acrisols in agricultural areas. The mineralogical assemblage of both soil types is similar. The soil clay mineral kaolinite (Kt), common in Southern Brazil, is predominant (Kämpf & Klamt, 1978). Besides, other clay minerals were identified, e.g., cristobalite, disordered kaolinite or halloysite or both, and hydroxy-Al interlayered 2:1 phyllosilicates (Kämpf & Klamt, 1978; Dalmolin et al., 2006). There is evidence for the occurrence of a small amount of hematite. The soils vary in texture and differ significantly in clay and iron oxide contents (Table 1). Nine disturbed and undisturbed soil samples were collected from two layers (0-10 and 10-20 cm) for each land use, to determine the intrinsic soil properties clay and iron oxide contents (Table 1). The clay content was assessed by the pipette method (Kilmer & Alexander, 1949), and total free iron oxide (FeDCB) extracted by hot dithionite-citrate-bicarbonate (Mehra & Jackson, 1960). The variable soil properties potential acidity, sum of bases (sum of Ca2+, Mg2+ and K+), organic carbon, bulk density, and porosity were analyzed. Potential acidity (H+Al) was determined by titration with NaOH after extraction with 0.5 mol L-1 Ca(CH3COO)2.H2O solution at pH 7.0-7.1; Ca2+ and Mg2+ were quantified by atomic absorption after 1.0 mol L-1 KCl extraction; K+ was quantified by flame atomic emission spectrometry after extraction with Mehlich-1 solution (Mehlich, 1953); and organic carbon was determined by titration with Fe(NH4)2(SO4)2.6H2O after wet digestion with 0.2 mol L-1 K2Cr2O7 solution (Yeomans & Bremner, 1988). All methods were described by Tedesco et al. (1995) and Donagema et al. (2011). Statistical analysis was performed in the R environment (R Core Team, 2013) using frequency histograms, cumulative frequency distribution graphs, and the Cullen and Frey graph to test data normality (Delignette-Muller et al., 2010; R Core Team, 2013). The homogeneity of variances among land uses was tested using box plots, standard deviations and Bartlett's Test (Bartlett, 1954; Webster, 2001). The distribution of the intrinsic soil properties was normal in both depth ranges. But for some variable soil properties (potential acidity and organic carbon in the upper, sum of bases in the lower layer) skewed distribution and heteroscedasticity were observed. Logarithmic transformation was used to achieve normality and stabilize the variances. The transformations altered the correlations between soil properties. The differences between dependent correlations (before and after transformation) were not significant at 5 % probability of a type I error using Williams's Test (Steiger, 1980; Revelle, 2012). To test the hypothesis that the intrinsic soil properties were the same for all three land uses, we used ANOVA (analysis of variance) and analyzed the results by Fisher's LSDs (least significant differences) for means separation at a probability of 5 % of committing a type I error (Fisher, 1951; Webster, 2007). To test the effect of clay and iron oxide content on variable soil properties, we used ANCOVA (analysis of covariance) when the correlation between soil properties was significant at 5 % probability of a type I error in the Student's t-tests (Revelle, 2012). The relation between intrinsic and variable soil properties was also analyzed using partial correlations (Revelle, 2012) and scatterplots with a 95 % confidence interval of a linear fit (Wickham, 2009; R Core Team, 2013). Tables and plots are presented according to the latest international recommendations for statistical results in soil science (Webster, 2001; Webster, 2007). Soil data and R scripts are freely available at <http://www.soil-scientist.net>. RESULTS The intrinsic soil properties were correlated with the variable soil properties (Table 2). Potential acidity is significantly correlated with clay and FeDCB in both layers, but the correlation was strongest in the subsurface. However, the partial correlation (r) between potential acidity and FeDCB, when the clay constant was fixed, r was only 0.252 in the surface and 0.084 in the subsurface layer. Besides, the data points were not evenly distributed in the 95 % confidence interval of the mean in the scatterplots (Figure 1), that is, they were stratified according to the land use. Removing data of any land use significantly altered the correlation between soil properties. An exception is the correlation of clay with potential acidity for the subsurface. The ANCOVA table (Table 3) showed that the effect of FeDCB on potential acidity was not significant. Clay also played a minor role, accounting for no more than 5 % of the total variance, while land use accounted for more than 80 % in both layers. Sum of bases was negatively correlated with clay and FeDCB content (Table 2). However, the partial correlation of the sum of bases with FeDCB in the 0-10 cm layer at a fixed clay content was only r = -0.295. In addition, data points were not evenly distributed in the 95 % confidence interval of the mean in the scatterplots (Figure 1). Again, removing data from any land use significantly altered the correlation between soil properties. The ANCOVA results (Table 3) showed a minor effect of the clay content on the sum of bases in the surface layer (10 % of the total variance), and no effect in the subsurface layers. Organic carbon correlated negatively with clay content in the 0-10 cm layer and with FeDCB in the 10-20 cm layer (Table 2). However, data points were stratified in the scatterplots according to land use (Figure 1). The ANCOVA (Table 3) showed a minor effect of clay on soil organic carbon, accounting for less than 5 % of the total variance. Moreover, the effect of FeDCB was not significant compared to the effect of the land uses, which accounted for 80 % of the total variance in the 0-10 cm layer and for 36 % in the 10-20 cm layer. Bulk density was positively correlated with clay content in the surface layer (Table 2). However, some degree of stratification of data points according to the land use was identified in the scatterplot (Figure 1). Besides, few data points were within the 95 % confidence interval of the mean. The ANCOVA showed that clay had no significant effect on soil bulk density (Table 4). Land use accounted for more than 80 % of the total variance. Porosity was negatively correlated with clay content (Table 2). In the scatterplot, we observed a marked stratification of data points according to the land use (Figure 1). Thus, removing data from any land use significantly altered the correlations between soil properties. The ANCOVA showed that clay did not affect soil porosity and that land use accounted for more than 85 % of the total variance (Table 4). DISCUSSION Scatterplots and ANCOVA were essential to identify the true correlations between soil properties. The use of correlation coefficients alone would have resulted in a misinterpretation of the relations between soil properties and led to an overestimation of the effect of intrinsic soil properties. In this study, all correlations were statistically significant, however not meaningful, e.g., the correlation between clay content and sum of bases. Phyllosilicate rich-soils, such as the Leptosols and Acrisols studied here, are known to have a net negative charge (Essington, 2004). In this case, a positive correlation would be expected between sum of bases and clay content. However, conversely, we observed a spurious correlation. An accidental correlation is one to which no "natural" causal interpretation can be applied (Haig, 2003). In this study, the spurious correlation observed was a result of the accidental occurrence of the highest sum of bases and lowest clay content in soils under no-tillage on the one hand, and the lowest sum of bases and highest clay content in soils under conventional agriculture on the other hand. According to the land-use history, soils under no-tillage are periodically treated with high fertilization rates, while in soils under conventional cultivation, only a small part of the nutrients exported by crops is replaced. On the other hand, in soils under native forest the sum of bases is slightly higher than under conventional agriculture because there is no export of nutrients. Besides, soils under native forest are Leptosols, which means that the weathering was less intense than of Acrisols, resulting in less leaching of bases from the soil profile (Chesworth, 2008). Because of all these relations, the variation in the sum of bases in these soils is a result of the land-use history (and maybe soil type). Clay and iron oxide content have no major influence, suggesting that they are not relevant in the assessment of soil quality. The correlation of intrinsic soil properties with organic carbon, bulk density and porosity were accidental as well. It is known that clay and iron oxides can play a positive role in organic carbon accumulation in soils (Dalmolin et al., 2006). We also expected a negative correlation of clay content with bulk density, in the same way as lower clay contents reduce soil porosity (Michelon et al., 2007). However, the correlation in the studied soils was opposite to the one expected. Similarly to the trend of the sum of bases, the variation observed in organic carbon, bulk density and porosity was mainly affected by the land-use history. Only the correlation between intrinsic soil properties and potential acidity was as expected (Tan, 2010). However, despite the statistically significant correlations, ANCOVA showed that clay content had only a minor effect on potential acidity and that FeDCB had none. The land use had the greatest effect on the potential acidity in both layers, as evidenced by the stratification of the points in the scatterplots. Once again there was an accidental correlation, where soils with lower clay content were under a land use (no-tillage) that reduced potential acidity. The soils under no-tillage were periodically limed. On the other hand, soils with higher clay content were under a land use (conventional agriculture) that induced the degradation of soil fertility. Under conventional agriculture, the soil was limed more than four years before sampling. In turn, soils under native forest were not fertilized or limed and consequently organic matter is the component with the greatest impact on availability of exchangeable aluminium (Silva et al., 2008) and thus on potential acidity. Since these soils have higher organic carbon content than soils under conventional agriculture, the potential acidity will be lower. Our results showed that under the tested conditions, soil properties important for land use such as potential acidity, sum of bases, organic carbon, bulk density, and porosity are affected by the land use. They are efficient in reflecting land use and are little affected by intrinsic soil properties such as clay and iron oxide content, suggesting that they can be included in soil quality assessments. Our data also indicated that despite statistically significant differences in clay and iron oxide content among soils, these soil properties have no effect on the results of soil quality assessments. Webster (2001) previously stated that when a difference is considered statistically significant, it is not necessarily important or physically or biologically meaningful. Our findings and those of many other researchers agree with this statement. Mandal et al. (2010) studied Himalayan Entisols and Inceptisols with clay contents between 106 and 172 g kg-1 and found that soil management strongly influenced the soil quality indicators. Marzaioli et al. (2010) studied sandy loam and loam sandy Mediterranean soils (clay content 54 - 171 g kg-1) and also observed that any differences in soil quality could be attributed to differences in land use. Similar results were found by Schipper et al. (2010), who expected to find fewer changes in organic carbon after 40 years in Andisols in New Zealand than in other soil orders due to the remarkable capacity of allophane and imogolite to sequester and protect organic carbon (Parfitt, 2009). Miralles et al. (2009) found that clay and iron oxide contents in soil samples of Leptosols, Calcisols, Vertisols, Luvisols, Kastanozems, and Chernozems (clay content between 219 and 483 g kg-1 and FeDCB content between 4.4 and 14.1 g kg-1) were weakly correlated with organic carbon (r < 0.26, p = 0.05, n = 68). Percival et al. (2000) studied New Zealand soils (Andisols, Inceptisols, Ultisols, Mollisols, Oxisols, Alfisols, Spodosols, Entisols, Aridisols, and Ultisols) with clay contents from 38 g kg-1 to 329 g kg-1 and found that carbon contents correlated poorly with clay content across and within soil orders. According to the authors, the clay content accounted for only 1 % of the variation in carbon content in both cases. However, several studies have shown that intrinsic soil properties can control the expression of other soil properties. In the above-mentioned study, Miralles et al. (2009) observed a strong correlation of clay and iron oxide content with cation exchange capacity, porosity and bulk density. Likewise, Dieckow et al. (2009) compared a clay Ferralsol (clay: 711 g kg-1; FeDCB: 87.6 g kg-1), a sandy clay Ferralsol (clay: 360 g kg-1; FeDCB: 48.3 g kg-1) and a sandy clay loam Acrisol (clay: 220 g kg-1; FeDCB: 11.8 g kg-1) and observed that clay and iron oxide contents were strongly correlated with organic carbon. According to Dieckow et al. (2009), in terms of organic carbon losses, the conversion of native vegetation to cropland is less harmful in clayey than in sandy soils, possibly because of a chemical and physical protection. The inconsistency of the findings of different researchers shows the need for further research on the role of intrinsic soil properties in the control of the effect of land uses on soil properties included in soil quality assessments. Meanwhile, we can suggest that textural classes should not be used as isolated tools to describe the heterogeneous character of soils and to explain the effects of intrinsic soil properties on other chemical and physical soil properties. Our study shows that despite the soil classification in different texture classes, the clay content has no apparent effect on other chemical and physical properties. Whenever possible, the information on intrinsic soil properties should be used as a continuous (clay content) and not as a discrete variable (textural class). Therefore, ANCOVA becomes highly relevant in the assessment of the extent to which chemical and physical soil properties can be explained by intrinsic soil properties under different land uses and thus affect soil quality assessments. Future studies must focus on the identification and definition of clay and iron oxide content ranges and other oxides in soils that affect the behavior of soil chemical and physical properties used in soil quality assessments in different ways. Finally, a general recommendation for studies on little weathered 1:1 and 2:1 phyllosilicate rich-soils with a similar mineralogy as in the soils of this study. If the difference between the maximum and the minimum clay content under different land uses is less than about 200 g kg-1 and the iron oxide content less than approximately 15 g kg-1, the physico-chemical soil properties in the surface layer (0-20 cm) are determined mostly by the land use and history. The effect of intrinsic soil properties is negligible and only detectable in deeper layers (>20 cm), where management practices have less effect. CONCLUSIONS 1. Organic carbon content, sum of bases, potential acidity, bulk density and porosity are easy measureable soil properties, reflect the land use and are little affected by intrinsic soil properties (clay and iron oxides). This suggests that they can be included in the assessment of the soil quality on small family farms. 2. Despite the classification of soils in different texture classes, intrinsic soil properties (clay and iron oxide contents) may have no effect on other chemical and physical soil properties. Whenever possible, information on intrinsic soil properties in soil quality assessments should be used as continuous rather than as discrete variables. 3. For little weathered 1:1 and 2:1 phyllosilicate rich-soils, if the difference between the maximum and the minimum clay content under different land uses is less than about 200 g kg-1 and the iron oxide content less than about 15 g kg-1, the physico-chemical soil properties in the surface layer are determined mostly by the land use. ACKNOWLEDGEMENTS We are grateful to Dr. Richard Webster, at the Rothamsted Research Centre, whose comments were decisive to improve our work. The authors gratefully acknowledge the support of the National Council for Scientific and Technological Development (CNPq) and of the Research Support Foundation of the State of Rio Grande do Sul (FAPERGS) grants to the authors and financial support for this research. LITERATURE CITED Received for publication on January 20, 2012 Approved on May 8, 2013. BARTLETT, M.S. A note on the multiplying factors for various chi-squared approximations. J. Royal Statist. Soc. Series B., 16:296-298, 1954. A note on the multiplying factors for various chi-squared approximations J. Royal Statist. Soc. Series B. 1954 296 298 16 BARTLETT M.S. BASHER, L.R. Is pedology dead and buried? Aust. J. Soil Res., 35:979-994, 1997. Is pedology dead and buried? Aust. J. Soil Res. 1997 979 994 35 BASHER L.R. BEWKET, W. & STROOSNIJDERB, L. Effects of agroecological land use succession on soil properties in Chemoga watershed, Blue Nile basin, Ethiopia. Geoderma, 111:85-98, 2003. Effects of agroecological land use succession on soil properties in Chemoga watershed, Blue Nile basin, Ethiopia Geoderma 2003 85 98 111 BEWKET W. STROOSNIJDERB L. CHESWORTH, W. Encyclopedia of earth sciences series. Dordrecht, Springer, 2008. 902p. Encyclopedia of earth sciences series 2008 CHESWORTH W. CLEMENS, G.; FIEDLER, S.; CONG, N.D.; van DUNG, N.; SCHULER, U. & STAHR, K. Soil fertility affected by land use history, relief position, and parent material under a tropical climate in NW-Vietnam. Catena, 81:87-96, 2010. Soil fertility affected by land use history, relief position, and parent material under a tropical climate in NW-Vietnam Catena 2010 87 96 81 CLEMENS G. FIEDLER S. CONG N.D. van DUNG N. SCHULER U. STAHR K. DALMOLIN, R.S.; GONÇALVES, C.N.; DICK, D.P.; KNICKER, H.; KLAMT, E. & KÖGEL-KNABNER, I. Organic matter characteristics and distribution in Ferralsol profiles of a climosequence in southern Brazil. Eur. J. Soil Sci., 57:644-654, 2006. Organic matter characteristics and distribution in Ferralsol profiles of a climosequence in southern Brazil Eur. J. Soil Sci. 2006 644 654 57 DALMOLIN R.S. GONÇALVES C.N. DICK D.P. KNICKER H. KLAMT E. KÖGEL-KNABNER I. DELIGNETTE-MULLER, M.L.; POUILLOT, R.; DENIS, J.-B. & DUTANG, C. Fitdistrplus: Help to fit of a parametric distribution to non-censored or censored data. Viena, R Project, 2010. 29p. Fitdistrplus: Help to fit of a parametric distribution to non-censored or censored data 2010 DELIGNETTE-MULLER M.L. POUILLOT R. DENIS J.-B. DUTANG C. DIECKOW, J.; BAYER, C.; CONCEIÇÃO, P.C.; ZANATTA, J.A.; MARTIN-NETO, L.; MILORI, D.B.M.; SALTON, J.C.; MACEDO, M.M.; MIELNICZUK, J. & HERNANI, L.C. Land use, tillage, texture and organic matter stock and composition in tropical and subtropical Brazilian soils. Eur. J. Soil Sci., 60:240-249, 2009. Land use, tillage, texture and organic matter stock and composition in tropical and subtropical Brazilian soils Eur. J. Soil Sci. 2009 240 249 60 DIECKOW J. BAYER C. CONCEIÇÃO P.C. ZANATTA J.A. MARTIN-NETO L. MILORI D.B.M. SALTON J.C. MACEDO M.M. MIELNICZUK J. HERNANI L.C. DONAGEMA, G.K.; CAMPOS, D.V.B.; CALDERANO, S.B.; TEIXEIRA, W.G. & VIANA, J.H.M. Handbook of methods of soil analysis. Rio de Janeiro, Embrapa Solos, 2011. 230p. Handbook of methods of soil analysis 2011 DONAGEMA G.K. CAMPOS D.V.B. CALDERANO S.B. TEIXEIRA W.G. VIANA J.H.M. DORAN, J.W. & PARKIN, T.B. Defining and assessing soil quality. In: DORAN, J.W.; COLEMAN, D.C.; BEZDICEK, D.F. & STEWART, B.A., eds. Defining soil quality for a sustainable environment. Madison, Soil Sience Socienty of America, 1994. p.3-21. Defining soil quality for a sustainable environment 1994 3 21 DORAN J.W. PARKIN T.B. DORAN J.W. COLEMAN D.C. BEZDICEK D.F. STEWART B.A. ESSINGTON, M.E. Soil and water chemistry: An integrative approach. Boca Raton, CRC Press, 2004. 534p. Soil and water chemistry: An integrative approach 2004 ESSINGTON M.E. FISHER, R.A. The design of experiments. London, Oliver and Boyd, 1951. 248p. The design of experiments 1951 FISHER R.A. HAIG, B.D. What is a spurious correlation? Underst. Stat., 2:125-132, 2003. What is a spurious correlation? Underst. Stat. 2003 125 132 2 HAIG B.D. GEISSEN, V. & GUZMAN, G.M. Fertility of tropical soils under different land use systems - a case study of soils in Tabasco, Mexico. Appl. Soil Ecol., 31:169-178, 2006. Fertility of tropical soils under different land use systems: a case study of soils in Tabasco, Mexico Appl. Soil Ecol. 2006 169 178 31 GEISSEN V. GUZMAN G.M. ISLAM, K.R. & WEIL, R.R. Soil quality indicator properties in mid-atlantic soils as influenced by conservation management. J. Soil Water Conserv., 55:69-78, 2000. Soil quality indicator properties in mid-atlantic soils as influenced by conservation management J. Soil Water Conserv. 2000 69 78 55 ISLAM K.R. WEIL R.R. IUSS WORKING GROUP WRB - IUSS-WRB. World reference base for soil resources 2006. 2.ed. Roma, FAO, 2006. 132p. World reference base for soil resources 2006 2006 KÄMPF, N. & KLAMT, E. Mineralogy and genesis of Oxisols and Ultisols from the Northeastern region of the Sul-riograndense Plateau. R. Bras. Ci. Solo, 2:68-73, 1978. Mineralogy and genesis of Oxisols and Ultisols from the Northeastern region of the Sul-riograndense Plateau R. Bras. Ci. Solo 1978 68 73 2 KÄMPF N. KLAMT E. KILMER, V.J. & ALEXANDER, L. Methods of making mechanical analyses of soils. Soil Sci., 68:15-24, 1949. Methods of making mechanical analyses of soils Soil Sci. 1949 15 24 68 KILMER V.J. ALEXANDER L. LARSON, W.E. & PIERCE, F.J. Conservation and enhancement of soil quality. In: INTERNATIONAL WORKSHOP ON EVALUATION FOR SUSTAINABLE LAND MANAGEMENT IN THE DEVELOPING WORLD, Chiang Rai, 1991. Proceedings... Bangkok, International Board for Soil Research and Management, 1991. p.175-203. Proceedings... 1991 175 203 LARSON W.E. PIERCE F.J. MANDAL, D.; SINGH, R.; DHYANI, S. & DHYANI, B. Landscape and land use effects on soil resources in a Himalayan watershed. Catena, 81:203-208, 2010. Landscape and land use effects on soil resources in a Himalayan watershed Catena 2010 203 208 81 MANDAL D. SINGH R. DHYANI S. DHYANI B. MALUF, J. A new climatic classification for the state of Rio Grande do Sul, Brazil. R. Bras. Agrometeorol., 8:141-150, 2000. A new climatic classification for the state of Rio Grande do Sul, Brazil R. Bras. Agrometeorol. 2000 141 150 8 MALUF J. MARZAIOLI, R.; D'ASCOLI, R.; PASCALE, R.D. & RUTIGLIANO, F. Soil quality in a Mediterranean area of Southern Italy as related to different land use types. Appl. Soil Ecol., 44:205-212, 2010. Soil quality in a Mediterranean area of Southern Italy as related to different land use types Appl. Soil Ecol. 2010 205 212 44 MARZAIOLI R. D'ASCOLI R. PASCALE R.D. RUTIGLIANO F. MEHLICH, A. Determination of P, Ca, Mg, K, Na and NH4. Raleigh, North Carolina Soil Testing Division, 1953. p.195. (Não Publicado) Determination of P, Ca, Mg, K, Na and NH4 1953 195 MEHLICH A. MEHRA, O.P. & JACKSON, M.L. Iron oxide removal from soils and clays by a dithionite-citrate system buffered with sodium bicarbonate. In: NATIONAL CONFERENCE ON CLAYS AND CLAY MINERALS, New York, 1960. Proceedings... New York, Pergamon, 1960. p.317-327. Proceedings... 1960 317 327 MEHRA O.P. JACKSON M.L. MICHELON, C.J.; CARLESSO, R.; PETRY, M.T.; DAVID, G. & DALLA SANTA, C. Physical quality of irrigated soil of Rio Grande do Sul State-Brazil. Ci. Rural, 37:1308-1315, 2007. Physical quality of irrigated soil of Rio Grande do Sul State-Brazil Ci. Rural 2007 1308 1315 37 MICHELON C.J. CARLESSO R. PETRY M.T. DAVID G. DALLA SANTA C. MIRALLES, I.; ORTEGA, R.; ALMENDROS, G.; SÁNCHEZ-MARAÑÓN, M. & SORIANO, M. Soil quality and organic carbon ratios in mountain agroecosystems of South-east Spain. Geoderma, 150:120-128, 2009. Soil quality and organic carbon ratios in mountain agroecosystems of South-east Spain Geoderma 2009 120 128 150 MIRALLES I. ORTEGA R. ALMENDROS G. SÁNCHEZ-MARAÑÓN M. SORIANO M. MUNIZ, L.C.; MADARI, B.E.; TROVO, J.B.F.; CANTANHÊDE, I.S.L.; MACHADO, P.L.O.A.; COBUCCI, T. & FRANÇA, A.F.S. Soil biological attributes in pastures of different ages in a crop-livestock integrated system. Pesq. Agropec. Bras., 46:1262-1268, 2011. Soil biological attributes in pastures of different ages in a crop-livestock integrated system Pesq. Agropec. Bras. 2011 1262 1268 46 MUNIZ L.C. MADARI B.E. TROVO J.B.F. CANTANHÊDE I.S.L. MACHADO P.L.O.A. COBUCCI T. FRANÇA A.F.S. PARFITT, R.L. Allophane and imogolite: role in soil biogeochemical processes. Clay Miner., 44:135-155, 2009. Allophane and imogolite: role in soil biogeochemical processes Clay Miner. 2009 135 155 44 PARFITT R.L. PERCIVAL, H.J.; PARFITT, R.L. & SCOTT, N.A. Factors controlling soil carbon levels in New Zealand grasslands: is clay content important? Soil Sci. Soc. Am. J., 64:623-1630, 2000. Factors controlling soil carbon levels in New Zealand grasslands: is clay content important? Soil Sci. Soc. Am. J. 2000 623 1630 64 PERCIVAL H.J. PARFITT R.L. SCOTT N.A. R CORE TEAM. R: A language and environment for statistical computing. Viena, R Project, 2011. R: A language and environment for statistical computing 2011 REICHERT, J.M.; REINERT, D.J. & BRAIDA, J.A. Qualidade dos solos e sustentabilidade de sistemas agrícolas. R. Ci. Amb., 27:29-48, 2003. Qualidade dos solos e sustentabilidade de sistemas agrícolas R. Ci. Amb. 2003 29 48 27 REICHERT J.M. REINERT D.J. BRAIDA J.A. REVELLE, W. Psych: Procedures for psychological, psychometric, and personality research. Viena, R Project, 2012. 282p. Psych: Procedures for psychological, psychometric, and personality research 2012 REVELLE W. SAMUEL-ROSA, A.; MIGUEL, P.; DALMOLIN, R.S.D. & PEDRON, F.A. Uso da terra no Rebordo do Planalto do Rio Grande do Sul. Ci. Natura, 33:161-173, 2011. Uso da terra no Rebordo do Planalto do Rio Grande do Sul Ci. Natura 2011 161 173 33 SAMUEL-ROSA A. MIGUEL P. DALMOLIN R.S.D. PEDRON F.A. SANTOS, H.G.; JACOMINE, P.K.T.; ANJOS, L.H.C.; OLIVEIRA, V.A.; OLIVEIRA, J.B.; COELHO, M.R.; LUMBRERAS, J.F. & CUNHA, T.J.F., eds. Sistema brasileiro de classificação de solos. 2.ed. Rio de Janeiro, Embrapa Solos, 2006. 306p. Sistema brasileiro de classificação de solos 2006 SANTOS H.G. JACOMINE P.K.T. ANJOS L.H.C. OLIVEIRA V.A. OLIVEIRA J.B. COELHO M.R. LUMBRERAS J.F. CUNHA T.J.F. SCHIPPER, L.; PARFITT, R.L.; ROSS, C.; BAISDEN, W.; CLAYDON, J. & FRASER, S. Gains and losses in C and N stocks of New Zealand pasture soils depend on land use. Agric. Ecosyst. Environ., 139:611-617, 2010. Gains and losses in C and N stocks of New Zealand pasture soils depend on land use Agric. Ecosyst. Environ. 2010 611 617 139 SCHIPPER L. PARFITT R.L. ROSS C. BAISDEN W. CLAYDON J. FRASER S. SILVA, L.B.; DICK, D.P. & INDA, A.V. Highland subtropical soils: Chemical attributes, content of organic matter and its resistance to chemical oxidation. Ci. Rural, 38:1167-1171, 2008. Highland subtropical soils: Chemical attributes, content of organic matter and its resistance to chemical oxidation Ci. Rural 2008 1167 1171 38 SILVA L.B. DICK D.P. INDA A.V. SOIL SURVEY STAFF. Keys to soil taxonomy. Washington DC: USDA-NRCS, 2010. 346p. Keys to soil taxonomy 2010 STEIGER, J. Tests for comparing elements of a correlation matrix. Psychol. Bull., 87:245-251, 1980. Tests for comparing elements of a correlation matrix Psychol. Bull. 1980 245 251 87 STEIGER J. TAN, K.H. Principles of soil chemistry. 4.ed. Boca Raton, CRC Press, 2010. 362p. Principles of soil chemistry 2010 TAN K.H. TEDESCO, M.J.; GIANELLO, C.; BISSANI, C.A.; BOHNEN, H. & VOLKWEISS, S.J. Analysis of soil, plants and other materials. Porto Alegre, Universidade Federal do Rio Grande do Sul, 1995. 174p. Analysis of soil, plants and other materials 1995 TEDESCO M.J. GIANELLO C. BISSANI C.A. BOHNEN H. VOLKWEISS S.J. WEBSTER, R. Statistics to support soil research and their presentation. Eur. J. Soil Sci., 52:331-340, 2001. Statistics to support soil research and their presentation Eur. J. Soil Sci. 2001 331 340 52 WEBSTER R. WEBSTER, R. Analysis of variance, inference, multiple comparisons and sampling effects in soil research. Eur. J. Soil Sci., 58:74-82, 2007. Analysis of variance, inference, multiple comparisons and sampling effects in soil research Eur. J. Soil Sci. 2007 74 82 58 WEBSTER R. WICKHAM, H. ggplot2: Elegant graphics for data analysis. New York, Springer, 2009. 212p. ggplot2: Elegant graphics for data analysis 2009 WICKHAM H. YEOMANS, J.C. & BREMNER, J.M. A rapid and precise method for routine determination of organic carbon in soil. Commun. Soil Sci. Plant Anal., 19:1467-1476, 1988. A rapid and precise method for routine determination of organic carbon in soil Commun. Soil Sci. Plant Anal. 1988 1467 1476 19 YEOMANS J.C. BREMNER J.M.