Soil salinization is one of the significant concerns regarding irrigation with saline waters as an alternative resource for limited freshwater resources in arid and semi-arid regions. Thus, the investigation of proper management methods to control soil salinity for irrigation with saline waters is inevitable. The HYDRUS-1D model is a well-known numerical model that can facilitate the exploration of management scenarios to mitigate the consequences of irrigation with saline waters, especially soil salinization. However, before using the model as a decision support system, it is crucial to calibrate the model and analyze the model’s parameters and outputs’ uncertainty. Therefore, the generalized likelihood uncertainty estimation (GLUE) algorithm was implemented for the HYDRUS-1D model in the R environment to calibrate the model and assess the uncertainty aspects for simulating soil salinity of corn root zone under saline irrigation with linear move sprinkle irrigation system. The results of the study have detected a lower level of uncertainty in the α, n, and θs (saturated soil water content) parameters of water flow simulations, dispersivity (λ), and adsorption isotherm coefficient (Kd) parameters of solute transport simulations comparing to the other parameters. A higher level of uncertainty was found for the diffusion coefficient as its corresponding posterior distribution was not considerably changed from its prior distribution. The reason for this phenomenon could be the minor contribution of diffusion to the solute transport process in the soil compared with advection and hydrodynamic dispersion under saline water irrigation conditions. Predictive uncertainty results revealed a lower level of uncertainty in the model outputs for the initial growth stages of corn. The analysis of the predictive uncertainty band also declared that the uncertainty in the model parameters was the predominant source of uncertainty in the model outputs. In addition, the excellent performance of the calibrated model based on 50% quantiles of the posterior distributions of the model parameters was observed in terms of simulating soil water content (SWC) and electrical conductivity of soil water (ECsw) at the corn root zone. The ranges of NRMSE for SWC and ECsw simulations at different soil depths were 0.003 to 0.01 and 0.09 to 0.11, respectively. The results of this study have demonstrated the authenticity of the GLUE algorithm to seek uncertainty aspects and calibration of the HYDRUS-1D model to simulate the soil salinity at the corn root zone at field scale under a linear move irrigation system.

Utilizing degraded quality waters such as saline water as irrigation water with proper management methods such as leaching application is a potential answer to water scarcity in agricultural systems. Leaching application requires understanding the relationship between the amount of irrigation water and its quality with the dynamic of salts in the soil. The HYDRUS-1D model can simulate the dynamic of soil salinity under saline water irrigation conditions. However, these simulations are subject to uncertainty. A study was conducted to assess the uncertainty of the HYDRUS-1D model parameters and outputs to simulate the dynamic of salts under saline water irrigation conditions using the Markov Chain Monte Carlo (MCMC) based Metropolis-Hastings algorithm in the R-Studio environment. Results indicated a low level of uncertainty in parameters related to the advection term (water movement simulation) and water stress reduction function for root water uptake in the solute transport process. However, a higher level of uncertainty was detected for dispersivity and diffusivity parameters, possibly because of the study’s scale or some error in initial or boundary conditions. The model output (predictive) uncertainty showed a high uncertainty in dry periods compared to wet periods (under irrigation or rainfall). The uncertainty in model parameters was the primary source of total uncertainty in model predictions. The implementation of the Metropolis-Hastings algorithm for the HYDRUS-1D was able to conveniently estimate the residual water content (θr) value for the water simulation processes. The model’s performance in simulating soil water content and soil water electrical conductivity (ECsw) was good when tested with the 50% quantile of the posterior distribution of the parameters. Uncertainty assessment in this study revealed the effectiveness of the Metropolis-Hastings algorithm in exploring uncertainty aspects of the HYDRUS-1D model for reproducing soil salinity dynamics under saline water irrigation at a field scale.