The diversity of processes taking place in hillslope agro-ecosystems makes the estimation of vadose zone dynamics rather challenging. This paper presents the first insight into the research of ...volumetric water content, granulometric composition, meteorological data, precipitation and soil-water isotopic composition conducted within the SUPREHILL project at its vadose zone observatory. The main goals of this research are related to the evaluation of soil-water origin at the hillslope vineyard, but also to the estimation of depths until which precipitation infiltrates and where the occurrence of preferential flow is possible. For that purpose, hydrometeorological data, granulometric composition and stable isotopesof hydrogen (δ2H) and oxygen (δ18O) from precipitation and sampled soil water have been used. The results indicate the existence of a different isotopic signature in soil water, which suggests different infiltration patterns in the investigated area. Also, the results point out that surface runoff, subsurface runoff, and most of the passive wick lysimeters respond to precipitation, while the response of suction probes located at deeper depth is not that evident. This corresponds to the results related to the variation of water content at different depths. All the results indicate the possible existence of a low permeable layer at an approximate depth of 60 cm. Furthermore, preferential flow, if it exists, can be expected from the shallowest depths of the vineyard to a maximum depth of 80 cm. It is expected that an established long-term monitoring network at the SUPREHILL Observatory will give a more precise definition of soil-water behaviour and the existence of preferential flow.
Processes in hillslope soils present a particular challenge for agricultural production and soil management due to their hydropedological specifics and high soil erosion risk. Soil heterogeneities ...can cause preferential and/or lateral flow on the entire hillslope resulting in the off-site movement of water, fertilizers and chemicals used in crop production. A study was conducted under controlled conditions in a laboratory with undisturbed soil cores (250 cm3), which were used to estimate the soil hydraulic properties (SHP) using HYPROP and WP4C devices, while undisturbed soil columns (diameter = 16 cm, length = 25 cm) were used for the evaluation of preferential flow pathways using potassium bromide and Brilliant Blue. Samples were excavated in triplicate from the hilltop, backslope and footslope regions within the inter-rows of a vineyard from a critical zone observatory, SUPREHILL, in Croatia in Dystric Luvic Stagnosol. The aim of this study was to determine if the erosion-affected hillslope position affected the physical, chemical and hydraulic properties of soil and to identify water flow and possible preferential flow using dye and bromide tracers. The results of the sensor measurements and estimated SHPs were in agreement, showing a faster leaching of the irrigated rainwater in the footslope column. The tracer experiments showed variability even in the columns taken from the same position on the hillslope, which can be linked to plant roots and soil fauna activity. Altogether, the results showed a deeper loose layer at the footslope as a consequence of the soil erosion, which then resulted in higher hydraulic conductivity and the leached mass of the bromide due to better soil structure and pore connectivity. Thus, due to significant differences in the leached mass of bromide, this research should be later expanded in field experiments to reveal the impact of surface runoff, subsurface preferential and lateral flow on a larger scale.
This paper presents research related to the estimation of the precipitation fraction in the soil water of a sloped vineyard at the SUPREHILL Critical Zone Observatory (CZO) in Zagreb, Croatia. ...Numerous investigations have shown that exploration of hillslope soils can be very challenging due to the existence of heterogeneity and different soil properties, as well as due to anthropogenically induced processes, which can affect precipitation infiltration and soil water flow. Within this research, physicochemical soil properties, soil water content (SWC), and isotopic composition of soil water and precipitation (δ2H and δ18O) have been examined. The isotopic signature of soil water was monitored in 24 points, at 4 depths, throughout the hillslope vineyard. Soil water isotopic composition from all monitoring points coincided with the Local Meteoric Water Line (LMWL), with almost no variability at 100 cm depth, which was consistent with the smallest variation of SWC at 80 cm depth and indicated that most of water mixing takes place in the shallower part of the hillslope. Results suggested the existence of heterogeneity, uneven erosion processes in the footslope of the observed vineyard, and different infiltration patterns. Fractions of precipitation varied significantly depending on the depth and position in the vineyard, from approximately 1% up to 98%, where more precipitation fraction has been determined in the surface and subsurface runoff. Additionally, statistical analysis and a more detailed evaluation of precipitation fractions at the 40 cm depth, where wick lysimeters are installed, have shown that Corg content is related to the silt fraction, while the first results indicate that the infiltration patterns were dependent on the common influence of all observed physicochemical properties.
The negative environmental impact of conventional agriculture threatens agroecosystem stability and food security. Therefore, searching for optimal soil management practices is crucial for ...maintaining and improving soil functions. This work aims to determine the impact of conventional and conservation tillage on Stagnosols in a semi-humid environment in Marija Magdalena (Croatia) during 2021 and 2022. Under each tillage treatment, subplots were biochar, cattle manure, and control (split-plot design). The conservation tillage exhibits lower compaction in addition to conventional tillage. In 2021, at 0–15 cm and 15–30 cm depths, control plots had the highest bulk density (BD), while biochar plots had the lowest. In 2022, biochar and manure treatments under conventional tillage had significantly higher BD than those under conservation tillage. Penetration resistance did not exceed 2 MPa in all treatments. Soil water content was high in conservation treatments at 0–15 cm. Water-stable aggregates were higher in biochar and manure plots under both tillage treatments. Maize yield was higher in conservation treatments in 2021 and in conventional during 2022. Manure and biochar in the conventional system showed a better impact on grain yields than under conservation. Conservation tillage in rain-fed farming maintains crop yields and reduces soil compaction.
A quantitative understanding of actual evapotranspiration (ETa) and soil–water dynamics in a hillslope agroecosystem is vital for sustainable water resource management and soil conservation; however, ...the complexity of processes and conditions involving lateral subsurface flow (LSF) can be a limiting factor in the full comprehension of hillslope soil–water dynamics. The research was carried out at SUPREHILL CZO located on a hillslope agroecosystem (vineyard) over a period of two years (2021–2022) by combining soil characterization and field hydrological measurements, including weighing lysimeters, sensor measurements, and LSF collection system measurements. Lysimeters were placed on the hilltop and the footslope, both having a dynamic controlled bottom boundary, which corresponded to field pressure head measurements, to mimic field soil–water dynamics. Water balance components between the two positions on the slope were compared with the goal of identifying differences that might reveal hydrologically driven differences due to LSF paths across the hillslope. The usually considered limitations of these lysimeters, or the borders preventing LSF through the domain, acted as an aid within this installation setup, as the lack of LSF was compensated for through the pumping system at the footslope. The findings from lysimeters were compared with LSF collection system measurements. Weighing lysimeter data indicated that LSF controlled ETa rates. The results suggest that the onset of LSF contributes to the spatial crop productivity distribution in hillslopes. The present approach may be useful for investigating the impact of LSF on water balance components for similar hillslope sites and crops or other soil surface covers.
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.
The solubility and mobility of copper (Cu) in soil is strongly influenced by the presence of dissolved organic carbon (DOC); however, the interactions between Cu and DOC are complex and not yet fully ...understood. In this study, Cu and DOC concentrations were measured monthly for two years in leachates from self-constructed lysimeters installed at inter- and intra-row vineyard hilltop, backslope, and footslope areas at the SUPREHILL Critical Zone Observatory, Croatia. The aim was to quantify Cu and DOC leaching from the hilltop towards the backslope and the footslope. The assumed strong relationship between Cu and DOC in the leachates was statistically analyzed and explained using chemical equilibrium software. Leachates were analyzed for pH, EC, DOC, Cu, and major ion concentrations. The highest Cu concentrations found in leachates from the intra-row footslope suggested Cu downhill transport. Although not strong, a significant positive correlation between Cu and DOC in footslope leachates confirmed the relevance of Cu complexation by DOC. Speciation confirmed that more than 99.9% of total Cu in leachates was found as a Cu-DOC complex. Data implied the role of soil water flow pathways in explaining Cu downhill transport. Critical timing for applying Cu fungicides at sloped vineyards was highlighted.
Soil and water loss due to traditional intensive types of agricultural management is widespread and unsustainable in Croatian croplands. In order to mitigate the accelerated land degradation, we ...studied different cropland soil management strategies to obtain feasible and sustainable agro-technical practices. A rainfall simulation experiment was conducted at 58 mm h–1 over 30 min on 10 paired plots (0.785 m2), bare and straw covered (2 t ha−1). The experiment was carried out in maize cultivation (Blagorodovac, Croatia) established on Stagnosols on slopes. Measurements were conducted during April (bare soil, after seeding), May (five-leaves stage), and June (intensive vegetative growth) making 60 rainfall simulations in total. Straw reduced soil and water losses significantly. The highest water, sediment loss, and sediment concentrations were identified in tillage plots during May. Straw addition resulted in delayed ponding (for 7%, 63%, and 50% during April, May and June, respectively) and runoff generation (for 37%, 32%, and 18% during April, May and June, respectively). Compared with the straw-mulched plot, tillage and bare soil increased water loss by 349%. Maize development reduced the difference between bare and straw-mulched plots. During May and June, bare plots increase water loss by 92% and 95%, respectively. The straw mulch reduced raindrop kinetic energy and sediment detachment from 9, 6, and 5 magnitude orders in April, May, and June, respectively. Overall, the straw mulch was revealed to be a highly efficient nature-based solution for soil conservation and maize cultivation protection.
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.
Nutrient leaching from agricultural soils presents an economic loss for farmers and can degrade the quality of the surrounding environment. Thus, leachates from 18 in situ wick lysimeters, installed ...at 40 cm soil depth at the vineyard hilltop, backslope, and footslope intra- and inter-row area (SUPREHILL Critical Zone Observatory, Croatia) were collected monthly over two years and analyzed for major plant nutrient ions. Our objectives were to quantify nutrient losses via leaching from the hilltop towards the backslope and to the footslope, and to compare leaching from vine plant rows (intra-row) with grassed areas between vine rows (inter-row). We found that the concentrations of nitrate, orthophosphate, and potassium were significantly higher in leachates collected at the footslope as compared to the hilltop and backslope only at intra- and not at inter-row positions, while ammonium was independent of the slope and row positions. The vineyard intra-row is identified as the probable spatial origin of nutrient leaching along the slope, thus confirming spatially different contributions of overall hillslope to major plant nutrients leaching. The experimental field scheme used in this study, which separately analyses vineyard intra- and inter-row, was confirmed to be an adequate approach for optimizing vineyard management practices.