Root Biomechanical Traits in a Montane Mediterranean Forest Watershed : Variations with Species Diversity and Soil Depth

Plant roots play a key role in stabilizing slopes, particularly in the Mediterranean region, characterized by rough and unstable terrain. However, forest species differ in their stabilizing capacities. The purpose of this study is to fill the gap of knowledge on root biomechanical properties of relevant Mediterranean trees and shrubs in relation to slope stability. Root specimens of typical montane Mediterranean tree and shrub species were sampled in Southern Italy. Root characteristics, such as tensile strength (Tr) and root area ratio (RAR), were assessed from live roots sampled in trenches, while root cohesion was calculated. Power law functions yielded the best fit for the relationship of Tr versus root diameter; however, no significant relationship was found between root strength and root moisture content. RAR varied amongst different tree and shrub species. Roots of Quercus cerris L. were the most resistant to breaking under tension, while roots of Ilex aquifolium L. had the highest tensile strength among all shrub species. Results provide quantitative information on the role of root systems of montane Mediterranean forest species in stabilizing soils and will improve modeling of landslide susceptibility to the prevention and mitigation of natural hazards in mountain environments.


Introduction
Slope stability and soil erosion are significant environmental and societal issues worldwide, particularly in the Mediterranean region [1,2].The importance of these geomorphic processes is growing in relation to the intensification of precipitation, as well as in response to the increasing occurrence of extreme rain events in recent years [3], which are predicted to increase even more in the future [4].In addition to rainfall, earthquakes, storm surges, shallow groundwater-level changes, and streambank erosion can all induce instability in hillslopes, triggering hydrogeomorphic hazards, such as mass movements of rock, debris, and soil [5].Moreover, factors such as poor land management, forest degradation and conversion, forest fires, road construction, urban sprawl, and climate change can all affect slope stability [6].In this context, forest conservation and appropriate forest management play significant roles in preventing and mitigating hydrogeomorphic hazards induced or exacerbated by human activities.In fact, the presence of plant roots enhances soil properties (e.g., organic matter content, soil structure, pore size distribution) and significantly increases the shear strength of soils [7][8][9][10].Analyzing landslide inventories clearly shows an increase in landslide frequency and an accelerated displacement of existing landslides after forest harvesting and vegetation removal [11,12].These studies attribute enhanced landslide occurrence to the reduction in root reinforcement following the removal of forest vegetation.Forest cover has also been suggested to ameliorate slope stability via rainfall interception [13], although other studies have shown that this influence is very limited [14].Indeed, forests improve the infiltration capacity of the soil, thus reducing the potential for surface erosion [6].
The anchorage of roots into stable substrate and the consequent improvement of slope stability is related to increased soil shear strength [7,15], which also depends on specific properties of the root systems, such as root density and distribution [16] and root tensile strength.The role of vegetation in slope stability has been clearly established [9], whereby root reinforcement mechanisms are particularly dominant in forest soils [17].In fact, on steep soil slopes, techniques such as restoring woody vegetation cover are applied to prevent shallow failures [18].Plant roots increase the strength of the soil-root mass by enhancing the confining pressure and the resistance to sliding [19].Root tensile strength together with root distribution and root area ratio (RAR) are important factors for determining soil reinforcement and contributing to the safety factor of hillslopes [20].RAR provides a measure of root density within the soil and is defined as the fraction of the soil cross-sectional area occupied by roots [21].RAR is strongly influenced by local soil and climate characteristics, land management, and associated vegetation.RAR decreases with depth and distance from the tree trunk [22].Another essential parameter that influences slope stability is root cohesion, which varies according to the species and is highly dependent on root morphology [23].Zhang et al. [24] showed that root cohesion is produced by the transfer of shear stress in the soil due to tensile resistance in the roots.
Many studies have evaluated mechanical properties of plant roots and their role in soil stability [17,[25][26][27], but there is a lack of such information in Mediterranean environments and vegetation.Furthermore, assessing root cohesion values for different tree and shrub species could improve our understanding of how vegetation affects slope stability and how best to apply phytotechnology for landscape restoration.One aspect which has been rarely taken into account in such investigations is the moisture content in the roots for different species and site characteristics, and its possible effect on mechanical properties [24].The basic assumption is derived from wood technology studies [28], which report a remarkable effect of moisture content on mechanical properties (including tensile strength) in the wood hygroscopic range from 0 to 30%, where 30% is roughly assumed to be the wood fiber saturation point (FSP) and values above this usually do not influence mechanical properties [28].The fiber saturation point represents the moisture content at which the maximum water content bonds to cell walls by means of hydrogen bonds and the assumed value of 30% only significantly changes in very degraded wood, such as waterlogged wood [29].Root biomechanical traits are greatly affected by species and site conditions [30], while they represent relevant properties to be considered for modeling slope stability and landslide initiation.This study aims to investigate the role that important montane Mediterranean forest trees and shrubs play in soil stability and erosion control in largely forested regions of Southern Italy, such as the Sila National Park.In this context, soil erosion and slope stability are critical issues due to geomorphology, ecosystem management, forest fires, and rainfall regimes with irregular but intense events.The purpose of this study was to validate the following objectives:

•
Soil-root reinforcement in mountains and hillslopes is highly dependent on forest species biodiversity, particularly in the Mediterranean region; • Root properties, such as tensile strength and root cohesion, are important biomechanical traits that could be efficiently utilized in nature-based solutions to improve soil stability in the region; • Root moisture conditions affect root tensile strength.

Study Area and Plant Material
The investigation focuses on a mountain catchment in Southern Italy within the Calabria Region.This region is characterized by more than 600,000 ha of forest, corresponding to 40% of the total land surface of the region [31].According to the Italian National Forest and the Carbon Inventory [31], the main forest types in Calabria are: Chemonstnut (Castanea sativa Mill.) (11.3%),Calabrian pine (Pinus nigra subsp.laricio (Poir.)Maire) (12.2%), and oaks (Quercus spp.) (7.6%).The study catchment is located in the Sila Greca mountain range (39 • 28 50 N, 16 • 30 12 E) in proximity to Sila National Park and is named "Bonis" after the stream that crosses through it (Figure 1). Root moisture conditions affect root tensile strength.

Study Area and Plant Material
The investigation focuses on a mountain catchment in Southern Italy within the Calabria Region.This region is characterized by more than 600,000 ha of forest, corresponding to 40% of the total land surface of the region [31].According to the Italian National Forest and the Carbon Inventory [31], the main forest types in Calabria are: Chemonstnut (Castanea sativa Mill.) (11.3%),Calabrian pine (Pinus nigra subsp.laricio (Poir.)Maire) (12.2%), and oaks (Quercus spp.) (7.6%).The study catchment is located in the Sila Greca mountain range (39°28′50" N, 16°30′12" E) in proximity to Sila National Park and is named "Bonis" after the stream that crosses through it (Figure 1).The catchment has a total area of 139 ha, most of which is forested.The area was largely reforested 50 years ago with Calabrian pine plantations, and includes small areas of mixed stands with chestnut and alder (Alnus glutinosa (L.) Gaertn) riparian forests.The climate of the Calabria region is Mediterranean, characterized by mild winters and hot summers with little precipitation.However, in mountainous and inland areas, colder winters with snowfall and cooler summers with some precipitation occur [32].Mean annual precipitation at the Bonis site is 1179 mm and mean annual air temperature is 9 °C [33].Acidic plutonic rocks [34] underlie two types of soils (Typic Xerumbrepts and Ulpic Haploxeralfs) in the catchment.The Typic Xerumbrepts are variable in thickness, from several centimeters to over a meter, with a dominant sandy texture; the presence of the granite bedrock is common in the most eroded areas.The Ultic Haploxeralfs are moderately thick with a silty clay loam texture [34].The study area is strictly uniform from a geological point of view.The soils have mainly a sandy-silty and sandy-silty clayey texture [33].The strong sandy component of the soil suggests a low degree of pedogenetic evolution in the area.The measured pH values allow the soils to be classified, regardless of the sampling depth, from ultra-acids (<4.5) to moderate acids (5.1-6.0)[33].These values are also related to the nature of the granite lithological matrix emerging in the study area [33].The catchment has a total area of 139 ha, most of which is forested.The area was largely reforested 50 years ago with Calabrian pine plantations, and includes small areas of mixed stands with chestnut and alder (Alnus glutinosa (L.) Gaertn) riparian forests.The climate of the Calabria region is Mediterranean, characterized by mild winters and hot summers with little precipitation.However, in mountainous and inland areas, colder winters with snowfall and cooler summers with some precipitation occur [32].Mean annual precipitation at the Bonis site is 1179 mm and mean annual air temperature is 9 • C [33].Acidic plutonic rocks [34] underlie two types of soils (Typic Xerumbrepts and Ulpic Haploxeralfs) in the catchment.The Typic Xerumbrepts are variable in thickness, from several centimeters to over a meter, with a dominant sandy texture; the presence of the granite bedrock is common in the most eroded areas.The Ultic Haploxeralfs are moderately thick with a silty clay loam texture [34].The study area is strictly uniform from a geological point of view.The soils have mainly a sandy-silty and sandy-silty clayey texture [33].The strong sandy component of the soil suggests a low degree of pedogenetic evolution in the area.The measured pH values allow the soils to be classified, regardless of the sampling depth, from ultra-acids (<4.5) to moderate acids (5.1-6.0)[33].These values are also related to the nature of the granite lithological matrix emerging in the study area [33].
The basin extends from an altitude of 975 to 1300 m a.s.l.(above sea level), with an average elevation of 1131 m a.s.l.The length of the main stream channel is 2.2 km with a median gradient of 43.4%.The catchment was instrumented in 1986 to investigate the hydrological cycle of the area [35]; in 2004, monitoring of carbon and water vapor exchange at the canopy level (eddy covariance technique [36]) was initiated.The catchment is now associated with the ALForLab public-private laboratory (www.alforlab.it).Roots of typical montane Mediterranean tree and shrub species, Castanea sativa Mill, Pinus nigra subsp.laricio (Poir.)Maire, Alnus glutinosa (L.) Gaertn., Quercus cerris L., Crataegus monogyna Jacq., Cytisus scoparius L., Ilex aquifolium L., and Spartium junceum L. (Table 1), were sampled to assess root properties of biomechanical interest, including tensile strength, root area ratio, and root cohesion for each species.The root samples of the study plants were collected in the Bonis catchment at different locations in summer and fall, 2017, by sampling two trees or shrubs for each species.Sample trees were selected according to the mean diameter breast height (DBH) of the population of each selected tree in the forest catchment, while sampled shrubs had a height and size representative of the shrub population of the forest (Table 2).For each woody shrub, 30 × 30 cm trenches were dug 50 cm from the stem.For sample trees, 50 cm wide × 50 cm deep trenches were dug 100 cm from the stem [25,35].Trenches were dug along the mountain slope contours [37] and root samples were collected along the trench section.To preserve the moisture content of excavated roots, they were stored in sealed plastic bags and transported to the laboratory in a refrigerated box; later, the root moisture content was measured from the root samples.The protocol was adopted because one purpose of the investigation was to determine the original water content of the roots and the related mechanical behavior.It is well known that wood and, in general, lignocellulosic materials lose water in drier conditions because of capillarity and evaporation.Within its hygroscopic field (0-30%), wood exchanges OH − hydroxyl groups with the environment in the gas phase [28,29].In the laboratory, the bags were stored in a refrigerated room at a constant temperature of 4 • C for about a week after sampling.Afterwards, 346 total samples of all species (20 cm long) were selected, consisting of about 40 samples for each species.The diameters were measured with a digital caliper (Mitutoyo 500, Lainate, Milano, Italy) (accuracy: ±0.03 mm) at three different positions along their length (at the ends and center).Average diameter was assigned to the entire root and used in data analysis.Roots were grouped into four diameter classes (0-1 mm, 1-2 mm, 2-5 mm, 5-10 mm) [36].The sampled trees were considered representative of the average conditions of the forest stand.

Root Tensile Strength
Root terminals of all samples were reinforced with synthetic resin, according to a method proposed by Nilaweera and Nutalaya [38].After resin application, roots were left for 2 days in a 4 • C cold room for resin thickening.This resin application method enables a better transfer of force from the clamping devices for tensile strength to the root without slipping.The 346 root samples were randomly divided into two subgroups of equal numbers, including each diameter class and each plant species; the first subgroup was utilized for biomechanical measurements with roots at field moisture conditions, whereas the samples of the second subgroup were saturated before biomechanical measurement.Fresh water content (FWC) is defined as the weight of water in plant tissues, expressed as a percentage of the oven-dry weight of the plant sample.Usually, moisture content of roots ranges from about 30% to more than 200% of the weight of an oven-dry plant sample [39].Root water saturation was reached after soaking the roots in water for 48 h.After 48 h, the samples reached a constant weight.For all steps, weight was carefully measured at an accuracy of 0.001 g.After saturation, the samples were oven-dried for 24 h at 103 • C and reweighed.The temperature of 103 • C is normally used in wood technological analysis because it allows water to evaporate, but it does not cause wood combustion [28,40].Root fragments used to measure moisture content were discarded and not used for further analysis.The FWC and the maximum water content (MWC) of the root samples was calculated as: where w f c is the weight of the root at field moisture condition of the sample collected from the soil, w wc is the weight of the root at water saturated state, and w 0 is the weight of the oven-dry root.
Values obtained by Equations ( 1) and ( 2) can be related to the expected behavior of roots at the fiber saturation point (FSP), which is a well-known threshold for changing mechanical performance.Each root was tested for tensile stress until failure.When breakage occurred at one of the clamping sites, the result was discarded, and only samples that broke near the middle of the root were accepted for data analysis [41].A universal testing machine ("ProLine"-Zwick Roell Company, Ulma, Germany, cell load capacity 50 kN) was used to measure root tensile strength (T r ).The machine combines three functions: (1) Traction force generation; (2) measuring load and displacement; and (3) data acquisition (Figure 2).The samples were subjected to movement at a constant test speed of 10 mm/min.Tensile strength (T r ) was calculated according to Bischetti et al. [42]: where F max is the maximum force needed to break the root and D is the mean root diameter before stretching.The samples were subjected to movement at a constant test speed of 10 mm/min.Tensile strength (Tr) was calculated according to Bischetti et al. [42]: where Fmax is the maximum force needed to break the root and D is the mean root diameter before stretching.

Root Area Ratio
Vertical profiles of 50 × 50 cm for trees and 30 × 30 cm for shrubs were used to calculate RAR [43].Within each trench, root diameter was measured using an electronic caliper; at a later stage, the number of roots in the trench was estimated via high-definition photographs [43].Before analysis, photographs were corrected to eliminate geometric distortion (Figure 3).For RAR calculation, we subdivided the excavation area into individual 10 cm depth intervals, yielding five depth class intervals for trees and three intervals for shrubs [44].

Root Area Ratio
Vertical profiles of 50 × 50 cm for trees and 30 × 30 cm for shrubs were used to calculate RAR [43].Within each trench, root diameter was measured using an electronic caliper; at a later stage, the number of roots in the trench was estimated via high-definition photographs [43].Before analysis, photographs were corrected to eliminate geometric distortion (Figure 3).For RAR calculation, we subdivided the excavation area into individual 10 cm depth intervals, yielding five depth class intervals for trees and three intervals for shrubs [44].RAR values at depth increments of 10 cm were calculated by the following equation: where, πr 2 is the total root cross-sectional area in each layer (mm 2 ) and A is the soil trench area (50 cm × 50 cm) measured in each layer (10 cm in depth and 50 cm wide).

Root Cohesion
Waldron and Dakessian [45] and Wu et al. [10] were pioneers in modeling relevant soil-root interactions and quantifying the root reinforcement effect on soil.The model developed by these authors assumes that all roots cross the shear plane during slope failure.According to Waldron [45], the primary influence of root reinforcement can be expressed as a cohesion term in the Mohr-Coulomb equation.The Mohr-Coulomb break criterion is a linear equation in main stress space, characterizing the conditions for which an isotropic material will fail [46,47].To include root reinforcement, the Coulomb equation, which quantifies the resistance to soil sliding, is extended to root-permeated soils (Equation ( 5)) by considering an increased shear strength as a function also of root cohesion [7]: where S sr is the soil shear strength against slippage (including root reinforcement), c s is the effective cohesion of the soil, c r is the apparent cohesion provided by roots, σ is the normal stress due to the weight of the soil within the sliding mass, u is the soil pore water pressure, and ϕ is the effective internal friction angle of the soil.In this model, c r (kPa) depends on root tensile strength (T r (kPa)) and on the cross-sectional area of roots at the shear plane (RAR): where i indicates the diameter class and N is the number of classes, k represents the random orientations of the roots with respect to the failure plane and depends on the friction angle of the soil and the distortion angle of sheared roots, and k is a factor for quantifying non-simultaneous breaking of roots.For a friction angle ϕ = 30 • , k varies between 1.03 and 1.13, whereas 0.88 < k < 1.02 is used for ϕ = 20 • and k = 1.2 for ϕ = 40.The values of k' should be selected according to site conditions, considering a soil friction angle of 30 • .Several k values have been defined [45,48,49], and for forest areas, this value was set to 0.56, according to Hammond et al. [49].
The generally accepted form of the relationship between T r and d is a simple power function [7,9]: where α and β are empirical constants.

Statistical Analysis
The Kruskal-Wallis nonparametric multiple comparisons test was used to test for significant differences in tensile strength of FWC and MWC roots of each species and for each of the four root classes.The relationship between root tensile strength (T r ) and root diameter (D) was defined by a power-law regression.Goodness-of-fit models used to assess these two parameters were R 2 and p-values, utilizing the log-transformed values of T r and D; Spearman's rank-order correlation (r s ) between predicted and observed values was also utilized to test statistical significance.The variation of tensile strength values among trees and shrubs was evaluated by ANCOVA analysis, accounting for the diameter as a covariate factor in the linear regression (logT r -logD) (Supplementary Materials Figure S1).Least significant difference (LSD) post hoc tests were used to check differences among samples.This ANCOVA procedure was also used to compare the RAR values among tree and shrub species and finally in trees and shrubs in the first three depth classes, considering depth classes as a covariate.In this context, an LSD post hoc test was performed to determine which species in the sample differ.Lastly, differences in RAR values between species in soil layer classes were verified by a two-way ANOVA test.Statistical analyses were conducted using the R software program (www.r-project.org,R version 3.5.3,University of Auckland, Auckland, New Zealand).
ANCOVA analysis indicated significant differences in root tensile strength among species (F (7,730) = 33.96,p-value < 0.001) (root diameter as covariate).Indeed, analysing T r (D) log-transformed values, we found that the angular coefficients were different among the species, resulting in non-parallel curves (Figure 4).
The ANCOVA analysis confirms that the tensile strength depends on diameter and that the slope of the regression between Tr and D is species-specific.The LSD test divided the tested species into four significantly different groups based on tensile strength: 1-Quercus cerris L.; 2-Ilex aquifolium L. and Cytisus scoparius L.; 3-Spartium junceum L.; 4-Crataegus monogyna Jacq., Castanea sativa Mill., Pinus nigra subsp.laricio (Poir.)Maire, and Alnus glutinosa (L.) Gaertn.Quercus cerris L. had the largest values of Tr and therefore greater root resistance.Moreover, considering plant type (trees versus shrubs) the two series are parallel (ANCOVA, F (1,342) = 0.11, p-value = 0.73) and significantly different (ANCOVA, F (1,342) = 2.72, p-value < 0.1).The LSD post hoc test highlighted that shrubs have a significantly higher root tensile strength than trees and, therefore, the shrubs have a higher resistance to tension stress.The LSD post hoc test highlighted that shrubs have a significantly higher root tensile strength than trees and, therefore, the shrubs have a higher resistance to tension stress.
Young modulus could add interesting information to predicting root reinforcement [50].Nevertheless, the stress strain curves of many samples show a shape quite tortuous, as in agreement Young modulus could add interesting information to predicting root reinforcement [50].Nevertheless, the stress strain curves of many samples show a shape quite tortuous, as in agreement with Commandeur and Pyles [50] too.Indeed, in order to measure the Young modulus, several traits of stress strain curve must be selected, applying for each trait a different model to determine the root's real behavior.This part of the research needs a focused investigation.The results of the modulus of elasticity, processed by testing machine "ProLine"-Zwick Roell" for each species, were: Castanea sativa Mill.
For tree species, results of the ANCOVA showed that RAR is significantly different among species F (4,35) = 4.26, p-value < 0.0065 (depth as covariate).The post hoc LSD test showed that Pinus nigra subsp.laricio (Poir.)Maire and Quercus cerris L. had significantly higher RAR values than the other tree species.This result is confirmed by the ANOVA test, in fact, the main effect of species was significant, F (3,20) = 17.31, p value < 0.0001, as was the main effect of layer depth classes, F (4,20) = 14.06, p-value < 0.0001.Furthermore, the interaction of these two factors was significant, F (12,20) = 2.59, p-value = 0.029.The post hoc LSD test showed that the RAR values in the 10-20 cm, 20-30 cm, and 30-40 cm depth classes were significantly different with respect to the other classes (0-10 cm and 40-50 cm).For shrubs, results of the ANCOVA confirmed that the RAR is not significantly different among species F( 4,19 ) = 0.075, p-value = 0.99 (depth as covariate).This result is confirmed by the ANOVA test, in which the main effect of species was not significant, F( 3,12 ) = 0.09, p value = 0.96, nor was the main effect of layer depth classes, F( 2,12 ) = 3.08, p-value = 0.08.The interaction of these two factors was also not significant, F( 6,12 ) = 0.11, p-value = 0.99.The ANCOVA for tree and shrub categories in the first three depth classes showed that the RAR is not significantly different among categories (F (2,45) = 2.34, p-value = 0.10).

Root Cohesion
Equation ( 5) for estimating soil reinforcement by roots was applied and the results were highly variable for both species and depth.Values of c r ranged from 1.90 kPa for Alnus glutinosa (L.) Gaertn to 77.53 kPa for Quercus cerris L., while for shrubs, the highest c r was for Ilex aquifolium L. (35.80 kPa) and the lowest for Spartium junceum L. (13.86 kPa) (Table 7 and Figure 6).Root cohesion in the 10-20 cm depth was highest in all species.For trees, the second and third depth layer had elevated c r values.Below this depth, c r declined similarly to other studies [9,30,41,43,52,53] Table 7. Mean and standard deviation values of root cohesion for each species.effect of layer depth classes, F(2,12) = 3.08, p-value = 0.08.The interaction of these two factors was also not significant, F(6,12) = 0.11, p-value = 0.99.The ANCOVA for tree and shrub categories in the first three depth classes showed that the RAR is not significantly different among categories (F(2,45) = 2.34, p-value = 0.10).

Root Cohesion
Equation ( 5) for estimating soil reinforcement by roots was applied and the results were highly variable for both species and depth.Values of cr ranged from 1.90 kPa for Alnus glutinosa (L.) Gaertn to 77.53 kPa for Quercus cerris L., while for shrubs, the highest cr was for Ilex aquifolium L. (35.80 kPa) and the lowest for Spartium junceum L. (13.86 kPa) (Table 7 and Figure 6).Root cohesion in the 10-20 cm depth was highest in all species.For trees, the second and third depth layer had elevated cr values.Below this depth, cr declined similarly to other studies [9,30,41,43,52,53] Table 7. Mean and standard deviation values of root cohesion for each species.

Discussion
The analysis of the root systems of important montane Mediterranean tree and shrub species indicates significant differences in their strength characteristics.A total of 346 undamaged roots with diameters smaller than 1 cm and a constant root length of 20 cm were tested for tensile strength.In

Discussion
The analysis of the root systems of important montane Mediterranean tree and shrub species indicates significant differences in their strength characteristics.A total of 346 undamaged roots with diameters smaller than 1 cm and a constant root length of 20 cm were tested for tensile strength.In general, the root tensile strength (T r ) decreased with increasing root diameter (D).The generally accepted form of the relationship between T r and D is a simple power function [7,9].Many studies have confirmed the validity of this relationship [41,45], while the influence of root diameter on tensile properties has been studied from the perspective of chemical components [23].ANCOVA tests showed that species are the most important factor determining root tensile strength, as also shown by previous studies [41].Furthermore, root diameter acts together with species characteristics to affect root tensile strength.Differences in tensile strengths between trees and shrubs were not statistically significant.Tensile strength values of Quercus cerris L. roots were the highest among the species tested; Alnus glutinosa (L.) Gaertn and Pinus nigra subsp.laricio (Poir.)Maire had the lowest tensile strength values, while other species showed an intermediate behavior.The tensile strength values of Pinus nigra subsp.laricio (Poir.)Maire are comparable with the other Pinus species [24,54].Generally, roots of deciduous trees have higher tensile strengths than conifers [54] and this is confirmed in our study, except for Alnus glutinosa (L.) Gaertn, which could be explained by its wood density being similar to the density of pine trees.The values of tensile strength obtained in all the hardwoods and shrubs fall within ranges reported by other studies [24,27,41,55], especially with regard to Quercus cerris L. [26].However, in many studies, most of the species showed large variability in tensile strength [56].This variability can arise from different conditions in the growing environment [26], including soil moisture content and physical characteristics, and other factors, such as the methods used for preserving roots prior to measurements, the time elapsed before testing, and differences in machines used to conduct tensile strength tests [41].
The comparison between tensile strength measured in fresh roots (FWC), as collected in the field, and water-saturated roots (MWC) requires more discussion.Apparently, by choosing the wettest conditions, the test should represent the minimum strength value of a particular root and, more likely, the strength of the root in conditions of root failure during landslide initiation caused by soil water saturation [54].The role of moisture content in roots is controversial because an increase in moisture has been shown to induce a decrease in mechanical strength [54,57].Yang et al. [54] showed that the quadratic equation was the most applicable model for explaining the relationship between root moisture content and tensile strength.A comparable result was obtained by [58], who also took into account decreasing root diameter due to water content loss [58,59] and, consequently, root strength.However, in our study, the variation of water content in the roots, from ambient soil moisture conditions to maximum water content, did not significantly affect tensile strength performance.This can be related to the fact that moisture content values were generally higher than the wood fiber saturation point (about 30%), above which the physical-mechanical properties usually did not change significantly [56].In any case, in relation to the high variability of soil moisture content, soil properties, climate conditions, root structure (especially diameter [54]), species (hardwood or conifers), and tree physiological status [57] (vegetative season or dormant period), fresh root water content needs to be further investigated, as it would be preferable to make experiments on the real fiber saturation point at different conditions.
Our findings confirm the trend of decreasing RAR with soil depth that has been associated with reduced availability of nutrients, decrease in soil aeration, and the presence of more compact layers [60].RAR patterns were investigated separately for trees and shrubs and assessed to depths of 0.50 m and 0.30 m, respectively.Our method only considered roots typically smaller than those of second-and third-order roots.Furthermore, the effects of different root system structures for various species were not considered.For example, we found that Pinus nigra subsp.laricio (Poir.)Maire, Quercus cerris L., and Alnus glutinosa (L.) Gaertn have a similar RAR distribution trend with very deep taproots and well-developed radical structures, whereas Castanea sativa Mill.had the highest RAR value in the second and fourth depth class layer and roots expanded more laterally than in depth [61].For woody shrubs, RAR values were very similar for all four species, while the largest number of roots were present in the 10 to 20 cm soil depth interval.Significantly, because of our mountainous site, more lateral roots were located below the upper 10 cm soil depth; plants growing on relatively flat sites typically have more abundant lateral roots in the upper 10 cm of soil compared to hillslopes [62].RAR trends for shrub and tree profiles in our study can be compared with previous results [63].Lombardi et al. [62] found that the number of fine roots of Spartium junceum L. was highest in the upper 10 cm and decreased with soil depth.In our study, RAR values were generally higher compared to the values obtained in other Mediterranean and alpine sites [27,41], except for values of chestnut, which are comparable with the results obtained by Bischetti et al. [41].
The model reported by Wu et al. [10] was used to calculate root cohesion.The choice of the coefficients k and k , as reported in Equations ( 5) and ( 6), play a crucial role in estimating root cohesion.Many authors have demonstrated that k can assume values much smaller than 1.0, dramatically affecting the value of root reinforcement [45,48].Hammond et al. [49] proposed a k factor of 0.56 for forest vegetation.Moreover, the value of root cohesion varies as a function of the RAR distribution patterns, although some discrepancies can be associated with the dependency on root tensile strength and root diameter distribution.Usually, estimated root cohesion values are significantly higher in the upper soil layers (10-40 cm) and decline with soil depth.Furthermore, our results show that the investigated plant species have comparable reinforcement to previously studied tree and shrub species [22,60].Maximum root cohesion occurred from 10 to 20 cm depth in all species examined in our study.The type of vegetation affects the value of root cohesion as well as the land use.Schmidt et al. [30] noted lower root cohesion values in intensively managed clear-cut forests compared to natural forests.The root cohesion values obtained in the Bonis basin were very similar to those found by Schmidt et al. [30] for unharvested natural forests and those reported by Greenway [22] (average root cohesion up to 40 kPa for trees).

Conclusions
Our findings contribute to knowledge about the root tensile strength, RAR, and the cohesion model of some typical trees and wood shrubs in the Mediterranean region.The results could significantly increase knowledge to improve natural hazard prevention and mitigation in Mediterranean forest environments.The results related to tensile strength, RAR, and cohesion are not all comparable with values obtained in previous studies, indicating the significant role of site features, such as soil type and management history of forest stands.While Quercus cerris L. maintained the strongest role in soil reinforcement, pioneer shrub species like Spartium junceum L., with root tensile strengths comparable or even higher than some tree species, have important stabilizing roles.Shrubs can increase soil shear strength considerably without having the more complex effects attributed to wind stress or static surcharge of the trees.Furthermore, retaining woody shrubs with high root strength via forest cable logging operations could benefit soil stability, compared to harvesting operations that damage the understory layer [64].Increases in soil shear strength are mostly attributed to roots (1 < D < 10 mm) that often exert their maximum tensile strength during soil displacement [43].Our results confirmed the very important role of small roots, which increase soil shear strength by anchoring the soil mantle and by forming a binding network within the soil profile.The role played by Pinus nigra subsp.laricio (Poir.)Maire, which has been used extensively in reforestation throughout Europe and in other continents, was quite important.Although the tensile strength of Pinus nigra subsp.laricio (Poir.)Maire is modest, root cohesion is quite high, as displayed by maximum values of RAR.Thus, Pinus nigra subsp.laricio (Poir.)Maire can significantly contribute to soil slope stabilization and, providing that appropriate forest management practices are used, could also allow for recolonization by other late-successional trees and woody shrubs, further stabilizing the soil, considering that diverse tree root types greatly contribute to slope stability [65].Knowledge of the behavior of root systems can be helpful for designing and using bioengineering techniques [16,56,66] and can be useful for mapping slope stability in wooded areas, particularly in Mediterranean environments.Moreover, the results of our experimental data can be used in slope stability models, opening up new research scenarios, especially by including root cohesion in models of mountain stability [6,67].Root properties such as tensile strength and root morphology are important biomechanical traits that could be efficiently utilized in nature-based solutions to improve soil stability in the region.In this context, forest management in sustainable ways could play a key role in improving slope stability [68] and enhance the ecosystem services associated with forests [69].Therefore, managing forests against mass wasting hazards is one of the most important ecosystem services provided in mountain regions [68,70].

Figure 1 .
Figure 1.Study site located in South Italy (Calabria Region).Red line marks the "Bonis Basin".

Figure 1 .
Figure 1.Study site located in South Italy (Calabria Region).Red line marks the "Bonis Basin".

Figure 2 .
Figure 2. Prepared root samples and tensile strength measurements.(a) Root terminals cast in synthetic resin.Tensile force on roots was measured in the z-direction; d = diameter measurement; L0 = length before test.(b) Root secured between clamping jaws before tension test.1: Wedge tensioner, lower part; 2: Wedge tensioner, upper part; and 3: Vice clamps.

Figure 2 .
Figure 2. Prepared root samples and tensile strength measurements.(a) Root terminals cast in synthetic resin.Tensile force on roots was measured in the z-direction; d = diameter measurement; L 0 = length before test.(b) Root secured between clamping jaws before tension test.1: Wedge tensioner, lower part; 2: Wedge tensioner, upper part; and 3: Vice clamps.

Figure 4 .
Figure 4. Tr and D relationships for eight different species fitted to a power law function: (a) Trees, (b) shrubs.

Figure 4 .
Figure 4. Tr and D relationships for eight different species fitted to a power law function: (a) Trees, (b) shrubs.

Figure 6 .
Figure 6.Root cohesion distributions of the tree (a) and shrubs species (b).

Figure 6 .
Figure 6.Root cohesion distributions of the tree (a) and shrubs species (b).

Table 1 .
Main characteristics of the study tree and shrub species.

Table 2 .
Mean diameter and height of study tree species.

Table 3 .
Moisture content for fresh water content (FWC) and maximum water content (MWC) in the studied tree and shrub species.

Table 4 .
Descriptive statistics for diameter and tensile strength of the tested roots of different species.

Table 5 .
Coefficients of the power law equation, with statistical parameters and coefficient of correlation.

Table 6 .
Mean and standard deviation values of Root Area Ratio (RAR) for each species.

Table 6 .
Mean and standard deviation values of Root Area Ratio (RAR) for each species.