Water Related Properties to Assess Soil Quality in Two Olive Orchards of South Spain under Different Management Strategies

: Soil quality is usually assessed through the measurement of selected soil properties. However, in spite of the diversity of the chosen properties, use of the soil water retention curve, like the pressure head or the speciﬁc water capacity at the inﬂection point, provides relevant information of degradation or improvement of soil. The main aim of this study was to evaluate the methods based on these indices in the evaluation of short-term changes of olive cropped soils under typical Mediterranean agricultural conditions. For this reason, soil properties (bulk density, hydraulic conductivity, aggregates stability, and organic matter content) were measured in a short-term trial settled in two olive orchards under different soil managements: tillage and cover crop. Several sampling areas were also distinguished: (i) along the inter tree row and under the canopies’ projection and (ii) at 0–10 cm and 10–20 cm depth. In addition, water retention curves were determined and ﬁtted using two models (van Genutchen’s and Kosugi’s) in order to obtain the inﬂection point and therefore the S index. This index is the maximum value of the slope of the soil water retention curve and is related to soil quality. At the two sites, changes in soil management, even after a brief period of two years, had a relatively quick effect, especially in organic matter content along the inter tree row. The use of indices based on soil water retention curves helps to detect soil degradation or improvement changes. Future research, including the inclusion of more soil types and longer time periods, might lead to the development of more reﬁned tools for the assessment of soil health.


Introduction
Soil erosion is a severe environmental threat for the future of our society [1], especially on steep slopes and under a Mediterranean climate where the rainy season starts after a hot and dry summer season without a protective plant cover [2]. It is precisely in those areas where an important part of traditional olive groves is located. Andalusia, in southern Spain, is the first oil-producing region of the world. In the region, about 25% of the olive orchards grow on steep terrain with average slopes of over 20%, while another 59% of the area is cultivated in slopes ranging between 5% and 20%. Therefore, the food supply and conservation of the Earth pose a challenge for the society; in the words of Sposito [3], to sustain the "genius of soil", which means the optimization of the natural resources like the soil, respecting the complexity of its ecosystem.

Studied Area
The experimental trials were carried out in two representative olive orchard farms with different soil characteristics. The first one was called Pedrera, the name of the closest village, and was located in the southeast of Seville province (Spain), 37 • 13 N; 4 • 53 W. It is a rainfed olive orchard that was planted in 1999 with a frame of 6 × 7 m 2 with the variety "Hojiblanca". Test plots were established in an area with an average slope of 4%. The soil, a clay loam, belongs to subgroup Typic Calcixerept [24]. It is a well-drained soil with abundant stoniness, about 36% of carbonate content, with low organic matter content at 1.6%. The climate of the area is Mediterranean, Csa, within the Köppen-Geiger climate classification [25], temperate with dry and hot summers. The mean annual temperature is 18 • C and it has an annual rainfall of 550 mm. The second farm, called Benacazón, is the neighboring village, and is located west from the capital of the Seville province (Spain), 37 • 21 N; 6 • 11 W, in a small spur of reduced elevation roughly parallel to the Guadalquivir river course. It is an irrigated farm of green table olives planted with the variety "Gordal" in 1980 in a frame of 8 × 7 m 2 . Test plots were established on an average slope of 11%. The soil, a sandy loam, belongs to the subgroup Petrocalcic Palexeralf [24]. It is a well-drained soil, with no appreciable stoniness, an average content of 1.3% organic matter, and 28% of carbonate content. The climate is Mediterranean, similar to the other farm, Csa, with a mean annual temperature of 18.6 • C and annual precipitation of 534 mm [26].
In both farms, two runoff plots parallel to the maximum slope gradient were installed using collector tanks connected by primers. Two different soil management systems were chosen (one per plot). In Pedrera, settled in May 2002, were 12 × 65 m 2 , while in Benacazón, installed in June 2003, the plots' dimensions were 8 × 60 m 2 . One management system was conventional tillage (hereafter referred as tillage) consisting of weeds control by 2-4 cultivator passes at a depth of 15 cm. The number of passes per year varied depending on weather conditions and, therefore, growth of adventitious vegetation. The usual distribution is one pass at the beginning of the year (after harvesting in autumn), followed by several passes (between one and three) at the end of winter and spring.
The second soil management system was temporary cover crop, hereafter called cover. This management system consisted of a sown cover crop of ryegrass (Lolium multiflorum) along the inter tree rows every year, with the first autumn rains. The cover was sown by hand at a seed density of 40 kg ha −1 and fertilized during the winter with 14-14-16 at a dose of 270 kg ha −1 . In early spring (between mid-March and mid-April) depending on annual rainfall, the cover was chemically killed with paraquat 12% + diquat 8% to avoid water and nutrient competition with the olive trees according to the local mowing date recommendations [27].
These two managements were set up in Pedrera and Benacazón before closing the runoff plots. Thus, soil management began in fall 2001 and spring 2003, respectively. The previous soil management systems in the whole farms were conventional tillage in Pedrera and spontaneous cover in Benacazón.

Field Measurements
Field measurements and soil sampling for laboratory analysis were performed during January and April 2005 along the inter tree rows (X) and under olive trees' canopies (C) at two different depths, 0-10 cm (S) and 10-20 cm (P). In the case of the tillage plots, measurements and sampling coincided when soil compaction was maximum, just before the first plough pass in early spring.
Selected soil chemical and physical properties were measured: organic matter (OM), bulk density (BD), saturated hydraulic conductivity (Ks), and macroaggregate stability (AE). In addition, water retention curves (WRC) were determined in undisturbed soil samples (samples in 50 × 50 mm 2 cores) during a drying cycle using a sand box and a sand/kaolin box for pF determination (Eijkelkamp Giesbeek, The Netherlands), which range to 1 wcm and from 1 to 5 wcm, respectively. For lower water potentials of higher than 5 wcm of suction, disturbed soil samples were analyzed using a dewpoint potentiameter (WP4C, Decagon, Pullman, WA, USA). The number of samples taken, or measurements per management, area, and depth was eight for BD, AE, and WRC, and four for the rest of the parameters. Detailed information about the protocols can be found in [28].

Water Retention Curves and S Index
The soil water retention curve has been characterized by several equations, what Raats denominated subclasses of a van Genuchten type [29,30] relating the saturation degree, S e , a normalized moisture content, θ, using the maximum, θ 0 , and minimum value, θ r with the pressure head, h, normalized with the parameter α, using some auxiliary parameters n and m.

The water retention equation is
The value of the pressure head at the inflection point of the soil water retention curve, h IP , is Finally, the value of the Dexter S index is Even though the two indices h IP and S are basically similar, use of the logarithmic transformation in the pressure head implies some changes in their respective final equations. Inclusion of the logarithmic transformation could not be justified in some cases [31], but in this report it is kept.
In some soils, the simplest Brutsaert equation [32], a particular case of the van Genuchten subclass with m = 1, fits very well the experimental data of the water retention curve. Another possible alternative is the Kosugi equation [33] based on the lognormal probability distribution function with three parameters [34], the pressure head at the inflection point, h 0 , the pressure head at the air-entry state, h c , and the shape parameter σ, for the pressure heads greater than the air-entry value, using the complementary of the error function, erfc [35]: Water 2019, 11, 367

of 12
In this equation the value of the Dexter S index is with ∆θ The two soil water retention equations were fitted to the measured data using the Rosenbrock optimization algorithm [36] by minimization with constraints of the sum of the squares of the differences between measured and computed values as the objective function.

Statistical Analysis
Basic statistical analysis (ANOVA and mean comparisons with a Tukey significance level of p < 0.05) were carried out to evaluate the possible differences on soil properties regarding management (tillage, cover), areas (inter tree row, under canopies), and depths (0-10 cm, 10-20 cm). The data was analyzed using Statistix software v.8.1 (Analytical Software, Tallahassee, FL, USA). Table 1 shows a summary of the soil properties measured at the two farms (Pedrera and Bencazón), managements (tillage or cover), depths (0-10 cm or 10-20 cm) and areas (inter tree row or under canopies).

Field Measurements
At the Pedrera farm, OM was slightly greater in the cover plot than in the tillage one, especially along the inter tree row, which denotes a quick response of the soil to the establishment of the cover. On the contrary, areas and depths presented similar content. BD was significantly different in all cases (managements, areas and depths) with the highest values at the top soil (0-10 cm) of the inter tree row of the tilled plot, due to the sampling date, before the first tillage pass. Differences in AE were detected regarding depths, denoting that an interaction between managements and areas exists.
At Benacazón site, the OM was significantly different between the soil management systems. The tilled plot showed in general terms greater values than the cover plot, suggesting the short-term effects of tillage in a plot subject to a traditional system, which included a spontaneous cover crop. These effects can be attributed to the homogenization of the profile by mixing different soil layers. Clear differences between areas were noticed, due to the small BD measured under the canopies in all cases. Some interactions between depths and areas were observed.
Concerning the Ks, Pedrera plots showed lower values than Benacazón due to differences in the clay content at both sites. Moreover, this property presented a great variability due to the occasional appearance of surface cracks during the infiltration tests at Pedrera. However, the highest values were found under the canopies, possibly due to the presence of olive tree roots and the associated porosity. Table 1. Average values and standard deviation (in brackets) at the two sites (Pedrera and Benacazón) and managements (tillage and cover) for the different areas: SX (0-10 cm inter tree row), PX (10-20 cm inter tree row), SC (0-10 cm under tree canopies) and PC (10-20 cm under tree canopies). Soil properties are bulk density (BD), saturated hydraulic conductivity (Ks), macroaggregate stability (AE) and organic matter (OM).

Soil Quality and S Index
The mode of the effective pore radius of the different treatments, the ordinate of the inflection point in the fitted soil water retention curves, are depicted in the box and whiskers plots of Figure 1. The data of this figure were estimated with the Kosugi water retention curve, whose results were similar to those computed with the van Genuchten equation, except for some soil samples of the Pedrera farm, whose retention curve did not show a clear air entry state. This property, also found in the SWRC of swelling soils of the zone [37], posed some problems to other authors [38]; Pierson and Mulla [13] suggested the addition of a second-degree polynomial to the van Genuchten water retention equation.
In any case, the plots of Figure 1 indicate the benefit of the cover crop treatment, which increased the range of greater effective pore volume fraction, as compared to the analog volume fraction of the tilled plot soils. This advantage is not as clear for the data of Benacazón farm. At this site, Benacazón, the plots of the subsoil, depth interval 10-20 cm, of the tilled and cover treatment of the soil samples collected under the canopy and in the alley between the tree lines, are greater than the corresponding values of the surface layer, 0-10 cm, revealing alleviation of the compaction caused by the former cover crop in the soil layers by surface tillage in one case and by the root development in the other one. The observed differences of the mean values of samples collected under the canopy and in the inter tree area suggest the localized compaction of the olive fruit harvesting operation where the workers stand by the tree milking the lower branches bearing the fruits.   Figure 2 shows the S index obtained with the mean WRC of Pedrera (X axis) and Benacazón (Y axis). As Dexter [16][17][18] indicated, S indexes over 0.035 are usually found in non-degraded soils. These results confirm the good quality of soils at both sites, and, in particular, Benacazón detected by farmers and agronomists.
However, it is worthy to comment the S index values of the SX plots (0-10 cm along the inter tree row). In Pedrera (S = 0.041), the beneficial effects of the cover can be clearly appreciated even in the brief period of implementation. On the other hand, in Benacazón the relatively small value of the S index, (S = 0.031) could be attributed to the change from a spontaneous cover with autochthonous vegetation to a sown, monovarietal cover of a gramineous. the good quality of soils at both sites, and, in particular, Benacazón detected by farmers and agronomists.
However, it is worthy to comment the S index values of the SX plots (0-10 cm along the inter tree row). In Pedrera (S = 0.041), the beneficial effects of the cover can be clearly appreciated even in the brief period of implementation. On the other hand, in Benacazón the relatively small value of the S index, (S = 0.031) could be attributed to the change from a spontaneous cover with autochthonous vegetation to a sown, monovarietal cover of a gramineous. Comparing the S index values from the average WRC with some of the mean measured soil properties, no clear trend is evident; Figure 3. The lack of significant correlation might be the result of the confluence of several factors, such as the spatial variability of the measured properties, and the interaction among management Comparing the S index values from the average WRC with some of the mean measured soil properties, no clear trend is evident; Figure 3. The lack of significant correlation might be the result of the confluence of several factors, such as the spatial variability of the measured properties, and the interaction among management systems, areas, and depths. Futhermore, for most of the cases, a change in the management systems was not always reflected by a significant variation of soil properties or S index values. As was previously discussed, this effect was clearer at SX at both sites, which can be explained due to the fact that top soil at the inter tree rows is the most influenced zone within the two managements.
systems, areas, and depths. Futhermore, for most of the cases, a change in the management systems was not always reflected by a significant variation of soil properties or S index values. As was previously discussed, this effect was clearer at SX at both sites, which can be explained due to the fact that top soil at the inter tree rows is the most influenced zone within the two managements. These results suggest that application of the S index for the evaluation of soil quality in Mediterranean soils of woody plantations, such as olive groves, does not seem free of uncertainties. They also indicate the need of more detailed research (e.g., including a larger set of soil properties and expanding the temporal scale) to address whether the aspects discussed in this study, such as the lack of correlation of the S indexes derived from different models, can be satisfactorily resolved for their application in soil characterization and evaluation. In any case, it is remarkable that such a simple index can describe relevant physical properties of the soil.  These results suggest that application of the S index for the evaluation of soil quality in Mediterranean soils of woody plantations, such as olive groves, does not seem free of uncertainties. They also indicate the need of more detailed research (e.g., including a larger set of soil properties and expanding the temporal scale) to address whether the aspects discussed in this study, such as the lack of correlation of the S indexes derived from different models, can be satisfactorily resolved for their application in soil characterization and evaluation. In any case, it is remarkable that such a simple index can describe relevant physical properties of the soil.

Conclusions
The two soils studied showed how the change of soil management, from cover to tillage in Benacazón, or the reverse, from tillage to cover, in Pedrera, have generated a relatively rapid change of some soil properties. This suggests that, for the studied soils, soil degradation processes associated with tillage or with a simplified cover crop occur relatively quickly (Benacazón), but they demonstrate also that the improvement of soil properties associated with the introduction of a cover crop (Pedrera) follows the same trend. Future work expanding the number of soils studied and measuring the evolution of soil properties over a longer time interval should improve these preliminary results.
The observed correlation between the stability of macro-aggregates in water, a classical indicator, and some porosity characteristics of soils, detected by both indices of soil quality, recommends their adoption in future studies of soil health.
Given the limited nature of the dataset, our results suggest the usefulness of the inclusion of this combined approach in future studies due to the relative simplicity of field and laboratory work involved in the determination of selected properties and WRCs. Deeper research advances must be carried out to fully assess its potential with evaluating soil degradation due to different agricultural managements in Mediterranean areas. Funding: This work was supported by project P12-AGR-931 (Andalusian Government), projects ProTerra I-II and Biosuelo (Syngenta) and EU-FEDER funds.