Acessibilidade / Reportar erro

BEYOND THE “LEAST LIMITING WATER RANGE”: RETHINKING SOIL PHYSICS RESEARCH IN BRAZIL(1) (1) A version of this text in Portuguese may be obtained with the corresponding author.

ALÉM DO “INTERVALO HÍDRICO ÓTIMO”: REPENSANDO A PESQUISA EM FÍSICA DO SOLO NO BRASIL

Abstracts

As opposed to objective definitions in soil physics, the subjective term “soil physical quality” is increasingly found in publications in the soil physics area. A supposed indicator of soil physical quality that has been the focus of attention, especially in the Brazilian literature, is the Least Limiting Water Range (RLL), translated in Portuguese as "Intervalo Hídrico Ótimo" or IHO. In this paper the four limiting water contents that define RLLare discussed in the light of objectively determinable soil physical properties, pointing to inconsistencies in the RLLdefinition and calculation. It also discusses the interpretation of RLL as an indicator of crop productivity or soil physical quality, showing its inability to consider common phenological and pedological boundary conditions. It is shown that so-called “critical densities” found by the RLL through a commonly applied calculation method are questionable. Considering the availability of robust models for agronomy, ecology, hydrology, meteorology and other related areas, the attractiveness of RLL as an indicator to Brazilian soil physicists is not related to its (never proven) effectiveness, but rather to the simplicity with which it is dealt. Determining the respective limiting contents in a simplified manner, relegating the study or concern on the actual functioning of the system to a lower priority, goes against scientific construction and systemic understanding. This study suggests a realignment of the research in soil physics in Brazil with scientific precepts, towards mechanistic soil physics, to replace the currently predominant search for empirical correlations below the state of the art of soil physics.

modeling; field capacity; drought stress; mechanical stress; anoxic stress


Em oposição a definições objetivas da física do solo, o termo subjetivo “qualidade física do solo” aparece com frequência cada vez maior em trabalhos publicados da área. Um suposto indicador de qualidade física do solo que recebe muita atenção, especialmente na literatura brasileira, é o Intervalo Hídrico Ótimo (IHO), definido por quatro teores-limite de água. Nesse texto, discutem-se seus quatro teores-limite à luz de propriedades físicas do solo determináveis objetivamente, apontando-se incoerências na definição e no cálculo do IHO. Discute-se a interpretação do IHO como indicador da produtividade das culturas ou da qualidade física do solo, demonstrando-se sua incapacidade de abranger condições de contorno fenológicas e pedológicas comuns. Demonstra-se, também, que as densidades críticas encontradas pelo IHO por método comumente empregado são questionáveis. Considerando a disponibilidade de modelos comprovadamente robustos para as áreas de agronomia, ecologia, hidrologia, meteorologia e outras conexas, a popularidade do IHO como indicador na física do solo brasileira não pode ser entendido pela sua eficácia, pois essa nunca foi comprovada, mas deve ser explicado pela simplicidade com que ele é tratado. A determinação dos respectivos teores-limite de forma simplificada, deixando o estudo ou a preocupação com o real funcionamento do sistema para o segundo plano, vai à contramão da construção científica e do entendimento sistêmico. Sugere-se um realinhamento da pesquisa em física do solo no Brasil com os preceitos da ciência, na direção da física do solo mecanística em detrimento da busca por correlações empíricas aquém do estado da arte em física do solo.

modelagem; capacidade de campo; estresse hídrico; estresse mecânico; estresse anóxico


INTRODUCTION

Soil physics is the branch of soil science that deals with soil physical properties and processes. The resulting knowledge is used to predict the behavior of natural or managed ecosystems. To this effect, soil physics focuses on measuring and modeling the dynamics of physical soil components, comprised of mineral and organic solids, water, solutes and soil air, as well as heat flows in the soil system. A dissociation of measurement and modeling generates a limited and fragmented understanding of the systems of interest to soil physics and it is common practice that soil physicists develop instruments to measure the properties of interest, construing them in the light of hypotheses on the operation of the studied system. Soil physicists translate the evolution of hypotheses and the increased associated knowledge into models or algorithms, allowing a quantitative prediction of system variables in time and space.

As opposed to objective definitions of soil physics, the subjective term “soil physical quality” has appeared more and more often in studies published in the area of soil physics. Whereas “soil physics” is a clearly defined field of science, “soil physical quality” has no absolute definition. In contrast with objectively defined soil properties like conductivities and diffusivities with exact dimensions following from their physical definition, “quality” has no defined dimension. From the agricultural point of view, the inherent subjectivity of the term “soil physical quality” can be translated objectively by soil physics into properties of transfer and storage of mass and energy that correspond to contents of water, solutes, air and heat. These contents will then be appropriate to maximize the development of crops, minimize environmental degradation, and ensure soil structural stability to maintain its biological health and allow root growth.

Specific models have been developed aiming to predict and assess such processes. These models are primarily based on the physical and physicochemical description of mass and energy transfer processes in the soil and in the soil-plant-atmosphere system and allow the quantitative assessment of the adequacy of soil conditions to crop requirements. The development, calibration, validation, sensitivity analysis and the use of such models have engaged soil scientists worldwide. After defining a scenario, it is translated into input parameters for the respective model, allowing prediction of the system behavior subject to the established boundary conditions. Examples of well-known models of this type are Hydrus (Simunek and Van Genuchten, 2008Simunek J, van Genuchten MT. Modeling nonequilibrium flow and transport processes using HYDRUS. Vadose Zone J. 2008;7:782-97.),SWAP (Kroes et al., 2008Kroes JG, Van Dam JC, Groenendijk P, Hendriks RFA, Jacobs CMJ. SWAP version 3.2. Theory description and user manual. Wageningen: 2008. (Alterra Report, 1649). Disponível em: http://www.swap.alterra.nl.
http://www.swap.alterra.nl...
) and SWAT (Arnold et al., 2012Arnold JG, Moriasi DN, Gassman PW, Abbaspour KC, White MJ, Srinivasan R, Santhi C, Harmel RD, Van Griensven A, Van Liew MW, Kannan N, Jha MK. SWAT: Model use, calibration, and validation. Trans Am Soc Agric Biol Eng. 2012;55:1491-508.), among many others. The models differ in process description and the dimensions of the simulated system (1-D, 2-D or 3-D), the method of solving numerical problems (implicitly or explicitly) or the emphasis given to one or another sub-process. Under this systemic approach, soil physicists have advanced in objectively understanding how the system works, something that is intangible under the subjective approach of “soil physical quality”.

In an attempt to avoid dealing with the complex relations that govern the system and the assessment of the input parameters required by the respective models, indicators of the system status have been suggested to address more specific or simpler manifestations of a system. An indicator cannot only provide information on the progress of a system towards a given goal, but can also be understood as a tool that helps perceive a tendency or a phenomenon that, otherwise, would not be easily detected (Hammond et al., 1995Hammond A, Adriaanse A, Rodenburg E, Bryant D, Woodward R. Environmental indicators: A systematic approach to measuring and reporting on environmental policy performance in the context of sustainable development. Washington, D.C.: World Resources Institute; 1995.). Indicators sometimes just represent a characteristic of the system or are calculated as a function of some characteristics by an empirical formula. By name, “indicator” means something that indicates, suggests, and one cannot expect that an indicator will determine or provide predictions of system behavior. Thus, several indicators can be used to understand something as complex as “soil physical quality”. Soil aggregate stability (Jastrow et al., 1998Jastrow JD, Miller RM, Lussenhop J. Contributions of interacting biological mechanisms to soil aggregate stabilization in restored prairie. Soil Biol Biochem. 1998;30:905-16.; Saygin et al., 2012Saygin SD, Cornelis WM, Erpul G, Gabriels D. Comparison of different aggregate stability approaches for loamy sand soils. Appl Soil Ecol. 2012;54:1-6.), bulk density (Reynolds et al., 2009Reynolds WD, Drury CF, Tan CS, Fox CA, Yang XM. Use of indicators and pore volume-function characteristics to quantify soil physical quality. Geoderma. 2009;152:252-63.), organic matter content (Gilley et al., 2001Gilley JE, Doran JW, Eghball B. Tillage and fallow effects on selected soil quality characteristics of former conservation reserve program sites. J Soil Water Conserv. 2001;56:126-32.), available water capacity (Allen et al., 1998Allen RG, Pereira LS, Raes DE, Smith M. Crop evapotranspiration: Guidelines for computing crop water requirements. Rome: FAO; 1998. (Irrigation and Drainage Paper, 56).),S index (Dexter, 2004Dexter AR. Soil physical quality: Part I. Theory, effects of soil texture, density, and organic matter, and effects on root growth. Geoderma. 2004;120:201-14.), among many others, have all been proposed as indicators of soil physical quality in some context.

The use of these properties as indicators of soil physical quality is implicitly empirical and has little to do with the description of soil physics given in the first paragraph of this introduction. Clearly, the “quality of the indicator” depends on its effectiveness in representing the attributes to which the robust model (read as “state of the art”) indicates sensitivity.

An indicator that has received much attention, especially in Brazil, is the Least Limiting Water Range (Silva et al., 1994Silva AP, Kay BD, Perfect E. Characterization of the least limiting water range. Soil Sci Soc Am J. 1994;58:1775-81.). It was introduced in Brazilian literature by Tormena et al. (1998)Tormena CA, Silva AP, Libardi PL. Caracterização do intervalo hídrico ótimo de um Latossolo Roxo sob plantio direto. R Bras Ci Solo. 1998;22:573-81. and translated to Portuguese as “Intervalo Hídrico Ótimo” (IHO). The Least Limiting Water Range (RLL) is a mathematical proposition for the Non Limiting Water Range (NLWR), introduced by Letey (1985)Letey J. Relationship between soil physical properties and crop production. Adv Soil Sci. 1985;1:277-94., who defined it as the range of soil water contents for which neither the matric potential nor aeration nor mechanical impedance would be limiting factors to the plants growth. In its conception and in the way it was introduced by Letey (1985)Letey J. Relationship between soil physical properties and crop production. Adv Soil Sci. 1985;1:277-94., NLWR is a complex, comprehensive concept, to the extent that it comprises several aspects related to the soil physical conditions necessary for plant development. It does not, however, include hydraulic conductivity, closely associated with water availability (Gardner, 1960Gardner WR. Dynamic aspects of water availability to plants. Soil Sci. 1960;89:63-7.; Cowan, 1965Cowan IR. Transport of water in the soil-plant-atmosphere system. J Appl Ecol. 1965;2:221-39.), nor soil thermal properties, also influenced by the water content, which can determine plant growth (Abdelhafeez et al., 1975Abdelhafeez AT, Harssema H, Verkerk K. Effects of air temperature, soil temperature and soil moisture on growth and development of tomato itself and grafted on its own and egg-plant rootstock. Sci Hortic. 1975;3:65-73.; Reddell et al., 1985Reddell P, Bowen GD, Robson AD. The effects of soil temperature on plant growth, nodulation and nitrogen fixation in Casuarina cunninghamiana Miq. New Phytol. 1985;101:441-50.). Silva et al. (1994)Silva AP, Kay BD, Perfect E. Characterization of the least limiting water range. Soil Sci Soc Am J. 1994;58:1775-81. replaced the fundamentals of the NLWR soil-root interface processes by empirical limits of soil physical properties to establish a mathematical proposition for RLL.

Similarly to what is seen for the S index proposed by Dexter (2004)Dexter AR. Soil physical quality: Part I. Theory, effects of soil texture, density, and organic matter, and effects on root growth. Geoderma. 2004;120:201-14. and discussed by De Jong van Lier (2014)De Jong van Lier Q. Revisiting the S-index for soil physical quality and its use in Brazil. R Bras Ci Solo. 2014;38:1-10., Brazil leads the world ranking in number of published research papers onRLL (Gubiani et al., 2013a)Gubiani PI, Reichert JM, Reinert DJ. Indicadores hídrico-mecânicos de compactação do solo e crescimento de plantas. R Bras Ci Solo. 2013a;37:1-10.. This indicates a large number of human resources involved in determining theRLL, and therefore not engaged in studying fundamental physical properties, how to model them, and, based on that, how to make sound predictions on the functioning of the soil-plant-atmosphere system. Aiming to investigate the scope of the information contained in theRLL and its reliability to indicate “quality, upon which it was based, this work aims to analyze its limiting contents and compare them with basic soil physical properties. With this analysis, the authors invite students and researchers of soil physics to become acquainted with an objective interpretation of theRLL and suggest a realignment of the soil physics research in Brazil with scientific precepts.

DEVELOPMENT

Definition of the indicator

The Least Limiting Water Range (RLL, m3 m-3) is defined as:

where “max” represents the “maximum” function (max(i,j) = i ifij; max(i,j) = j ifi<j), and “min” is the “minimum” function (min(i,j) = i ifij; min(i,j) = j ifi>j). The definition ofRLL includes four limiting soil water contents, all on a volume base (m3 m-3):θfc is the water content at field capacity; θair, is the limiting water content for proper aeration;θpwp is the water content at the permanent wilting point; andθpr, is the soil water content corresponding to the limiting penetration resistance. TheRLL indicator can be compared to the quantity available water A (m3m-3), defined as:

Figure 1 shows the graphical representation of RLL and A,according to equations 1 and 2 for the four possible combinations of the relative position of θpwp andθpr and ofθair andθfc. Awill be equal to RLL in the case ofθpwpθprandθfcθair(case 1, Figure 1). In all other cases,RLL will be smaller thanA.

Figure 1
Least Limiting Water Range (RLL) and available water (A) on a scale of volumetric water contents from completely dry (0) to saturation (θs), according to equations 1 and 2 for the four possible combinations of the relative position of water contents of permanent wilting and limiting root penetration resistance (θpwp andθpr) and of water contents of limiting air-filled porosity and field capacity (θair andθfc).

The limiting water contents, their coherence and parametrization

The limiting contents that comprise theRLL definition (Equation 1) are:

  • The field capacity limiting water content (θfc), defined byVeihmeyer and Hendrickson (1931)Veihmeyer FJ, Hendrickson AH. The moisture equivalent as a measure of the field capacity of soils. Soil Sci. 1931:32:181-193. as “the amount of water held in the soil after the excess gravitational water has drained away and after the rate of downward movement of water has materially decreased, which usually takes place within 2 or 3 days in pervious soils of uniform structure and texture”.

  • The limiting water content for proper aeration (θair), which is the maximum water content corresponding to an air-filled porosity that allows an adequate exchange of gases between soil and atmosphere, ensuring that in any rooted portion of the soil there is no excess of CO2 or lack of O2. The crop stress associated with this parameter (ifθ>θair) is the anoxic stress.

  • The limiting water content corresponding to the permanent wilting point (θpwp), which delimits soil conditions in which a plant will permanently (that is, without recovery) wilt in short time. The stress associated with this parameter (ifθ<θpwp) is the drought stress.

  • The limiting water content related to excess root penetration resistance (θpr), representing the water content below which the plant roots cannot penetrate the soil for mechanical restrictions. The stress associated with this parameter (ifθ<θpr) is the mechanical stress.

A cursory analysis of the limiting values that are used to defineRLL reveals some key aspects. Two of the limiting water contents,θpwp andθpr, represent conditions that correspond to the “fatal” limit of the respective stress - drought and mechanical (permanent wilting and zero root growth), while a third water content, θair, corresponds to the onset of the anoxic stress phase, when aeration conditions have just begun to cause damage to the plant, but below a fatal stress level. The three types of stress regarding their respective limiting contents are schematically represented in figure 2. In the example illustrated in this figure, there is a near-fatal stress condition (water content “1”, where the water content is slightly aboveθpwp), which, nevertheless, is within the range considered asRLL.

Figure 2
The Least Limiting Water Range (RLL) corresponding to case 1 of Figure 4c, with indication of drought, mechanical and anoxic stresses. The water content indicated with 1 represents a water content withinRLL but with near-fatal stress, and water content 2 refers to a water content outside the RLL, but without stress.

An opposite situation, i.e., water content outside theRLL, but without causing stress to the plant, may occur whenθfc<θair(cases 1 and 3 of Figure 1) and whenθfc<θ<θair(situation 2, Figure 2). In such conditions, the plant is not suffering any kind of stress, however, the corresponding water content is not included in theRLL. Considering the definition of field capacity, water contents higher thanθfc will rapidly decrease in time due to drainage, and θfcwill be reached in a matter of days. Depending on rainfall and other boundary conditions, the amount of water taken up by the plant in the range ofθfc andθair, whenθfc<θair,may, nevertheless, be considerable.

In virtually all studies published in Brazil on the subject, the parametrization of the four limiting contents of theRLL is fake and just consists of using a set of fixed values. Forθfc a water content is used that corresponds to a given matric potential (-10 or -33 kPa, in general);θair is chosen to match the air content of 0.10 m3 m-3;θpwp is considered to correspond to the matric potential of -1,500 kPa; andθpr is chosen equivalent to the water content at which the soil has a predefined resistance to penetration (often 2 MPa). Except in some studies (Klein and Camara, 2007Klein VA, Camara RK. Rendimento da soja e intervalo hídrico ótimo em Latossolo Vermelho sob plantio direto escarificado. R Bras Ci Solo. 2007;31:221-7.; Kaiser et al., 2009Kaiser DR, Reinert DJ, Reichert JM, Collares GL, Kunz M. Intervalo hídrico ótimo no perfil explorado pelas raízes de feijoeiro em um Latossolo sob diferentes níveis de compactação. R Bras Ci Solo. 2009;33:845-55.; Gubiani et al., 2013bGubiani PI, Goulart RZ, Reichert JM, Reinert DJ. Crescimento e produção de milho associados com o intervalo hídrico ótimo. R Bras Ci Solo. 2013b;37:1502-11.), the values in general are not tested and compared with biological measurements, and a sensitivity analysis of the outcomes of studies regarding these values is rarely performed. In some cases, only the effect of the value chosen for penetration resistance in θpr is assessed (Betioli Júnior et al., 2012Betioli Júnior E, Moreira WH, Tormena CA, Ferreira CJB, Silva AP, Giarola NFB. Intervalo hídrico ótimo e grau de compactação de um Latossolo Vermelho após 30 anos sob plantio direto. R Bras Ci Solo. 2012;36:971-82.; Moreira et al., 2014bMoreira WH, Tormena CA, Betioli Junior E, Petean LP, Alves SJ. Influência da altura de pastejo de azevém e aveia em atributos físicos de um Latossolo Vermelho distroférrico, após sete anos sob integração lavoura-pecuária. R Bras Ci Solo. 2014b;38:1315-26.). The justification for using one or another set of values, if any, has been limited to quoting other publications in which the same values were used. The scientific flaw of doing so becomes clear when one realizes that the cited papers did, on their turn, not present sustainable justification; the values employed were chosen by someone at a certain moment, but their validity is much more the product of self-assertion than experimental scientific confirmation.

Moreover, few studies focus on the association ofRLL and crop response. Those who do show low correlation, as found in a wide literature review (Gubiani et al., 2013aGubiani PI, Reichert JM, Reinert DJ. Indicadores hídrico-mecânicos de compactação do solo e crescimento de plantas. R Bras Ci Solo. 2013a;37:1-10.), as well as in an experimental corroboration with eight maize crop cycles (Gubiani et al., 2013bGubiani PI, Goulart RZ, Reichert JM, Reinert DJ. Crescimento e produção de milho associados com o intervalo hídrico ótimo. R Bras Ci Solo. 2013b;37:1502-11.).

Such apparent contradiction - it would be reasonable to assume that an indicator of soil physical quality applied to agronomy would have correlation with yield - may have two main explanations. There may be a structural problem in the adopted model; in other words, the parameters considered or the way they are related to calculate the RLL (Equation 1) are not sufficient to represent the processes that determine the crop yield. Alternatively, the calibration of the model may be wrong, i.e. the values attributed to the limiting contents are incorrect.

Taking into account that the parameters included in the calculation ofRLL can represent drought stress and anoxic stress, and also include mechanical root growth stress, the most likely explanation for the observed low correlation is the simple way by which the parameters are related inRLL as well as inaccuracies in their calibration.

Limiting contents versus physical properties

In the following, we describe some relationships between theRLL limiting water contents and primary soil physical properties that can be determined objectively.

The field capacity limiting water content versus hydraulic retention and conductivity

Field capacity, by definition, is related to the vertical water draining movement and, therefore, depends on the water flow through the soil profile. The water flow in porous media is described by the Darcy-Buckingham Law, which establishes that the water flux density is equal to the product of hydraulic conductivity and the total hydraulic gradient. This law, when combined with the mass conservation principle, results in the Richards Equation, allowing predictions on soil water conditions in time and space. In addition to information on the specific geometry of the problem, in order to apply the Richards Equation it is necessary to know the relationship between hydraulic conductivity (K), matric potential (h) and water content (θ). Such relationships are highly nonlinear and hysteretic, and describing them in a mathematical format allowing a relatively simple solution is still one of the challenges of soil physics.

The relation between soil water content or matric potential at field capacity and soil hydraulic properties may be assessed simulating an internal drainage experiment using a hydrologic model based on the Richards Equation. Following this approach, Twarakavi et al. (2009)Twarakavi NKC, Sakai M, Simunek J. An objective analysis of the dynamic nature of field capacity. Water Resour Res. 2009;45:1-9.used the Hydrus model (Simunek and van Genuchten, 2008Simunek J, van Genuchten MT. Modeling nonequilibrium flow and transport processes using HYDRUS. Vadose Zone J. 2008;7:782-97.) and estimated field capacity for a large number of soils, discussing the sensitivity of the determination to the soil hydraulic properties and the criterion of negligibility of the drainage rate to be adopted. These authors proposed the following empirical equation to estimate the water content at field capacity as a function of other soil physical parameters:

where θr,θs and nare parameters of the water retention equation by van Genuchten (1980)van Genuchten MT. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci Soc Am J. 1980;44, 892-8.;Ks is the saturated hydraulic conductivity and qfc is the flux considered negligible for the purposes of field capacity estimation. In this equation, Ks andqfc should be expressed in the same unit of time per length (TL-1).

By combining equation 3 with the van Genuchten (1980)van Genuchten MT. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci Soc Am J. 1980;44, 892-8. water retention equation, the following expression for the matric potential at field capacity is obtained:

In this equation, hfc assumes the inverse unit of the α parameter from the van Genuchten (1980)van Genuchten MT. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci Soc Am J. 1980;44, 892-8. equation. For example, ifα is chosen in cm-1,hfc will be in cm. Twarakavi et al. (2009)Twarakavi NKC, Sakai M, Simunek J. An objective analysis of the dynamic nature of field capacity. Water Resour Res. 2009;45:1-9. concluded that the use of a “negligible” flux of qfc = 0.01 cm d-1 resulted in good estimations of field capacity for a wide range of soils; and for this qfcvalue, equations 3 and 4, forKs in cm d-1, become, respectively:

Based on their studies, Twarakavi et al. (2009)Twarakavi NKC, Sakai M, Simunek J. An objective analysis of the dynamic nature of field capacity. Water Resour Res. 2009;45:1-9., similarly to Souza and Reichardt (1996)Souza LD, Reichardt K. Estimativas da capacidade de campo. R Bras Ci Solo. 1996;20:183-9., concluded that the use of a fixed value for the matric potential to estimate field capacity does not correspond to reality.Figure 3 shows thehfc values as a function of the saturated hydraulic conductivityKs, as calculated by equation 6 for soils with α = 0.015 cm-1 and for five values of parameter n. The figure illustrates the tendency that is intuitively perceived: higher Ksvalues result in quicker drainage, causing the soil to dry more (more negative matric potential) until reaching the draining rate that was pre-established as negligible.

Figure 3
Matric potential corresponding to field capacity (hfc) as a function of saturated hydraulic conductivity (Ks), according to Twarakavi et al. (2009) (Equation 6) considering a negligible bottom flux for the purpose of field capacity (qfc) of 0.01 cm d-1, for five van Genuchten (1980) soils (five values of shape factor n; shape factor α equal to 0.015 cm-1 in all cases).

The equation developed by Twarakavi et al. (2009)Twarakavi NKC, Sakai M, Simunek J. An objective analysis of the dynamic nature of field capacity. Water Resour Res. 2009;45:1-9. was included here just as an example. For its development a soil database was used that contains predominantly soils from temperate regions. Moreover, equation 4 does not allow the inclusion of vertical heterogeneity (stratification) of the soil, which, when it occurs, may have a major effect on the water content and on the matric potential corresponding to field capacity. It would not be possible to include such heterogeneity in a simple empirical equation like equation 4, because it would increase the number of parameters in such a way that a solution by regression, as obtained by Twarakavi et al. (2009)Twarakavi NKC, Sakai M, Simunek J. An objective analysis of the dynamic nature of field capacity. Water Resour Res. 2009;45:1-9., would become inviable. Moreover, applications to solve for the problem of multi-layer systems already exists, in the form of numerical hydrologic models based on the Richards Equation. One of these models (Hydrus) was used by Twarakavi et al. (2009)Twarakavi NKC, Sakai M, Simunek J. An objective analysis of the dynamic nature of field capacity. Water Resour Res. 2009;45:1-9. themselves and depends only on the hydraulic properties as a function of depth to calculate the flows.

The limiting water content corresponding to the permanent wilting pointversus retention and hydraulic conductivity

The water content corresponding to the permanent wilting point (θpwp), by definition, refers to the lower limit of extractable soil water by plant roots. Its experimental determination is time-consuming and obtained values are discussed, e.g., byHaise et al. (1955)Haise HR, Haas HJ, Jensen LR. Soil mixture studies of some great plain soils. II - Field capacity as related to 1/3-atmosphere percentage and ‘minimum point’ as related to 15- and 26-atmosphere percentage. Soil Sci Soc Am Proc. 1955;34:20-5., Veihmeyer and Hendrickson (1955)Veihmeyer FJ, Hendrickson AH. Does transpiration decrease as the soil moisture decreases? Trans Am Geophys Union. 1955;36:425-48. and Cutford et al. (1991)Cutford HW, Jefferson PG, Campbell CA. Lower limit of available water for three plant species grown on a medium-textured soil in southwestern Saskatchewan. Can J Soil Sci. 1991;71:247-52.. The soil water content can become lower than θpwponly by evaporation or drainage. Therefore, usingθpwp to calculate crop available water A seems plausible, but its use in theRLL definition, with the adjective least limiting, does not seem reasonable, as will be explained in the following.

The inclusion of θpwp inRLL as the lower limit corresponding to drought stress represents a gross simplification of the existing knowledge on transpiration reduction due to soil drying, as observed in experiments like those reported by Meyer and Green (1980)Meyer WS, Green GC. Water use by wheat and plant indicators of available soil water. Agron J. 1980;72:253-7., Wright and Smith (1983)Wright GC, Smith RGC. Differences between two grain sorghum genotypes in adaptation to drought stress. II - Root water uptake and water use. Aust J Agric Res. 1983;34:627-35., Casaroli et al. (2010)Casaroli D, De Jong van Lier Q, Dourado Neto D. Validation of a root water uptake model to estimate transpiration constraints. Agric Water Manage. 2010;97:1382-8.and many others. To improve the empirical modeling of this phenomenon, Doorenbos and Kassam (1986)Doorenbos J, Kassam AH. Yield response to water. Rome: FAO; 1986. (FAO Irrigation and Drainage Paper, 33)., Feddes et al. (1978Feddes RA, Kowalik PL, Zaradny H, editors. Simulation of field water use and crop yield. New York: John Wiley & Sons; 1978., 1988Feddes RA, Kabat P, Van Bakel PJT, Bronswijk JJB, Halbertsma J. Modelling soil water dynamics in the unsaturated zone - state of the art. J Hydrol. 1988;100:69-111.) and van Genuchten (1987)van Genuchten MT. A numerical model for water and solute movement in and below the root zone. Riverside: USDA, ARS; 1987. (U.S. Salinity Laboratory Research Report). developed mathematically simple proposals to describe the phenomenon, used in several hydrological models, such as Hydrusand SWAP. Lascano and van Bavel (1986)Lascano RJ, Van Bavel CHM. Root water uptake and soil water distribution: test of an availability concept. Soil Sci Soc Am J. 1986;48:233-7., Jarvis (1989Jarvis NJ. Simple empirical model of root water uptake. J Hydrol. 1989;107:57-72., 2011Jarvis NJ. Simple physics-based models of compensatory plant water uptake: Concepts and eco-hydrological consequences. Hydrol Earth Syst Sci. 2011;15:3431-46.) and Li et al. (2001)Li KY, De Jong R, Boisvert JB. An exponential root-water-uptake model with water stress compensation. J Hydrol. 2001;252:189-204. offered relatively simple proposals that allow considering the soil vertical heterogeneity as well.

The permanent wilting water contentθpwp is defined based on the process of water uptake by plant roots; to predict it, one should model the water flow towards the plant roots, which, in the same way as forθfc, can be done by means of the Richards equation with appropriate boundary conditions. Classical publications on the item are those by Gardner (1960)Gardner WR. Dynamic aspects of water availability to plants. Soil Sci. 1960;89:63-7. and Cowan (1965)Cowan IR. Transport of water in the soil-plant-atmosphere system. J Appl Ecol. 1965;2:221-39. who developed theories based on the Richards equation and basic theory of water extraction from the soil by the plants.

Numerical algorithms based on the theory by Gardner (1960)Gardner WR. Dynamic aspects of water availability to plants. Soil Sci. 1960;89:63-7. were presented in several recent publications (Javaux et al., 2008Javaux M, Schröder T, Vanderborght J, Vereecken H. Use of a three-dimensional detailed modeling approach for predicting root water uptake. Vadose Zone J. 2008;7:1079-88., 2013Javaux M, Couvreur V, Vanderborght J, Vereecken H. Root water uptake: from 3-D biophysical processes to macroscopic modeling approaches, Vadose Zone J. 2013. doi:10.2136/vzj2013.02.0042, 2013; De Jong van Lier et al., 2006De Jong van Lier Q, Metselaar K, Van Dam JC. Root water extraction and limiting soil hydraulic conditions estimated by numerical simulation. Vadose Zone J. 2006;5:1264-77., 2013De Jong van Lier Q, Van Dam JC, Durigon A, Santos MA, Metselaar K. Modeling water potentials and flows in the soil-plant system comparing hydraulic resistances and transpiration reduction functions. Vadose Zone J. 2013;12:1-20.; Couvreur et al., 2012Couvreur V, Vanderborght J, Javaux M. A simple three-dimensional macroscopic root water uptake model based on the hydraulic architecture approach. Hydrol Earth Syst Sci. 2012;16:2957-71.). The plant root geometry and properties, specifically the limiting root (or leaf) water potential and internal plant resistances, together with the soil water potential and the hydraulic conductivity of the soil near the roots, determine whether a plant can or cannot withdraw water at rates compatible with the atmospheric demand, allowing the assessment of drought stress.

Consequently, regarding the lower limit for crop water availability, retention and hydraulic conductivity are the physical properties of major interest. De Jong van Lier et al. (2006De Jong van Lier Q, Metselaar K, Van Dam JC. Root water extraction and limiting soil hydraulic conditions estimated by numerical simulation. Vadose Zone J. 2006;5:1264-77., 2008De Jong van Lier Q, Van Dam JC, Metselaar K, De Jong R, Duijnisveld WHM. Macroscopic root water uptake distribution using a matric flux potential approach. Vadose Zone J. 2008;7:1065-78.) showed that the matric flux potential (M, m2 d-1) (Gardner, 1958Gardner WR. Some steady-state solutions of the unsaturated moisture flow equation with application to evaporation from a water table. Soil Sci. 1958;85:228-32.; Raats, 1970Raats PAC. Steady infiltration from line sources and furrows. Soil Sci Soc Am Proc. 1970;34:709-14.; Pullan, 1990Pullan AJ. The quasilinear approach for unsaturated porous media flow. Water Resour Res. 1990;26:1219-34.), a soil physical property that integrates retention and conductivity, is the physical parameter most directly linked to crop available water. M is defined as:

where hl is the limiting matric potential, the most negative value that the water in the roots can assume. By writing the Darcy equation in terms of M it is shown that the spatial gradient of M (dM/dx) is equal to the water flux density in the absence of gravitational flow, which, approximately, is the case for water uptake from the soil by plant roots.

When the K(h) function is known, the determination ofM(h) is a matter of analytical or numerical integration.Figure 4 illustrates, for two differently textured soils, hydraulic conductivity K (Figure 4a) and its integral, the matric flux potential M (Figure 4b) as calculated by equation 7, as a function of the matric potential. For a given matric potential, the M value corresponds to the ease with which a plant withdraws water from the soil. Figure 4c shows the ratio of the M values of the two soils as a function of the matric potential. It can be seen that for not too negative matric potentials (wet soil), the relative difference between theM values for both soils is small, but for more negative matric potentials (drier soil) the M of the clay loam soil becomes orders of magnitude higher than in the sandy loam soil. This indicates a much higher availability of water in the clay loam soil.

Figure 4
Hydraulic conductivity K (a), the matric flux potential M (b) (integral as calculated by equation 7, with hl = -15,000 cm) as a function of the matric potential for two different soil textures and the ratio of the M values of both soils (c) as a function of the matric potential.

The limiting water content for proper aeration versus gas diffusion

Soil aeration is the process through which the gases produced or consumed in soil are exchanged with the atmosphere, primarily by diffusion. The main gas consumed in the soil is oxygen, whereas the main gas produced is carbon dioxide. Thus, soil aerobic processes depend on the flow of O2 from the atmosphere into the soil, and of CO2 in the opposite direction. Diffusion of these gases is much higher in air than in water. So, the higher the soil water content the greater the difficulty for an adequate aeration, which can hinder plant growth and should be considered for the modeling of plant production.

By integrating various factors that control the aeration process in order to meet the biological demand, as described by De Jong van Lier (2001)De Jong van Lier Q. Oxigenação do sistema radicular: uma abordagem física. R Bras Ci Solo. 2001;25:233-8., De Jong van Lier and Cichota (2004)De Jong van Lier Q, Cichota R. Modelagem da aerabilidade do solo. In: Anais da 15ª Reunião Brasileira de Manejo e Conservação do Solo e da Água; 2001; Santa Maria. Santa Maria: Sociedade Brasileira de Ciência do Solo; 2004. v.1. presented the following equation to estimate the minimum required aeration porosity as a function of the O2consumption in the soil, the depth and total porosity of soil, the O2content in the atmosphere and O2 diffusion in the air:

In this equation, βmin(m3 m-3) is the minimum aeration porosity necessary to ensure that no part of the plant root system will suffer from lack of oxygen;Smax (kg m-3s-1) is the oxygen consumption rate per unit of soil volume at the soil surface (Smax was considered as 28.2∙10-8 kg m-3 s-1);Z(m) is the rooting depth; χ(m3 m-3) is total porosity (χ was considered equal to 0.5 m3 m-3); [O2]atm (kg m-3), the oxygen concentration in atmosphere, considered as 0.269 kg m-3; [O2]min (kg m-3) is the minimum oxygen concentration to sustain life at a given soil depth (considered [O2]min equal to [O2]atm/4);Da (m2s-1) is air diffusivity for oxygen (Da = 1.78∙10-5m2 s-1); and p is a factor determining the shape of the function of O2 consumption decrease with depth (SO2, kg m-3 s-1) (De Jong van Lier and Cichota, 2004De Jong van Lier Q, Cichota R. Modelagem da aerabilidade do solo. In: Anais da 15ª Reunião Brasileira de Manejo e Conservação do Solo e da Água; 2001; Santa Maria. Santa Maria: Sociedade Brasileira de Ciência do Solo; 2004. v.1.), whereSO2 is equal toSmax at the surface (z = 0) and zero when z =Z, according to:

Equation 8 was applied for rooting depths up to 1.0 m and forp-values of equation 9 that represent different oxygen consumption distributions with depth, from consumption concentrated at the soil surface (p=0.1) to a practically equal consumption at all depths occupied by the root system (p=10) (Figure 5). As a result,βmin ranged from 0.02 to 0.24 m3 m-3. This range could be enlarged considering other values for χ andSmax. It is worth noting that the βmin values, as calculated by equation 8 and represented in figure 5refer to the minimum air content at the soil surface, through which all O2 consumed and CO2 produced in the soil profile is flowing. The greater the depth, the smaller the flow of O2 and CO2 and the smaller will be the content of air required for an adequate aeration. Clearly, the application of any fixed value forβmin, such as the value of 0.1 m3 m-3 used for calculation ofθair ofRLL,means failing to consider the complexity of the aeration process and believing that it is possible to infer biological consequences from a fixed fraction of soil porosity not occupied by water.

Figure 5
Minimum air-filled porosity for proper aerationβmin as a function of the empirical parameter p and depth of the root system Z, as calculated by equation 8.

A more complete approach of the subject is presented by Bartholomeus et al. (2008)Bartholomeus RP, Witte JM, Van Bodegom PM, Van Dam JC, Aerts R. Critical soil conditions for oxygen stress to plant roots: Substituting the Feddes-function by a process-based model. J Hydrol. 2008;360:147-65., who modeled, besides the soil macroscopic diffusion, the microscopic oxygenation process of individual roots, and found βmin values varying from 0.02 and 0.08 m3 m-3 for diverse scenarios, together with a high sensitivity for soil temperature and rooting depth. According to these authors, fixed values for βmin(such as the 0.1 m3 m-3 value adopted in virtually all works involving RLL) very unlikely will represent an approximate value of the actual value.

The limiting water content related to penetration resistanceversus the physics of root elongation

The resistance that solid particles in the soil apply against a motion or displacement results from complex relations of cohesive forces of solid-solid interactions, adhesive forces of solid-liquid interactions and frictional forces generated during displacement (Keller et al., 2013Keller T, Lamandé M, Peth S, Berli M, Delenne J, Baumgarten W, Rabbel W, Radjaï F, Rajchenbach J, Selvadurai A, Or D. An interdisciplinary approach towards improved understanding of soil deformation during compaction. Soil Till Res. 2013;128:61-80.). When a root penetrates the pores of a soil, its growth will cause a significant displacement of the particles when it does not change its geometry to fit the irregular geometry of the pore space, when it forces its way through the pore cavity that it had filled or when pore space must be created (Bengough et al., 1997Bengough AG, Croser C, Pritchard J. A biophysical analysis of root growth under mechanical stress. Plant Soil. 1997;189:155-64., 2006Bengough AG, Bransby MF, Hans J, Mckenna SJ, Roberts TJ, Valentine TA. Root responses to soil physical conditions; growth dynamics from field to cell. J Exp Bot. 2006;57:437-47., 2011Bengough AG, Mckenzie BM, Hallett PD, Valentine TA. Root elongation, water stress, and mechanical impedance: a review of limiting stresses and beneficial root tip traits. J Exp Bot. 2011;62:59-68.; Clark et al., 2003Clark LJ, Whalley WR, Barraclough PB. How do roots penetrate strong soil? Plant Soil. 2003;255:93-104.). Root growth can be considered a cell mass flow into the soil (cell division and elongation). Directing the cell flow to less resistant areas is a biological strategy to save the energy captured in the photosynthesis process.

The Lockhart’s model (Lockhart, 1965Betioli Júnior E, Moreira WH, Tormena CA, Ferreira CJB, Silva AP, Giarola NFB. Intervalo hídrico ótimo e grau de compactação de um Latossolo Vermelho após 30 anos sob plantio direto. R Bras Ci Solo. 2012;36:971-82.; Jordan and Dumais, 2010Jordan BM, Dumais J. Biomechanics of plant cell growth. In: Encyclopedia of life sciences (ELS). Chichester: John Wiley & Sons; 2010.), employed in plant physiology to study the force relations involved in cell elongation, considers that cell walls behave as a Bingham fluid. A Bingham fluid, as opposed to a Newtonian fluid, requires a minimum pressureσY>0 to cause its deformation. The relative deformation rate over time will be proportional to the applied pressure σ (Pa) reduced byσY (Pa):

where L (m) represents the root length, t (s) is time, and the property μ (Pa-1 s-1) is called extensibility. In the case of cell elongation,σY can be understood as the cell wall resistance (R, Pa), i.e. the pressure at which the cell walls resist deformation. The cell turgor (T, Pa) is the internal cell pressure and is equivalent to the σ of equation 10 (Bengough et al., 1997Bengough AG, Croser C, Pritchard J. A biophysical analysis of root growth under mechanical stress. Plant Soil. 1997;189:155-64.). The difference between T and R is called net cell pressure (Plc, Pa):

By comparing equations 10 and 11, we see thatPlc represents the rate of deformation over time divided by its extensibility.

The soil mechanical resistance to root growth effectively occurs only when turgid roots exert pressure on the surrounding environment. In these cases a reaction pressure (Prs) is produced by the soil, its value being limited to the soil mechanical impedancePrs,max. The direction of the resulting force from such pressure in the cell wall is opposite to the force provided by the cell turgor pressure (Plc). The roots continue growing and moving soil particles whilePlc>Prs,max. Note that the presence of roots is indispensable for the existence ofPrs, but their presence alone is not sufficient, because those roots withPlc=0 do not provide the required physical conditions for the occurrence ofPrs.

Quantification of the soil mechanical impedance can be achieved by using available types of penetrometers and measurement techniques, as discussed byStolf (1991)Stolf R. Teoria e teste experimental de fórmula de transformação dos dados de penetrômetro de impacto em resistência do solo. R Bras Ci Solo. 1991;5:229-35.. In general, the soil mechanical impedance measured by penetrometers (RP, Pa) is the quantification of the soil response to pressure, resulting from the reaction to the pressure caused by the penetration of a metal cone positioned at the end of a rod into the soil. Although RP is useful to detect differences in soil resistance as a function of cohesive, adhesive or frictional forces and to distinguish pedological and tillage factors that affect such forces, it will not represent the manifestation ofPrs, given the physical conditions required for the occurrence ofPrs and because of the different scales and geometries of the cell expansion when compared to the penetrometer cone.

Because of the difficulty in determiningPrs,max, some researchers attempted to establish empirical relationships betweenRP and the cell elongation rate dL/dt. In a recent literature review, some relationships were presented and discussed by Bengough et al. (2011)Bengough AG, Mckenzie BM, Hallett PD, Valentine TA. Root elongation, water stress, and mechanical impedance: a review of limiting stresses and beneficial root tip traits. J Exp Bot. 2011;62:59-68.. From their experiments it became evident that in order to obtain empirical relations a homogeneous soil matrix (samples prepared in laboratory) was a prerequisite, as well as matric potentials less negative than -0.5 to -0.1 MPa. In the same moist conditions, the expectation is that the relation between dL/dt andRP will not be confirmed in cases where cracks and biopores are present, occupying a large portion of the soil volume (the plant directs the cell flow to less resistant spaces) which is common in conservation-tillage crop areas (Ehlers et al., 1983Ehlers W, Köpke U, Hesse F, Böhm W. Penetration resistance and growth root of oats in tilled and untilled loess soil. Soil Till Res. 1983;3:261-75.). For such areas and based on the criterion that there has not been a close relationship between some biological variable withRP, some researchers suggest that critical RP values should be higher than 2 MPa (Klein and Camara, 2007Klein VA, Camara RK. Rendimento da soja e intervalo hídrico ótimo em Latossolo Vermelho sob plantio direto escarificado. R Bras Ci Solo. 2007;31:221-7.; Reichert et al., 2009Reichert JM, Suzuki LEAS, Reinert DJ, Horn R, Hakansson I. Reference bulk density and critical degree-of-compactness for no-till crop production in subtropical highly weathered soils. Soil Till Res. 2009;102:242-54.; Moraes et al., 2014Moraes MT, Debiasi H, Carlesso R, Franchini JC, Silva VR. Critical limits of soil penetration resistance in a rhodic Eutrudox. R Bras Ci Solo. 2014;38:288-98.). However, the use of values above 2 MPa (e.g. 3.5 MPa as suggested by Moraes et al., 2014Moraes MT, Debiasi H, Carlesso R, Franchini JC, Silva VR. Critical limits of soil penetration resistance in a rhodic Eutrudox. R Bras Ci Solo. 2014;38:288-98.) often corresponds to matric potentials close to the permanent wilting point. In some cases of the study conducted by Moraes et al. (2014)Moraes MT, Debiasi H, Carlesso R, Franchini JC, Silva VR. Critical limits of soil penetration resistance in a rhodic Eutrudox. R Bras Ci Solo. 2014;38:288-98.,RP would reach 3.5 MPa only for matric potentials below the wilting point. In such conditions,Plc and, consequently,Prs, can be so small that the biological stress associated with the mechanical factor is negligible if compared to the stress associated with the drought factor. In these cases, interpreting the occurrence of loss in biological variables as caused by a highRP means the incorrect attribution of the cause to the mechanical factor when it is essentially a lack of water.

The empirical relationships established betweenRP and other soil properties, such as those used to estimateRP as a function of bulk density (ρ) and water content [RP = f (ρ, θ)], and which constitute the mathematical sub-model of the limiting mechanical impedance often used in combination with the calculation ofRLL, are useful to represent the relationship between the soil properties included in the sub-model. However, such models of relationships between soil properties do not contain any information about the occurrence and magnitude ofPrs. Its use in the context of RLL aiming to estimateθpr corresponding to a fixedRP only informs that the mechanical reaction of soil to a metal rod, in terms of pressure, would be equal to the fixed RP value. Assuming that the roots face this same resistance means accepting all interpretation errors of the cause factor and the conceptual errors previously mentioned.

Problems and mistakes in interpreting theRLL

Interpretation of RLL values is achieved by comparison. It is often concluded that the soil physical quality worsened when the RLL decreased as a function of some agricultural tillage system or other factor. Such comparison usually includes an assessment of different scenarios. Changes in soil management, for example, represent changes in the plant root density, rooting depth, and soil hydraulic properties. At the same time, theRLL in its original conception has limitations when the scenario presents more complex boundary conditions. Is it plausible to use fixed values for the limiting contents regardless of the phenological and pedological boundary conditions? The previous sections raised several arguments for a negative answer. So, would the calculated RLL that disregards a change in the boundary conditions be valuable as a comparative indicator? Probably not, as previously discussed.

Then, how could the RLL be used as an indicator? According to the previous discussion, it would be necessary to estimate the limiting contents considering the changing boundary conditions. This process would be as complex, or even more complex, than the use of a detailed mechanistic model and, therefore, the advantage ofRLL being a “simple” indicator would disappear.

Regarding the upper limit of theRLL, the inconsistencies in the choice of θair orθfc as limiting values have already been discussed. For the lower limit, it was observed thatθpwp corresponds to a condition of total drought stress, which would better be replaced by a value corresponding to a considerable stress, but not total, or, in analogy to what is done in the case of θair, to the water content corresponding to the onset of drought stress. The relationship between drought stress and matric potential is discussed in detail by Metselaar and De Jong van Lier (2007)Metselaar K, De Jong van Lier Q. The shape of the transpiration reduction function under plant water stress. Vadose Zone J. 2007;6:124-39. and by De Jong van Lier et al. (2009)De Jong van Lier Q, Dourado Neto D, Metselaar K. Modeling of transpiration reduction in Van Genuchten–Mualem type soils. Water Res Res. 2009;45:W02422. doi:10.1029/2008WR006938. The theory presented in these publications shows that drought stress, quantified by relative transpiration, increases in a nonlinear way when the matric potential decreases from a certain value, generally in the range of -30 to -100 kPa (Kukal et al., 2005Kukal SS, Hira GS, Sidhu AS. Soil matric potential-based irrigation scheduling to rice (Oryza sativa). Irrig Sci. 2005;23:153-9.; Liu et al., 2012Liu H, Yang H, Zheng J, Jia D, Wang J, Li Y, Huang G. Irrigation scheduling strategies based on soil matric potential on yield and fruit quality of mulched-drip irrigated chili pepper in Northwest China. Agric Water Manage. 2012;115:232-41.), much less negative than the permanent wilting point. Thus, the use of a water content corresponding to this order of magnitude of potential (critical water content,θcr) to replaceθpwp in the calculation of theRLL sounds reasonable and was proposed yet recently by Silva et al. (2015)Silva BM, Oliveira GC, Serafim ME, Silva EA, Ferreira MM, Norton LD, Curi N. Critical soil moisture range for a coffee crop in an oxidic Latosol as affected by soil management. Soil Till Res. 2015:154:103-13.. Their data show that, as a consequence, the permanence ofθpr in theRLL calculation would be superfluous, once its value, in almost all situations, would be lower thanθcr. This finding corroborates with the fact that the major models of plant growth do not use, among their input data, information on the penetration resistance, and yet they result in good predictions. Apparently, the process (plant growth) is not sensitive to parameters that relate penetration resistance to water content.

Two boundary conditions, common in numerous production systems, will be discussed in the following sections. The first refers to the phenology of the crop under observation; the second, to variations in soil properties with depth. For neither of them RLL can be easily introduced as an indicator.

Phenological boundary conditions and their consequences for theRLL

During crop development, plants undergo a series of phenological changes that affect the transfer of mass and energy in the soil-plant-atmosphere system. The most obvious parameters include changes in the rooting depth, distribution and growth, leaf area index, carbohydrate partitioning among plant organs and other physiological parameters. Hydrological models, when coupled to plants growth models, usually treat these parameters as a function of time, development stage or accumulated heat sum. For the purpose ofRLL, interpretation would be achieved based on the limiting contents that, except forθfc, which is exclusively related to soil, should be considered as variable during crop development.

As the rooting depth and root length density increase,βmin increases and the corresponding θair decreases. Thus, for the first phenological stage, the required air-filled porosity is very low and θair becomes virtually equal toθs; at the peak of root growth, according to the theory presented,βmin can represent half of the total porosity (Figure 5). The great variation of these values during the cycle of an annual crop shows that the use of a RLL with fixed limiting values is not consistent with the real conditions in this context either.

Regarding θpr, in the reproductive stage root growth becomes irrelevant. Hence, in this phenological stage, the mechanical resistance is no longer biologically important,θpr equals zero, and the expression to calculate RLL(Equation 1) could be reduced to:

This fact is important, for example, when assessing the correlation between grain productivity and RLL, which may show a significant reduction of RLLif θpr is the lower limit ofRLL andθpwp<θ<θpr, but this reduction is not correlated to dry matter accumulation during the reproductive phase.

The value of θpwp (or the critical water content θcr that might replace it) should also be dependent on crop phenology. As previously discussed,θpwp is influenced by the root system's ability to withdraw a sufficient amount of water from the soil to meet the transpiration demand. The transpiration demand is determined, among other factors, by the leaf area index, which is dependent on development stage. Simultaneous with an increase in leaf area, the root system also grows in depth and density. The combined effect of these factors should be considered for an accurate determination of the value ofθpwp orθcr.

Pedological boundary conditions and their consequences for theRLL

Most soils present morphological stratification. The distribution of the root system in the soil, with a density that normally decreases with depth, superposes this pedological stratification. Consequently, the rate at which plant roots take up water from the soil is not uniform with depth, and as soil water redistribution is a slow process, water content variations between depths will arise. Thus, one can expect that the values of the several limiting water contents, regardless of how they are estimated, vary with depth. Taking theRLL indicator, its definition is applicable to only one soil horizon, however, scenarios in which the conditions of one horizon are outside theRLL, while another horizon is within RLL may be supposed to occur frequently. In these cases, the indicator does not allow a unique prediction. To what extent would the favorable conditions in one horizon compensate the unfavorable conditions in another one?

Considering available water A, the integration of its value along the rooted soil depth ze (m) results in the available water capacityCAW (m):

A similar integration could be proposed for theRLL, but it would not make much sense once the purpose of RLLis its use as an indicator of soil quality. A simple integration, as in the case of CAW, would disregard any compensatory effect that may exist when one soil horizon is outside the optimum range and another one is within it.

Critical density

Very common in publications addressing the least limiting water range is its estimation for several soil densities or degrees of compaction, aiming to predict the effects of a possible compaction on the soil physical quality and crop yields (Imhoff et al., 2001Imhoff S, Silva AP, Dias Junior MS, Tormena CA. Quantificação de pressões críticas para o crescimento das plantas. R Bras Ci Solo. 2001;25:11-8.; Leão et al., 2004Leão TP, Silva AP, Macedo MCM, Imhoff S, Euclides VPB. Intervalo hídrico ótimo na avaliação de sistemas de pastejo contínuo e rotacionado. R Bras Ci Solo. 2004;28:415-23.; Beutler et al., 2008Beutler AN, Centurion JF, Silva AP, Centurion MAPC, Leonel CL, Freddi OS. Soil compaction by machine traffic and least limiting water range related to soybean yield. Pesq Agropec Bras. 2008;43:1591-600.; Blainski et al., 2009Blainski E, Gonçalves ACA; Tormena CA, Folegatti MV, Guimarães RML. Intervalo hídrico ótimo num Nitossolo Vermelho distroférrico irrigado. R Bras Ci Solo. 2009;33:273-81.; Kaiser et al., 2009Kaiser DR, Reinert DJ, Reichert JM, Collares GL, Kunz M. Intervalo hídrico ótimo no perfil explorado pelas raízes de feijoeiro em um Latossolo sob diferentes níveis de compactação. R Bras Ci Solo. 2009;33:845-55.; Petean et al., 2010Petean LP, Tormena CA, Alves SJ. Intervalo hídrico ótimo de um Latossolo Vermelho distroférrico sob plantio direto em sistema de integração lavoura-pecuária. R Bras Ci Solo. 2010;34:1515-26.; De Lima et al., 2012De Lima CLR, Miola ECC, Timm LC, Pauletto EA, Silva AP. Soil compressibility and least limiting water range of a constructed soil under cover crops after coal mining in Southern Brazil. Soil Till Res. 2012;124:190-5.; Farias et al., 2013Farias IL, Pacheco EP, Viégas PRA. Characterisation of the optimal hydric interval for a Yellow Argisol cultivated with sugarcane on the coastal plains of Alagoas, Brazil. R Ci Agron. 2013;44:669-75.; Guimarães et al., 2013Guimarães RML, Tormena CA, Blainski E, Fidalski J. Intervalo hídrico ótimo para avaliação da degradação física do solo. R Bras Ci Solo. 2013;37:1512-21.; Moreira et al., 2014aMoreira FR, Dechen SCF, Silva AP, Figueiredo GC, Maria IC, Pessoni PT. Intervalo hídrico ótimo em um Latossolo Vermelho cultivado em sistema semeadura direta por 25 anos. R Bras Ci Solo. 2014a;38:118-27.; Seben Junior et al., 2014Seben Junior GF, Corá JE, Lal R. Effect of cropping systems in no-till farming on the quality of a Brazilian Oxisol. R Bras Ci Solo. 2014;38:1268-80.). In this case,θfc andθpwp are usually estimated as a function of an empirical relationship between the water content, bulk density and matric potential proposed by Silva et al. (1994)Silva AP, Kay BD, Perfect E. Characterization of the least limiting water range. Soil Sci Soc Am J. 1994;58:1775-81. or other also empirical equations (Guimarães et al., 2013Guimarães RML, Tormena CA, Blainski E, Fidalski J. Intervalo hídrico ótimo para avaliação da degradação física do solo. R Bras Ci Solo. 2013;37:1512-21.), using fixed values for the matric potential; θpr is calculated as a function of soil bulk density using the empirical model proposed by Busscher (1990)Busscher WJ. Adjustment of flat-tipped penetrometer resistance data to a common water content. Trans Am Soc Agron Eng. 1990;33:519-24., assessing it for a fixed impedance value, commonly 2 MPa; andθair is calculated by subtracting 0.1 m3 m-3 from total porosity, which is a function of the bulk and particle density.

The results are usually presented in the form of a graph of the four limiting water content values, where the abscissa represents the soil bulk density and the ordinate represents the water content, as exemplified in figure 6. In this context, “critical bulk density” is defined as the lowest soil density for which theRLL equals zero. For such bulk density, or higher values, irrespective of the water content in the soil, it would be outside the RLL, and, therefore, productivities in such soils would always be reduced.

Figure 6
Representation ofRLL limiting water contents versus bulk density. Solid lines represent the limiting water contents ofRLL estimated based on fixed parameters; dotted lines show the expected behavior of the limiting water contents based on physical knowledge. The gray colored area represents theRLL integral over the range of represented densities.θpwpθfc,θair,θrp: water water contents corresponding to permanent wilting, field capacity, limiting air-filled porosity and limiting penetration resistance, ρ: bulk density, ρs: particle density.

The described procedure is emblematic for the shift of attention that resulted from the introduction of the RLL in Brazilian scientific circles. The processes and mechanisms that govern the system are not discussed. Instead, conclusions are based on the perfection attributed to the indicator. To its value, obtained by extrapolation of experimental observations using empirical relationships, predictive qualities and recommendations on soil management to establish required bulk densities are attributed. However, it is clear that, in case of soil compaction (a decrease in total porosity, particularly in macroporosity), there are predictable changes for almost all soil physical properties. Such changes will influence theRLL limiting contents in an equally predictable manner, as exemplified by the dotted lines in figure 6 and detailed as follows:

  1. The saturated or near-saturated hydraulic conductivity decreases as macroporosity decreases. Hence, the draining process slows down, and the water content at field capacity (dependent on the negligibility criterion) will correspond to a less negative matric potential (wetter soil) than for the same non-compacted soil, as also discussed in the context of the equation 6 developed by Twarakavi et al. (2009)Twarakavi NKC, Sakai M, Simunek J. An objective analysis of the dynamic nature of field capacity. Water Resour Res. 2009;45:1-9..

  2. By reducing total porosity, according to De Jong van Lier (2001)De Jong van Lier Q. Oxigenação do sistema radicular: uma abordagem física. R Bras Ci Solo. 2001;25:233-8. and Bartholomeus et al. (2008)Bartholomeus RP, Witte JM, Van Bodegom PM, Van Dam JC, Aerts R. Critical soil conditions for oxygen stress to plant roots: Substituting the Feddes-function by a process-based model. J Hydrol. 2008;360:147-65., the minimum air-filled porosity is expected to decrease. Taking into account that compaction is likely to cause a reduction of the root growth as well as of the depth of the soil explored by the root system, the minimum air-filled porosity would be even more reduced. For example, using the values presented by De Jong van Lier (2010)De Jong van Lier Q. Disponibilidade de água às plantas. In: De Jong van Lier Q, editor. Física do solo. Viçosa, MG: Sociedade Brasileira de Ciência do Solo; 2010. p.283-98. for high oxygen consumption conditions, if total porosity decreases from 0.6 to 0.5 m3 m-3 and the depth of the root system decreases from 0.5 to 0.4 m, the minimum air-filled porosity would decrease from 0.17 to 0.13 m3 m-3.

  3. Compaction reduces the soil hydraulic conductivity at high water contents, but increases the hydraulic conductivity under drier conditions (Hillel, 2003Hillel D. Introduction to environmental soil physics. London: Academic Press; 2003.). Thus, in drier soils, access to water will be easier in compressed soil, and the matric potential corresponding toθpwp orθcr is expected to decrease (become more negative).

These predictable facts show the non-validity of fixed values of air content and matric potentials for calculation of theRLL limiting contents when the objective is to compare its value before and after a compaction. AsRLL assessment is always made with fixed values to assess the limiting levels, whereas the existing knowledge indicates that such levels should change because of compaction, the reliability of the critical densities found by this estimation method is, at least, questionable.

FINAL CONSIDERATIONS

Considering the availability of robust models for agronomy, ecology, hydrology, meteorology and other related areas, the appeal of theRLL as an indicator in Brazilian soil physics cannot be assumed for its effectiveness, never proven, but should be explained by the simplicity with which it is dealt. The exercise of estimating the limiting water contents in the light of the state of the art of soil physics, instead of considering them a function of simple constants, proves to be as or more complex than the existing mechanistic models. The amount and extent of objective information contained in theRLL limiting contents, when based on constants, is minimal when compared to the amount of objective information contained in continuous functions of properties of biological and physical processes. Thus, the great disadvantage of using theRLL as an indicator compared to the use of models that integrate such functions based on mechanistic knowledge becomes evident. The RLL is a mathematical model, but a sensitivity analysis ofRLL to the limiting values was only performed in one of the first publications on the subject (Silva et al., 1994Silva AP, Kay BD, Perfect E. Characterization of the least limiting water range. Soil Sci Soc Am J. 1994;58:1775-81.). Since then, what has been done as an approximation of calibration is not much more than verifying the dependence between RLL and the soil bulk density. The few tests of the correlation betweenRLL and crop yields, as described by Gubiani et al. (2013a)Gubiani PI, Reichert JM, Reinert DJ. Indicadores hídrico-mecânicos de compactação do solo e crescimento de plantas. R Bras Ci Solo. 2013a;37:1-10., in most cases refuted it, as also reported by Gubiani et al. (2013b)Gubiani PI, Goulart RZ, Reichert JM, Reinert DJ. Crescimento e produção de milho associados com o intervalo hídrico ótimo. R Bras Ci Solo. 2013b;37:1502-11.. Further research conducted with the purpose of revealing the relationship between RLL and soil bulk density or compaction, without questioning the limiting values or testing the sensitivity of RLL to these values will not add knowledge and should be discouraged. Research onRLL can provide scientific contribution, but should inform in which circumstanceRLL is a reliable indicator of the soil physical quality for plant production by comparing plant variables with water content in the context ofRLL.

CONCLUSIONS

The four limiting contents of the RLLare related to crop yield. However, their determination using fixed values of air content, matric potentials and mechanical impedance of soil marginalizes the existing knowledge in the field of soil physics, plant physiology and agrometeorology.

When calculated using the limiting water content for the onset of drought stress instead of the permanent wilting point,RLL shows insensitive to parameters that relate penetration resistance to water content. This corroborates to the fact that major crop growth and yield prediction models do not include soil mechanical information among their input parameters, and suggests that the determination of these parameters in the context of soil water availability and crop growth is counterproductive.

The search for quantifying concepts such as the least limiting water rangeRLL goes against the systemic understanding, being restricted to the simple determination of the respective limiting values, relegating the study or concern with the actual functioning of the system to a secondary role. Contrary to more complex models, indicators such as theRLL, with fixed limiting values as most researchers use them, show no correlation with crop productivity.

In the light of such findings, the attention given to theRLL in Brazilian soil physics is not justified. The authors of this study suggest a realignment of the research in Brazilian soil physics with the scientific principles towards emphasizing mechanistic soil physics rather than the search for empirical correlations far below the state of the art of soil physics.

REFERENCES

  • Abdelhafeez AT, Harssema H, Verkerk K. Effects of air temperature, soil temperature and soil moisture on growth and development of tomato itself and grafted on its own and egg-plant rootstock. Sci Hortic. 1975;3:65-73.
  • Allen RG, Pereira LS, Raes DE, Smith M. Crop evapotranspiration: Guidelines for computing crop water requirements. Rome: FAO; 1998. (Irrigation and Drainage Paper, 56).
  • Arnold JG, Moriasi DN, Gassman PW, Abbaspour KC, White MJ, Srinivasan R, Santhi C, Harmel RD, Van Griensven A, Van Liew MW, Kannan N, Jha MK. SWAT: Model use, calibration, and validation. Trans Am Soc Agric Biol Eng. 2012;55:1491-508.
  • Bartholomeus RP, Witte JM, Van Bodegom PM, Van Dam JC, Aerts R. Critical soil conditions for oxygen stress to plant roots: Substituting the Feddes-function by a process-based model. J Hydrol. 2008;360:147-65.
  • Bengough AG, Bransby MF, Hans J, Mckenna SJ, Roberts TJ, Valentine TA. Root responses to soil physical conditions; growth dynamics from field to cell. J Exp Bot. 2006;57:437-47.
  • Bengough AG, Croser C, Pritchard J. A biophysical analysis of root growth under mechanical stress. Plant Soil. 1997;189:155-64.
  • Bengough AG, Mckenzie BM, Hallett PD, Valentine TA. Root elongation, water stress, and mechanical impedance: a review of limiting stresses and beneficial root tip traits. J Exp Bot. 2011;62:59-68.
  • Betioli Júnior E, Moreira WH, Tormena CA, Ferreira CJB, Silva AP, Giarola NFB. Intervalo hídrico ótimo e grau de compactação de um Latossolo Vermelho após 30 anos sob plantio direto. R Bras Ci Solo. 2012;36:971-82.
  • Beutler AN, Centurion JF, Silva AP, Centurion MAPC, Leonel CL, Freddi OS. Soil compaction by machine traffic and least limiting water range related to soybean yield. Pesq Agropec Bras. 2008;43:1591-600.
  • Blainski E, Gonçalves ACA; Tormena CA, Folegatti MV, Guimarães RML. Intervalo hídrico ótimo num Nitossolo Vermelho distroférrico irrigado. R Bras Ci Solo. 2009;33:273-81.
  • Busscher WJ. Adjustment of flat-tipped penetrometer resistance data to a common water content. Trans Am Soc Agron Eng. 1990;33:519-24.
  • Casaroli D, De Jong van Lier Q, Dourado Neto D. Validation of a root water uptake model to estimate transpiration constraints. Agric Water Manage. 2010;97:1382-8.
  • Clark LJ, Whalley WR, Barraclough PB. How do roots penetrate strong soil? Plant Soil. 2003;255:93-104.
  • Couvreur V, Vanderborght J, Javaux M. A simple three-dimensional macroscopic root water uptake model based on the hydraulic architecture approach. Hydrol Earth Syst Sci. 2012;16:2957-71.
  • Cowan IR. Transport of water in the soil-plant-atmosphere system. J Appl Ecol. 1965;2:221-39.
  • Cutford HW, Jefferson PG, Campbell CA. Lower limit of available water for three plant species grown on a medium-textured soil in southwestern Saskatchewan. Can J Soil Sci. 1991;71:247-52.
  • De Jong van Lier Q, Cichota R. Modelagem da aerabilidade do solo. In: Anais da 15ª Reunião Brasileira de Manejo e Conservação do Solo e da Água; 2001; Santa Maria. Santa Maria: Sociedade Brasileira de Ciência do Solo; 2004. v.1.
  • De Jong van Lier Q. Disponibilidade de água às plantas. In: De Jong van Lier Q, editor. Física do solo. Viçosa, MG: Sociedade Brasileira de Ciência do Solo; 2010. p.283-98.
  • De Jong van Lier Q. Oxigenação do sistema radicular: uma abordagem física. R Bras Ci Solo. 2001;25:233-8.
  • De Jong van Lier Q. Revisiting the S-index for soil physical quality and its use in Brazil. R Bras Ci Solo. 2014;38:1-10.
  • De Jong van Lier Q, Dourado Neto D, Metselaar K. Modeling of transpiration reduction in Van Genuchten–Mualem type soils. Water Res Res. 2009;45:W02422. doi:10.1029/2008WR006938
  • De Jong van Lier Q, Metselaar K, Van Dam JC. Root water extraction and limiting soil hydraulic conditions estimated by numerical simulation. Vadose Zone J. 2006;5:1264-77.
  • De Jong van Lier Q, Van Dam JC, Durigon A, Santos MA, Metselaar K. Modeling water potentials and flows in the soil-plant system comparing hydraulic resistances and transpiration reduction functions. Vadose Zone J. 2013;12:1-20.
  • De Jong van Lier Q, Van Dam JC, Metselaar K, De Jong R, Duijnisveld WHM. Macroscopic root water uptake distribution using a matric flux potential approach. Vadose Zone J. 2008;7:1065-78.
  • De Lima CLR, Miola ECC, Timm LC, Pauletto EA, Silva AP. Soil compressibility and least limiting water range of a constructed soil under cover crops after coal mining in Southern Brazil. Soil Till Res. 2012;124:190-5.
  • Dexter AR. Soil physical quality: Part I. Theory, effects of soil texture, density, and organic matter, and effects on root growth. Geoderma. 2004;120:201-14.
  • Doorenbos J, Kassam AH. Yield response to water. Rome: FAO; 1986. (FAO Irrigation and Drainage Paper, 33).
  • Ehlers W, Köpke U, Hesse F, Böhm W. Penetration resistance and growth root of oats in tilled and untilled loess soil. Soil Till Res. 1983;3:261-75.
  • Farias IL, Pacheco EP, Viégas PRA. Characterisation of the optimal hydric interval for a Yellow Argisol cultivated with sugarcane on the coastal plains of Alagoas, Brazil. R Ci Agron. 2013;44:669-75.
  • Feddes RA, Kowalik PL, Zaradny H, editors. Simulation of field water use and crop yield. New York: John Wiley & Sons; 1978.
  • Feddes RA, Kabat P, Van Bakel PJT, Bronswijk JJB, Halbertsma J. Modelling soil water dynamics in the unsaturated zone - state of the art. J Hydrol. 1988;100:69-111.
  • Gardner WR. Dynamic aspects of water availability to plants. Soil Sci. 1960;89:63-7.
  • Gardner WR. Some steady-state solutions of the unsaturated moisture flow equation with application to evaporation from a water table. Soil Sci. 1958;85:228-32.
  • Gilley JE, Doran JW, Eghball B. Tillage and fallow effects on selected soil quality characteristics of former conservation reserve program sites. J Soil Water Conserv. 2001;56:126-32.
  • Gubiani PI, Goulart RZ, Reichert JM, Reinert DJ. Crescimento e produção de milho associados com o intervalo hídrico ótimo. R Bras Ci Solo. 2013b;37:1502-11.
  • Gubiani PI, Reichert JM, Reinert DJ. Indicadores hídrico-mecânicos de compactação do solo e crescimento de plantas. R Bras Ci Solo. 2013a;37:1-10.
  • Guimarães RML, Tormena CA, Blainski E, Fidalski J. Intervalo hídrico ótimo para avaliação da degradação física do solo. R Bras Ci Solo. 2013;37:1512-21.
  • Haise HR, Haas HJ, Jensen LR. Soil mixture studies of some great plain soils. II - Field capacity as related to 1/3-atmosphere percentage and ‘minimum point’ as related to 15- and 26-atmosphere percentage. Soil Sci Soc Am Proc. 1955;34:20-5.
  • Hammond A, Adriaanse A, Rodenburg E, Bryant D, Woodward R. Environmental indicators: A systematic approach to measuring and reporting on environmental policy performance in the context of sustainable development. Washington, D.C.: World Resources Institute; 1995.
  • Hillel D. Introduction to environmental soil physics. London: Academic Press; 2003.
  • Imhoff S, Silva AP, Dias Junior MS, Tormena CA. Quantificação de pressões críticas para o crescimento das plantas. R Bras Ci Solo. 2001;25:11-8.
  • Jarvis NJ. Simple empirical model of root water uptake. J Hydrol. 1989;107:57-72.
  • Jarvis NJ. Simple physics-based models of compensatory plant water uptake: Concepts and eco-hydrological consequences. Hydrol Earth Syst Sci. 2011;15:3431-46.
  • Jastrow JD, Miller RM, Lussenhop J. Contributions of interacting biological mechanisms to soil aggregate stabilization in restored prairie. Soil Biol Biochem. 1998;30:905-16.
  • Javaux M, Couvreur V, Vanderborght J, Vereecken H. Root water uptake: from 3-D biophysical processes to macroscopic modeling approaches, Vadose Zone J. 2013. doi:10.2136/vzj2013.02.0042, 2013
  • Javaux M, Schröder T, Vanderborght J, Vereecken H. Use of a three-dimensional detailed modeling approach for predicting root water uptake. Vadose Zone J. 2008;7:1079-88.
  • Jordan BM, Dumais J. Biomechanics of plant cell growth. In: Encyclopedia of life sciences (ELS). Chichester: John Wiley & Sons; 2010.
  • Kaiser DR, Reinert DJ, Reichert JM, Collares GL, Kunz M. Intervalo hídrico ótimo no perfil explorado pelas raízes de feijoeiro em um Latossolo sob diferentes níveis de compactação. R Bras Ci Solo. 2009;33:845-55.
  • Keller T, Lamandé M, Peth S, Berli M, Delenne J, Baumgarten W, Rabbel W, Radjaï F, Rajchenbach J, Selvadurai A, Or D. An interdisciplinary approach towards improved understanding of soil deformation during compaction. Soil Till Res. 2013;128:61-80.
  • Klein VA, Camara RK. Rendimento da soja e intervalo hídrico ótimo em Latossolo Vermelho sob plantio direto escarificado. R Bras Ci Solo. 2007;31:221-7.
  • Kroes JG, Van Dam JC, Groenendijk P, Hendriks RFA, Jacobs CMJ. SWAP version 3.2. Theory description and user manual. Wageningen: 2008. (Alterra Report, 1649). Disponível em: http://www.swap.alterra.nl.
    » http://www.swap.alterra.nl
  • Kukal SS, Hira GS, Sidhu AS. Soil matric potential-based irrigation scheduling to rice (Oryza sativa). Irrig Sci. 2005;23:153-9.
  • Lascano RJ, Van Bavel CHM. Root water uptake and soil water distribution: test of an availability concept. Soil Sci Soc Am J. 1986;48:233-7.
  • Leão TP, Silva AP, Macedo MCM, Imhoff S, Euclides VPB. Intervalo hídrico ótimo na avaliação de sistemas de pastejo contínuo e rotacionado. R Bras Ci Solo. 2004;28:415-23.
  • Letey J. Relationship between soil physical properties and crop production. Adv Soil Sci. 1985;1:277-94.
  • Li KY, De Jong R, Boisvert JB. An exponential root-water-uptake model with water stress compensation. J Hydrol. 2001;252:189-204.
  • Liu H, Yang H, Zheng J, Jia D, Wang J, Li Y, Huang G. Irrigation scheduling strategies based on soil matric potential on yield and fruit quality of mulched-drip irrigated chili pepper in Northwest China. Agric Water Manage. 2012;115:232-41.
  • Lockhart JA. An analysis of irreversible plant cell elongation. J Theoret Biol. 1965;8:264-75.
  • Metselaar K, De Jong van Lier Q. The shape of the transpiration reduction function under plant water stress. Vadose Zone J. 2007;6:124-39.
  • Meyer WS, Green GC. Water use by wheat and plant indicators of available soil water. Agron J. 1980;72:253-7.
  • Moraes MT, Debiasi H, Carlesso R, Franchini JC, Silva VR. Critical limits of soil penetration resistance in a rhodic Eutrudox. R Bras Ci Solo. 2014;38:288-98.
  • Moreira WH, Tormena CA, Betioli Junior E, Petean LP, Alves SJ. Influência da altura de pastejo de azevém e aveia em atributos físicos de um Latossolo Vermelho distroférrico, após sete anos sob integração lavoura-pecuária. R Bras Ci Solo. 2014b;38:1315-26.
  • Moreira FR, Dechen SCF, Silva AP, Figueiredo GC, Maria IC, Pessoni PT. Intervalo hídrico ótimo em um Latossolo Vermelho cultivado em sistema semeadura direta por 25 anos. R Bras Ci Solo. 2014a;38:118-27.
  • Petean LP, Tormena CA, Alves SJ. Intervalo hídrico ótimo de um Latossolo Vermelho distroférrico sob plantio direto em sistema de integração lavoura-pecuária. R Bras Ci Solo. 2010;34:1515-26.
  • Pullan AJ. The quasilinear approach for unsaturated porous media flow. Water Resour Res. 1990;26:1219-34.
  • Raats PAC. Steady infiltration from line sources and furrows. Soil Sci Soc Am Proc. 1970;34:709-14.
  • Reddell P, Bowen GD, Robson AD. The effects of soil temperature on plant growth, nodulation and nitrogen fixation in Casuarina cunninghamiana Miq. New Phytol. 1985;101:441-50.
  • Reichert JM, Suzuki LEAS, Reinert DJ, Horn R, Hakansson I. Reference bulk density and critical degree-of-compactness for no-till crop production in subtropical highly weathered soils. Soil Till Res. 2009;102:242-54.
  • Reynolds WD, Drury CF, Tan CS, Fox CA, Yang XM. Use of indicators and pore volume-function characteristics to quantify soil physical quality. Geoderma. 2009;152:252-63.
  • Saygin SD, Cornelis WM, Erpul G, Gabriels D. Comparison of different aggregate stability approaches for loamy sand soils. Appl Soil Ecol. 2012;54:1-6.
  • Seben Junior GF, Corá JE, Lal R. Effect of cropping systems in no-till farming on the quality of a Brazilian Oxisol. R Bras Ci Solo. 2014;38:1268-80.
  • Silva AP, Kay BD, Perfect E. Characterization of the least limiting water range. Soil Sci Soc Am J. 1994;58:1775-81.
  • Silva BM, Oliveira GC, Serafim ME, Silva EA, Ferreira MM, Norton LD, Curi N. Critical soil moisture range for a coffee crop in an oxidic Latosol as affected by soil management. Soil Till Res. 2015:154:103-13.
  • Simunek J, van Genuchten MT. Modeling nonequilibrium flow and transport processes using HYDRUS. Vadose Zone J. 2008;7:782-97.
  • Souza LD, Reichardt K. Estimativas da capacidade de campo. R Bras Ci Solo. 1996;20:183-9.
  • Stolf R. Teoria e teste experimental de fórmula de transformação dos dados de penetrômetro de impacto em resistência do solo. R Bras Ci Solo. 1991;5:229-35.
  • Tormena CA, Silva AP, Libardi PL. Caracterização do intervalo hídrico ótimo de um Latossolo Roxo sob plantio direto. R Bras Ci Solo. 1998;22:573-81.
  • Twarakavi NKC, Sakai M, Simunek J. An objective analysis of the dynamic nature of field capacity. Water Resour Res. 2009;45:1-9.
  • van Genuchten MT. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Sci Soc Am J. 1980;44, 892-8.
  • van Genuchten MT. A numerical model for water and solute movement in and below the root zone. Riverside: USDA, ARS; 1987. (U.S. Salinity Laboratory Research Report).
  • Veihmeyer FJ, Hendrickson AH. Does transpiration decrease as the soil moisture decreases? Trans Am Geophys Union. 1955;36:425-48.
  • Veihmeyer FJ, Hendrickson AH. The moisture equivalent as a measure of the field capacity of soils. Soil Sci. 1931:32:181-193.
  • Wright GC, Smith RGC. Differences between two grain sorghum genotypes in adaptation to drought stress. II - Root water uptake and water use. Aust J Agric Res. 1983;34:627-35.

Publication Dates

  • Publication in this collection
    Jul-Aug 2015

History

  • Received
    16 Sept 2014
  • Accepted
    19 Feb 2015
Sociedade Brasileira de Ciência do Solo Sociedade Brasileira de Ciência do Solo, Departamento de Solos - Edifício Silvio Brandão, s/n, Caixa Postal 231 - Campus da UFV, CEP 36570-900 - Viçosa-MG, Tel.: (31) 3612-4542 - Viçosa - MG - Brazil
E-mail: sbcs@sbcs.org.br