Dissolved organic carbon in soil solution of the Jena Experiment (Main Experiment, year 2008)

This data set contains measurements of dissolved organic carbon in samples of soil water collected from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. In April 2002 glass suction plates with a diameter of 12 cm, 1 cm thickness and a pore size of 1-1.6 mm (UMS GmbH, Munich, Germany) were installed in depths of 10, 20, 30 and 60 cm to collect soil solution. The sampling bottles were continuously evacuated to a negative pressure between 50 and 350 mbar, such that the suction pressure was about 50 mbar above the actual soil water tension. Thus, only the soil leachate was collected. Cumulative soil solution was sampled biweekly and analyzed for dissolved organic carbon concentration by a high TOC elemental analyzer (Elementar Analysensysteme GmbH, Hanau, Germany). Samples were analyzed as soon as possible and stored at 4°C if necessary. Often in summer, no free soil solution was available for collection, especially in the upper soil layers. Annual mean values of measured biweekly concentrations of dissolved organic carbon are provided.

(Supplemental Table 2) Plant characteristics, GHG fluxes, and soil chemical and microbial properties in experimental treatments in Alexandra Fiord

We evaluated above- and belowground ecosystem changes in a 16 year, combined fertilization and warming experiment in a High Arctic tundra deciduous shrub heath (Alexandra Fiord, Ellesmere Island, NU, Canada). Soil emissions of the three key greenhouse gases (GHGs) (carbon dioxide, methane, and nitrous oxide) were measured in mid-July 2009 using soil respiration chambers attached to a FTIR system. Soil chemical and biochemical properties including Q10 values for CO2, CH4, and N2O, Bacteria and Archaea assemblage composition, and the diversity and prevalence of key nitrogen cycling genes including bacterial amoA, crenarchaeal amoA, and nosZ were measured. Warming and fertilization caused strong increases in plant community cover and height but had limited effects on GHG fluxes and no substantial effect on soil chemistry or biochemistry. Similarly, there was a surprising lack of directional shifts in the soil microbial community as a whole or any change at all in microbial functional groups associated with CH4 consumption or N2O cycling in any treatment. Thus, it appears that while warming and increased nutrient availability have strongly affected the plant community over the last 16 years, the belowground ecosystem has not yet responded. This resistance of the soil ecosystem has resulted in limited changes in GHG fluxes in response to the experimental treatments.

Aboveground community and species-specific plant biomass from the Jena Experiment (Dominance Experiment, time series since 2002)

This data set comprises a time series of aboveground community plant biomass (Sown plant community, Weed plant community, Dead plant material, and Unidentified plant material; all measured in biomass as dry weight) and species-specific biomass from the sown species of the dominance experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the dominance experiment, 206 grassland plots of 3.5 x 3.5 m were established from a pool of 9 species that can be dominant in semi-natural grassland communities of the study region. In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 3, 4, 6, and 9 species). Plots were maintained by bi-annual weeding and mowing. Aboveground community biomass was harvested twice a year, generally in May and August (in 2002 only once in September) on all experimental plots of the dominance experiment. This was done by clipping the vegetation at 3 cm above ground in two rectangles of 0.2 x 0.5 m per experimental plot. The location of these rectangles was assigned by random selection of new coordinates every year within the central area of the plots (excluding an outer edge of 50cm). The positions of the rectangles within plots were identical for all plots. The harvested biomass was sorted into categories: individual species for the sown plant species, weed plant species (species not sown at the particular plot), detached dead plant material, and remaining plant material that could not be assigned to any category. Biomass was dried to constant weight (70°C, >= 48 h) and weighed. Sown plant community biomass was calculated as the sum of the biomass of the individual sown species. The mean of both samples per plot and the individual measurements are provided in the data file. Overall, analyses of the community biomass data have identified species richness and the presence of particular species as an important driver of a positive biodiversity-productivity relationship.

Mineral nitrogen from KCl extractions of soil from the Jena Experiment (Main Experiment up to 15cm depth, year 2008)

As an estimate of plant-available N, this data set contains measurements of inorganic nitrogen (NO3-N and NH4-N, the sum of which is termed mineral N or Nmin) determined by extraction with 1 M KCl solution of soil samples from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Soil sampling and analysis: Five soil cores (diameter 0.01 m) were taken at a depth of 0 to 0.15 m of the mineral soil from each of the experimental plots in March and October 2008. In October 2008, also the plots of the management experiment, that altered mowing frequency and fertilized subplots (see further details below) were sampled. Samples of the soil cores per plot (subplots in case of the management experiment) were pooled during each sampling campaign. NO3-N and NH4-N concentrations were determined by extraction of soil samples with 1 M KCl solution and were measured in the soil extract with a Continuous Flow Analyzer (CFA, AutoAnalyzer, Seal, Burgess Hill, United Kingdom).

Mineral nitrogen from KCl extractions of soil from the Jena Experiment (Main Experiment up to 15cm depth, year 2006)

As an estimate of plant-available N, this data set contains measurements of inorganic nitrogen (NO3-N and NH4-N, the sum of which is termed mineral N or Nmin) determined by extraction with 1 M KCl solution of soil samples from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Soil sampling and analysis: Five soil cores (diameter 0.01 m) were taken at a depth of 0 to 0.15 m of the mineral soil from each of the experimental plots in March 2006. In October 2006 also the plots of the management experiment, that altered mowing frequency and fertilized subplots (see further details below) were sampled. Measurements from the management experiment are separated into 0 to 0.08 m and 0.08 to 0.15 m. Samples of the soil cores per plot (subplots in case of the management experiment) were pooled during each sampling campaign. NO3-N and NH4-N concentrations were determined by extraction of soil samples with 1 M KCl solution and were measured in the soil extract with a Continuous Flow Analyzer (CFA, AutoAnalyzer, Seal, Burgess Hill, United Kingdom).

Seasonally aggregated values of dissolved inorganic phosphorus in soil solution of the Jena Experiment (Main Experiment, year 2003)

This data set contains measurements of inorganic phosphorus in samples of soil solution collected in 2003 from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below) that have been aggregated to seasonal values. In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Glass suction plates with a diameter of 12 cm, 1 cm thickness and a pore size of 1-1.6 µm (UMS GmbH, Munich, Germany) were installed in April 2002 in depths of 10, 20, 30 and 60 cm to collect soil solution. Manual soil matric potential measurements were used to regulate the vacuum system. Manual soil matric potential measurements were used to regulate the vacuum system. The sampling bottles were continuously evacuated to a negative pressure between 50 and 350 mbar, such that the suction pressure was about 50 mbar above the actual soil water tension. Thus, only the soil leachate was collected. Cumulative soil solution was sampled biweekly and analyzed for dissolved inorganic P (PO4P). Here volume-weighted mean values are provided as aggregated seasonal values (spring = March to May, summer = June to August, fall = September to November, winter = December to February) for 2003 in spring, fall, and winter. To calculate these values, the sampled volume of soil solution is used as weight for P concentrations of the respective sampling date. Inorganic phosphorus concentrations in the soil solution were measured photometrically with a continuous flow analyzer (CFA SAN++, Skalar [Breda, The Netherlands]). Ammonium molybdate catalyzed by antimony tartrate reacts in an acidic medium with phosphate and forms a phospho-molybdic acid complex. Ascorbic acid reduces this complex to an intensely blue-colored complex. As the molybdic complex forms under strongly acidic conditions, we could not exclude the hydrolysis of labile organic P compounds in our samples. Furthermore, the molybdate reaction is not sensitive for condensed phosphates. The detection limits of both TDP and PO4P were 0.02 mg P l-1 (CFA, Skalar).

Seasonally aggregated values of dissolved inorganic phosphorus in soil solution of the Jena Experiment (Main Experiment, year 2005)

This data set contains measurements of inorganic phosphorus in samples of soil solution collected in 2005 from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below) that have been aggregated to seasonal values. In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Glass suction plates with a diameter of 12 cm, 1 cm thickness and a pore size of 1-1.6 µm (UMS GmbH, Munich, Germany) were installed in April 2002 in depths of 10, 20, 30 and 60 cm to collect soil solution. Manual soil matric potential measurements were used to regulate the vacuum system. Manual soil matric potential measurements were used to regulate the vacuum system. The sampling bottles were continuously evacuated to a negative pressure between 50 and 350 mbar, such that the suction pressure was about 50 mbar above the actual soil water tension. Thus, only the soil leachate was collected. Cumulative soil solution was sampled biweekly and analyzed for dissolved inorganic P (PO4P). Here volume-weighted mean values are provided as aggregated seasonal values (spring = March to May, summer = June to August, fall = September to November, winter = December to February) for 2005 in spring, and winter. To calculate these values, the sampled volume of soil solution is used as weight for P concentrations of the respective sampling date. Inorganic phosphorus concentrations in the soil solution were measured photometrically with a continuous flow analyzer (CFA Autoanalyzer [Bran&Luebbe, Norderstedt, Germany]). Ammonium molybdate catalyzed by antimony tartrate reacts in an acidic medium with phosphate and forms a phospho-molybdic acid complex. Ascorbic acid reduces this complex to an intensely blue-colored complex. As the molybdic complex forms under strongly acidic conditions, we could not exclude the hydrolysis of labile organic P compounds in our samples. Furthermore, the molybdate reaction is not sensitive for condensed phosphates. The detection limits of both TDP and PO4P were 0.04 mg P l-1 (Autoanalyzer, Bran&Luebbe).

Mineral nitrogen from KCl extractions of soil from the Jena Experiment (Main Experiment up to 15cm depth, year 2007)

As an estimate of plant-available N, this data set contains measurements of inorganic nitrogen (NO3-N and NH4-N, the sum of which is termed mineral N or Nmin) determined by extraction with 1 M KCl solution of soil samples from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Soil sampling and analysis: Five soil cores (diameter 0.01 m) were taken at a depth of 0 to 0.15 m of the mineral soil from each of the experimental plots in March and October 2007. In March and in October 2007 also the plots of the management experiment, that altered mowing frequency and fertilized subplots (see further details below) were sampled. Samples of the soil cores per plot (subplots in case of the management experiment) were pooled during each sampling campaign. NO3-N and NH4-N concentrations were determined by extraction of soil samples with 1 M KCl solution and were measured in the soil extract with a Continuous Flow Analyzer (CFA, AutoAnalyzer, Seal, Burgess Hill, United Kingdom).

Dissolved phosphorus in soil solution of the Jena Experiment (Main Experiment, year 2002)

This data set contains measurements of dissolved phosphorus (total dissolved nitrogen: TDP, dissolved inorganic phosphorus: PO4P and dissolved organic phosphorus: DOP) in samples of soil water collected in 2002 from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Glass suction plates with a diameter of 12 cm, 1 cm thickness and a pore size of 1-1.6 µm (UMS GmbH, Munich, Germany) were installed in April 2002 in depths of 10, 20, 30 and 60 cm to collect soil solution. Manual soil matric potential measurements were used to regulate the vacuum system. The sampling bottles were continuously evacuated to a negative pressure between 50 and 350 mbar, such that the suction pressure was about 50 mbar above the actual soil water tension. Thus, only the soil leachate was collected. Cumulative soil solution was sampled bi-weekly, in 2002 at the 23.10.2002; 05.11.2002; 20.11.2002; 05.12.2002; and 28.12.2002, and analyzed for dissolved inorganic P (PO4P) and total dissolved phosphorus (TDP). Inorganic phosphorus concentrations in the soil solution were measured photometrically with a continuous flow analyzer (CFA SAN++, Skalar [Breda, The Netherlands]). Ammonium molybdate catalyzed by antimony tartrate reacts in an acidic medium with phosphate and forms a phospho-molybdic acid complex. Ascorbic acid reduces this complex to an intensely blue-colored complex. Total dissolved P in soil solution was analyzed by irradiation with UV and oxidation with K2S2O8 followed by reaction with ammonium molybdate (Skalar catnr. 503-553w/r). As the molybdic complex forms under strongly acidic conditions, we could not exclude the hydrolysis of labile organic P compounds in our samples. Furthermore, the molybdate reaction is not sensitive for condensed phosphates. The detection limits of both TDP and PO4P were 0.02 mg P l-1 (CFA, Skalar). Dissolved organic P (DOP) in soil solution was calculated as the difference between TDP and PO4P. In a low number of samples, TDP was equal to or smaller than PO4P; in these cases, DOP was assumed to be zero.

Dissolved organic carbon in soil solution of the Jena Experiment (Main Experiment, year 2004)

This data set contains measurements of dissolved organic carbon in samples of soil water collected from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. In April 2002 glass suction plates with a diameter of 12 cm, 1 cm thickness and a pore size of 1-1.6 mm (UMS GmbH, Munich, Germany) were installed in depths of 10, 20, 30 and 60 cm to collect soil solution. The sampling bottles were continuously evacuated to a negative pressure between 50 and 350 mbar, such that the suction pressure was about 50 mbar above the actual soil water tension. Thus, only the soil leachate was collected. Cumulative soil solution was sampled biweekly and analyzed for dissolved organic carbon concentration by a high TOC elemental analyzer (Elementar Analysensysteme GmbH, Hanau, Germany). Samples were analyzed as soon as possible and stored at 4°C if necessary. Often in summer, no free soil solution was available for collection, especially in the upper soil layers. Annual mean values of measured biweekly concentrations of dissolved organic carbon are provided.

Plant communities and soil chemistry of shrub-steppe and pine savannah grasslands in southern British Columbia, Canada

Habitat fragmentation alters the edges of remnant habitat patches. We examined changes in the plant community and soil in relation to distance from edge and edge type for shrub-steppe and pine savannah grasslands in southern British Columbia, Canada. Community composition showed significant nonlinear relationships with distance-to-edge more frequently at paved roads and fruit crops than at dirt roads or control sites (i.e., in the interior of grassland patches), with changes typically extending 25-30 m. More exotic species and fewer native species were found near edges, and edges showed decreased cryptogam cover and increased bare ground, especially near paved roads. The soil factors that best predicted compositional changes were soil pH and Cu/Mn at paved roads, soil pH and nitrogen at fruit crops, and soil resistance at dirt roads. Variation partitioning suggested that both direct (e.g., propagule pressure) and indirect (environmental change) factors mediated edge-related community changes, and provided evidence that nonlinear responses at developed edges were not due to natural gradients. Given the range of grassland patch sizes in this region (many patches 1-100 ha), the edge effects we observed represent a considerable loss of "core" habitat, which must be accounted for in conservation planning and site restoration.

Dissolved nitrogen in soil solution of the Jena Experiment (Main Experiment, year 2002)

This data set contains measurements of dissolved nitrogen (total dissolved nitrogen: TDN, dissolved organic nitrogen: DON, dissolved ammonium: NH4+, and dissolved nitrate: NO3-) in samples of soil water collected from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. In April 2002 glass suction plates with a diameter of 12 cm, 1 cm thickness and a pore size of 1-1.6 µm (UMS GmbH, Munich, Germany) were installed in depths of 10, 20, 30 and 60 cm to collect soil solution. The sampling bottles were continuously evacuated to a negative pressure between 50 and 350 mbar, such that the suction pressure was about 50 mbar above the actual soil water tension. Thus, only the soil leachate was collected. Cumulative soil solution was sampled biweekly and analyzed for nitrate (NO3-) and ammonium (NH4+) concentrations with a continuous flow analyzer (CFA, Skalar, Breda, The Netherlands). Nitrate was analyzed photometrically after reduction to NO2- and reaction with sulfanilamide and naphthylethylenediamine-dihydrochloride to an azo-dye. Our NO3- concentrations contained an unknown contribution of NO2- that is expected to be small. Simultaneously to the NO3- analysis, NH4+ was determined photometrically as 5-aminosalicylate after a modified Berthelot reaction. The detection limits of NO3- and NH4+ were 0.02 and 0.03 mg N L-1, respectively. Total dissolved N in soil solution was analyzed by oxidation with K2S2O8 followed by reduction to NO2- as described above for NO3-. Dissolved organic N (DON) concentrations in soil solution were calculated as the difference between TDN and the sum of mineral N (NO3- + NH4+). In 5% of the samples, TDN was equal to or smaller than mineral N. In these cases, DON was assumed to be zero.

Dissolved nitrogen in soil solution of the Jena Experiment (Main Experiment, year 2003)

This data set contains measurements of dissolved nitrogen (total dissolved nitrogen: TDN, dissolved organic nitrogen: DON, dissolved ammonium: NH4+, and dissolved nitrate: NO3-) in samples of soil water collected from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. In April 2002 glass suction plates with a diameter of 12 cm, 1 cm thickness and a pore size of 1-1.6 µm (UMS GmbH, Munich, Germany) were installed in depths of 10, 20, 30 and 60 cm to collect soil solution. The sampling bottles were continuously evacuated to a negative pressure between 50 and 350 mbar, such that the suction pressure was about 50 mbar above the actual soil water tension. Thus, only the soil leachate was collected. Cumulative soil solution was sampled biweekly and analyzed for nitrate (NO3-) and ammonium (NH4+) concentrations with a continuous flow analyzer (CFA, Skalar, Breda, The Netherlands). Nitrate was analyzed photometrically after reduction to NO2- and reaction with sulfanilamide and naphthylethylenediamine-dihydrochloride to an azo-dye. Our NO3- concentrations contained an unknown contribution of NO2- that is expected to be small. Simultaneously to the NO3- analysis, NH4+ was determined photometrically as 5-aminosalicylate after a modified Berthelot reaction. The detection limits of NO3- and NH4+ were 0.02 and 0.03 mg N L-1, respectively. Total dissolved N in soil solution was analyzed by oxidation with K2S2O8 followed by reduction to NO2- as described above for NO3-. Dissolved organic N (DON) concentrations in soil solution were calculated as the difference between TDN and the sum of mineral N (NO3- + NH4+). In 5% of the samples, TDN was equal to or smaller than mineral N. In these cases, DON was assumed to be zero.

Mineral nitrogen from KCl extractions of soil from the Jena Experiment (Main Experiment up to 30cm depth, year 2003)

As an estimate of plant-available N, this data set contains measurements of inorganic nitrogen (NO3-N and NH4-N, the sum of which is termed mineral N or Nmin) determined by extraction with 1 M KCl solution of soil samples from the main experiment plots of a large grassland biodiversity experiment (the Jena Experiment; see further details below). In the main experiment, 82 grassland plots of 20 x 20 m were established from a pool of 60 species belonging to four functional groups (grasses, legumes, tall and small herbs). In May 2002, varying numbers of plant species from this species pool were sown into the plots to create a gradient of plant species richness (1, 2, 4, 8, 16 and 60 species) and functional richness (1, 2, 3, 4 functional groups). Plots were maintained by bi-annual weeding and mowing. Soil sampling and analysis: Five soil cores (diameter 0.01 m) were taken at a depth of 0 to 0.15 m and 0.15 to 0.3 m of the mineral soil from each of the experimental plots in March, June, and October 2003. Samples of the soil cores per plot were pooled during each sampling campaign. NO3-N and NH4-N concentrations were determined by extraction of soil samples with 1 M KCl solution and were measured in the soil extract with a Continuous Flow Analyzer (CFA, Skalar, Breda, Netherlands).

Effects of plant functional traits on soil stability: intraspecific variability matters

This data set describes different vegetation, soil and plant functional traits (PFTs) of 15 plant species in 30 sampling plots of an agricultural landscape in the Haean-myun catchment in South Korea. We divided the data set into two main tables, the first one includes the PFTs data of the 15 studied plant species, and the second one includes the soil and vegetation characteristics of the 30 sampling plots. For a total of 150 individuals, we measures the maximum plant height (cm) and leaf size (cm**2), which means the leaf surface area for the aboveground compartment of each individual. For the belowground compartment, we measured root horizontal width, which is the maximum horizontal spread of the root, rooting length, which is the maximum rooting depth, root diameter, which is the average root diameter of a the whole root, specific root length (SRL), which is the root length divided by the root dry mass, and root/shoot ratio, which is the root dry mass divided by the shoot dry mass. At each of the 30 studied plots, we estimated three different variables describing the vegetation characteristics: vegetation cover (i.e. the percentage of ground covered by vegetation), species richness (i.e. the number of observed species) and root density (estimated using a 30 cm x 30 cm metallic frame divided into nine 10 cm x 10 cm grids placed on the soil profile), as we calculated the total number of roots that appear in each of the nine grids and then we converted it into percentage based on the root count, following. Moreover, in each plot we estimated six different soil variables: Bulk density (g/cm**3), clay % (i.e. percentage of clay), silt % (i.e. percentage of silt), soil aggregate stability, using mean weight diameter (MWD), penetration resistance (kg/cm**2), using pocket penetrometer and soil shear vane strength (kPa).