Penetration Resistance of Laacher SeE-tephra Andosols—Evaluating Rooting Conditions of Undisturbed and Excavated Forest Soils in SW-Germany

: This study discusses penetration resistance (PR) of forested Pumice-Andosol sites. PR, a key soil property inﬂuencing root growth and elongation, exerts a substantial inﬂuence on ecological site quality and tree growth. Andosols were expected to show low PR because of their unique characteristics (low bulk density, loose soil matrix). Five sites, two undisturbed and three backﬁlled, were sampled. The latter result from pumice excavation and were examined to quantify potential PR alterations in the aftermath of backﬁlling and pumice removal. Penetrologger sampling on undisturbed sites showed mean PR not exceeding 3 MPa, a literature-based, critical threshold restricting root growth, in the upper 0.80 m, indicating conditions fostering tree rooting. Backﬁlled sites mostly exhibited increased ( > 3 MPa) PR, leading to rooting restrictions even beginning at − 0.21 m. Deviations from undisturbed soils range from − 15.6 to + 109.3% depending on depth and age of the backﬁlled site. Furthermore, GIS-based data interpolation helped to identify spatial PR patterns and allowed a direct comparison before / after backﬁlling at one site. Statistical analysis revealed signiﬁcantly altered PR after backﬁlling, while a concluding ANOVA provided at least signiﬁcant governing factors (depth, area, clay + silt content, soil organic matter), albeit with only small e ﬀ ect sizes.

Pumice excavation takes place, inter alia, at the NW-exposed slope of the investigation area. As described in Section 1, local concrete factories only excavate small lots at once, leading to a progressing excavation transversally following the pumice depositions across the slope. That is why excavated areas of different ages lie side by side along with areas that have not been excavated yet. The latter function as reference areas for undisturbed conditions, the former, however, provide the opportunity to evaluate different post-excavation states regarding temporal distance from excavation. While younger backfilled soils may have not consolidated yet, leading to presumably lower PR, older sites may feature higher PR due to soil matrix consolidation or eventually loosening of the soil matrix caused by reestablished vegetation and rooting, therefore again leading to presumably lower PR in the topsoil.
Overall, five different sub-areas were observed in this study ( Figure 2). Overall, five different sub-areas were observed in this study ( Figure 2).
Pumice excavation takes place, inter alia, at the NW-exposed slope of the investigation area. As described in Section 1, local concrete factories only excavate small lots at once, leading to a progressing excavation transversally following the pumice depositions across the slope. That is why excavated areas of different ages lie side by side along with areas that have not been excavated yet. The latter function as reference areas for undisturbed conditions, the former, however, provide the opportunity to evaluate different post-excavation states regarding temporal distance from excavation. While younger backfilled soils may have not consolidated yet, leading to presumably lower PR, older sites may feature higher PR due to soil matrix consolidation or eventually loosening of the soil matrix caused by reestablished vegetation and rooting, therefore again leading to presumably lower PR in the topsoil.
Overall, five different sub-areas were observed in this study ( Figure 2).  FOR represents the totally undisturbed state with a forest stand consisting of 70-80 a old Fagus sylvatica. CUT represents an undisturbed state with no excavation but already cut-down forest (Fagus sylvatica) in preparation for excavation in winter 2018/2019. CUT EX is the same area,  Figure 2 shows the state as in summer 2019 with an already backfilled area CUT EX. The transitional stage between CUT and CUT EX, showing excavation works, is depicted in Figure 3. FOR represents the totally undisturbed state with a forest stand consisting of 70-80 a old Fagus sylvatica. CUT represents an undisturbed state with no excavation but already cut-down forest (Fagus sylvatica) in preparation for excavation in winter 2018/2019. CUT EX is the same area, but observed in summer 2019, four months after excavation and backfilling, planted with solitary Quercus petrea and Fagus sylvatica seedlings. SUC was excavated in 2015 and already shows extensive succession vegetation in addition to Fagus sylvatica and Quercus petrea seedlings planted within. REF was excavated and reforested with dense Betula pendula stands and solitary Alnus glutinosa and Quercus rubra in 2002. The aerial photography in Figure 2 shows the state as in summer 2019 with an already backfilled area CUT EX. The transitional stage between CUT and CUT EX, showing excavation works, is depicted in Figure 3.

Sampling
Soil data was sampled in the fall of 2018 prior to this study for FOR, CUT, SUC and REF, using a variety of parameters in a soil erosion study [48]. This dataset consists of 485 individual soil samples, taken at 117 sampling points (3 catenae per area, each consisting of 10 sampling points; 3 sampling points were excluded because of corrupted data). On undisturbed sites (FOR and CUT), a horizon-wise sampling took place, while five depth-slices of 0.2 m were sampled on disturbed/backfilled sites, covering the upper first meter of the soil, as they did not have pedogenetic horizons anymore. This dataset is expanded by soil moisture measurements, which were conducted simultaneously but were not included in [48]. On CUT EX, an additional 205 samples were taken in early summer of 2019, consisting of 51 individual sampling points and again five depth classes covering 0.2 meters each. These samples were also not included in [48]. Parameters used and discussed in this study are soil texture, soil moisture (SM) and soil organic matter (SOM). In order to

Sampling
Soil data was sampled in the fall of 2018 prior to this study for FOR, CUT, SUC and REF, using a variety of parameters in a soil erosion study [48]. This dataset consists of 485 individual soil samples, taken at 117 sampling points (3 catenae per area, each consisting of 10 sampling points; 3 sampling points were excluded because of corrupted data). On undisturbed sites (FOR and CUT), a horizon-wise sampling took place, while five depth-slices of 0.2 m were sampled on disturbed/backfilled sites, covering the upper first meter of the soil, as they did not have pedogenetic horizons anymore. This dataset is expanded by soil moisture measurements, which were conducted simultaneously but were not included in [48]. On CUT EX, an additional 205 samples were taken in early summer of 2019, consisting of 51 individual sampling points and again five depth classes covering 0.2 m each. These samples were also not included in [48]. Parameters used and discussed in this study are soil texture, soil moisture (SM) and soil organic matter (SOM). In order to achieve the best possible comparability between the 2018 and 2019 datasets, similar boundary conditions were sought. While soil texture and soil organic matter were defined as practically constant, soil moisture was treated as a potentially volatile parameter. Countering this possible dataset inconsistency, climate data from weather station Grenzau (situated 2 km North of the investigation area and operated by the DLR Rhineland-Palatinate) was used to calculate an antecedent precipitation index (API) according to [49], incorporating precipitation input 21 days prior to sampling and accounting seasonality using an empirical adjustment factor. API was 0.22 mm in 2018 and 0.50 mm in 2019, showing comparable and low antecedent precipitation.
Penetration resistance was measured using an Eijkelkamp Penetrologger [28]. This device samples PR in a vertical resolution of 0.01 m from 0-−0.80 m, while simultaneously locating sampling positions via built-in GPS. GPS-inaccuracies were diminished by conducting a parallel measurement with a Garmin eTrex handheld GPS using the built in point averaging method. On FOR, CUT, SUC and REF, PR was measured at 150 individual sampling points, and on CUT EX, 153 points were sampled ( Figure 4).
Forests 2016, 7, x 5 of 5 achieve the best possible comparability between the 2018 and 2019 datasets, similar boundary conditions were sought. While soil texture and soil organic matter were defined as practically constant, soil moisture was treated as a potentially volatile parameter. Countering this possible dataset inconsistency, climate data from weather station Grenzau (situated 2 km North of the investigation area and operated by the DLR Rhineland-Palatinate) was used to calculate an antecedent precipitation index (API) according to [49], incorporating precipitation input 21 days prior to sampling and accounting seasonality using an empirical adjustment factor. API was 0.22 mm in 2018 and 0.50 mm in 2019, showing comparable and low antecedent precipitation. Penetration resistance was measured using an Eijkelkamp Penetrologger [28]. This device samples PR in a vertical resolution of 0.01 m from 0-−0.80 m, while simultaneously locating sampling positions via built-in GPS. GPS-inaccuracies were diminished by conducting a parallel measurement with a Garmin eTrex handheld GPS using the built in point averaging method. On FOR, CUT, SUC and REF, PR was measured at 150 individual sampling points, and on CUT EX, 153 points were sampled ( Figure 4).

Data Analysis and Further Statistical Analysis
Soil characteristics and PR datasets are discussed on the basis of different statistical approaches. As described, pre-collected data used in a soil erosion study in the same investigation area was combined with completely new datasets of subarea CUT EX and previously unpublished data of FOR, CUT, SUC and REF. In particular, soil texture, explicitly the sum of clay and silt (C+Si), soil moisture (SM) and soil organic matter (SOM) are discussed, as they exert an influence on PR and provide a good insight into soil structure alterations made by excavation and backfilling.  Figure 4 shows all individual penetrologger sampling points, exact soil sampling positions can be derived from [48].

Data Analysis and Further Statistical Analysis
Soil characteristics and PR datasets are discussed on the basis of different statistical approaches. As described, pre-collected data used in a soil erosion study in the same investigation area was combined with completely new datasets of subarea CUT EX and previously unpublished data of FOR, CUT, SUC and REF. In particular, soil texture, explicitly the sum of clay and silt (C + Si), soil moisture (SM) and soil organic matter (SOM) are discussed, as they exert an influence on PR and provide a good insight into soil structure alterations made by excavation and backfilling.
PR is the main focus of this study, which is why this dataset is discussed in-depth, beginning with a descriptive analysis against the background of PR threshold values concerning root growth.
Following up, a spatial calculation and presentation is carried out, using sampling points as input parameters for a Kriging-interpolation. Calculations were made with Golden Software Surfer (version 8), and ESRI ArcGIS 10.1 was used for visualization. Root mean square error (RMSE) was used to compare interpolated and in-situ values for every sampling point, allowing an evaluation of interpolated data quality. All spatial interpolations were carried out for every depth-slice (0-−0.20, −0.21-−0.40, −0.41-−0.60 and −0.61-−0.80 m) and with a spatial resolution of 1 × 1 m. This data analysis should provide information, if there is any sort of spatial gradient or distribution of distinctively high and/or low PR values. Additionally, area CUT EX was interpolated separately. With the help of this dataset, a direct comparison between CUT and CUT EX was possible, viewing at both absolute and relative PR changes in the aftermath of excavation and backfilling.
In addition, a deeper statistical analysis was performed, comparing sub-areas in view of their specific PR. As a first step, the Kolmogorov-Smirnov-test was conducted to check if the datasets are normally distributed. Following this, the non-parametrical Kruskal-Wallis test was used both simplified and pairwise to evaluate if there are statistical similarities between the sites. In a further step, both undisturbed sites (FOR and CUT) were aggregated to one simplified class ("Undisturbed"), allowing a better and clearer analysis of excavation-induced changes, explicitly viewing at (a) relative deviations from the undisturbed state and (b) again testing similarities using both the simplified and pairwise Kruskal-Wallis test.
Further on, a correlation analysis was carried out, testing if PR is associated with governing factors cited in existent studies. For this purpose, measuring points of PR were associated with interpolated datasets of (clay + silt), SM and SOM, as the latter provided spatially corresponding point data for each PR sampling point. As the existent literature assumes non-linear correlations between governing factors and PR [20,27,28], the Spearman correlation coefficient was calculated. Categorization of correlation strength was based on a five step rule of thumb as proposed in [50,51] (Table 1). Table 1. Strength of correlation coefficients according to [50,51].

Correlation Coefficient
Strength of Correlation  Table 2. All statistical analyses were performed using IBM SPSS Version 25, statistical plots were made with QtiPlot version 0.9.8.3. All featured box-plot diagrams show Q 10 and Q 90 as whisker, Q 25 and Q 75 as box boundaries and Q 50 as a straight line within the box. Minimum, maximum and arithmetic mean appear as scatter points.

Soil Characteristics: Descriptive Analysis
Depth-aggregated soil texture data (Table 3) clearly indicates disturbance after excavation. Undisturbed sites (FOR, CUT) mostly show a silty loam texture in the upper 0.2 m, followed by more clayish textures. They indicate the B-Horizon, consisting of weathered tephra between −0.21 and −0.60 m. From −0.61 m on, sandier textures were found, representing slightly weathered and unweathered tephra, addressed as (weathered) C-Horizon. In contrast, disturbed and backfilled sites mostly exhibit a noticeable majority (up to 86%) of sandy loam textures throughout all depths, indicating a thoroughly mixed soil due to backfilling, which rather represents the mean texture of the formerly separated horizons.
Besides analyzing soil texture classes, a graphical analysis was conducted to illustrate mean percentages of clay, silt and sand for each site and depth category ( Figure 5). Besides analyzing soil texture classes, a graphical analysis was conducted to illustrate mean percentages of clay, silt and sand for each site and depth category ( Figure 5).     Table 3, backfilled sites exhibit a more or less homogeneous distribution with a dominant part of sand in all depths, which even shows a decreasing trend from the upper to the lower parts of the soil profiles. In the case of CUT EX, this actually leads to silt being the dominant grain size. This pattern is also observable in the box-plot diagrams ( Figure 6), where undisturbed sites (FOR, CUT) show decreasing clay and silt content with increasing depth, while backfilled sites show the opposite trend ( Figure 6, Table 4).

CL
Forests 2016, 7, x 9 of 5 only slightly weathered or unweathered LST. Analogous to soil texture percentages presented in Table 3, backfilled sites exhibit a more or less homogeneous distribution with a dominant part of sand in all depths, which even shows a decreasing trend from the upper to the lower parts of the soil profiles. In the case of CUT EX, this actually leads to silt being the dominant grain size. This pattern is also observable in the box-plot diagrams (Figure 6), where undisturbed sites (FOR, CUT) show decreasing clay and silt content with increasing depth, while backfilled sites show the opposite trend ( Figure 6, Table 4).  Concerning soil moisture, there is also a noticeable difference between undisturbed and disturbed sites (Figure 7, Table 5).  Concerning soil moisture, there is also a noticeable difference between undisturbed and disturbed sites ( Figure 7, Table 5).   Antecedent precipitation input is still preserved in undisturbed topsoil layers, as FOR and CUT clearly show increased soil moisture content in the upper parts of the soil profile. CUT throughout shows higher soil moisture values, most likely caused by the already cut-down forest stand, which eventually leads to less/no rainfall interception and/or water abstraction by trees. SUC and REF both show rather homogeneous depth profiles, indicating a faster water transport from top-to bottomlayers. This pattern is very likely linked to comparable SOM content patterns ( Figure 8, Table 6), as SOM possesses good water binding properties. Both 2018 and 2019 datasets show comparable SM contents, therefore potentially exerting a similar influence on PR. Antecedent precipitation input is still preserved in undisturbed topsoil layers, as FOR and CUT clearly show increased soil moisture content in the upper parts of the soil profile. CUT throughout shows higher soil moisture values, most likely caused by the already cut-down forest stand, which eventually leads to less/no rainfall interception and/or water abstraction by trees. SUC and REF both show rather homogeneous depth profiles, indicating a faster water transport from top-to bottom-layers. This pattern is very likely linked to comparable SOM content patterns ( Figure 8, Table 6), as SOM possesses good water binding properties. Both 2018 and 2019 datasets show comparable SM contents, therefore potentially exerting a similar influence on PR.   Concerning SOM, again, a clear difference between undisturbed and disturbed sites was found. While humic A-Horizons were/are existent on undisturbed sites-indicated by remarkably high SOM contents (SOM 16.0-19.4%) in the upper 0.2 m-a noteworthy drop of SOM-percentage is observable on disturbed sites. Here, the mean content is reduced more than 10% in the upper depth-class and it remains lower in all depth classes (Table 6). This is also a strong reference to proper backfilling practices: in the present case, thoroughly mixed material was backfilled. A separation of former topsoil during excavation and backfilling would lead to a less prominent loss of SOM and would also provide the opportunity to remodel a SOM-rich topsoil layer. Apart from that, the expectable SOM pattern with a distinctly higher SOM content in topsoil layers is only existent on undisturbed areas and all three disturbed areas show more or less homogeneous SOM distributions.

Penetration Resistance: Descriptive Analysis
Mean PR datasets are depicted in Figure 9 and Table 7, using depth-aggregated datasets.  (Table 6). This is also a strong reference to proper backfilling practices: in the present case, thoroughly mixed material was backfilled. A separation of former topsoil during excavation and backfilling would lead to a less prominent loss of SOM and would also provide the opportunity to remodel a SOM-rich topsoil layer. Apart from that, the expectable SOM pattern with a distinctly higher SOM content in topsoil layers is only existent on undisturbed areas and all three disturbed areas show more or less homogeneous SOM distributions.

Penetration Resistance: Descriptive Analysis
Mean PR datasets are depicted in Figure 9 and Table 7, using depth-aggregated datasets. Undisturbed soils (FOR, CUT) show comparable depth-patterns, which is also the case for old backfilled sites (SUC, REF); therefore, a clear distinction can be made (Figure 9): while undisturbed sites exhibit a continuously and constantly rising PR with depth, older disturbed sites show a distinct and sharp PR-increase in the upper 0.4 m, with PR increasing further in the two lower depth classes. CUT EX, the most recently backfilled site, is more comparable to undisturbed sites and additionally features rather low PR with increasing depth.
Even though there are PR max-values >2 MPa, mean PR clearly stays <2 MPa in the upper 0.2 m of all soil profiles, most likely indicating no noteworthy rooting restrictions (Table 7). Yet, viewing at   Undisturbed soils (FOR, CUT) show comparable depth-patterns, which is also the case for old backfilled sites (SUC, REF); therefore, a clear distinction can be made (Figure 9): while undisturbed sites exhibit a continuously and constantly rising PR with depth, older disturbed sites show a distinct and sharp PR-increase in the upper 0.4 m, with PR increasing further in the two lower depth classes. CUT EX, the most recently backfilled site, is more comparable to undisturbed sites and additionally features rather low PR with increasing depth.

Soil Characteristics and Penetration Resistance: Spatial Distribution
Spatial interpolations of clay and silt content ( Figure 10) show a patchy pattern for all sites in the uppermost depth category with SUC exhibiting the lowest (C + Si) contents. Going deeper, a less patchy pattern is observable at FOR and CUT, indicating comparable and ubiquitous pedogenetic processes, which led to differentiated spatial patterns. For example, a band of higher (C + Si) content is visible at −0.21-−0.40 m, retracing the slope geometry. At undisturbed sites, there is no distinctive superior pattern observable but rather abrupt spatial changes, indicating thoroughly mixed, backfilled soil. Interpolation quality was found to be sufficient, as RMSE showed no noteworthy deviations (Table 8).    Spatial SM datasets ( Figure 11) clearly depict trends described in Section 3.1, with area CUT showing by far the highest soil water content. Overall, the highest SM values were found in the upper soil layer. Here, a spatial pattern is clearly visible, as highest values in both undisturbed sites are located at the foot of the slope, indicating a surface-near water movement that follows the main slope. On the backfilled sites, which have a smoother terrain, no comparable trend is visible. Totally undisturbed soil (FOR) shows a distinctive drop of SM in deeper zones, even beginning in depth category −0.21-−0.40 m, this tendency is less prominent on disturbed sites. The RMSE for all raster datasets was satisfying (Table 9).  soil layer. Here, a spatial pattern is clearly visible, as highest values in both undisturbed sites are located at the foot of the slope, indicating a surface-near water movement that follows the main slope. On the backfilled sites, which have a smoother terrain, no comparable trend is visible. Totally undisturbed soil (FOR) shows a distinctive drop of SM in deeper zones, even beginning in depth category −0.21-−0.40 m, this tendency is less prominent on disturbed sites. The RMSE for all raster datasets was satisfying (Table 9).  Analogous to soil moisture datasets, interpolated soil organic matter data generally retraces findings explained in Section 3.1, as the highest contents are located in the upper soil layer of undisturbed sites, which is an expectable pattern for forest soils. Yet, singular areas on CUT show comparably high SOM content, even exceeding 30%, and therefore indicating organic soil horizons ( Figure 12). This might be caused by accumulation of eroded topsoil, as these areas are situated at the foot of the slope and CUT was already deforested, leaving topsoil prone to water erosion. Several heavy rainfall events that occurred in late summer could have triggered erosion processes; yet, no distinct channeling was observable, therefore sheet erosion would be the most likely accountable process. Disturbed sites show constantly lower SOM content, indicating the removal of former organic and/or humic A-Horizons and thoroughly mixed, backfilled soil material. Again, RMSE shows satisfactory results (Table 10).     Table 10. Root mean square error (RMSE), interpolated raster soil organic matter. Spatial patterns of soil characteristics on CUT EX ( Figure 13) show indications for certain backfilling actions: Beginning in the second depth class, increased clay and silt content is visible at the east/northeast half of the area. Here, soil moisture and soil organic matter content is in contrast relatively decreased. Low SOM and high C + Si are most likely an indicator for backfilled former subsoil material, originating from weathered B-Horizons. The RMSE of all datasets was again satisfying (Table 11). Table 11. Root mean square error (RMSE), interpolated raster datasets of CUT EX.

Depth [m]
Clay  (Table 12). Because of the larger number and denser spatial distribution of singular input data points, the PR RMSE is higher than those of general soil characteristics (Tables 8-11).
Spatial patterns of soil characteristics on CUT EX (Figure 13) show indications for certain backfilling actions: Beginning in the second depth class, increased clay and silt content is visible at the east/northeast half of the area. Here, soil moisture and soil organic matter content is in contrast relatively decreased. Low SOM and high C+Si are most likely an indicator for backfilled former subsoil material, originating from weathered B-Horizons. The RMSE of all datasets was again satisfying (Table 11).     (Table 12). Because of the larger number and denser spatial distribution of singular input data points, the PR RMSE is higher than those of general soil characteristics (Tables 8-11).   Figure 15 shows the PR of CUT EX. As this area provided the possibility to compare before/after excavation states, absolute and relative deviation of the PR from CUT is also depicted.  Figure 15 shows the PR of CUT EX. As this area provided the possibility to compare before/after excavation states, absolute and relative deviation of the PR from CUT is also depicted.  CUT EX shows two prominent corridors of increased PR, both starting in the SE corner and continuing to the NE and NW respectively. The SE−NW corridor shows by far the highest PR, located in the very bottom soil layer. Yet, there a still parts of CUT EX that show PR < 2 MPa even at −0.61-−0.80 m. Changes in PR are distinct areas that formally featured low PR (NW corner) and now show high PR and vice versa (SE corner). Therefore, absolute changes between −4.77 and 6.49 MPa and relative changes between −79% and 906% are observable. However, looking at mean changes, a minimal reduction of absolute PR is observable, with topsoil being the only exception ( Figure 16, Table 13). in the very bottom soil layer. Yet, there a still parts of CUT EX that show PR < 2 MPa even at −0.61-−0.80 m. Changes in PR are distinct areas that formally featured low PR (NW corner) and now show high PR and vice versa (SE corner). Therefore, absolute changes between −4.77 and 6.49 MPa and relative changes between −79% and 906% are observable. However, looking at mean changes, a minimal reduction of absolute PR is observable, with topsoil being the only exception ( Figure 16, Table 13).

Statistical Comparison of Undisturbed and Disturbed Sites
Further on, a statistical comparison of undisturbed and disturbed sites was carried out. An initial Kolmogorov-Smirnov test showed mixed results, as only six datasets were normally distributed (Table 14).

Statistical Comparison of Undisturbed and Disturbed Sites
Further on, a statistical comparison of undisturbed and disturbed sites was carried out. An initial Kolmogorov-Smirnov test showed mixed results, as only six datasets were normally distributed (Table 14). As the majority of the observed data was not normally distributed, the Kruskal-Wallis test was chosen to determine statistical similarities (Table 15).  The simplified test, only comparing all sub-areas within one analysis, shows p = 0.000 for all depths, indicating significantly dissimilar (p < 0.05) datasets. In order to obtain a more detailed analysis, a pairwise Kruskal-Wallis test was carried out for every sub-area and depth (Table 16). With only one exception, namely the topsoil layer, undisturbed sites (FOR, CUT) show statistically similar results. This is rather expectable, considering the generally more heterogeneous dataset distribution in this depth. Apart from that, both REF and SUC exhibit mostly similar penetration resistances, showing a consistent similarity with increasing depth. Recently backfilled soil at CUT EX is often comparable to undisturbed soil, but there are observable deviations indicating altered PR, especially in deeper soil zones (cf. Sections 3.2 and 3.3).
Concluding, a simplified approach was used, comparing undisturbed with disturbed sites: Both FOR and CUT were aggregated to the new category "Undisturbed", as the Kruskal-Wallis test showed overall statistically similar datasets. On the basis of this established benchmark-category, deviations from an undisturbed state could be calculated ( Figure 17).
It is clear to see, that the youngest backfilled site (CUT EX)-while still showing areas with increased PR (cf. Section 3.3)-is the only site with mainly reduced PR, contrasting with the older areas SUC and REF, which show notable PR increases. The general exception is the topsoil layer, where all areas show only small deviations.
Again, simplified and pairwise Kruskal-Wallis tests were performed, testing backfilled areas and the combined, undisturbed dataset regarding their similarities. As expected, the simplified Kruskal-Wallis tests showed overall significantly differing datasets, except for the topsoil layer (Table 17). soil at CUT EX is often comparable to undisturbed soil, but there are observable deviations indicating altered PR, especially in deeper soil zones (cf. Sections 3.2, 3.3).
Concluding, a simplified approach was used, comparing undisturbed with disturbed sites: Both FOR and CUT were aggregated to the new category "Undisturbed", as the Kruskal-Wallis test showed overall statistically similar datasets. On the basis of this established benchmark-category, deviations from an undisturbed state could be calculated ( Figure 17). Again, simplified and pairwise Kruskal-Wallis tests were performed, testing backfilled areas and the combined, undisturbed dataset regarding their similarities. As expected, the simplified Kruskal-Wallis tests showed overall significantly differing datasets, except for the topsoil layer (Table 16).  Figure  17). Significantly different areas (p < 0.05) are indicated bold.

Correlation Analysis and ANOVA
Concluding correlation analysis and ANOVA showed ambiguous results. Spearman-correlations for PR and its supposedly governing factors show a consistent trend for SM, SOM and depth (Table 19). Table 19. Spearman-correlation for penetration resistance (PR), clay + silt content (C + Si), soil moisture (SM), soil organic matter (SOM). Strong (< −0.700, > 0.700) and very strong (< −0.900, > 0.900) correlations are indicated bold. Concluding, an ANOVA was carried out to determine multi-factorial dependencies between PR and its most likely governing factors. Class breaks, basing on predefined quantiles (cf. Section 2) are shown in Table 20.  Figure 18 shows relative shares of (C + Si)-, SM-and SOM-classes for each area and depth. Viewing at rather stable parameters ((C + Si) and SOM), undisturbed soils tend to show less wide dataset distributions per depth, indicating evolved pedogenetic horizons with rather constant and characteristic properties over all sampling points. On the other hand, backfilled sites show more fragmented class distributions.

PR FOR
Resulting significant factors, calculated via ANOVA, are listed in Table 21. Overall, no medium (0.25 or higher) or large (0.40 or higher) effect sizes were identified. The largest effect size was achieved for depth, representing rising PR with increasing depth. The second largest effect size was calculated for the combined factor of (C + Si) and area. The remaining three significant factors show comparable, smaller effect sizes and involve factor "Area", or "Area" in combination with Depth, (C + Si) or SOM.  Table 20. Overall, no medium (0.25 or higher) or large (0.40 or higher) effect sizes were identified.

Discussion
It is difficult to universally classify the PR results found in this study. The two main reasons-besides the fact that there are no consistent datasets for Andosols apart from this study-are: (a) differing literature values for PRs that affect root growth and elongation, and (b) the underlying consideration, what kind of vegetation is the subject of the analysis. The latter is often the reason for differing literature values for critical PR. Fundamentally, the most consistent range of critical PR lies between 1.8 and 2.0 MPa, with [53] reporting the lowest values (rooting affected > 1.4 MPa and restricted > 1.8 MPa) and other studies stating 2.0 MPa as a critical threshold value [2,13]. Viewing at the mean PR values, undisturbed LST-Andosols exhibit noticeably low PRs, especially in the upper 0.4 m, compared to other forest soils [54].
With Fagus sylvatica forests being the potential natural vegetation in the investigation area [55,56] and given the fact that Fagus sylvatica is/was the main tree species planted on FOR, CUT and CUT EX, this seems to be a promising result, as studies point out, that root systems of Fagus sylvatica show highest densities in the upper 0.2-0.5 m [57,58]. Yet, a noteworthy number of roots was reported in depths up to 1.5 m in sandy soils [58], that is why not only the topmost soil layer has to be considered when classifying the results.
Overall, Ref. [13] most likely provides the best benchmark values for root growth, as this study was conducted at forested and sandy, backfilled sites just like the ones observed in this study. Yet, one has to consider that pumice particles may provide lower PR due to their low specific particle density. Besides that, Ref. [13] broke PR down into three major classes, <2 MPa, 2-3 MPa and >3 MPa, representing low, medium and high root growth restriction. This gave the possibility to visualize the datasets on the basis of an established categorization (Figures 19 and 20).  Figure 19. Interpolated penetration resistance, root growth restriction classes according to [13]. Right column: Area CUT EX. Figure 19. Interpolated penetration resistance, root growth restriction classes according to [13]. Right column: Area CUT EX.
Using these categories, it is clear to see, that backfilled sites show PRs that most likely will affect root growth and elongation. Critical values even begin at −0.21 m and show, almost exclusively, PR > 3 MPa. The only exception is CUT EX, which was recently backfilled. Here, particularly low PR was found. These findings were also supported by the statistical comparison between undisturbed and disturbed sites, which showed the most satisfying and clear results within the whole statistical analysis.
Here, especially the pairwise Kruskal-Wallis test showed similarities between both undisturbed sites and between disturbed sites, with only CUT EX representing a disturbed site that is more comparable to undisturbed conditions. A possible explanation for this is that the backfilled soil at CUT EX has not consolidated yet, as backfilling took place only a few months before PR measurements. Useful studies concerned with self-weight soil consolidation are scarce, as the majority of consolidation analysis has studied construction sites and their soil stability or self-weight consolidation of suspended sediment in a water column (e.g., [59][60][61][62][63]). However, consolidation timespans described may hint towards an ongoing self-weight induced compaction taking place after sampling.
Another successive process eventually leading to higher PR in the future might be internal erosion, more precisely the process of suffusion by seepage leading to clogged pores as described in [64][65][66][67][68]. Within the thoroughly mixed, backfilled material, an interflow-or seepage-driven suffusion of fine particles may most likely occur, as the still present, coarse pumice particles provide wide pores allowing seepage particle transport [69]. In [64], the authors even explicitly discuss the effect of suffusion on PR of Andosols, but they describe the opposite phenomenon, where the upper soil layers lose fine particles through suffusion. However, they found clear evidence that loss of fine particles Figure 20. Relative shares of root growth restriction classes according to [13], interpolated datasets.
Using these categories, it is clear to see, that backfilled sites show PRs that most likely will affect root growth and elongation. Critical values even begin at −0.21 m and show, almost exclusively, PR > 3 MPa. The only exception is CUT EX, which was recently backfilled. Here, particularly low PR was found. These findings were also supported by the statistical comparison between undisturbed and disturbed sites, which showed the most satisfying and clear results within the whole statistical analysis.
Here, especially the pairwise Kruskal-Wallis test showed similarities between both undisturbed sites and between disturbed sites, with only CUT EX representing a disturbed site that is more comparable to undisturbed conditions. A possible explanation for this is that the backfilled soil at CUT EX has not consolidated yet, as backfilling took place only a few months before PR measurements. Useful studies concerned with self-weight soil consolidation are scarce, as the majority of consolidation analysis has studied construction sites and their soil stability or self-weight consolidation of suspended sediment in a water column (e.g., [59][60][61][62][63]). However, consolidation timespans described may hint towards an ongoing self-weight induced compaction taking place after sampling.
Another successive process eventually leading to higher PR in the future might be internal erosion, more precisely the process of suffusion by seepage leading to clogged pores as described in [64][65][66][67][68]. Within the thoroughly mixed, backfilled material, an interflow-or seepage-driven suffusion of fine particles may most likely occur, as the still present, coarse pumice particles provide wide pores allowing seepage particle transport [69]. In [64], the authors even explicitly discuss the effect of suffusion on PR of Andosols, but they describe the opposite phenomenon, where the upper soil layers lose fine particles through suffusion. However, they found clear evidence that loss of fine particles reduces PR by destabilizing the soil matrix. Therefore, it seems plausible, that enrichment of fine particles in deeper soil layers leads to the opposite effect. Subject to the condition that a successive downward transport of fine particles will occur in the aftermath of backfilling, CUT EX most likely has not fully consolidated yet and may show higher PR in the future. Ongoing penetrologger measurements of this area are taking place, and further investigations have to clarify if a further consolidation takes place that would explain PR values differing from older backfilled sites.
Trying to isolate major correlations and/or governing factors within the dataset was not as successful as hoped. Concerning correlations, at least general trends were observable. Some of them were expectable, like a positive correlation of depth and PR, others hinted towards the beneficial effects of SOM and SM, as they showed negative correlations.
ANOVA analysis revealed effect sizes merely exceeding small effects according to [52]. Even depth-a factor that should, on first sight, clearly govern PR-showed an effect size of only 0.167. However, this rather surprisingly small effect size is most likely caused by the specific properties of Andosols (undisturbed soils): High sand content and loose soil matrix in unweathered and undisturbed C-Horizons led to less PR-increase with depth, thus not generating a PR increase that would have led to a larger effect size.
Interestingly, ANOVA showed only small effect sizes (0.095 and 0.090) for the factor "area", even though it was proven that they feature statistically dissimilar datasets. Yet the effect is still significant, albeit small. Larger effect sizes were achieved for the combined factor "(C + Si) and Area", indicating that disturbance (backfilling) is still one of the main verifiable factors governing PR, as (C + Si) is certainly affected by backfilling, again indicated by homogenized soil textures.

Conclusions
As one of the first studies published, PR benchmarks for forested Andosol sites were collected. They clearly show undisturbed conditions that may foster root growth and elongation against the background of published critical threshold values. Additionally, effects of pumice excavation and soil backfilling were quantified. On backfilled sites, PR was significantly altered, most likely affecting root growth and elongation. Further studies, partly already conducted and in preparation for publication, have to address two main topics. First of all, a comprehensive analysis of root systems has to be carried out in order to draw conclusions on correlations and dependencies between PR and rooting. Secondly, repeated PR measurements on CUT EX have to be conducted to verify the hypothesis that low PR values are caused by a soil matrix that has not consolidated yet, potentially leading to higher PR in the longer aftermath of backfilling actions. Within this study, the soil physical fundament has been laid for further investigations.