Runoff and Soil Erosion Assessment on Forest Roads Using a Small Scale Rainfall Simulator

Forestry operations can significantly alter hydrological and erosional processes in a catchment. In the course of developing timberland, a network of persistent roads and skid trails causing soil compaction is usually established. Hereby, the infiltration rate of the soil is distinctly reduced, which leads to the generation of overland flow—this may also cause soil erosion. In this study, a small-scale rainfall simulator is used to investigate hydrological and erosional processes on forest roads and skid trails. The results show increased runoff rates on forest roads, up to 25 times higher than on undisturbed forest topsoil. On skid trails, the runoff rates were altered especially in rutted areas (16 times higher) while unrutted parts showed a lesser change (four times higher). With sufficient overland flow, soil erosion rates also rose, particularly when the vegetation cover of the surface was removed: bare road surfaces featured higher mean erosion rates (195 g ̈ m ́2) than partly or completely vegetated skid trails (13 g ̈ m ́2) and undisturbed sites (5 g ̈ m ́2). The findings presented in this study indicate the need for the use of compaction reducing technology during forestry operations and a revegetation of road surfaces in order to minimize the detrimental factor of roads and skid trails on water retention and soil conservation.


Introduction
Forestry operations usually go along with the construction of a more or less dense forest road network. These forest roads feature soil properties that differ significantly from adjacent undisturbed topsoil. In the course of establishing a sufficient trafficability, different actions are taken to ensure a stable and durable road surface. The most commonly realized method of forest road construction is a simple grading, which produces a solid surface out of the in-situ soil without the addition of other construction material. If there is a need for a more stable road surface, for instance because of heavy traffic or on steep road segments, an additional fortification is often used. In most cases, a subgrade made out of chippings or gravel is built to increase the surface stability [1][2][3].
In both cases, the soil bulk density is increased and the pore system is heavily disturbed. Because of this soil compaction, infiltration rates and capacities are significantly reduced, causing a much higher tendency towards the generation of Hortonian overland flow (HOF), when rainfall intensities exceed the infiltration rate of the road surface [4][5][6][7].
This process-related interaction also occurs on skid trails, where no grading or fortification takes place. Here, the soil compaction is caused by the direct overpassing of previously undisturbed topsoil with forestry equipment. Therefore, especially the rutted areas of the skid trail show a severe soil compaction, whereas the unrutted areas are less affected [8,9].
Beyond that, the influence of traffic induced soil compaction on skid trails is a comparatively long-lasting effect. An increased soil bulk density was measured in former rutted parts of skid trails The rainfall simulations were carried out in two different forested headwater catchments in Rhineland-Palatinate, Western Germany, located in the nature reserve Laacher See (50°25′47.4″N; 7°15′49.1″E) and the national park Hunsrück-Hochwald (49°43′47.0″N; 7°5′23.0″E) (Figure 1).
Both catchments are situated within the Rhenish Massif and feature mostly slate and quartzite bedrock with dominant quartzite in the national park Hunsrück-Hochwald. As the Laacher See catchment is located in the East-Eifel volcanic field and in direct vicinity of the Laacher See caldera and numerous quaternary basaltic volcanoes, Laacher See tephra (LST) and basaltic vulcanites dominate its local geology. Nevertheless, both study areas show similar USDA soil textures (Table  1), basing upon a total of 447 soil samples, which were taken at the rainfall simulation sites and during several soil sampling campaigns in both catchments. Regarding their land use, the catchments are characterized by forested land, which, as a consequence, involves frequent traffic on dense road and skid trail networks.

The Portable Rainfall Simulator
A nozzle type rainfall simulator was designed based on the structural design of [27] for measuring interrill erosion ( Figure 2). The nozzle was attached to a 1 m² square aluminum frame at a height of 2 m. Four adjustable telescoping aluminum legs allowed a horizontal alignment of the frame and a windshield was attached to the legs to minimize the influence of wind during the experiments. The water flow was driven by a battery powered electric bilge pump, while the flow rate was controlled using a small gate valve. Artificial raindrops were generated using a Lechler

The Portable Rainfall Simulator
A nozzle type rainfall simulator was designed based on the structural design of [27] for measuring interrill erosion ( Figure 2). The nozzle was attached to a 1 m 2 square aluminum frame at a height of 2 m. Four adjustable telescoping aluminum legs allowed a horizontal alignment of the frame and a windshield was attached to the legs to minimize the influence of wind during the experiments. The water flow was driven by a battery powered electric bilge pump, while the flow rate was controlled using a small gate valve. Artificial raindrops were generated using a Lechler 406.608 full cone nozzle, which is stated to produce a good replication of natural rainfall characteristics in numerous publications [27][28][29]. As a second flow control mechanism, a manometer was installed near the nozzle to ensure a constant water pressure at the drop height of 2 m. Hydrology 2016, 3, 25 4 of 22 406.608 full cone nozzle, which is stated to produce a good replication of natural rainfall characteristics in numerous publications [27][28][29]. As a second flow control mechanism, a manometer was installed near the nozzle to ensure a constant water pressure at the drop height of 2 m. The observed plot surface was limited by a 0.64 m² square frame made of stainless steel, which was driven into the topsoil. It also featured a runoff apron that allowed a safe sampling of generated overland flow and eroded soil by attaching wide neck bottles to the outlet.

Calibration of the Rainfall Simulator
To ensure that the simulated rainfall replicates natural rainfall characteristics as closely as possible, a wide range of calibration methods was used. First of all, the spatial rainfall distribution and the intensity were measured to identify an operating pressure with an acceptable trade-off of both spatial homogeneity and desired intensity. The low drop height of 2 m conditioned a diminished terminal fall velocity of the artificial raindrops compared to natural rainfall. Therefore, a comparatively high intensity >40 mm•h −1 was planned, to ensure kinetic energy values high enough to mobilize soil particles, as the rainfall kinetic energy is a product of mass (respectively intensity) and fall velocity [27,60].
An increase of intensity was obtained by reducing the operation pressure which caused on the other hand a more and more inhomogeneous rainfall distribution because of the nozzle specifications: Lechler specifies a minimum pressure of 0.5 bar for operating the type-460 nozzles with a symmetrical spray cone. Thus, only a narrow range of viable pressures was given to produce a homogenous and erosive rainfall. To measure both intensity and spatial homogeneity, 263 collectors were set up beneath the rainfall simulator, covering the base area of 1 m 2 . After a rainfall exposure of 15 min, each collector was weighed to determine the quantity of the collected rainfall at the respective The observed plot surface was limited by a 0.64 m 2 square frame made of stainless steel, which was driven into the topsoil. It also featured a runoff apron that allowed a safe sampling of generated overland flow and eroded soil by attaching wide neck bottles to the outlet.

Calibration of the Rainfall Simulator
To ensure that the simulated rainfall replicates natural rainfall characteristics as closely as possible, a wide range of calibration methods was used. First of all, the spatial rainfall distribution and the intensity were measured to identify an operating pressure with an acceptable trade-off of both spatial homogeneity and desired intensity. The low drop height of 2 m conditioned a diminished terminal fall velocity of the artificial raindrops compared to natural rainfall. Therefore, a comparatively high intensity >40 mm¨h´1 was planned, to ensure kinetic energy values high enough to mobilize soil particles, as the rainfall kinetic energy is a product of mass (respectively intensity) and fall velocity [27,60].
An increase of intensity was obtained by reducing the operation pressure which caused on the other hand a more and more inhomogeneous rainfall distribution because of the nozzle specifications: Lechler specifies a minimum pressure of 0.5 bar for operating the type-460 nozzles with a symmetrical spray cone. Thus, only a narrow range of viable pressures was given to produce a homogenous and erosive rainfall. To measure both intensity and spatial homogeneity, 263 collectors were set up beneath the rainfall simulator, covering the base area of 1 m 2 . After a rainfall exposure of 15 min, each collector was weighed to determine the quantity of the collected rainfall at the respective plot position. With the aid of the given collector area, the rainfall intensity was calculated. The obtained dataset was used to interpolate the rainfall intensity across the plot and beyond that to calculate the spatial homogeneity using the Christiansen-Coefficient (CU) presented in [61]. Interpolations were conducted with the help of a variogram-based kriging, using Golden Software Surfer 8.
An operation pressure of 0.09 bar, resulting in a mean intensity of 45 mm¨h´1 (equals 0.48 L¨min´1 applied to the plot surface of 0.64 m 2 ) showed the best results of both CU and desired intensity ( Table 2). Using the determined operation pressure, the drop size distribution and fall velocity was measured using a Laser Precipitation Monitor (LPM) by Thies© (Göttingen, Germany) which allowed a contactless real-time measurement at single-drop resolution. The LPM was placed at 9 different positions beneath the rainfall simulator, with a measurement time of 15 min per position. Amongst others, the most important rainfall characteristics, median drop diameter (d50) and kinetic rainfall energy (KE) were derived from the LPM-measurements. The latter were calculated based on [27,60].

Experimental Setup
The rainfall characteristics during the simulation runs are presented in Table 2.
The chosen intensity represents heavy rainfall events of 30 min with a return period of approximately 5 years for the NPHH catchment [5]. Due to a lack of reliable data for the LS catchment, the return period can only be assumed. The catchment shows basically comparable precipitation patterns, with the only difference that convective precipitation events during the summer months occur more frequently [62]. As those events usually exhibit high rainfall intensities, a slightly shorter return period than 5 years seems to be probable.
Each run consisted of a sequence of three rainfall events with a duration of 30 min separated by a pause of 5 min, similar to the rainfall pattern presented in [4], where a comparable experimental setup regarding plot characteristics and research question was at hand. The separations allowed evaluating the changing infiltration patterns at the beginning of each event caused by the gradually higher antecedent soil moisture. Infiltration volumes were derived indirectly by subtracting measured runoff volumes from the amount of rainfall input in the same time interval.
The rainfall intensity for each simulation run was quantified at the beginning and the end of every 30-min event by measuring the output of the simulator directly at the nozzle. Additionally, the operating pressure of 0.09 bar was constantly checked throughout the rainfall events. Because of the fact that in-situ experiments are difficult to control perfectly, a certain degree of inconstant intensity was expected. To overcome those uncertainties, experiments on forest roads were used to calculate possible deviations. A complete saturation of the compacted topsoil-and therefore saturation flow conditions with P « Q-was assumed, if: (a) the coefficient of variation (CV) of the runoff volumes in the last 5 min of the experiment was <5%, indicating stable runoff conditions; and (b) the mean runoff sum in the last 5 min of the experiment was 100%˘10% of the desired intensity. A similar approach was used by [4] to address saturated conditions during rainfall simulations. After this, the mean runoff sum of the same time interval was calculated to obtain the range of potential deviation from the targeted rainfall intensity.
Overland flow and sediment yield were collected in a high temporal resolution of one minute. Therefore, it was possible to generate precise datasets of hydrological characteristics and soil erodibility over a rainfall event. Foremost, runoff volumes for each minute of the simulated rainfall events were measured directly. With known rainfall intensity, it was also possible to calculate the runoff-coefficient (C), the percentage of runoff related to the amount of rainfall for any desired time-interval greater or equal 1 min. Additionally, the suspended sediment load was measured for each minute of the simulations. Detailed sediment graphs with a high temporal resolution of one minute can be found in [63]. As eroded particles often settled out of the flowing water at the plot apron and were collected after each simulation run, only the single event based values with a temporal resolution of 30 min are presented within this article. All of the obtained results are displayed normalized for an area of 1 m 2 in order to allow a better comparability with other studies using the same normalization [4,5,[33][34][35].
To elaborate the differences between each road type, a road categorization was established. A first distinction was made between forest roads and skid trails, as there are major differences concerning their constructional layout and persistence ( Figure 3). minute of the simulations. Detailed sediment graphs with a high temporal resolution of one minute can be found in [63]. As eroded particles often settled out of the flowing water at the plot apron and were collected after each simulation run, only the single event based values with a temporal resolution of 30 min are presented within this article. All of the obtained results are displayed normalized for an area of 1 m 2 in order to allow a better comparability with other studies using the same normalization [4,5,[33][34][35].
To elaborate the differences between each road type, a road categorization was established. A first distinction was made between forest roads and skid trails, as there are major differences concerning their constructional layout and persistence ( Figure 3). Forest roads were subdivided into fortified (graveled) and unfortified (graded) roads, as the absence of a base course in case of unfortified roads was believed to have a noticeable influence on hydrological and erosional processes. Because of the small plot size, even rutted and unrutted areas on the road surface could be distinguished. Especially skid trail surfaces were divided into rutted and unrutted areas as the main traffic induced soil compaction was assumed to be present in forestry equipment tracks. This assumption was proved right, as the results obtained on persistent forest roads featured only minor differences between rutted and unrutted areas. Thus, the respective results are displayed aggregated for both rutted and unrutted parts of the road. In contrast, soil compaction conditioned a clear pattern on skid trails: Here, especially the rutted parts showed higher rates of runoff generation than unrutted areas. As a consequence, both are displayed separately.
Regarding the plot characteristics (Table 3), the soil moisture of the plot was measured before and after the simulation run by collecting sample ring probes with a volume of 0.1 dm 3 , which were taken in direct vicinity of (before) and on the plot surface (after). In the course of this sampling the bulk density (ρ) and gravel content could also be derived. Additionally, the vegetation cover of the plot surface was determined, as it is a major control factor for soil erosion processes [30,35,[64][65][66]. For this purpose, digital photographs of the plot surfaces were analyzed by marking covered areas with a mask tool and calculating the percentage of vegetation covered pixels. At last, hydrophobicity of the topsoil was measured using a classification basing on the Water Drop Penetration Time (WDPT) method as presented in [5,67], since soil hydrophobicity may additionally trigger overland flow generation [5]. The WDPT method evaluates the hydrophobicity on the basis of how long it takes to infiltrate a single drop of water that is applied to the surface. WDPT classes range from 0 (no noticeable water repellency, infiltration time <5 s) to 6 (severe water repellency, infiltration time >300 s). Forest roads were subdivided into fortified (graveled) and unfortified (graded) roads, as the absence of a base course in case of unfortified roads was believed to have a noticeable influence on hydrological and erosional processes. Because of the small plot size, even rutted and unrutted areas on the road surface could be distinguished. Especially skid trail surfaces were divided into rutted and unrutted areas as the main traffic induced soil compaction was assumed to be present in forestry equipment tracks. This assumption was proved right, as the results obtained on persistent forest roads featured only minor differences between rutted and unrutted areas. Thus, the respective results are displayed aggregated for both rutted and unrutted parts of the road. In contrast, soil compaction conditioned a clear pattern on skid trails: Here, especially the rutted parts showed higher rates of runoff generation than unrutted areas. As a consequence, both are displayed separately.
Regarding the plot characteristics (Table 3), the soil moisture of the plot was measured before and after the simulation run by collecting sample ring probes with a volume of 0.1 dm 3 , which were taken in direct vicinity of (before) and on the plot surface (after). In the course of this sampling the bulk density (ρ) and gravel content could also be derived. Additionally, the vegetation cover of the plot surface was determined, as it is a major control factor for soil erosion processes [30,35,[64][65][66]. For this purpose, digital photographs of the plot surfaces were analyzed by marking covered areas with a mask tool and calculating the percentage of vegetation covered pixels. At last, hydrophobicity of the topsoil was measured using a classification basing on the Water Drop Penetration Time (WDPT) method as presented in [5,67], since soil hydrophobicity may additionally trigger overland flow generation [5]. The WDPT method evaluates the hydrophobicity on the basis of how long it takes to infiltrate a single drop of water that is applied to the surface. WDPT classes range from 0 (no noticeable water repellency, infiltration time <5 s) to 6 (severe water repellency, infiltration time >300 s).

Statistical Analysis
The plot characteristics and the category based mean values of runoff and eroded sediment were described statistically according to [5], as this publication represents one of the few existent studies that provide datasets for runoff generation and soil erosion on skid trails in a European catchment. IBM SPSS 23 was used to calculate the required statistical parameters. With the help of the Kolmogorov-Smirnov test, it was confirmed, that the datasets were not normally distributed. Therefore, only non-parametric tests were performed. In respect thereof, the same statistical basic conditions as in [5] were on hand. Then, boxplots for each surface category were plotted and the Kruskal-Wallis test was performed for independent samples to verify if the categories showed significant differences. Finally, using the surface categories, a paired Mann-Whitney test was conducted to show significant differences between them. As the Mann-Whitney test focuses mainly on the median of the dataset, it was suitable for the rather small number of simulations for each surface category [68]. To test if runoff and erosion results and certain plot characteristics correlate, Spearman-Rho correlation matrices were also calculated.  Table 4 shows the runoff values of each simulation, separated by the single rainfall events. Additionally, the calculated runoff-coefficients are specified.

Overland Flow
Detailed hydrographs were plotted as a result of the high-resolution sampling method. Exemplary hydrographs basing upon simulation runs on different surfaces (simulations #2, #6, #10, #14 and #19) give an insight in unprocessed datasets ( Figure 4). They show typical runoff patterns of the respective surface category.
Both fortified and unfortified forest roads show runoff rates that are significantly different compared to undisturbed forest topsoil. After only a few minutes into the first rainfall event, both surfaces show an almost instantaneous shift from no measurable runoff to flow rates in the range of the rainfall input. By comparison, the unfortified road tends to generate high runoff amounts earlier during the first minutes of each rainfall event.
The rutted skid trail shows the most apparent pattern of successive topsoil saturation. During rainfall Event 1, the first runoff is already measured after 3 min-presumably because of the hydrophobicity of the soil (WDPT class 2)-and an increased amount of runoff generation starts at 7 min. Towards the end of the second rainfall event, the rutted skid trail already shows runoff comparable to the rainfall input. A fast saturation at the beginning of Event 3 is clearly visible.
After this response time, runoff volumes level off at rainfall input rates, representing an assumed 100% saturation of the topsoil leading to saturation flow conditions.   Both fortified and unfortified forest roads show runoff rates that are significantly different compared to undisturbed forest topsoil. After only a few minutes into the first rainfall event, both surfaces show an almost instantaneous shift from no measurable runoff to flow rates in the range of the rainfall input. By comparison, the unfortified road tends to generate high runoff amounts earlier during the first minutes of each rainfall event.
The rutted skid trail shows the most apparent pattern of successive topsoil saturation. During rainfall Event 1, the first runoff is already measured after 3 min-presumably because of the hydrophobicity of the soil (WDPT class 2)-and an increased amount of runoff generation starts at 7 min. Towards the end of the second rainfall event, the rutted skid trail already shows runoff On the unrutted skid trail, no runoff was measured until Minute 18, which indicates a comparatively high infiltration rate of the topsoil. Even the following rainfall events show a distinct saturation stage, as there are 4 (Event 2) and 6 (Event 3) min without any runoff at the beginning of each event. The undisturbed forest soil shows a stable but almost negligible overland flow. There is no noticeable pattern of rising runoff rates, which shows that even after 90 min rainfall with a rainfall intensity of 45 mm¨h´1 and 43.6 L rainfall applied to the plot surface, the infiltration capacity is not exceeded. The measured runoff volumes represent only a small amount of HOF that occurs because of the high rainfall intensity.
The hydrographs also show an example for calculating the possible deviation of rainfall intensity during the simulations: There are data points that indicate more runoff than applied rainfall, especially when looking at the third event on fortified and unfortified roads. For example, the deviation for simulation #2, conducted on a fortified road is displayed in Figure 4 as grey area above the indicated amount of rainfall. The CV of the runoff in the last 5 min was 1.2%, implying stable runoff conditions with a mean C of 107%-therefore, a deviation of 7% was assumed.
The largest deviation in the whole dataset was observed during simulation #5 with a mean C of 109% and a CV of 2.2% during Minutes 85-90. Because of that, a possible general error of 10% was accounted for all of the experiments and integrated in Figure 6 which shows the aggregated volumes of runoff and infiltration for each surface category. Nevertheless, all of the simulations showed a total C of <100% (Table 4). Besides the possibility of deviating rainfall intensities, a partial retention of the superficially flowing water caused by irregularities on the plot surface, for example imprinted tire treads, may also have caused higher C-values towards the end of the simulations. In those imprints, small pools of water were built up, which were drained when the structures were eroded by splash and flowing water. Figure 5 provides an example observed on a rutted forest road. treads, may also have caused higher C-values towards the end of the simulations. In those imprints, small pools of water were built up, which were drained when the structures were eroded by splash and flowing water. Figure 5 provides an example observed on a rutted forest road. The described runoff patterns also apply to the mean values of each category. These datasets are displayed as processed datasets, the mean percentages of runoff and infiltration for each surface category with a time interval of 10 min (Figure 6). The described runoff patterns also apply to the mean values of each category. These datasets are displayed as processed datasets, the mean percentages of runoff and infiltration for each surface category with a time interval of 10 min ( Figure 6).
Again, forest roads produce the highest runoff of all surface types. Only during the first 10 min of rainfall Event 1 both surfaces show higher infiltration than runoff. Hereafter, the mean runoff volumes are permanently higher. When looking at the mean runoff volumes of all simulation runs (Table 5), the fortified surfaces show lower runoff volumes than unfortified roads. Both surfaces are compacted but there is a constructional difference between them that causes the deviating runoff patterns. Most likely it can be assumed that the base course featured on fortified roads led to, on average, larger soil pores compared to graded unfortified roads, as the coarse material stabilizes the matrix and therefore the pore system to a certain degree. Fine material between the crushed stones is not compacted that much compared to graded surfaces where the whole topsoil is affected thoroughly and an almost complete loss of pore volume has to be stated. Additionally, fortified roads feature higher percentages of sand, presumably because of the insertion of extraneous material, which again leads to a higher porosity of the topsoil. The described runoff patterns also apply to the mean values of each category. These datasets are displayed as processed datasets, the mean percentages of runoff and infiltration for each surface category with a time interval of 10 min (Figure 6).  The skid trails showed comparatively great differences between rutted and unrutted areas. According to the findings presented on the basis of Figure 4, the saturation of the soil successively proceeds over the three rainfall events leading to only minor infiltration volumes in case of rutted surfaces.
Here, the mean runoff percentage of rainfall Event 3 exceeds the values of rainfall Event 1 on forest roads. It is assumable, that the saturation of the topsoil would proceed and eventually cause a runoff pattern similar to forest roads if the duration of the simulated rainfall was prolonged. As the road surfaces represent a certain reference point for the generation of HOF over time, an almost completely saturated soil has to be expected within additional 60 min. This tendency is not expectable for unrutted surfaces, as they merely reach runoff percentages of 25% within rainfall Event 3.
Undisturbed topsoil shows only minimal runoff percentages. Also, the variation between each rainfall event is minimal. This indicates a stable infiltration rate, which leads to almost no notable HOF; therefore, the topsoil was not saturated.
In addition to the study of the temporal changes in runoff and infiltration, a statistical analysis of the mean values for each category provides a better insight and allows addressing significant differences between surface types. Figure 7 shows the boxplots for each surface category, Table 5 the mean total runoff volumes and resulting runoff coefficients.
As expected, both forest road types show high C-values because of soil compaction resulting from constructional measures. Comparatively fast topsoil saturation in case of graded road surfaces is reflected in a slightly higher runoff-coefficient. The most noticeable result was achieved on rutted skid trails. Even though no fortification was existent and only the weight of forestry equipment caused a soil compaction, the mean runoff-coefficient lies clearly above the mean C-values of undisturbed forest topsoil, showing a 16-times higher percentage of overland flow. This increase does not only occur on rutted skid trails but also on unrutted parts. Here, the mean runoff-coefficient is noticeably lower than on rutted surfaces but it is still four times higher than the mean values of undisturbed soil, indicating a compacted soil.
forest roads. It is assumable, that the saturation of the topsoil would proceed and eventually cause a runoff pattern similar to forest roads if the duration of the simulated rainfall was prolonged. As the road surfaces represent a certain reference point for the generation of HOF over time, an almost completely saturated soil has to be expected within additional 60 min. This tendency is not expectable for unrutted surfaces, as they merely reach runoff percentages of 25% within rainfall Event 3.
Undisturbed topsoil shows only minimal runoff percentages. Also, the variation between each rainfall event is minimal. This indicates a stable infiltration rate, which leads to almost no notable HOF; therefore, the topsoil was not saturated.
In addition to the study of the temporal changes in runoff and infiltration, a statistical analysis of the mean values for each category provides a better insight and allows addressing significant differences between surface types. Figure 7 shows the boxplots for each surface category, Table 5 the mean total runoff volumes and resulting runoff coefficients.   In an additional step, differences between surface categories were tested regarding their significance ( Table 6). As expected, the Mann-Whitney test showed significant differences between among road categories, unrutted skid trails and undisturbed topsoil. Apart from that, rutted skid trails showed no significant difference from forest roads, indicating the increased runoff volumes due to topsoil compaction. In contrast, unrutted skid trails show a strong similarity with undisturbed topsoil. To investigate, if plot characteristics influenced the runoff generation, Spearman-Rho coefficients were calculated (Table 7).
There are identifiable significant correlations between the mean runoff sum and cover, gravel content and WDPT class when observing all conducted experiments. These strong correlations have to be relativized, as vegetation cover, gravel content and WDPT class differ depending on the constructional layout of the road surface. To avoid a spurious relationship, it is helpful to consider category-based values. Because of the rather small population in each surface category, only few significant correlations were achieved. The most noticeable correlations were found for rutted skid trails; here, a strong negative correlation between vegetation cover and mean runoff volume is existent. A strong positive correlation between WPDT class and mean runoff volume was also confirmed, indicating higher runoff volumes with severer water repellency. Compared to existent studies discussing the influence of roads and skid trails in Central Europe, the findings presented in this study show similar results. Rainfall simulations with a comparable intensity conducted by [5] showed a median runoff-coefficient of 78.6% on skid trails and harvester tracks. In contrast, a median of only 0.2%-1.3% was measured on undisturbed topsoil. Therefore, an even higher influence of skid trails was determined, making the presented results appear plausible-related results were also confirmed by [15,32]. Along with skid trails, forest roads are characterized as a potential source of overland flow, endangering the water retention potential in forested catchments [69]. This finding could be confirmed by the conducted rainfall experiments. Reported increases of overland flow due to forest roads and soil compaction on skid trails for catchments outside Central Europe and obtained by rainfall simulations with divergent intensities also confirm these findings [4,6,17,24,26,[33][34][35]55]: Here, for example, a faster generation of overland flow on forest roads was proven, most likely caused by soil compaction [4,[33][34][35]. Despite covering only a small percentage of a catchment, road surfaces were accounted for a runoff increase of 10% in selected catchments [6]. Especially the parts of the road lacking a vegetation cover were addressed as areas showing noticeably higher runoff coefficients [33][34][35]. One of the few studies that explicitly compare undisturbed sites with skid trails showed with the help of comparably large runoff plots (2ˆ15 m and 2ˆ10 m) that skid trails generated a mean of 11.99 L¨m´2 over a period of two months. In contrast, undisturbed sites only generated 0.13 L¨m´2 [55].

Soil Erosion
As a consequence of the generated overland flows, the expected soil erosion processes developed. Table 8 shows the eroded sediment sums for every simulation and rainfall event.
Regarding the mean total eroded sediment for each rainfall event, it is clearly visible, that forest roads showed the highest erosion rates during the simulation runs, whereas on skid trail surfaces and undisturbed topsoil virtually negligible amounts of sediment were eroded (Figure 8).
Forest roads generated the highest sums of eroded sediment because of two main reasons: On the one hand, an overland flow sufficient for erosional processes was generated in the first place. On the other hand, the vegetation cover was removed due to grading and compaction of the road surface. Beyond that, pioneer flora that could trigger an ecological succession was repeatedly removed because of the continuous use of the road by forestry vehicles.
The bare surfaces are much more vulnerable to erosion processes induced by rainfall and overland flow than skid trail surfaces and forest topsoil, where interception by vegetation and forest floor occurs.
Thus, not only the incoming rainfall is retained; since the first impact of the raindrops occurs on understory vegetation or litter, the effect of splash erosion is significantly diminished. Regarding the mean total eroded sediment for each rainfall event, it is clearly visible, that forest roads showed the highest erosion rates during the simulation runs, whereas on skid trail surfaces and undisturbed topsoil virtually negligible amounts of sediment were eroded (Figure 8).
Forest roads generated the highest sums of eroded sediment because of two main reasons: On the one hand, an overland flow sufficient for erosional processes was generated in the first place. On the other hand, the vegetation cover was removed due to grading and compaction of the road surface. Beyond that, pioneer flora that could trigger an ecological succession was repeatedly removed because of the continuous use of the road by forestry vehicles.
The bare surfaces are much more vulnerable to erosion processes induced by rainfall and overland flow than skid trail surfaces and forest topsoil, where interception by vegetation and forest floor occurs. Thus, not only the incoming rainfall is retained; since the first impact of the raindrops occurs on understory vegetation or litter, the effect of splash erosion is significantly diminished.
Erosion rates on forest roads were significantly higher than those measured on undisturbed soil and skid trails. Mean sediment sums of fortified and unfortified roads were 36-and 68-times higher than on forest topsoil, respectivly. In the course of the simulation, forest roads showed erosion rates that even exceeded literature values for agricultural land in Central Europe [70].
A closer look at the different surface types reveals that fortified forest roads produced an almost constant sediment yield during all rainfall events, while the unfortified road surfaces yielded decreasing amounts of sediment (Table 9).  Erosion rates on forest roads were significantly higher than those measured on undisturbed soil and skid trails. Mean sediment sums of fortified and unfortified roads were 36-and 68-times higher than on forest topsoil, respectivly. In the course of the simulation, forest roads showed erosion rates that even exceeded literature values for agricultural land in Central Europe [70].
A closer look at the different surface types reveals that fortified forest roads produced an almost constant sediment yield during all rainfall events, while the unfortified road surfaces yielded decreasing amounts of sediment (Table 9). Table 9. Mean sums of sediment eroded from different surface categories during three consecutive rainfall events.

Surface Category
Mean Σ Sediment (g¨m´2)

Event 1 Event 2 Event 3 Total
Forest road, fortified 37 This temporal pattern is directly linked to the missing fortification on graded roads. Here, the surface is compacted to ensure the trafficability, but the lack of a base course that stabilizes the whole roadway allows tire treads to dig deep into the surface and to detach entire soil aggregates. This material is now prone to soil erosion processes as it is no longer embedded in the compacted matrix. During the three rainfall events, recently detached particles were gradually eroded, leading to a more and more decreasing erosion rate as the percentage of already transportable soil was diminished and particles at first had to be detached by splash and flowing water.
In contrast, fortified roads showed a comparatively high soil erosion rate because of the bare surface, but because of the stabilized roadbed there was only a certain kind of base value that could be eroded from the topmost soil layer. In this case, constant rainfall kinetic energy and intensity limited the amount of detached and eroded sediment. It is assumable, that unfortified roads would eventually show similar erosion rates, once all of the previously traffic-detached material is washed off the plot surface.
Another simultaneously observable process that potentially limits the soil erosion rates over time is selective soil erosion. Splash detachment and particle transport by flowing water particularly affects high erodible grain sizes like fine sand and silt. Therefore, these grain fractions are transported away from the surface leaving mainly coarse fractions like larger sand and gravel. Over time, these coarse grains function as an armoring layer, a stone mulch preventing the erosion of finer soil particles beneath them or in their immediate vicinity. If the kinetic energy of the applied rainfall is insufficient for detaching the armoring layer, a decreasing soil erosion rate has to be expected. After 90 min of artificial rainfall, however, an enrichment of coarse material was observed (Figure 9).
A concluding analysis of the particle-size-distribution using sedimentation was conducted to specify the soil properties of the irrigated plots and to evaluate the possible differences between the initial particle-size-distribution and the grain size composition of the eroded material. As skid trails and forest soils did not generate sufficient erosion to obtain enough material for sedimentation analysis, only samples of forest roads were examined. Figure 10 features the results of both plot surfaces and mean values for different surfaces. On the basis of the grain-size-distribution, it is clear to see that the unfortified, graded roads show soil textures akin to skid trails and undisturbed topsoil. Both were classified as sandy loam according to USDA soil taxonomy. Again, this indicates that no alteration of the soil texture by any fortification was existent. In contrast to that, fortified roads showed higher percentages of sand, even leading to a classification as another soil type, loamy sand.
Compared to these results, fortified and unfortified roads showed almost identical grain-size-distributions of the eroded material. Both mean soil textures are distinctly finer than their reference specimen. This result is related to the constant rainfall characteristics, which limited the maximum grain size erodible by splash and flowing water. Therefore, both particle-size distributions approximate each other. In this connection, the mean results provide a strong evidence for constant rainfall characteristics during all simulated rainfall events. Beyond that, a selective erosion of fine particles is evident.
eventually show similar erosion rates, once all of the previously traffic-detached material is washed off the plot surface.
Another simultaneously observable process that potentially limits the soil erosion rates over time is selective soil erosion. Splash detachment and particle transport by flowing water particularly affects high erodible grain sizes like fine sand and silt. Therefore, these grain fractions are transported away from the surface leaving mainly coarse fractions like larger sand and gravel. Over time, these coarse grains function as an armoring layer, a stone mulch preventing the erosion of finer soil particles beneath them or in their immediate vicinity. If the kinetic energy of the applied rainfall is insufficient for detaching the armoring layer, a decreasing soil erosion rate has to be expected. After 90 min of artificial rainfall, however, an enrichment of coarse material was observed (Figure 9). A concluding analysis of the particle-size-distribution using sedimentation was conducted to specify the soil properties of the irrigated plots and to evaluate the possible differences between the initial particle-size-distribution and the grain size composition of the eroded material. As skid trails and forest soils did not generate sufficient erosion to obtain enough material for sedimentation analysis, only samples of forest roads were examined. Figure 10 features the results of both plot surfaces and mean values for different surfaces. On the basis of the grain-size-distribution, it is clear to see that the unfortified, graded roads show soil textures akin to skid trails and undisturbed topsoil. Both were classified as sandy loam according to USDA soil taxonomy. Again, this indicates that no alteration of the soil texture by any fortification was existent. In contrast to that, fortified roads showed higher percentages of sand, even leading to a classification as another soil type, loamy sand.
Compared to these results, fortified and unfortified roads showed almost identical grain-sizedistributions of the eroded material. Both mean soil textures are distinctly finer than their reference specimen. This result is related to the constant rainfall characteristics, which limited the maximum grain size erodible by splash and flowing water. Therefore, both particle-size distributions approximate each other. In this connection, the mean results provide a strong evidence for constant rainfall characteristics during all simulated rainfall events. Beyond that, a selective erosion of fine particles is evident. Additionally, Table 10 shows the percentage of gravel (particles >2 mm). Mixed samples of fortified roads were intentionally taken exclusive of base course material to examine only the topmost layer that was affected by erosional processes. Again, the limited kinetic rainfall energy caused a distinct difference between eroded soil and the mixed samples of the plot surfaces. While the roads featured almost 50% gravel, only 0.5% and 1.1% of gravel, respectively, was eroded. As explained with the aid of Figure 9, this shows the selective soil erosion observed during the simulation runs.
Still, the results show another allusion to selective erosion considering the results of the mixed samples. The increased percentage of gravel of unfortified forest roads compared to skid trails und forest topsoil most likely indicates an already lacking part of finer particles that was already eroded Additionally, Table 10 shows the percentage of gravel (particles >2 mm). Mixed samples of fortified roads were intentionally taken exclusive of base course material to examine only the topmost layer that was affected by erosional processes. Again, the limited kinetic rainfall energy caused a distinct difference between eroded soil and the mixed samples of the plot surfaces. While the roads featured almost 50% gravel, only 0.5% and 1.1% of gravel, respectively, was eroded. As explained with the aid of Figure 9, this shows the selective soil erosion observed during the simulation runs.
Still, the results show another allusion to selective erosion considering the results of the mixed samples. The increased percentage of gravel of unfortified forest roads compared to skid trails und forest topsoil most likely indicates an already lacking part of finer particles that was already eroded when the sample was taken. Soil erosion processes were also investigated on a statistical basis. With the help of Boxplot-Diagrams of the mean erosion sums for each surface category, the distinct difference between mostly bare road surfaces and predominantly vegetated skid trails and undisturbed sites is clearly visible (Figure 11).  Soil erosion processes were also investigated on a statistical basis. With the help of Boxplot-Diagrams of the mean erosion sums for each surface category, the distinct difference between mostly bare road surfaces and predominantly vegetated skid trails and undisturbed sites is clearly visible ( Figure 11). The Mann-Whitney test (Table 11) shows significant similarities between fortified and unfortified forest roads as well as between unrutted skid trails and undisturbed sites. Rutted skid trails are only similar to fortified roads, indicating a transitional stage between heavily disturbed road surfaces and mostly undisturbed sites. Again, Spearman-Rho correlations were calculated to investigate, if there were any significant correlations between plot characteristics and soil erosion processes. Only vegetation cover, plot sloping, soil bulk density and gravel content were examined, as soil moisture and WDPT class were not believed to have a direct influence. In exchange, the mean runoff sum was included into the correlation matrix (Table 12). The Mann-Whitney test (Table 11) shows significant similarities between fortified and unfortified forest roads as well as between unrutted skid trails and undisturbed sites. Rutted skid trails are only similar to fortified roads, indicating a transitional stage between heavily disturbed road surfaces and mostly undisturbed sites. Again, Spearman-Rho correlations were calculated to investigate, if there were any significant correlations between plot characteristics and soil erosion processes. Only vegetation cover, plot sloping, soil bulk density and gravel content were examined, as soil moisture and WDPT class were not believed to have a direct influence. In exchange, the mean runoff sum was included into the correlation matrix (Table 12). Throughout the surface categories, there is a strong correlation between mean runoff volume and eroded sediment sum. There is only one exception, unfortified road surfaces. Here, lacking fortification and bare surfaces condition high erosion rates, irrespective of the amount of superficially flowing water. There is also a strong correlation between vegetation cover and sediment yield that tends to be less distinctive with decreasing anthropogenic disturbance. On undisturbed sites, there is still a negative correlation (´0.101), but it is only a fraction of the correlation observed on fortified road surfaces (´0.949).
Studies investigating soil erosion on forest roads and skid trails concordantly state, that the vegetation cover is one of the most important control factors for soil erosion processes. There are numerous studies showing that erosion rates increase significantly with the absence of ground covering vegetation or an artificially applied mulch cover [30,35,40,[64][65][66]71,72]. These findings were also confirmed during the presented rainfall simulations. As described, the only exception to this rule is mainly uncompacted sites, were the high infiltration rates impede runoff generation. Thus, no eroding or transporting agent was present, limiting erosion processes only to splash effects.
Regarding the quantities of eroded sediment, [5] again gives a reference for comparable surrounding conditions. Here, a median soil erosion of 41.1 g¨m´2 on skid trails and 0.3-0.4 g¨m´2 on undisturbed sites is stated, indicating up to 10 times higher erosion rates on skid trails. Converted, the findings presented in this study, show a median of 20.02 g¨m´2 on rutted skid trails and 4.64 g¨m´2 on undisturbed sites. Given the fact, that the runoff volumes on skid trails presented in [5] were higher than in this study, it seems plausible that the erosion rates are also higher. Nevertheless, there is still a comparable pattern when comparing skid trail surfaces and undisturbed sites, as both studies show a clear increase of soil erosion rates on skid trails. Forest roads were not included in [5]; concerning this study at hand, they featured a median of 96.05 g¨m´2 (fortified) and 263.51 g¨m´2 (unfortified), showing distinctly higher values than skid trails presented in [5].

Conclusions
As there are only very few studies on this subject in Europe, the objective of this study was to investigate the effects of roads and skid trails on runoff generation, water retention and soil erosion processes in forested catchments. Regardless of the type of construction, persistent forest roads caused noticeably increased runoff volumes and soil erosion rates.
Compared to undisturbed sites, where the mean runoff coefficient during 43.6 l of rainfall applied to the plot surface (I = 45 mm¨h´1) was 3.7%, forest roads showed mean runoff coefficients up to 92.1%. While undisturbed forest soil showed only negligible erosion rates (4.7 g¨m´2 of eroded soil after 90 min of simulated rainfall), distinctly higher rates were measured on forest roads (118.5-272.2 g¨m´2).
Skid trails on the other hand also showed increased runoff coefficients, but mainly in the compacted, rutted parts. Here the mean runoff coefficient was 59.9%, while it was only 14.8% in the unrutted parts. Despite the increased runoff coefficients, there was no clearly detectable rise of soil erosion rates, as rutted skid trails only showed a mean sum of eroded sediment of 21.4 g¨m´2 and unrutted skid trails a mean sum of 5.1 g¨m´2 after 90 min of artificial rainfall.
Considering the dependency between plot characteristics, runoff generation and soil erosion processes, soil compaction most likely triggers overland flow processes. This could be observed especially on skid trails. As usual in soil erosion processes, a vegetation cover seems to have the most important influence as it prevents particles from being detached by splash erosion and retards the generation of overland flow by interception.
In summary, the obtained results indicate that not only forest roads but also skid trails have a massive influence on hydrological properties. This finding leads to potentially severe consequences: while forest road networks are clearly identifiable, skid trails tend to become less clearly visible when they are not used over a longer period. Nevertheless, soil compaction and therefore altered hydrological properties still remain effective. Consequently, there might be an underestimated percentage of forest soils on which significantly altered runoff processes take place when extreme rainfall events occur. Runoff generated on forest roads can be channeled quite easily by lateral ditches, which allow a certain regulation of flowing water. On skid trails, in contrast, such measures are very rarely realized.
Another aspect of this set of problems is the high density of forest road networks. Often, road network densities greater than 100 m¨ha´1 can be observed. With an increasing network density, water retention in forested catchments is more and more endangered. Therefore it has to be ensured that road networks are planned in a preferably low density and that soil compaction due to traffic is minimized.
Concerning the findings of the soil erosion assessment, high erosion rates only occurred on bare surfaces and especially on surfaces that were not sufficiently fortified. There are two possibilities to minimize soil erosion on roads. First and foremost, a vegetation cover has to be established, for example on the median strip, to reduce the erosional effects of splash and flowing water. Secondly, if there is a need for the construction of a new forest road, a fortification has to be considered, as it eventually may cause less environmental impact in the long run.