Open Access
How to translate text using browser tools
12 June 2020 Experimental Comparison of Runoff Generation and Initial Soil Erosion Between Vineyards and Croplands of Eastern Croatia: A Case Study
Igor Bogunovic, Leon Josip Telak, Paulo Pereira
Author Affiliations +

Soil and water loss in agricultural fields is a global problem. Although studies about soil erosion in croplands and vineyards exist, the direct comparison between these land uses is missing, especially under continental climates in Europe. Therefore, it is needed to find control measures to the impacts of these land-use management strategies on soil properties and hydrological response. The objective of this work is to estimate and compare the impacts of croplands and vineyards under conventional management croplands and vineyards on soil properties (water holding capacity—WHC; bulk density—BD; soil water content—SWC; water stable aggregates—WSA; mean weight diameter—MWD; soil organic matter—SOM; available phosphorus—AP; total nitrogen—TN) and hydrological response (runoff—Run; sediment content—SC; sediment loss—SL; carbon loss—C loss; phosphorus loss—P loss; nitrogen loss—N loss) in Eastern Croatia. To achieve these goals, a study was set up using rainfall simulation tests at 58 mm h-1 over 30 minutes on 2 locations (Zmajevac: 45°48′N; 18°46′E; Erdut: 45°30′N; 19°01′E). In total, 32 rainfall simulations were carried out, 8 repetitions in vineyards and 8 in cropland plots of 0.876 m2, per location. Bulk density was significantly higher in cropland plots compared with the vineyard. Soil water content was significantly higher in Zmajevac cropland compared with Erdut plots. Also, SWC was significantly lower in Zmajevac vineyard than in the cropland located in the same area. Water stable aggregates and MWD were significantly higher in vineyard plots than in the cropland. Also, SOM and TN were significantly lower in Zmajevac cropland compared with the vineyard located in the same area. Available phosphorus was significantly high in Zmajevac plots than in Erdut. The rainfall simulations showed that Run was significantly higher in Erdut vineyard (8.2 L m-2) compared with Zmajevac (3.8 L m-2). Also, the Run in Erdut Cropland was significantly lower than in the vineyard. Sediment content did not show significant differences among locations. In Erdut, vineyard plots had a significantly lower SL (28.0 g m-2) than the cropland ones (39.1 g m-2). C loss was significantly higher in Zmajevac cropland than in Erdut. Also, C loss was significantly lower in Zmajevac vineyard compared with the cropland. We did not observe significant differences in P loss, and N loss also did not show significant differences. The principal component analysis showed that SOM was associated with WSA, AP, and TN. These variables were negatively related to slope, SWC, and C loss (factor 1). Also, MWD was inversely related to SL, P, and N loss (factor 2). Bulk density and SC were negatively related to Run. Overall, we conclude that noninvertive tillage practices in vineyards preserve soil structure, enhance soil quality, and reduce the extent of soil degradation.


Erosion by water is considered one of the major threats to soil resources.1 Although it is directly connected to the climatic, pedological, edaphic, and geomorphological conditions, the agriculture practices (tillage) used increase exponentially soil erosion rates2,3 In this context, conventional tillage, involving a sequence of plowing, disking, and harrowing, and agrochemicals are recognized as primary drivers of unsustainable soil erosion rates.4,5 Soil erosion rates can be 10 to 40 times higher than the rate of soil formation.6

Studies focused on soil erosion in croplands dominate in the international literature.7 In the last 2 decades, there has been a high increase in studies conducted in vineyards, highlighting the importance of this land-use type in soil erosion works.8,9 Several studies were carried out in the Mediterranean and reported soil erosion rates up to 14 t ha-1 y-1.10,11 Numerous studies also reported a positive impact of conservation strategies to mitigate soil erosion problems.8,9,12 Nowadays, despite a few exceptions, most of the studies conducted to study soil erosion were focused on 1 land-use type.13,14 Therefore, it is necessary to determine the impact of land use and the respective managements on soil erosion and identify which land uses and practices are more detrimental to the soil. This could be beneficial to identify areas and propose measures to achieve better sustainable management and land degradation neutrality.15

Some authors confirmed that in croplands, soil erosion by water is influenced by tillage intensity, slope direction, crop, planting direction, and/or orientation.16,17 In plowed fields, unsustainable soil erosion rates were observed in per-humid,18 semi-humid,19 arid,20 or semi-arid21 environments on sandy,22 silty,23 or clay soils.24 Similarly, soils in vineyards are highly sensitive to management, which affects their hydrological response dramatically. Vineyard management involves frequent tillage and tractor traffic. These practices reduce soil structure and hydraulic properties25,26 and increase soil compaction, runoff, and soil erosion above the tolerant levels.2 Previous works revealed that soil erosion in vineyards increases with the increasing slope,8,10 tillage intensity,9 and soil compaction.26 On the contrary, a decrease in vegetation cover increases it.2 In croplands, soil erosion increases with tillage intensity21,27 and slope20 and decreases with a plant density of cover crops.16 Several works compared the impacts of different land uses and management on soil properties or hydrological response in Europe,7,11,28-38 Africa,39-41 Asia,20,42-50 North America,51-54 and South America.55-57 Some studied the impact of land use on soil properties and hydrological response.12,28-30,35,37,48,50,52,53,55,57 From these works, none has compared croplands and vineyards. To our knowledge, 1 study investigated cropland and vineyard soil management on soil properties and hydrological response, but in the Mediterranean environment and a different type of soil (terra rossa).28 Therefore, studies at the pedon-scale are needed in croplands and vineyards to understand the impacts of land uses and the respective practices on soil properties and hydrological response in different climate zones such as temperate continental. Nowadays, standard methodology for determination soil erodability (rainfall simulator type, rainfall intensities, experimental time, plot sizes) often differs, which enables to make a qualitative comparison of the different land-use impact on soil erosion.58,59 This work aims to study the impact of croplands and vineyards on soil properties and hydrological response in the Podunavlje region (Croatia). Both land uses are traditional in this area. The specific objectives are analyzing the effects on different soil properties (bulk density [BD], soil water content [SWC], mean weight diameter [MWD], water stable aggregates [WSA], soil organic matter [SOM], available phosphorus [AP] and total nitrogen [TN]) and their impact on the hydrological response: runoff (Run), sediment content (SC), sediment loss (SL), carbon loss (C loss), phosphorus loss (P loss), and nitrogen loss (N loss).

Materials and Methods

Study site

The study was carried out in 2 sites in Eastern Croatia (Zmajevac: 45°21′N; 13°26′E, 296 m a.s.l and Erdut: 45°48′N; 18°46′E, 155 m a.s.l), at an average elevation of 296 and 155 m above sea level, respectively (Figure 1). The landscape of the studied area was mainly flat. The parent material is loess, and both soils are loamy textured, classified as anthrosols created from chernozems.60 Basic soil properties are presented in Table 1. The climate is moderate continental with Cfwbx description, according to Köppen climate classification.61 The average annual precipitation (1980-2018) is 677.9 mm, ranging from a minimum of 317.0 mm (2000) to a maximum of 1038.2 mm (2010) (Hydrological and Meteorological Service of Croatia). The rainfall amount is the highest during June and May and the lowest in the period from January to March. However, high inter-annual and annual variations are usual too. The mean annual temperature is 11.5°C, where January is the coldest (0.1°C) and July the warmest (22.2°C) (Hydrological and Meteorological Service of Croatia). Cropland is the dominant land use in this region. However, vineyards and orchards cover the slopes in the study area.

Figure 1.

Study area.


Table 1.

General soil properties in investigated areas.


Cropland and vineyard were used as treatments in Erdut and Zmajevac. Before initiating the experiment, we collected information from landowners about land-use practices. At Erdut site, vineyards are subsoiled every second season to 50 cm depth during spring, followed by cultivation (≈10 cm) or harrowing (≈10 cm), depending on soil moisture conditions. During the season, vegetation covers the inter-row position, whereas between the vines, herbicides are used to control the weeds. Following the season, inter-rows are not tilled, only mulched 2 to 3 times. Each year cropland soil is tilled using a mouldboard plow, followed by disking and roto-harrow and seedbed preparation before sowing. Herbicides and insecticides are used annually. Primary tillage for summer crops is implemented in October or November in the previous autumn, and supplementary tillage followed in the spring before planting. Tillage practices for winter crops (primary and secondary) were carried out in September or October. The crops that are grown followed a typical rotation for this area, which included maize, soybean, winter wheat, sunflower, and barley. During the measurements, winter wheat was a crop, whereas the preceding crop was sunflower. At Zmajevac location, annual vineyard management consists of chiseling (≈20 cm) in autumn, followed by tine tillage (≈12 cm) in the spring. Cropland management at Zmajevac is similar to the cropland at the Erdut site. Winter wheat was the crop, whereas the preceding crop was oil seed rape used as green manure. Vineyard and cropland treatments at both locations were similarly managed for at least 10 years.

Field experiments, soil sampling, and rainfall simulation experiments

In each treatment, a transect was established to carry out the experiments. We selected 8 sampling points, separated by 6 m. At each point, we sampled soils (disturbed and undisturbed—10 cm ring) before carrying out rainfall simulation experiments (8 per treatment, 16 per location, 32 in total). Plots were established in nontraffic areas on each plot. The rainfall simulation experiments were performed during April 2019 under relatively dry soil conditions (15% SWC). Rainfall simulation experiments were carried out in cropland in the stage of stem extension, whereas vineyard was tilled and generally bare (Figure 1). Rainfall simulator (UGT Rainmaker, Munich, Germany) was used in this experiment. Rainfall intensity was adjusted by the time that the nozzle (VeeJet 80/100 nozzle; pressure at 0.5 bar) remains at the reversal points and nozzle turning speed.62 Plastic collectors (n = 144) were placed under the rainfall simulator to collect the drops. After 30 minutes of the experiment at an intensity of 58 mm h-1, we observed an average of 64 g (CV 4.66). The mean drop size was 0.7 mm, and the mean falling velocity was 6.263 m s-1. Plots used in this work are circular plots of 0.785 m2 (metal ring of 100 cm diameter). This intensity was chosen because 93% of the annual soil loss was detected in a single rainstorm event at a rainfall intensity of 59 mm h-1.2 Before the simulations, the simulator was calibrated using the plastic vessel of known dimensions. The slope was measured inside the ring area.

Undisturbed soil samples were used to measure SWC and BD. These properties were determined by drying core samples in an oven at 105°C for 48 hours.63 Additional undisturbed samples were collected to determine MWD, according to Diaz-Zorita et al.64 The percentage of WSA was determined and calculated by the procedure described in Kemper and Rosenau.65 Disturbed samples were dried in the laboratory for 5 days at room temperature (20°C-23°C). Subsequently, samples were sieved with <2 mm mesh to analyze chemical properties. The SOM content was calculated according to the digestion method.66 The AP content was determined by extraction with the ammonium lactate (AL) method67 using a spectrophotometer (model DR/2000; Hach, Dusseldorf, Germany). Total nitrogen in soils and carbon (C) and nitrogen (N) in sediments were obtained by a dry combustion method using Vario MACRO CHNS analyzer.

During rainfall simulation experiments, the overland flow was stored in plastic canisters. Canisters were weighed and filtered to obtain overland flow. Sediment yield was determined after air-drying at room temperature (20°C-23°C) and weighting of the filter paper. Mass of the sediment was deduced from the mass of overland flow to obtain the Run. Sediment content was calculated by dividing the mass of the sediment with the mass of the overland flow. Dried sediments were milled and passed through 2 mm mesh as a preparation for C, N, and P determination.

Statistical data processing

Prior to the statistical analysis, the Shapiro-Wilk (S-W) and Levene tests were applied to test the data normality and homogeneity of the variances (P > .05). Most of the variables did not respect the Gaussian distribution and heteroscedasticity. Therefore, several data normalization methods were performed to achieve the data normality, including natural logarithm, logarithm with base 10, and Box-Cox transformation. From all the tests, Box-Cox transformed data followed data normality and homogeneity of the variances. Therefore, it was used to apply the 2-way analysis of variance. The site and land use were used as factors. In cases where significant differences were found, the Tukey honestly significant difference post hoc test was applied. Significant differences were considered at P < .05. The data presented in the tables are the original one. A principal component analysis (PCA) was performed using the Box-Cox transformed data and was based on the correlation matrix to identify correlations among both soil and overland flow variables. No rotation procedure was applied. Statistical analyses were carried out using Statistica 12.0 for Windows (StatSoft, Tulsa, USA). Graphics were done using Plotly


Environmental plot characteristics

The slope did not differ significantly between treatments in both locations. The slope was significantly higher in Zmajevac vineyard compared with the Erdut one (Table 2). Bulk density of the soils ranged from 1.20 to 1.44 g cm-3. Bulk density was significantly higher in the cropland than in the vineyard at Erdut, whereas at Zmajevac it did not differ between land uses. However, BD was significantly higher in Erdut cropland than in Zmajevac. The SWC and MWD varied between 13.4% and 19.9% and 2.04 and 2.93 mm, respectively. In Zmajevac, SWC was significantly lower in the vineyard than in the cropland, whereas MWD was significantly higher in the vineyard compared with the cropland. Mean weight diameter was significantly lower in the cropland than in the vineyard at both sites. Water stable aggregates and SOM ranged from 58.9% to 88.8% and from 0.7% to 1.7%, respectively. Water stable aggregates were significantly higher in both land uses in Erdut compared with Zmajevac. In Erdut’s location, WSA was significantly higher in the vineyard than in the cropland. The SOM content was significantly lower in the cropland than in the vineyard at Zmajevac, whereas the cropland in Erdut had significantly higher SOM than the cropland in Zmajevac. Finally, AP and TN ranged from 165.5 to 417.2 mg kg-1 and from 0.04% to 0.1%, respectively. The AP content was significantly lower in cropland and vineyard in the Zmajevac site, compared with Erdut. The TN content was significantly higher in Zmajevac cropland compared with the Erdut one. In the Zmajevac site, TN was significantly lower in the cropland compared with the vineyard (Table 2).

Table 2.

Results of 2-way analysis of variance (n = 32).


Hydrological response

The effects of soil management on the overland flow properties are summarized in Table 3. Run and SC varied from 28.0 to 82.0 m3 ha-1 and from 3.6 to 12.3 g kg-1, respectively. In Erdut, vineyard plots had a significantly higher Run than at vineyard in the Zmajevac site. In the Erdut site, vineyard Run was significantly higher compared with cropland. The SC was significantly higher at cropland than in vineyard plots in both sites. Also, SL ranged from 2.8 to 8.2 L m-2. The SL content was significantly lower in the vineyard than in the cropland. Finally, C loss ranged from 0.86 to 2.23 kg m-2. It was significantly lower in Erdut cropland compared with Zmajevac. The vineyard located in Zmajevac had significantly lower C losses than cropland. Finally, N loss and P loss ranged from 0.004 to 0.013 g m-2 and 0.604 and 0.570 g m-2, respectively. In both cases, no differences were observed between sites and land use.

Table 3.

Results of 1-way analysis of variance (n = 32).


PCA analysis

The first 3 factors explained 87.4% of the total variance. Factor 1 explained 32.8% of all variance, whereas factors 2 and 3 explained 27.6% and 18.0%, respectively. Also, factor 1 had high positive loadings in WSA, SOM, AP, and TN and high negative loadings in slope, SWC, and C loss. Factor 2 had high negative loading in WSA and high positive loading in SL, P loss, and N loss. Finally, factor 3 had high positive loading in BD and SC and high negative loading in Run (Table 4). The intersection between factor 1 and factor 2 showed that WSA, SOM, AP, and TN are inversely related to slope, SWC, C loss, SC, and SL (Figure 2A). The impact of rainfall on soil properties and the hydrological response was more different in Zmajevac than in the Erdut site. The highest variability was observed in Zmajevac plots (Figure 2B).

Table 4.

Loading matrix with the first 3 factors extracted from the principal component analysis.


Figure 2.

Relation between factor 1 and factor 2: (A) variables and (B) cases. AP indicates available phosphorus; BD, bulk density; C loss, carbon loss; MWD, mean weight diameter; N loss, nitrogen loss; P loss, phosphorus loss; Run, runoff; SC, sediment content; SL, sediment loss; SOM, soil organic matter; SWC, soil water content; TN, Total nitrogen; WSA, water stable aggregates.



To understand the impact of land use on soils, it is crucial to recognize the impacts of soil structure. This is key to understand the hydrological response. The results obtained indicate that the cropland treatment had implications for BD in both study sites. In this study, BD was higher under cropland treatment than under the vineyard. Previous research generally studies the difference between forests, grasslands (or pastures), and croplands and reveals that BD is the highest in the croplands.37,42,69 Direct comparison of soil properties in vineyards and croplands is scarce. Bogunovic et al28 found that the soil BD was greater in croplands compared with tine-tilled vineyard treatments on a clay loam soil, whereas on sandy loam soils in South Africa Materechera70 found higher BD at cultivated vegetable croplands in addition to untilled vineyards. High topsoil BD is usually attributed to compaction, structural damage, and destruction of macro-pores of topsoil by use of machinery71 or intensive tillage operations.72 Our results are thus consistent because croplands were managed annually with heavy machinery when compared with light tractor–used management in vineyards.

The intensity of land-use practices has an impact on soil structure as it has been observed in other studies.25,72 Our study results reveal that croplands have lower MWD and WSA in addition to vineyards. Such findings show that in the area of research, vineyard management favoring aggregates stability and enhances soil resistance to disaggregation. Small and unstable aggregates were found on soils that are heavily plowed, whereas more persistent soil structure was found on soils with noninvertive tillage operations.73 Although direct studies involving vineyards and cropland are missing, we can support this statement with tillage management comparison. In different environments on loam,74 sandy,75 and clay76 soils, chiseling and/or tine tillage management enhance soil aggregation in addition to plowed soils. The higher SOM concentrations provide higher aggregate stability, as has been observed in previous studies.77 Intensive tillage distorts soil aggregates and exposes them to air, enhancing microbial activity and accelerating decomposition and mineralization of SOM.44 Although our study results for SOM support the previous statement only in Erdut, no significant difference in SOM concentrations between treatments in Zmajevac can be explained by the occurrence of green manure (oil rapeseed) buried before winter wheat sowing. Therefore, the soil compaction observed for cropland soils is likely the result of tillage intensity and field traffic at a particular site, which modifies soil structure.71

Mitigation of soil erosion is obligatory in agriculture land to achieve sustainable management because of the high soil erosion rates that occur under intensive agro-technical management. This is especially important in agricultural soils in the Continental climate of Croatia, where bare topsoil prevails due to intense tillage or use of herbicides. In the case of Croatia, most of the farmers still perform conventional tillage,21,27 indicating that current agricultural soil management is unsustainable. In this context, the adoption of environmental-friendly soil management is required because low soil quality and increased soil erosion act as a consequence of unsustainable soil management.78 Our measurements in 16 cropland plots on 2 locations reveal high soil erosion rates with an average value of 43.6 g m-2 per 1 rainstorm. Vineyard losses were in the average value of 22.2 g m-2. Despite the differences in catchment size plots and dropping systems, other rainfall simulation experiments in vineyards report losses from 2.6 to 54.7 g m-2,12 and from 31 to 186 g m-2,79 at a rainfall intensity from 51 to 70 mm h-1 and at 60 mm h-1, respectively. This study demonstrates that the use of intensive tillage (cropland) increases Run, SC, and soil erosion rates. Although measured in different cropping systems, our study agrees with others. In the Swiss midlands, plowing seems to be responsible for 9.8 times higher soil erosion rates in addition to noninvertive tillage.80 In clay soil, Bertol et al81 reported 34% higher soil erosion rates on plowed plots in addition to harrowed plots. Similar results were reported by Dickey et al13 for clay loam soil when they compared plowing and chisel tillage. Such results reveal that noninvertive, less-intensive tillage practice conserves soil better than plowing. Soil conserving tillage practices have been recognized as effective methods for controlling soil erosion.8,15,17

The PCA revealed that WSA, SOM, AP, and TN are positively associated. This agrees with other findings identified in the literature because usually high SOM concentration in soils increases WSA and acts as a source of nutrients.73,82 This dynamic was observed in the Zmajevac location with a high SOM (Table 1). As mentioned previously, intensive tillage like plowing affects soil aggregate size, exposes them to air and fastens the decomposition of SOM, and reduces the aggregate stability.19,28,75 Despite the increase in SOM, the studied soils record the decrease in C loss. Soil with high SOM content has a stabile structure where organic carbon was entrapped and unexposed to air, making aggregates resistant to particle detachment during rainfall.44 Factor 1 relates slope and SWC. Usually, SWC decreases with an increasing slope.53 In our case, SWC was high, where the slope was high as well. Nevertheless, the slope of these plots was small, and therefore this relationship is very likely an artifact. Factor 2 shows that the relation between MWD and SL, P loss, and N loss is negative. The MWD had the highest values in vineyard plots, whereas SL, P loss, and N loss had the lowest values (Table 1 and 2). The high values of MWD were attributed to favorable soil structure, which helps to reduce the SL,78 as discussed above. Factor 3 relates BD and SC, whereas Run was negatively associated with the first 2 properties. This agrees with other studies that confirm the lower SC on plots with higher Run.83 High SC in Run may be attributed to the effect of tillage practices that generally reduce the size of MWD51 and increase the vulnerability of soil particles to be detached by raindrop impact. This is especially evident in plowed soils,14 as we observed in our work (Tables 1 and 2). However, low BD usually decreases Run,83 although our results showed differently. This was observed in both locations (Table 1 and 2). To our knowledge, no studies were performed on BD critical limit for water infiltration that was carried out in chernozems. However, Wilson et al84 observed that this threshold was 1.44 g cm-3 for mollisol. We hypothesize that the limits for chernozems can be higher than 1.44 g cm-3 in mollisol. More research is necessary to reveal unanswered questions in this process.

Overall, the fact that soil has a low BD does not necessarily mean that it decreases Run. The type of tillage management trough disturbance is a crucial factor in controlling soil erosion. Noninvertive tilled soils (vineyards) have better structure, lover compaction, and are more resistant to sediment detachment and soil loss.


Different land-use practices affect soil properties and hydrological response. Noninvertive tillage practices reduced soil BD and increased WSA and MWD in vineyards. In these plots, the better preservation of soil structure enhanced soil quality. Also, intensive invert tillage like plowing increased the SC and SL. Vineyard soils with loosening-type tillage had a high infiltration and resistance to soil detachment, and this shows that vineyard management in the studied area on loamy soils represents the less endangered practice and land use with higher resistance to rainfall impact. From a hydrological response perspective, the practices applied in croplands need to be either reconsidered or coupled with conservation strategies to mitigate soil erosion problems. Outcomes support the fact that management practices on the studied area should use reduced tilling frequency to have persistent structure and to avoid soil erosion. This study contributes to better land-use management in Eastern Croatia on Continental-type climate.


We are grateful to firms Erdutski vinogradi d.o.o. and Vinarija Josić for their help during field work. Aleksandra Percin is acknowledged for her cooperation during the laboratory work.

Author Contributions

Conceptualization and methodology: IB; Formal analysis: LJT and PP; Investigation: IB; Writing—original draft preparation: IB; Writing—review and editing: PP and LJT; Visualizatio: PP; Supervision: IB and PP.



Vanwalleghem T , Gómez JA , Infante Amate J , et al. Impact of historical land use and soil management change on soil erosion and agricultural sustainability during the Anthropocene. Anthropocene. 2017;17:13–29. Google Scholar


Biddoccu M , Ferraris S , Opsi F , Cavallo E. Long-term monitoring of soil management effects on runoff and soil erosion in sloping vineyards in Alto Monferrato (North–West Italy). Soil Till Res. 2016;155:176–189. Google Scholar


Nearing MA , Xie Y , Liu B , Ye Y. Natural and anthropogenic rates of soil erosion. Int Soil Water Conserv Res. 2017;5:77–84. Google Scholar


Skinner RJ , Chambers BJ. A survey to assess the extent of soil water erosion in lowland England and Wales. Soil Use Manage. 1996;12:214–220. Google Scholar


Mhazo M , Chivenge P , Chaplot V. Tillage impact on soil erosion by water: discrepancies due to climate and soil characteristics. Agri Ecosyst Environ. 2016;230:231–241. Google Scholar


Pimentel D , Burgess M. Soil erosion threatens food production. Agriculture. 2013;3:443–463. Google Scholar


Cerdan O , Govers G , Le Bissonnais Y , et al. Rates and spatial variations of soil erosion in Europe: a study based on erosion plot data. Geomorphology. 2010;122:167–177. Google Scholar


Prosdocimi M , Cerdà A , Tarolli P. Soil water erosion on Mediterranean vineyards: a review. Catena. 2016;141:1–21. Google Scholar


Rodrigo-Comino J. Five decades of soil erosion research in “terroir.” The State-of-the-Art. Earth-Sci Rev. 2018;179:436–447. Google Scholar


Verheijen FG , Jones RJ , Rickson RJ , Smith CJ. Tolerable versus actual soil erosion rates in Europe. Earth-Sci Rev. 2009;94:23–38. Google Scholar


Blavet D , De Noni G , Le Bissonnais Y , et al. Effect of land use and management on the early stages of soil water erosion in French Mediterranean vineyards. Soil Till Res. 2009;106:124–136. Google Scholar


Arnaez J , Lasanta T , Ruiz-Flaño P , Ortigosa L. Factors affecting runoff and erosion under simulated rainfall in Mediterranean vineyards. Soil Till Res. 2007;93:324–334. Google Scholar


Dickey EC , Shelton DP , Jasa PJ , Peterson TR. Soil erosion from tillage systems used in soybean and corn residues. Trans ASAE. 1985;28:1124–1130. Google Scholar


Ramos MC , Pareja-Sánchez E , Plaza-Bonilla D , Cantero-Martínez C , Lampurlanés J. Soil sealing and soil water content under no-tillage and conventional tillage in irrigated corn: effects on grain yield. Hydrol Process. 2019;33:2095–2109. Google Scholar


Cowie AL , Orr BJ , Castillo Sanchez VM , et al. Land in balance: the scientific conceptual framework for Land Degradation Neutrality. Environ Sci Policy. 2018;79:25–35. Google Scholar


Basic F , Kisic I , Mesic M , Nestroy O , Butorac A. Tillage and crop management effects on soil erosion in central Croatia. Soil Till Res. 2004;78:197–206. Google Scholar


Lal R. Mechanized tillage systems effects on soil erosion from an Alfisol in watersheds cropped to maize. Soil Tillage Res. 1984;4:349–360. Google Scholar


Lal R. , Soil Erosion in the Tropics: Principles and Management. New York, NY: McGraw Hill; 1990. Google Scholar


Kisic I , Bogunovic I , Birkás M , Jurisic A , Spalevic V. The role of tillage and crops on a soil loss of an arable Stagnic Luvisol. Arch Agron Soil Sci. 2017;63:403–413. Google Scholar


Ziadat FM , Taimeh AY. Effect of rainfall intensity, slope, land use and antecedent soil moisture on soil erosion in an arid environment. Land Degrad Dev. 2013;24:582–590. Google Scholar


Williams JD , Wilkins DE , Douglas JrCL , Rickman RW. Mow-plow crop residue management influence on soil erosion in north-central Oregon. Soil Till Res. 2000;55:71–78. Google Scholar


Cantón Y , Solé-Benet A , Asensio C , Chamizo S , Puigdefábregas J. Aggregate stability in range sandy loam soils relationships with runoff and erosion. Catena. 2009;77:192–199. Google Scholar


Johnson CB , Moldenhauer WC. Effect of chisel versus moldboard plowing on soil erosion by water. Soil Sci Soc Am J. 1979;43:177–179. Google Scholar


Mermut AR , Acton DF , Eilers WD. Estimation of soil erosion and deposition by a landscape analysis technique on clay soils in southwestern Saskatchewan. Can J Soil Sci. 1983;63:727–739. Google Scholar


Coulouma G , Boizard H , Trotoux G , Lagacherie P , Richard G. Effect of deep tillage for vineyard establishment on soil structure: a case study in Southern France. Soil Till Res. 2006;88:132–143. Google Scholar


Capello G , Biddoccu M , Ferraris S , Cavallo E. Effects of tractor passes on hydrological and soil erosion processes in tilled and grassed vineyards. Water. 2019;11:2118. Google Scholar


Bogunovic I , Pereira P , Kisic I , Sajko K , Sraka M. Tillage management impacts on soil compaction, erosion and crop yield in Stagnosols (Croatia). Catena. 2018;160:376–384. Google Scholar


Bogunovic I , Telak LJ , Pereira P. Agriculture management impacts on soil properties and hydrological response in Istria (Croatia). Agronomy. 2020;10:282. Google Scholar


Dunjó G , Pardini G , Gispert M. The role of land use–land cover on runoff generation and sediment yield at a microplot scale, in a small Mediterranean catchment. J Arid Environ. 2004;57:239–256. Google Scholar


Arhonditsis G , Giourga C , Loumou A. Ecological patterns and comparative nutrient dynamics of natural and agricultural Mediterranean-type ecosystems. Environ Manage. 2000;26:527–537. Google Scholar


Bini C , Gemignani S , Zilocchi L. Effect of different land use on soil erosion in the pre-alpine fringe (North-East Italy): ion budget and sediment yield. Sci Total Environ. 2006;369:433–440. Google Scholar


García-Ruiz JM , Lasanta T , Ortigosa L , Ruiz-Flaño P , Martí C , González C. Sediment yield under different land uses in the Spanish Pyrenees. Mt Res Dev. 1995;15:229–240. Google Scholar


Romero-Díaz A , Cammeraat LH , Vacca A , Kosmas C. Soil erosion at three experimental sites in the Mediterranean. Earth Surf Proc Land. 1999;24:1243–1256. Google Scholar


Kosmas C , Danalatos N , Cammeraat LH , et al. The effect of land use on runoff and soil erosion rates under Mediterranean conditions. Catena. 1997; 29:45–59. Google Scholar


Nunes AN , De Almeida AC , Coelho CO. Impacts of land use and cover type on runoff and soil erosion in a marginal area of Portugal. Appl Geogr. 2011;31:687–699. Google Scholar


Zuazo VD , Martínez JF , Raya AM. Impact of vegetative cover on runoff and soil erosion at hillslope scale in Lanjaron, Spain. Environmentalist. 2004;24:39–48. Google Scholar


Bormann H , Klaassen K. Seasonal and land use dependent variability of soil hydraulic and soil hydrological properties of two Northern German soils. Geoderma. 2008;145:295–302. Google Scholar


Maetens W , Poesen J , Vanmaercke M. How effective are soil conservation techniques in reducing plot runoff and soil loss in Europe and the Mediterranean? Earth Sci Rev. 2012;115:21–36. Google Scholar


Bewket W , Stroosnijder L. Effects of agroecological land use succession on soil properties in Chemoga watershed, Blue Nile basin, Ethiopia. Geoderma. 2003;111:85–98. Google Scholar


Taye G , Poesen J , Wesemael BV , et al. Effects of land use, slope gradient, and soil and water conservation structures on runoff and soil loss in semi-arid Northern Ethiopia. Phys Geogr. 2013;34:236–259. Google Scholar


El-Hassanin AS , Labib TM , Gaber EI. Effect of vegetation cover and land slope on runoff and soil losses from the watersheds of Burundi. Agric Ecosyst Environ. 1993;43:301–308. Google Scholar


Celik I. Land-use effects on organic matter and physical properties of soil in a southern Mediterranean highland of Turkey. Soil Till Res. 2005;83:270–277. Google Scholar


Gao X , Wu P , Zhao X , Wang J , Shi Y. Effects of land use on soil moisture variations in a semi-arid catchment: implications for land and agricultural water management. Land Degrad Dev. 2014;25:163–172. Google Scholar


Liu M , Han G , Zhang Q. Effects of soil aggregate stability on soil organic carbon and nitrogen under land use change in an erodible region in Southwest China. Int J Environ Res Publ Health. 2019;16:3809. Google Scholar


Peng T , Wang SJ. Effects of land use, land cover and rainfall regimes on the surface runoff and soil loss on karst slopes in southwest China. Catena. 2012;90:53–62. Google Scholar


Wei W , Chen L , Fu B , Huang Z , Wu D , Gui L. The effect of land uses and rainfall regimes on runoff and soil erosion in the semi-arid loess hilly area, China. J Hydrol. 2007;335:247–258. Google Scholar


El Kateb H , Zhang H , Zhang P , Mosandl R . Soil erosion and surface runoff on different vegetation covers and slope gradients: a field experiment in Southern Shaanxi Province, China. Catena. 2013;105:1–10. Google Scholar


Feng Q , Zhao W , Wang J , et al. Effects of different land-use types on soil erosion under natural rainfall in the Loess Plateau, China. Pedosphere. 2016;26:243–256. Google Scholar


Kothyari BP , Verma PK , Joshi BK , Kothyari UC. Rainfall–runoff-soil and nutrient loss relationships for plot size areas of Bhetagad watershed in Central Himalaya, India. J Hydrol. 2004;293:137–150. Google Scholar


Mohammad AG , Adam MA. The impact of vegetative cover type on runoff and soil erosion under different land uses. Catena. 2010;81:97–103. Google Scholar


Nath AJ , Lal R. Effects of tillage practices and land use management on soil aggregates and soil organic carbon in the North Appalachian Region, USA. Pedosphere. 2017;27:72–176. Google Scholar


Jacinthe PA , Lal R , Owens LB , Hothem DL. Transport of labile carbon in runoff as affected by land use and rainfall characteristics. Soil Till Res. 2004;77:111–123. Google Scholar


Gilley JE , Doran JW , Karlen DL , Kaspar TC. Runoff, erosion, and soil quality characteristics of a former Conservation Reserve Program site. J. Soil Water Conserv. 1997;52:189–193. Google Scholar


Rimal BK , Lal R. Soil and carbon losses from five different land management areas under simulated rainfall. Soil Till Res. 2009;106:62–70. Google Scholar


Návar J , Synnott TJ. Surface runoff, soil erosion, and land use in Northeastern Mexico. Terra Latinoamericana. 2000;18:247–253. Google Scholar


Dos Santos JCN , de Andrade EM , Medeiros PHA , Guerreiro MJS , de Queiroz Palácio HA. Land use impact on soil erosion at different scales in the Brazilian semi-arid. Rev Cienc Agron. 2017;48:251–260. Google Scholar


Molina A , Govers G , Vanacker V , Poesen J , Zeelmaekers E , Cisneros F. Runoff generation in a degraded Andean ecosystem: interaction of vegetation cover and land use. Catena. 2007;71:357–370. Google Scholar


Mayerhofer C , Meißl G , Klebinder K , Kohl B , Markart G. Comparison of the results of a small—plot and a large—plot rainfall simulator—Effects of land use and land cover on surface runoff in Alpine catchments. Catena. 2017;156:184–196. Google Scholar


Dunkerley D. How is overland flow produced under intermittent rain? An analysis using plot-scale rainfall simulation on dryland soils. J Hydrol. 2018;556:119–130. Google Scholar


IUSS Working Group WRB. World Reference Base for Soil Resources 2014, Update 2015: International Soil Classification System for Naming Soils and Creating Legends for Soil Maps (World Soil Resources Reports No. 106, 192); 2015, Google Scholar


Kottek M , Grieser J , Beck C , Rudolf B , Rubel F. World map of the Köppen-Geiger climate classification updated. Meteorol Z. 2006;15:259–263. Google Scholar


Schindewolf M , Schmidt J. Parameterization of the EROSION 2D/3D soil erosion model using a small-scale rainfall simulator and upstream runoff simulation. Catena. 2012;91:47–55. Google Scholar


Black CA . Method of Soil Analysis. Part 2. Chemical and Microbiological Properties. 1st ed. Madison, WI: American Society of Agronomy; 1965:771–1569. Google Scholar


Diaz-Zorita M , Perfect E , Grove JH. Disruptive methods for assessing soil structure. Soil Till Res. 2002;64:3–22. Google Scholar


Kemper WD , Rosenau RC . Aggregate stability and size distribution. In: Klute A , , ed. Methods of Soil Analysis, Part 1. Physical and Mineralogical Methods. Vol. 1, 2nd ed.Madison, WI: American Society of Agronomy; 1986:425–442. Google Scholar


Walkly A , Black IA. An examination of digestion methods for determining soil organic matter and a proposed modification of the chromic and titration. Soil Sci Soc Am J. 1934;37:29–38. Google Scholar


Egnér HANS , Riehm H , Domingo WR . Untersuchungen über die chemische Bodenanalyse als Grundlage für die Beurteilung des Nährstoffzustandes der Böden. II. Chemische Extraktionsmethoden zur Phosphor-und Kaliumbestimmung. Lantbrukshögskolans Ann. 1960;26:199–215. Google Scholar


Plotly. (accessed on 17 March 2020) Google Scholar


Murty D , Kirschbaum MU , Mcmurtrie RE , Mcgilvray H. Does conversion of forest to agricultural land change soil carbon and nitrogen? A review of the literature. Global Change Biol. 2002;8:105–123. Google Scholar


Materechera SA. Tillage and tractor traffic effects on soil compaction in horticultural fields used for peri-urban agriculture in a semi-arid environment of the North West Province, South Africa. Soil Till Res. 2009;103:11–15. Google Scholar


Hamza MA , Anderson WK. Soil compaction in cropping systems: a review of the nature, causes and possible solutions. Soil Till Res. 2005;82:121–145. Google Scholar


Dekemati I , Simon B , Vinogradov S , Birkás M. The effects of various tillage treatments on soil physical properties, earthworm abundance and crop yield in Hungary. Soil Till Res. 2019;194:104334. Google Scholar


Birkás M , Szemők A , Antos G , Neményi M. , Environmentally-Sound Adaptable Tillage. Amsterdam, The Netherlands: Akadémiai Kiadó; 2008. Google Scholar


Pagliai M , Vignozzi N , Pellegrini S. Soil structure and the effect of management practices. Soil Till Res. 2004;79:131–143. Google Scholar


Carter MR. Influence of reduced tillage systems on organic matter, microbial biomass, macro-aggregate distribution and structural stability of the surface soil in a humid climate. Soil Till Res. 1992;23:361–372. Google Scholar


Pagliai M , Raglione M , Panini T , Maletta M , La Marca M. The structure of two alluvial soils in Italy after 10 years of conventional and minimum tillage. Soil Till Res. 1995;34:209–223. Google Scholar


Mujdeci M , Isildar AA , Uygur V , Alaboz P , Unlu H , Senol H. Cooperative effects of field traffic and organic matter treatments on some compaction-related soil properties. Solid Earth. 2017;8:189–198. Google Scholar


Pereira P , Bogunovic I , Muñoz-Rojas M , Brevik EC. Soil ecosystem services, sustainability, valuation and management. Curr Opin Environ Sci Health. 2018;5:7–13. Google Scholar


Battany MC , Grismer ME. Rainfall runoff and erosion in Napa Valley vineyards: effects of slope, cover and surface roughness. Hydrol Process. 2000;14:1289–1304. Google Scholar


Prasuhn V. On-farm effects of tillage and crops on soil erosion measured over 10 years in Switzerland. Soil Till Res. 2012;120:137–146. Google Scholar


Bertol I , Mello EL , Guadagnin JC , Zaparolli ALV , Carrafa MR. Nutrient losses by water erosion. Sci Agr. 2003;60:581–586. Google Scholar


Bot A , Benites J. , The Importance of Soil Organic Matter: Key to Drought-Resistant Soil and Sustained Food Production (No. 80). Rome, Italy: Food and Agriculture Organization; 2005. Google Scholar


Keesstra S , Pereira P , Novara A , et al. Effects of soil management techniques on soil water erosion in apricot orchards. Sci Total Environ. 2016;551:357–366. Google Scholar


Wilson MG , Sasal MC , Caviglia OP. Critical bulk density for a Mollisol and a Vertisol using least limiting water range: effect on early wheat growth. Geoderma. 2013;192:354–361. Google Scholar
© The Author(s) 2020 This article is distributed under the terms of the Creative Commons Attribution-NonCommercial 4.0 License ( which permits non-commercial use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (
Igor Bogunovic, Leon Josip Telak, and Paulo Pereira "Experimental Comparison of Runoff Generation and Initial Soil Erosion Between Vineyards and Croplands of Eastern Croatia: A Case Study," Air, Soil and Water Research 13(1), (12 June 2020).
Received: 27 April 2020; Accepted: 28 April 2020; Published: 12 June 2020
land use
rainfall simulation experiments
soil erosion
Back to Top