Evaluation of the E ﬀ ect of Dolomitic Lime Application on Pastures—Case Study in the Montado Mediterranean Ecosystem

: The Montado ecosystem, predominant in the Mediterranean region, consists of poor soils, a sparse cover of cork and holm with an understory of natural biodiverse pastures, grazed by animals in extensive regime. The recommended procedure for increasing productivity of these pastures is based on the application of phosphate fertilizer. One of the main productivity-limiting factors is, however, associated with soil acidity. The objective of this work was to evaluate the simultaneous e ﬀ ect of the holm oak canopy and the application of dolomitic lime on the productivity and quality of a permanent biodiverse pasture, grazed by sheep, in an acid soil (pH = 5.4 ± 0.3). Pasture was monitored at the end of autumn 2018 and winter and spring 2019. The results show that amendment of soil acidity is a slow and gradual process that improves soil Mg / Mn ratio and has a positive impact on pasture productivity and quality. Pasture crude protein availability (CP, kg · ha − 1 ), which is based on both pasture dry matter yield (kg · ha − 1 ) and quality (CP, %), proved to be a very practical indicator of the contributions of tree canopy and soil acidity correction to the holistic management of the Montado ecosystem.


Introduction
The silvo-pastoral ecosystem, predominant in the Alentejo region, is based on Cambisols whose genesis derives from granitic bedrock [1]. As a result of their origin, these soils have characteristics that make them not very fertile, due to their low organic matter content, the very high acidity, or even the reduced cation exchange capacity [2]. In the case of grazed soils, the tendency towards acidification is accentuated, not only by the accumulation of manure, but also by the extraction of nutrients and the leaching of nitrates [3]. On this soil substrate, natural biodiverse pastures develop under a canopy of cork and holm trees, while grazed by animals in extensive regime. The recommended procedure for the recovery of these soils is the installation of permanent pastures and increasing soil fertility through chemical fertilizer applications [4]. One of the main limiting factors of greater productivity is, however, associated with soil acidity, which raises the obvious question: Does the effectiveness of fertilizer application increase if, in advance, the soil pH is corrected?
In addition to soils of low fertility, the Mediterranean climate of these regions is characterized by great irregularity and seasonality, which imposes hot and dry summers and a high concentration of precipitation events in a relatively short winter [5]. Typically, dryland pastures in the Mediterranean climate are composed of annual plants, whose cycle begins with the first rains, in early autumn. Their growth takes place during the autumn, winter, and spring period [4]. After the reproductive phase, at the end of spring, the vegetative cycle ends, with the production of flowers and seeds; this allows these plants to survive the demands of the summer period [4]. The development and productivity of dryland pastures depend mainly on the distribution of precipitation throughout the year and on its combination with the air temperature. Additionally, the fall in temperatures in winter mainly affects legumes, which are less resistant to cold [4].
In addition to soil and climate, the effect of trees is also relevant to pasture growth and results from the deposition of biomass on the soil and the ensuing biological activity. This process has a positive impact on fertility and structure of the soil under the tree canopy, with an increase in organic matter, greater availability of nutrients, greater water retention capacity, and a slightly less acidic pH [2]. Pasture productivity under tree canopy (UTC) can be negatively affected as a result of competition for light, moisture, and nutrients [6]; however, the quality of pasture tends to improve, although this trend varies according to the level of soil fertility [7]. The improvement in crude protein (CP) levels UTC was justified by Pullanagari et al. [8] as a result of a reduction in the photosynthetic rate, and consequently a delay in the physiological development of the pasture UTC, as well as variations in the floristic composition [3]. The combination of these factors results in the plants remaining physiologically younger in the understory and, therefore, with higher metabolic levels for a longer period of time [6].
One of the most important decisions in the Montado management, prompted by the existence of usually poor and acidic soils, is related to fertilization and soil amendment [9]. The basic strategy for soil conservation, reversal of the degradation of this ecosystem, improvement of pastures productivity, and, consequently, of animal production begins with the correction of soil acidity [4]. Carvalho et al. [2] demonstrated the importance of ionic antagonism between soil concentrations of magnesium (Mg) and manganese (Mn) in crop productivity. They called attention to the interest of applying dolomitic lime to correct soil acidity since it simultaneously improves the Mg/Mn ratio in the soil solution. This is followed by the correction of other possible deficiencies, such as the reduced availability of some nutrients, through the application of fertilizers. Legumes in particular can benefit from the acidity correction [10] but also from the application of phosphate fertilizers, which result in an increase in the total biomass production of the pasture [11].
The objective of this work, consisting of regular and systematic monitoring of the evolution of the pasture, was to evaluate the simultaneous effect of the holm oak canopy and the application of dolomitic lime on the productivity and quality of a permanent biodiverse pasture, grazed by sheep, in an acid soil.

Material and Methods
In this work, an experimental field of Montado (4 ha), located at Mitra Farm (38 • 53.1 N; 8 • 01.1 W), was monitored. The predominant soil in this parcel is classified as Cambisol, the major reference soil group in Portugal [2], derived from granite and characterized by low fertility, cultivated under a mixed, agro-forestry production systems [1]. These soils are generally characterized by their coarse texture (mostly sandy loam), small cation exchange capacity (CEC), organic matter content (≤1%), and water holding capacity, and strong acidity (pH ≤ 5.5 in water) [2]. This ecosystem consists of dryland biodiverse pastures, Quercus ilex ssp. rotundifolia Lam. trees, and is grazed by sheep. During the experimental period, the whole field was occupied by adult sheep in maintenance, none of which were pregnant or lactating. The animals were kept on each of the plots in continuous grazing, at an annual stocking rate of eight animals per hectare. The sheep flock consisted of a mixture of Merino White and Merino Black breeds ( Figure 1). Herd management consisted only of vaccinations and deworming according to the health management calendar, a hay supplementation during the autumn period, and a commercial protein supplement in the summer. The animals had permanent access to blocks of mineral salt with macro-and micronutrients.
Sustainability 2020, 12, x FOR PEER REVIEW  3 of 19 vaccinations and deworming according to the health management calendar, a hay supplementation during the autumn period, and a commercial protein supplement in the summer. The animals had permanent access to blocks of mineral salt with macro-and micronutrients.

Characterization of the Climate
The Mediterranean climate is characterized by most of the rainfall concentrated in the winter and very dry, hot summers. The average monthly temperature varies between about 8 °C in January and around 25 °C in August, with minimum temperatures close to 0 °C between December and February and maximum temperatures above 40 °C in summer. The annual accumulated precipitation in the region varies between 400 and 650 mm, distributed mainly between October and March, and practically nonexistent in summer. Figure 2 illustrates the thermo-pluviometric diagram of the Meteorological Station of Mitra (Évora, Portugal). This figure shows the evolution of the monthly mean temperature and monthly rainfall between September 2018 and August 2019, and the average monthly rainfall in a 30-year period (between 1981 and 2020). In the agricultural year of 2018/2019 and relative to the normal average for the region, the figure shows a reduction of 60% in the monthly rainfall between December and June (148 mm versus 370 mm).

Characterization of the Climate
The Mediterranean climate is characterized by most of the rainfall concentrated in the winter and very dry, hot summers. The average monthly temperature varies between about 8 • C in January and around 25 • C in August, with minimum temperatures close to 0 • C between December and February and maximum temperatures above 40 • C in summer. The annual accumulated precipitation in the region varies between 400 and 650 mm, distributed mainly between October and March, and practically nonexistent in summer. Figure 2 illustrates the thermo-pluviometric diagram of the Meteorological Station of Mitra (Évora, Portugal). This figure shows the evolution of the monthly mean temperature and monthly rainfall between September 2018 and August 2019, and the average monthly rainfall in a 30-year period (between 1981 and 2020). In the agricultural year of 2018/2019 and relative to the normal average for the region, the figure shows a reduction of 60% in the monthly rainfall between December and June (148 mm versus 370 mm). vaccinations and deworming according to the health management calendar, a hay supplementation during the autumn period, and a commercial protein supplement in the summer. The animals had permanent access to blocks of mineral salt with macro-and micronutrients.

Characterization of the Climate
The Mediterranean climate is characterized by most of the rainfall concentrated in the winter and very dry, hot summers. The average monthly temperature varies between about 8 °C in January and around 25 °C in August, with minimum temperatures close to 0 °C between December and February and maximum temperatures above 40 °C in summer. The annual accumulated precipitation in the region varies between 400 and 650 mm, distributed mainly between October and March, and practically nonexistent in summer. Figure 2 illustrates the thermo-pluviometric diagram of the Meteorological Station of Mitra (Évora, Portugal). This figure shows the evolution of the monthly mean temperature and monthly rainfall between September 2018 and August 2019, and the average monthly rainfall in a 30-year period (between 1981 and 2020). In the agricultural year of 2018/2019 and relative to the normal average for the region, the figure shows a reduction of 60% in the monthly rainfall between December and June (148 mm versus 370 mm).

Soil Sampling, Interventions and Measurements
The original experimental field (2 ha) was previously characterized by soil samples, collected in October 2015 [12], and pasture floristic composition, carried out in May 2016.
In November 2017, an extra area (2 ha) was added to the experimental field because manganese toxicity was identified and the effect of dolomitic lime and its interaction with the trees was intended to be tested. Dolomitic lime (42% calcium oxide (CaO) and 10% magnesium oxide (MgO); 2000 kg ha −1 ) was applied to the surface of the original area of the experimental field (COR). In each area of 2 ha (amended area (COR) and unamended area (UCOR)), twelve grazing exclusion cages (dimensions of 0.50 m × 0.50 m) were installed as geo-referenced sampling points. Half of each of these sampling points was installed under the tree canopy (UTC) and the other half outside the tree canopy (OTC; Figure 4).
Soil samples were kept in plastic bags, air-dried, and sieved. The fine components (< 2 mm) were analyzed using the following methods [13]: (i) pH in 1:2.5 (soil:water) suspension, using the potenciometric method; (ii) organic matter was measured by combustion and CO2 measurement, using an infrared detection cell; (iii) Ntotal with the Kjeldahl method; (iv) P and K were extracted by the Egner-Riehm method, P was measured using colorimetric method, and K content was measured

Soil Sampling, Interventions and Measurements
The original experimental field (2 ha) was previously characterized by soil samples, collected in October 2015 [12], and pasture floristic composition, carried out in May 2016.
In November 2017, an extra area (2 ha) was added to the experimental field because manganese toxicity was identified and the effect of dolomitic lime and its interaction with the trees was intended to be tested. Dolomitic lime (42% calcium oxide (CaO) and 10% magnesium oxide (MgO); 2000 kg·ha −1 ) was applied to the surface of the original area of the experimental field (COR). In each area of 2 ha (amended area (COR) and unamended area (UCOR)), twelve grazing exclusion cages (dimensions of 0.50 m × 0.50 m) were installed as geo-referenced sampling points. Half of each of these sampling points was installed under the tree canopy (UTC) and the other half outside the tree canopy (OTC; Figure 4).
In October 2018, soil samples in these 24 geo-referenced points were collected using a gouge auger and a hammer, in a depth range of 0-0.30 m. This soil was characterized in terms of pH, organic matter (OM), total nitrogen (N total ), phosphorus (P), potassium (K), magnesium (Mg), manganese (Mn), and exchange basis (calcium (Ca 2+ ), magnesium (Mg 2+ ), and potassium (K + )).
Soil samples were kept in plastic bags, air-dried, and sieved. The fine components (<2 mm) were analyzed using the following methods [13]: (i) pH in 1:2.5 (soil:water) suspension, using the potenciometric method; (ii) organic matter was measured by combustion and CO 2 measurement, using an infrared detection cell; (iii) N total with the Kjeldahl method; (iv) P and K were extracted by the Egner-Riehm method, P was measured using colorimetric method, and K content was measured with a flame photometer; (v) Mg and Mn were measured using atomic absorption spectrometry; and (vi) exchange basis were extracted with ammonium acetate.
In December 2018 and March and June 2019, cone index (CI, in kPa) measurements were carried out in the 24 sampling points with a cone penetrometer electronic FieldScout SC 900 (Spectrum Technologies, Aurora, IL, USA). In this study, the CI values between 0 and 0.20 m, as the average of three measurements dates, were used. On the same dates, soil samples for determining the soil moisture content (SMC) in the top 0-0.20 m soil layer were taken at all sampling points using a gouge auger and a hammer. These soil samples were transported to the lab in metallic boxes, weighed, and then dried at 105 °C for 24 h. Once cooled, they were weighed again to establish SMC. The volumetric SMC was then obtained through multiplying the values by the bulk density.

Pasture Sample Collection and Analysis
During the pasture flowering period (May 2016), a floristic inventory of species present in each of the sampling points was carried out. In each sampling area (0.25 m 2 ), the percentage of coverage by each species was recorded.
In the vegetative cycle of 2018/2019, at the end of autumn (December) 2018, of winter (March), and of spring (June) 2019, pasture productivity and quality parameters were monitored at 24 sampling points. Inside each sampling point, pasture was harvested with a portable electric grass shears at 1 cm above ground level and stored in marked plastic bags. The collected pasture samples were then subjected to the laboratory analysis: (i) weighing (to obtain fresh mass or green matter In November 2018, topographic and apparent soil electrical conductivity (EC a ) surveys were completed. Topographic survey of the area was carried out using an all-terrain vehicle equipped with a real-time kinematics (RTK) GPS instrument (Trimble RTK/PP-4700 GPS, Trimble Navigation Limited, Sunnyvale, CA 94085, USA). The relative field elevation (RFE) map was created using the linear interpolation TIN tool from ArcGIS 9.3 and converted to a grid surface with a 1-m grid resolution. To measure the EC a in the experimental field, the Veris 2000 XA contact-type sensor (Veris Technologies, Salina, KS, USA), equipped with a global positioning system (GPS) antenna, was used.
In December 2018, ammonium phosphate (18% of N and 46% P; 100 kg·ha −1 ) was applied to the surface throughout the experimental area.
In December 2018 and March and June 2019, cone index (CI, in kPa) measurements were carried out in the 24 sampling points with a cone penetrometer electronic FieldScout SC 900 (Spectrum Technologies, Aurora, IL, USA). In this study, the CI values between 0 and 0.20 m, as the average of three measurements dates, were used. On the same dates, soil samples for determining the soil moisture content (SMC) in the top 0-0.20 m soil layer were taken at all sampling points using a gouge auger and a hammer. These soil samples were transported to the lab in metallic boxes, weighed, and then dried at 105 • C for 24 h. Once cooled, they were weighed again to establish SMC. The volumetric SMC was then obtained through multiplying the values by the bulk density.

Pasture Sample Collection and Analysis
During the pasture flowering period (May 2016), a floristic inventory of species present in each of the sampling points was carried out. In each sampling area (0.25 m 2 ), the percentage of coverage by each species was recorded.
In the vegetative cycle of 2018/2019, at the end of autumn (December) 2018, of winter (March), and of spring (June) 2019, pasture productivity and quality parameters were monitored at 24 sampling points. Inside each sampling point, pasture was harvested with a portable electric grass shears at Sustainability 2020, 12, 3758 6 of 19 1 cm above ground level and stored in marked plastic bags. The collected pasture samples were then subjected to the laboratory analysis: (i) weighing (to obtain fresh mass or green matter (GM) in kg·ha −1 ); (ii) dehydrated (72 h at 65 • C); (iii) weighed again (to obtain pasture moisture content (PMC) in percent and dry matter (DM) in kg·ha −1 ); and (iv) the dehydrated samples were subjected to standard analysis of wet chemistry according to the Association of Official Analytical Chemists [14] to determine key components of pasture quality, CP and fiber (NDF), both expressed as a percentage of DM. In the "Discussion", CP is expressed graphically in kg·ha −1 , calculated from Equation (1). The "Discussion" also presents a pasture quality index (PQI), calculated according to Equation (2).
where CP is crude protein (kg·ha −1 and %, respectively) and DMis dry matter (kg·ha −1 ); where PQI is the pasture quality index (%), CP is the crude protein (%), PMC is the pasture moisture content (%), and NDF is the neutral detergent fiber (%).

Statistical Analysis
Descriptive statistical analysis (mean, standard variation, and range) was performed for soil and pasture parameters. Then, ANOVA of the data was carried out considering a two factor experiment (soil correction, COR vs. UCOR, and influence of tree canopy, UTC vs. OTC), using "MSTAT-C" software with a 95% significance level (p < 0.05). Because soil correction was not repeated in space, the interactions between fields and replicas were used to generate the error to compare the two fields. The "Fisher" ("Fisher's least significant difference, LSD") test was applied whenever the ANOVA results presented significant differences between factors.

Variability Pattern of the Soil Parameters
The starting point for this study was the evaluation of the soil of the experimental field carried out in October 2015, which showed that this is an acid soil (mean pH H2O = 5.4 ± 0.3) [12].
In addition to acidity, the soils of this region of Portugal have some other characteristics that affect crop productivity. Figure 5a shows the relative field elevation (RFE) map of the experimental field, with undulated relief typical of the southern region of Portugal, and an 8-m difference (220-228 m) between the highest and the lowest altimetry. Altimetric information is important inasmuch as it affects soil drainage and fertility [15], the major drivers of the productivity in dryland pastures [16]. Figure 5b shows the soil apparent electrical conductivity (EC a ) map, with a mean value of 1.7 ± 0.2 mS·m −1 , which reflects the sandy loam texture (sand = 80.6 ± 2.3%; silt = 10.1 ± 1.7%; and clay = 9.3 ± 1.4%) and other structural limitations of these soils [17]. The spatial information that results from the combination of these surveys (altimetry and EC a ) has a great potential for principal component analysis (PCA) in order to identify and define soil management zones) and select the factors that determine crop yield [15,18,19].
In November 2018, a year after soil amendment, the superficial soil layer was evaluated. Tables 1  and 2 show, respectively, the descriptive and inferential statistics of the soil parameters at the experimental field on this date. The main aspects to be highlighted are: (i) The spatial variability of all soil parameters, except for pH H2O is high; the high CV (30-70%) is a result of the simultaneous effect of trees and animal grazing [12]. (ii) The soil amendment was significant in terms of pH H2O (5.58 ± 0.15 in COR areas and 5.30 ± 0.19 in UCOR areas; P = 0.0193) but not in terms of macronutrients availability; only the UCOR areas showed significantly more Mn (66.8 ± 31.4 mg·kg −1 in UCOR areas Sustainability 2020, 12, 3758 7 of 19 and 33.6 ± 16.1 mg·kg −1 in COR areas; P = 0.0000); as the dolomitic lime is rich in calcium and Mg, the Mg/Mn ratio increased from 1.1 in UCOR areas to 2.5 in COR areas (reaching 4.2 UTC), which is a relevant aspect due to the toxicity problems of Mn in acid Cambisols, recognized as the major limitation for pasture and forage production in the Montado ecosystem [2]. (iii) Tree canopy significantly improved all the evaluated parameters (pH, OM, N total , P, and K), except for Mg, Mn, and exchange basis. Mg and exchange basis showed no significant differences, while Mn showed significantly higher values OTC. This positive effect of trees on soil fertility has also been demonstrated by several other studies [6,20].
Sustainability 2020, 12, x FOR PEER REVIEW 7 of 19 relevant aspect due to the toxicity problems of Mn in acid Cambisols, recognized as the major limitation for pasture and forage production in the Montado ecosystem [2]. (iii) Tree canopy significantly improved all the evaluated parameters (pH, OM, Ntotal, P, and K), except for Mg, Mn, and exchange basis. Mg and exchange basis showed no significant differences, while Mn showed significantly higher values OTC. This positive effect of trees on soil fertility has also been demonstrated by several other studies [6,20].

Variability Pattern of the Pasture Parameters
The objective of this work was to evaluate the simultaneous effect of the application of dolomitic lime and the tree canopy on the productivity and quality of a permanent biodiverse pasture in an acid soil. Figure 6 shows the predominant botanical species of the biodiverse pasture of the studied field evaluated in spring 2016, UTC and OTC. These nine species represent approximately 80% of the total mean cover. Most of these species are annuals, with pronounced patchiness of vegetation communities and marked seasonality of plant cycles, which germinate according to yearly environmental conditions (e.g., temperature and rainfall); consequently, the species composition of grassland in Montado can have a high inter-annual variability [9]. The most representative species are the Erodium moschatum (UTC) and the Chamamelum mixtum (OTC). This spatial variability of the botanical composition is caused by changes in the microclimate (mainly shade tolerance), soil properties, and livestock grazing [9].

Variability Pattern of the Pasture Parameters
The objective of this work was to evaluate the simultaneous effect of the application of dolomitic lime and the tree canopy on the productivity and quality of a permanent biodiverse pasture in an acid soil. Figure 6 shows the predominant botanical species of the biodiverse pasture of the studied field evaluated in spring 2016, UTC and OTC. These nine species represent approximately 80% of the total mean cover. Most of these species are annuals, with pronounced patchiness of vegetation communities and marked seasonality of plant cycles, which germinate according to yearly environmental conditions (e.g., temperature and rainfall); consequently, the species composition of grassland in Montado can have a high inter-annual variability [9]. The most representative species are the Erodium moschatum (UTC) and the Chamamelum mixtum (OTC). This spatial variability of the botanical composition is caused by changes in the microclimate (mainly shade tolerance), soil properties, and livestock grazing [9].   Tables 3 and 4 show, respectively, the descriptive and inferential statistics of the pasture parameters at the experimental field on three dates: end of autumn (December 2018), end of winter (March 2019), and end of spring (June 2019). The main aspects to be highlighted are: (i) Similar to what was observed in terms of the soil parameters, there is a high spatial variability in pasture productivity (GM and DM; CV = 40-70%) and CP (CV = 30-40%) and a low to moderate variability of PMC and NDF (CV = 5-40%). (ii) The soil amendment showed a positive effect on pasture productivity (GM and DM), which was significant only in the end of the winter. (iii) Tree canopy effect is positive on CP, and significant in the autumn and in the spring, and is negative on pasture productivity (GM and DM), and significant in the winter and spring. This effect of tree canopy on pasture productivity and quality has also been demonstrated by several other studies [6,7,12].

Variability Pattern of the Soil Parameters
The 21st century began with two major challenges for agriculture: (i) to increase productivity; and (ii) to attain greater sustainability of the productive systems. One of the great challenges in modern agriculture is, therefore, creating systems that combine lower levels of input with higher food production and minimum environmental impacts [21]. The projections indicate a 30% increase in the world population between 2015 and 2050 (from 7.3 to 9.5 billion), thus the corresponding increase in food demand requires an increase in agricultural production of~70% [22]. On the other hand, soil is a non-renewable resource on human time scales with its vulnerability to degradation depending on complex interactions among processes, factors, and causes occurring at a range of spatial and temporal scales [22]. It is in this context that it is important to evaluate the possibility of intensifying extensive livestock production by increasing pasture productivity in Mediterranean Montado ecosystem, and simultaneously contributing to soil preservation. Furthermore, grasslands are among the most widespread vegetation types worldwide, covering between 14% and 26% of the Earth's surface, and they play an important role in the carbon cycle and in mitigating climate changes [23].
The sustainability of these large cork and holm oak areas is, nevertheless, currently threatened by forecasted climate changes and inappropriate management practices [24]. In the Southern region of Portugal, it is common practice to apply phosphate fertilizers at intervals of two or three years with the aim of improving the fertility of these poor soils [4]. However, the impact of this improvement on pasture productivity is conditioned by an objective limitation: soil acidity. The yield of crops in these Cambisols developed on granitic formations is generally poor as a consequence of the Mn toxicity [2]. In this study, the surface application of dolomitic lime (2000 kg·ha −1 ) in November 2017 had a direct impact on soil pH H2O (from 5.4 ± 0.3 in 2015 to 5.6 ± 0.2 in 2018), with a clearer effect UTC areas (mean pH H2O = 5.7) than in OTC areas (mean pH H2O = 5.5) (Figure 7). These results show that surface lime application was an efficient practice for increasing pH and reducing the exchangeable acidity [25]. However, the application of amendments does not increase much the soil pH H2O immediately, but rather gradually over time [26]. This is in line with other studies [27], which also reflects the way in which lime was applied in the field (surface application, not incorporated into the soil through mobilization). The application also had a positive and significant effect on the Mg/Mn ratio, increasing from close to 1 in UCOR areas to more than 4 in COR × UTC areas (Figure 7). This aspect is particularly important since Mn toxicity has long been recognized as an important factor limiting plant growth on acid soils [2]. The result of this toxicity is a reduction in the soil productive potential, conditioning, in the case of biodiverse pastures, its floristic diversity, which makes them less interesting from an agronomic and nutritional point of view [2]. According to Halim et al. [26], soil amendment is suitable for two primary categories of problems, including reducing contaminant bioavailability/phytoavailability and improving poor soil health and ecosystem function, since it improves cation exchange capacity (CEC), increases water holding capacity and re-establishes microbial communities as well as soil structure by alleviating compaction.
Appropriate indicators of soil chemical quality include pH H2O , CEC, nutrient availability, favorable elemental balance, and absence of any toxicity or deficiency [22]. The improvement of soil chemical attributes with the superficial application of lime observed in this study follows the trend observed in other works [28]. In the context of the mixed Montado ecosystem, the evaluation of the effect of soil pH H2O correction cannot be dissociated from the tree canopy effect on the soil. Areas corrected and under tree canopies (COR × UTC) showed significantly higher levels of OM (2.13% in COR × UTC areas vs. 1.01% in UCOR × OTC areas), nitrogen (0.12% in COR × UTC areas vs. 0.08% in UCOR × OTC areas), phosphorus (72.3 mg·kg −1 in COR × UTC areas vs. 51.2 mg·kg −1 in UCOR × OTC areas), and potassium (212.3 mg·kg −1 in COR × UTC areas vs. 102.0 mg·kg −1 in UCOR × OTC areas). The positive effect of tree canopy on soil fertility has been documented in several studies [6,12,29]. Improving soil quality (i.e., increasing soil organic carbon (SOC) pool, improving soil structure, and enhancing soil fertility) can reduce risks of soil degradation (physical, chemical, biological, and ecological) while improving the environment [22]. Developing strategies to ensure that SOC pool (key indicator of soil quality and an important driver of agricultural sustainability) stays above the threshold or critical level of 1.0-1.5%, especially in the European semi-arid Mediterranean regions, is essential for reducing soil degradation risks and reversing degradation trends [22].
Sustainability 2020, 12, x FOR PEER REVIEW 13 of 19 ecological) while improving the environment [22]. Developing strategies to ensure that SOC pool (key indicator of soil quality and an important driver of agricultural sustainability) stays above the threshold or critical level of 1.0-1.5%, especially in the European semi-arid Mediterranean regions, is essential for reducing soil degradation risks and reversing degradation trends [22]. In terms of topsoil compaction (0-0.20 m surface layer, where grass roots normally develop [4]), Figure 8a shows a trend towards greater soil penetration resistance in UCOR areas relatively to COR areas in the soil layer between 0.10 and 0.15 m, which will justify in the future a more detailed study to identify possible causes of this behavior and its impact on pasture development. Figure 8b shows, on the other hand, very similar soil penetration resistance patterns in UTC and OTC areas, which do not indicate a possible greater compaction that would result from the animals choosing the tree shade as preferred resting areas. This aspect might be explained by the fact that the animals use these places especially in the summer season, looking for shade, when the soil is relatively dry and, therefore, less susceptible to compaction. On the other hand, increased soil compaction due to the effect of tree roots is expected at greater depths [9].

Variability Pattern of the Pasture Parameters
In the previous section, the contribution of lime application to improve soil fertility and the sustainability of the Montado ecosystem is presented, highlighting the interactions that result from the presence of trees and grazing animals. In this section, the impact of these interventions and interactions on productivity and quality of pasture is evaluated, as final indicators of their success and the possibility of intensifying livestock production, necessary to respond to the growing global food requirements. In terms of topsoil compaction (0-0.20 m surface layer, where grass roots normally develop [4]), Figure 8a shows a trend towards greater soil penetration resistance in UCOR areas relatively to COR areas in the soil layer between 0.10 and 0.15 m, which will justify in the future a more detailed study to identify possible causes of this behavior and its impact on pasture development. Figure 8b shows, on the other hand, very similar soil penetration resistance patterns in UTC and OTC areas, which do not indicate a possible greater compaction that would result from the animals choosing the tree shade as preferred resting areas. This aspect might be explained by the fact that the animals use these places especially in the summer season, looking for shade, when the soil is relatively dry and, therefore, less susceptible to compaction. On the other hand, increased soil compaction due to the effect of tree roots is expected at greater depths [9].

Variability Pattern of the Pasture Parameters
In the previous section, the contribution of lime application to improve soil fertility and the sustainability of the Montado ecosystem is presented, highlighting the interactions that result from the presence of trees and grazing animals. In this section, the impact of these interventions and interactions on productivity and quality of pasture is evaluated, as final indicators of their success and the possibility of intensifying livestock production, necessary to respond to the growing global food requirements.
Grassland productivity and quality depends on the interaction between soil, plant, animal, and climate factors, where the correction of soil reaction and fertility [28] plays a crucial role. In extreme cases, of very acidic soils, soil amendments are elements added to the soil to improve its capacity to support plant life [26]. In this study, the topsoil shows medium acidity [4] and, therefore, what is intended is to evaluate the possibility of optimizing pasture productivity and quality. Grassland productivity and quality depends on the interaction between soil, plant, animal, and climate factors, where the correction of soil reaction and fertility [28] plays a crucial role. In extreme cases, of very acidic soils, soil amendments are elements added to the soil to improve its capacity to support plant life [26]. In this study, the topsoil shows medium acidity [4] and, therefore, what is intended is to evaluate the possibility of optimizing pasture productivity and quality.
In all seasons (autumn, winter, and spring), soil amendment showed a tendency to improve pasture productivity (GM and DM) and quality (higher CP and PMC and lower NDF) (Tables 3 and  4). On the other hand, tree canopy showed a tendency to improve pastures CP, presenting an irregular season-dependent behavior with regard to PMC and NDF and a tendency for a negative effect on pasture productivity (GM and DM). This effect of tree canopy on pasture productivity and quality has also been demonstrated by several other studies [6,7,12]. Dryland pasture development is dependent on the agricultural year in terms of rainfall distribution [4]. In this case, the effects of a relatively normal autumn, with a peak precipitation of about 100 mm in November 2018, a dry winter (accumulated rainfall of about 75 mm), and a late spring (May and June) without precipitation ( Figure  2), were reflected in the SMC of the surface layer (0-0.20 m; Figure 9), with average values of 14.7%, 11.3%, and 3.5%, respectively, in the end of autumn, winter, and spring. This distribution of precipitation tends to favor the autumnal growth of the pasture, but compromises the rest of the pasture cycle and, consequently, the spring peak of production [4]. Indeed, the dry winter imposed lower SMC levels than its equivalent in autumn, which was reflected in the positive and significant effect of lime application on winter pasture productivity (Table 4). According to Efe Serrano [4], the distribution of precipitation between March and June is one of the factors that determine the annual pasture production, mainly due to its influence on the extension of the vegetative cycle. Competition for water is usually the main limiting factor for pasture growth, particularly in regions subjected to summer droughts with high temperatures and incident radiation. Soil water deficit results in lower forage DM yields primarily by limiting leaf area development and reducing photosynthesis due to stomatal closure [6]. In all seasons (autumn, winter, and spring), soil amendment showed a tendency to improve pasture productivity (GM and DM) and quality (higher CP and PMC and lower NDF) (Tables 3  and 4). On the other hand, tree canopy showed a tendency to improve pastures CP, presenting an irregular season-dependent behavior with regard to PMC and NDF and a tendency for a negative effect on pasture productivity (GM and DM). This effect of tree canopy on pasture productivity and quality has also been demonstrated by several other studies [6,7,12]. Dryland pasture development is dependent on the agricultural year in terms of rainfall distribution [4]. In this case, the effects of a relatively normal autumn, with a peak precipitation of about 100 mm in November 2018, a dry winter (accumulated rainfall of about 75 mm), and a late spring (May and June) without precipitation (Figure 2), were reflected in the SMC of the surface layer (0-0.20 m; Figure 9), with average values of 14.7%, 11.3%, and 3.5%, respectively, in the end of autumn, winter, and spring. This distribution of precipitation tends to favor the autumnal growth of the pasture, but compromises the rest of the pasture cycle and, consequently, the spring peak of production [4]. Indeed, the dry winter imposed lower SMC levels than its equivalent in autumn, which was reflected in the positive and significant effect of lime application on winter pasture productivity (Table 4). According to Efe Serrano [4], the distribution of precipitation between March and June is one of the factors that determine the annual pasture production, mainly due to its influence on the extension of the vegetative cycle. Competition for water is usually the main limiting factor for pasture growth, particularly in regions subjected to summer droughts with high temperatures and incident radiation. Soil water deficit results in lower forage DM yields primarily by limiting leaf area development and reducing photosynthesis due to stomatal closure [6].
The abrupt drop in moisture in the topsoil from March (end of winter) to June (end of spring) resulted in a strong decrease in the quality of pasture (with very significant breaks in PMC and CP and a substantial increase in NDF, all expressed in percent of DM), reflected in the abrupt PQI break in the whole area of the experimental field ( Figure 10). Sustainability 2020, 12, x FOR PEER REVIEW 15 of 19  Monitoring of pasture quality over time is critical to define the nutritional value of pastures and design balanced diets for grazing animals [30]. However, these qualitative parameters do not take  The abrupt drop in moisture in the topsoil from March (end of winter) to June (end of spring) resulted in a strong decrease in the quality of pasture (with very significant breaks in PMC and CP and a substantial increase in NDF, all expressed in percent of DM), reflected in the abrupt PQI break in the whole area of the experimental field ( Figure 10). PQI, PMC, CP, and NDF, expressed in percent of DM, are pasture qualitative indicators. Monitoring of pasture quality over time is critical to define the nutritional value of pastures and design balanced diets for grazing animals [30]. However, these qualitative parameters do not take PQI, PMC, CP, and NDF, expressed in percent of DM, are pasture qualitative indicators. Monitoring of pasture quality over time is critical to define the nutritional value of pastures and design balanced diets for grazing animals [30]. However, these qualitative parameters do not take into account the evolution of DM. Calculation of the CP of the experimental field, in absolute terms (in kg·ha −1 based on Equation (1)), quantifies the CP that is truly available at each moment (in this case at the end of each season: autumn, winter and spring). Figure 11 shows the DM (Figure 11a) and CP (Figure 11b) accumulated in each season (both in kg·ha −1 and as a percentage of accumulated total). While the DM production peak takes place in the spring (between 38.8% and 58.9% of accumulated total, depending on the specific area of the field), the CP availability is equally distributed among the three seasons, with a slight preponderance in winter. This approach is more appropriate since it allows the agricultural manager to have an idea of the effective CP availability per hectare. This indicator is much more stable throughout the pasture vegetative cycle ( Figure 12) because it is based on two complementary parameters (while CP in percent tends to decrease during the vegetative cycle, DM in kg·ha −1 tends to increase) and, therefore, responds better to the food need of grazing animals (without resorting to supplementation). into account the evolution of DM. Calculation of the CP of the experimental field, in absolute terms (in kg ha −1 based on Equation (1)), quantifies the CP that is truly available at each moment (in this case at the end of each season: autumn, winter and spring). Figure 11 shows the DM ( Figure 11a) and CP (Figure 11b) accumulated in each season (both in kg ha −1 and as a percentage of accumulated total). While the DM production peak takes place in the spring (between 38.8% and 58.9% of accumulated total, depending on the specific area of the field), the CP availability is equally distributed among the three seasons, with a slight preponderance in winter. This approach is more appropriate since it allows the agricultural manager to have an idea of the effective CP availability per hectare. This indicator is much more stable throughout the pasture vegetative cycle ( Figure 12) because it is based on two complementary parameters (while CP in percent tends to decrease during the vegetative cycle, DM in kg ha −1 tends to increase) and, therefore, responds better to the food need of grazing animals (without resorting to supplementation).   The tree canopy effect on pasture is a direct consequence of the extent to which they modify the microclimate and soil properties [6]. In particular, this effect in pasture CP (Figure 11) was only advantageous in autumn, tending to have greater absolute availability of CP in OTC areas in winter and spring due to the marked growth in productivity (DM) in these areas. A reduction in the quantity and quality of light UTC affects directly the physiological processes of plants, decreasing pasture carbohydrate manufacture and net DM production [6].

Conclusions
This study showed that amendment of soil acidity based on dolomitic lime application is a slow and gradual process that improves soil Mg/Mn ratio and has a positive impact on pasture productivity and quality. Improving soil fertility and, consequently, pasture productivity/quality is a fundamental strategy that promotes greater sustainability of livestock production in two ways: (i) incorporation of organic matter, reducing the risks of soil degradation processes; and (ii) intensification of animal production systems, allowing higher number of animals per hectare and reducing the animals' need for supplementary feeding during critical periods. Pasture crude protein availability (CP, in kg ha −1 ), since it is based on both pasture productivity (DM, in kg ha −1 ) and quality (CP, in %), proved to be a very practical indicator of the contributions of tree canopy and soil acidity correction to the holistic management of the Montado ecosystem. In the future, the common practice of phosphate fertilizer application in these acid soils of Mediterranean region to improve pasture productivity and quality must be complemented with previous and repeated application of dolomitic lime. In the context of climate change, the evaluation of the pasture biodiversity through floristic In Figure 11, it is possible to observe the increased CP (kg·ha −1 ) production in COR areas (275.9, 454.1, and 317.6 kg·ha −1 , respectively, in autumn, winter, and spring) relative to UCOR areas (201.2, 257.9, and 276.6 kg·ha −1 , respectively, in autumn, winter, and spring). According to Crusciol et al. [25] and Halim et al. [26], this yield improvement provided by the surface lime application results from a positive effect on the soil nutritional condition, which is reflected in plant nutrition and chlorophyll content.
The tree canopy effect on pasture is a direct consequence of the extent to which they modify the microclimate and soil properties [6]. In particular, this effect in pasture CP ( Figure 11) was only advantageous in autumn, tending to have greater absolute availability of CP in OTC areas in winter and spring due to the marked growth in productivity (DM) in these areas. A reduction in the quantity and quality of light UTC affects directly the physiological processes of plants, decreasing pasture carbohydrate manufacture and net DM production [6].

Conclusions
This study showed that amendment of soil acidity based on dolomitic lime application is a slow and gradual process that improves soil Mg/Mn ratio and has a positive impact on pasture productivity and quality. Improving soil fertility and, consequently, pasture productivity/quality is a fundamental strategy that promotes greater sustainability of livestock production in two ways: (i) incorporation of organic matter, reducing the risks of soil degradation processes; and (ii) intensification of animal production systems, allowing higher number of animals per hectare and reducing the animals' need for supplementary feeding during critical periods. Pasture crude protein availability (CP, in kg·ha −1 ), since it is based on both pasture productivity (DM, in kg·ha −1 ) and quality (CP, in %), proved to be a very practical indicator of the contributions of tree canopy and soil acidity correction to the holistic management of the Montado ecosystem. In the future, the common practice of phosphate fertilizer application in these acid soils of Mediterranean region to improve pasture productivity and quality must be complemented with previous and repeated application of dolomitic lime. In the context of climate change, the evaluation of the pasture biodiversity through floristic composition patterns will be fundamental for monitoring the adoption of the recommended management practices. Funding: This work was funded by National Funds through FCT (Foundation for Science and Technology) under the Project UIDB/05183/2020 and by the projects PDR2020−101-030693 and PDR2020−101-031244 ("Programa 1.0.1-Grupos Operacionais").