BAYESIAN ESTIMATION OF HYDRAULIC AND SOLUTE TRANSPORT PARAMETERS FROM LABORATORY SOIL COLUMN EXPERIMENTS

Author(s):  
Helcio R. B. Orlande ◽  
Martinus Th. van Genuchten ◽  
Renato M. Cotta ◽  
Paulo H. Moreira
2015 ◽  
Vol 10 (1) ◽  
pp. 139
Author(s):  
Panom Chaiyasit ◽  
Piya Duangpatra ◽  
Visoot Verasan ◽  
Varawoot Vudhivanich

<p class="zhengwen"><span lang="EN-GB">An experiment was conducted on the purpose to study movement of water and salt through soil column. Salt-affected paddy soil was assessed for its relevant transport parameters consisting of the hydraulic and the solute transport parameters. The hydraulic parameters included soil hydraulic conductivity (K<sub>s</sub>) and the van Genuchten’s parameters (θ<sub>s</sub>, θ<sub>r</sub>, α, and n). In this experiment the solute transport parameters was referred to the coefficient of Langmuir’s isotherm which consisted of k<sub>d</sub> and η. Experience showed that hydraulic parameters were sensitive to changes of soil bulk density (ρ<sub>b</sub>). Therefore pedotransfer functions describing the relations between these parameters with ρ<sub>b</sub> were established. Straight line functions were found for θ<sub>s</sub> and n, exponential functions were found for α and K<sub>s</sub>, and logarithmic function was found for θ<sub>r</sub>. Packing the soil in the physical column inevitably caused horizontal differentiation of different ρ<sub>b</sub>. Bulk density of each layer was estimated by analysis of water flow through soil column at steady-state. Then ρ<sub>b</sub> of each layer was calculated from the relation K<sub>s</sub> (ρ<sub>b</sub>). After the ρ<sub>b</sub> was known the van Genuchten’s parameters were calculated from the pedotransfer functions. A physical column of 4 inches diameter and 50 cm length was constructed. Sodium chloride solution EC 6 dS/m was fed on soil surface during the process of salinization and the feeding solution was changed to fresh water during the process of desalinization. Breakthrough solution was analyzed for Na concentration and the breakthrough curves were constructed. The relevant parameters as well as initial and boundary conditions were fed into Hydrus-1D on the purpose to simulate the breakthrough curves. Statistical comparison results using t-test and RMSE suggested that Hydrus-1D could be used successfully to monitor transport of water and salt through soil column.</span></p><p class="zhengwen"><span lang="EN-GB">Five scenarios concerning water and solute transport through soil profile under rice and mung bean cropping were simulated using Hydrus-1D. Simulation results suggested that continuous flooding was the most efficient way to leach soluble salts down to ground water. Wet/dry irrigation scheme for rice production could drain salts only when rice crop was in the first period of growth where crop water uptake was small. During later stages of growth concentration profile of Na remained almost unchange indicating negligible downward movement of salts. Leaving the soil to dry naturally during the dry season caused upward accumulation of salt to the degree smaller than cultivating mung bean since water content and hence the diffusion coefficient of the soil were higher and more favorable for upward salt diffusion than in the former case. Inserting the capillary rise hindering soil layer underneath mung bean root zone was found to retard upward diffusion of salt to the degree comparable to leaving the soil to dry naturally.</span></p>


2012 ◽  
Vol 45 (7) ◽  
pp. 675-683 ◽  
Author(s):  
E. Raymundo-Raymundo ◽  
Yu. N. Nikolskii ◽  
A. K. Guber ◽  
C. Landeros-Sanchez

2016 ◽  
Vol 64 (1) ◽  
pp. 30-44 ◽  
Author(s):  
Paulo H. S. Moreira ◽  
Martinus Th. van Genuchten ◽  
Helcio R. B. Orlande ◽  
Renato M. Cotta

Abstract In this study the hydraulic and solute transport properties of an unsaturated soil were estimated simultaneously from a relatively simple small-scale laboratory column infiltration/outflow experiment. As governing equations we used the Richards equation for variably saturated flow and a physical non-equilibrium dual-porosity type formulation for solute transport. A Bayesian parameter estimation approach was used in which the unknown parameters were estimated with the Markov Chain Monte Carlo (MCMC) method through implementation of the Metropolis-Hastings algorithm. Sensitivity coefficients were examined in order to determine the most meaningful measurements for identifying the unknown hydraulic and transport parameters. Results obtained using the measured pressure head and solute concentration data collected during the unsaturated soil column experiment revealed the robustness of the proposed approach.


Water ◽  
2019 ◽  
Vol 11 (8) ◽  
pp. 1660 ◽  
Author(s):  
Anis Younes ◽  
Jabran Zaouali ◽  
Sabri Kanzari ◽  
Francois Lehmann ◽  
Marwan Fahs

Numerical modeling has become an irreplaceable tool for the investigation of water flow and solute transport in the unsaturated zone. The use of this tool for real situations is often faced with lack of knowledge of hydraulic and soil transport parameters. In this study, advanced experimental and numerical techniques are developed for an accurate estimation of the soil parameters. A laboratory unsaturated flow and solute transport experiment is conducted on a large undisturbed soil column of around 40 cm length. Bromide, used as a nonreactive contaminant, is injected at the surface of the undisturbed soil, followed by a leaching phase. The pressure measurements at different locations along the soil column as well as the outflow bromide concentration are collected during the experiment and used for the statistical calibration of flow and solute transport. The Richards equation, combined with constitutive relations for water content and permeability, is used to describe unsaturated flow. Both linear and non-equilibrium mobile–immobile transport models are investigated for the solute transport. All hydraulic and mass transport parameters are inferred using a one-step Bayesian estimation with the Markov chain Monte Carlo sampler. The results prove that the pressure and concentration measurements are able to identify almost all hydraulic and mass transport parameters. The mobile–immobile transport model better reproduces the infiltration experiment. It produces narrower uncertainty intervals for soil parameters and predictive output concentrations.


2008 ◽  
Vol 43 (8) ◽  
pp. 732-741 ◽  
Author(s):  
Zoi Magga ◽  
Dimitra N. Tzovolou ◽  
Maria A. Theodoropoulou ◽  
Theodora Dalkarani ◽  
Konstantinos Pikios ◽  
...  

2015 ◽  
Vol 19 (6) ◽  
pp. 2617-2635 ◽  
Author(s):  
M. Sprenger ◽  
T. H. M. Volkmann ◽  
T. Blume ◽  
M. Weiler

Abstract. Determining the soil hydraulic properties is a prerequisite to physically model transient water flow and solute transport in the vadose zone. Estimating these properties by inverse modelling techniques has become more common within the last 2 decades. While these inverse approaches usually fit simulations to hydrometric data, we expanded the methodology by using independent information about the stable isotope composition of the soil pore water depth profile as a single or additional optimization target. To demonstrate the potential and limits of this approach, we compared the results of three inverse modelling strategies where the fitting targets were (a) pore water isotope concentrations, (b) a combination of pore water isotope concentrations and soil moisture time series, and (c) a two-step approach using first soil moisture data to determine water flow parameters and then the pore water stable isotope concentrations to estimate the solute transport parameters. The analyses were conducted at three study sites with different soil properties and vegetation. The transient unsaturated water flow was simulated by solving the Richards equation numerically with the finite-element code of HYDRUS-1D. The transport of deuterium was simulated with the advection-dispersion equation, and a modified version of HYDRUS was used, allowing deuterium loss during evaporation. The Mualem–van Genuchten and the longitudinal dispersivity parameters were determined for two major soil horizons at each site. The results show that approach (a), using only the pore water isotope content, cannot substitute hydrometric information to derive parameter sets that reflect the observed soil moisture dynamics but gives comparable results when the parameter space is constrained by pedotransfer functions. Approaches (b) and (c), using both the isotope profiles and the soil moisture time series, resulted in good simulation results with regard to the Kling–Gupta efficiency and good parameter identifiability. However, approach (b) has the advantage that it considers the isotope data not only for the solute transport parameters but also for water flow and root water uptake, and thus increases parameter realism. Approaches (b) and (c) both outcompeted simulations run with parameters derived from pedotransfer functions, which did not result in an acceptable representation of the soil moisture dynamics and pore water stable isotope composition. Overall, parameters based on this new approach that includes isotope data lead to similar model performances regarding the water balance and soil moisture dynamics and better parameter identifiability than the conventional inverse model approaches limited to hydrometric fitting targets. If only data from isotope profiles in combination with textural information is available, the results are still satisfactory. This method has the additional advantage that it will not only allow us to estimate water balance and response times but also site-specific time variant transit times or solute breakthrough within the soil profile.


Sign in / Sign up

Export Citation Format

Share Document