scholarly journals Study on Movement of Water and Salt through Soil Column and Utilization of Hydrus-1D Program to Simulate Five Scenarios of Crop Production in Salt Affected Paddy Soil

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>

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.


2020 ◽  
Author(s):  
Chaoqi Wang ◽  
Xiaoguang Wang ◽  
Samer Majdalani ◽  
Vincent Guinot ◽  
Hervé Jourde

&lt;p&gt;An important phenomenon often encountered when interpreting tracer tests in karst aquifers is the occurrence of dual-peaked breakthrough curves (BTCs). The dual-peaked BTCs are usually attributed to tracer transport through a conduit system consisting of a dual-conduit structure: an auxiliary conduit that deviates from the main conduit at the upstream and converges back at the downstream. In order to understand how the geometric configuration of the dual-conduit structure influences the BTCs, laboratory experiments utilizing plastic tubes were conducted. The physical models were constructed by varying: 1) the total length of the conduits, while fixing the length ratio; 2) length ratio between the two conduits, while fixing the length of the main conduit; and 3) conduits connection angle. The tracer experiments are then fitted by a Multi-Region Advection Dispersion model and a Transfer Function model to derive effective transport parameters. This allows us to quantitatively compare the experimental results, and thus to analyse the conduit geometry effects on solute transport and to compare the performance of the two models.&lt;/p&gt;&lt;p&gt;Results show that the dual-conduit structure causes the double peaks of BTCs. Keeping the length ratio of the two conduits and increasing their total length leads to a larger separation of the two peaks of the BTCs. Keeping the length of main conduit while increasing the length of the secondary conduit causes similar effects. As (&amp;#952;&lt;sub&gt;1&lt;/sub&gt;-&amp;#952;&lt;sub&gt;2&lt;/sub&gt;) increases, the first peak concentration value decreases, the second peak concentration value increases.&lt;/p&gt;&lt;p&gt;&lt;strong&gt;Keywords&lt;/strong&gt;: karst, lab experiment, dual-peaked BTCs, modelling&lt;/p&gt;


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.


Soil Research ◽  
2006 ◽  
Vol 44 (2) ◽  
pp. 175 ◽  
Author(s):  
M. A. Mojid ◽  
D. A. Rose ◽  
G. C. L. Wyseure

A transfer-function method has been applied to determine solute-transport parameters from earlier sections of complete breakthrough data. Time-domain reflectometry allows the measurement of breakthrough data in unsaturated soil. In fine-textured soils, the flow of water must be kept low to maintain unsaturated conditions, and so experiments for a complete breakthrough of solute may last a very long time. Substantial savings of time and computer memory might be achieved if data could be analysed from an earlier section of breakthrough data. Data at 2 vertical positions (input at upper and response at lower position) from a complete breakthrough of calcium chloride applied as a pulse input to 4 unsaturated soils (coarse sand, sandy loam, clay loam, clay) were divided into 4 sets of increasing duration. Transport parameters of calcium chloride were determined by a transfer function, which results in similar values of the parameters from the last 3 datasets in all 4 soils. In the clay soil, however, because of erroneous breakthrough data the fit between the measured and estimated breakthrough curves (BTCs) was poor, but the transport parameters were consistent among different segments of data. We show that it is possible to determine successfully solute-transport parameters from partial breakthrough data, which include the peak of the response BTC. This transfer-function method is thus a powerful tool to shorten breakthrough experiments.


2002 ◽  
Vol 33 (2-3) ◽  
pp. 145-164 ◽  
Author(s):  
Magnus Persson ◽  
Ronny Berndtsson

Two quasi steady-state solute transport experiments were carried out in a loamy sand under field conditions. The flux was 40 mm/d in experiment 1 and 18.7 mm/d in experiment 2. Both water content (θ) and resident concentration (Cr) measurements were taken using 64 time domain reflectometry probes at depths ranging from 0.05 to 0.90 m. The Cr measurement was calibrated in situ for each probe location in the field. The convective dispersive equation (CDE) and convective lognormal transfer function (CLT) models were fitted to the breakthrough curves (BTCs). The results indicated fingered flow, which has been shown to exist in previous studies of this soil. The finger width was larger in experiment 1 leading to smaller horizontal heterogeneity and a relatively smaller solute transport velocity. The location of the fingers was consistent between the two experiments resulting in a high correlation between the velocity and mass balance fields. Mass balance calculations showed that the solute mass integrated over depth one day after the solute application was larger than the mass balance for the entire experiment (integrated over time). The probable reason being that solutes were transported out of the measurement volume by horizontal flow across the Ap/E horizon boundary. The investigation of the transport parameters revealed that both the CDE and CLT models could be successfully used to predict most individual BTCs. Horizontally averaged global CDE and CLT models were also fitted to the data. Global solute transport was better modeled with the CDE model in experiment 1, while in experiment 2, the CLT model was better. This study clearly shows the applicability of using TDR with the in situ calibration technique in field experiments with varying water content.


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.


2014 ◽  
Vol 11 (10) ◽  
pp. 11203-11245 ◽  
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 two 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 numerically solving the Richards equation with the finite-element code of Hydrus-1D. The transport of deuterium was simulated with the advection-dispersion equation, and the Hydrus code was modified to allow for 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 satisfying model performances 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 thus increases parameter realism. Approaches (b) and (c) both outcompeted simulations run with parameters derived from pedotransfer functions, which did not result in satisfying model efficiencies. 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.


Author(s):  
Muskan Mayank ◽  
Pramod Kumar Sharma

Abstract Environmental concerns have drawn much research interest in solute transport through porous media. Thus, contaminants of groundwater permeate through pores in the ground, and adsorption attenuates the pollution concentration as the pollutants adhere to the solid surface. Mathematical models based on certain simplifying assumptions have been used to predict solute transport. The transport of solutes in porous media is governed by a partial differential equation known as the advection-dispersion equation. In this study, a two-dimensional numerical model has been developed for solute transport through porous media. Results of spatial moments have been predicted and analysed in the presence of both constant and time-dependent dispersion coefficients. Afterward, a numerical model is used to simulate experimentally observed breakthrough curves for both conservative and non-conservative solutes. Thus, transport parameters are estimated through numerical simulation of observed breakthrough curves. Finally, this model gives the best simulation of observed breakthrough curves, and it can also be used in the field scale.


Sign in / Sign up

Export Citation Format

Share Document