Estimation of water retention curve from mercury intrusion porosimetry and van Genuchten model

The water retention curve (WRC) is a hydraulic characteristic of concrete required for advanced modeling of water (and thus solute) transport in variably saturated, heterogeneous concrete. Unfortunately, determination by a direct experimental method (for example, measuring equilibrium moisture levels of large samples stored in constant humidity cells) is a lengthy process, taking over 2 years for large samples. A surrogate approach is presented in which the WRC is conveniently estimated from mercury intrusion porosimetry (MIP) and validated by water sorption isotherms: The well-known Barrett, Joyner and Halenda (BJH) method of estimating the pore size distribution (PSD) from the water sorption isotherm is shown to complement the PSD derived from conventional MIP. This provides a basis for predicting the complete WRC from MIP data alone. The van Genuchten equation is used to model the combined water sorption and MIP results. It is a convenient tool for describing water retention characteristics over the full moisture content range. The van Genuchten parameter estimation based solely on MIP is shown to give a satisfactory approximation to the WRC, with a simple restriction on one. of the parameters.

Alternative Analytical Expressions for the General van Genuchten-Mualem and van Genuchten-Burdine Hydraulic Conductivity Models

The van Genuchten expressions for the unsaturated soil hydraulic properties, first published in 1980, are used frequently in various vadose zone flow and transport applications assuming a specific relationship between the m and n soil hydraulic parameters. By comparison, probably because of the complexity of the hydraulic conductivity equations, the more general solutions with independent m and n values are rarely used. We expressed the general van Genuchten-Mualem and van Genuchten-Burdine hydraulic conductivity equations in terms of hypergeometric functions, which can be approximated by infinite series that converge rapidly for relatively large values of the van Genuchten-Mualem parameter n but only very slowly when n is close to one. Alternative equations were derived that provide very close approximations of the analytical results. The newly proposed equations allow the use of independent values of the parameters m and n in the soil water retention model of van Genuchten for subsequent prediction of the van Genuchten-Mualem and van Genuchten-Burdine hydraulic conductivity models, thus providing more flexibility in fitting experimental pressure-head-dependent water content, theta(h), and hydraulic conductivity, K(h), or K(theta) data.

Modeling of transpiration reduction in van Genuchten-Mualem type soils

We derive an analytic expression for the matric flux potential (M) for van Genuchten-Mualem (VGM) type soils which can also be written in terms of a converging infinite series. Considering the first four terms of this series, the accuracy of the approximation was verified by comparing it to values of M estimated by numerical finite difference integration. Using values of the parameters for three soils from different texture classes, the proposed four-term approximation showed an almost perfect match with the numerical solution, except for effective saturations higher than 0.9. Including more terms reduced the discrepancy but also increased the complexity of the equation. The four-term equation can be used for most applications. Cases with special interest in nearly saturated soils should include more terms from the infinite series. A transpiration reduction function for use with the VGM equations is derived by combining the derived expression for M with a root water extraction model. The shape of the resulting reduction function and its dependency on the derivative of the soil hydraulic diffusivity D with respect to the soil water content theta is discussed. Positive and negative values of dD/d theta yield concave and convex or S-shaped reduction functions, respectively. On the basis of three data sets, the hydraulic properties of virtually all soils yield concave reduction curves. Such curves based solely on soil hydraulic properties do not account for the complex interactions between shoot growth, root growth, and water availability.

Modelagem da curva de retenção de água de Latossolos utilizando a Equação Duplo Van Genuchten

Os Latossolos da região do Cerrado brasileiro, em função da mineralogia da fração argila e estrutura granular, apresentam elevado volume de poros grandes e de poros extremamente pequenos, não mostrando significativo volume de poros intermediários, o que faz com que neles a disponibilidade de água para as plantas seja baixa. Este trabalho objetivou caracterizar e modelar o comportamento da retenção de água em Latossolos oxídicos da região do Cerrado, pertencentes a diferentes classes texturais. Foram coletadas amostras do horizonte Bw de 10 Latossolos oxídicos sob vegetação nativa. A retenção de água nos potenciais matriciais de 1, 2, 4, 6, 8 e 10 kPa foi obtida em unidade de sucção nos potenciais de 33, 60, 100, 500 e 1.500 kPa no extrator de Richards, e a água retida sob potenciais de 1.500 a 300.000 kPa foi quantificada utilizando o psicrômetro de termopar WP4-T. O modelo duplo van Genuchten foi proposto para ajustar os dados experimentais de retenção de água por meio de procedimentos de ajuste de modelos não lineares do software R 2.10.1; também foi avaliada a relação entre as estimativas dos parâmetros do modelo, bem como a inclinação do ponto de inflexão com as propriedades texturais dos solos, aplicando o teste de correlação de Pearson. Os resultados mostraram o bom ajuste do modelo e alto poder de predição, sendo observada correlação do conteúdo de argila do solo com os parâmetros da equação (Usat, Upmp, Ures), assim como, com a inclinação do segundo ponto de inflexão (Itex). A textura dos Latossolos influenciou o comportamento das curvas de retenção de água. As curvas de retenção de água dos Latossolos em estudo apresentaram comportamento bimodal no intervalo de potencial matricial estudado.

Pedotransfer functions to estimate water retention parameters of soils in northeastern Brazil

Pedotransfer functions (PTF) were developed to estimate the parameters (α, n, θr and θs) of the van Genuchten model (1980) to describe soil water retention curves. The data came from various sources, mainly from studies conducted by universities in Northeast Brazil, by the Brazilian Agricultural Research Corporation (Embrapa) and by a corporation for the development of the São Francisco and Parnaíba river basins (Codevasf), totaling 786 retention curves, which were divided into two data sets: 85 % for the development of PTFs, and 15 % for testing and validation, considered independent data. Aside from the development of general PTFs for all soils together, specific PTFs were developed for the soil classes Ultisols, Oxisols, Entisols, and Alfisols by multiple regression techniques, using a stepwise procedure (forward and backward) to select the best predictors. Two types of PTFs were developed: the first included all predictors (soil density, proportions of sand, silt, clay, and organic matter), and the second only the proportions of sand, silt and clay. The evaluation of adequacy of the PTFs was based on the correlation coefficient (R) and Willmott index (d). To evaluate the PTF for the moisture content at specific pressure heads, we used the root mean square error (RMSE). The PTF-predicted retention curve is relatively poor, except for the residual water content. The inclusion of organic matter as a PTF predictor improved the prediction of parameter a of van Genuchten. The performance of soil-class-specific PTFs was not better than of the general PTF. Except for the water content of saturated soil estimated by particle size distribution, the tested models for water content prediction at specific pressure heads proved satisfactory. Predictions of water content at pressure heads more negative than -0.6 m, using a PTF considering particle size distribution, are only slightly lower than those obtained by PTFs including bulk density and organic matter content.

Filter Paper Method for the Determination of the Soil Water Retention Curve

ABSTRACT High cost and long time required to determine a retention curve by the conventional methods of the Richards Chamber and Haines Funnel limit its use; therefore, alternative methods to facilitate this routine are needed. The filter paper method to determine the soil water retention curve was evaluated and compared to the conventional method. Undisturbed samples were collected from five different soils. Using a Haines Funnel and Richards Chamber, moisture content was obtained for tensions of 2; 4; 6; 8; 10; 33; 100; 300; 700; and 1,500 kPa. In the filter paper test, the soil matric potential was obtained from the filter-paper calibration equation, and the moisture subsequently determined based on the gravimetric difference. The van Genuchten model was fitted to the observed data of soil matric potential versus moisture. Moisture values of the conventional and the filter paper methods, estimated by the van Genuchten model, were compared. The filter paper method, with R2 of 0.99, can be used to determine water retention curves of agricultural soils as an alternative to the conventional method.

Bimodal pore distribution on soils under conservation management system for coffee crop

This study aims at detailing bimodal pore distribution by means of water retention curve in an oxidic-gibbsitic Latosol and in a kaolinitic cambisol Latossol under conservation management system of coffee crop. Samples were collected at depths of 20; 40; 80; 120 and 160 cm on coffee trees rows and between rows under oxidic-gibbsitic Latosol (LVd) and kaolinitic cambisol Latossol (LVAd). Water retention curve was determined at matrix potentials (Ψm) -1; -2; -4; -6; -10 kPa obtained from the suction unit; the Ψm of -33; -100; -500; -1,500 kPa were obtained by the Richards extractor, and WP4-T psychrometer was used to determine Ψm -1,500 to -300,000 kPa. The water retention data were adjusted to the double van Genuchten model by nonlinear model procedures of the R 2.12.1 software. Was estimated the model parameter and inflection point slope. The system promoted changes in soil structure and water retention for the conditions evaluated, and both showed bimodal pores distribution, which were stronger in LVd. There was a strong influence of mineralogy gibbsitic in the water retention more negative than Ψm -1500 kPa, reflected in the values of the residual water content.

Número mínimo de tensões para determinação da curva característica de retenção de água de um latossolo vermelho eutrófico sob sistema de semeadura direta

Para o conhecimento do comportamento físico-hídrico do solo, é fundamental a determinação da curva de retenção de água (CRA). A aquisição de dados para obtenção da CRA envolve processos demorados e custos elevados. A hipótese deste estudo foi a de que é possível determinar a CRA com menor número de tensões, reduzindo seu tempo de obtenção, sem, contudo, comprometer a acurácia dos resultados. Assim, o objetivo deste trabalho foi determinar quais tensões podem ser utilizadas para determinação da CRA de Latossolo Vermelho eutrófico de maneira acurada e em menor tempo. Foram determinados os conteúdos de água retidos em amostras de um Latossolo Vermelho eutrófico sob sistema de semeadura direta (SSD). As tensões utilizadas na determinação dos conteúdos de água foram: 0, 10, 20, 40, 60, 80, 100, 200, 300, 500, 700, 1.000, 3.000, 5.000 e 15.000 hPa. Os dados foram combinados e determinaram-se 450 CRAs ajustadas pelo modelo de van Genuchten. Os parâmetros α, m, n e Ug res, gerados pelo modelo, foram submetidos à análise de variância (teste F) e as médias comparadas pelo teste de Scott-Knott (p = 0,05). A hipótese deste estudo foi confirmada, ou seja, é possível determinar a CRA com menor número de tensões, reduzindo seu tempo de obtenção em até cinco vezes, sem, contudo, comprometer a acurácia dos resultados. A tensão de 15.000 hPa deve estar contida na combinação de tensões a serem utilizadas para determinação da curva característica de retenção de água, quando se emprega o modelo matemático de van Genuchten para ajuste dos dados. As curvas características de retenção de água do solo, definidas com as combinações de tensões (0, 60, 700 e 15.000 hPa), (0, 80, 700 e 15.000 hPa), (0, 100, 1.000 e 15.000 hPa), (0, 20, 60, 100, 700 e 15.000 hPa), (0, 80, 300, 700 e 15.000 hPa) e (0, 100, 300, 1.000 e 15.000 hPa), sendo os dados ajustados ao modelo de van Genuchten, proporcionaram estimativas dos conteúdos de água, retidos no Latossolo Vermelho eutrófico sob sistema de semeadura direta, com a mesma acurácia, quando comparada à curva de referência partindo-se de um conjunto de 15 valores de tensões.

Comparison among monitoring strategies to assess water flow dynamic and soil hydraulic properties in agricultural soils

Irrigated agriculture is usually performed in semi-arid regions despite scarcity of water resources. Therefore, optimal irrigation management by monitoring the soil is essential, and assessing soil hydraulic properties and water flow dynamics is presented as a first measure. For this purpose, the control of volumetric water content, θ, and pressure head, h, is required. This study adopted two types of monitoring strategies in the same experimental plot to control θ and h in the vadose zone: i) non-automatic and more time-consuming; ii) automatic connected to a datalogger. Water flux was modelled with Hydrus-1D using the data collected from both acquisition strategies independently (3820 daily values for the automatic; less than 1000 for the non-automatic). Goodness-of-fit results reported a better adjustment in case of automatic sensors. Both model outputs adequately predicted the general trend of θ and h, but with slight differences in computed annual drainage (711 mm and 774 mm). Soil hydraulic properties were inversely estimated from both data acquisition systems. Major differences were obtained in the saturated volumetric water content, θs, and the n and α van Genuchten model shape parameters. Saturated hydraulic conductivity, Ks, shown lower variability with a coefficient of variation range from 0.13 to 0.24 for the soil layers defined. Soil hydraulic properties were better assessed through automatic data acquisition as data variability was lower and accuracy was higher.

Examining the information content of time-lapse crosshole GPR data collected under different infiltration conditions to estimate unsaturated soil hydraulic properties

Time-lapse geophysical data acquired during transient hydrological experiments are being increasingly employed to estimate subsurface hydraulic properties at the field scale. In particular, crosshole ground-penetrating radar (GPR) data, collected while water infiltrates into the subsurface either by natural or artificial means, have been demonstrated in a number of studies to contain valuable information concerning the hydraulic properties of the unsaturated zone. Previous work in this domain has considered a variety of infiltration conditions and different amounts of time-lapse GPR data in the estimation procedure. However, the particular benefits and drawbacks of these different strategies as well as the impact of a variety of key and common assumptions remain unclear. Using a Bayesian Markov-chain-Monte-Carlo stochastic inversion methodology, we examine in this paper the information content of time-lapse zero-offset-profile (ZOP) GPR traveltime data, collected under three different infiltration conditions, for the estimation of van Genuchten-Mualem (VGM) parameters in a layered subsurface medium. Specifically, we systematically analyze synthetic and field GPR data acquired under natural loading and two rates of forced infiltration, and we consider the value of incorporating different amounts of time-lapse measurements into the estimation procedure. Our results confirm that, for all infiltration scenarios considered, the ZOP GPR traveltime data contain important information about subsurface hydraulic properties as a function of depth, with forced infiltration offering the greatest potential for VGM parameter refinement because of the higher stressing of the hydrological system. Considering greater amounts of time-lapse data in the inversion procedure is also found to help refine VGM parameter estimates. Quite importantly, however, inconsistencies observed in the field results point to the strong possibility that posterior uncertainties are being influenced by model structural errors, which in turn underlines the fundamental importance of a systematic analysis of such errors in future related studies.

Bayesian Markov chain Monte Carlo inversion of geophysical data for the purpose of estimating the hydraulic properties of the unsaturated zone

Ground-penetrating radar (GPR) has the potential to provide valuable information on hydrological properties of the vadose zone because of their strong sensitivity to soil water content. In particular, recent evidence has suggested that the stochastic inversion of crosshole GPR data within a coupled geophysical-hydrological framework may allow for effective estimation of subsurface van-Genuchten-Mualem (VGM) parameters and their corresponding uncertainties. An important and still unresolved issue, however, is how to best integrate GPR data into a stochastic inversion in order to estimate the VGM parameters and their uncertainties, thus improving hydrological predictions. Recognizing the importance of this issue, the aim of the research presented in this thesis was to first introduce a fully Bayesian inversion called Markov-chain-Monte-carlo (MCMC) strategy to perform the stochastic inversion of steady-state GPR data to estimate the VGM parameters and their uncertainties. Within this study, the choice of the prior parameter probability distributions from which potential model configurations are drawn and tested against observed data was also investigated. Analysis of both synthetic and field data collected at the Eggborough (UK) site indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when these data are combined with a realistic, informative prior. A subsequent study explore in detail the dynamic infiltration case, specifically to what extent time-lapse ZOP GPR data, collected during a forced infiltration experiment at the Arrenaes field site (Denmark), can help to quantify VGM parameters and their uncertainties using the MCMC inversion strategy. The findings indicate that the stochastic inversion of time-lapse GPR data does indeed allow for a substantial refinement in the inferred posterior VGM parameter distributions. In turn, this significantly improves knowledge of the hydraulic properties, which are required to predict hydraulic behaviour. Finally, another aspect that needed to be addressed involved the comparison of time-lapse GPR data collected under different infiltration conditions (i.e., natural loading and forced infiltration conditions) to estimate the VGM parameters using the MCMC inversion strategy. The results show that for the synthetic example, considering data collected during a forced infiltration test helps to better refine soil hydraulic properties compared to data collected under natural infiltration conditions. When investigating data collected at the Arrenaes field site, further complications arised due to model error and showed the importance of also including a rigorous analysis of the propagation of model error with time and depth when considering time-lapse data. Although the efforts in this thesis were focused on GPR data, the corresponding findings are likely to have general applicability to other types of geophysical data and field environments. Moreover, the obtained results allow to have confidence for future developments in integration of geophysical data with stochastic inversions to improve the characterization of the unsaturated zone but also reveal important issues linked with stochastic inversions, namely model errors, that should definitely be addressed in future research.

Plant-available soil water capacity: estimation methods and implications

The plant-available water capacity of the soil is defined as the water content between field capacity and wilting point, and has wide practical application in planning the land use. In a representative profile of the Cerrado Oxisol, methods for estimating the wilting point were studied and compared, using a WP4-T psychrometer and Richards chamber for undisturbed and disturbed samples. In addition, the field capacity was estimated by the water content at 6, 10, 33 kPa and by the inflection point of the water retention curve, calculated by the van Genuchten and cubic polynomial models. We found that the field capacity moisture determined at the inflection point was higher than by the other methods, and that even at the inflection point the estimates differed, according to the model used. By the WP4-T psychrometer, the water content was significantly lower found the estimate of the permanent wilting point. We concluded that the estimation of the available water holding capacity is markedly influenced by the estimation methods, which has to be taken into consideration because of the practical importance of this parameter.

Assessment of pedotransfer functions for estimating soil water retention curves for the amazon region

Knowledge of the soil water retention curve (SWRC) is essential for understanding and modeling hydraulic processes in the soil. However, direct determination of the SWRC is time consuming and costly. In addition, it requires a large number of samples, due to the high spatial and temporal variability of soil hydraulic properties. An alternative is the use of models, called pedotransfer functions (PTFs), which estimate the SWRC from easy-to-measure properties. The aim of this paper was to test the accuracy of 16 point or parametric PTFs reported in the literature on different soils from the south and southeast of the State of Pará, Brazil. The PTFs tested were proposed by Pidgeon (1972), Lal (1979), Aina & Periaswamy (1985), Arruda et al. (1987), Dijkerman (1988), Vereecken et al. (1989), Batjes (1996), van den Berg et al. (1997), Tomasella et al. (2000), Hodnett & Tomasella (2002), Oliveira et al. (2002), and Barros (2010). We used a database that includes soil texture (sand, silt, and clay), bulk density, soil organic carbon, soil pH, cation exchange capacity, and the SWRC. Most of the PTFs tested did not show good performance in estimating the SWRC. The parametric PTFs, however, performed better than the point PTFs in assessing the SWRC in the tested region. Among the parametric PTFs, those proposed by Tomasella et al. (2000) achieved the best accuracy in estimating the empirical parameters of the van Genuchten (1980) model, especially when tested in the top soil layer.

Stabilization of soil hydraulic properties under a long term no-till system

The area under the no-tillage system (NT) has been increasing over the last few years. Some authors indicate that stabilization of soil physical properties is reached after some years under NT while other authors debate this. The objective of this study was to determine the effect of the last crop in the rotation sequence (1st year: maize, 2nd year: soybean, 3rd year: wheat/soybean) on soil pore configuration and hydraulic properties in two different soils (site 1: loam, site 2: sandy loam) from the Argentinean Pampas region under long-term NT treatments in order to determine if stabilization of soil physical properties is reached apart from a specific time in the crop sequence. In addition, we compared two procedures for evaluating water-conducting macroporosities, and evaluated the efficiency of the pedotransfer function ROSETTA in estimating the parameters of the van Genuchten-Mualem (VGM) model in these soils. Soil pore configuration and hydraulic properties were not stable and changed according to the crop sequence and the last crop grown in both sites. For both sites, saturated hydraulic conductivity, K0, water-conducting macroporosity, εma, and flow-weighted mean pore radius, R0ma, increased from the 1st to the 2nd year of the crop sequence, and this was attributed to the creation of water-conducting macropores by the maize roots. The VGM model adequately described the water retention curve (WRC) for these soils, but not the hydraulic conductivity (K) vs tension (h) curve. The ROSETTA function failed in the estimation of these parameters. In summary, mean values of K0 ranged from 0.74 to 3.88 cm h-1. In studies on NT effects on soil physical properties, the crop effect must be considered.

Estimating vadose zone hydraulic properties using ground penetrating radar: The impact of prior information

A number of geophysical methods, such as ground-penetrating radar (GPR), have the potential to provide valuable information on hydrological properties in the unsaturated zone. In particular, the stochastic inversion of such data within a coupled geophysical-hydrological framework may allow for the effective estimation of vadose zone hydraulic parameters and their corresponding uncertainties. A critical issue in stochastic inversion is choosing prior parameter probability distributions from which potential model configurations are drawn and tested against observed data. A well chosen prior should reflect as honestly as possible the initial state of knowledge regarding the parameters and be neither overly specific nor too conservative. In a Bayesian context, combining the prior with available data yields a posterior state of knowledge about the parameters, which can then be used statistically for predictions and risk assessment. Here we investigate the influence of prior information regarding the van Genuchten-Mualem (VGM) parameters, which describe vadose zone hydraulic properties, on the stochastic inversion of crosshole GPR data collected under steady state, natural-loading conditions. We do this using a Bayesian Markov chain Monte Carlo (MCMC) inversion approach, considering first noninformative uniform prior distributions and then more informative priors derived from soil property databases. For the informative priors, we further explore the effect of including information regarding parameter correlation. Analysis of both synthetic and field data indicates that the geophysical data alone contain valuable information regarding the VGM parameters. However, significantly better results are obtained when we combine these data with a realistic, informative prior.