scholarly journals HYDRUS-2D Simulations of nitrate nitrogen and potassium transport characteristics under fertilizer solution infiltration of furrow irrigation

Author(s):  
Wei-Bo Nie ◽  
Kun-Kun Nie ◽  
Yi-Bo Li ◽  
Xiao-Yi Ma

Abstract Understanding the characteristics of soil solute transport is fundamental to the design and management of furrow irrigation systems. This study determined the soil hydraulic and solute transport parameters by inverse solution with HYDRUS-2D and then verified them. The experimental data were obtained from the infiltration of clay loam and sandy loam of different potassium nitrate (KNO3) concentrations under furrow irrigation. Then, the initial soil water content (θ0), KNO3 concentration, and water depth (h0) affecting the transport characteristics of nitrate nitrogen (NO3−-N) and potassium (K+) were analyzed. The results indicated that the soil hydraulic and solute transport parameters determined from the inversion solution with HYDRUS-2D were reliable. The soil saturated water content, saturated hydraulic conductivity, and empirical parameter n in the van Genuchten–Mualem model increase with the increase of KNO3 concentrations, whereas the empirical parameter a shows a decreasing tendency. The distribution range of NO3−-N increased with the increases of θ0 and the KNO3 concentration, which had barely any effect on the range of K+ distribution. The horizontal distribution range of NO3−-N and K+ increased with the increase of h0, but it had no obvious influence on the vertical range.

1997 ◽  
Vol 1 (4) ◽  
pp. 801-811 ◽  
Author(s):  
B. Lennartz ◽  
S. K. Kamra ◽  
S. Meyer-Windel

Abstract. The spatial variability of transport parameters has to be taken into account for a reliable assessment of solute behaviour in natural field soils. Two field sites were studied by collecting 24 and 36 small undisturbed soil columns at an uniform grid of 15 m spacing. Displacement experiments were conducted in these columns with bromide traced water under unsaturated steady state transport conditions. Measured breakthrough curves (BTCs) were evaluated with the simple convective-dispersive equation (CDE). The solute mobility index (MI) calculated as the ratio of measured to fitted pore water velocity and the dispersion coefficient (D) were used to classify bromide breakthrough behaviour. Experimental BTCs were classified into two groups: type I curves expressed classical solute behaviour while type II curves were characterised by the occurrence of a bromide concentration maximum before 0.35 pore volumes of effluent (MI<0.35) resulting from preferential flow conditions. Six columns from site A and 8 from site B were identified as preferential. Frequency distributions of the transport parameters (MI and D) of both sites were either extremely skewed or bimodal. Log-transformation did not lead to a normal distribution in any case. Contour maps of bromide mass flux at certain time steps indicated the clustering of preferential flow regions at both sites. Differences in the extent of preferential flow between sites seemed to be governed by soil structure. Linear cross correlations among transport parameters and independently measured soil properties revealed relations between solute mobility and volumetric soil water content at time of sampling, texture and organic carbon content. The volumetric field soil water content, a simple measure characterising the soil hydraulic behaviour at the sampling location, was found to be a highly sensitive parameter with respect to solute mobility and preferential flow situations. Almost no relation was found between solute transport parameters and independently determined soil properties when non-preferential and preferential samples were considered separately in regression analyses. Future work should concentrate to relate integrated parameters such as the infiltration rate or the soil hydraulic functions to solute mobility under different flow situations.


1998 ◽  
Vol 62 (5) ◽  
pp. 1172-1178 ◽  
Author(s):  
Francis X. M. Casey ◽  
Robert Horton ◽  
Sally D. Logsdon ◽  
Dan B. Jaynes

2020 ◽  
Vol 34 (3) ◽  
pp. 310-324
Author(s):  
Leonardo Ezequiel Scherger ◽  
Victoria Zanello ◽  
Claudio Lexow

The aim of this work is to compare the use of the inverse solution approach in the estimation of soil hydraulic properties with traditional tension disk infiltrometer (TDI) data analysis, field retention data and commonly used pedotransfer functions (PTFs). Field data were collected in an experimental plot located at Bahía Blanca, Argentina. Field infiltration under saturated conditions was measured by the inverse auger hole method and infiltration under unsaturated conditions were carried out with TDI. Field retention data (θ(h)) were also collected periodically. The HYDRUS 2D/3D software was used to optimize soil hydraulic parameters by inverse solution according to TDI data. The saturated hydraulic conductivity measured by inverse auger hole method (5.53 cm.h-1) and calculated by Wooding analytical approach (5.35 cm.h-1) and inverse numerical simulations (5.36 cm.h-1) showed very close values. According to macroporosity estimates infiltrated water is mainly conducted through soils micro and mesopores.  Macropores only channeled 15.9% of total infiltrated flow.  Soil water retention curves (SWRC) predicted by PTFs did not represented correctly field retention data. The best adjustment between water content at specific pressure heads predicted by SWRCs and field measured water content was reached by the TDI inverse solution approach (RMSE: 0.050 cm3.cm-3). The inverse solution approach probed to be a simple and practical method to obtain an accurate estimate of both, SWRC and hydraulic conductivity curve.


2007 ◽  
Vol 46 (8) ◽  
pp. 1275-1289 ◽  
Author(s):  
Gerd Schädler

Abstract Continuous time series of soil water content over a period of more than 9 months for a midlatitude sandy loam soil covered by grass are calculated with the Campbell and the van Genuchten soil hydraulic functions and the Clapp–Hornberger, Cosby et al., and Rawls–Brakensiek parameter sets. The results are compared with soil water content observed at several soil depths, and the water balance components are evaluated. The Campbell soil hydraulic functions are often used by meteorologists, whereas the van Genuchten functions are widespread among hydrologists. The simulations are performed with the “VEG3D” soil–vegetation model in stand-alone mode forced by on-site meteorological observations. The soil water content and meteorological observations were obtained within the Regional Climate Project (REKLIP) at a site in the Rhine valley in southern Germany with 10-min temporal resolution. Apart from the different soil hydraulic functions and parameter sets, the effects of different lower boundary conditions and initializations on the simulations are compared in terms of statistical quantities like mean error, bias, correlation coefficient, and least squares fit. Large differences between the various combinations are found. For the situation considered in this paper, the van Genuchten–Clapp–Hornberger, the Campbell–Cosby et al., and the van Genuchten–Rawls–Brakensiek combinations give the best overall agreement with the observations.


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.


1991 ◽  
Vol 22 (5) ◽  
pp. 275-302 ◽  
Author(s):  
K. Høgh Jensen ◽  
J. C. Refsgaard

Natural field systems exhibit a large degree of soil heterogeneity which affects the movement of water and solutes and thus leads to highly varying observations of water content and solute concentration. To investigate this problem comprehensive field investigation programs were carried out at two field sites in Denmark representing two different soil types, a coarse sand and a sandy loam, respectively. The field investigations included collection of soil samples for analysis of textural composition, retention, and hydraulic conductivity, measurements of water content and suction, and measurements of radioactive tracer concentration, all carried out at a number of positions within the two field sites. Models for one-dimensional vertical unsaturated flow and solute transport were applied to the two field sites, and the simulation results were compared to field measurements of water content, suction and solute concentration. This paper describes results from model simulations in individual soil profiles, while the variability issues at field scale are described in the two accompanying papers. The modelling approach was based on numerical solutions to Richards' equation for water flow and the convection-dispersion equation (CDE) for solute transport. The model results from the coarse sand field site compared relatively well to measurements of water content, suction, and concentration except for the upper soil layer (∼ 10 cm depth) where the measured water contents appeared to be somewhat uncertain. Due to the neglecting of hysteresis and macropore flow (by-pass) in the model the measured retention curves (drainage based) and the hydraulic conductivity functions at the sandy loam field site had to he modified empirically through the calibration procedure in order to match the measurements.


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.


Sign in / Sign up

Export Citation Format

Share Document