Hillslope hydrology in agricultural landscapes is complex due to a variety of hydropedological processes and field management possibilities. The aim was to test if there are any differences in soil properties and water regime along the hillslope and to compare vineyard rows (vine) with inter-rows (grass) area for those properties. The study determined that there are significant differences in the contents of soil particle fractions, pH, and humus content along the slope (P < 0.0001), with lower confidence level in bulk density (P < 0.05). Differences between row and inter-row space were significant for the pH, humus, and silt content, but for sand and clay content, and bulk density differences were not determined. The study determined differences in soil water content among five slope positions (P < 0.0001), and between row and inter-row vineyard space (all with P < 0.05). Where in the upper slope positions (e. g., P1) soil water content was higher than on lower slope positions. Higher soil water content was observed at higher slope positions, associated with clay content. However, it can be concluded that the retention of moisture on the slope is more influenced by local-scale soil properties (primarily soil texture) and variability of the crop (row/inter-row) than the position on the slope.

In urban areas, land use usually increases soil degradation. However, there are areas occupied by agriculture and woodlands with an essential role in provisioning food and other services such as water and climate regulation. The objective of this work was to assess the effect of long-term land use and soil management practices on peri-urban soils in Zagreb (Croatia). Samples were collected at depth 0–10 cm within intensively tilled cropland (CROP) and vineyard (VINE), traditional grass-covered orchard (ORCH), and forest (FOR). The results showed that bulk density was significantly higher in VINE and CROP than in ORCH and FOR. The opposite dynamic was observed in water-holding capacity, air-filled porosity, aggregate stability, organic matter, and soil organic matter stocks (SOCS). Soil water infiltration was higher in FOR plot compared to the other plots. Overall, land-use change had a substantial impact on soil properties and SOCS, especially in CROP and VINE soils. Tillage, pesticides, and fertilizer applications were presumably the reasons for altered soil quality properties. Intensively used areas (VINE and CROPS) may reduce soil ecosystems services such as the capacity for flood retention and C sequestration.

Soil heterogeneities can impact hillslope hydropedological processes (e.g., portioning between infiltration and runoff), creating a need for in-depth knowledge of processes governing water dynamics and redistribution. The presented study was conducted at the SUPREHILL Critical Zone Observatory (CZO) (hillslope vineyard) in 2021. A combination of field investigation (soil sampling and monitoring campaign) and numerical modeling with hydrological simulator HYDRUS-1D was used to explore the water dynamics in conjunction with data from a sensor network (soil water content (SWC) and soil-water potential (SWP) sensors), along the hillslope (hilltop, backslope, and footslope). Soil hydraulic properties (SHP) were estimated based on (i) pedotransfer functions (PTFs), (ii) undisturbed soil cores, and (iii) sensor network data, and tested in HYDRUS. Additionally, a model ensemble mean from HYDRUS simulations was calculated with PTFs. The highest agreement of simulated with observed SWC for 40 cm soil depth was found with the combination of laboratory and field data, with the lowest average MAE, RMSE and MAPE (0.02, 0.02, and 5.34%, respectively), and highest average R² (0.93), while at 80 cm soil depth, PTF model ensemble performed better (MAE = 0.03, RMSE = 0.03, MAPE = 7.55%, R² = 0.81) than other datasets. Field observations indicated that heterogeneity and spatial variability regarding soil parameters were present at the site. Over the hillslope, SWC acted in a heterogeneous manner, which was most pronounced during soil rewetting. Model results suggested that the incorporation of field data expands model performance and that the PTF model ensemble is a feasible option in the absence of laboratory data.

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 cm³ cm⁻³ while the hydraulic conductivity varied from 3 cm day⁻¹ to 990.9 cm day⁻¹. 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.

Erosion has been reported as one of the top degradation processes that negatively affect agricultural soils. The study objective was to identify hydropedological factors controlling soil water dynamics in erosion-affected hillslope vineyard soils. The hydropedological study was conducted at identically-managed Jastrebarsko (location I), and Jazbina (II) and (III) sites with Stagnosol soils. Soil Hydraulic Properties (SHP) were estimated on intact soil cores using Evaporation and WP4C methods; soil hydraulic functions were fitted using HYPROP-FIT software. For Apg and Bg/Btg horizons, uni- and bimodal soil hydraulic models could be well fitted to data; although, the bimodal model performed better in particular cases where data indicated non-uniform pore size distribution. With these SHP estimations, a one-year (2020) water flow scenario was simulated using HYDRUS-1D to compare water balance results obtained with uni- and bimodal hydraulic functions. Simulation results revealed relatively similar flux distribution at each hillslope position between the water balance components infiltration, surface runoff, and drainage. However, at the bottom profile at Jastrebarsko, bimodality of the hydraulic functions led to increased drainage. Soil water storage was reduced, and the vertical movement increased due to modified soil water retention curve shapes. Adequate parameterization of SHP is required to capture the hydropedological response of heterogenous erosion-affected soil systems.

Variations in soil pH have been shown to affect mesotrione adsorption, which in turn, may have an impact on crop susceptibility. Therefore, a greenhouse experiment was conducted to evaluate the effect of simulated mesotrione residues on pea crop grown in the typical agricultural soil (gleysol) of north-western Croatia. The soil pH was manipulated to obtain neutral (pH 7.0) and acidic (pH 5.0) values. Simulated mesotrione residues were 1.1, 2.3, 4.5, 9.0, 18, 36 and 72 g a.i. ha⁻¹. Crop visual injuries as well as reductions in chlorophyll fluorescence and aboveground dry biomass were higher at pH 7.0 than at pH 5.0. With increasing mesotrione residues, the reductions in chlorophyll fluorescence ranged from 38.8% to 89.7% at pH 5.0 and from 63.7% to 99.3% at pH 7.0. Compared to chlorophyll fluorescence, the reductions in dry biomass were smaller and ranged from 49.2% to 96.8% at pH 7.0 and from 32.0% to 82.6% at pH 5.0 for the mesotrione residues from 1.1 to 72 g a.i. ha⁻¹. These results indicate that soil pH is an important factor determining the susceptibility of pea crop to simulated mesotrione residues.