Impacts of Extensive Sheep Grazing on Soil Physical and Chemical Quality in Open Mountain Forests, NE Portugal

Grazing and mechanical clearing are common techniques for vegetation management in open Mediterranean forests. Despite its recognized bene ﬁ ts in the prevention of high-intensity and severity forest ﬁ res, it is essential to consider its impacts on the physical and chemical soil properties. In an open mountain forest located in the NE of Portugal, soil samples were analyzed at depths 0 – 5, 5 – 10, and 10 – 20 cm collected at two moments: before mechanical clearing, (Control) and after 18 months of extensive sheep grazing, in areas without grazing, only mechanical clearing (MC) and in areas with both mechanical clearing and grazing (MCG). The results indicate that vegetation cutting has induced a signi ﬁ cant decrease in extractable potassium, and an increase in the soil organic matter and total nitrogen. The exchangeable bases and the exchangeable acidity did not undergo expressive changes, as indicated by the pH values and the cation exchangeable capacity. After grazing, extractable phosphorus and potassium, organic matter, total nitrogen, exchangeable bases, and cation exchangeable capacity have increased signi ﬁ cantly in the topsoil (0 – 5 cm), reducing soil acidity. Regarding physical properties, only soil permeability has been negatively affected by grazing. Mediterranean mountain open forests management with the combination of vegetation clearing and extensive sheep grazing proved to be effective in reducing vegetable fuel availability and improving soil quality.


INTRODUCTION
Grazing effects on soil quality and ecosystems services across the world are often controversial (Lai and Kumar, 2020;Fielding, 2022).However, extensive pastoral systems are an effective tool for vegetation management in open mountain forests in the Mediterranean region.It contributes to vegetation maintenance (herbaceous and shrubby) and reduces hazardous fuels, rising the prevention of wildfire occurrence (Pardini et al., 2007;Fonseca et al., 2021;Castro et al., 2022).Pastoralism should be used as a complementary technique to the prescribed fire or mechanical clearing (Torres-Manso, 2011;Castro et al., 2022), which can reduce the costs associated with these traditional practices by around 75%, following the increase in intervals of time between interventions (12-year intervals with grazing versus 3-year intervals without grazing) (Varela et al., 2007).In addition, environmental gains must be added, as the recurrent use of prescribed fire causes substantial damage to ecosystems, such as the degradation of soil and water quality and greenhouse gases emission (Fonseca et al., 2017).Grazing by herbivores can be a climate mitigation strategy as it influences the quantity and stability of the soil carbon pool (Naidu et al., 2022).
In Portugal, herding activity has been practiced mainly in mountain areas, where soil resource is scarce (Figueiredo, 2013) and dominant land use is composed of forests, bushes, and sparse pastures.In this sense, extensive grazing with species and breeds highly adapted to these constraints is, in many rural zones, one of the few economic activities capable of guaranteeing the local inhabitants' subsistence (Pinto et al., 2023).The Portuguese historical past, where extensive livestock played a crucial role in the economy of rural areas, is changing due to the scarcity and demanding workforce (shepherds and breeders), rural exodus, population aging, and lack of rural spaces organization (Fonseca et al., 2021).Effectively, in the last 100 years, there has been a sharp reduction in herding activity and in the number of shepherds and raisers (INE, 2021).The National Strategy for the Portuguese Forest indicates extensive grazing as a fuel management technique, but its use for this purpose is still very limited (Castro et al., 2018), since the most important areas for the practice of livestock are located in places with low population density.
Adequate animal load and a short residence time in a given space are essential factors to avoid excessive vegetation disappearance and consequent exposure of the soil to erosion processes and deterioration of the quality of habitats for wildlife (Díaz-Solís et al., 2016;Fielding, 2022).A well-managed pasture in space and time can provide wider biodiversity, nutrients distribution in the soil, and soil quality improvement.Soil quality is a concept centered on the dynamics of soil properties and processes (Karlen et al., 1994;Fonseca and Figueiredo, 2016), constituting an essential factor in the global sustainability of the biosphere and, in particular, in the sustainability of forest systems (Fonseca et al., 2019).In this sense, despite the recognized benefits of extensive grazing in reducing fuel and preventing forest fires, it is important to know the impacts caused on physical and chemical soil properties.Soil is a strategic resource for the sustainable development of humanity, and its preservation is of unquestionable importance.The present work aims to study the effects of mechanical clearing and extensive sheep grazing in an open Mediterranean mountain forest, on soil physical and chemical quality.Soils were sampled before the application of mechanical clearing (Control, C) and after 18 months the effects of grazing were studied.Grazed areas (mechanical clearing plus grazing, MCG) were compared with non-grazing areas (mechanical clearing, MC) and original soil (C).

MATERIALS AND METHODS
The study area, with forestry aptitude, was located in Romeu, Mirandela municipality, NE Portugal (41 °32′25″N 7 °02′15″W), in the range between 474 and 534 m altitude (Figure 1).The annual average temperature is 14 °C and annual average precipitation is 508 mm, concentrated from October to April and with a marked dry season among June and August (Agroconsultores and Coba, 1991).Soils are classified as Distric Leptosols derived from schist (FAO, 2015).It is an open mountain forest where the arboreal vegetation is essentially composed of Quercus rotundifolia, Quercus faginea and Quercus suber.Associated with these species there is a relatively dense understory of Cytisus multiflorus, Cytisus scoparius, Cystus ladanifer and Lavandula stoechas, and some herbaceous vegetation species.
The experimental area occupies 4,300 m 2 and comprises three treatments: 1) plot without mechanical clearing with extensive grazing, which represents the traditional treatment, with an area of 906 m 2 (Control, C); 2) plot of 400 m 2 with mechanical clearing without grazing (MC); 3) plot of 41,694 m 2 with mechanical clearing plus grazing (MCG) (Figure 2).Extensive grazing by a flock of 150 sheep took place over an 18-month period, from December 2019 to May 2021.
In 2019, the soil was characterized before the application of mechanical clearing (C).Disturbed soil samples were collected from four randomly distributed sites at depths 0-5, 5-10, and 10-20 cm.Also, in the surface soil layer (0-5 cm) undisturbed soil samples with 100 cm 3 of volume were taken to determine physical soil properties.In 2021, in the area subjected to mechanical clearing, the effects of grazing on soil physical and chemical properties were evaluated.Disturbed soil samples were collected at depths 0-5, 5-10, and 10-20 cm in four sites randomly distributed in the area preferentially grazed by sheep (MCG) (Figure 3), and in the non-pastured area (MC).Permeability was measured at the same points in undisturbed soil samples taken in 100 cm 3 cylinders from the topsoil (0-5 cm).In the traditional treatment (C) no evaluations were carried out in 2021, as the flock of sheep remained in this area for a very short time, and significant changes in soil properties were not foreseeable in the period considered (18 months).
Disturbed soil samples were air-dried and sieved with a 2 mm mesh sieve, in order to assess the organic matter, nutrients concentration, soil pH, exchangeable acidity, cation exchangeable capacity, and soil particle size (soil texture).Methods applied in laboratory analysis comprised the Tinsley method for soil organic matter (Tinsley, 1950), Kjeldahl method for Total N (Bremner, 1996), the Egner-Riehm method for extractable phosphorus and potassium (Almeida and Balbino, 1960), and hydrometer method for soil texture determination.The pH was determined in a soilwater suspension (1: 2.5 soil water ratio).Exchangeable bases were evaluated by atomic adsorption (Ca 2+ and Mg 2+ ) and by flame emission spectrophotometry (K + and Na + ).Cation exchangeable capacity was calculated as the sum of exchangeable bases and exchangeable acidity.Soil permeability was measured, in undisturbed soil samples, in a constant head closed circuit laboratory permeameter.In the same undisturbed soil samples, other physical properties were also evaluated, which are related to permeability, such as bulk density, total porosity, microporosity, macroporosity, field capacity, and water saturation capacity.
The soil classification in the study area was carried out based on the description and observation of five soil profiles and the analysis of samples collected in the horizons identified in the profiles.A representative soil profile was selected, based on the greater differentiation of horizons and their clearest delimitation.
The statistical analysis of the information collected included descriptive statistics, analysis of variance (one-way), and tests of means comparison (Tukey, p < 0.05).

Soil Description and Classification
The morphological properties (Figure 4) and the physical and chemical soil properties evaluated (Table 1) show that the soils are very shallow (<30 cm), with high stoniness, sandy loam texture, low pH (<5.5), medium to low in organic matter content, low total N, low exchangeable bases and cation exchangeable capacity, low P and moderate to high K contents.There will be no problems associated with aluminum toxicity, as the acid saturation is less than 20% (Porta et al., 2003).However, along the slope in very restricted areas, the soil thickness may slightly exceed 30 cm in depth.According to FAO (2015) and soil properties evaluated, the soils are classified as Distric Leptosols derived from schist.
Leptosols comprise very shallow soils over continuous rock to a depth of ≤30 cm or with a high proportion of unconsolidated stony material with <20% (by volume) fine earth up to 75 cm in depth.Thus, Leptosols can be either very shallow soils directly lying on bedrock, or soils with excess rock fragments in a profile.No calcic, chernic, duric, gypsic, petrocalcic, petrogypsic, petroplinthic or spodic horizon (Agroconsultores and Coba, 1991;FAO, 2015).Leptosols are found in all climate zones (many of them in hot or cold dry regions), in particular in strongly eroding areas.Leptosols have a resource potential for wet-season grazing and as forest land.In NE Portugal, Leptosols constitute the dominant soil unit, occupying about 70% of the territory (Agroconsultores and Coba, 1991;Figueiredo, 2013).

Soil Organic Matter and Nitrogen
Extensive grazing (MCG) contributed to a significant increase in soil organic matter (SOM) concentration in the soil surface layer (0-5 cm).The organic matter content almost doubled in relation to the traditional treatment (C) and registered an increase of about 30% compared to the treatment with mechanical clearing (MC), presenting in the remaining layers (5-10 and 10-20 cm), values quite similar among treatments (Figure 5A).As the soil samples were collected in the preferred grazing sites by the sheep, the increase in organic matter in the topsoil may be due to animal manure (feces and urine) and also to the vegetal components humification that was partially incorporated into the surface layer during mechanical clearing (Herrero et al., 2016;Fonseca and Figueiredo, 2018;Fielding, 2022).Similarly, Conant et al. (2003) and Wang et al. (2015) found only SOM increases in the topsoil layer.This increase in carbon also corroborates Abdalla et al. (2018) in a study where different grazing intensities were compared, showing increments of soil carbon under lowintensity grazing.The translocation of organic components in the soil may explain the slight increase in organic matter content in the 5-10 cm layer (Patrício et al., 2018;Fonseca et al., 2022).
The non-significant changes observed in SOM between the areas with mechanical vegetation cutting with grazing (MCG) and without grazing (MC) agree with the low soil disturbance by the mechanized operation to cut the vegetation and its actual low impact.However, even slight soil disturbances might increase the SOM mineralization, which explains the negligible decrease in the SOM content in the surface soil layer (0-5 cm) (Ma et al., 2015;Fonseca and Figueiredo, 2018;Fonseca et al., 2019).
SOM plays a fundamental role in sustaining soil quality (Percival et al., 2000), and is often referred to as one of the main indicators of soil quality (Karlen et al., 1994;Sharma et al., 2005;Lal, 2015;Fonseca and Figueiredo, 2016).The increase in SOM due to grazing may represent an important contribution to soil fertility, water retention, aggregation stability, and carbon sequestration (C = 58% SOM), contributing to soil quality improvement and climate change mitigation (Piñeiro et al., 2010;Lal, 2015;Fonseca et al., 2019).Data obtained in this work are in line with the climate-positive strategy, as this grazing model increases the stocks of the soil organic carbon pool (Lai and Kumar, 2020).In the extensive grazing, the animals spend a short time in the same area, moving and leaving behind trampled vegetation covered by manure and urine, which turn out to be incorporated into the soil, improving physical and chemical soil properties (Abdalla et al., 2018;Fielding, 2022).
Nitrogen is found in the soil essentially in organic form (about 98%).In this sense, it presents a strict relationship with the organic matter content (N = 5% SOM), tending the carbon and nitrogen fluxes to be connected within the ecosystems (Piñeiro et al., 2010).Total nitrogen values are very low and present statistical differences (p < 0.05) only in the surface layer of the MC and MCG treatments in relation to the traditional treatment (C) (Figure 5B).Compared to mechanical clearing (MC), nitrogen concentration was not significantly affected by grazing (MCG), and the highest values were recorded on topsoil (0-5 cm).As expected, variations across the soil profile closely follow variations in organic matter.
Likewise, Lai and Kumar (2020), based on a meta-analysis of livestock grazing impacts on soil properties, resultant from the   (Ellsworth, 2004).

Extractable Phosphorus and Potassium
Extractable phosphorus and potassium showed different behavior by grazing effect.In grazing treatment (MCG), phosphorus increased significantly in the topsoil (0-5 cm), maintaining similar values in the deeper layers (5-10, 10-20 cm), but with a tendency to decrease in relation to mechanical clearing (MC) (Figure 6A).The effects of mechanical clearing without grazing (MC) are not very expressive compared to the other treatments (C and MCG).Values are considered low (26-50 mg kg −1 ) to medium (51-100 mg kg −1 ) (Santos, 2015).
With regard to extractable potassium, grazing (MCG) contributed to significant increases in all soil layers (0-5, 5-10, and 10-20 cm).The remaining treatments (MC and C) maintained similar values, with the exception of the topsoil layer (0-5 cm), where the vegetation cutting (MC) showed significantly lower values.Extensive grazing contributed to the values go by from medium (51-100 mg kg −1 ) to very high (>200 mg kg −1 ) (Figure 6B).
Phosphorus low values may be related to retention and fixation phenomena.Portuguese soils, in general, are very poor in this nutrient (Arrobas and Coutinho, 2002).Potassium has high mobility in the soil, undergoing redistribution along the soil profile (Fonseca et al., 2017).These parameter variations in the soil seem to be directly dependent on the effect of the manure (feces and urine) produced by the animals.

Exchangeable Bases, Exchangeable Acidity, Cation Exchangeable Capacity and pH
Table 2 presents the data concerning soil exchangeable complex (exchangeable bases, sum exchangeable bases, exchangeable acidity, cation exchangeable capacity) and pH values.As a result of grazing (MCG), the sum of exchangeable bases (SEB) almost doubled in the surface layer (0-5 cm), recording much lower values in the deeper layers (5-10 and 10-20 cm), but always higher in the treatment with grazing (MCG) when compared to the treatment without grazing (MC).Mechanical clearing (MC) contributed to the loss of exchangeable bases at all soil depths.This is a result of the basic cations having high mobility in the soil and the losses are associated with leaching and erosion processes (Fonseca et al., 2017), promoted by the greater exposure of the soil as a consequence of the vegetation removal.On the opposite, the exchangeable acidity (EA) decreased in the three soil layers owing to the grazing effect (MCG), leading to a reduction of the soil acidity.In general, at all soil layers and treatments, the values of base cations follow the sequence Ca 2+ > Mg 2+ > K + > Na + , a common pattern observed in the soils.Calcium is the most important base cation, contributing more than 50% of the sum of exchangeable bases in the upper soil layers (0-5 and 5-10 cm).
Grazing produced important changes in the sum of exchangeable bases (SEB) and in the exchangeable acidity (EA), which are reflected in the cation exchangeable capacity (CEC) and in the soil pH values.The great increase in CEC in the surface layer (0-5 cm) is related to the corresponding increase in organic matter.Soil pH values were affected by the increase in base cations mainly in the most surface layer, and by the decrease in exchangeable acidity (EA).Due to the effect of grazing, in the topsoil (0-5 cm), soil changed from acid (5.3) to sub-acid (5.8) (Costa, 2011).
Among other factors, soil permeability depends on soil porosity and bulk density (Leite et al., 2013;Fonseca et al., 2017).However, Table 3 shows that grazing did not significantly affect the soil physical properties evaluated in the topsoil (0-5 cm).In this sense, the reduction in permeability may be related to the production of hydrophobic substances on the soil surface (Keizer et al., 2008;Loss et al., 2015).On the other hand, since organic matter represents the key factor in soil quality, its significant increase in the surface layer (Figure 5A) may have contributed to nullifying the effect of trampling on physical soil properties.Also, Shrestha et al. (2020) did not find significant differences in soil bulk density, soil moisture, and soil compaction in extensive grazing systems of small ruminants.
In general, based on the soil physical and chemical properties evaluated, the presence of sheep flocks in the field contributed to increasing soil fertility and carbon sequestration.Identical results were obtained by Silva et al. (2020), in extensive grazing with cattle.Over a grazing period, the degree of soil physical degradation changes with factors such as animal load, hoof pressure, grazing duration, and the previous grazing history (Donovan and Monaghan, 2021).However, the existing literature and examples of light grazing effects on soil properties are limited and not conclusive.

CONCLUSION
The main effects of mechanical clearing without grazing (MC) and mechanical clearing with grazing (MCG) were recorded in the physical and chemical properties of the topsoil layer (0-5 cm), keeping the deeper layers (5-10 and 10-20 cm), in most cases, identical values to the original soil.Mechanical clearing has contributed to a significant increase in soil organic matter and total nitrogen.On the opposite, extractable potassium decreased significantly.Exchangeable bases and the exchangeable acidity did not record expressive changes, as indicated by the pH values and the   Extensive Grazing and Soil Quality cation exchangeable capacity.The soil physical properties did not change due to the cutting vegetation effect.Extensive sheep grazing contributed overall to improved soil fertility.This improvement in fertility was mainly visible in the surface layer (0-5 cm) and resulted in increases in organic matter, total nitrogen, extractable phosphorus, and potassium, the sum of exchangeable bases, cation exchangeable capacity, and soil pH values.In general, grazing did not significantly affect the physical soil properties, only a significant reduction in soil permeability was observed.
Based on the results obtained, in this kind of mountain environment, grazing systems for small ruminants, in conjunction with vegetation clearing, can play an important role in reducing fuel availability (prevention of large forest fires) and in improving productive (increased fertility) and environmental functions of the soil (increased carbon sequestration).
Copyright © 2023 Fonseca, Castro, Alves, Castro and de Figueiredo.This is an openaccess article distributed under the terms of the Creative Commons Attribution License (CC BY).The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice.No use, distribution or reproduction is permitted which does not comply with these terms.

FIGURE 2 |
FIGURE 2 | Experimental design and spatial distribution of the sampling points and soil profiles in the study area.

FIGURE 4 |
FIGURE 4 | Representative profile of the original soil in the study area.

TABLE 1 |
Physical and chemical soil properties of the original soil in the study area.