Infiltration and Soil Loss Changes during the Growing Season under Ploughing and Conservation Tillage

Decreased water retention and increased runoff and soil loss are of special importance concerning soil degradation of hilly crop fields. In this study, plots under ploughing (conventional) tillage (PT) and conservation tillage (CT; 15 years) were compared. Rainfall simulation on 6 m2 plots was applied to determine infiltration and soil loss during the growing season. Results were compared with those measured from 1200 m2 plots exposed to natural rainfalls in 2016. Infiltration was always higher under CT than PT, whereas the highest infiltration was measured under the cover crop condition. Infiltration under seedbed and stubble resulted in uncertainties, which suggests that natural pore formation can be more effective at improving soil drainage potential than can temporary improvements created by soil tillage operations. Soil erodibility was higher under PT for each soil status; however, the seedbed condition triggered the highest values. For CT, soil loss volume was only a function of runoff volume at both scales. Contrarily, on PT plots, some extreme precipitation events triggered extremely high soil loss owing to linear erosion, which meant no direct connection existed between the scales. Improved soil conditions due to conservation practice are more important for decreasing soil loss than the better surface conditions.


Introduction
Soil water retention capacity determines most of the physical, chemical, and biological properties of soil and is responsible for soil fertility.In many parts of the world, mainly in semiarid areas, available water content limits crop production.In such places, it is essential to store as much precipitation in the soil as is possible, for crop production that is more secure [1,2].Infiltration, one of the most important soil properties, is responsible for water entering the soil; however, it is not easy to measure and predict, as it is highly variable both spatially and temporally [3].Some major driving forces, such as soil surface crusting due to the high kinetic energy of falling raindrops [4,5]; soil sealing on the surface or at deeper levels (e.g., plough pan [6]); and rill and gully formation, which decreases water contact to the soil surface, decrease infiltration [7], whereas other forces, such as increasing plant growth, surface roughness, and surface coverage [8], increase the volume of infiltration.Ploughing tillage (PT; conventional cultivation system based on annual moldboard ploughing tillage) plays an ambivalent role.Some researchers emphasize that it increases the infiltration volume right after ploughing [9], while others highlight how ploughing tillage destroys the soil's structure [10] and compact the soil [11].Therefore, new or revised practices in order to maximize water moisture during the growing season have been developed and tested worldwide.Most focus on tillage mitigation and surface cover increase using dead (mulch) or living (cover crops) biomass and are called "conservation agriculture" [12,13].In the last two decades, the application of cover crops came into the spotlight [14,15].Although there are many advantages of cover crops, such as physical protection of the surface, organic carbon source as green manure, and soil structure improvement, there are some disadvantages as well (potential source of plant diseases and concurrent soil moisture and nitrogen loss) [16].
Soil has a higher infiltration capacity in dry conditions.As moisture content increases, the infiltration value decreases, until it reaches its theoretical minimum, which is dependent on the soil texture and structure.In reality, this minimum value is not a constant but rather a function of rainfall intensity [17].Higher rainfall intensity triggers a higher final infiltration rate, as described by Rose [18].
where i is the infiltration rate (mm h −1 ), p is the irrigation rate (mm h −1 ), and i s is the soil permeability (mm h −1 ).Moreover, precipitation that is unable to infiltrate creates runoff and triggers soil erosion and flash flooding.The on-site costs due to soil erosion is estimated at $165-410 year −1 ha −1 in the EU states [19].Higher runoff volume does not necessarily mean higher soil loss.Theoretically, a compacted, smooth, crusted soil surface, which is the reason for higher runoff volume, could protect against soil detachment via runoff [20].Accordingly, defense against soil erosion should be second only to water retention as the aim of cultivation.
As most physical properties of the soil are highly variable during the year, it is crucial to take these changes into account in order to create and apply the best management practice [21].To do so, one has to gather measured infiltration and soil erodibility values from in situ plot-scale crop fields.As runoff plots measuring the effects of natural precipitation are not suitable for investigation over a wide spectrum of intensity under certain soil conditions, rainfall simulation was applied [22].This method is highly effective, reproducible, and suitable for the purposes of this study [23][24][25].
The aims of this study were to (i) quantify instantaneous infiltration values as a function of precipitation intensity; (ii) compare these functions under various soil conditions, such as under cover crop, seedbed, stubble, and tillage systems (conventional and conservation); and (iii) compare soil loss results.

Study Site
Measurements were carried out at the Szentgyörgyvár soil erosion experimental site in Hungary (46 • 44 54.1"N; 17 • 08 48.4"E) (Figure 1).PT and conservation tillage (CT) (with no inversion, using a reduced number of tillage operations and leaving minimum 30% crop residues on the soil surface) were compared on four 24 × 50 m plots, two per tillage type.Since 2004, runoff and soil loss of the plots have been collected continuously until now [26].The recent tillage techniques used on the plots are listed in Table 1.The plots' soil is Luvisol on sandy loess, with a particle size distribution of 4% clay (<0.002 mm), 60% silt (0.002-0.02), and 36% sand [27].The mean annual temperature is 11 • C, and precipitation is 628 mm.The slope steepness is 9-10% [28,29].

Rainfall Simulation and Cultivation
A Shower Power 02 in situ rainfall simulator, designed by the Hungarian Geographical Institute, was used for infiltration measurement [30].The size of the irrigated plot was 3 m × 4 m = 12 m 2 .To exclude border effects, only the inner 2 m × 3 m area was investigated.For drop formation, 80,100 VeeJet nozzles were used [31].Two nozzles were placed 2 m apart from each other.The nozzles were exactly in the same plane with each other and the alternating axis.This overlap of the two nozzles ensured homogeneous drop distribution 3 m under the nozzles.The nozzles were elevated above 3 m height, so the biggest drops could reach their final velocity before hitting the ground.The drop spectrum at 0.41 bars of pressure is corresponded well with that of natural rainstorms [32].Intensity varied 30-130 mm h −1 based on the axis alternation frequency, therefore drop spectrum was the same for each precipitation.Intensities were calibrated for certain frequencies using laboratory measurements on rain volume.One alternation theoretically equaled 0.07 mm precipitation.
Rainfall simulation was carried out three times in 2016.The first occurred when a green cover crop was present in April; the second occurred a month later (May), after the green cover crop had

Rainfall Simulation and Cultivation
A Shower Power 02 in situ rainfall simulator, designed by the Hungarian Geographical Institute, was used for infiltration measurement [30].The size of the irrigated plot was 3 m × 4 m = 12 m 2 .To exclude border effects, only the inner 2 m × 3 m area was investigated.For drop formation, 80,100 VeeJet nozzles were used [31].Two nozzles were placed 2 m apart from each other.The nozzles were exactly in the same plane with each other and the alternating axis.This overlap of the two nozzles ensured homogeneous drop distribution 3 m under the nozzles.The nozzles were elevated above 3 m height, so the biggest drops could reach their final velocity before hitting the ground.The drop spectrum at 0.41 bars of pressure is corresponded well with that of natural rainstorms [32].Intensity varied 30-130 mm h −1 based on the axis alternation frequency, therefore drop spectrum was the same for each precipitation.Intensities were calibrated for certain frequencies using laboratory measurements on rain volume.One alternation theoretically equaled 0.07 mm precipitation.
Rainfall simulation was carried out three times in 2016.The first occurred when a green cover crop was present in April; the second occurred a month later (May), after the green cover crop had been destroyed and the soil was disturbed to create a seedbed for planting maize using different soil cultivation methods; and, finally, the third occurred in October just after harvesting.
Runoff and erosion data were generated for CT and PT.Each investigation included five runs at various constant intensities (30,40,60,90, or 120 mm h −1 ), with 9% slope steepness (which is widespread for crop fields in Hungary).Natural precipitation events with 30-40 mm h −1 intensities are casual, higher values are rare in the studied area.The investigated plot was fenced by metal sheets dug into the soil in order to inhibit surface runoff and run-on flow.Right after the set-up process, a 40 mm h −1 rainfall was simulated, yet runoff measurement was not performed.The aim of this first artificial precipitation event was to fill the soil to field capacity in order to ensure standard circumstances for the measurements.
After this, the first pretreatment, measurements were carried out in an increasing order of precipitation intensity.For each measurement, the total amount of runoff was collected and measured.Partial runoff and soil losses were measured in separate units, with special emphasis on time.

Calculation
During each measurement changes in infiltration within the precipitation event were recorded.Based on runoff dynamics, the apparent infiltration intensity was calculated (Figure 2).Sustainability 2017, 9, 1726 4 of 13 been destroyed and the soil was disturbed to create a seedbed for planting maize using different soil cultivation methods; and, finally, the third occurred in October just after harvesting.Runoff and erosion data were generated for CT and PT.Each investigation included five runs at various constant intensities (30,40,60,90, or 120 mm h −1 ), with 9% slope steepness (which is widespread for crop fields in Hungary).Natural precipitation events with 30-40 mm h −1 intensities are casual, higher values are rare in the studied area.The investigated plot was fenced by metal sheets dug into the soil in order to inhibit surface runoff and run-on flow.Right after the set-up process, a 40 mm h −1 rainfall was simulated, yet runoff measurement was not performed.The aim of this first artificial precipitation event was to fill the soil to field capacity in order to ensure standard circumstances for the measurements.
After this, the first pretreatment, measurements were carried out in an increasing order of precipitation intensity.For each measurement, the total amount of runoff was collected and measured.Partial runoff and soil losses were measured in separate units, with special emphasis on time.

Calculation
During each measurement changes in infiltration within the precipitation event were recorded.Based on runoff dynamics, the apparent infiltration intensity was calculated (Figure 2).For each measurement, linear functions were fitted to the last six-ten measured runoff volume data points against time.In this last period, runoff intensity became constant; accordingly, this last section of the data points showed a linear linkage.Estimation of the apparent infiltration intensity was based on the steepness of this linear function (mm s −1 ).These values are listed as apparent infiltration in Table 1.
The infiltration rate for each event was fitted to Equation (1).Hence, the value of is needs to be determined by fitting Equation (1) to the irrigation and infiltration data.This was done manually by adjusting the value of is until the normalized root mean squared error (NRMSE) was minimized, given by the equation: For each measurement, linear functions were fitted to the last six-ten measured runoff volume data points against time.In this last period, runoff intensity became constant; accordingly, this last section of the data points showed a linear linkage.Estimation of the apparent infiltration intensity was based on the steepness of this linear function (mm s −1 ).These values are listed as apparent infiltration in Table 1.
The infiltration rate for each event was fitted to Equation (1).Hence, the value of i s needs to be determined by fitting Equation (1) to the irrigation and infiltration data.This was done manually by adjusting the value of i s until the normalized root mean squared error (NRMSE) was minimized, given by the equation: where P j and O j are the jth predicted and observed values, respectively, of N values, and O av is the mean of the observed values.From this fit, because evapotranspiration is negligible on the time scale of these experiments, it follows that the runoff rate (q) is simply: The amount of eroded sediment, C (kg m −1 ), for each event was fitted to the equation, as given by the Pan-European Soil Erosion Risk Assessment (PESERA) model [33]: where k is the soil erodibility, V is the volume of runoff per unit width (L m −1 ), and S is the local slope gradient (dimensionless).Moreover, soil erodibility K factor values were also calculated according to the USLE model [34] for the entire year of the field scale, because this model is still widely used [35][36][37].
The USLE K factor was created for annual temporal resolution; therefore, for single simulated rainfalls, the modified USLE (MUSLE) [38] was applied.For comparisons, boxplots and one-way ANOVA were applied using IBM SPSS software.

Earthworm Investigation
In order to estimate natural aggregation and porosity formation effects earthworms were surveyed at the end of October by taking samples 10 cm in diameter and 10 cm in depth at nine points per plot with a soil auger.The number and total weight of earthworms per soil core were recorded by applying the method of Harper Adams University College [39].

Results and Discussion
All 15 artificial rainfalls triggered runoff and soil loss from the investigated plots.The main measured and calculated rainfall and soil properties are presented in Table 2.
Table 2. Main parameters of the artificial rainfalls and corresponding infiltration and soil loss values created on the 6 m 2 plots, with 9% slope steepness.The shaded rows were under conventional tillage, whereas the unshaded rows were under conservation tillage.MUSLE = Modified soil loss equation [38]; K = soil erodibility factor.

Hydrological Results for the Small Plots
Based on the calculated final apparent infiltration values, an infiltration function was determined for both tillage systems and each soil status of the growing season (Figure 3).

Hydrological Results for the Small Plots
Based on the calculated final apparent infiltration values, an infiltration function was determined for both tillage systems and each soil status of the growing season (Figure 3).In April, under crop cover, both soils had the highest annual infiltration capacity; however, infiltration was higher under CT.Both the highest infiltration and the highest increase in infiltration owing to precipitation intensity increase were found on the CT plot under the cover crop condition.
Under the seedbed condition in May, both functions indicated a slightly lower infiltration intensity, although one of the main purposes of cultivation is to increase infiltration.This is in accordance with the results of many studies [40][41][42], though infiltration and saturated and unsaturated hydraulic conductivity values measured on undisturbed soil samples showed high variability because of spatial heterogeneity of the soil [21].On the other hand, rainfall simulation was carried out after 2 weeks of tillage, Freese et al. [43] reported results of significant infiltration rate decrease for this period in regard to three tillage systems.
Infiltration was still higher under CT; however, irrigation intensity dependence on apparent infiltration intensity decreased (curve steepness on Figure 3).This means the infiltration rate reduction compared to the cover crop status was higher as the precipitation intensity increased.Accordingly, the advantage of cover crops in increasing infiltration should be emphasized for the heaviest rainstorms.According to climate projections, considerable change in the distribution of precipitation is expected in the Carpathian Basin (wider surrounding of this study).This means lessfrequent but higher-intensity precipitation events, with increased precipitation during winter [44,45].Therefore, it is necessary to sow the cover crop as soon after harvest as is possible because this period is frequently endangered by extreme precipitation events.
In October, under stubble, infiltration was still higher under CT, although it decreased in both tillage systems, compared to the initial values.The function fitted on the PT data was less steep, which means a very limited increase in infiltration response owing to precipitation intensity increase.Over 40 mm h −1 the effect of precipitation intensity was practically negligible.Moreover, the stubble infiltration function under CT was very similar to that of PT under the seedbed condition (Figure 3, Table 2).Accordingly, at the end of the growing season, settled and crusted surface soil can obtain and retain as much water as the recently tilled seedbed under PT.Therefore, in the present case, infiltration In April, under crop cover, both soils had the highest annual infiltration capacity; however, infiltration was higher under CT.Both the highest infiltration and the highest increase in infiltration owing to precipitation intensity increase were found on the CT plot under the cover crop condition.
Under the seedbed condition in May, both functions indicated a slightly lower infiltration intensity, although one of the main purposes of cultivation is to increase infiltration.This is in accordance with the results of many studies [40][41][42], though infiltration and saturated and unsaturated hydraulic conductivity values measured on undisturbed soil samples showed high variability because of spatial heterogeneity of the soil [21].On the other hand, rainfall simulation was carried out after 2 weeks of tillage, Freese et al. [43] reported results of significant infiltration rate decrease for this period in regard to three tillage systems.
Infiltration was still higher under CT; however, irrigation intensity dependence on apparent infiltration intensity decreased (curve steepness on Figure 3).This means the infiltration rate reduction compared to the cover crop status was higher as the precipitation intensity increased.Accordingly, the advantage of cover crops in increasing infiltration should be emphasized for the heaviest rainstorms.According to climate projections, considerable change in the distribution of precipitation is expected in the Carpathian Basin (wider surrounding of this study).This means less-frequent but higher-intensity precipitation events, with increased precipitation during winter [44,45].Therefore, it is necessary to sow the cover crop as soon after harvest as is possible because this period is frequently endangered by extreme precipitation events.
In October, under stubble, infiltration was still higher under CT, although it decreased in both tillage systems, compared to the initial values.The function fitted on the PT data was less steep, which means a very limited increase in infiltration response owing to precipitation intensity increase.Over 40 mm h −1 the effect of precipitation intensity was practically negligible.Moreover, the stubble infiltration function under CT was very similar to that of PT under the seedbed condition (Figure 3, Table 2).Accordingly, at the end of the growing season, settled and crusted surface soil can obtain and retain as much water as the recently tilled seedbed under PT.Therefore, in the present case, infiltration was not limited by the surface conditions but by the conditions in deeper layers.Accordingly, soil porosity formed by natural processes can be more effective in drainage than that formed by tillage induced, temporary, vulnerable processes.The same conclusion was made based on the earthworm survey.A significant difference in both the number and weight of earthworms was revealed following the October rainfall simulation, in favor of CT.This resulted in 2.5 times more and 5.3 times larger earthworms, on average (p < 0.01) (Figure 4).The role of soil fauna is essential in establishing the optimal structure and porosity of the soil.CT does not damage these layer-specific organisms, and a macropore system is established, which is capable of absorbing rainwater quickly [46][47][48].
Sustainability 2017, 9, 1726 7 of 13 was not limited by the surface conditions but by the conditions in deeper layers.Accordingly, soil porosity formed by natural processes can be more effective in drainage than that formed by tillage induced, temporary, vulnerable processes.The same conclusion was made based on the earthworm survey.A significant difference in both the number and weight of earthworms was revealed following the October rainfall simulation, in favor of CT.This resulted in 2.5 times more and 5.3 times larger earthworms, on average (p < 0.01) (Figure 4).The role of soil fauna is essential in establishing the optimal structure and porosity of the soil.CT does not damage these layer-specific organisms, and a macropore system is established, which is capable of absorbing rainwater quickly [46][47][48].In the lower precipitation intensity range, all three soil conditions reacted more or less the same to precipitation intensity increase.This is the case up to 10 mm h −1 regarding CT, whereas under PT, this value is 18 mm h −1 .This suggests that independently from annual variation of single cultivation runs or temporary soil hydraulic conditions, PT has a general long-term effect on infiltration increase, and the earthworms and other soil fauna have a key role in this effect.The calculated soil permeability values also demonstrate these trends (Table 3).CT had higher infiltration values throughout the year, even though most studies have emphasized that infiltration to CT soils increases over the long run, compared to PT [31,49].Interestingly, right after seedbed preparation, CT retained 30% water surplus, although the preparation technically was the same.Therefore, plant residuals left on the field, even in a partly buried status, can trigger water retention.Plant residues protect the soil surface, keep the macropores open, and provide natural nutrition for the macrofauna.Earthworms' life conditions are better under CT, owing to the increased amount of organic matter, more favorable moisture conditions, and reduced disturbance (Figure 4).The highest difference between the two tillage systems was found under the stubble condition.In Central Europe, this period is the most vulnerable from soil erosion aspects, as this is the season of the heaviest extreme thunderstorms while the soil is mostly uncovered [50].An infiltration potential increase of 75% under these circumstances can lead to significant benefits regarding both water retention and soil loss mitigation.In the lower precipitation intensity range, all three soil conditions reacted more or less the same to precipitation intensity increase.This is the case up to 10 mm h −1 regarding CT, whereas under PT, this value is 18 mm h −1 .This suggests that independently from annual variation of single cultivation runs or temporary soil hydraulic conditions, PT has a general long-term effect on infiltration increase, and the earthworms and other soil fauna have a key role in this effect.The calculated soil permeability values also demonstrate these trends (Table 3).CT had higher infiltration values throughout the year, even though most studies have emphasized that infiltration to CT soils increases over the long run, compared to PT [31,49].Interestingly, right after seedbed preparation, CT retained 30% water surplus, although the preparation technically was the same.Therefore, plant residuals left on the field, even in a partly buried status, can trigger water retention.Plant residues protect the soil surface, keep the macropores open, and provide natural nutrition for the macrofauna.Earthworms' life conditions are better under CT, owing to the increased amount of organic matter, more favorable moisture conditions, and reduced disturbance (Figure 4).The highest difference between the two tillage systems was found under the stubble condition.In Central Europe, this period is the most vulnerable from soil erosion aspects, as this is the season of the heaviest extreme thunderstorms while the soil is mostly uncovered [50].An infiltration potential increase of 75% under these circumstances can lead to significant benefits regarding both water retention and soil loss mitigation.

Soil Loss and Erodibility Results for the Small Plots
In general, soil was found to be resistant against erosion.The highest soil loss value measured, due to the most extreme precipitation event, even under the seedbed condition, did not exceed the value of tolerable soil loss [51].The single soil loss results owing to the simulated rainfalls are listed in Table 1.The highest difference was measured under the seedbed soil condition, regardless of the tillage system.This is likely the result of high heterogeneity of soil porosity [52] and an increased value of surface roughness [53].The most endangered soil status was the seedbed condition, soil covered by crops had the highest resistance against soil erosion, and the stubble condition fell between the two.
Theoretically, soil erodibility is a certain value that refers to a soil layer independent of its status, coverage, or porosity.Therefore, predicting soil erodibility is a function of only the chemical and physical properties of the soil layer.The USLE method instructs one to measure erodibility under the seedbed condition [34].The PESERA model (Equation ( 4)) does not take into account soil status; therefore, in this case, the seedbed condition was also selected to be the standard.Thus, the soil erodibility (K factor) was determined for each single precipitation event using both the MUSLE and PESERA methods (Figure 5).

Soil Loss and Erodibility Results for the Small Plots
In general, soil was found to be resistant against erosion.The highest soil loss value measured, due to the most extreme precipitation event, even under the seedbed condition, did not exceed the value of tolerable soil loss [51].The single soil loss results owing to the simulated rainfalls are listed in Table 1.The highest difference was measured under the seedbed soil condition, regardless of the tillage system.This is likely the result of high heterogeneity of soil porosity [52] and an increased value of surface roughness [53].The most endangered soil status was the seedbed condition, soil covered by crops had the highest resistance against soil erosion, and the stubble condition fell between the two.
Theoretically, soil erodibility is a certain value that refers to a soil layer independent of its status, coverage, or porosity.Therefore, predicting soil erodibility is a function of only the chemical and physical properties of the soil layer.The USLE method instructs one to measure erodibility under the seedbed condition [34].The PESERA model (Equation ( 4)) does not take into account soil status; therefore, in this case, the seedbed condition was also selected to be the standard.Thus, the soil erodibility (K factor) was determined for each single precipitation event using both the MUSLE and PESERA methods (Figure 5).However, although both methods calculate soil erodibility, they are not the same in terms of interpretation.Therefore, the values are comparable within methods but not between models.The PESERA method indicated no significant difference between PT and CT erodibility, whereas the soil layer under PT had significantly higher erodibility (p < 0.05) based on the MUSLE method.Both methods calculate soil erodibility based on runoff volume; however, the MUSLE method considers runoff intensity, whereas the PESERA method operates only with slope steepness.Therefore, in the present case, the MUSLE method is more suitable for indicating changes in erodibility due to variation in properties of the same soil.Accordingly, differences, for instance, in aggregate stability can hardly affect erodibility predictions [54].
Fitting a function onto erodibility factor of Equation ( 4), all treatments have different theoretical erodibility results (Figure 6).It is true again that long-term CT decreased the theoretical soil erodibility throughout the year independent of the actual soil status; however, the lowest difference was found under the cover crop condition.This was expected, as crop cover protects the soil surface and therefore mitigates the role of the current physical soil properties in erosion control.On the other hand, the difference highlights the filtering effect of plant residuals left on the surface via CT, as was However, although both methods calculate soil erodibility, they are not the same in terms of interpretation.Therefore, the values are comparable within methods but not between models.The PESERA method indicated no significant difference between PT and CT erodibility, whereas the soil layer under PT had significantly higher erodibility (p < 0.05) based on the MUSLE method.Both methods calculate soil erodibility based on runoff volume; however, the MUSLE method considers runoff intensity, whereas the PESERA method operates only with slope steepness.Therefore, in the present case, the MUSLE method is more suitable for indicating changes in erodibility due to variation in properties of the same soil.Accordingly, differences, for instance, in aggregate stability can hardly affect erodibility predictions [54].
Fitting a function onto erodibility factor of Equation ( 4), all treatments have different theoretical erodibility results (Figure 6).It is true again that long-term CT decreased the theoretical soil erodibility throughout the year independent of the actual soil status; however, the lowest difference was found under the cover crop condition.This was expected, as crop cover protects the soil surface and therefore mitigates the role of the current physical soil properties in erosion control.On the other hand, the difference highlights the filtering effect of plant residuals left on the surface via CT, as was also reported by Alliaume et al. [14].The highest difference was found under the seedbed condition, which likely from the changes in the physical soil properties after cultivation.
Sustainability 2017, 9, 1726 9 of 13 also reported by Alliaume et al. [14].The highest difference was found under the seedbed condition, which likely from the changes in the physical soil properties after cultivation.

Field-Scale Results
At the field-scale, nine natural precipitation events occurred, which triggered soil loss from at least one plot (Table 4).In general, CT plots had lower runoff and soil loss values compared to those of PT plots.At the field scale, CT also had the higher infiltration average; however, the current runoff rate strongly depends on precipitation properties, soil status, and coverage.In general, CT had a one order of magnitude lower runoff rate than PT did during the growing season.
The PESERA k values were also predicted based on the field-scale results (Figure 7).For CT, the predicted k value was smaller than the smallest value at the plot scale (cover crop) and had a fairly good fit, because each precipitation event occurred under the same soil and canopy cover circumstances, i.e., high maize canopy coverage (Table 2).This suggests that, under high canopy

Field-Scale Results
At the field-scale, nine natural precipitation events occurred, which triggered soil loss from at least one plot (Table 4).In general, CT plots had lower runoff and soil loss values compared to those of PT plots.At the field scale, CT also had the higher infiltration average; however, the current runoff rate strongly depends on precipitation properties, soil status, and coverage.In general, CT had a one order of magnitude lower runoff rate than PT did during the growing season.
The PESERA k values were also predicted based on the field-scale results (Figure 7).For CT, the predicted k value was smaller than the smallest value at the plot scale (cover crop) and had a fairly good fit, because each precipitation event occurred under the same soil and canopy cover circumstances, i.e., high maize canopy coverage (Table 2).This suggests that, under high canopy coverage and CT circumstances, net soil loss is a function of runoff volume and-at least within the investigated scales-independent from spatial scales.On the other hand, under PT conditions, some runoff events triggered extremely high soil loss; therefore, the function did not fit (Figure 7).This accords with the statements of Stroosnijder [22], who reported the scaling issue to be one of the biggest difficulties of erosion measurements.The reason is presumably the contribution of rill and gully erosion-as it was observed on PT plots-which can multiply soil loss and greatly increase the sediment concentration in runoff [55].Accordingly, no adequate k value was predictable for this period at field scale under PT.
Sustainability 2017, 9, 1726 10 of 13 investigated scales-independent from spatial scales.On the other hand, under PT conditions, some runoff events triggered extremely high soil loss; therefore, the function did not fit (Figure 7).This accords with the statements of Stroosnijder [22], who reported the scaling issue to be one of the biggest difficulties of erosion measurements.The reason is presumably the contribution of rill and gully erosion-as it was observed on PT plots-which can multiply soil loss and greatly increase the sediment concentration in runoff [55].Accordingly, no adequate k value was predictable for this period at field scale under PT.Both the USLE and MUSLE methods consider the crop and cover management factor (C) and conservation practice factor (P) for calculations [34,38].C and P are determined in a more generalized way-measuring soil losses from plots with various managements and conservation practice-for much coarser spatial resolution; therefore, no adequate values are available for making erodibility (K) predictions in the present case.On the other hand, if formerly calculated (plot scale) MUSLE K values are taken into the USLE model, C and P can be calculated.Therefore, C was calculated for the PT plots at 0.39, which is in line with database values [56].Using this C value for CT plots, P was found to be 0.07, which means 7% soil loss mitigation due to the CT techniques applied.This reduction is in accordance with the results of Vogel et al. [57], who found CT the most effective erosion control tool compared to grassed waterways or contoured buffer strips under maize.Moreover, Panagos et al. [58] could only estimate the effect of stone walls, grassed waterways, and contour tillage on the P factor at the continental scale, because the effect of CT is still unpredictable.Accordingly, the main conservation effect manifests as improvement of soil properties (differences in K) and the role of conservation practice is secondary.

Conclusions
The determination and application of apparent soil infiltration as a function of precipitation intensity seems to be a very useful tool for in situ hydrological modelling.Instead of a single value, such as saturated hydraulic conductivity, it is able to describe process reaction due to changes in the current environmental conditions.It is particularly suitable for making comparisons among tillage systems.Although CT was constructed originally for water retention purposes in arid and semi-arid environments, it has considerable benefits in sub-humid climates as well, both in terms of infiltration improvement and soil erosion mitigation.Theoretically, seedbed preparation is for increased water infiltration and improved porosity conditions, however no benefits of seedbed preparation were found under either tillage system.Contrarily, cover crop application was found to be effective against runoff and soil erosion, even in the second half of the growing season.Conservation practice is very valuable particularly by improving physical soil conditions and so decreasing soil loss.The role of Both the USLE and MUSLE methods consider the crop and cover management factor (C) and conservation practice factor (P) for calculations [34,38].C and P are determined in a more generalized way-measuring soil losses from plots with various managements and conservation practice-for much coarser spatial resolution; therefore, no adequate values are available for making erodibility (K) predictions in the present case.On the other hand, if formerly calculated (plot scale) MUSLE K values are taken into the USLE model, C and P can be calculated.Therefore, C was calculated for the PT plots at 0.39, which is in line with database values [56].Using this C value for CT plots, P was found to be 0.07, which means 7% soil loss mitigation due to the CT techniques applied.This reduction is in accordance with the results of Vogel et al. [57], who found CT the most effective erosion control tool compared to grassed waterways or contoured buffer strips under maize.Moreover, Panagos et al. [58] could only estimate the effect of stone walls, grassed waterways, and contour tillage on the P factor at the continental scale, because the effect of CT is still unpredictable.Accordingly, the main conservation effect manifests as improvement of soil properties (differences in K) and the role of conservation practice is secondary.

Conclusions
The determination and application of apparent soil infiltration as a function of precipitation intensity seems to be a very useful tool for in situ hydrological modelling.Instead of a single value, such as saturated hydraulic conductivity, it is able to describe process reaction due to changes in the current environmental conditions.It is particularly suitable for making comparisons among tillage systems.Although CT was constructed originally for water retention purposes in arid and semi-arid environments, it has considerable benefits in sub-humid climates as well, both in terms of infiltration improvement and soil erosion mitigation.Theoretically, seedbed preparation is for increased water infiltration and improved porosity conditions, however no benefits of seedbed preparation were found under either tillage system.Contrarily, cover crop application was found to be effective against runoff and soil erosion, even in the second half of the growing season.Conservation practice is very valuable particularly by improving physical soil conditions and so decreasing soil loss.The role of better surface conditions is less important from the aspect of soil erosion.To extrapolate the results to various soil types and slope steepness, additional measurements among spatial scales are needed.

Figure 1 .
Figure 1.Location of the study area (solid red square) (A), aerial photo of the Szentgyörgyvár soil erosion experimental site, with the four plots outlined in red (B).CT: conservation tillage, PT: ploughing tillage.

Figure 1 .
Figure 1.Location of the study area (solid red square) (A), aerial photo of the Szentgyörgyvár soil erosion experimental site, with the four plots outlined in red (B).CT: conservation tillage, PT: ploughing tillage.

Figure 2 .
Figure 2. Fitted function on measured runoff values with apparent infiltration intensity (measurements for the seedbed condition under ploughing tillage (a) and conservation tillage (b);Measurement codes refer to the properties in Table1).The calculated apparent infiltration values are highlighted.(y = cumulated runoff; x = irrigation time; R 2 = coefficient of determination).

Figure 2 .
Figure 2. Fitted function on measured runoff values with apparent infiltration intensity (measurements for the seedbed condition under ploughing tillage (a) and conservation tillage (b); Measurement codes refer to the properties in Table 1).The calculated apparent infiltration values are highlighted.(y = cumulated runoff; x = irrigation time; R 2 = coefficient of determination).

Figure 3 .
Figure 3. Final apparent infiltration rate changes as a function of irrigation intensity fitted using Equation (2); (a) conservation tillage; (b) ploughing tillage.(Legend is for both subfigures).

Figure 3 .
Figure 3. Final apparent infiltration rate changes as a function of irrigation intensity fitted using Equation (2); (a) conservation tillage; (b) ploughing tillage.(Legend is for both subfigures).

Figure 4 .
Figure 4. Effect of tillage on earthworm number (a) and weight (b) at Szentgyörgyvár in October 2016 (CT: conservation tillage, PT: ploughing tillage) for n = 18, Differences between the tillage systems are significant at the p < 0.01 level.

Figure 4 .
Figure 4. Effect of tillage on earthworm number (a) and weight (b) at Szentgyörgyvár in October 2016 (CT: conservation tillage, PT: ploughing tillage) for n = 18, Differences between the tillage systems are significant at the p < 0.01 level.

Figure 5 .
Figure 5. Soil erodibility values of the seedbed condition soils determined by rainfall simulation, calculated according to the Pan-European Soil Erosion Risk Assessment[33] and Modified Universal Soil Loss Equation[38] methods; CT: conservation tillage, PT: ploughing tillage.

Figure 5 .
Figure 5. Soil erodibility values of the seedbed condition soils determined by rainfall simulation, calculated according to the Pan-European Soil Erosion Risk Assessment[33] and Modified Universal Soil Loss Equation[38] methods; CT: conservation tillage, PT: ploughing tillage.

Figure 6 .
Figure 6.Rainfall simulation soil loss data as a function of runoff volume (V) and slope steepness (S).Based on Equation (4), the steepness of the fitted function represents soil erodibility (k).(a) conservation tillage and (b) ploughing tillage.(y = soil loss; x = V 2 S 1.5 ; R 2 = coefficient of determination).

Figure 6 .
Figure 6.Rainfall simulation soil loss data as a function of runoff volume (V) and slope steepness (S).Based on Equation (4), the steepness of the fitted function represents soil erodibility (k).(a) conservation tillage and (b) ploughing tillage.(y = soil loss; x = V 2 S 1.5 ; R 2 = coefficient of determination).

Figure 7 .
Figure 7. Field-scale soil loss data owing to natural rainfalls as a function of runoff volume (V) and slope steepness (S).Based on Equation (4), the steepness of the fitted function represents soil erodibility (k).Scales are logarithmic.(y = soil loss; x = V 2 S 1.5 ; R 2 = coefficient of determination).

Figure 7 .
Figure 7. Field-scale soil loss data owing to natural rainfalls as a function of runoff volume (V) and slope steepness (S).Based on Equation (4), the steepness of the fitted function represents soil erodibility (k).Scales are logarithmic.(y = soil loss; x = V 2 S 1.5 ; R 2 = coefficient of determination).

Table 1 .
Cultivation activities used at the studied experimental sites.

Table 1 .
Cultivation activities used at the studied experimental sites.

Table 4 .
Precipitation, runoff, and soil loss properties measured on 1200 m 2 field plots at the Szentgyörgyvár experimental site during 2016.I30 is the highest 30 min rainfall intensity during the precipitation; CT: conservation tillage; PT: ploughing tillage; cc = concentration.

Table 4 .
Precipitation, runoff, and soil loss properties measured on 1200 m 2 field plots at the Szentgyörgyvár experimental site during 2016.I 30 is the highest 30 min rainfall intensity during the precipitation; CT: conservation tillage; PT: ploughing tillage; cc = concentration.