Introduction

Land degradation is caused by natural and anthropogenic processes threatening the capacity of natural resources to carry out their ecological functions1. Soil degradation is related to a decline in the multiple ecosystem services relying on it, as a result of adverse changes in its biological, chemical, physical and hydrological properties2, affecting the soil capacity to perform a range of ecological functions3.

In the north-western sectors of Murge plateau, an inland part of Apulia region, southern Italy (Fig. 1), land degradation is mainly attributed to the intensive deforestation that started since the roman time due to the large amount of timber removed, which was necessary for the construction of warships4. The original floristic composition of the north-western sectors of Murge plateau was dominated by Quercus dalechampii Ten. and Quercus virgiliana Ten., with the contribution of species belonging to the Quercetalia pubescenti—petraeae order. The Stipo Bromoidis—Quercetum Dalechampii phytocoenosis was completed by Stipa bromoides L., Crataegus laevigata (Poir.) DC., Lonicera etrusca Santi, Carex hallerana Asso, and Iris collina N. Terracc. The ancient forest is today reduced to remnant strips because of the intense exploitation for coppicing, and sheep and goat pasturing5. The pastures not protected by trees were affected by intense erosion and “evolved” into degraded pastures, shallow soils or even outcropping rocks6.

Figure 1
figure 1

Typical landscape of Murge plateau.

In order to solve such critical situation, several reforestations were realized over time utilizing pioneer species, such as Pinus halepensis L., Pinus pinea L. and Cupressus spp., and the nitrogen fixing black locust (Robinia pseudoacacia L.), due to their ability to grow up on shallow and eroded soils. Although the black locust does not enjoy a good reputation as a forest tree because of its invasiveness7,8, it has been utilized extensively for reforestation. These pioneer trees prevent erosion by covering the soil, reduce the impact of rainstorms erosion and produce, over time, litter and root exudates that increase the depth and the organic matter content of soils and can modify the biogeochemical cycles of N and C. These circumstances restore the soil ecological functions to reestablish native forestry species, such as several oaks, walnut, linden, etc., that need higher level of soil quality in comparison to pioneer forestry species.

Since the reforestations have not been conducted in a homogeneous way, in the Murge plateau one can find reforestations that differ in age of plantation and pioneer trees used. Therefore, we selected two sites that diverge for the aforementioned characteristics and investigated their effects on chemical, microbial, and enzymatic soil properties involved in the biogeochemical nutrients cycling, such as microbial biomass, organic matter, soil enzyme activities and soil microbial entropy that are some of the most well-known soil quality indicators9,10.

Materials and methods

Site description and samples collection

The Murgia plateau decreases its elevation from about 680 m a.s.l down to the Adriatic Sea, is characterized by an annual rainfall of about 650 mm11, concentrated mainly in fall and winter, and the rainwater reaches the sea through several incised karst valleys, locally called “lame”, that act as ephemeral streams only during the rainstorms12. From 2004 this area belongs to the “Alta Murgia National Park”.

The first sampling site was Lama Giulia (40°57′12.21′′ N; 16°20′14.75′′ E), belonging to the municipality of Gravina in Puglia, characterized by Calcic Haploxeralf fine loamy mixed thermic soils13. It was reforested with Cupressus arizonica Greene (CYG) during the years 1976–1980 (about 40 years of reforestation at the time of soil sampling), and with Robinia pseudoacacia L. (ROG), whose plantation dates to the end of 1990’s—beginning of 2000’s (about 17 years of reforestation at the time of soil sampling)14. Cupressus arizonica Greene is a conifer belonging to the family of Cupressaceae native to Arizona and New Mexico. It is also used as an ornamental tree, for timber and as a windbreak15. Robinia pseudoacacia L. is a broad-leaf tree belonging to the family of Leguminosae native to the North America. It is a promising fast-growing species for biomass production in short rotation coppice system16. It is also used as ornamental tree and windbreak, and for honey production8. The pasture near CYG and ROG, consisting in degraded but not reforested fields, was considered as control (COG) to check the effects of reforestation. It is mainly composed of shrubs and herbs belonging to four vegetation classes such as Artemisietea Vulgaris, Lygeo-Stipetea, Festuco-Brometea, and Stellarietea Mediae17. This site is close to an abandoned enclosure built with dry stone walls, commonly used in the past for sheltering sheeps and goats.

The other sampling site was located near the Locone Lake (41°04′52.88′′ N; 16°00′54.45′′ E), belonging to the municipality of Minervino Murge, characterized by Typic Calcixerept sandy mixed thermic soils13. It was mainly reforested with Robinia pseudoacacia L. (ROM) in 1985 (about 30 years of reforestation at the time of soil sampling)18. As per Lama Giulia, the pasture near ROM was considered as control (COM). It is currently characterized by arid mediterranean meadows comprising Brachipodium resutum (Pers.) P.Beauv., Brachipodium ramosum Roem. & Schult., Trachynia distachya (L.) P.Beauv., Bromus madritensis L., and Lagurus ovatus L.19. This site shows no signs of farming, even in the past, so grazing is presumably due only to wild herbivores.

Each site was sampled in spring 2017 using an auger, removing the litter and the herbaceous layer, and coring the soil until the abundant root apparatus depth (about 20–25 cm). In details, three composite samples (each made of five subsamples) were cored from CYG, ROG, COG, ROM, and COM, air dried, gently crushed and passed through a 2-mm sieve.

Physical and chemical soil analyses

Soil samples were analyzed according to standard methods20. The soil organic carbon (SOC) content was measured by the Walkley–Black method, and the SOC enrichment was calculated as the difference between the reforested SOC content and the SOC content of the corresponding control soil (pasture). The total nitrogen (TKN) was determined by the Kjeldahl method and the TKN enrichment was calculated as per SOC enrichment. The available phosphorus (Pava) was determined by UV–Vis spectrophotometry according to Olsen method. The exchangeable Ca and K were extracted with a BaCl2 and triethanolamine solution buffered at pH 8.2 and determined using the inductively coupled plasma optical emission spectrometer (ICP-OES) iCAP 6300 (Thermo Electron). The available fractions of Fe, Zn, Cu, Mn and Mo were extracted with diethylenetriaminepentaacetic acid (DTPA) and quantified using the ICP-OES.

Particle-size distribution was determined by the pipette method, and the soil texture was identified using the USDA soil textural classification system13.

Dissolved organic carbon (DOC), soluble carbohydrates and phenolic compounds of soil were analyzed on water extracts obtained at a soil to water ratio 1:10. DOC was determined with a Shimadzu TOC-VCSH analyzer (Kyoto, Japan). Soluble carbohydrates were determined spectrophotometrically by the anthrone method21 and expressed as glucose C equivalents (mg C kg-1). Phenolic compounds were determined by the method of Lowry et al.22 based on the Folin-Ciocalteu reagent and expressed as p-coumaric acid equivalents (mg kg-1).

Soil biochemical analyses

Dehydrogenase activity was measured as the reduction of 2-p-iodophenyl-3-p-nitrophenyl-5-phenyltetrazolium chloride (INT) to iodonitrophenyl-formazan (INTF)23 and expressed as micrograms of INTF produced per soil gram (dry mass) per h (μg INTF g−1 h−1). Microbial biomass C (MBC) content was determined by fumigation of the soil sample with ethanol-free CHCl3 and extraction with 0.5 M K2SO4, according to the method of Vance et al.24 modified by Gregorich et al.25, measured in a Shimadzu TOC-VCSH analyzer (Kyoto, Japan) and expressed in mg C per soil kg.

Alkaline phosphatase activity was determined following the method of Nannipieri et al.26, using p-nitrophenyl phosphate disodium (pNPP, 0.115 M) as substrate for this enzyme. This assay is based on the release and detection of p-nitrophenol (pNP) measured spectrophotometrically. β-glucosidase activity was analyzed according to the method of Tabatabai27, using p-nitrophenyl β-d-glucopyranoside (pNG, 0.05 M) as substrate for this enzyme. Phosphatase and β-glucosidase activities were expressed in micromoles of p-nitrophenol (pNP) produced per soil gram (dry weigh) per h (μmol pNP g−1 h−1).

Urease and protease activities were determined in 0.1 M phosphate buffer at pH 7; 1 M urea and 0.03 M N-α-benzoyl-argininamide (BAA) were used as substrates, respectively. Two mL of buffer and 0.5 mL of substrate were added to 0.5 g of the soil sample, which was incubated at 30 °C (urease) or 39 °C (protease) for 90 min. Both enzymatic activities were determined by the amount of NH4+ released26, expressed in units of micromoles of ammonium-N produced per soil gram (dry mass) per h (μmol NH4+-N g−1 h−1).

Statistical analysis

All experimental data were tested against the normal distribution of variables (Shapiro–Wilk’s test) and the homogeneity of variance (Bartlett’s test) using the RStudio software. The variables normally distributed and with homogeneity of variances were subjected to one-way ANOVA and HSD test, while results not normally distributed were subjected to Levene’s test to check the homogeneity of variance: data with heterogeneous variance were subjected to the nonparametric Friedman test.

The principal component analysis (PCA) was performed to compare the two sites that had different chemical, microbial, enzymatic and botanical characteristics. It was achieved using XLSTAT software and shown as biplots of scores (sites) and loadings (variables).

The regression analyses were performed using the RStudio software and the goodness of fit was evaluated through R2.

Results

All soils from Lama Giulia showed a silty loam texture, with an average composition of sand, silt and clay of 22.23, 54.04 and 23.73%, respectively, whereas all soils from Locone Lake were sandy (84% sand, 9.1% silt and 6.8% clay, on average).

Table 1 shows the average values of the main chemical parameters of the soil samples from Lama Giulia. Noteworthy, the SOC content was significantly (p = 0.0003) the highest in the soil samples of cypress with an average value of 94.3 g kg−1, followed by Robinia pseudoacacia L. (58.0 g kg−1) and pasture soils (45.1 g kg−1). The SOC enrichment of each reforested soil was linearly correlated with the age of each reforestation (Fig. 2). TKN resembled the SOC trend, with a significantly (p = 0.005) higher value recorded in soil samples from cypress (7.0 g kg−1) than black locust and pasture, that showed 4.9 and 3.9 g TKN kg−1, respectively. As per SOC, the linear correlation between the age of each reforestation and the TKN enrichment of each reforested soil showed that the differences of TKN content between the tree species was highly dependent by the time of planting (Fig. 2). In contrast, the Pava was the significantly highest (p = 0.002) in the soil samples from black locust (390.33 mg kg−1), while in cypress and pasture soils was 14.1 and 32.66 mg kg−1, respectively. Finally, the black locust soil showed numerically and significantly higher content of exchangeable Zn and Cu (p = 0.03), respectively, than the other two soils.

Table 1 Chemical properties (mean ± standard deviation of three replicates) of Lama Giulia soils.
Figure 2
figure 2

Linear regression between the age of each reforestation and the soil organic carbon enrichment. ROG black locust at Lama Giulia), ROM black locust at Locone Lake, CYG cypress at Lama Giulia.

The main chemical characteristics of soils from Locone Lake are reported in Table 2. The reforested soil had significantly higher SOC (p = 0.012) and TKN (p = 0.006) content than pasture, while all other parameters did not show any significant differences.

Table 2 Chemical properties (mean ± standard deviation of three replicates) of Locone Lake soils.

Table 3 reports the microbial and enzymatic properties of Lama Giulia soils. The cypress soil showed significantly higher MBC (p = 0.0001) and DOC (p = 0.003) content than black locust and control soil. In addition, also the MBC/SOC ratio was significantly higher (p = 0.011) in the cypress soil than the other two ones. Finally, the protease activity was significantly higher (p = 0.02) in reforested soils than pasture.

Table 3 Enzymatic activities, microbial biomass C and dissolved organic C of Lama Giulia soils (mean ± standard deviation of three replicates).

The microbial and enzymatic characteristics of Locone Lake soils are reported in Table 4. The soil from black locust showed significantly higher phosphatase (p = 0.016), protease (p = 0.048) and urease (p = 0.03) activities than pasture, together to significantly higher DOC (p = 0.0074) and phenolic compounds (p = 0.043) content. In contrast, the MBC/SOC ratio was significantly lower (p = 0.0061) in reforested soil than pasture.

Table 4 Enzymatic activities, microbial biomass C and dissolved organic C of Locone Lake soils (mean ± standard deviation of three replicates).

The PCA analysis isolated five principal components, but only the first two were considered because they explained 74.2% of the variability (Fig. 3). Along the first principal component F1 there was a clear segregation between the points belonging to Lama Giulia (on the right) and those belonging to Locone Lake (on the left) (Fig. 3a). This result explained 50.15% of variability and it was mainly supported by SOC (load 0.940), TKN (load 0.952), soluble carbohydrates (load 0.813), phenolic compounds (load 0.783) and DOC (load 0.771) content, pH value (load −0.928 and −0.909 for pHH2O and pHKCl, respectively) and phosphatase activity (load 0.886) and, to a lesser extent, by Ca, Fe, Zn, Cu and Mn content (Fig. 3b; Table 5). Along the second principal component F2, the black locus from Lama Giulia (ROG) was clearly segregated from the black locust of the other site (ROM) and from the cypress of the same site (CYG) (Fig. 3a). This result explained 20.06% of variability and it was mainly supported by K content (load −0.740), and by glucosidase (load 0.722), urease (load 0.701) and dehydrogenase (load 0.768) activities (Table 5).

Figure 3
figure 3

Principal component analysis shown as biplots of scores. (a) Sites, (b) variables loadings. ROM black locust at Locone Lake, COM degraded but not reforested soils at Locone Lake, CYG cypress at Lama Giulia, ROG black locust at Lama Giulia, COG degrade but not reforested soils at Lama Giulia.

Table 5 PCA loadings.

Discussion

The two sites were segregated for their SOC content, with Lama Giulia soils showing a general higher SOC content than Locone Lake site. The larger SOC content of the cypress soil than black locust soils could be the result of about forty years of litter deposition occurred in the former soil with respect to the about thirty and seventeen years of black locust litter deposition at Locone Lake and Lama Giulia sites, respectively. In addition, the quality of the plant residues could have also played a role in C sequestration since the litter contribution is related to the type of plant28 and, in particular, black locust litter has a relatively higher input of aliphatic carbon, while the conifer litter composition is more aromatic and shows a slower decomposition29.

The Lama Giulia site differed from Locone Lake one for the TKN content too. Although a higher TKN content was expected in black locust soils because of the biological N fixation (BNF), the concentration of TKN was the highest in the cypress soil and mainly related to the litter deposition over time, since 90–95% of soil N is organic. Certainly, cypress also establishes relationships with beneficial microorganisms of the rhizosphere, while the BNF can be influenced by environmental conditions. Wei et al.30 reported that R. pseudoacacia L. forms nodules with many different bacterial species and this can explain the environmental success of black locust, but Schulze et al.31 and Zahran32 reported that the symbiotic association and the BNF are sensitive to soil water availability and drought stress. For example, Veste and Kriebitzsch33 demonstrated that, under long-term drought stress, the leaf N content of black locust reduced due to the negative influence of extreme drought on BNF. Even though Robinia pseudoacacia L. can grow in semi-arid environments, such as the Apulia region ones, it originated from climatic regions with annual rainfall of 1020 to 1830 mm34. Therefore, since in Apulia the balance between rainfall and evapotranspiration is frequently negative during late spring and summer, the drought may have affected the BNF.

The pH values of soils from Lama Giulia were lower than those from Locone Lake and inversely proportional to their SOC content since the soil organic matter offers a large source of acidic protons in the form of carboxyl groups that can determine a decrease of pH values35.

According to the PCA analysis, Locone Lake site showed a general lower phosphatase activity than Lama Giulia and, within each site, no significant difference was found in Lama Giulia (Table 3), while the black locust of Locone Lake showed significantly larger phosphatase activity than the corresponding pasture (Table 4). Since the phosphatase activity is mostly promoted by higher content of soil organic P and total N, and lower pH, the results are in accordance with the trend of those variables between the two sites and within each site. However, Pava has been demonstrated that did not correlate with phosphatase activity36. The highest content of Pava recorded in ROG almost certainly was related to the grazing cattle present in that area which enriched the site with their excreta, since black locust is also a food source for livestock and wild herbivores. This result is apparently confirmed also by the significantly higher content of Zn and Cu in the black locust soil within the Lama Giulia site (ROG), since those elements are often related to the animal droppings.

The two sites were segregated also by the soluble carbohydrates, phenolic compounds and DOC that represent labile fractions of the soil organic matter directly affected by the aboveground and belowground litter inputs. In general, the significantly largest DOC content and the numerically higher carbohydrates and phenolic compounds contents of the cypress soil suggest that the quality of the cypress litter has a better effect than black locust and pastures litters to maintain a pool of labile organic C that provide a source of energy for soil microbial activities. In fact, the chemical composition of the litter species such as lignin, polyphenols, and the C-to-N and lignin-to-N ratios control the litter decomposition and the environmental conditions for soil microbiota37. Indeed, the majority of enzymatic activities are correlated with SOC and the easily degradable organic compounds in soil38,39, which stimulate the growth of soil microbial community.

The PCA analysis indicated, even to a lesser extent, that CYG and ROM soils shared a slightly higher dehydrogenase activity than ROG. Dehydrogenase activity is an index of the global soil microbial metabolic activity and health status of the soil ecosystem40. Dehydrogenase is exclusively an intracellular enzyme that is only present in living microbial cells and therefore is a well-founded marker for monitoring changes in soil metabolism41. As the ecosystems develop with time, there is a reduction of metabolic activity and tendency to reach energy equilibrium in the system42,43, therefore the results pointed out that all soils showed different environmental conditions that modulated the changes in the overall soil metabolism. This result was confirmed also by the activities of the extracellular enzymes β-glucosidase and urease, involved in organic matter turnover and in soil N metabolism, respectively44. In fact, they showed lower activities in ROG soil than the other tree soils.

The protease showed a significantly tree-related stimulation of its activity over the pasture soils within each site (Tables 3 and 4), in response to the higher TKN content of such soils than corresponding pastures that stimulated the synthesis of this enzyme by soil microbial community45.

The MBC is a reliable indicator of changes in soil quality reflecting the long-term effects of soil-vegetation system and it was positively correlated to the SOC contents (R2 = 0.84). Even if the MBC did not segregate the two sites, within the Lama Giulia one it reached the significantly highest value in the cypress soil because of the longer inputs of litter that increased the food resources for microbes and the amount of soil microorganisms. In addition, De Marco et al.46 found that black locust soils had low MBC and respiratory response because of the presence of inhibitory substances, such as 4-hydroxyacetophenon, in black locust litter extracts, that negatively affects the abundance and activity of the soil microbial community and the microbial-derived organic matter formation.

Soil microbial entropy, defined as the ratio of soil MBC to SOC, is a sensitive indicator of the soil quality, assessing the SOC quality and its turnover rate47. It reflects the percentage of the active SOC in the soil ecosystem and the changes of soil environment48. MBC generally account for 1–4% of the total SOC49. Comparing the two pastures, the higher soil microbial entropy of COM with respect to COG was possibly due to the sandy texture of the former soil, that ensured a higher SOC turnover. For the same reason, ROM had higher MBC/SOC ratio than ROG. Within Lama Giulia site, CYG showed the slightly but significantly highest soil microbial entropy because of the amount of litter deposition and, to a lesser extent, the quality of the litter.

Conclusions

The effects of Robinia pseudocacia L. on Locone Lake were remarkable, with an increase of SOC and TKN content in comparison to the corresponding pasture. Comparing the two black locust plantations, the different soil texture between the two sites determined higher values of such parameters in Lama Giulia than Locone Lake, despite the latter was reforested earlier than the former site. In fact, the Locone Lake’s sandy soils showed a higher microbial entropy and turnover, especially in the not reforested soil, suggesting a less conservative soil state than Lama Giulia’s silty loam soil.

Within the Lama Giulia site, CYG soil showed the highest values of some chemical and microbial parameters due to the longest time of reforestation. In contrast, almost all soil parameters studied were not statistically different between the ROG and the corresponding pasture, suggesting that about 17 years of reforestation with black locust apparently are not sufficient to establish general better soil quality.