Effects of Different Land Uses (Abandoned Farmland, Intensive Agriculture and Forest) on Soil Hydrological Properties in Southern Spain

A detailed knowledge of soil water repellency (SWR) and water infiltration capacity of soils under different land uses is of fundamental importance in Mediterranean areas, since these areas are prone to soil degradation risks (e.g., erosion, runoff of polluting compounds) as a response to different hydrological processes. The present study evaluates the effects of land uses on SWR and soil hydraulic conductivity (SHC) by direct measurements at the plot scale in three areas representing (1) intensive agricultural use, (2) abandoned farmland, and (3) a forest ecosystem in Southern Spain under Mediterranean climatic conditions. The physico-chemical properties and water content of the experimental soils were also measured. Significant SWR and SHC differences were found among the analyzed land uses. Forest soils showed high SWR and low SHC, while the reverse effects (that is, low SWR and high SHC) were detected in soils subjected to intensive agriculture. Organic matter and bulk density were important soil properties influencing SWR and SHC. The study, demonstrating how land uses can have important effects on the hydrological characteristics of soils, give land managers insights into the choice of the most suitable land use planning strategies in view of facing the high runoff and erosion rates typical of the Mediterranean areas.


Introduction
Mediterranean areas are very prone to soil degradation risks (e.g., surface runoff, erosion, transport of nutrients and other polluting compounds): the soils are generally shallow with low levels of organic matter, low aggregate stability, and nutrient content [1], and the climate is characterized by frequent and intense rainstorms producing a high magnitude of flash floods with high erosive power [2].This combination of soil type and climate leads to a peculiar hydrological response with high runoff rates that have high erosive power.This response is also affected by both land use and soil cover and their temporal and spatial variability, which are considered the most important factors affecting the intensity and frequency of surface runoff and soil erosion [3,4].Inappropriate land use or a poor soil cover may accelerate water runoff and soil erosion dynamics [4,5], leading to unsustainable land degradation processes.The main causes of such negative environmental impacts are agricultural practices, deforestation, overgrazing, land abandonment, wildfires, and civil works [6,7].For instance, agriculture is thought to generate high erosion rates [4,8,9], and the use of land for intensive agriculture, in general, may cause soil damage [10,11].As well known, agriculture, particularly when intensive and subject to frequent tillage operations, may strongly affect the physico-chemical properties of soil, making it more prone to erosion and quality decay.
Erosion problems have also been associated with land abandonment of both agricultural and marginal farmlands [1,12] with particular reference to the Mediterranean of Western Europe [9,13,14].For instance, in SE Spain, where drastic human impacts were recorded in the second half of the past century [15], anthropogenic land use changes have triggered soil erosion and led to severe land degradation [1].Here, since the mid-twentieth century, the use of land for intensive agriculture and the gradual process of vegetation recovery of abandoned farmlands and marginal areas have been the main causes of land degradation [10,11].These land use changes, whose environmental impacts are worsened by the specific climate and the soil fragility, make this region very prone to runoff generation and soil erosion [3] with consequent pollution of water bodies, soil organic matter decline, devastating floods, reservoir siltation, and mass failures [1,15].
In Mediterranean areas, the hydrological processes generating water runoff, soil erosion, and transport of polluting compounds are dominated by the infiltration-excess mechanism [16].In soils of the semi-arid Mediterranean climate, exhibiting low hydraulic conductivity, surface runoff, and soil erosion, can be high.Moreover, such soils could be expected to also be affected by water repellency [17,18], further decreasing infiltration rates [19], which in turn leads to increased runoff and erosion [20,21], to accelerated leaching of agrochemicals [22], and to a reduction in the vegetal cover of soils, leaving the latter bare and thus prone to erosion [23,24].Thus, depending on the water repellency level, Mediterranean soils may have an infiltration rate of up to several orders of magnitude lower than would be expected (e.g., [20,21,25,26]).It is thus evident how a detailed knowledge of the soil water repellency (hereinafter SWR) and water infiltration capacity (in terms of soil hydraulic conductivity, hereinafter SHC) under different land uses is of fundamental importance in Mediterranean areas to control hydrological risks and other environmental impacts linked to these risks.
SWR has been studied worldwide [20,21], in both forest [27] and agricultural soils [24,28].The agricultural soils are usually considered wettable [29] and thus a little subject to SWR, while other studies have demonstrated that some management practices may induce SWR in cultivated soils [30][31][32].The attention paid to surface runoff and soil erosion rates in marginal and abandoned lands (e.g., [33]) has highlighted sometimes contrasting results in the Mediterranean [34] and, consequently, the difficulty to fully understand the effects of land abandonment on hydrological response and soil erosion [9], as modified by both SWR and SHC.Therefore, it is important to know about repellency and infiltration under typical land uses in the drier parts of the Mediterranean basin [18].A better comprehension of these fundamental soil parameters is important for both agricultural production and protection of abandoned farmland from hydro-geological risks [24].
This study evaluates the effects of land use on SWR and SHC by direct measurements at the plot scale in three areas representing (1) intensive agricultural use, (2) abandoned farmland, and (3) a forest ecosystem in Southern Spain.In addition, the physico-chemical properties and water content, which, as it is well known, can influence SWR and SHC, were also measured on these representative Mediterranean soils.The objective of the study is to evaluate which of the analyzed land uses (agriculture, abandoned farmland and forest) shows the highest SWR and lowest SHC in the experimental conditions, also linking these properties to important soil characteristics.By demonstrating how land uses can have important effects on the hydrological characteristics of soils, we want to give land managers insights on the choice of the most suitable land use planning in view of facing the high runoff and erosion rates typical of the Mediterranean areas.

Study Area
The study area is located in the municipality of Villamalea (39.36422 • N, −1.59689 • E, Albacete, SE Spain) (Figure 1).The climate is typically Mediterranean, "Csa", according to the Köppen-Geiger classification [35].The average annual rainfall and temperature are 407 mm and 14 • C, respectively (Figure 2).The rainfall is mainly distributed in spring and autumn, with a long drought in summer that usually lasts from June-September.
Elevation of the studied plots ranges between 760 and 770 m a.s.l., a flat terrain, typical for the Iberian Plateau (Meseta).The plateau is characterized by high elevated (>600 m a.s.l.) region with an undulating (hilly) landscape that is prone to high erosion rates due to the dry climate conditions that reduce the vegetation cover.
The study area is a traditional Mediterranean forest and agricultural land, where pine plantations, natural forest, almonds, olive, cereals, and vineyards are widespread (Figure 1).The natural tree vegetation of the study area mainly consists of Pinus pinea L. and Pinus halepensis M. The main shrubs and herbaceous species found at the study site are Rosmarinus officinalis L., Brachypodium retusum (Pers.)Beauv., Lavandula latifolia Medik., Thymus vulgaris L., Stipa tenacissima (L.), Quercus coccifera L. and Plantago albicans L. Agriculture consists mainly of vineyards, olive orchards and maize crops, with vine production providing the main agricultural income in the study region.
Water 2019, 11 FOR PEER REVIEW 3 of 14

Study Area
The study area is located in the municipality of Villamalea (39,36422° N, −1,59689° E, Albacete, SE Spain) (Figure 1).The climate is typically Mediterranean, "Csa", according to the Köppen-Geiger classification [35].The average annual rainfall and temperature are 407 mm and 14 °C, respectively (Figure 2).The rainfall is mainly distributed in spring and autumn, with a long drought in summer that usually lasts from June-September.
Elevation of the studied plots ranges between 760 and 770 m a.s.l., a flat terrain, typical for the Iberian Plateau (Meseta).The plateau is characterized by high elevated (>600 m a.s.l.) region with an undulating (hilly) landscape that is prone to high erosion rates due to the dry climate conditions that reduce the vegetation cover.
The study area is a traditional Mediterranean forest and agricultural land, where pine plantations, natural forest, almonds, olive, cereals, and vineyards are widespread (Figure 1).The natural tree vegetation of the study area mainly consists of Pinus pinea L. and Pinus halepensis M. The main shrubs and herbaceous species found at the study site are Rosmarinus officinalis L., Brachypodium retusum (Pers.)Beauv., Lavandula latifolia Medik., Thymus vulgaris L., Stipa tenacissima (L.), Quercus coccifera L. and Plantago albicans L. Agriculture consists mainly of vineyards, olive orchards and maize crops, with vine production providing the main agricultural income in the study region.According to FAO-UNESCO [36] and the IUSS Working Group WRB [37], the soils of the study area are classified as Calcic cambisols.Organic matter content is on average 1.1 g/m 3 , varying from 1.05-1.12g/m 3 .The texture is sandy clay loam (60% sand, 10% silt, and 27% clay) over a limestone parent material.

Plot Description
In autumn 2018, nine 10 × 10 m 2 plots were set up in each of the three land uses (27 plots in total) of the same study area with the typical climate and soil characteristics described in Section 2.1.The analyzed land uses were (1) intensive agriculture; (2) abandoned farmland; (3) natural forest.The plots with similar slopes (1-5%) were distributed selecting certain site characteristics, slopes, and aspects to ensure comparability among the 27 plots.Distances between plots were always over 500 m.
The land use related to intensive agriculture (hereinafter "IA") consisted of an olive grove about 20 years old.Cropping operations follow the usual standards of the local farmers.The orchard is tilled three times per year, scheduled depending on the soil moisture and weed removal needs.Tillage is carried out just before or on the occasion of weed germination in order to keep the soil surface clean and tidy and is refined by eliminating herbs near the trunk of each tree (in June or July) by hoes.No fertilizers are used in the olive grove.Residues of pruning, operated in February or March, are concentrated and burned at the margins of the grove, but not chipped, as usually made for other crops (e.g., citrus and apricots).
Forest plots (hereinafter "FO") are covered by Pinus halepenis M., about 50 years old, with a density of 560 trees per ha.Shrubs found underneath the trees are Quercus coccifera L. and Macrochloa tenacissima L. According to FAO-UNESCO [36] and the IUSS Working Group WRB [37], the soils of the study area are classified as Calcic cambisols.Organic matter content is on average 1.1 g/m 3 , varying from 1.05-1.12g/m 3 .The texture is sandy clay loam (60% sand, 10% silt, and 27% clay) over a limestone parent material.

Plot Description
In autumn 2018, nine 10 × 10 m 2 plots were set up in each of the three land uses (27 plots in total) of the same study area with the typical climate and soil characteristics described in Section 2.1.The analyzed land uses were (1) intensive agriculture; (2) abandoned farmland; (3) natural forest.The plots with similar slopes (1-5%) were distributed selecting certain site characteristics, slopes, and aspects to ensure comparability among the 27 plots.Distances between plots were always over 500 m.
The land use related to intensive agriculture (hereinafter "IA") consisted of an olive grove about 20 years old.Cropping operations follow the usual standards of the local farmers.The orchard is tilled three times per year, scheduled depending on the soil moisture and weed removal needs.Tillage is carried out just before or on the occasion of weed germination in order to keep the soil surface clean and tidy and is refined by eliminating herbs near the trunk of each tree (in June or July) by hoes.No fertilizers are used in the olive grove.Residues of pruning, operated in February or March, are concentrated and burned at the margins of the grove, but not chipped, as usually made for other crops (e.g., citrus and apricots).
Forest plots (hereinafter "FO") are covered by Pinus halepenis M., about 50 years old, with a density of 560 trees per ha.Shrubs found underneath the trees are Quercus coccifera L. and Macrochloa tenacissima L.

Soil Property Measurements
Herbal cover, rock fragments, dead woody matter, and bare soil covers were measured at three 5-m transects in each plot, measuring the percent cover of each property over a grid of 1 m × 1 m.
Water 2019, 11, 503 5 of 14 Bulk density was calculated on triple samples per plot as the weight of soil in a given volume of the core extracted by a small cylinder at a depth of 20-30 cm.Soil water content (SWC) was estimated using a HOBOnet Soil Moisture Sensor, which integrates the field-proven ECH2O™ EC5 Sensor and provides readings directly in volumetric water content.SWC was measured hourly for 24 h on the occasion of the SWR and SHC surveys and the measures were then averaged.
Based on the parameters above, the carbon-to-nitrogen ratio (C/N) and cation exchange capacity (CEC) were calculated.
The water infiltration capacity of soils was estimated measuring SHC by a MiniDisk infiltrometer (MDI).In more detail, first the measured cumulative infiltration values (I, [m]) were fitted against the measurement intervals (t, [s]), both given by MDI, using Equation (1): and the coefficients C 1 [m/s] and C 2 [m/s 1/2 ] were estimated by interpolation.Coefficient C 1 is related to SHC, and C 2 is the sorptivity [44].Then, SHC (k, [mm/h]) was calculated using the following equation: where coefficient A is a value relating the Van Genuchten parameters (n and α) for a certain soil type to the suction rate (h 0 ) and the infiltrometer disk radius (2.25 cm).Entering the values of n, α, and h 0 (assumed in this study to be equal to −2 cm) of the experimental soils in the table reported in the MDI manual [44], a value of 2.8 was derived for A. Equations ( 1) and (2) were proposed by Reference [45].SWR was measured as follows: 15 drops of distilled water were released, using a pipette, on the soil surface of a 1-m transect, to homogenize the changing soil conditions; the time necessary for drops to infiltrate completely into the soil was measured by a stopwatch.Before measurement, litter cover was removed, and the soil surface was cleaned using a brush.The high number of replicates (ten points per plot) assured the best reliability for this measurement.The method used is recognized as one of the most appropriate for evaluating the SWR degree in field measurements [46].

Statistical Analysis
A one-way ANOVA (using land use as the independent factor) was applied to evaluate the statistical significance of the variations in soil cover, physico-chemical properties, SWR and SHC.All the plots were considered spatially independent.An independent Fisher's minimum significant difference test (LSD) was used for the post hoc analysis comparisons.A p < 0.05 level of significance was adopted.It was not necessary to perform data transformations for the analysis.ANOVA assumes normality and this assumption was checked using QQ-plots.All measured variables were used to perform the principal component analysis (PCA), the latter being based on a Spearman rank correlation matrix, to reduce the dimensionality of the data set.The statistical analysis was performed by version 3.24 of the R Project for Statistical Computing.

Results
The soil physical characteristics as measured for the three land uses are shown in Table 1.The texture of the experimental soils was statistically similar between AF and IA; the texture of the FO soil was significantly different from the other land uses.In general, the soil was finer in the AF and IA land uses (clay loam and loam, respectively) and coarser in FO.Soils subject to IA were more compacted (mean bulk density of 1.25 g cm −3 ), while the FO soils had a lower bulk density (0.83 g cm −3 ) (Table 1).The IA showed the lowest vegetation cover (on average only 1% of the total sampled area).This latter parameter was similar between the other land uses (about 30% and not significantly different between abandoned farmland and FO).On average, only 1% of the soil in the FO plots was bare against more than 70% in IA (Table 1).The percentage of dead woody matter (plant material less than 2 mm in diameter) was statistically significantly different for all land uses, being highest in FO plots, followed by IA and AF soils.
The soil chemical characteristics as measured in the three land uses are shown in Table 2.The values of pH (on average 8.5-8.7) and EC (0.21-0.24 mS cm −1 ) were not significantly different among the investigated land uses.The OM content, in tune with the soil total carbon, was similar between AF and FO (mean values between 1.3 and 1.5%), and significantly lower in IA (on average 0.6%).OM was not correlated to dead woody matter, presumably due to the different mineralization levels of the investigated soils, on which the previous and current tillage practices may have played a role.Also, the contents of P, K, and cations (Na + , Ca ++ and Mg ++ ) were significantly different between land uses.FO soils had the lowest concentrations (except for Ca ++ ) and the AF the highest (except for P).Conversely, there were no statistically significant differences in N content (on average 0.05-0.1%).Due to this, the C/N ratio was quite similar among all the studied land uses (mean values between 11.9 and 16.3), in spite of the differences recorded in C mean contents (from 0.6% of IA to 1.5% of AF).Finally, CEC was significantly higher in AF (on average 16.3 meq/100 g of soil), while FO soils showed the lowest value (6.8 meq/100 g, this latter not significantly different from IA, 9.9 meq/100 g) (Table 2).The analysis of the hydrological properties of the investigated soils showed much higher mean SWR values on soil surface compared to measurements made at 2 cm below ground (4-fold, AF, and 20-fold, FO).Only in soils under IA, the surface, and below ground SWR were similar.Moreover, FO soils had the highest SWR at both investigated depth (on average 74 s, on the soil surface and 3.5 s, 2 cm below the ground), while the lowest values were measured in IA (2.5 s) on soil surface, and AF (2.0 s) at 2 cm below the ground.Soil water content was 10.0% for AF and FO and 15.2% for IA) at the date of measurements (Figure 3a).
The analysis of the hydrological properties of the investigated soils showed much higher mean SWR values on soil surface compared to measurements made at 2 cm below ground (4-fold, AF, and 20-fold, FO).Only in soils under IA, the surface, and below ground SWR were similar.Moreover, FO soils had the highest SWR at both investigated depth (on average 74 s, on the soil surface and 3.5 s, 2 cm below the ground), while the lowest values were measured in IA (2.5 s) on soil surface, and AF (2.0 s) at 2 cm below the ground.Soil water content was 10.0% for AF and FO and 15.2% for IA) at the date of measurements (Figure 3a).
In soils under IA, the highest SHC was detected (on average 17 mm/h), while the FO soils showed the lowest values (5 mm/h), similar (but significantly different) from SHC measured in the AF (mean value of 7 mm/h) (Figure 3b and Table 3).In soils under IA, the highest SHC was detected (on average 17 mm/h), while the FO soils showed the lowest values (5 mm/h), similar (but significantly different) from SHC measured in the AF (mean value of 7 mm/h) (Figure 3b and Table 3).The PCA evidenced a clear clustering of the three investigated land uses with regard to the properties of soils, assumed as original variables (Figure 4).The PC 1 and PC 2 explained 51% and 37% (in total 88%), respectively, of the total variance of the original variables (Table 4).EC, pH, CEC, and nutrient contents were the variables with the higher loading factors on PC 1 , whereas plot % of vegetation cover, bare soil, dead woody matter, total carbon content, and bulk density showed the higher loadings on PC 2 .
Water 2019, 11 FOR PEER REVIEW 9 of 14

Discussion
Many studies have demonstrated how and by what extent land uses influence the physicochemical and hydrological characteristics of soil (e.g., [4,11,16,34,47]).However, these characteristics also depend on the specific soil properties, such as texture and vegetal cover.For the investigated soils (sandy to clayey texture with variable vegetation cover), it has been shown that some of the chemical characteristics (namely pH, EC, N content and the C/N ratio) were very similar among the investigated land uses.Conversely, significant gradients in the OM content and CEC, two very important soil properties regarding soil fertility and productivity, were noticed when comparing AF

Discussion
Many studies have demonstrated how and by what extent land uses influence the physico-chemical and hydrological characteristics of soil (e.g., [4,11,16,34,47]).However, these characteristics also depend on the specific soil properties, such as texture and vegetal cover.For the investigated soils (sandy to clayey texture with variable vegetation cover), it has been shown that some of the chemical characteristics (namely pH, EC, N content and the C/N ratio) were very similar among the investigated land uses.Conversely, significant gradients in the OM content and CEC, two very important soil properties regarding soil fertility and productivity, were noticed when comparing AF (showing the highest values) to soils subject to IA.This soil depletion in OM and CEC may be attributable to plant uptake, due to the crop growth [48].Moreover, soil under IA has lower OM and higher C/N ratio on average, although not optimal [49]; the other two land uses (AF and FO) have much less favorable C/N ratio for plants.
The most important variations among the investigated land uses were noticed in the soil hydrological properties: SWR and SHC varied significantly among the experimental land uses.More specifically, the topsoil of forested land had a much higher SWR compared to the surface soils subject to IA and AF.Even a slight SWR of FO soil may have noticeable effects on water infiltration rates, generating more runoff-particularly in summer, when the soil is drier-and thus enhancing soil erosion [27].However, although it has been reported that SWR is particularly common under rangeland or forestland such as permanent grassland and deciduous shrub and forest terrain [50][51][52], it can also occur on agricultural land [27].SWR affects both coarse and fine-textured soils [21,29,53] and occurs at low to moderate water contents [26,54].In particular, coarse-textured soils, as those of the forestland analyzed in this study, are more prone to water repellency, even when they are permanently vegetated [55].The level of SWR depends on the soil particle fraction with a hydrophobic surface coating [52] and is influenced by the surface area of the particles, which varies considerably with soil texture [56].Therefore, in the sandy soils of our forestland, which have the lowest surface area, a hydrophobic surface impacts a larger proportion of particles than for a loamy or clayey soil where the surface area is up to three orders of magnitude greater [56,57].Moreover, the plant species surveyed in the forest of this study (eucalyptus and pines) are perennial trees with a considerable concentration of resins, waxes or aromatic oils, which have commonly been associated with SWR [58][59][60].Under these conditions, it may be advisable to develop proper strategies to reduce SWR, such as a more effective soil management, the addition of clays to increase particle surface area, tillage to break-up and abrade hydrophobic surfaces and the use of chemical wetting agents [29,61].Moreover, regions with a Mediterranean climate with prolonged dry periods, such as southern Spain, could be particularly affected by SWR and its hydrological impacts, which bring soils within a water content range in which SWR is exhibited [18,62,63].
Conversely, SWR levels surveyed in the other two investigated land uses (IA and AF) are less severe.In the soils previously subjected to agricultural activities, plant natural succession after land abandonment helps to avoid or reduce SWR, contributing to water penetration into the deeper soil layers, thanks to preferential flow paths via plant roots and stem flow [64].Regarding the plots subjected to IA, agricultural soils are usually considered wettable [29], even though some studies have found that soil management practices can induce water repellency also in cropped areas [30][31][32].In general, tilled sites are virtually unaffected by SWR [55]; furthermore, cultivation promotes rainfall infiltration, and, as a consequence, the runoff and erosion rates are significantly reduced [65].However, although runoff and erosion are mostly reduced after plowing, these processes may increase again over time because of crusting, especially in silty soils [66].
Beside the low SWR levels, the reduced SHC of forestland detected in this study, compared to soils subject to IA, deserves much attention.As a matter of fact, depending on the severity of SWR, such soils may have water infiltration rates lower than would be expected on the basis of their pore size distribution and SWC (e.g., [20,25]).The combination of a reduced wettability (due to high SWR) and low SHC may enhance surface runoff (e.g., [26,67]).With this type of land use, the exposure of soils to the highest and most erosive rainfalls in the experiment area during autumn/winter may aggravate the erosive risks.However, this particular hydrological response of forest sandy soils can be prevented by a proper ground cover and by natural vegetation, which, having a strong influence on soil hydrological properties, reduces soil erosion rates [68].
Also, in the AF of this study, a lower SHC was detected (of the same order of magnitude as FO), which may lead to the degradation of the ecosystem at least during the first 3-5 years after abandonment [65].This reduced SHC could be mainly due to the finer texture of AF soil compared to IA, but also the land use may have had a role.However, it has been demonstrated that, after land abandonment, the vegetation recovery should improve soil water-retention capacity and hydraulic conductivity, thus increasing infiltration and decreasing runoff and erosion rates [65,69].In such a way, plant re-growth starts to control the hydrological and erosional soil response [64].Moreover, the C/N ratio of AF soils is less favorable than in agricultural soils (and this could be due to the absence of fertilizer applications), which subsequently influenced plant succession (depending on the duration of abandonment) [49].
Overall, the worse hydrological response of FO soils compared to IA and AF plots, detected in this study, suggest that, although the forest tree planting by reforestation has been adopted as a viable land management strategy in many parts of the Mediterranean basin, natural scrubland, similar as those of revegetation processes in AF, may be more appropriate with regard to water use efficiency and soil conservation measures, since they assure a lower SWR and a slightly higher SHC [18].
However, it must be understood that the SHC measurements of the tension infiltrometer carried out in this study relate to the unsaturated soil.This parameter may be initially more important than saturated hydraulic conductivity at the time scale of convective rainfall, which is typically short and only lasts 20-60 min, common in many Mediterranean areas [70].Since this investigation showed that soils of forestland and AF are affected by a lower SHC compared to agricultural areas, these land uses require caution in semi-arid or arid Mediterranean ecosystems, where runoff and soil erosion risks may be high.As a matter of fact, given that in the soils typical of this climatic context (in particular, those of prevalent sandy texture) the Hortonian (that is, infiltration-excess) overland flow type dominates over concentrated runoff, a reduction of SHC may worsen the hydrological behavior of soils subjected to these land uses, with a possible increase of water runoff and soil erosion [71].In fact, a lower water infiltration rate increase reduces the water storage of soil during heavy storms, thus increasing the share of precipitation, which generates surface runoff [16,47].

Conclusions
The effects of land use (abandoned farmland, intensive agriculture and natural forest) have been evaluated at the plot scale with particular reference to the physico-chemical properties, water content and repellency, and the hydraulic conductivity of soils.While most soil properties were not significantly different between the land uses, the hydraulic properties of the investigated soils showed specific responses to the different land uses or plant covers.Forest soils showed -high water repellency and low infiltrability, which worsens their hydrological behavior under heavy and frequent storms typical of the Mediterranean landscape.This behavior may increase the risks of soil erosion and pollutant runoff downstream in sloping areas.Abandoned soils previously subjected to agriculture showed a moderate water repellency, but their low hydraulic conductivity can cause serious problems in terms of runoff and soil erosion.However, shrub vegetation recovery, resulting from plant succession, can decrease this concern by increasing soil cover, which may reduce its erodibility.Compared to the forestland and the abandoned land, the agricultural soils were less affected by water repellency and low infiltrability, presumably due to the periodical tillage operations.
Overall, the main conclusion of this study is the important effect of land use on the hydrological characteristics of soils, and indirectly on their different susceptibility to surface runoff and erosion.This suggests paying attention to the specific land use and soil type under the Mediterranean climate (namely, steep sandy soils of forest ecosystems and also in the context of the expected climate changes), which may be affected by high runoff and erosion rates.

Figure 1 .
Figure 1.Location of the study area (A, Villamalea, Castilla La Mancha, SE Spain) with the typical land uses (B: Forest; C: Intensive agriculture; D: Abandoned farmland).

Figure 1 .
Figure 1.Location of the study area (A, Villamalea, Castilla La Mancha, SE Spain) with the typical land uses (B: Forest; C: Intensive agriculture; D: Abandoned farmland).

Figure 2 .
Figure 2. Annual precipitation (total P, mm) as well as mean maximum, average and minimum air temperatures (Tmax, Tmed and Tmin) in the study area (Villamalea, Castilla La Mancha, SE Spain).

Figure 2 .
Figure 2.Annual precipitation (total P, mm) as well as mean maximum, average and minimum air temperatures (T max , T med and T min ) in the study area (Villamalea, Castilla La Mancha, SE Spain).

Figures 3 .
Figures 3. Soil water repellency (SWR, measured on the soil surface and at 2 cm below ground, (a), water content (SWC, b) and hydraulic conductivity (SHC, b) in the study area (Villamalea, Castilla La Mancha, SE Spain).Different lower-case letters indicate statistically significant differences comparing land uses for SWC and SHC following LSD test (P-value < 0.05).

Figure 3 .
Figure 3. Soil water repellency (SWR, measured on the soil surface and at 2 cm below ground, (a), water content (SWC, b) and hydraulic conductivity (SHC, b) in the study area (Villamalea, Castilla La Mancha, SE Spain).Different lower-case letters indicate statistically significant differences comparing land uses for SWC and SHC following LSD test (P-value < 0.05).

Figure 4 .
Figure 4. Principal component analysis applied to properties of soils subject to different land uses (AF: Abandoned Farmland; IA: Intensive Agriculture; FO: Forest) in the study area (Villamalea, Castilla La Mancha, SE Spain) (the ellipses represent the evident clusters achieved by coupling ANOVA and PCA and using land uses as factors).

Figure 4 .
Figure 4. Principal component analysis applied to properties of soils subject to different land uses (AF: Abandoned Farmland; IA: Intensive Agriculture; FO: Forest) in the study area (Villamalea, Castilla La Mancha, SE Spain) (the ellipses represent the evident clusters achieved by coupling ANOVA and PCA and using land uses as factors).

Table 1 .
Soil physical characteristics (mean ± standard deviation) in the study area (Villamalea, Castilla La Mancha, SE Spain).Different lower-case letters indicate statistically significant differences following LSD test (P-value < 0.05).

Table 2 .
Soil chemical characteristics (mean ± standard deviation) in the study area (Villamalea, Castilla La Mancha, SE Spain).Different lower-case letters indicate statistically significant differences following the LSD test (P-value < 0.05).

Table 3 .
Results of ANOVA applied to land uses (Abandoned Farmland; Intensive Agriculture; Forest) for soil water content (SWC), repellency (SWR) and hydraulic conductivity (SHC) in the study area (Villamalea, Castilla La Mancha, SE Spain).

Table 4 .
Factor loadings on the first three principal components (PC) applied to properties of soils subjected to different land uses in the study area (Villamalea, Castilla La Mancha, SE Spain).: in parentheses the percentage of the total variance explained by each PC is reported. Note