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.

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.

Prostorno modeliranje radi opisivanja prostorne varijabilnosti fizikalnih svojstava tla u istočnoj Hrvatskoj Bogunović, I., Đekemati, I., Magdić, I., Vrbanić, M., Matošić, S., Mesić, M. Poljoprivreda

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−1. 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−1. These results indicate that soil pH is an important factor determining the susceptibility of pea crop to simulated mesotrione residues.

Quantification of soil plasticity is usually based on Atterberg limits or indices, which are then used for engineering and agricultural purposes on clay soils. Because these limits/indices are seldom analyzed during routine soil surveys, they are sometimes estimated from available soil properties using pedotransfer functions (PTFs). Main aim of this study was to test if two robust PTFs, previously created by other authors for other soils/areas, may be used for predicting plasticity of typical clay soils on fluvial deposits in Croatia. We analyzed mineralogy, plasticity, and related physicochemical properties of four representative soil profiles along Sava and Drava, two longest Croatian rivers. Particle size distribution patterns pointed to more uniform sedimentation along Sava, compared to Drava. Also, more clay was found within Sava profiles. Soil texture was finer farther away from the sources of both rivers. Soil cation exchange capacity (CEC) was almost fully positively correlated to clay content. On the other hand, C org content showed no correlation to CEC. This is attributed to the significant presence of smectite across all studied soils. Clay was recognized as the main factor influencing soil plasticity. When compared to the measured values, the predicted values of plasticity index and liquid limit were heavily underestimated. Hence, region-specific PTFs should be developed for more accurate prediction of plasticity in these soils. According to the WRB-2015, the soils were classified as Eutric Reductigleyic Stagnic Gleysols (Clayic, Humic, Protovertic). Because soils were not dry during field description, shrink-swell cracks were not prominent, and therefore soils were not classified as Vertisols. We suggest that field criteria for classification of Vertisols should not depend on actual soil moisture. Instead, plastic limits/indices could be used.

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 R2 (0.93), while at 80 cm soil depth, PTF model ensemble performed better (MAE = 0.03, RMSE = 0.03, MAPE = 7.55%, R2 = 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.