*
 

iForest - Biogeosciences and Forestry

iForest - Biogeosciences and Forestry
*

Impact of deforestation on the soil physical and chemical attributes, and humic fraction of organic matter in dry environments in Brazil

iForest - Biogeosciences and Forestry, Volume 15, Issue 6, Pages 465-475 (2022)
doi: https://doi.org/10.3832/ifor4016-015
Published: Nov 18, 2022 - Copyright © 2022 SISEF

Research Articles

Deforestation of Caatinga and inadequate land use of these dry environments have impacted soil quality in Northeastern Brazil. The objectives of this study were: (a) to evaluate the effect of deforestation and different agricultural uses on the physical and chemical properties of soil, and humic fractions of soil organic matter in dry environments; and (b) to detect the soil properties that were most affected by anthropic actions. We evaluated four dry areas in Chapada do Araripe, NE Brazil: preserved native vegetation; degraded native vegetation; cassava conventional cultivation; and eucalyptus agro-energy cultivation. Soil fertility, total organic carbon and humic fractions of soil organic matter were lower in the degraded native vegetation area. The best indicators for soil quality evaluation were: macroporosity; bulk density; soil resistance penetration; sum of bases (mainly Ca2+); available P; and saturation by Al3+. Total organic carbon and humic acid fractions of soil organic matter were important in improving soil quality. These properties were influenced by deforestation and agricultural uses, suggesting that the deforestation of native vegetation in dry environments has high capacity to degrade the soil, preventing its regeneration.

Degraded Soil, Forest Soil, Land Use Change, Soil Quality

  Introduction 

Deforestation and inadequate use of soil in semi-arid regions, along with the unfavorable climatic conditions of these regions, poorly developed soils and low biomass production, influence the regenerative capacity of the vegetation and soil degradation. Therefore, the use of more efficient and sustainable management systems is necessary for mitigate these negative effects on vegetation and soil in these regions ([39], [40]). Mekuria et al. ([26]) highlighted the importance of studies evaluating soil degradation in semi-arid regions, though research on this subject is still incipient.

The gypsum mining in the Araripe region of the Brazilian semi-arid region presents a high biomass demand for the calciners, which is responsible for the high rates of deforestation, making this region susceptible to desertification and reducing the productive capacity of soils ([14]). According to Araújo Filho et al. ([3]), soil degradation in the Brazilian semi-arid region has been influenced by the removal of vegetation cover. Studies by Vieira et al. ([48]) indicated that the semi-arid region of the Northeast of Brazil is subject to a high degree of desertification, and the inadequate removal of vegetation cover is one of its main causes.

Changes in the environmental landscape and land use in the semi-arid region are also associated with agricultural expansion. Several studies showed that the agricultural uses in this region have caused soil degradation ([5], [39], [24]). Cassava is the main subsistence crop in Chapada do Araripe (Pernambuco State, NE Brazil) due to its widespread use in family agriculture. This agricultural use of the soil, together with the harvest of the native vegetation as a source of energy for the calcination of gypsum from mining have been the main causes of deforestation of the Caatinga.

The impact on forest resources in the Araripe region has led to the use of strategic alternatives to contain deforestation ([32]). Afforestation with eucalyptus plantations for biomass production to replace the native Caatinga vegetation has been one of these alternatives. Despite these environmental benefits, there is still little scientific information on the impact of the use of this exotic species on the quality of the soil in the region. Indeed, assessing these impacts and establishing an adequate management plan for land uses in this region can directly contribute to soil recovery and conservation, and can indirectly improve the quality of life of local people relying on these natural resources. Monitoring the impact of different land uses and assessing soil quality are keys for early diagnosis of environmental disturbances ([7]).

In this study, we hypothesized that the deforestation of Caatinga for cultivation of cassava monoculture or eucalyptus plantations for use as biomass energy may be the cause of environmental degradation in different soil properties responsible for maintaining the quality and sustainability of the environment.

Recently, soil degradation in dry environment has expanded in the Chapada do Araripe region, mainly due to the intense exploitation of native vegetation, thus making it essential to study the physical and chemical aspects of soil quality. In this study, we meant to answer the following questions: (i) which soil properties are being most impacted by deforestation? (ii) What is the influence of subsistence family farming (cassava cultivation) on soil properties? (iii) Can the establishment of alternative energy crops, such as eucalyptus, affect soil properties and reduce its quality?

  Material and methods 

Study area

The study was carried out in top areas of Chapada do Araripe (07° 27′ 32″ S and 40° 24′ 55″ W), located in the extreme west of the state of Pernambuco, in the municipality of Araripina, in an area of gypsum mining and production. According to the classification of Koppen ([2]), the climate of the region is warm semi-arid, steppe-type with summer-autumn rainfall. Mean annual precipitation is 77 mm and mean annual temperature over the last 30 years is 24 °C. The topography is flat, with altitude varying from 828 to 837 m a.s.l. The soil of the study areas was classified as Dystrophic Yellow Latosol ([31]), corresponding to Oxisol ([38]).

The native vegetation of the top areas of Chapada do Araripe includes the following dominant species: Guapira opposita (Vell.) Reitz; Croton limae A.P.S. Gomes, M.F. Sales & P.E. Berry; Metrodorea mollis Taub.; Pilocarpus spicatus subsp. aracatensis Kaastra; Annona leptopetala (R.E.Fr.) H. Rainer; Senegalia langsdorffii (Benth.) Seigler & Ebinger; Swartzia psilonema Harms; Zanthoxylum hamadryadicum Pirani; Byrsonima gardneriana A. Juss.; Piptademia viridiflora (kunth) Benth. ([30]).

The study was carried out in four distinct though adjacent areas: (i) preserved native vegetation area (PNV); (ii) degraded native vegetation area, with low natural regeneration (DNV); (iii) cassava conventional cultivation area (CCC), with plantations established 46 years ago; and (iv) eucalyptus agro-energetic cultivation area (EAC), aged only 11 years (Fig. 1). Information on land use history and management used in each study area was obtained through interviews to old farmers and local research institutions (Tab. S1 in Supplementary material).

Fig. 1 - The studied sites. (a) Preserved native vegetation (PNV); (b) degraded native vegetation (DNV); (c) Cassava conventional cultivation (CCC); (d) Eucalyptus agro-energy cultivation (EAC).

  Enlarge/Shrink   Download   Full Width  Open in Viewer

Data collection

The number of samples needed to accurately represent each soil property (sample adequacy) was calculated as follows ([37] - eqn. 1):

\begin{equation} N = { \left (\frac{t\alpha \cdot CV} {f} \right )} ^ {2} \end{equation}

where N is the number of samples; tα is the t value from Student distribution for a probability α (p<0.05) and N - 1 degrees of freedom; f is the angular error in relation to the mean; and CV (=100 · standard deviation / mean) is the coefficient of variation of the soil property.

Physical properties

Sampling for determination of the physical properties was performed at the four vertices and in the center of a 100 × 100 m square plot located in the center of each study area, at depths of 0-10 cm and 10-20 cm, totaling ten samples per area (40 samples in total). Samples were collected in volumetric rings with 5 × 5 cm (height × diameter) without structure deformation. The value for the physical properties was 2.7764 for 4 degrees of freedom (5 samples) and probability α (p<0.05). The CV of each physical property was calculated and f values of up to 50% were used, as physical properties are less variable than chemical properties and SOM fractions.

The properties analyzed were total porosity (TP), macroporosity (MaP), microporosity (MiP), field capacity moisture content (FC), wilting point moisture content (WP), soil resistance penetration (SRP) and bulk density (BD). Available water (AW) content was calculated by the difference between FC and WP. A preliminary analysis revealed small and negligible variation in soil texture and particle size distribution among the different areas (Tab. S2 in Supplementary material).

Samples were saturated for 24 h for determination of TP, MaP and MiP. After the saturation period, samples were submitted to the tension of 6 KPa using a tension table for MiP determination. MaP was obtained by the difference between TP and MiP. To determine BD, samples were placed in an oven at 105 °C until constant weight. To obtain the FC, samples were placed on the tension table at 10 KPa, and to obtain the WP, samples were placed in a Richards extractor at 1.500 KPa ([12]). Samples were weighted after equilibration at 80 KPa. Then, determination of SRP was performed using an electronic penetrometer with a needle to simulate the penetration of the root at a speed of 1 cm min-1. The base of the cone was 4 mm thick ([9]).

Chemical properties

Sampling for determination of chemical properties, TOC and SOM fractions was performed within the 100 × 100 m square plot in the center of each study area. Twenty-five samples were collected within each plot, spaced 25 × 25 m at depths of 0-5 cm, 5-10 cm and 10-20 cm, totaling 75 samples per area (300 samples in total). Samples were collected with a Dutch type auger with a deformed structure. The value for these properties was 2.0639 for 24 degrees of freedom (25 samples) and probability α (p<0.05). The CV of each property was calculated and f values of up to 30% were considered ([18]).

The collection of undisturbed samples to evaluate the physical properties of the soil was carried out on the side of a small soil profile by introducing a volumetric ring of 5 cm in diameter. In this case, soil sampling was carried out at depths 0-10 cm (topsoil) and 10-20 cm (subsoil), as there was not enough space to collect samples at the same depths used for the evaluation of soil chemical properties (0-5 cm; 5-10 cm; and 10-20 cm - see above). Additionally, we sampled the entire A horizon and a fraction of the AB horizon, as the most effective absorption capacity of the root system is up to 20 cm deep.

The following chemical properties were determined on soil samples: pH in H2O; exchangeable Ca2+; Mg2+; K+; Al3+; potential acidity (H+Al); and available P. The pH was measured with a glass electrode in a 1:2.5 ratio of the soil solution in distilled water. The Ca2+, Mg2+ and Al3+ were extracted with KCl 1.0 mol L-1. Ions Ca2+ and Mg2+ were measured by atomic absorption spectroscopy and Al3+ was measured by titration in the presence of the blue bromothymol indicator and titrated with NaOH (0.025 mol L-1). Available P and K+ were estimation with Mehlich-1 extraction. P concentration was dosed by colorimetry and K+ content was dosed by flame photometry. The (H+Al) was extracted with 0.5 mol L-1 calcium acetate and titrated with NaOH and phenolphthalein as indicator. All chemical protocols were taken from Teixeira et al. ([46]).

The following indexes derived from the above chemical properties were obtained: sum of bases (SB); potential cation exchange capacity (CECp); saturation by bases (V) and saturation by aluminum (Al saturation).

Humic substances fractions of organic matter

TOC was determined by oxidation of the organic matter by potassium dichromate (K2Cr2O7) 0.020 mol L-1 and titrated with ammonium ferrous sulphate (Mohr’s salt) 0.005 mol L-1, according to Yeomans & Bremner ([50]).

Chemical fractionation of SOM was based on different solubilities of fulvic acid (FA-C), humic acid (HA-C) and humin (HU-C) fractions in acid and alkaline media, according to Swift ([42]). For this, 2.0 g of soil was immersed in 20 mL of 0.1 mol L-1 NaOH for 24 h. Separation between alkaline extract (FA-C + HA-C) and residue was performed by centrifugation at 3000×g for 30 min. The residue was washed again with 20 mL of the NaOH solution and centrifugation again at 3000×g for 30 min. This washing solution was mixed with the extract. The pH of the alkaline extract was adjusted to 1.0 (±0.1) using 20% H2SO4 before decantation for 18 h. The precipitate (HA-C) was separated from the soluble fraction (FA-C) by filtration. The material remaining in the centrifuge tubes was considered as the HU-C fraction. The fractions FA-C and HA-C was determined in 5 mL extract aliquots mixed with 1.0 mL of 0.042 mol L-1 potassium dichromate and 5 mL of concentrated H2SO4 in a block digester at 150 °C (30 min), followed by titration with 0.0125 mol L-1 ferrous ammonium sulfate. In the oven-dried residue, HU-C was determined by adding 5 mL of 0.1667 mol L-1 potassium dichromate and 10 mL of concentrated H2SO4 to a block digester at 150 °C (30 min) and titrating with 0.25 mol L-1 ferrous ammonium sulfate and using ferroin indicator solution. After determinations, the following humification indexes were derived: humic substances fractions (HS-C), ratio between HA-C/FA-C and ratio between HS-C/TOC. HS-C is the sum of HA-C, FA-C and HU-C.

Statistical procedures

Data from physical and chemical properties, TOC and SOM fractions were submitted to analysis of variance by the F test (p<0.05). When there was a significant effect of the different land uses on these properties, the Scott-Knott’s mean test (p<0.05) was applied.

Data were also analyzed by multivariate statistics, through principal component analysis (PCA) and cluster analysis (CA) to select properties that were most influenced by different land uses and to identify similar groups ([35]). The Ward’s algorithm was used as the agglomeration method and the Euclidean distance as measure of dissimilarity ([35]).

Data were centralized and normalized to mean zero and variance one to ensure that all properties equally contributed to the multivariate models used ([35]). In order to select the properties with greater variability, correlations with values of r ≥ 0.70 were considered as significant ([4]). This criterion considers that the number of principal components needed for the interpretation of the results is based on the explanation of at least 70% of the total variation.

  Results 

Soil physical properties

Different land uses had an effect on MaP, TP and BD in the soil surface layer (topsoil: 0-10 cm) and on SRP in both soil layers. The soil physical properties FC, WP AW were only influenced by different land uses in the subsurface layer (subsoil: 10-20 cm - Tab. 1).

Tab. 1 - Soil physical attributes of topsoil (0-10 cm) and subsoil (10-20 cm) under forest, deforested, cropland and agro-energy sites. (PNV): preserved native vegetation; (DNV): degraded native vegetation; (CCC): Cassava conventional cultivation; (EAC): Eucalyptus agro-energy cultivation; (MaP): macroporosity; (MiP): microporosity; (TP): total porosity; (BD): bulk density; (SRP): soil resistance penetration; (FC): field capacity moisture content; (WP): wilting point moisture content; (AW): available water. Mean values and standard deviation are displayed. Values in columns with similar letters are not significantly different (p>0.05); (ns): not significant.

Soil
Depth
Site MaP
(m-3 m-3)
MiP
(m-3 m-3)
TP
(m-3 m-3)
BD
(Mg m-3)
SRP
(Mpa)
FC
(m-3 m-3)
WP
(m-3 m-3)
AW
(m-3 m-3)
0-10 cm PNV 0.17 ± 0.03 a 0.31 ± 0.01 a 0.48 ± 0.02 a 1.27 ± 0.06 b 1.12 ± 0.22 b 0.30 ± 0.01 a 0.09 ± 0.01 a 0.21 ± 0.02 a
DNV 0.08 ± 0.01 b 0.33 ± 0.01 a 0.41 ± 0.02 b 1.51 ± 0.07 a 3.44 ± 0.69 a 0.33 ± 0.01 a 0.10 ± 0.01 a 0.22 ± 0.02 a
CCC 0.13 ± 0.02 a 0.36 ± 0.02 a 0.49 ± 0.02 a 1.34 ± 0.06 b 0.86 ± 0.17 b 0.36 ± 0.02 a 0.10 ± 0.01 a 0.25 ± 0.02 a
EAC 0.16 ± 0.03 a 0.33 ± 0.01 a 0.49 ± 0.02 a 1.30 ± 0.06 b 1.07 ± 0.21 b 0.32 ± 0.01 a 0.10 ± 0.01 a 0.22 ± 0.02 a
Mean 0.15 ± 0.02 0.33 ± 0.01 0.46 ± 0.02 1.35 ± 0.07 1.62 ± 0.32 0.32 ± 0.01 0.09 ± 0.01 0.22 ± 0.02
F-test 4.05* 3.37 ns 4.28* 3.78* 26.16* 3.20 ns 1.00 ns 1.33 ns
10-20 cm PNV 0.12 ± 0.02 a 0.33 ± 0.01 a 0.45 ± 0.01 a 1.34 ± 0.02 a 0.93 ± 0.14 b 0.33 ± 0.01 b 0.07 ± 0.003 b 0.25 ± 0.01 a
DNV 0.08 ± 0.01 a 0.35 ± 0.01 a 0.43 ± 0.01 a 1.43 ± 0.02 a 2.05 ± 0.30 a 0.35 ± 0.01 a 0.10 ± 0.004 a 0.25 ± 0.01 a
CCC 0.10 ± 0.02 a 0.36 ± 0.01 a 0.47 ± 0.01 a 1.41 ± 0.02 a 1.46 ± 0.22 b 0.36 ± 0.01 a 0.10 ± 0.004 a 0.27 ± 0.02 a
EAC 0.10 ± 0.02 a 0.33 ± 0.01 a 0.42 ± 0.01 a 1.42 ± 0.02 a 2.69 ± 0.40 a 0.32 ± 0.01 b 0.10 ± 0.004 a 0.21 ± 0.01 b
Mean 0.10 ± 0.02 0.34 ± 0.01 0.44 ± 0.01 1.40 ± 0.02 1.78 ± 0.26 0.34 ± 0.01 0.09 ± 0.004 0.24 ± 0.01
F-test 1.60 ns 14.28 ns 2.14 ns 3.49 ns 10.52* 16.03* 16.19* 17.73*

  Enlarge/Reduce  Open in Viewer

The DNV area in the topsoil had the lowest values of MaP and TP and the highest values of BD and SRP. The DNV area presented the highest SRP value, being 207, 300 and 221% higher than the PNV, CCC and EAC areas. The PNV, CCC and EAC areas presented MaP values above 0.10 m3 m-3. The DNV area at both depths and EAC at the subsurface layer showed SRP values greater than 2.0 MPa (Tab. 1).

The physical properties that better accounted for data variation were MaP, TP, BD and SRP in topsoil, which were significantly correlated with PC1, showing the influence of the different soil uses on these properties. Regarding PC2, MiP, FC and AW showed the highest and significant correlation values with this axis (Tab. 2).

Tab. 2 - Correlation between principal components (PC1, PC2) and soil physical attributes of topsoil (0-10 cm) and subsoil (10-20 cm). (MaP): macroporosity; (MiP): microporosity; (TP): total porosity; (BD): bulk density; (SRP): soil resistance penetration; (FC): field capacity moisture content; (WP): wilting point moisture content; (AW): available water; (*): p<0.05.

Attribute 0-10 cm 10-20 cm
PC1 PC2 PC1 PC2
MaP -0.983* 0.027 -0.647 -0.568
MiP 0.365 0.930* -0.258 0.966*
TP -0.834* 0.525 -0.735* 0.548
BD 0.991* -0.124 0.816* 0.427
SRP 0.851* -0.520 0.997* -0.074
FC 0.487 0.864* -0.252 0.959*
WP 0.692 0.530 0.935* 0.276
AW 0.234 0.955* -0.695 0.713*
Absolute variance (%) 54.80 40.90 54.10 37.80
Cumulated variance (%) 54.80 95.70 54.10 91.90

  Enlarge/Reduce  Open in Viewer

Regarding the 10-20 cm depth layer (subsoil), the first two PC axes explained 91.90% of the total variance in physical properties among the 4 different land uses. PC1 a 54.10% and presented a high correlation with TP, SRP, BD and WP, while PC2 accountd for 37.80% of the variability and showed a high correlation with MiP, FC and AW (Tab. 2).

In the derived PC biplots (Fig. 2), BD and SRP vectors were opposed to MaP and TP along the first PC axis in the top 0-10 cm layer (Fig. 2a). SRP was also inversely related to MiP, FC and AW in the 10-20 cm depth layer (Fig. 2b).

Fig. 2 - Biplots of the soil physical properties. (a) Topsoil (0-10 cm) and (b) subsoil (10-20 cm). (PNV): Preserved native vegetation; (DNV): Degraded native vegetation; (CCC): Cassava conventional cultivation; (EAC): Eucalyptus agro-energy cultivation. Soil physical attributes: (MaP) macroporosity; (MiP) microporosity; (TP) total porosity; (BD) bulk density; (SRP) soil resistance penetration; (FC) field capacity moisture content; (WP) wilting point moisture content; (AW) available water. The arrow represents the direction of high weighting (maximum variation) of soil physical properties along the principal components. PC1 and PC2 of topsoil (a) explain 95.7% of the variation, while PC1 and PC2 of subsoil (b) explain 91.9% of the total variation.

  Enlarge/Shrink   Download   Full Width  Open in Viewer

The different land uses sampled were distributed in different quadrants of the biplots in both layer of the soil, forming 4 distinct groups (Fig. 2). This shows that the different land uses have very different physical properties of the soil. For example, the DNV land use is characterized by extreme values of DB and SRP in the topsoil (0-10 cm - Fig. 2a) and by extreme values of BD and WP in the subsurface layer (10-20 cm - Fig. 2b). This means that soil bulk density (BD) is largely affected in the degraded native vegetation (DNV) as compared with the preserved native vegetation (PNV). Similarly, MaP and TP in the topsoil (0-10 cm) and SRP in the subsoil (10-20 cm) showed extreme values under the EAC land use (Fig. 2). Moreover, MiP, FC, WP and AW properties in the soil surface layer (Fig. 2a) had extreme values under the CCC land use, as well as TP, MiP, FC and AW in the subsurface layer (Fig. 2b). The preserved native vegetation (PNV) had the higher macroporosity values (MaP) in the subsurface layer (10-20 cm - Fig. 2b, Tab. 1), and in general it showed more balanced physical properties as compared with the other land use types.

Cluster analysis based on soil physical attributes highlighted the existence of three groups of land uses, according to the different types of land use and soil depth (Fig. 3). In the superficial layer, the first group was formed only by the DNV area, the second group was formed only by the CCC area and the third group joined the PNV and EAC areas (Fig. 3a). In the subsurface layer, DNV and CCC areas were grouped and PNV and EAC areas were separated and formed individual groups (Fig. 3b). The more intensive use of the soil in the CCC area showed degradation of its physical properties and evolution towards the levels of the degraded natural vegetation (DNV). The use of eucalyptus (EAC) was more conservative, preserving the favorable physical properties of the preserved area (PNV).

Fig. 3 - Dendrograms of the different land uses analyzed according to the soil physical properties. (a) Topsoil (0-10 cm) and (b) subsoil (10-20 cm).(PNV): preserved native vegetation; (DNV): degraded native vegetation; (CCC): Cassava conventional cultivation; (EAC): Eucalyptus agro-energy cultivation.

  Enlarge/Shrink   Download   Full Width  Open in Viewer

Soil chemical properties

Deforestation and different land uses altered soil chemical properties at all depths (Tab. 3). The pH values were generally higher in the CCC and EAC plots, except in the surface layer (0-5 cm depth), where the highest value was recorded in the EAC. The lowest pH values occurred in the DNV plot. The pH was also low in the PNV area, matching the DNV plot in the 10-20 cm depth layer (Tab. 3). The other components of soil acidity, such as Al3+ and Al saturation were higher in the DNV area at all depths (Tab. 3). The areas under agricultural cultivation had higher soil pH and lower Al saturation values, especially in superficial layers. However, the value of (H + Al) and CECp were higher in the PNV area (Tab. 3).

Tab. 3 - Soil chemical attributes of topsoil (0-5 cm) and subsoil (5-10 and 10-20 cm) under forest, deforested, cropland and agro-energy sites. (PNV): preserved native vegetation; (DNV): degraded native vegetation; (CCC): Cassava conventional cultivation; (EAC): Eucalyptus agro-energy cultivation. Soil chemical attibutes: (pH) hydrogen potential (H2O/1:2.5); (P) available phosphorus; (K+) exchangeable potassium; (Ca2+) exchangeable calcium; (Mg2+) exchangeable magnesium; (Al3+) exchangeable aluminum; (H+Al) potential acidity; (SB) sum of bases; (CECp) cations exchange capacity potential; (V) base saturation; (Al saturation) saturation by aluminum. Mean values and standard deviation are reported. Values in columns with similar letters are not significantly different (p>0.05).

Depth Site pH P
(mg dm-3)
K+
(cmolc dm-3)
Ca2+
(cmolc dm-3)
Mg2+
(cmolc dm-3)
Al3+
(cmolc
dm-3)
H+Al
(cmolc dm-3)
SB
(cmolc dm-3)
CECp
(cmolc dm-3)
V
(%)
Al
saturation
(%)
0-5 cm PNV 4.83±0.25b 1.26±0.63b 0.14±0.03a 0.65±0.02b 0.18±0.08b .0.51±0.13b 6.68±1.87a 0.97±0.29c 7.66±1.98a 13.36±4.80b 35.19±9.04b
DNV 4.49±0.20c 0.95±0.68c 0.06±0.01c 0.36±0.05c 0.07±0.02c 0.66±0.12a 4.15±1.05c 0.52±0.08d 4.67±1.04d 11.70±3.25b 55.70±6.03a
CCC 4.96±0.38b 2.10±0.63a 0.10±0.03b 0.98±0.19a 0.18±0.06b 0.38±0.09c 4.27±1.14c 1.29±0.25b 5.57±1.14c 24.17±6.67a 23.23±7.23c
EAC 5.42±0.51a 2.59±0.89a 0.14±0.07a 1.04±0.38a 0.32±0.18a 0.33±0.17c 4.78±0.79a 1.53±0.57a 6.31±0.79b 24.19±8.51a 19.94±13.43d
Mean 4.92±0.34 1.72±0.71 0.11±0.03 0.75±0.16 0.19±0.09 0.47±0.13 4.97±1.21 1.07±0.30 6.05±1.24 18.35±5.81 33.51±8.93
F-test 25.97* 25.47* 16.24* 40.91* 22.94* 37.72* 24.56* 36.49* 25.32* 30.75* 79.04*
5-10 cm PNV 5.00±0.21b 1.26±0.62b 0.07±0.02b 0.80±0.21b 0.15±0.05b 0.53±0.11b 6.60±1.59a 1.03±0.27b 7.63±1.33a 13.67±2.91b 34.71±7.01b
DNV 4.72±0.12c 0.81±0.66c 0.03±0.01d 0.45±0.06c 0.04±0.01c 0.60±0.10a 4.48±0.94b 0.53±0.07c 5.02±0.97c 10.93±1.99c 52.78±5.02a
CCC 5.23±0.32a 2.04±0.69a 0.05±0.01c 1.04±0.21a 0.18±0.08b 0.36±0.11c 4.98±1.29b 1.29±0.24a 6.27±1.33b 21.16±4.85a 22.37±8.01c
EAC 5.22±0.46a 1.78±0.91a 0.09±0.05a 1.00±0.37a 0.28±0.11a 0.36±0.16c 5.16±1.00b 1.39±0.56a 6.55±1.05b 21.37±7.76a 22.82±12.90c
Mean 5.04±0.28 1.47±0.72 0.06±0.02 0.82±0.21 0.16±0.06 0.46±0.12 5.30±1.21 1.06±0.29 6.37±1.17 16.78±4.38 33.17±8.24
F-test 15.70* 15.13* 14.89* 28.36* 20.95* 27.70* 14.60* 28.91* 17.44* 28.10* 63.27*
10-20 cm PNV 4.83±0.18b 0.56±0.59a 0.06±0.01b 0.28±0.12c 0.15±0.05b 0.59±0.09b 5.41±1.35a 0.49±0.18b 5.91±1.44a 8.39±2.85c 55.90±10.8b
DNV 4.70±0.15b 0.11±0.18b 0.02±0.01d 0.07±0.03d 0.04±0.01c 0.66±0.09a 4.36±1.25b 0.14±0.05c 4.50±1.24b 3.51±1.79d 82.51±5.32a
CCC 4.93±0.29a 0.90±0.59a 0.05±0.01c 0.64±0.33a 0.17±0.05b 0.48±0.10c 4.17±1.09b 0.86±0.38a 5.03±1.13b 17.61±8.96a 38.08±13.08c
EAC 5.04±0.43a 0.79±0.52a 0.08±0.01a 0.45±0.32b 0.25±0.05a 0.51±0.17c 5.15±1.37a 0.79±0.51a 5.95±1.63a 13.11±6.17b 43.17±18.09c
Mean 4.88±0.26 0.59±0.47 0.05±0.01 0.35±0.20 0.15±0.04 0.56±0.11 4.78±1.27 0.57±0.28 5.35±1.36 10.65±4.94 54.91±11.82
F-test 6.15* 11.94* 15.25* 23.97* 20.95* 14.62* 6.13* 22.45* 7.32* 25.41* 60.69*

  Enlarge/Reduce  Open in Viewer

Exchangeable bases (Ca2+, Mg2+ and K+) had the highest levels in the EAC area at all soil depths, with the exception of Ca2+ in the 10-20 cm depth layer, where the highest content was found in the CCC area. These results reflected in the increment of SB, CECp and V of the soil (Tab. 3). The soil of the DNV area presented high deficiency of exchangeable bases and the PNV area was naturally poor in nutrients. The areas under cultivation (CCC and EAC) had the highest P levels in the soil, but in subsoil they were not different from the PNV area (Tab. 3).

In general, we found low soil fertility in all the different land uses. The exchangeable bases, available P and pH were low, and Al saturation was high. CCC or EAC contributed to raise fertility levels and reduce acidity indicators, but the CECp of the PNV area was higher than the areas of cultivation (CCC and EAC). Further, the inappropriate exploitation of the DNV area degraded the soil, reducing its fertility and increasing acidity.

Total organic carbon and humic fractions of organic matter

Deforestation and different land use affected TOC (all depths) and the humic fractions of organic matter HA-C (all depths), HU-C at 0-5 cm depth and HS-C at 0-5 cm and 10-20 cm depths (Tab. 4). The PNV and EAC areas presented the highest levels of TOC in the superficial layer. Reduction of TOC in the CCC and DNV areas was of 10% and 20%, respectively, in relation to the PNV area. In the subsurface layers (5-10 and 10-20 cm), the CCC area accumulated more organic C and TOC contents reached the same levels of PNV and EAC areas (Tab. 4).

Tab. 4 - Total organic carbon and humic fractions of the organic matter of topsoil (0-5 cm) and subsoil (5-10 and 10-20 cm) under forest, deforested, cropland and agro-energy. (PNV): preserved native vegetation; (DNV): degraded native vegetation; (CCC): Cassava conventional cultivation; (EAC): Eucalyptus agro-energy cultivation. (TOC): total organic carbon; (HA-C): humic acid fraction; (FA-C): fulvic acid fraction; (HU-C): humin fraction; (HS-C): humic substances fractions; (HA-C/FA-C): humic acid fraction/fulvic acid fraction; (HS-C/TOC): humic substances fractions/total organic carbon. Mean values and standard deviation are reported. Values in columns with similar letter are not significantly different (p>0.05). (ns): not significant.

Soil
Depth
Site TOC
(g kg-1)
HA-C
(g kg-1)
FA-C
(g kg-1)
HU-C
(g kg-1)
HS-C
(g kg-1)
HA-C/FA-C
(g kg-1)
HS-C/TOC
(g kg-1)
0-5 cm PNV 14.16 ± 2.23 a 3.21 ± 1.73 a 2.82 ± 2.05 a 6.34 ± 1.73 a 12.37 ± 2.88 a 1.14 ± 0.51 a 0.87 ± 0.26 a
DNV 11.26 ± 3.00 c 1.97 ± 1.43 b 2.56 ± 1.59 a 4.81 ± 2.73 b 9.34 ± 3.55 b 0.77 ± 0.34 b 0.83 ± 0.39 a
CCC 12.78 ± 2.56 b 2.33 ± 1.75 b 2.88 ± 1.72 a 5.67 ± 1.80 a 10.88 ± 3.06 b 0.81 ± 0.36 b 0.85 ± 0.30 a
EAC 14.31 ± 2.72 a 3.67 ± 2.47 a 2.84 ± 1.28 a 5.75 ± 2.01 a 12.26 ± 3.48 a 1.29 ± 0.58 a 0.86 ± 0.28 a
Mean 13.13 ± 2.63 2.80 ± 1.85 2.78 ± 1.66 5.64 ± 2.15 11.21 ± 3.24 1.00 ± 0.45 0.85 ± 0.31
F-test 7.70* 3.93* 0.62 ns 2.83* 5.30* 1.67* 0.38 ns
5-10 cm PNV 12.17 ± 2.19 a 3.06 ± 0.82 a 2.92 ± 0.95 a 3.59 ± 2.23 a 9.57 ± 2.96 a 1.04 ± 0.68 a 0.79 ± 0.36 a
DNV 9.48 ± 1.70 b 1.93 ± 0.52 b 2.51 ± 0.73 a 3.11 ± 1.93 a 7.55 ± 2.84 a 0.77 ± 0.80 b 0.80 ± 0.41 a
CCC 11.16 ± 2.00 a 2.04 ± 0.55 b 2.83 ± 0.97 a 3.19 ± 1.79 a 8.06 ± 2.66 a 0.72 ± 0.54 b 0.72 ± 0.34 a
EAC 12.00 ± 2.16 a 3.07 ± 0.83 a 3.05 ± 1.13 a 3.23 ± 2.24 a 9.35 ± 4.62 a 1.01 ± 0.67 a 0.78 ± 0.36 a
Mean 11.20 ± 2.05 2.52 ± 0.68 2.83 ± 0.95 3.28 ± 2.05 8.63 ± 3.27 0.88 ± 0.67 0.77 ± 0.37
F-test 2.49* 2.97* 1.45 ns 0.26 ns 2.18 ns 2.31* 0.19 ns
10-20 cm PNV 7.73 ± 2.08 a 1.90 ± 1.25 a 2.85 ± 1.38 a 1.98 ± 0.95 a 6.73 ± 2.21 a 0.67 ± 0.51 a 0.87 ± 0.35 a
DNV 6.03 ± 1.28 b 0.66 ± 0.80 b 2.20 ± 1.73 a 1.81 ± 0.89 a 4.67 ± 2.01 b 0.30 ± 0.44 b 0.77 ± 0.36 a
CCC 7.49 ± 2.07 a 1.11 ± 0.75 b 2.61 ± 1.53 a 1.83 ± 0.84 a 5.55 ± 1.74 b 0.43 ± 0.86 b 0.74 ± 0.41 a
EAC 7.61 ± 1.73 a 1.24 ± 0.94 b 2.89 ± 2.22 a 1.86 ± 1.02 a 5.99 ± 2.40 a 0.43 ± 0.44 b 0.79 ± 0.33 a
Mean 7.21 ± 1.79 1.23 ± 0.94 2.64 ± 1.72 1.87 ± 0.93 5.74 ± 2.09 0.45 ± 0.56 0.79 ± 0.36
F-test 5.24* 7.89* 0.93 ns 0.19 ns 5.13* 1.82* 0.61 ns

  Enlarge/Reduce  Open in Viewer

The HA-C fractions in the PNV and EAC areas were also higher than in the CCC and DNV areas in the superficial layers (0-5 and 5-10 cm), but the PNV area had the highest C content in the HA-C fraction in the subsurface layer (10-20 cm). The C of the HU-C fraction in the CCC area showed a reduction of 24% in relation to the more preserved area (PNV) at 0-5 cm depth. In the CCC and DNV areas, there was a reduction of the C content on the HS-C fraction in layers of 0-5 cm and 10-20 cm depth (Tab. 4).

In general, a low level of humification (<45%) was observed in the soil of all land uses, as organic C of the HU-C fraction represented on average 43%, 29% and 26% of TOC at 0-5 cm, 5-10 cm and 10-20 cm depths, respectively (Tab. 4). However, the PNV plot showed the highest level of humification (45%) in the surface layer.

TOC reduction with depth ranged from 41% to 47% between the 0-5 cm and 10-20 cm depth layers, with the EAC area showing the greatest reduction (Tab. 4). The litter produced by eucalyptus, concentrating more organic C on the soil surface. The HS-C fraction was reduced in depth by 46%, 50%, 49% and 51% in the areas of PNV, DNV, CCC and EAC, respectively (Tab. 4). The HA-C fraction was reduced with depth approximately 1.5-fold more in the DNV area than in the PNV area (Tab. 4). The stability of the organic C in the DNV area was higher in the superficial layer and in the PNV area the distribution of the more stable fraction of the organic C (HA-C) was more homogeneous. The HA-C/FA-C ratio was higher than 1 in the superficial layers of soil of the PNV and EAC areas (Tab. 4), showing a predominance of humic acids over fulvic acids and confirming the greater stability of SOM.

The first two components of the PCA carried out on soil chemical properties accounted for 96% of the total variation in the 0-5 cm depth layer. PC1 was associated with most of the chemical properties (pH, P, K+, Ca2+, Mg2+, Al3+, SB and Al saturation) and all humic fractions of SOM (TOC, HA-C, FA-C, HU-C, HS-C and HA-C/FA-C), with the exception of HS-C/TOC, while PC2 was associated with properties (H+Al) and V (Tab. 5).

Tab. 5 - Correlation between each principal component (PC) and soil chemical attributes and humic fractions of the organic matter in topsoil (0-5 cm) and subsoil (5-10 cm and 10-20 cm). (pH): hydrogen potential (H2O/1:2.5); (P): available phosphorus; (K+): exchangeable potassium; (Ca2+): exchangeable calcium; (Mg2+): exchangeable magnesium; (Al3+): exchangeable aluminum; (H+Al): potential Acidity; (SB): sum of bases; (CECp): cations exchange capacity potential; (V): base saturation; (Al saturation): saturation by aluminum; (TOC): total organic carbon; (HA-C): humic acid fraction; (FA-C): fulvic acid fraction; (HU-C): humin fraction; (HS-C): humic substances fractions; (HA-C/FA-C): humic acid fraction/fulvic acid fraction; (HS-C/TOC): humic substances fractions/total organic carbon. (*): p<0.05.

Attribute 0-5 cm 5-10 cm 10-20 cm
PC1 PC2 PC1 PC2 PC1 PC2
pH 0.931* -0.325 0.847* -0.508 0.870* -0.393
P 0.794* -0.606 0.703* -0.671 0.887* -0.430
K+ 0.941* 0.332 0.949* 0.151 0.881* 0.192
Ca2+ 0.832* -0.509 0.852* -0.473 0.746* -0.614
Mg2+ 0.960* -0.195 0.934* -0.288 0.935* -0.204
Al3+ -0.879* 0.454 -0.727* 0.687 -0.798* 0.592
(H+Al) 0.410 0.884* 0.563 0.701* 0.553 0.766*
SB 0.904* -0.416 0.914* -0.400 0.858* -0.506
CECp 0.694 0.677 0.742* 0.606 0.875* 0.421
Al saturation -0.890* 0.370 -0.865* 0.469 -0.903* 0.402
V 0.659 -0.737* 0.718* -0.696 0.733* -0.633
TOC 0.954* 0.294 0.936* 0.340 0.981* 0.094
HA-C 0.924* 0.277 0.819* 0.532 0.717* 0.645
FA-C 0.952* 0.233 0.999* 0.010 0.966* 0.236
HU-C 0.798* 0.587 0.498 0.750* 0.513 0.826*
HS-C 0.921* 0.386 0.861* 0.506 0.826* 0.545
HA-C /FA-C 0.873* 0.294 0.679 0.683 0.639 0.684
HS-C /TOC 0.696 0.658 0.650 0.725* 0.253 0.964*
Absolute variance (%) 71.60 24.40 65.20 29.30 63.30 31.00
Cumulated variance (%) 71.60 96.00 65.20 94.50 63.30 94.30

  Enlarge/Reduce  Open in Viewer

Regarding the 5-10 cm depth layer, the first two components explained 94.5% (PC1: 65.20%, PC2: 29.30%) of the variance in soil chemical properties and humic fractions of SOM. With the exception of (H+Al), HU-C, HA-C/FA-C and HS-C/TOC, all other properties were correlated with PC1, while PC2 was associated with (H+Al), HU-C and HS-C/TOC (Tab. 5).

In the 10-20 cm depth layer, PC1 and PC2 were responsible for 94.30% of the total variance (63.30% and 31.00%, respectively). Soil chemical properties and the SOM fractions that most correlated with PC1 were: pH, P, K+, Ca2+, Mg2+, Al3+, SB, CECp, Al saturation, V, TOC, HA-C, FA-C and HS-C. The soil property (H + Al) and SOM fractions HU-C and HS-C/TOC were those showing the stronger association with PC2 (Tab. 5).

The acidity indicators Al3+ and Al saturation correlated negatively with pH, SB (except for K+), V and P in all soil layers (Fig. 4). In areas with a predominance of acidity, soil fertility was limited both by lower levels of exchangeable basic cations and by low soil P content. The CECp correlated positively with the contents of (H + Al), TOC and the humic fractions of SOM. The ability to adsorb cations in these environments was dependent on SOM, showing its influence on soil fertility (Fig. 4).

Fig. 4 - Biplots of the soil chemical properties and humic fractions of the organic matter. (a) topsoil (0-5 cm), (b) subsoil 1 (5-10 cm) and (c) subsoil 2 (10-20 cm). (PNV): preserved native vegetation; (DNV): degraded native vegetation; (CCC): Cassava conventional cultivation; (EAC): Eucalyptus agro-energy cultivation; (pH): hydrogen potential (H2O/1:2.5); (P): available phosphorus; (K): exchangeable potassium; (Ca): exchangeable calcium; (Mg): exchangeable magnesium; (Al): exchangeable aluminum; (H+Al): potential Acidity; (SB): sum of bases; (CECp): cations exchange capacity potential; (V): base saturation; (Al saturation): saturation by aluminum; (TOC): total organic carbon; (HA-C): humic acid fraction; (FA-C): fulvic acid fraction; (HU-C): humin fraction; (HS-C): humic substances fractions; (HA-C/FA-C): humic acid fraction/fulvic acid fraction; (HS-C/TOC): humic substances fractions/total organic carbon. The arrows represent the directions of the high weighting of soil chemical properties and humic fractions of the organic matter along the first (PC1) and second (PC2) principal components. PC1 and PC2 of topsoil explained 96.00% of the total variation, while they accounted for 94.50% and 94.30% of the variation is soil chemical parameters in subsoil 1 (5-10 cm) and subsoil 2 (10-20 cm), respectively.

  Enlarge/Shrink   Download   Full Width  Open in Viewer

The TOC and humic fractions of SOM were positively correlated. However, their correlation is stronger in the superficial layer and becomes weaker in depth, with TOC showing greater association with FA-C in the layer 10-20 cm. Moreover, the vector of the parameter HU-C in Fig. 4is less overlapping with that of TOC, indicating the low humification level of this TOC.

The chemical properties pH, P, Ca2+, Mg2+ and V reached their largest values in the soil of EAC and CCC plots. The PNV areas were associated with TOC and humic fractions of SOM. The DNV area was more influenced by acidic properties Al3+ and Al saturation (Fig. 4). The CCC soil properties approached those of the DNV area in the surface layer, while they are less similar as for deeper layers of the soil (Fig. 4).

The dendrogram of similarity between different land uses based on soil chemical properties (Fig. 5) showed that the DNV area differed from the other areas in all soil layers. The PNV area is similar to the cultivated areas in superficial layers (CCC, 0-5 cm and EAC, 5-10 cm).

Fig. 5 - Dendrograms of the different land uses according to the soil chemical properties and humic fractions of the organic matter. (a) topsoil (0-5 cm); (b) subsoil 1 (5-10 cm); (c) subsoil 2 (10-20 cm). (PNV): preserved native vegetation; (DNV): degraded native vegetation; (CCC): Cassava conventional cultivation; (EAC): Eucalyptus agro-energy cultivation.

  Enlarge/Shrink   Download   Full Width  Open in Viewer

  Discussion 

The soils of areas with different land uses showed comparable texture at the studied depths (Tab. 2). Many physical and chemical properties of the soil are directly affected by texture, such as bulk density, pore size distribution, CEC, SOM, and others. In view of this, the deforestation and land use changes significantly influenced their properties.

Soil physical properties

We found that the increase in BD in DNV was explained by reduction in soil porosity (mainly MaP), agreeing with the findings of Bottinelli et al. ([8]) and Chen et al. ([11]). According to Hebb et al. ([16]), the decrease of MaP, together with the increase of BD, decreases soil water permeability because drainage pores are main responsible for water percolation in the soil. Kuncoro et al. ([22]) found a positive correlation between water permeability and MaP, associating this result with the size and discontinuity of macropores.

The MaP little varied in the subsurface layer, indicating that different land uses weakly affect this physical property. Indeed, soil macroporosity may reflect the contribution of plant roots and plant residues, which occur mainly on the soil surface ([8]). In addition, we found a negative correlation between MaP and BD, which was generally higher in the subsurface layer, causing MaP reduction.

The degradation of soil physical properties in the DNV area can be favored by the reduced vegetation cover and the consequent higher exposure to rainfall ([3]) and animal trampling ([13]), which may cause the compaction of the surface layer. Taylor & Ashcroft ([45]) stated that values of MaP below 0.10 m3 m-3 limit gas exchange, impairing root respiration and development. In the PNV and EAC areas, the action of roots on the soil surface and larger increment of SOM may have promoted greater aggregation and structuring of the soil, as indicated by the higher values of MaP.

The highest values of BD and SRP in DNV were also related with the lower TOC content in this area in the 0-5 and 10-20 cm layers, promoting lower root content and microbial action as SOM improves structuring and porosity ([19]). SRP values greater than 2.0 MPa are above the critical levels reported by Taylor et al. ([44]), and this can hamper plant root growth and favor loss of water and air from soil pores.

The strong association between the MiP, FC and AW soil properties in the CCC area may be explained by the greater homogenization due to tillage applied for crop cultivation, thus reducing pores and increases water retention ([25]). However, the similarity between the EAC area in the subsurface layer and the DNV area may be due to a lower decomposition of plant residues in the subsurface ([51]), which results in low aggregation capacity of the soil in deeper layers and consequently in the degradation of the physical properties of the soil.

Soil chemical properties

The highest values of pH in cultivated areas (CCC and EAC) can be attributed to liming and fertilization practices that increased the number of exchangeable bases and soil pH ([1]). In these areas, also Al+3 complexation may occur due to SOM acid groups ([17]). In CCC area the plant residues are usually incorporated in the soil through tillage or hoeing, while in EAC area are maintained on the soil surface.

The high acidification of the DNV area may be related to the reduction of exchangeable bases due to surface erosion and leaching during periods of higher rainfall, as a consequence of the removal of native vegetation leaving the soil exposed ([27]). However, in areas of preserved native vegetation, soil pH is also usually acid ([6]). According to Silva et al. ([36]) in areas of preserved vegetation, the ecosystem in equilibrium tends to a natural acidification of the soil due to SOM mineralization and the release of acidic exudates by plant root systems.

The high amounts of exchangeable bases in the soil of planted forests such as eucalyptus cultivation may be associated with fertilization and liming practices in these areas and to the contribution of plant material to the soil surface. According to Monroe et al. ([28]), forested areas have the capacity to increase soil fertility by the addition of organic matter to soil from both the above- and belowground biomass.

The amount of exchangeable bases in areas of preserved vegetation are naturally low ([25]). Allocation of these nutrients to plant biomass and its return to the soil through the cycling process adequately nourishes the plants ([23]).

The high P availability in the soil of these dry ecosystems is due to chemical fertilization. However, in areas with preserved vegetation, it can be promoted by the organic material contribution, which reduces the adsorption and precipitation of P ([6]), besides being a source of release ([52]) and recycling of P in organic form ([28]).

Total organic carbon and humic fractions of organic matter

The similarity of TOC contents between PNV and EAC areas in the soil surface layer showed that increase in TOC content is favored by continuous deposition of litter and by the presence of more developed root systems on the surface ([10], [28]), which may contribute to maintain soil moisture and reduce soil erosion ([15]). In this study, the reduction of TOC in the DNV and CCC areas indicated that the soil preparation, which exposes SOM to weathering, and the lack of conservation practices in the CCC area have reduced physical protection of SOM and left the soil exposed to erosion ([27], [49]). These results corroborate with Sharma et al. ([34]) who found lower TOC levels in an environment of conventional agriculture and in degraded areas compared to those where the soil was kept covered.

The increase of HA-C in PNV and EAC areas reveals that the management that maintains soil coverage, adequate humidity and minimal or no anthropic action may promote microbial action and humic fraction production, and more developed root systems as well ([29], [43]). In fact, lower accumulation of HA-C in DNV and CCC areas may be related to the process of degradation of the areas due to soil exposure to erosive processes ([27]) and/or to the low increment of plant biomass.

The higher proportion of HU-C in relation to HA-C and FA-C in the topsoil of different land uses is due to its molecule size which determines a higher chemical affinity to colloids ([27]), making them more stable and resistant to weathering. However, the microbial activity is reduced in semi-arid regions ([39]), resulting in a low content of this fraction of SOM. Levels of HU-C similar to those of this study were found by Sousa et al. ([39]) and Melo et al. ([27]) who investigated the behavior of the different SOM fractions in soils of the Brazilian semi-arid region.

Low biomass production and low organic matter input to soil are common in semi-arid regions ([27]). However, the HA-C/FA-C ratio was > 1 in the PNV and EAC areas due to the greater contribution of organic material, reflecting the high degree of humification in these soils ([10]). Contrastingly, the HA-C/FA-C ratio < 1 in DNV and CCC areas suggests a lower contribution of organic residues and lower microbial activity, which hampers the formation of more condensed humic substances ([27]). The HA-C/FA-C ration indicates the quality of humus of the soil, with smaller ratios expected for soils subjected to more intense degradation processes ([41]).

The humification degree (HS-C/TOC) varying from 72 to 87% in all areas under different land uses indicates that most of TOC is in humic fractions, suggesting a high level of humification.

We found that P, Ca2+ and TOC, SB and Al saturation were the most important chemical properties to explain data variation, in addition to the HA-C fraction of SOM, which strongly discriminate areas under different land uses. Marinho et al. ([24]) and Sousa Neto et al. ([40]), studying different management systems in the semi-arid region of Brazil, found that the different sites could be discriminated based on soil chemical properties such as pH, P, K+, Ca2+, Mg2+, CEC and TOC.

The negative correlation of Al3+ with pH, nutrient availability and humic fractions of SOM suggests a complexation of Al3+ by humic substances ([47]), resulting in an increase of nutrient availability and soil pH. In fact, our results highlight the influence of organic matter on soil fertility as a source of nutrients for the plants, with K+ being the most mineralized nutrient in the study areas. The relationship between CEC and TOC showed that a higher organic content in the soil can be very important in the generation of electric charges, maximizing the soil’s ability to retain cations ([47]).

The close association of FA-C with TOC in the subsurface may be due to characteristics of the FA, which is highly soluble and mobile in the soil, and thus its content incraeses with depth ([33]). Contrastingly, we observed a greater distance between the HU-C and TOC vectors, suggesting a small influence of this SOM fraction on deeper layers of the soil.

We found that the differences in soil chemical properties between DNV and the other land uses become weaker in the subsurface layer, and this may reflect the low capacity of the different land uses to improve chemical quality and SOM at the lower soil depths in these dry habitats. The higher recalcitrance of the eucalyptus plant material ([51]), associated with the water conditions of the semi-arid region ([27]) may explain these results.

The clustering of CCC and EAC areas in soil subsurface indicated that cultivation practices adopted in these areas did not affect the soil chemical properties and SOM fractions. The low capacity of eucalyptus to improve soil conditions at lower depth ([51]) may explain the above similarity, mainly in terms of TOC and SOM fractions. Further, the independent clustering of DNV from the other land uses may be due to the strong degradation of the soil quality resulting from lower vegetation cover, greater erosion and leaching, and lower input of organic material.

SOM preservation and transformation processes in different land uses has the greatest importance and has to be studied in detail. However, it has been widely questioned whether the traditional chemical fractionation involving extraction with strong alkali can provide information on SOM dynamics and turnover in soils ([20]). The alkali extraction is often considered less suitable to obtain information on SOM preservation and turnover than physical fractionation methods, such as density fractionation procedures ([21]), which can distinguish light active fractions, occluded intermediate fractions and mineral-associated, more stable C pools. Nonetheless, the chemical fractionation method used in this study provided good indicators of the SOM preservation and transformation processes.

Finally, the improvement of soil properties in the surface layer by the EAC does not prove the hypothesis tested.

  Conclusions 

The physical and chemical properties that best described the soil quality of different land uses in Chapada do Araripe were MaP, BD, SRP, SB (mainly Ca2+), available P and Al saturation. The TOC and the HA-C fraction of SOM were also important in improving soil quality. The DNV area showed degradation of both chemicals and physical properties, resulting in the reduction of soil quality. Improvement of these properties in the EAC area was limited to the surface layer. Grouping obtained by cluster analysis showed that DNV (degraded native vegetation) and PNV (preserved native vegetation) areas are different from the cultivated areas (EAC and CCC) in terms of soil physical and chemical parameters, though EAC and CCC areas approach the PNV area in the soil surface layers, while DNV showed a strong reduction in soil quality.

Our study suggests that the disordered cutting of native vegetation causes land degradation in Chapada do Araripe, which results in the reduction of soil quality. Crops that allow the input of biomass and keep the soil covered, protecting it from weathering, can maintain or improve soil properties and functions, although this improvement has not been observed below the soil surface (>10 cm depth).

  Author contributions 

JSR: investigation, methodology, formal analysis and writing the original draft; FJF; conceptualization, funding acquisition, supervision, writing review and editing; JCAF: resources, methodology and data curation; MBGSF: resources, methodology and visualization; BGA: formal analysis, data curation and visualization; LRCS: formal analysis and data curation.

  Acknowledgements 

The authors thank the Experimental Station of Araripina of the Agronomic Institute of Pernambuco (IPA), Brazil for the logistic and operational support. The study was supported by Pernambuco State Science and Technology Support Foundation (FACEPE) and by the National Council for Scientific and Technological Development (CNPq), project no. APQ-0729-5.01/14 “Energia de biomassa em Pernambuco: impactos ambientais e alternativas sustentáveis (ProMassa)”.

  References

(1)
Abegaz A, Adugna A (2015). Effects of soil depth on the dynamics of selected soil properties among the highlands resources of Northeast Wollega, Ethiopia: are these sign of degradation? Solid Earth 7: 2011-2035.
CrossRef | Gscholar
(2)
Alvares CA, Stape JL, Sentelhas PC, Gonçalves JLM, Sparovek G (2013). Köppen’s climate classification map for Brazil. Meteorologische Zeitschrift 22: 711-728.
CrossRef | Gscholar
(3)
Araújo Filho RN, Freire MBGS, Wilcox BD, West JB, Freire FJ, Marques FA (2018). Recovery of carbon stocks in deforested Caatinga dry forest soils requires at least 60 years. Forest Ecology and Management 407: 210-220.
CrossRef | Gscholar
(4)
Arruda LEV, Portela JC, Dias NS, Porto VCN, Ferreira Neto M, Silva SYAM, Dias MCC, Oliveira Filho TJ (2015). Oxisolphysical properties under different agricultural uses in Brazil. African Journal Agricultural Research 10: 4331-4337.
CrossRef | Gscholar
(5)
Ashagrie Y, Zech W, Guggenberger G, Mamo T (2007). Soil aggregation, and total and particulate organic matter following conversion of native forests to continuous cultivation in Ethiopia. Soil and Tillage Research 94: 101-108.
CrossRef | Gscholar
(6)
Avila MA, Souza SR, Veloso MDM, Santos RM, Fernandes LA, Nunes YRF (2016). Structure of natural regeneration in relation to soil properties and disturbance in two swamp forests. Cerne 22: 1-10.
CrossRef | Gscholar
(7)
Bone J, Barraclough D, Eggleton P, Head M, Jones DT, Voulvoulis N (2012). Prioritising soil quality assessment through the screening of sites: the use of publicly collected data. Land Degradation and Development 25: 251-266.
CrossRef | Gscholar
(8)
Bottinelli N, Hallaire V, Goutal N, Bonnaud P, Ranger J (2014). Impact of heavy traffic on soil macroporosity of two silty forest soils: initial effect and short-term recovery. Geoderma 217- 218: 10-17.
CrossRef | Gscholar
(9)
Bradford JM (1986). Penetrability. In: “Methods of soil analysis: Physical and mineralogical methods” (Klute A eds). SSSA Book Series, American Society of Agronomy, Madison, WI, USA, pp. 463-478.
CrossRef | Gscholar
(10)
Cardoso JAF, Lima AMN, Cunha TJF, Rodrigues MS, Hernani LC, Amaral AJ, Oliveira Neto MBO (2015). Organic matter fractions in a quartzipsamment under cultivation of irrigated mango in the lower São Francisco valley region, Brazil. Revista Brasileira de Ciência do Solo 39: 1068-1078.
CrossRef | Gscholar
(11)
Chen G, Weil RR, Hill RL (2014). Effects of compaction and cover crops on soil least limiting water range and air permeability. Soil and Tillage Research 136: 61-69.
CrossRef | Gscholar
(12)
Dane JH, Topp CG (2002). Methods of soil analysis: Part 4. Physical methods. SSSA Book Series, Soil Science Society of American, Madison, WI, USA, pp. 1663.
Gscholar
(13)
Evans CRW, Krzic M, Broersma K, Thompson DJ (2012). Long-term grazing effects on grassland soil properties in southern British Columbia. Canadian Journal of Soil Science 92: 685-693.
CrossRef | Gscholar
(14)
Fontenele MN, Silva JAA, Ferreira RLC, Berger R, Gadelha FHL, Guera MOG (2018). Volumetric and economic evaluation of Eucalyptus spp. clones in high forest, coppice and native vegetation systems at the Gypsum Pole of Araripe - PE. Forest Sciences 46: 77-86.
CrossRef | Gscholar
(15)
Fultz LM, Moore-Kucera J, Zobeck TM, Acosta-Martínez V, Allen VG (2013). Aggregate-carbon pools after 13 years of integrated crop-livestock management in semiarid soils. Soil Science Society of America Journal 77: 1659-1666.
CrossRef | Gscholar
(16)
Hebb C, Schoderbek D, Hernandez-Ramirez G, Hewins D, Carlyle CN, Bork E (2017). Soil physical quality varies among contrasting land uses in Northern Prairie regions. Agriculture Ecosystems Environment 240: 14-23.
CrossRef | Gscholar
(17)
Inda AV, Torrent J, Barron V, Bayer C (2010). Aluminum hydroxy-interlayered minerals and chemical properties of a subtropical Brazilian Oxisol under no-tillage and conventional tillage. Revista Brasileira de Ciência do Solo 34: 33-41.
CrossRef | Gscholar
(18)
James DW, Wells KL (1990). Soil sample collection and handling: technique based on source and degree of field variability. In: “Soil testing and plant analysis” (Westerman RL eds). SSSA Book Series vol. 2, Soil Science Society of America, Madison, WI, USA, pp. 25-44.
CrossRef | Gscholar
(19)
Kassa H, Dondeynec S, Poesenc J, Franklb A, Nyssen J (2017). Impact of deforestation on soil fertility, soil carbon and nitrogen stocks: the case of the Gacheb catchment in the White Nile Basin, Ethiopia. Agriculture Ecosystems Environment 247: 273-282.
CrossRef | Gscholar
(20)
Kleber M, Lehmann J (2019). Humic substances extracted by alkali are invalid proxies for the dynamics and functions of organic matter in terrestrial and aquatic ecosystems. Journal of Environmental Quality 48: 207-216.
CrossRef | Gscholar
(21)
Kögel-Knabner I, Guggenberger G, Kleber M, Kandeler E, Kalbitz K, Scheu S, Eusterhues K, Leinweber P (2008). Organo-mineral associations in temperate soils: Integrating biology, mineralogy, and organic matter chemistry. Journal of Plant Nutrition and Soil Science 171: 61-82.
CrossRef | Gscholar
(22)
Kuncoro PH, Koga K, Satta N, Muto Y (2014). A study on the effect of compaction on transport properties of soil gas and water. II: Soil pore structure indices. Soil and Tillage Research 143: 180-187.
CrossRef | Gscholar
(23)
Leite FP, Silva IR, Novais RF, Barros NF, Neves JCL (2010). Alterations of soil chemical properties by eucalyptus cultivation in five regions in the Rio Doce Valley. Revista Brasileira de Ciência do Solo 34: 821-831.
CrossRef | Gscholar
(24)
Marinho ACCS, Portela JC, Silva EF, Dias NS, Sousa Júnior FS, Silva AC, Silva JF (2016). Organic matter and physicochemical properties of a Cambisol under different agricultural uses in a semi-arid region of Brazil. Australian Journal of Crop Science 10: 32-41. -
Online | Gscholar
(25)
Martins AF, Salcedo IG, Oliveira FP, Pereira WE (2018). Physical, chemical, and microbiological properties of soil under different plant covers in the Seridó desertification region in the Brazilian semiarid. Revista Brasileira de Ciência do Solo 43: 1-12.
CrossRef | Gscholar
(26)
Mekuria W, Veldkamp E, Haile M, Nyssen J, Muys B, Gebrehiwot K (2007). Effectiveness of exclosures to restore degraded soils as a result of overgrazing in Tigray, Ethiopia. Journal of Arid Environments 69: 270-284.
CrossRef | Gscholar
(27)
Melo AAB, Valladares GS, Ceddia MB, Pereira MG, Soares I (2016). Spatial distribution of organic carbon and humic substances in irrigated soils under different management systems in a semi-arid zone in Ceará, Brazil. Semina: Ciencias Agrarias 37: 1845-1856.
CrossRef | Gscholar
(28)
Monroe PHM, Gama-Rodrigues EF, Gama-Rodrigues AC, Marques JRB (2016). Soil carbon stocks and origin under different cacao agroforestry systems in Southern Bahia, Brazil. Agriculture Ecosystems Environment 221: 99-108.
CrossRef | Gscholar
(29)
Ramos DD, Silva EF, Ensinas SC, Souza NH, Potrich DC, Freitas ME, Formagio ASN, Vieira MC (2013). Stocks of carbon, total nitrogen and humicsubstances in soil under different cropping systems. Semina: Ciencias Agrarias 34: 2219-2228.
CrossRef | Gscholar
(30)
Raulino WNC, Freire FJ, Assunção EAA, Ataí De KMP, Silva HV, Silva ACF (2020). Nutrition of tree species in tropical dry forest and rainforest environments. Revista Ceres 67: 70-80.
CrossRef | Gscholar
(31)
Santos HG, Jacomine PKT, Anjos LHC, Oliveira VA, Lumbreras JF, Coelho MR, Almeida JA, Araújo Filho JC, Oliveira JB, Cunha TJF (2018). Sistema brasileiro de classificação de solos [Brazilian system of soil classification]. EMBRAPA, Brasília, Brazil, pp. 356. [in Portuguese]
Gscholar
(32)
Santos RL, Freire FJ, Rocha AT, Silva JAA, Tavares JA, Ferreira EGBS, Oliveira ECA (2015). Elephant grass (Pennisetum purpureum Schum.) biomass production as promising alternative source of energy in Brazil’s semiarid area using gypsum. Australian Journal of Crop Science 9: 1082-1088. -
Online | Gscholar
(33)
Schnitzer M (1986). Binding of humic substances by soil mineral colloids. In: “Interactions of soil minerals with natural organics and microbes”. (Huang PM, Schnitzel M eds). SSSA Special Publications vol. 23, Soil Science Society of American, Madison, WI, USA, pp. 83-87.
CrossRef | Gscholar
(34)
Sharma V, Hussain S, Sharma KR, Aryam VM (2014). Labile carbon pools and soil organic carbon stocks in the foothill Himalayas under different land use systems. Geoderma 232- 234: 81-87.
CrossRef | Gscholar
(35)
Shiker MAK (2012). Multivariate statistical analysis. British Journal of Science 6: 55-66.
Online | Gscholar
(36)
Silva GF, Santos D, Silva AP, Souza JM (2015). Soil quality indicators under different land use systems in the Agreste region of Paraiba, Brazil. Revista Caatinga 28: 25-35.
CrossRef | Gscholar
(37)
Snedecor GW, Cochran WG (1967). Statistical methods. Iowa State University Press, Ames, IO, USA, pp. 593.
Gscholar
(38)
Soil Survey Staff (2014). Keys to soil taxonomy. Natural Resources Conservation Service, USDA, Washington, DC, USA, pp. 372.
Gscholar
(39)
Sousa FP, Ferreira TO, Mendonça ES, Romero RE, Oliveira JGB (2012). Carbon and nitrogen in degraded Brazilian semi-arid soils undergoing desertification. Agriculture Ecosystems Environment 148: 11-21.
CrossRef | Gscholar
(40)
Sousa Neto ON, Dias NS, Lira RB, Silva EF, Ferreira ALL, Freitas JJR (2017). Chemical properties of traditional agriculture and Caatinga managed at different depths in an Inceptisol. Revista Brasileira de Engenharia Agrícola e Ambiental 21: 50-55.
CrossRef | Gscholar
(41)
Sousa RF, Brasil EPF, Figueiredo CC, Leandro WM (2015). Soil organic matter fractions in preserved and disturbed wetlands of the Cerrado Biome. Revista Brasileira de Ciência do Solo 39: 222-231.
CrossRef | Gscholar
(42)
Swift RS (1996). Organic matter characterization. In: “Methods of soil analysis. Part 3. Chemical methods” (Sparks DL eds). SSSA Book Series vol. 23, Soil Science Society of American, Madison, WI, USA, pp. 1018-1020.
CrossRef | Gscholar
(43)
Tavares RLM, Nahas E (2014). Humic fractions of forest, pasture and maize crop soils resulting from microbial activity. Brazilian Journal of Microbiology 45: 963-969.
CrossRef | Gscholar
(44)
Taylor HM, Roberson GM, Parker JR (1966). Soil strength-root penetration relations to medium to coarse-textured soil materials. Soil Science 102: 18-22.
CrossRef | Gscholar
(45)
Taylor SA, Ashcroft GL (1972). Physical edaphology: the physics of irrigated and nonirrigated soils. W.H. Freeman, San Francisco, CA, USA, pp. 533.
Online | Gscholar
(46)
Teixeira PC, Donagemma GK, Fontana A, Teixeira WG (2017). Manual de métodos de análise de solo [Manual of soil analysis methods]. EMBRAPA, Brasília, Brazil, pp. 573. [in Portuguese]
Gscholar
(47)
Valladares GS, Pereira MG, Benites VM, Anjos LHC, Ebeling AG, Guareschi RO (2016). Carbon and nitrogen stocks and humic fractions in Brazilian organosols. Revista Brasileira de Ciência do Solo 40: 1-16.
CrossRef | Gscholar
(48)
Vieira RMSP, Tomasella J, Alvalá RCS, Sestini MF (2015). Identifying areas susceptible to desertification in the Brazilian northeast. Solid Earth 6: 347-360.
CrossRef | Gscholar
(49)
Villarino SH, Studdert GA, Baldassini P, Cendoya MG, Ciuffoli L, Mastrángelo M, Piñeiro G (2017). Deforestation impacts on soil organic carbon stocks in the Semiarid Chaco Region, Argentina. Science of The Total Environment 575: 1056-1065.
CrossRef | Gscholar
(50)
Yeomans JC, Bremner JM (1988). A rapid and precise method for routine determination of organic carbon in soil. Communications in Soil Science and Plant Analysis 19: 1467-1476.
CrossRef | Gscholar
(51)
Zancada MC, Almendros G, Ballesta RJ (2003). Humus quality after eucalypt reforestations in Astúrias (Northern Spain). Science of The Total Environment 313: 245-258.
CrossRef | Gscholar
(52)
Zhang P, Chen X, Wei T, Yang Z, Jia Z, Yang B, Han Q, Ren X (2016). Effects of straw incorporation on the soil nutrient contents, enzyme activities, and crop yield in a semiarid region of China. Soil and Tillage Research 160: 65-72.
CrossRef | Gscholar

Authors’ Affiliation

(1)
Jefrejan Souza Rezende 0000-0002-2606-9386
Universidade Estadual do Piauí, Campus Prof. Barros Araújo, Rodovia BR-316, KM 299, Altamira 64602000 Picos, PI (Brazil)
(2)
Fernando José Freire 0000-0002-3264-712X
Maria Betânia Galvão Dos Santos Freire 0000-0002-0872-5909
Brivaldo Gomes de Almeida 0000-0003-3260-0767
Leandro Reis Costa Santos
Universidade Federal Rural de Pernambuco, Rua Dom Manoel de Medeiros, s/n, Dois Irmãos, 52171-900 Recife, PE (Brazil)
(3)
José Coelho de Araújo Filho 0000-0002-8318-7418
Empresa Brasileira de Pesquisa Agropecuária (EMBRAPA), Rua Antônio Falcão, 402, Boa Viagem 51020-240 Recife, PE (Brazil)

Corresponding author

 
Jefrejan Souza Rezende
jefrejansouza@pcs.uespi.br

Citation

Souza Rezende J, Freire FJ, Araújo Filho JC, Dos Santos Freire MBG, Gomes de Almeida B, Costa Santos LR (2022). Impact of deforestation on the soil physical and chemical attributes, and humic fraction of organic matter in dry environments in Brazil. iForest 15: 465-475. - doi: 10.3832/ifor4016-015

Academic Editor

Giorgio Alberti

Paper history

Received: Nov 10, 2021
Accepted: Sep 12, 2022

First online: Nov 18, 2022
Publication Date: Dec 31, 2022
Publication Time: 2.23 months

© SISEF - The Italian Society of Silviculture and Forest Ecology 2022

  Open Access

This article is distributed under the terms of the Creative Commons Attribution-Non Commercial 4.0 International (https://creativecommons.org/licenses/by-nc/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Creative Commons Licence

Breakdown by View Type

(Waiting for server response...)

Article Usage

Total Article Views: 403
(from publication date up to now)

Breakdown by View Type
HTML Page Views: 0
Abstract Page Views: 0
PDF Downloads: 333
Citation/Reference Downloads: 0
XML Downloads: 70

Web Metrics
Days since publication: 496
Overall contacts: 403
Avg. contacts per week: 5.69

Article citations are based on data periodically collected from the Clarivate Web of Science web site
(last update: Nov 2020)

(No citations were found up to date. Please come back later)


 

Publication Metrics

by Dimensions ©

List of the papers citing this article based on CrossRef Cited-by.

 

iForest Similar Articles

 

This website uses cookies to ensure you get the best experience on our website. More info