SciELO - Scientific Electronic Library Online

 
vol.51 issue9Physical quality of an Oxisol under no-tillage subjected to different cropping systemsSprinkler-irrigated rice yield in lowlands as a function of water availability and soil attributes author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

Share


Pesquisa Agropecuária Brasileira

Print version ISSN 0100-204XOn-line version ISSN 1678-3921

Pesq. agropec. bras. vol.51 no.9 Brasília Sept. 2016

http://dx.doi.org/10.1590/s0100-204x2016000900057 

SOIL SCIENCE AND PRODUCTION

Pore network of an Inceptisol under different uses and relativized S-index as an indicator of soil physical quality

Rede de poros de um Cambissolo sob diferentes usos e índice S relativizado como indicador da qualidade física do solo

Raimundo Nonato de Assis Júnior(1) 

Jaedson Cláudio Anunciato Mota(1) 

Alcione Guimarães Freire(1) 

Thiago Leite de Alencar(1) 

(1)Universidade Federal do Ceará, Centro de Ciências Agrárias, Departamento de Ciências do Solo, Campus do Pici, Avenida Mister Hull, no 2.977, CEP 60356-001 Fortaleza, CE, Brazil. E-mail: assisjr@ufc.br, jaedson.mota@ufc.br, alcionegf@hotmail.com, thiagoleitealencar@yahoo.com.br


Abstract

The objective of this work was to evaluate soil porous media quality under different uses, as well as to propose a new index to monitor changes in its structure due to management: the Srelative index. An Inceptisol under irrigated banana cultivation and secondary forest, in the Brazilian Apodi Plateau, was evaluated at three soil layers (0.0-0.1, 0.1-0.2, and 0.2-0.3 m) as to: particle size, particle and bulk density, organic matter content, intrinsic permeability to air, and water retention curve. Total porosity, index for pore continuity (defined by the ratio between soil intrinsic permeability to air and aeration porosity), volume of blocked pores, and the S and Srelative indices were determined. The Srelative index was calculated as the ratio between S-index values from undisturbed and disturbed samples of the same soil (Sundisturbed/Sdisturbed). Under the studied conditions, the Srelative index is a good quantifier of the changes imposed by soil use and management. Soil cultivation maintains or improves the quality of the evaluated attributes, except for blocked porosity, when compared with the soil under native secondary forest with ten years of regeneration after agricultural use.

Index terms: blocked porosity; irrigation management; pore continuity; soil management; soil permeability to air; Srelative index

Resumo

O objetivo deste trabalho foi avaliar a qualidade do meio poroso de solo submetido a diferentes usos, bem como propor um novo índice para monitorar mudanças em sua estrutura devidas ao manejo: o índice Srelativo. Um Cambissolo sob cultivo de banana irrigada e mata secundária, na Chapada do Apodi, foi avaliado em três camadas (0,0-0,1, 0,1-0,2 e 0,2-0,3 m), quanto a: granulometria, densidade de partículas e do solo, conteúdo de matéria orgânica, permeabilidade intrínseca ao ar e curva de retenção de água. Foram determinados porosidade total, índice de continuidade dos poros (definido pela razão entre a permeabilidade ao ar e a porosidade de aeração), volume de poros bloqueados, e índices S e Srelativo. O índice Srelativo foi calculado como a razão entre os valores do índice S em amostras indeformadas e deformadas de um mesmo solo (Sindeformada/Sdeformada). Nas condições estudadas, o índice Srelativo mostra-se como um bom quantificador das alterações impostas pelo uso e pelo manejo do solo. O cultivo do solo mantém ou melhora a qualidade dos atributos avaliados, exceto para porosidade bloqueada, em comparação ao solo sob mata secundária com dez anos de regeneração após uso agrícola.

Termos para indexação: porosidade bloqueada; manejo da irrigação; continuidade de poros; manejo do solo; permeabilidade ao ar; índice Srelativo

Introduction

Although it is largely accepted that soil quality is a consequence of biological, chemical, and physical aspects, the latter are often neglected. This, according to Dexter (2004), is a great misconception, since the decrease in soil physical quality, which is widespread in agricultural areas, affects both chemical and biological conditions. These soil conditions, moreover, require information on physical aspects in order to be fully interpreted. In this context, the inclusion of visual indicators is recommended for the assessment of soil quality (Giarola et al., 2010).

Therefore, the research effort in seeking indicators or reference values to analyze soil physical quality for agricultural purposes has been intense. The following critical values of soil physical attributes, for instance, have been suggested: 2 MPa for resistance to root penetration, 10% for aeration porosity (Beutler et al., 2006), and 1 µm2 for soil intrinsic permeability to air (Ball et al., 1988).

The S-index proposed by Dexter (2004) is commonly used as a soil quality indicator. Mathematically, it is defined as the slope of the tangent line to the inflection point of the soil water retention curve, traced according to the gravimetric water content (w) and to the natural logarithm of the soil water matric potential (ln h). Physically, it is associated with the soil structure, since it separates the effects of soil structure and texture in the soil water retention curve. Dexter & Czyz (2007) defined soil quality categories for the S-index, namely: S<0.020, very poor; 0.020≤S<0.035, poor; 0.035≤S<0.05, good; and S≥0.05, very good.

Implicitly, the S-index assumes unimodal pore radius distribution and, therefore, it cannot be considered for soils with bimodal or multimodal pore radius distribution (Carducci et al., 2011). For van Lier (2014), the theoretical basis of the S-index is conceptually interesting, albeit it should not create expectations that the S-index, like an absolute index, is able to predict soil physical quality by itself. The author also emphasizes that the commonly used limit value of S=0.035 has no basis under the various boundary conditions described in the literature.

It should be noted that small changes in the S-index may cause substantial change in the quality attributed to the soil, which may not correspond to reality. Therefore, it is reasonable to consider that the S-index and other indicators still need adjustments in order to become real indicators of changes in soil structure. Relativizing these indexes, as has been done for soil density (Zhao et al., 2010), may prove a good alternative to cope with these deficiencies. The relativized S-index, here called Srelative, is obtained from the ratio between the S-index, calculated from undisturbed and disturbed samples of the same soil - that is, with the same texture -, allowing for quantitative inference on the effects of management on soil structure. Among the advantages of the Srelative index, is the fact that it does not require indicative ranges of soil quality, as does the S-index, which makes it applicable to evaluate alterations in the structure of any soil, regardless of the uses and managements to which it is subjected to.

The objective of this work was to evaluate soil porous media quality under different uses, as well as to propose a new index to monitor changes in its structure due to management: the Srelative index.

Materials and Methods

The study area is located at the Apodi Plateau, in the state of Ceará, Brazil, and has been uninterruptedly cultivated with banana (Musa sp.) for eight years. A contiguous preserved area of secondary native forest, once used for agricultural purposes ten years before, was taken as check. The areas under banana cultivation and secondary native forest had 3 and 0.5 ha, respectively, and were 5 m apart from each other. The soil is classified as a Cambissolo (Inceptisol), and water and nutrient supply for the crop were delivered through drip fertigation.

Disturbed and undisturbed soil samples were taken in the 0.0-0.1, 0.1-0.2, and 0.2-0.3-m soil layers. Undisturbed soil samples were collected with an Uhland sampler, in steel rings (0.05-m height x 0.05-m diameter), with five replicates. Samplings were performed in the central part and in the four corners of the experimental areas, excluding the lateral borders.

The particle size analysis was performed with the pipette method, whereas total porosity was calculated according to the results of particle density and soil bulk density. Particle density was determined using a volumetric flask; and soil bulk density, using undisturbed soil samples. Organic carbon was also determined, according to the method proposed by Yeomans & Bremmer (1988). Soil permeability to air was obtained based on Silva et al. (2009): a volume of air with 1,000-Pa pressure was passed through the volumetric rings containing undisturbed soil samples, in equilibrium with the water tensions of 2, 6, 10, 33, and 100 kPa.

The soil water retention curve was obtained considering the water content at a saturation equal to soil porosity. At the low tension points (2, 4, 6, 8, and 10 kPa), water content was determined using a Haines funnel, and, at the higher ones (33, 100, 300, 700, 1,000, and 1,500 kPa), using a Richards pressure chamber with porous plate (Klute, 1986). The fitted curve was obtained using the mathematical model proposed by van Genuchten (1980). The curve of cumulative pore-size frequency, or pore distribution curve, was obtained by the ratio θ/α - in which θ is the volumetric water content and α is total porosity - for each matric potential applied in the construction of the soil water retention curve (Libardi, 2012).

The S-index was determined according to Dexter (2004), but the proposed Srelative indicator was used to assist in its interpretation. The S-index is calculated with the soil water retention curves obtained for undisturbed (regular method) and disturbed (for Srelative) soil samples, whereas the Srelative index is obtained from the ratio S-indexundisturbed/S-indexdisturbed.

The disturbed soil samples were air-dried after passing through a 10-mesh sieve, then they were accommodated in rings with the same dimensions as those used for collecting the undisturbed samples. The accommodation must be done in a way that the soil particles are naturally arranged in the rings, without using a predefined value for soil density. It should be pointed out that, since the material was sieved, the accommodation of the material inside the cylinder occurs in such a way that the pores in it resemble the textural pores, obtained from the natural arrangement of the individual particles of sand, silt, and clay.

Sieved soil was used in the present study, since this is the same material employed in most of the chemical and physical analyses, including the soil water retention curve at lower potentials, in which the effect of soil structure on water retention is negligible (Reichardt & Tim, 2012).

Values of soil intrinsic permeability to air (Kair) were related to aeration porosity values (εair) using the Kozeny-Carman equation, similarly to what was done in the procedure adopted by Ball et al. (1988): Kair = MεairN, which can be written as logKair = logM + Nlogεair, in which M and N are empirical fitting parameters. In this case, the N value was considered as a pore continuity index that reflects either the increase in Kair with increasing εair, or the decrease in pore tortuosity and surface area with increasing pore fraction available for flow. Aeration porosity corresponds to the difference between total porosity and the water content measured at each potential.

The limiting air permeability (intercept in the "x" axis or in the "log εair" axis) was defined as 1 µm2, a value suggested by Ball et al. (1988) as the volume of blocked pores in the soil (εb) that does not contribute to the convective transport of air. Therefore, εb is the value of εair below which the air flow through the soil ceases as a consequence of the discontinuity in the aeration pore network. The εb value is obtained with the following equation: εb = 10(-logM)/N.

Groenevelt et al. (1984) also proposed a pore continuity index that is based on the ratio between Kair and εair, i.e., K1 = Kairair.

As for the S-index, Dexter & Czyz (2007) recommended that it should be considered the slope of the tangent line to the inflection point of the van Genuchten equation, traced with the gravimetric water content (w) and the natural logarithm of the soil water matric potential (ln h), according to the equation:

The dataset was analyzed considering a completely randomized statistical design, in a split-plot arrangement (cultivation systems and layers), with five replicates. The following tests were applied: the Kolmogorov-Smirnov test to verify data normality, the F-test for the analysis of variance, and Tukey's test to compare means, all at 5% probability, using the Assistat program (Silva, 2013).

Results and Discussion

The 0.0-0.3-m soil layer showed 26% sand, 39% silt, and 35% clay under cultivation; however, under native secondary forest, the same layer showed values of 28, 31, and 41%, respectively. For the 0.0-0.1, 0.1-0.2, and 0.2-0.3-m soil layers, soil organic matter contents in the native forest were 62, 32, and 25 g kg-1, respectively, which were higher than the values found in the cultivated area, which were 30, 27, and 24 g kg-1, respectively.

For the cultivated soil, the relation between soil intrinsic permeability to air and matric potential was linear, whereas, for the secondary forest soil, it was logarithmic (Figure 1). Therefore, the same mathematical model cannot be applied to adjust the data for both areas. Costa et al. (2014), while evaluating a Cambissolo (Inceptisol) located in the same region of the present study, also observed that Kair data versus the matric potential fitted well to a logarithmic model. Indeed, Kair values increased with the reduction in soil water contents, in response to higher pore volumes and to the formation of preferential paths for air flow as smaller-sized pores were progressively drained (Silva et al., 2009).

Figure 1 Soil permeability to air (Kair) as a function of the water matric potential of an Inceptisol in the layers: A, 0.0-0.1 m; B, 0.1-0.2 m; and C, 0.2-0.3 m, under cultivation with banana (Musa sp.); and D, 0.0-0.1 m; E, 0.1-0.2 m; and F, 0.2-0.3 m, under native secondary forest. 

For the cultivated area, at least in the studied water range from the established matric potentials, the less wet the soil gets, the greater is the area available for gas diffusion. However, in the area under secondary native forest, soil drying does not increase soil permeability to air in the same proportion, showing increasingly smaller increments of available area for air flow, for the same variation in the soil water matric potential.

The linear ratio between soil intrinsic permeability to air (log Kair) and aeration porosity (logεair) showed good adjustment (R2>0.84) for both soils and all layers (Figure 2). Similar relationships were reported by Ball et al. (1988), Dörner & Horn (2006), and Berisso et al. (2013). It should be highlighted that the angular coefficient (N) from the log-log relationship between the Kair and εair variables is considered an indicator of soil physical quality. Physically, it represents the rate of increase or reduction of connected pores as aeration porosity increases or decreases (Ball et al., 1988). Comparatively, the slope was greater in the soil under banana, showing a better pore connection.

Figure 2 Relationship between soil intrinsic permeability to air (log Kair) and aeration porosity (log εair) for an Inceptisol under cultivation with banana (Musa sp.) and under native secondary forest, in the layers: A, 0.0-0.1 m; B, 0.1-0.2 m; and C, 0.2-0.3 m. 

The cultivated soil, despite its better pore network quality, showed a very small value for log Kair, on the order of 6 µm2 (Figure 2), which is close to the limit of blocked porosity, of 1 µm2 (Berisso et al., 2013), when log εair was 0.1 m3 m-3 - the limiting aeration porosity for plants (Beutler et al., 2006). This suggests the need for a careful monitoring of the porous fraction of this soil. Similarly, Rodrigues et al. (2011), when analyzing conventional tillage and no-tillage systems in a Latossolo Vermelho (Rhodic Hapludox), found very low values for log Kair, even closer to the blocking situation than the one reported here.

There was a significant interaction of log Kair values between treatments, at the matric potentials of -33 and -100 kPa (Table 1). At the matric potential of -2 kPa, the pore volume exists but is blocked (log Kair = 1 µm2) and does not participate in the convective air transport (although soil gas exchanges occur predominantly through diffusion), in both systems and in all layers analyzed.

Table 1 Soil intrinsic permeability to air (Kair) and pore continuity index (K1) for an Inceptisol under banana (Musa sp.) cultivation and under native secondary forest, at different depths(1)

(1)Means followed by equal letters, uppercase in the rows and lowercase in the columns, do not differ by Tukey’s test, at 5% probability. **Significant interaction between the factors, at 1% probability. nsNonsignificant interaction between the factors.

This result is specifically important for the cultivated soil, because, when it is very close to saturation, the pores, although empty, do not participate in gas exchanges since they are blocked. This finding indicates a need for an irrigation management that provides water in the shortest time possible or, preferentially, at a rate that does not increase soil water content to an energy level higher than -2 kPa.

The pore continuity index K1, obtained by the ratio between Kair and εair, allows identifying whether the differences in the first are only due to the differences in the second, or whether they are a consequence of the porous space geometry, which depends on diameter distribution, tortuosity, and continuity of pores. For Groenevelt et al. (1984), soils with different Kair, but similar K1 values, show similar diameter and pore continuity distribution.

In the matric potential interval from -2 to -10 kPa, pore continuity was similar for both soil use systems and even for the different soil layers. However, at the matric potentials of -33 and -100 kPa, except for the layer of 0.1-0.2 m, the pores were either more continuous or less tortuous in the soil cultivated with banana than in the one under secondary forest. This result is similar to those found by Alencar (2014).

It is reasonable to consider that soil intrinsic permeability to air precedes εair regarding limitations for plant growth, since a volume of 0.1 m3 m-3 of unconnected pores does not meet the minimal condition for proper gas exchange between soil and atmosphere. Actually, Ball et al. (1988) have already reported that the continuity of conductive pores is probably more important than their volume or number. Although the porosity block limit for air flow is already known (1 µm2), soil physicists must research further the soil intrinsic permeability to air that would allow gas exchange and favor, or at least not restrain, plant growth. This variable could be useful, for instance, to determine the least limiting water range (Beutler et al., 2006).

For the cultivated soil, N values, according to the second equation, were higher than the ones observed in the secondary forest (Table 2), in all layers, besides showing greater pore continuity, as previously discussed. It was evident that the better pore network quality of the soil under cultivation was certainly due to wetting (when irrigated) and to drying (between irrigations), which contribute to the formation of cracks between aggregates, creating preferential paths for air flow. Furthermore, one should not disregard the role of the root system in the formation of root pores. Reports on root effects on soil permeability to air are given by Groenevelt et al. (1984), Uteau et al. (2013), and Chen et al. (2014).

Table 2 Regression parameters for the relationship between soil intrinsic permeability to air and soil air porosity(1) and the volume of pores blocked for air flow (εb), for an Inceptisol under banana (Musa sp.) cultivation and under native secondary forest, at different depths. 

(1)Described by the model: log Kair = log M + N log εair, in which Kair is the soil air intrinsic permeability; εair is the soil air porosity; and M and N are empirical fitting parameters.

The soil under native secondary forest showed lower volume of blocked pores in all soil layers, when compared with the soil cultivated with banana, which corroborates the results obtained by Alencar (2014) and by Zúñiga et al. (2015). Particularly in the 0.0-0.1-m soil layer, the soil under secondary native forest showed almost the totality of its pores available for gas flow, with only 0.39% of blocked pores. Despite being out of the scope of the present study, laboratorial analyses showed that, in this layer, the organic carbon content was statistically higher in soil under native secondary forest, which surely contributed to a better porous media quality, especially regarding pore continuity.

Soil water retention and pore distribution frequency are shown in Figure 3. Considering the pores with diameter ≥50 µm (Libardi, 2012), measured at a tension of 6 kPa, the macroporosity of the cultivated soil was: 0.036 m3 m-3 (7.4% of total porosity), in the 0.0-0.1-m layer; 0.044 m3 m-3 (9% of total porosity) in the 0.1-0.2-m layer; and 0.034 m3 m-3 (7.1% of total porosity) in the 0.2-0.3-m layer. For the soil under secondary forest, these values were: 0.055 m3 m-3 (10.9% of total porosity); 0.038 m3 m-3 (8.1% of total porosity); and 0.050 m3 m-3 (10.7% of total porosity), respectively.

Figure 3 Soil water retention (A) and pore frequency distribution (B) for an Inceptisol under cultivation with banana (Musa sp.) and under native secondary forest, in the layers: 0.0-0.1, 0.1-0.2, and 0.2-0.3 m. 

The soil under native secondary forest had a greater volume of pores, except for the layer of 0.1-0.2 m. However, as previously mentioned, this greater volume did not contribute to a higher intrinsic permeability to air at the tension of 6 kPa, because the pore continuity was very similar in both cases. Bouma (1982) observed that large pores may not contribute to the flow when they are discontinuous, while small pores can be good conductors when they have continuity.

The S-index did not show significant differences in soil quality between management systems, either for the same soil layer or between layers (Figure 4). Considering the classes proposed by Dexter & Czyz (2007), soil quality was characterized as very good, in all situations of soil use and depth. Regarding Srelative, both areas were statistically similar in the 0.0-0.1, 0.1-0.2, and 0.0-0.3-m soil layers, but different in the 0.2-0.3-m soil layer, in which the soil under cultivation showed improved structure when compared to the secondary forest soil. In addition, in the comparison between layers, the Srelative index detected that the layer of 0.0-0.1 m had significantly better structure than the other two deeper layers. Therefore, the Srelative index was more sensitive than the S-index to distinguish changes in soil physical quality.

Figure 4 S-index and Srelative index for an Inceptisol under cultivation with banana (Musa sp.) and under native secondary forest. Means followed by equal letters, uppercase between layer means and lowercase between soil use, do not differ by Tukey's test, at 5% of probability. 

Besides having a higher sensitivity than the S-index, the Srelative index was also able to quantify the imposed changes to soil structure, because it relativizes them to a standard (textural attributes) that is kept almost unchanged along time. Therefore, once the S value is obtained for the disturbed reference sample of a given soil, it is possible to obtain the Srelative index at any time in the future, after determining the S value of different management conditions.

The Srelative index is promising for the evaluation of soil structure quality, and further studies should include local observation in other soil types, managements, and climates, in order to assess the universality of its application.

Conclusions

  1. The Srelative index is a sensitive indicator of the effects of soil management on its physical quality.

  2. Soil cultivation with banana (Musa sp.) improves, or at least maintains, soil quality attributes, except for blocked porosity, when compared with the soil under native secondary forest with ten years of regeneration after agricultural use.

Acknowledgements

To Fundação Cearense de Apoio ao Desenvolvimento Científico e Tecnológico (Funcap, project No. 0024-00013.01.00/10), for financial support.

References

ALENCAR, T.L. Alterações físicas em um Cambissolo tratado com biofertilizante: indicadores de qualidade e refinamento do método de avaliação pelo Srelativo2014. 92p. Dissertação (Mestrado) - Universidade Federal do Ceará, Fortaleza. [ Links ]

BALL, B.C.; O'SULLIVAN, M.F.; HUNTER, R. Gas diffusion, fluid flow and derived pore continuity indices in relation to vehicle traffic and tillage. Journal of Soil Science, v.39, p.327-339, 1988. DOI: 10.1111/j.1365-2389.1988.tb01219.x. [ Links ]

BERISSO, F.E.; SCHJØNNING, P.; KELLER, T.; LAMANDÉ, M.; SIMOJOKI, A.; IVERSEN, B.V.; ALAKUKKU, L.; FORKMAN, J. Gas transport and subsoil pore characteristics: anisotropy and long-term effects of compaction. Geoderma, v.195-196, p.184-191, 2013. DOI: 10.1016/j.geoderma.2012.12.002. [ Links ]

BEUTLER, A.N.; CENTURION, J.F.; SILVA, A.P. da; BARBOSA, J.C. Intervalo hídrico ótimo e produtividade de cultivares de soja. Revista Brasileira de Engenharia Agrícola e Ambiental, v.10, p.639-645, 2006. DOI: 10.1590/S1415-43662006000300015. [ Links ]

BOUMA, J. Measuring the hydraulic conductivity of soil horizons with continuous macropores. Soil Science Society of America Journal, v.46, p.438-441, 1982. DOI: 10.2136/sssaj1982.03615995004600020047x. [ Links ]

CARDUCCI, C.E.; OLIVEIRA, G.C. de; SEVERIANO, E. da C.; ZEVIANI, W.M. Modelagem da curva de retenção de água de Latossolos utilizando a equação duplo van Genuchten. Revista Brasileira de Ciência do Solo, v.35, p.77-86, 2011. DOI: 10.1590/S0100-06832011000100007. [ Links ]

CHEN, G.; WEIL, R.R.; HILL, R.L. Effects of compaction and cover crops on soil least limiting water range and air permeability. Soil and Tillage Research, v.136, p.61-69, 2014. DOI: 10.1016/j.still.2013.09.004. [ Links ]

COSTA, P.A. da; MOTA, J.C.A.; ROMERO, R.E.; FREIRE, A.G.; FERREIRA, T.O. Changes in soil pore network in response to twenty-three years of irrigation in a tropical semiarid pasture from northeast Brazil. Soil and Tillage Research, v.137, p.23-32, 2014. DOI: 10.1016/j.still.2013.11.004. [ Links ]

DEXTER, A.R. Soil physical quality. Part I. Theory, effects of soil texture, density, and organic matter, and effects on root growth. Geoderma,v.120, p.201-214, 2004. DOI: 10.1016/j.geoderma.2003.09.004. [ Links ]

DEXTER, A.R.; CZYZ, E.A. Applications of S-theory in the study of soil physical degradation and its consequences. Land Degradation and Development, v.18, p.369-381, 2007. DOI: 10.1002/ldr.779. [ Links ]

DÖRNER, J.; HORN, R. Anisotropy of pore functions in structured Stagnic Luvisols in the Weichselian moraine region in N Germany. Journal of Plant Nutrition and Soil Science, v.169, p.213-220, 2006. DOI: 10.1002/jpln.200521844. [ Links ]

GIAROLA, N.F.B.; SILVA, A.P. da; TORMENA, C.A.; BALL, B.; ROSA, J.A. Visual soil structure quality assessment on Oxisols under no-tillage system. Scientia Agricola, v.67, p.479-482, 2010. DOI: 10.1590/S0103-90162010000400016. [ 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. DOI: 10.1016/0016-7061(84)90016-8. [ Links ]

KLUTE, A. Water retention: laboratory methods. In: KLUTE, A. (Ed.). Methods of soil analysis. 2nd ed. Madison: American Society of Agronomy, Soil Science Society of America, 1986. p.635-662. (SSSA. Agronomy monography, 9). DOI: 10.2136/sssabookser5.1.2ed.c26. [ Links ]

LIBARDI, P.L. Dinâmica da água no solo. 2.ed. Piracicaba: Edusp, 2012. 352p. [ Links ]

REICHARDT, K.; TIMM, L.C. Solo, planta e atmosfera: conceitos, processos e aplicações. 2.ed. Barueri: Manole, 2012. 524p. [ Links ]

RODRIGUES, S.; SILVA, A.P. da; GIAROLA, N.F.B.; ROSA, J.A. Permeabilidade ao ar em Latossolo Vermelho sob diferentes sistemas de manejo. Revista Brasileira de Ciência do Solo, v.35, p.105-114, 2011. DOI: 10.1590/S0100-06832011000100010. [ Links ]

SILVA, A.P. da; LEÃO, T.P.; TORMENA, C.A.; GONÇALVES, A.C.A. Determinação da permeabilidade ao ar em amostras indeformadas de solo pelo método da pressão decrescente. Revista Brasileira de Ciência do Solo, v.33, p.1535-1545, 2009. DOI: 10.1590/S0100-06832009000600003. [ Links ]

SILVA, F.A.S. Assistat software. 2013. Available at: <http://www.assistat.com/indexp.html>. Accessed on: Dec. 8 2015. [ Links ]

UTEAU, D.; PAGENKEMPER, S.K.; PETH, S.; HORN, R. Root and time dependent soil structure formation and its influence on gas transport in the subsoil. Soil and Tillage Research, v.132, p.69-76, 2013. DOI: 10.1016/j.still.2013.05.001. [ Links ]

VAN GENUCHTEN, M. Th. A closed-form equation for predicting the conductivity of unsaturated soils. Soil Science Society of America Journal, v.44, p.892-897, 1980. DOI: 10.2136/sssaj1980.03615995004400050002x. [ Links ]

VAN LIER, Q. de J. Revisiting the S-index for soil physical quality and its use in Brazil. Revista Brasileira de Ciência do Solo, v.38, p.1-10, 2014. DOI: 10.1590/S0100-06832014000100001. [ Links ]

YEOMANS, J.C.; BREMMER, J.M. A rapid and precise method for routine determination of organic carbon in soil. Communications in Soil Science and Plant Analysis, v.19, p.1467-1476, 1988. DOI: 10.1080/00103628809368027. [ Links ]

ZHAO, Y.; KRZIC, M.; BULMER, C.E.; SCHMIDT, M.G.; SIMARD, S.W. Relative bulk density as a measure of compaction and its influence on tree height. Canadian Journal of Forest Research, v.40, p.1724-173, 2010. DOI: 10.1139/X10-115. [ Links ]

ZÚÑIGA, F.; IVELIC-SÁEZ, J.; LÓPEZ, I.; HUYGENS, D.; DÖRNER, F.J. Temporal dynamics of the physical quality of an Andisol under a grazing system subjected to different pasture improvement strategies. Soil and Tillage Research, v.145, p.233-241, 2015. DOI: 10.1016/j.still.2014.09.014. [ Links ]

Received: June 25, 2015; Accepted: February 12, 2016

Creative Commons License This is an Open Access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.