Determination of soil hydraulic parameters and evaluation of water dynamics and nitrate leaching in the unsaturated layered zone: A modeling case study in Central Croatia
- Defterdarović, Jasmina, Filipović, Lana, Kranjčec, Filip, Ondrašek, Gabrijel, Kikić, Diana, Novosel, Alen, Mustać, Ivan, Krevh, Vedran, Magdić, Ivan, Rubinić, Vedran, Bogunović, Igor, Dugan, Ivan, Čopec, Krešimir, He, Hailong, Filipović, Vilim
- Authors: Defterdarović, Jasmina , Filipović, Lana , Kranjčec, Filip , Ondrašek, Gabrijel , Kikić, Diana , Novosel, Alen , Mustać, Ivan , Krevh, Vedran , Magdić, Ivan , Rubinić, Vedran , Bogunović, Igor , Dugan, Ivan , Čopec, Krešimir , He, Hailong , Filipović, Vilim
- Date: 2021
- Type: Journal article
- Relation: Sustainability (Basel, Switzerland) Vol. 13, no. 12 (2021), p. 6688
- Full Text:
- Reviewed:
- Description: 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.
- Authors: Defterdarović, Jasmina , Filipović, Lana , Kranjčec, Filip , Ondrašek, Gabrijel , Kikić, Diana , Novosel, Alen , Mustać, Ivan , Krevh, Vedran , Magdić, Ivan , Rubinić, Vedran , Bogunović, Igor , Dugan, Ivan , Čopec, Krešimir , He, Hailong , Filipović, Vilim
- Date: 2021
- Type: Journal article
- Relation: Sustainability (Basel, Switzerland) Vol. 13, no. 12 (2021), p. 6688
- Full Text:
- Reviewed:
- Description: 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.
Effect of slope position on soil properties and soil moisture regime of Stagnosol in the vineyard
- Magdić, Ivan, Safner, Toni, Rubinić, Vedran, Rutić, Filip, Husnjak, Stjepan, Filipović, Vilim
- Authors: Magdić, Ivan , Safner, Toni , Rubinić, Vedran , Rutić, Filip , Husnjak, Stjepan , Filipović, Vilim
- Date: 2022
- Type: Text , Journal article
- Relation: Journal of Hydrology and Hydromechanics Vol. 70, no. 1 (2022), p. 62-73
- Full Text:
- Reviewed:
- Description: 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 ( < 0.0001), with lower confidence level in bulk density < 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 ( < 0.0001), and between row and inter-row vineyard space (all with < 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.
- Authors: Magdić, Ivan , Safner, Toni , Rubinić, Vedran , Rutić, Filip , Husnjak, Stjepan , Filipović, Vilim
- Date: 2022
- Type: Text , Journal article
- Relation: Journal of Hydrology and Hydromechanics Vol. 70, no. 1 (2022), p. 62-73
- Full Text:
- Reviewed:
- Description: 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 ( < 0.0001), with lower confidence level in bulk density < 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 ( < 0.0001), and between row and inter-row vineyard space (all with < 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.
- Filipović, Vilim, Defterdarović, Jasmina, Krevh, Vedran, Filipović, Lana, Ondrašek, Gabrijel, Kranjčec, Filip, Magdić, Ivan, Rubinić, Vedran, Stipičević, Sanja, Mustać, Ivan, Bubalo Kovačić, Marina, He, Hailong, Haghverdi, Amir, Gerke, Horst
- Authors: Filipović, Vilim , Defterdarović, Jasmina , Krevh, Vedran , Filipović, Lana , Ondrašek, Gabrijel , Kranjčec, Filip , Magdić, Ivan , Rubinić, Vedran , Stipičević, Sanja , Mustać, Ivan , Bubalo Kovačić, Marina , He, Hailong , Haghverdi, Amir , Gerke, Horst
- Date: 2021
- Type: Text , Journal article
- Relation: Agronomy (Basel) Vol. 12, no. 1 (2021), p. 33
- Full Text: false
- Reviewed:
- Description: 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 methodssoil 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.
Investigation of hillslope vineyard soil water dynamics using field measurements and numerical modeling
- Krevh, Vedran, Groh, Jannis, Weihermüller, Lutz, Filipović, Lana, Defterdarović, Jasmina, Kovač, Zoran, Magdić, Ivan, Lazarević, Boris, Baumgartl, Thomas, Filipović, Vilim
- Authors: Krevh, Vedran , Groh, Jannis , Weihermüller, Lutz , Filipović, Lana , Defterdarović, Jasmina , Kovač, Zoran , Magdić, Ivan , Lazarević, Boris , Baumgartl, Thomas , Filipović, Vilim
- Date: 2023
- Type: Text , Journal article
- Relation: Water (Switzerland) Vol. 15, no. 4 (2023), p.
- Full Text:
- Reviewed:
- Description: 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. © 2023 by the authors.
- Authors: Krevh, Vedran , Groh, Jannis , Weihermüller, Lutz , Filipović, Lana , Defterdarović, Jasmina , Kovač, Zoran , Magdić, Ivan , Lazarević, Boris , Baumgartl, Thomas , Filipović, Vilim
- Date: 2023
- Type: Text , Journal article
- Relation: Water (Switzerland) Vol. 15, no. 4 (2023), p.
- Full Text:
- Reviewed:
- Description: 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. © 2023 by the authors.
- «
- ‹
- 1
- ›
- »