Modeling seasonal soil moisture dynamics in gley soils in relation to groundwater table oscillations in eastern Croatia

This study aims to explain complex vadose zone hydrology of fine-textured (gley) agricultural soils influenced by a shallow and dynamic groundwater (GW) levels. The field site was located in the Biđ field (Eastern Croatia), where a detailed soil survey was performed. The simulations included a three-year period (2016–2018) at four locations. Soil hydraulic parameters (SHP) were estimated based on variables determined in the laboratory, while soil water flow was monitored using in-field zero-tension lysimeters. Piezometers were installed and used to monitor daily oscillations of groundwater levels (average depth to GW 2.2 m), while data from nearby Sava River was monitored. Unsaturated flow and water regime assessment was performed using HYDRUS-1D numerical modeling. Additional SHP optimization of van Genuchten-Mualem parameters (α and n) was performed using Shuffled Complex Evolution algorithm (SCE). The autocorrelation analysis was used to detect patterns in the precipitation, GW, and river level time series, while the Mutual Information (MI) was used to estimate the codependence of the processes in unsaturated zone and the main hydrological events. The model successfully (R2 0.72 – 0.94) reproduced measured lysimeters outflows. The outflows from lysimeters were connected to precipitation patterns, transpiration intensity, and soil moisture content influenced by the shallow water table. Comparable MI values obtained for precipitation, GW, and river level suggest a concurrent role of these parameters in the unsaturated flow dynamics. The relationship between upward flux/water storage change into the domain, and transpiration/growth stages, suggests a strong connection between the water fluxes and the root water uptake. Results confirm the importance of GW for the agricultural production due to the major influence on upper soil layer moisture.