Search

Uslijed sve izraženijih klimatskih promjena, na području istočne Hrvatske dolazi do sve češćih i duljih sušnih razdoblja što rezultira promjenama u vodnom režimu područja te poteškoćama u gospodarenju tlom i vodom za potrebe poljoprivrede. Iz tog razloga javlja se sve češća potreba za izgradnjom sustava navodnjavanja u svrhu stabilne i isplative poljoprivredne proizvodnje. Stoga ovaj rad ima za cilj odrediti pogodnost za navodnjavanje dreniranog poljoprivrednog zemljišta (tala) na području istočne Hrvatske. Za analizu osnovnih klimatskih pokazatelja korišteni su podaci s meteorološke postaje (m.p) Osijek-Čepin i m.p Gradište kod Županje, za vremensko razdoblje od 2003. do 2022. godine. Utvrđen je blagi trend sniženja godišnjih količina oborina te porast mjesečnih temperatura zraka od 0,09 °C godišnje. Također utvrđena je visoka varijabilnost oborina te kontinuirani porast nedostatka vode u tlu i smanjenje viška vode u tlu. Na istraživanom dreniranom poljoprivrednom zemljištu istočne Hrvatske (ukupne površine 74.871,2 ha), obzirom na klase pogodnosti za potrebe navodnjavanja, najviše je zastupljena klasa P-2 (63,2 %), a zatim slijede klasa P-3 (26,0 %) i klasa P-1 (9,8 %). Svega 1 % istraživanih površina svrstan je u klasu nepogodnosti N-1, a klasa trajno nepogodnih tala N-2 nije utvrđena.

Major losses of agricultural production and soils are caused by erosion, which is especially pronounced on hillslopes due to specific hydrological processes and heterogeneity. Therefore, the aim of this study was to assess the impact of agricultural management on the compaction, infiltration, and seasonal water content dynamics of the hillslope. Measurements were made at the hilltop and footslope, i.e., soil water content and potential were measured using sensors, wick lysimeters were used to quantify water flux, while a mini-disk infiltrometer was used to measure the infiltration rate and calculate the unsaturated hydraulic conductivity (K_unsat). Soil texture showed differences between hillslope positions, i.e., at the hilltop after 50 cm depth, the soil is classified as silty clay loam, and from 75 cm onward, the soil is silty clay, while at the footslope, the soil is silt loam even at the deeper depths. The results show a higher K_unsat at the footslope as well as higher average water volumes collected in wick lysimeters compared to the hilltop. Average water volumes showed a statistically significant difference at p < 0.01 between the hilltop and the footslope. The soil water content and water potential sensors showed higher values at the footslope at all depths, i.e., 8.0% at 15 cm, 8.4% at 30 cm, and 27.3% at 45 cm. The results show that, even though the vineyard is located in a relatively small area, soil heterogeneity is present, affecting the water flow along the hillslope. This suggests the importance of observing water movement in the soil, especially today when facing extreme weather (e.g., short-term high-intensity rainfall events) in order to protect soil and water resources.

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.

The intensification of climate change in the area of Eastern Croatia is leading to more frequent and longer drought periods, resulting in changes in the water regime of the area and difficulties in managing soil and water for agricultural purposes. For this reason, there is an increasing need for the construction of irrigation systems for stable and profitable agricultural production. Therefore, this study aims to determine the suitability for irrigation of drained agricultural land in the Eastern Croatia area. Data from the meteorological stations Osijek-Čepin and Gradište kod Županje were used to analyze basic climate indicators for the period from 2003 to 2022. A slight trend of decreasing annual precipitation and increasing monthly air temperatures by 0.09°C per year was identified. High variability of precipitation and continuous increase in water deficit and decrease in water surplus in soil were also determined. Regarding the suitability classes for irrigation needs, on the researched drained agricultural land in Eastern Croatia (with a total area of 74,871.2 ha), the P-2 class is the most represented (63.2%), followed by the P-3 class (26.0%) and the P-1 class (9.8%). Only 1 % of the researched area were classified as unsuitable (N-1 class), and the permanently unsuitable N-2 class was not identified.

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.

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.

Organic agriculture (OA) is a continuously growing global concept that emphasizes the use of sustainable and environmentally-friendly practices. By adopting OA, it is possible to improve ecosystems services, increase biodiversity, decrease environmental pollution, reduce carbon footprints and mitigate greenhouse gas emissions, generating food that is free from harmful residues of agrochemicals, thereby enhancing food safety and security. This study provides a comprehensive review of the latest insights on the global utilization of land resources in OA, focusing particularly on some EU countries that experienced a notable and rapid progress in organic farming during the past two decades. With negligible 1.6% (75 Mha) of global cropland is currently dedicated to OA, there is ample opportunity to expand the adoption of OA and realize its multi-beneficial potential for farmers (by premium prices), and consumers (by healthier and nutritious food). The importance of OA has been recognized by the most recent EU agro-environmental policies and green strategies, with an ambitious goal to have at least 25% of agroecosystems under organic management by 2030. Despite numerous financial supports and a multifold increase in OA land area, many member states are unlikely to achieve this goal, including Croatia, which currently has a share of only 8% (~109,000 ha) of lands in OA. Furthermore, converting conventional land to organic farming has not always led to an increase in value-added final OA products. EU policies related to OA have been focused on the area of land cultivated organically, rather than overall production performance, and financial subsidies have been essential to achieve this policy. Therefore, some of critical obstacles and challenges for OA under rising pressures due to global climate change, public health and geopolitical crises need to be managed by specifically designed policies and regulations, which would contribute to more sustainable OA, i.e., food safety and security.

The interplay of surface and shallow subsurface fluxes plays a critical role in controlling water movement in hillslope agroecosystems and impacting soil and plant health during prolonged dry periods, demonstrating a need for in-field monitoring. This study was conducted for two years (2021–2022) by combining field monitoring of the grass-covered inter-row area (passive wick lysimeter, surface runoff, and meteorological data), laboratory determination of soil hydraulic properties (SHPs), and numerical modeling with the aim to explore near-surface fluxes at the SUPREHILL Critical Zone Observatory (CZO) located on a hillslope vineyard. Additionally, sensitivity analysis for basic root water uptake (RWU) parameters was conducted. The model was evaluated (R2, RMSE, and NSE) with lysimeter (hillslope) and runoff (footslope) data, producing good agreement, but only after the inverse optimization of laboratory estimated hydraulic conductivity was conducted, demonstrating that adequate parameterization is required to capture the hydropedological response of erosion-affected soil systems. Results exhibit the dependence of runoff generation on hydraulic conductivity, rainfall, and soil moisture conditions. The data suggest different soil-rewetting scenarios based on temporal rainfall variability. Sensitivity analysis demonstrated that Leaf Area Index (LAI) was the most responsive parameter determining the RWU. The study offers an approach for the investigation of fluxes in the topsoil for similar sites and/or crops (and covers), presenting the methodology of self-constructed soil–water collection instruments.

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.

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.

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.

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 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.

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.

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.

Straw mulching is a sustainable practice used to control soil erosion. However, different doses of mulch affect the efficiency of straw conservation. This study presents detailed research on how soil physicochemical properties and the hydrological response react to different types of vineyard soil management (Tilled, Grass, Low Straw, High Straw) and seasons (spring, summer, autumn) under conventional management on Anthrosols in Mediterranean conditions. To assess soil properties, core samples and disturbed samples were taken from the topsoil layer (0–10 cm). To evaluate erosion rates, a rainfall simulation experiment was conducted (58 mm h−1 for 30 min) with 10 replicates per treatment and season (120 in total). The results show higher water-stable aggregates (WSA) and soil organic matter (SOM) and lower bulk density (BD) in the mulch and grass treatment groups compared with the Tilled treatment group. High Straw treatment successfully mitigated runoff, while other treatments had significantly higher runoff that triggered sediment loss (SL) and translocation of P, K, Zn and Ni down the slope. There were 254% and 520% higher K losses with Tilled treatment in autumn compared with Low Straw and Grass treatments, respectively. Statistical analysis showed a strong association between element loss and SL, which indicates an ecological threat in degraded and endangered vineyards. Mulch application and grass cover reduce the vulnerability of vineyards, reduce evaporation, act as insulation against high temperatures, reduce erosion and suppress weed growth. The mulch dosage varies depending on the goals and conditions of the vineyard; thus, lower mulch dosage (2 t/ha) is appropriate when soil conditions are favourable and there is no significant need for moisture retention, while higher mulch dosage is necessary in dry regions to maintain soil moisture during high-temperature periods, as well as in sloped areas subjected to erosion.