A Comparative Analysis of Soil Loss Tolerance and Productivity of the Olive Groves in the Protected Designation of Origin (PDO) Areas Norte Alentejano (Portugal) and Estepa (Andalusia, Spain)

: Olive groves are Mediterranean systems that occupy more than 2.5 M ha in Spain and 0.352 M ha in Portugal. Assuming the differences between both countries in terms of olive grove regulation and considering their multifunctionality, it is useful to implement agronomic indices to estimate their sustainability. The Soil Loss Tolerance Index (SLTI) and the Soil Productivity Index (SPI) are two such indices. We calculated both indices in the Protected Designation of Origin (PDO) Norte Alentejano (Portugal). The SLTI index was adapted considering speciﬁc variables of the analysed olive groves (i.e., SLTIog). The values obtained were compared with those previously estimated for PDO Estepa (Spain). The negative impacts of erosion and the underlying agricultural practices on the sustainability of olive groves became evident, resulting in decreased soil productivity at the regional level. The SLTIog index showed higher values for crops, being a more realistic tool to analyse sustainability. A higher soil loss tolerance was detected for integrated groves in the PDO Norte Alentejano than for PDO Estepa due to the shorter age of olive cultivation in Portugal, with incipient soil impacts. These indices provide information on the degree of soil erosion, allowing farmers and decision-makers to apply practices to maximise the sustainability of olive


Introduction
Across the European Union (EU), agricultural crops currently occupy 175 million hectares (M ha), of which 5 M ha correspond to olive groves [1]. Olive groves, as socioecological and cultural landscapes, are widely represented in southern Europe, especially in the Iberian Peninsula where they occupy 2.82 M ha [2]. Of these, more than 2.5 M ha are in Spain, resulting from long-term landscape histories with multiple cultural values. In Portugal there are currently 352,000 ha, considered an incipient crop in terms of agricultural expansion and intensification [3,4]. These agricultural socio-ecological landscapes are inherently multifunctional, with a primary economic function, contributing 6% and 1.36% to the national farm income in Spain and Portugal. In addition, they also perform an important societal function by generating 10% and 1.60% of rural employment in Spain and Portugal, respectively. Lastly, olive groves play a key societal role by delivering multiple ecosystem services (ES), and they also shape key environmental a cultural values of mitigate the consequences of erosion in agriculture, maintaining the productive, socioeconomic and environmental dimensions of olive groves in a balanced way [18,24]. A very useful approach to analyse the consequences of erosion on agricultural sustainability is the implementation of generic agronomic indices such as the Soil Loss Tolerance Index (SLTI) and its specific modifications for certain crops, such as olive groves (i.e., SLTIog). These indices stablish the maximum threshold of soil loss tolerated to guarantee the durability and multifunctionality of a crop in the face of different timescale projections, enabling technical staff and decision-makers to implement agricultural conservation measures to increase farming sustainability [25,26]. The Soil Productivity Index (SPI) is equally relevant, quantifying the edaphic productive potential of any agricultural system, providing information on its profitability [27,28].
The scarce existing agricultural regulation in Portugal at the legislative level combined with the intrinsic vulnerability and instability of olive groves, whose development and exploitation is in its early stages in Portugal, demands, as the general objective of this research, the implementation of tools for estimating their productive and environmental sustainability with the aim of guaranteeing a stable food supply to society while maintaining their multifunctionality [10,11,17,24]. Considering the socio-cultural and economic relevance of olive groves in the Iberian Peninsula, and the differences in their management models and legislation between Spain and Portugal, a comparative study was carried out which compares a Portuguese olive grove with European certification and previous studies carried out in an olive grove of a Protected Designation of Origin (PDO) in Southern Spain. In Portugal, the study area was the PDO Norte Alentejano, with a surface of approximately 1350 ha of olive groves largely managed in either traditional, integrated, or intensive ways, and an annual total production of 170 tons (t) of olive oil [13,29]. In Spain, a comparison was made with PDO Estepa (Seville, Andalusia, Spain), which hosts almost 40,000 ha of olive groves and a production of 30,000 t of oil per year [30]. The specific objectives of this study were (a) estimating the SLTI and SPI indices in the PDO Norte Alentejano as a proxy of agronomic sustainability; (b) calibrating and implementing specific modifications of the SLTI for the olive grove management models present in the study area, considering the influence of edaphic variables as potential mitigating agents of erosive processes (i.e., SLTIog); and (c) comparatively evaluating the results obtained with similar studies carried out in the PDO Estepa [31].

Study Area
The main study area of this research corresponded to the PDO Norte Alentejano (Portugal). However, a comparison with previous results postulated for the PDO Estepa (Seville, Andalusia, Spain) was also carried out ( Figure 1).
PDO Norte Alentejano was registered in 2005 at the European level. The area has a hotsummer Mediterranean climate, with an average temperature of 16.3 • C and an average annual rainfall of 582 mm, with gentle slopes and an altitude between 200 and 600 m above sea level (m.a.s.l.). Olive groves are spatially scattered within a landscape mosaic dominated by cork oak and holm oak wood pastures [13,29,[32][33][34]. The soils in which olive trees are grown have a depth of up to 150 cm, with loam or sandy-loam textures, being mainly calcareous and highly susceptible to erosion [35,36]. Specifically, olive farming systems are largely located on land with a slope of less than 15% where laminar/sheet erosion predominates [1,27]. The characteristic soils in this region have a high content of limestone, intrusive igneous rocks (i.e., granite) with quartz, feldspar, plagioclase and mica, and laminar metamorphic and sedimentary schists [14,32]. In Alentejo, water availability is being rapidly transformed due to the presence of the Alqueva reservoir, with public irrigation permissions. Access to these permissions allows deficit/drip irrigation to be applied in certain farms (i.e., the incorporation of up to 1500 m 3 of water ha −1 during water stress), except for intensive management, where water addition responds directly to crop demand with no regulated maximum limit [19,37]. In consequence, olive grove management models implemented in the PDO Norte Alentejano include the following: (a) Traditional mechanised groves in soils with a slope of <20% covering 46% of the olive grove area in this region. The plant density of these crops is <140 trees ha −1 , and the use of chemical fertilisers is allowed [38,39]. (b) Integrated management, with crops representing 45% of the PDO, with a planting density between 140 and 399 trees ha −1 , where the addition of synthetic fertilisers is regulated by technical control agencies [30,39]. (c) Intensive management (9% of the olive grove land within the PDO), with a plant density higher to 400 trees ha −1 , and with use of agricultural machinery, irrigation, and fertilisers indiscriminately permitted [20,37]. The PDO Estepa was registered with the European Union in 2010, covering a total area of 39,694 ha of olive groves scattered across another 20,000 ha of relict vegetation [6,24]. This region has an altitude between 200 and 800 m.a.s.l., and a continental Mediterranean climate characterised by an average temperature of 17.5 • C and 477 mm of annual rainfall [30]. Olive trees are hereby gown on loam and calcareous soils of great depth (i.e., 30-150 cm), with a planting density of 100-500 trees ha −1 [24,31]. From a geological point of view, calcareous soils (i.e., Albariza type) predominate, with a low water absorption capacity together with a medium-high mechanical resistance, being highly susceptible to laminar erosion (i.e., losses of over 35 t ha −1 year −1 of soil) [6,29]. This fact results in a uniform soil loss where materials are dragged from the edaphic surface, generating gullies in the territory, threatening the sustainability of the olive grove [20,28]. Although the majority of these systems are rainfed, 10% of the olive area depends on deficit/drip irrigation, with farms vastly managed under integrated models. Additionally, this region presents an incipient development of organic rainfed olive groves, where only the use of organic fertilisers is allowed (i.e., 500 ha) [6,20,40] PDO Norte Alentejano currently includes three types of olive grove management models, depending on the density of planting and the energy inputs applied to the crop: (a) traditional management, (b) integrated management, and (c) intensive management. The olive groves belonging to this region are located in areas with slopes of up to 15% and are considered to be plain olive groves where the use of agricultural machinery is allowed [19]. Although the majority of olive groves are managed under rainfed systems (i.e., 77% of olive groves), it should be noted that the high water availability in the study area means that all farms can generally implement the addition of deficit/drip irrigation during times of summer drought [41,42].
For the sampling and selection of the farming plots, the study area was stratified according to the three olive grove management models identified and also according to the soil erosion ranks that can be classified based on both the slope and the pattern of land use and land covers, using cadastral and land use information from different geo-spatial data and sources and repositories such as CORINE Land Cover and Epic WebSIG Portugal (i.e., including the use of Digital Elevation Models, DEM) [43,44]. Based on such information and assuming the general character for estimating soil losses from the Universal Soil Loss Equation (USLE) [45,46], the study area was classified into different erosive states following the classification postulated by Moreira-Madueño for agricultural crops according to the slope of the territory [47]. Bearing in mind that in the PDO Norte Alentejano, laminar erosion phenomena predominate, which can lead to edaphic gullies due to the entrainment of soil particles at the surface. The USLE was estimated specifically for the study area as follows: where A: annual soil losses (t ha −1 year −1 ); R: rain erosivity (J ha −1 ); K: soil erodibility (Mg J −1 ); LS: length and degree of slopes (dimensionless, also in %); C: ground cover (dimensionless); P: agricultural conservation practices (dimensionless). Key references in the recent literature addressing the overall Alentejo region and for the case study area were used to calibrate the R and LS factors of the USLE model (Equation (1)) [47,48]. Factor K was alternatively estimated, following the experimental criteria set by Gisbert Blanquer et al. [49] and Auerswald et al. [50]. Factor C was calibrated following the criteria by Gómez et al. [51]. This last factor varies with the type of land management practices implemented based on (a) olive tree density (considered minimum in traditional olive groves, medium in integrated and organic management types, and higher in intensive agriculture); (b) canopy diameter (maximum in traditional, integrated, and organic management olive grove types and minimum in intensive olive groves); and (c) extent (width) of ground cover (partial in traditional and integrated olive groves, maximum in organic olive groves, and null in intensive ones).
Following these criteria, factor C was assigned a value of 0.25 for traditional olive groves in response to the presence of partial vegetation cover and adult olive groves (i.e., >50 years) of approximately a 2.5 m radius and also of inter-rows at least 2 m wide. In the case of integrated olive groves in the study area, considering that the structural characteristics of the traditional olive groves are maintained unchanged, but the inter-rows are widened to 4 m and the olive trees are between 5-49 years old, factor C was assigned a value of 0.16. Lastly, in intensive olive groves this factor was assigned a value of 0.41, due to the absence of plant covers, to the presence young trees (i.e., <5 years) of a 1 m radius, and inter-rows of 2 m. Finally, factor P was assigned an equal value of 1 for all management models, assuming similar characteristics to those considered for the PDO Estepa (i.e., farms are subject to tillage practices without any control of erosion) [31,40]. Table 1 shows the parameters of the USLE model as calibrated for the PDO Norte Alentejano. Table 1. Estimation of soil loss rates (A, in t ha −1 year −1 ), according the Universal Soil Loss Equation (USLE), where R: rain erosivity (J ha −1 ); K: soil erodibility (Mg J −1 ); LS: length and degree of slopes (dimensionless and in %); C: ground cover (dimensionless); and P: agricultural conservation practices (dimensionless). USLE has been calibrated for traditional, integrated, and intensive olive groves, considering all soil erosion levels, in the PDO Norte Alentejano. The erosion categories/levels vary depending essentially on soil erodibility and ground cover (K and C) factors. Based on cadastral information [43,44], all olive farms presented representative crop plots for each possible erosion level considered. Thus, according to the combination of the various classes defined in this PDO for olive grove management models and soil erosion levels, 12 treatments were identified. Within each treatment, four plots were randomly sampled, resulting in a final sample size of n = 48 plots ( Figure 2). In each sampling plot, a 1 km long and 5 m wide transect was delimited, on each of which three soil samples with a core of 112.40 g (weight) and 141.37 cm 3 (volume) were collected. Additionally, soil depth was obtained from the DEM used [43,44]. Samples were dried at 105 • C over 24 h. For the calculation of the three agricultural indices proposed, the main physical-chemical soil variables analysed shown in Table 2 were measured. This task was based on previous research published by Sousa et al. [31], where a similar methodology to the one hereby defined for the PDO Norte Alentejano (Portugal) had been implemented. In the PDO Estepa there are currently two olive grove management models: (a) integrated with rainfed or with an irrigation regime and (b) rainfed organic on isolated plots. Following the USLE (Equation (1)) to estimate laminar erosion in this area, due to the plain relief of the olive groves (i.e., slopes up to 15%), these crops were classified accordingly based on the same two criteria of soil erosion and land management practices (Figure 3), resulting in eight treatments taking into account the representativeness of the integrated rainfed plots in all erosion levels (i.e., null with a 0% slope, slight with a 3% slope, moderate with a 7% slope, and severe with a 15% slope), the existence of integrated irrigated groves only in moderate and severe erosion levels, and the representativeness of organic plots only in lands with null and moderate erosion: Based on this categorisation of olive groves, four plots in each category were sampled, with a final sample size of n = 32 plots. In each of these plots, soil samples were collected from which the soil physical-chemical parameters required to compute the Soil Loss Tolerance Indices applied to agriculture and olive groves (i.e., SLTI and SLTIog) and the Soil Productivity Index (i.e., SPI) were then calculated. my 2021, 11, x FOR PEER REVIEW 8 of 24 Figure 3. Sampling design carried out in previous studies in the PDO Estepa. The water regimes of the study area are shown (i.e., rainfed and irrigated), together with the erosive states (i.e., null, slight, moderate, and severe). The olive grove of the whole study area is managed in an integrated way, with the exception of isolated plots managed in an organic way.

Soil Loss Tolerance Index for Crops (SLTI) and for Olive Groves (SLTIog)
To correctly calibrate the SLTI, soil formation and degradation rates over time must be considered [26,58], along with the depth of degradable soils (i.e., the difference between soil depth and the useful soil depth of 20 cm, which is assumed as a tillage layer and the dry bulk density [53]). Firstly, the weight of soil per unit area was estimated according to the Moreira-Madueño criterion [47]: where W: soil weight (t ha −1 ); Depth: soil depth (cm); DBD: dry bulk density (gr cm −3 ). Based on the SLTI defined by Lombardi-Neto and Bertoni [25] and following the criteria set by Rodríguez Sousa et al. [31], some assumptions/limitations were made to adapt this index to crops. These assumptions are the following: (a) the time period in which a crop yield should be considered as constant is 100 years and (b) that soil erosion rates will be assumed to remain constant during the overall simulation experiment. The simulation times defined were: 1, 10, 25, 50, and 150 years, spanning from the immediate to the distant future. The SLTI was, therefore, calculated according to the following equation: where SLTI: soil loss tolerance index for crops (t ha −1 year −1 ); W: weight of soil (t ha −1 ); E: soil erosion rate (t ha −1 year −1 ); R: soil regeneration rate (t ha −1 year −1 ); t: simulation time (years); 100: time interval in years, in which yields must remain similar to current yields. Soil regeneration, a variable that minimises soil loss due to erosion, was transformed from mm year −1 to t ha −1 year −1 through the following equation: Figure 3. Sampling design carried out in previous studies in the PDO Estepa. The water regimes of the study area are shown (i.e., rainfed and irrigated), together with the erosive states (i.e., null, slight, moderate, and severe). The olive grove of the whole study area is managed in an integrated way, with the exception of isolated plots managed in an organic way.

Soil Loss Tolerance Index for Crops (SLTI) and for Olive Groves (SLTIog)
To correctly calibrate the SLTI, soil formation and degradation rates over time must be considered [26,58], along with the depth of degradable soils (i.e., the difference between soil depth and the useful soil depth of 20 cm, which is assumed as a tillage layer and the dry bulk density [53]). Firstly, the weight of soil per unit area was estimated according to the Moreira-Madueño criterion [47]: where W: soil weight (t ha −1 ); Depth: soil depth (cm); DBD: dry bulk density (gr cm −3 ). Based on the SLTI defined by Lombardi-Neto and Bertoni [25] and following the criteria set by Rodríguez Sousa et al. [31], some assumptions/limitations were made to adapt this index to crops. These assumptions are the following: (a) the time period in which a crop yield should be considered as constant is 100 years and (b) that soil erosion rates will be assumed to remain constant during the overall simulation experiment. The simulation times defined were: 1, 10, 25, 50, and 150 years, spanning from the immediate to the distant future. The SLTI was, therefore, calculated according to the following equation: where SLTI: soil loss tolerance index for crops (t ha −1 year −1 ); W: weight of soil (t ha −1 ); E: soil erosion rate (t ha −1 year −1 ); R: soil regeneration rate (t ha −1 year −1 ); t: simulation time (years); 100: time interval in years, in which yields must remain similar to current yields. Soil regeneration, a variable that minimises soil loss due to erosion, was transformed from mm year −1 to t ha −1 year −1 through the following equation: where R: soil regeneration (t ha −1 year −1 ); R mm : soil regeneration in mm year −1 ; DBD: dry bulk density (gr cm −3 ). Despite the relevance of SLTI for agriculture, there is a need to devise more specific indices adapted to each type of crop. Thus, in the case of olive groves, a specific index exists that was published by Sousa et al. [31], defined as the Soil Loss Tolerance in olive groves (i.e., SLTIog). In order to establish a comparative study between the results in such research and the experiments described in this study, SLTIog was also computed for the Portuguese study area.
SLTIog incorporates the influence of key soil variables as mitigating agents for soil erosion processes. With the aim of excluding variables already considered as USLE factors, the main edaphic parameters that can mitigate soil erosion processes are the percentage of gravels by creating greater resistance to soil loss, porosity, contributing to the generation of edaphic aggregates, and organic matter, whose percentage is directly correlated to the increase in soil fertility [59,60]. Thus, firstly, it is necessary to standardise the values of each variable into a range between 0 and 1 in order to facilitate the comparison between the soil erosion mitigating potential across all management modes using a standardised erosion retardation factor (sERF). In order to correctly implement the SLTIog, the described edaphic variables were standardised according to the following MinMax scaler methodology [61]: where sX: standardised variable (dimensionless, from 0 to 1); X: original value of the variable; Xmin and Xmax: minimum and maximum values of the original variable.
To obtain the final standardised values, a linear regression model was applied, including, as described in Equation (5), a standard soil erosion retardation factor, which will become 0 in the case that no soil erosion exists: where sERF i : standardised erosion retardation factor of erosion level i (t ha −1 year −1 ); α: model intercept, dimensionless; β, γ, δ: coefficients of each independent variable, dimensionless; sG i , sPor i and sOM i : standardised gravels, porosity, and organic matter values of erosion level i. For the incorporation of this factor in the final SLTIog equation, the standardisation of the sERF i factor was removed, resulting in the soil erosion retardation factor (ERF i ): where ERF i : erosion retardation factor of erosion level i (t ha −1 year −1 ); sERF i : standardised erosion retardation factor of erosion level i (t ha −1 year −1 ); sERFmax and sERFmin: maximum and minimum value of the standardised erosion retardation factor (t ha −1 year −1 ). In line with the criteria defined by Sousa et al. [31], it must be considered as a limitation of this model that the ERF i factor is only reliable for crops where soil variables are defined by the following non-standard value ranges: (a) gravels, between 0 and 40%; (b) porosity, between 30 nd 90%; and (c) a range of 0-4% for organic matter. According to these criteria, the SLTIog modified formulation is detailed as follows: where SLTIog: soil loss tolerance index for olive groves (t ha −1 year −1 ); W i : weight of soil for erosion level i (t ha −1 ); R: regeneration soil rate (t ha −1 year −1 ); t: simulation time (years); E i : erosion of level i for the study area (t ha −1 year −1 ); ERF i : erosion retardation factor of erosion level i (t ha −1 year −1 ); 100: time period over which agricultural yields must be constant.
All statistical analyses were carried out using RStudio software, using the car library and the agricolae, dplyr and PMCMRplus packages, considering a level of significance of α = 0.05 in all analyses [62].

Soil Productivity Index (SPI)
The SPI, as formulated by Duan et al. [27], was calibrated and adapted for each treatment and erosion level of the study area: (9) where SPI i : soil productivity index in erosion level i (dimensionless, values between 0 and 1); SM i , SpH i , SOM i , SK i , SCl i : sufficiency values, in erosion level i, of moisture, pH, organic matter, potassium, and clays, respectively; W i : weight of soil in erosion level i (t ha −1 year −1 ); E i : erosion rate i (t ha −1 year −1 ); R: regeneration soil rate (t ha −1 year −1 ); t: simulation time (years); Wmaximum: maximum soil weight among all treatments (t ha −1 year −1 ).
To estimate this index (Equation (9)), the influence of the main variables linked to edaphic productivity (i.e., moisture, pH, organic matter, potassium, and clays) must be considered [40,63]. To properly calibrate the SPI, the sufficiency values obtained for each variable in each sampled treatment must be calculated (i.e., dimensionless values rescaled between 0 and 1). These values were estimated from the empirical data obtained according to Duan et al. [27,28,63]. Thus, the higher the SPI value, the higher the soil productivity. In this sense, SPI values < 0.4 correspond to a low soil productivity, values between 0.4 and 0.8 indicate an average soil productivity, and values > 0.8 correspond to high soil productivity [28]. Table 3 shows the results of the soil physicochemical variables estimated in the olive grove of the study area, with computed soil regeneration rates of 0.01 mm year −1 . These results, including the standardised variables and sufficiency values, are of direct use in calibrating the agronomic indices.

Soil Properties of the PDO Norte Alentejano
In line with these results and assuming that deficit/drip irrigation is incorporated into the crop for all olive groves studied, a detrimental effect of erosion on soil conservation could be observed. In general, dry bulk density presents higher values in intensively managed olive groves, whereas lower soil densities are measured for traditional and integrated olive groves. This last parameter increased in the most eroded plots by 13.99%, 17.27%, and 12.21% in the case of traditional, integrated, and intensive olive groves compared to plots where null soil erosion was detected. Soil depth decreased as soil erosion level increased, by 45.21% in the case of traditional olive groves, by 40.76% for integrated olive groves, and by 44.14% for intensive olive groves, with the latter occupying the shallowest soils across the PDO Norte Alentejano. These trends influence soil weight, with integrated olive groves generally presenting the highest values and with the maximum weight having been calculated at 16,870.20 t ha −1 in plots with null erosion. In contrast, the lowest values for soil weight were found in intensively managed olive grove plots.
Regarding soil texture, the study area was predominantly loam or sandy-loam, except for integrated plots with null or slight erosion, which showed a sandy-clay-loam texture. Taking into account that water runoff took place in fine soils, gravel content decreased in plots with high erosion, with soils having null gravels in intensive olive groves with moderate and severe erosive levels. Additionally, in this sense, the results indicated a slight decrease in the soil content of sand as the level of soil erosion increased due to a reduction of permeability and development of runoff processes. Similarly, the soil clay content also decreased with higher erosion processes that are aligned with agricultural intensification.
Structurally, the highest values of soil moisture were obtained for intensive olive groves. A decrease in porosity as soil erosion increased was observed, corresponding to values of 20.03%, 25.39%, and 26.54%, respectively, for traditional, integrated, and intensive olive groves, with the latter presenting the lowest values. Table 3. Mean and standard error (x ± SE) values for each of the edaphic variables analysed (DBD: dry bulk density; D: soil depth; W: soil weight; Eq: equivalence 1 t soil; G: gravel; S: sands; Sl: silts; Cl: clays; T: texture, where L: loam, SL: sandy-loam, SCL: sandy-clay-loam; P: porosity; M: moisture; OM: organic matter; K: potassium). Their units are specified according to the different management options (i.e., traditional, integrated, and intensive) and degrees of erosion (i.e., null, slight, moderate, and severe) defined for the olive grove in the study area. The standardised values for gravels, porosity, and organic matter (sG, sP, sOM), and the sufficiency values for moisture, pH, organic matter, potassium, and clays (i.e., SM, SpH, SOM, SK, SCl), used for the estimation of different agronomic indices, are also attached. Soil pH in Portuguese farms was generally slightly basic, with values above 7.00 for all types of olive groves examined. However, a mild acidification was observed in the most eroded and intensive olive plots. Organic matter losses of up to 39.95%, 32.56%, and 65.52% were also measured in traditional, integrated, and intensive olive groves, respectively. This was largely driven by soil erosion processes, with intensive agriculture showing the lowest overall carbon content. Edaphic potassium concentration, although decreasing with soil erosion (i.e., up to 47.30%), showed lower values for traditional and integrated olive groves than in intensive ones. Regarding to the estimated sufficiency values, high levels of soil moisture content are the norm at the sub-regional level, along with high pH sufficiency values. In contrast, the sufficiency of organic matter and clays in soils presented mediumhigh values in traditional and integrated olive groves, the latter being the treatment with the highest levels measured under null erosion, whilst this value decreased notably in olive groves under intensive management and with moderate to severe erosion. Finally, sufficiency for potassium was higher in intensively managed olive plots.

Soil Loss Tolerance and Soil Productivity in the PDO Norte Alentejano
In advance to showing the results of the two agronomic indices for the PDO Norte Alentejano, Table 4 shows the statistical results obtained for the multiple linear regression calculation on the basis of which the sERF factor was developed. These results show an R 2 = 51.5%, which can be considered as highly significant (i.e., p-value < 0.001 ***). Table 4. Coefficient, standard error, and significance for each variable considered in the multiple linear regression to define the standardised erosion retardation factor (sERF). Significance values are as follows: not significant (>0.05), significant (<0.05 *), very significant (<0.01 **), highly significant (<0.001 ***).

Intercept
Normalised Variables Based on this calibration, the independent variables showed a significant influence as mitigating agents of soil erosion. Agronomic indices, including the specific model designed for olive groves, were then executed for all types of crops and treatments (i.e., SLTI, SLTIog, SPI). The results are shown in Table 5. Table 5. Soil Loss Tolerance Index (SLTI), Soil Loss Tolerance Index for olive groves (SLTIog), and Soil Productivity Index (SPI) for the PDO Norte Alentejano, taking into account the implementation of drip irrigation in all treatments. If in any of the simulations the value of any of the indices is zero (0), it is considered to be an unsustainable situation (US).  A negative effect of erosion on soil loss tolerance became evident for all types of olive groves evaluated, with the latter decreasing with higher soil erosion rates and longer simulation timescales. It should be noted that the increase in soil loss tolerance detected in plots without any corrective measures implemented is reflective of a simple mathematical trend, whereby in the absence of significant long-term erosion processes, soil regeneration triggers the formation of a greater amount of soil over time. The general results for SLTI showed decreases of up to 46.40%, 36.31%, and 54.23% for traditional, integrated, and intensive olive groves for 1-year forecasts and of 77.04% and 56.96% for traditional and integrated olive groves in 150-year projections. In line with the SLTI, intensive olive groves were detected to become not sustainable from soil and agronomical perspectives in the long term. The results for SLTIog were very similar, also showing decreases in soil loss tolerance as a function of erosion processes. In general, integrated olive groves presented the highest values across all the estimated indices, while the intensive ones showed the lowest long-term sustainability rates. Although there were no significant differences encountered between the results of the SLTI and SLTIog indices (i.e., p-value > 0.05), the index implemented specifically for olive groves showed slightly higher tolerance values as plot-scale soil erosion rates and projected timescales increased (i.e., a maximum difference of 2.02% and 0.45% between SLTIog and SLTI was detected for 150-year projections on severely eroded soils in traditional and integrated olive groves, respectively).

Olive Grove Management
According to the SLTI and SLTIog indices, in the absence of soil erosion processes, integrated olive groves showed 10.23% and 40.51% more soil loss tolerance than traditional and intensive ones for the case of imminent timescale projections (i.e., 1 year), decreasing, respectively, to 10.22% and 40.40% for 150-year simulations. However, under severe erosion conditions, integrated olive groves continued to show (24.44% and 57.25%) higher results than both traditional and intensive plots for 1-year simulations, and 52.11% in the case of traditional olive groves for distant timescale simulations (i.e., 150 years). This makes the intensive plots with severe erosion unsustainable under such timescale projections. In addition, soil erosion processes were also found to have a negative influence on soil productivity, causing respective declines of 91.43%, 76.92%, and 93.75% in the cases of traditional, integrated, and intensive olive groves under 1-year timescale projections. These decreases went up to 97.14% and 84.21% for the traditional and integrated olive groves when they were considered under a 150-year timescale simulation. Thus, soil environment was unproductive for longer-term projections in the intensive olive groves. Finally, it is important to highlight how all of the olive grove management typologies considered showed low soil productivity rates (i.e., <0.4), except for the traditional and integrated olive grove plots without soil erosion processes, the latter of which were found to be the most productive with an SPI of 0.53.

Previous Results for the PDO Estepa
Results published for the PDO Estepa are briefly summarised in Figure 4. In this area, the linear regression model showed an R 2 > 90%, being a highly significant model (p-value < 0.001 ***). Thus, the generally sustainable character of the olive grove in this region becomes largely evident, except for the case of irrigated crops with severe erosion rates in long-term timescale simulations (i.e., 150 years), calculated according to the SLTI.
The level of soil erosion and of irrigation were thus the key factors contributing to minimising soil loss tolerance. The results suggest an inverse relationship between soil erosion and irrigation with respect to the level of tolerance to soil loss in all olive groves. In general terms, the results showed the sustainability of olive groves in Estepa region, and through the SLTIog, even the olive groves with severe erosion and irrigation may be viable at 150 years. The soil loss tolerance decreased from 135.09 t ha −1 year −1 to 65.56-65.57 t ha −1 year −1 in the integrated rainfed crops due to the erosion processes. The SLTI and SLTIog showed decreases of 20-25% in irrigated farms compared to the rainfed olive groves with equivalent erosion levels. For the organic olive groves, a slight increase was observed in soil loss tolerance regarding the integrated management (i.e., from 135.09 t ha −1 year −1 to 140.73 t ha −1 year −1 in plots with null erosion, and from 107.33 t ha −1 year −1 to 113.63-113.64 t ha −1 year −1 in plots with moderate erosion).
This study area is characterised by low-medium values of soil productivity, the most productive plots being the organic ones with null erosion rates, which presented values around 0.45. The influence of soil erosion resulted in decreases in the productivity ranging between 0.19 in the case of integrated rainfed olive groves with null erosion, and up to 0.02-0.01 in plots with severe erosion. In the case of plots with moderate erosion rates, the presence of irrigation led to a decrease of 40% in the SPI. . SLTI (t ha −1 year −1 (a1) and mm year −1 (a2)), SLTIog (t ha −1 year −1 (b1) and mm year −1 (b2)), and SPI (c, dimensionless) indices for the olive grove types and related soil erosion levels calculated in the PDO Estepa. Each mark on the trendlines represents the timescale projections considered, respectively, corresponding to 1, 10, 25, 50, and 150 years.

Influence of Erosion on the Soil Environment of Portuguese Olive Groves
Although the magnitude of soil erosion depends on the agricultural practices carried out on crops and can be mitigated through the correct use and management of the soil environment (i.e., multifunctional agricultural practices), it also depends intrinsically on geological and morphological factors of the territory [30,51]. In this sense, the slope of the crops is a key factor in determining soil loss, as the steeper the slope, the greater the risk of runoff and entrainment of soil materials, which can lead to the formation of gullies [59,64]. Several studies have shown how integrated olive groves combining the controlled supply of chemical fertilisers and deficit/drip irrigation, with the implementation of partial plant covers, minimises the magnitude of erosion processes, thereby mitigating surface runoff and reducing the decline in agricultural productivity due to loss of soil fertility [30,51,65,66]. On the other hand, intensively management models clearly trigger an increase in short-term production rates that are, in the case of olive oil, aligned with the current increase in market demand. Nonetheless, they also result in increased environmental impacts that potentially hamper the sustainability of crops over longer timescales [21,67]. Despite the relatively short history of olive intensification processes in Portugal, the results show a slight edaphic degradation. In this sense, an increase in dry bulk density was observed for intensive olive groves affected by higher erosion rates due to the characteristic tillage practices of this management, such as the indiscriminate use of agricultural machinery [20,[68][69][70]. On the other hand, integrated crops, where plant covers contributes to increasing soil fertility, showed less compacted soils than traditional and intensive ones [66,68]. A decrease in soil depth and soil weight in plots with higher erosion levels was also apparent. This is likely due to the dragging of materials derived from erosion processes and to the fact that the shallower soils are located in intensive groves [71]. An inverse trend with increasing erosion rates was observed for soils with a high content of gravels, which is in accordance with the mitigating effect of such soils according to the criteria set by Lal [59,60]. The integrated olive groves showed a higher level of stoniness, thus contributing to the maintenance of a stable soil structure [72]. Although the traditional olive groves showed an intermediate percentage of gravels, this was sharply reduced in intensively managed plots, which are generally affected by soil homogenisation works necessary to facilitate the use of heavy machinery [70,73]. Texturally, the highest representation of clays was found in the integrated farms, contributing to the generation of edaphic aggregates, increasing soil stability, and preventing to a large extent the runoff of materials [31,74]. Structurally, less edaphic pores at the higher erosion levels were detected, which can be attributed to soil compaction [75]. Considering that porosity is a mitigating factor for soil erosion, the intensively managed olive groves in the study area were placed over highly compacted soils with few pores, generating highly unstable environments where there are no physical mechanisms to promote the retention of soil materials [31,59,60]. Given that the plots in Portugal were characterised by deficit/drip irrigation [19,41], no clear differences related to soil moisture could be observed between the analysed olive managements due to the equitable water supply in all sampled crops.
Chemically, it was observed how soil erosion drives a slight acidification of soils, with the intensive olive groves showing lower pH values than the integrated and traditional ones. The underlying cause of the observed drop in pH is likely to be the non-controlled use of synthetic fertilisers, mainly nitrogen, phosphorus, and potassium (i.e., NPK fertilisers) in intensive farms [20,37]. Regarding the organic matter content, a clear correlation was detected between its concentration, the type of olive grove, and level of soil erosion. In this sense, the most eroded soils had a lower organic matter content, a situation that is particularly aggravated in intensive olive groves, where the decrease in soil fertility can ultimately impact agricultural yields in the medium to long term [76,77]. Another factor closely linked to olive grove production is the soil content in potassium cation [27,28,63]. While higher soil erosion rates showed a negative relationship with potassium content, because of the loss of edaphic nutrients [71], the highest concentrations of this ion were encountered in intensive groves due to the indiscriminate addition of chemical fertilisers to these crops [21,37].

Agronomic Synthesis for the PDO Norte Alentejano
Despite the great usefulness of the employment of indices based on soil loss tolerance and soil productivity, the applicability of these tools is still incipient [47]. In this sense, there are studies aimed at developing these indices from a methodological point of view, with the research of Li et al. [26] and Duan et al. [28] standing out, while other studies focus their applicability on how erosive processes affect the degradation of calcareous soil environments [58,63]. However, their use in agricultural environments and specifically in olive groves is novel, showing in general terms how erosive phenomena lead, in the long term, to the loss of fertility and soil horizons, minimising the sustainability of the crops [31,[78][79][80]. The present research consolidates, together with other similar studies, the necessary bases for the development of specific indices aimed at quantifying the estimation of the economic and environmental sustainability of these agricultural systems [40,[63][64][65]. Taking into account the intrinsic limitations of the indices developed in this research, in which only the existence of linear relationships between edaphic variables and erosion processes was explored, the SLTIog is a more specific and realistic index than the SLTI as it takes into account the influence of specific soil parameters as mitigating agents for the loss of soil materials in order to estimate the sustainability of olive groves [31].
In general terms, irrigation is considered an agricultural intensification measure that contributes to soil degradation, negatively affecting the long-term sustainability of agriculture [12,40]. However, in the study area, where a deficit/drip irrigation regime predominates in olive groves, the impacts derived from the addition of water to farms as a measure of intensification were minimised, improving water use efficiency and vegetative growth control [79][80][81]. Specifically, the addition of drip water in areas with water stress can nonetheless enhance the magnitude of soil erosion, driving the need for more efficient and integrated soil-water management practices [31,82]. In this sense, certain soil conservation measures, especially those oriented towards increasing the roughness of soil surface or the implementation of herbaceous plant covers, contribute to avoid soil runoff and nutrient loss [83][84][85]. Through these practices, a higher soil organic matter content is facilitated, whilst also improving the formation of soil aggregates and reducing compaction [59,60]. The results of this study made it possible to verify for the PDO Norte Alentejano that it is possible to move towards more sustainable management practices in olive cultivation in traditional, integrated, and intensive management by using the SLTI and SLTIog as monitoring tools, except for intensively managed and severely eroded plots. The higher sustainability levels were detected in integrated agriculture, a fact that is influenced by the conservation practices characteristic in this management model, highlighting the regulated use of synthetic fertilisers and machinery by technical control agencies and the implementation of partial herbaceous plant covers, measures contributing to increased soil fertility and stability and reducing soil degradation [20,24,66]. By contrast, the lower sustainability/durability rates appeared in the intensive olive groves. This is in line with the fact that agricultural intensification responds to a mere demand for increased levels of crop productivity, too often at the expense of key environmental dimensions [37,70].
Finally, it is important to remark how olive groves in the PDO Norte Alentejano showed low-medium values of soil productivity, with the integrated olive grove plots without soil erosion being those with the highest productive potential. Such low productivity rates are due to the low soil potassium contents of the sampled plots and to the equally low percentages of organic matter and clays in the most eroded treatments, which are variables that can be considered key factors in agricultural yields [27,28,59,63,86].

Soil and Agronomic Comparison between the PDO Norte Alentejano and PDO Estepa
After calibrating the USLE model [45][46][47][48][49][50][51], it was observed that both the traditional and integrated olive managements in Portugal showed lower soil erosion rates than the integrated ones in Spain. This is probably due to the longer history of olive growing in Southern Spain [6,14]. Olive groves in Spain have evolved over millenia and expanded early, resulting in highly eroded soils. However, the organic olive groves in the Estepa region showed the lowest erosion rates out of them all due to the employment of organic fertilisers and the implementation of total plant covers [68,83,84]. On the other hand, intensively managed groves in the PDO Norte Alentejano were those with the highest soil erosion rates. In these intensive farms, higher soil erosion rates are correlated with the environmental impacts from the use of agricultural machinery and the addition of fertilisers. These impacts with the absence of herbaceous vegetation covers that increase the retention of soil materials in intensive olive groves make these systems more susceptible to severe erosion processes [20,37,70,84,87].
In comparison with the results reported by Sousa et al. [31], Portuguese olive groves were characterised by higher dry bulk densities, the organic groves in the PDO Estepa being where the least compacted soils were found. This finding results from the lower magnitude of erosion processes and the presence of living plant covers in these crops [88]. The integrated management olive groves in Portugal contained higher values of soil depth and weight, with the traditional and intensive olive groves showing the lowest values. Similarly, soil texture was predominantly loam or sandy-loam in both study areas, with a high limestone content, making them highly susceptible to soil erosion [52][53][54][55]. Soil porosity was also similar in both regions, with the integrated olive groves of Portugal showing slightly higher values than the Spanish ones. In contrast, the olive groves in the PDO Norte Alentejano showed soil moisture contents equivalent to the irrigated groves surveyed in the PDO Estepa. This is explained by the use of deficit/drip irrigation in both areas [17][18][19][20]. Portuguese olive groves had slightly lower pH values than those of the PDO Estepa, constituting neutral or slightly basic soils in accordance with the Fuentes Yagüe classification [89]. The percentage of soil organic matter in the integrated olive groves in Portugal was higher than that of the organic olive groves in Spain, showing a higher soil fertility rates and thus positively influencing their productive yields [63]. However, the intensive managed plots showed lower fertility rates due to the loss of organic matter, which results in higher levels of soil erosion [37,90]. The highest potassium concentration was found in organic crops in Spain, with the Portuguese olive orchards having higher concentrations than those under integrated management practices in the PDO Estepa, but still lower than 170 mg kg −1 , thus limiting soil productivity [27,28].
In agronomic terms, the main limitation of the SLTIog model calibrated for the PDO Norte Alentejano consisted of the fact that the same soil parameters were used as in the model developed by Sousa et al. for the PDO Estepa [31]. Although this model presented an R 2 of 95.5% in Spain, in Portugal only an R 2 of 51.5% was obtained. This decrease in the R 2 is due to the fact that throughout this research the same SLTIog has been implemented for very different types of olive groves in Portugal, resulting in the extreme variability of the data that jointly make up this index. In contrast, data from PDO Estepa were more uniform, which is consistent with the fact that 95% of the land surveyed being managed under an integrated model [6,24,31]. Comparatively, according to the SLTI and SLTIog, the integrated olive groves in Portugal showed the highest soil loss tolerance values, with maximum values of 145.94 t ha −1 year −1 under null erosion conditions. This indicates a degree of sustainability even higher than that in the organic olive groves of the PDO Estepa, which are characterised by values of 140.89 t ha −1 year −1 . If soil loss tolerance is defined as the maximum annual soil degradation value permissible to guarantee the sustainability of an agricultural system over time [58,63,[78][79][80], the higher soil loss tolerance detected in the Portuguese integrated olive groves may be due to their shorter history, a fact that would explain that the soil degradation processes hampering the sustainability of these systems may have not yet appeared [3,4,13]. In addition, both the traditional and intensive olive groves in Portugal presented lower values of soil loss tolerance than the integrated and organic olive groves of the PDO Estepa, the intensive olive groves being those with the lowest results across the analysed timescale projections. Such lower rates of the intensive olive groves in the PDO Norte Alentejano are in line with the expected environmental impacts generally assumed for intensive management models and also the absence of semi-natural herbaceous plant covers that are capable of mitigating the magnitude of soil erosion processes encountered [20,37,70,87]. Although without erosion, it was the Portuguese integrated olive growing that presented the highest soil productivity; with erosive processes, the most productive scenario was the Spanish organic olive groves, given that the environmental practices of this agricultural management (i.e., organic fertilisers and plant covers) contribute to maximising soil conservation and fertility [66,88,90]. By contrast, the lowest productivity scenario was the intensive olive grove under severe erosion. This is explained by soil degradation and loss of nutrients (i.e., potassium), organic matter, and the generation of more compacted soils as a consequence of erosive processes [27,56,76,77].
Taking into account the results obtained, which show how the presence of plant covers in integrated and organic farming has a positive influence on the durability of olive groves, it is highly recommended that the implementation of farming practices characteristic of multifunctional agriculture be promoted from a legislative and political point of view, thus contributing to increasing agricultural sustainability [6,15,83,84,88]. In this regard, it should be noted that the current Common Agricultural Policy (i.e., CAP post-2020) includes within its structure a specific pillar for non-production-based aid, aimed at promoting rural development and soil conservation measures such as those commonly used in organic farming (i.e., minimal use of machinery, non-synthetic fertilisers, and avoiding bare soil), with a commitment to enhancing the value of productive and non-productive ES in agriculture in accordance with the current demands of society, mainly related to stable food supply and food security [24,91,92].

Conclusions
Considering the recent boom of olive grove expansion and intensification in Southern and Central Portugal, we hypothesise that not enough time may have yet elapsed to adequately assess the different soil impacts of different farm and olive grove management models in the study area. In this sense, it is important to highlight how erosive processes are influenced by the topological and edaphological characteristics of the territory, in addition to the agricultural practices carried out. Thus, erosion processes in olive groves are highly influenced by the slope of the territory, which favours runoff phenomena and the appearance of gullies. In general terms, integrated management showed the lowest erosion rates, most probably explained by the implementation of agricultural conservation practices that are generally absent from traditional and intensive management models. The SLTI and SPI indices prove to be relevant and innovative tools that allow us to know the conservation status of olive grove plots and, consequently, to adopt measures aimed at minimising the environmental impacts of agricultural practices. Specifically, the SLTIog, which is specifically tailored for olive groves, provides a more realistic indicator than the SLTI, which is generic for any other agricultural system. By incorporating the influence of different soil erosion mitigating factors, it is possible to estimate the admissible annual soil losses to ensure the maintenance of agricultural production across scales ranging from the short to the long-terms. The results showed a higher SLTI and SLTIog for the integrated olive crops in the PDO Norte Alentejano compared to the traditional and intensive ones, the latter being the management model with the lowest long-term sustainability standards.
Compared to the PDO Estepa in Spain, where olive groves are already the matrix of a monocultural landscape with a deep-rooted history, Portuguese integrated olive groves were detected to be those characterised by higher soil loss tolerance rates, even higher than the integrated and organic olive groves in Spain. However, the results obtained for traditional and intensive olive groves in the PDO Norte Alentejano suggested a lower potential for agronomic sustainability than the Spanish olive groves considering any of the various soil erosion levels. Similarly, the integrated olive groves in Portugal were also those that showed the highest soil productivity (with an intermediate level).
Finally, it is important to remark that although this comparative study provides relevant information on the state and future prospects for soil conservation of Portuguese and Spanish olive groves, it has the limitation of not having considered other possible soil parameters as erosion mitigating agents. Future studies which explore the relationship between the soil environment and agricultural sustainability should lay the foundations for the development of policies that promote the implementation of agri-environmental measures aimed at maximising the productive and environmental sustainability of agriculture. In this sense, it is necessary to carry out research aimed at the generation of ad hoc models adapted to each olive grove management considered in their respective contexts. This would permit assessing the potential consideration of soil variables that become apparent at the regional or local levels, potentially acting as mitigating agents of soil erosion processes, whilst also incorporating statistical models that are important to elucidate possible nonlinear relationships and interactions between each of these variables and final soil erosion rates.