scholarly journals Determination of Soil Hydraulic Parameters and Evaluation of Water Dynamics and Nitrate Leaching in the Unsaturated Layered Zone: A Modeling Case Study in Central Croatia

2021 ◽  
Vol 13 (12) ◽  
pp. 6688
Author(s):  
Jasmina Defterdarović ◽  
Lana Filipović ◽  
Filip Kranjčec ◽  
Gabrijel Ondrašek ◽  
Diana Kikić ◽  
...  

Nitrate leaching through soil layers to groundwater may cause significant degradation of natural resources. The aims of this study were: (i) to estimate soil hydraulic properties (SHPs) of the similar soil type with same management on various locations; (ii) to determine annual water dynamics; and (iii) to estimate the impact of subsoil horizon properties on nitrate leaching. The final goal was to compare the influence of different SHPs and layering on water dynamics and nitrate leaching. The study was conducted in central Croatia (Zagreb), at four locations on Calcaric Phaeozem, Calcaric Regosol, and Calcaric Fluvic Phaeozem soil types. Soil hydraulic parameters were estimated using the HYPROP system and HYPROP-FIT software. Water dynamics and nitrate leaching were evaluated using HYDRUS 2D/3D during a period of 365 days. The amount of water in the soil under saturated conditions varied from 0.422 to 0.535 cm3 cm–3 while the hydraulic conductivity varied from 3 cm day−1 to 990.9 cm day−1. Even though all locations have the same land use and climatic conditions with similar physical properties, hydraulic parameters varied substantially. The amount and velocity of transported nitrate (HYDRUS 2D/3D) were affected by reduced hydraulic conductivity of the subsoil as nitrates are primarily transported via advective flux. Despite the large differences in SHPs of the topsoil layers, the deeper soil layers, having similar SHPs, imposed a buffering effect preventing faster nitrate downward transport. This contributed to a very similar distribution of nitrates through the soil profile at the end of simulation period. This case study indicated the importance of carefully selecting relevant parameters in multilayered soil systems when evaluating groundwater pollution risk.

Water ◽  
2020 ◽  
Vol 12 (7) ◽  
pp. 1858 ◽  
Author(s):  
Jesús María Domínguez-Niño ◽  
Gerard Arbat ◽  
Iael Raij-Hoffman ◽  
Isaya Kisekka ◽  
Joan Girona ◽  
...  

Although surface drip irrigation allows an efficient use of water in agriculture, the heterogeneous distribution of soil water complicates its optimal usage. Mathematical models can be used to simulate the dynamics of water in the soil below a dripper and promote: a better understanding, and optimization, of the design of drip irrigation systems, their improved management and their monitoring with soil moisture sensors. The aim of this paper was to find the most appropriate configuration of HYDRUS-3D for simulating the soil water dynamics in a drip-irrigated orchard. Special emphasis was placed on the source of the soil hydraulic parameters. Simulations parameterized using the Rosetta approach were therefore compared with others parameterized using that of HYPROP + WP4C. The simulations were validated on a seasonal scale, against measurements made using a neutron probe, and on the time course of several days, against tensiometers. The results showed that the best agreement with soil moisture measurements was achieved with simulations parameterized from HYPROP + WP4C. It further improved when the shape parameter n was empirically calibrated from a subset of neutron probe measurements. The fit of the simulations with measurements was best at positions near the dripper and worsened at positions outside its wetting pattern and at depths of 80 cm or more.


2020 ◽  
Author(s):  
Jesús Fernández-Gálvez ◽  
Joseph Pollacco ◽  
Laurent Lassabatere ◽  
Rafael Angulo-Jaramillo ◽  
Sam Carrick

<p>Soil hydraulic characterization is crucial to describe the retention and transport of water in soil, but current methodologies limit its spatial applicability. This work presents a cost-effective general Beerkan Estimation of Soil Transfer parameters (BEST) methodology using single ring infiltration experiments to derive soil hydraulic parameters for any type of unimodal water retention and hydraulic conductivity functions. The proposed method relies on the BEST approach. The novelty lies in the use of Kosugi hydraulic parameters without need for textural information. Kosugi functions were chosen because they are based on physical principles (log-normal distribution for pore size distributions). A link between the Kosugi parameters (i.e., relationship between <em>σ</em> and <em>h</em><sub>kg</sub>) was introduced to reduce the number of parameters estimated and to avoid the need for information on the soil texture. This simplifies the procedures and avoids sources of errors related to the use of pedotransfer functions as for the previous BEST methods. Lastly, the method uses a quasi-exact formulation that is valid for all times, instead of the approximate expansions previously used, avoiding related inaccuracy and allowing the use of any infiltration data encompassing or not both transient and steady states. The new BEST methods were tested against numerically generated data for several contrasting synthetic soils, and the results show that these methods provide consistent hydraulic functions close to the target functions. The new BEST method is accurate and can use any type of water retention and hydraulic conductivity functions (Fernández-Gálvez et al., 2019).</p><p> </p><p> </p><p><strong>Reference</strong></p><p>Fernández-Gálvez, J., Pollacco, J.A.P., Lassabatere, L., Angulo-Jaramillo, R., Carrick, S., 2019. A general Beerkan Estimation of Soil Transfer parameters method predicting hydraulic parameters of any unimodal water retention and hydraulic conductivity curves: Application to the Kosugi soil hydraulic model without using particle size distribution data. Adv. Water Resour. 129, 118–130. https://doi.org/10.1016/j.advwatres.2019.05.005</p>


2011 ◽  
Vol 15 (10) ◽  
pp. 3043-3059 ◽  
Author(s):  
B. Scharnagl ◽  
J. A. Vrugt ◽  
H. Vereecken ◽  
M. Herbst

Abstract. In situ observations of soil water state variables under natural boundary conditions are often used to estimate the soil hydraulic properties. However, many contributions to the soil hydrological literature have demonstrated that the information content of such data is insufficient to accurately and precisely estimate all the soil hydraulic parameters. In this case study, we explored to which degree prior information about the soil hydraulic parameters can help improve parameter identifiability in inverse modelling of in situ soil water dynamics under natural boundary conditions. We used percentages of sand, silt, and clay as input variables to the ROSETTA pedotransfer function that predicts the parameters in the van Genuchten-Mualem (VGM) model of the soil hydraulic functions. To derive additional information about the correlation structure of the predicted parameters, which is not readily provided by ROSETTA, we employed a Monte Carlo approach. We formulated three prior distributions that incorporate to different extents the prior information about the VGM parameters derived with ROSETTA. The inverse problem was posed in a formal Bayesian framework and solved using Markov chain Monte Carlo (MCMC) simulation with the DiffeRential Evolution Adaptive Metropolis (DREAM) algorithm. Synthetic and real-world soil water content data were used to illustrate the approach. The results of this study demonstrated that prior information about the soil hydraulic parameters significantly improved parameter identifiability and that this approach was effective and robust, even in case of biased prior information. To be effective and robust, however, it was essential to use a prior distribution that incorporates information about parameter correlation.


2010 ◽  
Vol 14 (2) ◽  
pp. 251-270 ◽  
Author(s):  
G. Baroni ◽  
A. Facchi ◽  
C. Gandolfi ◽  
B. Ortuani ◽  
D. Horeschi ◽  
...  

Abstract. Data of soil hydraulic properties forms often a limiting factor in unsaturated zone modelling, especially at the larger scales. Investigations for the hydraulic characterization of soils are time-consuming and costly, and the accuracy of the results obtained by the different methodologies is still debated. However, we may wonder how the uncertainty in soil hydraulic parameters relates to the uncertainty of the selected modelling approach. We performed an intensive monitoring study during the cropping season of a 10 ha maize field in Northern Italy. The data were used to: i) compare different methods for determining soil hydraulic parameters and ii) evaluate the effect of the uncertainty in these parameters on different variables (i.e. evapotranspiration, average water content in the root zone, flux at the bottom boundary of the root zone) simulated by two hydrological models of different complexity: SWAP, a widely used model of soil moisture dynamics in unsaturated soils based on Richards equation, and ALHyMUS, a conceptual model of the same dynamics based on a reservoir cascade scheme. We employed five direct and indirect methods to determine soil hydraulic parameters for each horizon of the experimental profile. Two methods were based on a parameter optimization of: a) laboratory measured retention and hydraulic conductivity data and b) field measured retention and hydraulic conductivity data. The remaining three methods were based on the application of widely used Pedo-Transfer Functions: c) Rawls and Brakensiek, d) HYPRES, and e) ROSETTA. Simulations were performed using meteorological, irrigation and crop data measured at the experimental site during the period June – October 2006. Results showed a wide range of soil hydraulic parameter values generated with the different methods, especially for the saturated hydraulic conductivity Ksat and the shape parameter α of the van Genuchten curve. This is reflected in a variability of the modeling results which is, as expected, different for each model and each variable analysed. The variability of the simulated water content in the root zone and of the bottom flux for different soil hydraulic parameter sets is found to be often larger than the difference between modeling results of the two models using the same soil hydraulic parameter set. Also we found that a good agreement in simulated soil moisture patterns may occur even if evapotranspiration and percolation fluxes are significantly different. Therefore multiple output variables should be considered to test the performances of methods and models.


2009 ◽  
Vol 6 (3) ◽  
pp. 4065-4105 ◽  
Author(s):  
G. Baroni ◽  
A. Facchi ◽  
C. Gandolfi ◽  
B. Ortuani ◽  
D. Horeschi ◽  
...  

Abstract. Data of soil hydraulic properties forms often a limiting factor in unsaturated zone modelling, especially at the larger scales. Investigations for the hydraulic characterization of soils are time-consuming and costly, and the accuracy of the results obtained by the different methodologies is still debated. However, we may wonder how the uncertainty in soil hydraulic parameters relates to the uncertainty of the selected modelling approach. We performed an intensive monitoring study during the cropping season of a 10 ha maize field in Northern Italy. These data were used to: i) compare different methods for determining soil hydraulic parameters and ii) evaluate the effect of the uncertainty in these parameters on different outputs (i.e. evapotranspiration, water content in the root zone, fluxes through the bottom boundary of the root zone) of two hydrological models with different complexity: SWAP, a widely used model of soil moisture dynamics in unsaturated soils based on Richards equation, and ALHyMUS, a conceptual model of the same dynamics based on a reservoir cascade scheme. We employed five direct and indirect methods to determine soil hydraulic parameters for each horizon of the experimental field. Two methods were based on a parameter optimization of: a) laboratory measured retention and hydraulic conductivity data and b) field measured retention and hydraulic conductivity data. Three methods were based on the application of widely used Pedo-Transfer Functions: c) Rawls and Brakensiek; d) HYPRES; and e) ROSETTA. Simulations were performed using meteorological, irrigation and crop data measured at the experimental site during the period June–October 2006. Results showed a wide range of soil hydraulic parameter values evaluated with the different methods, especially for the saturated hydraulic conductivity Ksat and the shape parameter α of the Van Genuchten curve. This is reflected in a variability of the modeling results which is, as expected, different for each model. The variability of the simulated water content in the root zone and of the fluxes at the root zone bottom for different soil hydraulic parameter sets is found to be often larger than the difference between modeling results of the two models using the same soil hydraulic parameter set. Also we found that a good agreement in simulated soil moisture patterns may occur even if evapotranspiration and percolation fluxes are significantly different. Therefore multiple output variables should be considered to test the performances of methods and models.


Water ◽  
2015 ◽  
Vol 7 (10) ◽  
pp. 5676-5688 ◽  
Author(s):  
Ken Okamoto ◽  
Kazuhito Sakai ◽  
Shinya Nakamura ◽  
Hiroyuki Cho ◽  
Tamotsu Nakandakari ◽  
...  

2007 ◽  
Vol 8 (4) ◽  
pp. 715-729 ◽  
Author(s):  
Binayak P. Mohanty ◽  
Jianting Zhu

Abstract In this study, the authors investigate effective soil hydraulic parameter averaging schemes for steady-state flow in heterogeneous shallow subsurfaces useful to land–atmosphere interaction modeling. “Effective” soil hydraulic parameters of the heterogeneous shallow subsurface are obtained by conceptualizing the soil as an equivalent homogeneous medium. It requires that the effective homogeneous soil discharges the same mean surface moisture flux (evaporation or infiltration) as the heterogeneous media. Using the simple Gardner unsaturated hydraulic conductivity function, the authors derive the effective value for the saturated hydraulic conductivity Ks or the shape factor α under various hydrologic scenarios and input hydraulic parameter statistics. Assuming one-dimensional vertical moisture movement in the shallow unsaturated soils, both scenarios of horizontal (across the surface landscape) and vertical (across the soil profile) heterogeneities are investigated. The effects of hydraulic parameter statistics, surface boundary conditions, domain scales, and fractal dimensions in case of nested soil hydraulic property structure are addressed. Results show that the effective parameters are dictated more by the α heterogeneity for the evaporation scenario and mainly by Ks variability for the infiltration scenario. Also, heterogeneity orientation (horizontal or vertical) of soil hydraulic parameters impacts the effective parameters. In general, an increase in both the fractal dimension and the domain scale enhances the heterogeneous effects of the parameter fields on the effective parameters. The impact of the domain scale on the effective hydraulic parameters is more significant as the fractal dimension increases.


Geoderma ◽  
2017 ◽  
Vol 285 ◽  
pp. 247-259 ◽  
Author(s):  
Andrea Sz. Kishné ◽  
Yohannes Tadesse Yimam ◽  
Cristine L.S. Morgan ◽  
Bright C. Dornblaser

2021 ◽  
Author(s):  
Tailin Li ◽  
Nina Noreika ◽  
Jakub Jeřábek ◽  
Tomáš Dostál ◽  
David Zumr

<p>A better understanding of hydrological processes in agricultural catchments is not only crucial to hydrologists but also helpful for local farmers. Therefore, we have built the freely-available web-based WALNUD dataset (Water in Agricultural Landscape – NUčice Database) for our experimental catchment Nučice (0.53 km<sup>2</sup>), the Czech Republic. We have included observed precipitation, air temperature, stream discharge, and soil moisture in the dataset. Furthermore, we have applied numerical modelling techniques to investigate the hydrological processes (e.g. soil moisture variability, water balance) at the experimental catchment using the dataset.</p><p>The Nučice catchment, established in 2011, serves for the observation of rainfall-runoff processes, soil erosion and water balance of the cultivated landscape. The average altitude is 401 m a.s.l., the mean land slope is 3.9 %, and the climate is humid continental (mean annual temperature 7.9 °C, average annual precipitation 630 mm). The catchment consists of three fields covering over 95 % of the area. There is a narrow stream which begins as a subsurface drainage pipe in the uppermost field draining the water at catchment. The typical crops are winter wheat, rapeseed, mustard and alfalfa. The installed equipment includes a standard meteorological station, several rain gauges distributed in the area of the basin, and an H flume to monitor the stream discharge, water turbidity and basic water quality indicators. The soil water content (at point scale) and groundwater level are also recorded. Recently, we have installed two cosmic-ray soil moisture sensors (StyX Neutronica) to estimate large-scale topsoil water content at the catchment.</p><p>Even though the soil management and soil properties in the fields of Nučice seem to be nearly homogeneous, we have observed variability in the topsoil moisture pattern. The method for the explanation of the soil water regime was the combination of the connectivity indices and numerical modelling. The soil moisture profiles from the point-scale sensors were processed in a 1-D physically-based soil water model (HYDRUS-1D) to optimize the soil hydraulic parameters. Further, the soil hydraulic parameters were used as input into a 3D spatially-distributed model, MIKE-SHE. The MIKE-SHE simulation has been mainly calibrated with rainfall-runoff observations. Meanwhile, the spatial patterns of the soil moisture were assessed from the simulation for both dry and wet catchment conditions. From the MIKE-SHE simulation, the optimized soil hydraulic parameters have improved the estimation of soil moisture dynamics and runoff generation. Also, the correlation between the observed and simulated soil moisture spatial patterns showed different behaviors during the dry and wet catchment conditions.</p><p>This study has been supported by the Grant Agency of the Czech Technical University in Prague, grant No. SGS20/156/OHK1/3T/11 and the Project SHui which is co-funded by the European Union Project: 773903 and the Chinese MOST.</p>


Sign in / Sign up

Export Citation Format

Share Document