•Proper assessment of the soil structure of arable soils requires a combination of various methods.•Results of evaporation and leaching experiments were fitted using a single-porosity flow ...model.•Data collected during flow experiments did not provide evidence of non-equilibrium flow.•Additional dye-staining and X-ray imaging revealed a complex biopore network.•The bi-modal flow behavior was revealed using the dual-permeability modeling of dye transport.
The appearance and distribution of soil pores have a significant influence on water flow and solute transport in the soil vadose zone. The pore system is highly variable in arable soils where crop rotation, tillage, trafficking, soil amendments, and various management practices are commonly implemented. The aim of this study was to assess the porous system and preferential flow pathways in a vineyard soil using undisturbed soil columns, and by combining laboratory and numerical methods with dye staining and X-ray imaging. It was hypothesized that the integration of various methods could reveal more information about soil structure, and flow and transport behavior of structured arable soil. Soil water retention and hydraulic conductivity curves were obtained using the evaporation method, while water flow was assessed using intermittent leaching experiments. Water flow and the transport of Brilliant Blue were simulated using HYDRUS-1D. A single-porosity model of soil hydraulic properties provided a good description of data collected during the evaporation experiments. Data collected during leaching experiments did not provide enough experimental evidence for the occurrence of nonequilibrium flow patterns and the differentiation between the single- and dual-permeability models of soil hydraulic properties. However, dye staining and X-ray imaging revealed a complex pore-architecture network with large vertical and horizontal biopores. The staining patterns (Brilliant Blue FCF) within the vertical column sections documented the extent of preferential flow. The study showed that the bi-modal character of pore structure could often be hidden when a limited number or non-adequate methods are applied for its quantification from water flow behavior. The impact of preferential pathways on dye transport can be investigated with observations and simulations. A combination of various methods enabled us to adequately assess vineyard soil structure and fine-tune the description and extent of preferential water 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.
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.
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 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.
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.
Intensive agricultural practices increase agrochemical pollution, particularly nitrogen (N) based fertilizers, which present an environmental risk. This study aims to evaluate long-term (2009–2020) ...data on soil water regime and nitrate dynamics at an agricultural experimental site on fine-textured soils and to better understand the implications of N management in relation to groundwater pollution. The field site is located in the Biđ field (eastern Croatia), in the proximity of the Sava river. Zero-tension lysimeters were installed at six selected locations. Lysimeters were used to monitor the water regime, i.e., outflows in which nitrate concentration was measured, while additional soil-water samples were collected via 4 and 15-meter-deep monitoring wells. Soil hydraulic parameters were estimated by combining the laboratory measurements, and estimation in RETC software. Water regime and nitrate leaching in lysimeters were simulated using HYDRUS-1D for each year to allow crop rotation and to evaluate their effects individually. The HYDRUS-1D model successfully reproduced lysimeter outflows and nitrate dynamics, which was confirmed with high R2 values (water: 93% above 0.7, and nitrate: 73% above 0.7) indicating the good performance of the model simulating nitrification chain reactions. Principal component analysis (PCA) was performed to identify the relationships among all soil properties and environmental characteristics. The results showed the complex interaction of soil hydraulic properties, precipitation patterns, plant uptake, and N application. All locations have a decreasing trend of nitrate leaching over the investigation period. Most of the lysimeter outflows and elevated nitrate concentrations were connected to the wet period of the year when the soil was saturated, and evapotranspiration was low. The results of this study show that it is important to optimize N fertilizer applications for each particular environmental condition to reduce nitrate loss. The study indicates the importance of long-term field studies, key for agro-hydrological modeling and the improvement of agricultural practices.
Display omitted
•Interaction of soil hydraulic parameters, rainfall, crops, and nitrogen is presented.•HYDRUS-1D successfully reproduced long-term lysimeter outflows and nitrate dynamics.•Field conditions and evapotranspiration affected outflows and nitrate concentrations.•Long-term field data is key for modeling and improvement of agricultural 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.
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.