Estimating the Soil Hydraulic Functions of Some Olive Orchards: Soil Management Implications for Water Saving in Soils of Salento Peninsula (Southern Italy)

: Saving water resources in agriculture is a topic of current research in Mediterranean environments, and rational soil management can allow such purposes. The Beerkan Estimation of Soil Transfer parameters (BEST) procedure was applied in ﬁve olive orchards of Salento peninsula (southern Italy) to estimate the soil physical and hydraulic properties under alternative soil management (i.e


Introduction
Sustainable soil management in dry-farmed olive orchards represents a current topic in agricultural research because it is expected that rainfalls patterns (i.e., distributions and intensity) could vary due to climate changes [1,2]. As a consequence, coupled with the expected temperature variations, the aforementioned changes could seriously affect the water availability of Mediterranean orchards [3], and specific investigations aimed at quantifying the impact of soil management on the water availability of specific geographical areas are necessary.
The olive tree represents the main crop in the Salento peninsula (Puglia region, southern Italy). This spread mainly was due to the relatively high environmental aridity conditions of that area (i.e., relatively shallow soils with fine or medium texture, coupled with typical rainfall patterns of semiarid Mediterranean climate) [4]. Also, some olive orchards are not managed following the sustainable and advanced cultivation criteria because owners are often not full-time farmers. For these main reasons, low-cost management techniques (i.e., the no-tillage of the soil) coupled with chemical weed control (or with spontaneous grassing covers) is a widely applied approach by Salento landowners.
It is well known that soil tillage is an agronomic option that can appreciably modify the physical and hydraulic properties of the soil [5]. Tillage type, depth, and machinery characteristics, in fact, represent some major factors able to affect the aforementioned soil properties [6][7][8]. However, strategic minimum tillage without soil layers inversion was suggested to implement "options of good agricultural practices" for soil management (i) to increase the soil water storage, (ii) to improve soil water infiltration, (iii) to reduce evaporative losses from Mediterranean soils or, most recently, (iv) to ensure the control of disease vectors [4,9]. Consequently, because a lack of information for Salento soils exists on this topic, to provide specific findings is necessary.
There is experimental evidence that soil property assessments can represent an effective approach to investigate the influence of soil management on water storage [10], or to evaluate the impact of different soil conservation practices on soil functions [11,12]. Such studies have increased our knowledge on the sustainable management of orchards, but further investigations are needed for specific arid areas of the Mediterranean basin.
Knowledge of the soil water retention curve (WRC; i.e., the relationship between the volumetric soil water content and the soil water pressure head) and the hydraulic conductivity function (HCF; i.e., the relationship between the volumetric soil water content (or soil water pressure head) and the soil hydraulic conductivity) is crucial for adopting appropriate site-specific water management strategies (water volumes, irrigation scheduling, etc.). However, both WRC and HCF are variable in space and time, also because of the adopted soil tillage [13], and the choice of alternative soil management such as minimum tillage (MT) or no-tillage (NT) in orchards may represent a crucial factor because they also allow sustainable weed management [10]. Moreover, WRC-HCF determination is costly, time consuming, and/or practically inapplicable at the territorial scale (i.e., extensive farms, irrigation districts, or at larger scales) [14], and new experimental procedures that are cheap and relatively simplified should be selected and tested as a practical alternative to standard methods [15].
The Beerkan Estimation of Soil Transfer (BEST) procedure by Lassabatère et al. [16] allows to estimate both WRC and HCF starting from relatively easily measured or routinely surveyed soil data, such as the soil texture (i.e., soil particles distribution or, alternatively, clay, silt, and sand contents), the soil water content at the time of experiment, soil bulk density, and the cumulative infiltration by a very simple infiltrometric experiment in the field. Briefly, BEST is based on the analytical infiltration model proposed by Haverkamp et al. [17] and, assumed that that soil hydraulic properties (WRC and HCF) are expressed by given empirical relationships, it estimates the unknown parameters from a pedotransfer function and a simple field single-ring infiltration experiment conducted under saturated soil conditions. The van Genuchten [18] relationship with the condition of Burdine [19] is used for WRC and the Brooks and Corey [20] relationship for HCF [21]. More details on the BEST theory and on the three calculation algorithms can be found in Bagarello et al. [22]. Moreover, to allow easy application of the procedure, free tools were shared (i) to analyze the BEST data using an automatic analysis of multiple infiltration experiments using Scilab coded by Lassabatére, or using Excel and Visual Basic coded by Di Prima [23], and (ii) to build an automated infiltrometer for field experiments (i.e., infiltrometer constructive schemes for the automation of infiltration measures) [24]. Consequently, Agronomy 2020, 10, 177 3 of 17 new approaches for the acquisition of several infiltration runs at the same time and the improvement of the statistical treatment of the data were proposed [25].
To date, several BEST applications have proven the method's soundness in different agro-environmental investigations. In southern China, for example, the method has been applied to investigate the effects of Napier grass management on soil hydrologic functions in a karst landscape and to identify reasonable strategies for maintaining soil hydrologic function [26], or to establish the impact of different vegetation restoration types on saturated hydraulic conductivity, within a large global ecological restoration engineering project [27]. Khaledian et al. [28] applied BEST to investigate the impact of soil tillage on spatio-temporal variation of soil properties under drip irrigation. Di Prima et al. [29] assessed the soil physical quality of a Spanish orchard under three different types of soil management (i.e., no-tillage using herbicides, conventional tillage under chemical farming, and no-tillage under organic farming), and common indicators such as soil bulk density, organic carbon content, or structural stability index were considered in conjunction with capacitive indicators estimated by BEST. The results showed that independent and BEST-derived indicators yielded similar information, suggesting their ability to distinguish soil quality among contrasting soil management types. Therefore, cited references suggest that BEST is a suitable method for spatially distributed investigations, carried out under very varied operating conditions, but other field tests are necessary in different agro-environments to further verify their robustness.
The main objective of this study was to investigate the impact of soil tillage management (minimum tillage and no-tillage) on physical and hydraulic properties of porous media in typical rainfed olive orchards of the Salento peninsula. Specific goals were (i) to estimate the soil physical and hydraulic properties of some olive orchards, (ii) to test the BEST procedure's reliability to estimate the soil hydraulic functions, and (iii) to use the experimental information obtained to quantify the effects of soil management on soil water conservation.

Field Sites
This study was performed in five olive orchards, located in the Salento peninsula (Apulia region, southern Italy), and the inter-row has been subject to intense soil sampling. Overall, the investigated soils were heterogeneous in terms of soil management (i.e., minimum tillage (MT) and no-tillage (NT)) and soil textures (from relatively coarse to fine soil textures). However, given the prevailing soil management for olive orchards in that area, four of the five sites were NT soils. This has prevented to establish comparisons on all selected sites. For each orchard, sampling areas of 1.0-1.6 Ha were selected, in accordance with the purposes reported by Manici et al. [4]. Details on the experimental sites are shown in Table 1, and examples of soil tillage conditions at the time of experiments are depicted in Figure 1. Information on past soil management, that is to say, on the soil tillage in the previous years to the present investigation, was obtained only for the sites of NE and GC. As shown in Figure 1 (ST-MT), the soil at the Sternatia site was subjected to minimum tillage a few days before field measurements. Therefore, for this site, the comparison provides information on the extreme soil tillage conditions that were experimentally observable.

Field and Lab Measurements
A total of 34 BEST experiments were carried out in the spring-summer season of 2017 to obtain a complete hydraulic characterization (i.e., WRC and HCF) of investigated soils. Specifically, the BEST procedure was applied in 5 to 6 randomly selected sampling points at Gagliano del Capo (GC; five sampling points), Neviano (NE), Sternatia (ST), and Andrano (AN) sites, while 11 locations were sampled at Soleto (SO) where a water supply was available. For each sampling point, the application of BEST required performing a Beerkan infiltration and sampling a soil core (10 cm in height by 5 cm in diameter; 196 cm 3 ) to obtain both the soil water content at the time of the experiment, θ0, and the soil bulk density (BD). As common for the application of the method, a ring with a diameter of 15 cm was superficially inserted into the soil (1 cm), and 15 water volumes of 200 cm 3 each were repeatedly poured into the ring. The experimental cumulative infiltration, I(t), was then deduced. Therefore, because each experimental site required about 20 L of water, the lack of water supply in almost all the sites prevented sampling a greater number of points. Saturated soil water content, θs, was estimated from BD assuming a particle density of 2.65 g cm -3 . Following a standard experimental protocol for BEST application, the auxiliary determinations (θ0 and BD) concerned the 0-10 cm soil layer. A disturbed soil sample was also taken at each sampling point to determine both the soil particle size distribution (PSD) and the total organic carbon of the soil (TOC). The PSD was determined with the standard pipette method [30], while TOC was quantified through dry

Field and Lab Measurements
A total of 34 BEST experiments were carried out in the spring-summer season of 2017 to obtain a complete hydraulic characterization (i.e., WRC and HCF) of investigated soils. Specifically, the BEST procedure was applied in 5 to 6 randomly selected sampling points at Gagliano del Capo (GC; five sampling points), Neviano (NE), Sternatia (ST), and Andrano (AN) sites, while 11 locations were sampled at Soleto (SO) where a water supply was available. For each sampling point, the application of BEST required performing a Beerkan infiltration and sampling a soil core (10 cm in height by 5 cm in diameter; 196 cm 3 ) to obtain both the soil water content at the time of the experiment, θ 0 , and the soil bulk density (BD). As common for the application of the method, a ring with a diameter of 15 cm was superficially inserted into the soil (1 cm), and 15 water volumes of 200 cm 3 each were repeatedly poured into the ring. The experimental cumulative infiltration, I(t), was then deduced. Therefore, because each experimental site required about 20 L of water, the lack of water supply in almost all the sites prevented sampling a greater number of points. Saturated soil water content, θ s , was estimated from BD assuming a particle density of 2.65 g cm −3 . Following a standard experimental protocol for BEST application, the auxiliary determinations (θ 0 and BD) concerned the 0-10 cm soil layer. A disturbed soil sample was also taken at each sampling point to determine both the soil particle size distribution (PSD) and the total organic carbon of the soil (TOC). The PSD was determined with the standard pipette method [30], while TOC was quantified through dry combustion using a TOC Vario Select analyzer (Elementar, Langenselbold, Germany). TOC measurements were evaluated because, together with the finest fraction of the soil (clay and silt), the soil structure may have been affected. The clay, silt, and sand percentages, which were determined according to the USDA classification, were used to run BEST [31]. Therefore, the BEST-steady algorithm allowed estimating also the saturated hydraulic conductivity of the soil, Ks. Briefly, the BEST-steady algorithm by Bagarello et al. [22] makes use of both the intercept and the slope (b s and i s , respectively) of the straight line fitted to the data describing steady-state conditions on the I vs. t plot. According to Di Prima et al. [32], Ks can be directly calculated by the following equation: where A and C are constants [32]. More details on the applied methodology can be found in Castellini et al. [31].

Capacitive-Based Soil Indicators Determination
To account for the optimal proportion between water and air into the soil, three capacitive-based soil indicators were considered in this investigation, namely air capacity (AC), plant available water capacity (PAWC), and relative field capacity (RFC). Specifically, the aforementioned soil indicators were calculated by differences (i.e., AC and PAWC) or ratios (RFC) between the main points of the water retention curve (i.e., saturated soil water content, field capacity, and wilting point). The schematic example of Figure 2 shows the formulas for their calculation. Therefore, AC and PAWC were selected to account for the optimal (or not optimal) availability of air and water into the soil. On the other hand, RFC was also considered because it was applied in many investigations (among others, [33,34]). RFC is an effective soil indicator because it partially combines the air capacity and plant available water capacity indicators by expressing soil ability to store air and water relative to the soil's total pore volume [35]. Consequently, since RFC was also selected as the main variable to assess the effect of tillage [35], it seems suitable for the specific purposes of this investigation.
where A and C are constants [32]. More details on the applied methodology can be found in Castellini et al. [31].

Capacitive-Based Soil Indicators Determination
To account for the optimal proportion between water and air into the soil, three capacitive-based soil indicators were considered in this investigation, namely air capacity (AC), plant available water capacity (PAWC), and relative field capacity (RFC). Specifically, the aforementioned soil indicators were calculated by differences (i.e., AC and PAWC) or ratios (RFC) between the main points of the water retention curve (i.e., saturated soil water content, field capacity, and wilting point). The schematic example of Figure 2 shows the formulas for their calculation. Therefore, AC and PAWC were selected to account for the optimal (or not optimal) availability of air and water into the soil. On the other hand, RFC was also considered because it was applied in many investigations (among others, [33,34]). RFC is an effective soil indicator because it partially combines the air capacity and plant available water capacity indicators by expressing soil ability to store air and water relative to the soil's total pore volume [35]. Consequently, since RFC was also selected as the main variable to assess the effect of tillage [35], it seems suitable for the specific purposes of this investigation. For each of the three soil indicators considered, we selected optimal values or critical limits, according to the reference of literature [35,36]. Therefore, the following classes have been considered for AC (cm 3  For each of the three soil indicators considered, we selected optimal values or critical limits, according to the reference of literature [35,36]. Therefore, the following classes have been considered for AC (cm 3 cm −3 ), PAWC (cm 3 cm −3 ), and RFC (-): • AC < 0.10 poor, 0.10 ≤ AC ≤ 0.14 good, 0.14 < AC ≤ 0.26 optimal, AC > 0.26 poor; • PAWC < 0.10 poor, 0.10 ≤ PAWC < 0.15 limited, 0.15 ≤ PAWC < 0.20 good, PAWC ≥ 0.20 ideal; • RFC < 0.6 poor (water-limited soil), 0.6 ≤ RFC ≤0.7 optimal, RFC > 0.7 poor (aeration limited soil).

Data Analysis
WRC and HCF functions were estimated using the BEST-steady algorithm [22], and the obtained soil properties were used to establish comparisons among soils. WRC was used to estimate the selected capacitive-based soil indicators.
In order to verify the reliability of both soil functions and capacitive-based soil indicators obtained by BEST, a first analysis was carried out by subdividing the available dataset into classes of 1) soil texture, 2) soil bulk density, and 3) soil capillarity (λ c ). Specifically, after estimating soil functions (i.e., WRC) and calculating capacitive-based soil indicators (AC, PAWC, and RFC), they were discriminated against for the following: 1.
The macroscopic capillary length, λ c (mm), was calculated as where b is a dimensionless constant commonly set equal to 0.55 for field soils, and ∆θ is the difference between the soil water contents at saturation and at the time of the infiltration experiment (i.e., θ s -θ 0 ) [21]. The first classification is rather applied when is necessary to discretize the soil textures according to the prevailing hydrological behavior, while the others were developed both using literature references or adequately considering the available BD measurements [37][38][39][40]. Consequently, a qualitative assessment of the expected value for a given indicator (i.e., increasing or decreasing), according to the considered classes, was carried out. A further analysis was carried out to compare the water retention estimations obtained by BEST and the corresponding values independently obtained by standard lab methods. Specifically, 9 points of the soil water retention curve of Neviano-SL and Sternatia-CT were obtained by the hanging water column apparatus (for h values ranging from 5 to 100 cm) and pressure cells (for h values ranging from 1030 to 15300 cm) [14]. Moreover, some correlations between variables estimated by BEST (i.e., AC or RFC) and variables independently obtained (BD) were verified.
Finally, two comparisons were performed to establish the impact of soil management (i.e., soil tillage) on soil physical and hydraulic properties for the sites of Neviano and Sternatia.
For each variable considered in this investigation (clay, silt, sand, θ 0 , θ s , BD, TOC, λ c , AC, PAWC, RFC, and Ks), a given dataset was summarized by calculating the mean (M) and the associated coefficient of variation (CV). Arithmetic means were considered for all variables except Ks, which was assumed to be log-normally distributed, as commonly suggested in the literature [41]. Consequently, for the latter, geometric means and associated CVs were calculated using log-normal equations [41].

Field Sites
Results of soils texture characterization classified the no-tilled soil of Andrano (AN) as clay-loam, and the soils of Gagliano del Capo (GC) and Soleto (SO) as sandy-clay-loam (Table 1). Moreover, the two fields at the Neviano (NE) site, although adjacent, showed different soil textures, (i.e., sandy-loam or clay), while Sternatia soils were sandy-clay-loam or clay, according to the soil tillage (i.e., ST-MT or ST-NT, respectively, under minimum tillage or no-tillage). Therefore, the impact of both soil texture and soil tillage on physical and hydraulic properties of the soil were investigated at Neviano and Sternatia sites, respectively.
Overall, TOC values were usual for Mediterranean soils [4,31], with TOC values that ranged between 10 and 15 g kg −1 (Table 1). However, since previous investigations on these same soils (except for GC) have shown no significant differences in soil stability structure (i.e., obtained as a function of the contents of clay, silt, and soil organic carbon) [4], possible effects of soil organic carbon could be indirect (i.e., for the improvement of the physical and hydraulic properties of investigated soils). Figure 3 shows the differences in terms of WRC and HCF when investigated soils were grouped by soil texture classes, bulk density, and soil capillarity. Overall, soil texture gave more obvious differences. Specifically, WRCs estimated by BEST returned higher soil water contents for fine soils and lower for coarse soils (h < 100 mm), as expected. However, relatively lower soil water contents were shown close to water saturation for medium soils because of a higher soil bulk density of that class. Conversely, relatively higher discrepancies were detected for HCF because the saturated hydraulic conductivity function changed less abruptly as a function of soil water content for coarse soils than that for medium or fine textured soils ( Figure 3). On the other hand, grouping soils by bulk density or capillarity provided equally plausible findings. In fact, higher θ values were shown close to water saturation (h > 100 mm) for lower BD values (i.e., <1.2 g cm −3 ) and lower θ values for higher BD values (i.e., >1.2 g cm −3 ); this is expected from a physical point of view because more soil porosity determines greater water availability into the soil. However, although the differences in observed BD values were within the range 1.26-1.45 g cm −3 , BEST returned curves that overlapped for h values lower than about 100 mm (Figure 3). Similar results have been obtained for HCF. Intermediate discrepancies in θ values were obtained when investigated soils were grouped by soil capillarity. Also for this comparison, the obtained findings were plausible because (i) for a given h value, the corresponding soil water content was lower for a weak soil capillarity and higher for a strong one; and (ii) a less evident s-shaped curve was obtained as compared with those characterized by moderate-strong λ c values. Negligible effects were detected for HCF (Figure 3).

Evaluation of BEST for Soil Properties Classes
Tables 2-4 summarize the results of capacitive-based soil indicators grouped by texture classes (Table 2), bulk density (Table 3), and soil capillarity (Table 4). Overall, regardless of the classification adopted, results were in agreement with the theoretical assumptions. For instance, AC increased passing from finer to coarser soil textures, while it decreased as the bulk density of the soil, or its capillarity, increased. Conversely, RFC showed the opposite trend because lower, or higher, RFC values accounted for a greater, or lower, air presence into the soil. Although PAWC results appeared to be consistent with the expectations, for example PAWC increased with increasing capillarity of the soil (higher soil capillarity is expected for finer soils, such as clayey soils), it showed relatively small variations, and a clear schematization was not obtained in terms of texture and BD classes. This was attributed to the relatively low heterogeneous dataset available or to the different sample sizes within each class. However, despite the mentioned approximations, the preliminary analysis provided further evidence of the applicability of BEST to assess the effects of soil management on soil property modifications.  Tables 2-4 summarize the results of capacitive-based soil indicators grouped by texture classes (Table 2), bulk density (Table 3), and soil capillarity (Table 4). Overall, regardless of the classification adopted, results were in agreement with the theoretical assumptions. For instance, AC increased passing from finer to coarser soil textures, while it decreased as the bulk density of the soil, or its capillarity, increased. Conversely, RFC showed the opposite trend because lower, or higher, RFC values accounted for a greater, or lower, air presence into the soil. Although PAWC results appeared to be consistent with the expectations, for example PAWC increased with increasing capillarity of the soil (higher soil capillarity is expected for finer soils, such as clayey soils), it showed relatively small variations, and a clear schematization was not obtained in terms of texture and BD classes. This was attributed to the relatively low heterogeneous dataset available or to the different sample weak moderate strong f) Figure 3. Soil water retention curve and hydraulic conductivity function (left and right, respectively) obtained for a given class of soil texture (fine, medium, and coarse; (a,b)), soil bulk density range (lower than 1.2 g cm −3 , between 1.2 and 1.4 g cm −3 , and higher than 1.4 g cm −3 ; (c,d)), or soil capillarity class (10 < λ c < 42 mm, 42 ≤ λ c ≤ 250 mm, and 250 ≤ λ c ≤ 1000 mm representing weak, moderate, and strong, respectively; (e,f)).  The comparison between θ values predicted by BEST and the corresponding measurements obtained in lab are depicted in Figure 4. Overall, BEST's ability to estimate the soil water retention was comparable between the two sites because the root-mean-square deviations were equal to 0.068 and 0.072 cm 3 cm −3 , respectively, for Neviano and Sternatia. This suggests that, regardless of the soil texture (sandy loam or sandy clay loam) or the soil management (no tillage or minimum tillage), the accuracy degree of the θ estimations was comparable. However, some discrepancies were detected as a function of the applied soil pressure head (h) because BEST seems to have been more accurate near to water saturation for ST and more accurate under unsaturated soil conditions for NE ( Figure 4). In this regard, in order to estimate the capacitive-based soil indicators, relatively small differences have been identified at saturation or at field capacity, while higher discrepancies were detected at the wilting point. As a consequence, a relatively higher reliability is expected for AC and RFC rather than PAWC.  Finally, for a given site, the correlation between the BEST estimations (i.e., AC and RFC) and the measurements independently obtained (BD) were consistent with expectations [4], and they provide further evidence on the good reliability of BEST in providing consistent information on the relationships between soil properties ( Figure 5).

Impact of Soil Tillage on Physical and Hydraulic Properties of the Soil
Overall, the lack of alternative soil tillage at Andrano, Gagliano del Capo, and Soleto sites prevented to establish comparisons. Consequently, only the sites of Neviano and Sternatia were  Finally, for a given site, the correlation between the BEST estimations (i.e., AC and RFC) and the measurements independently obtained (BD) were consistent with expectations [4], and they provide further evidence on the good reliability of BEST in providing consistent information on the relationships between soil properties ( Figure 5). relationships between soil properties ( Figure 5).

Impact of Soil Tillage on Physical and Hydraulic Properties of the Soil
Overall, the lack of alternative soil tillage at Andrano, Gagliano del Capo, and Soleto sites prevented to establish comparisons. Consequently, only the sites of Neviano and Sternatia were considered for further investigation. Results of Neviano are reported in Figure 6. Although the no-tilled sites of Neviano were very close (about 50 m), so we assumed they could have the same soil texture, discrepancies were detected in the particle size distribution, as they were sandy-loam and clay soils. Based on available information, we established that only one of the two sites (i.e., clay soil) was undisturbed for a long time, while the second one (sandy-loam soil) was periodically tilled until sometime before (4-5 years or more). This finding is plausible because, considering that the bedrock was quite superficial in the investigated soils, continuous soil tillage (with or without inversion of the soil layers) may have brought back relatively coarser soil particles to the surface, or determined the downward migration of finer ones. This would explain the highest sand (about 69%) and the lowest clay (about 10%) contents observed at the NE site (Table 1). However, NE sites were comparable in terms of BD, that is, they were not statistically different according to Manici et al. [4], and the probability that the soil

Impact of Soil Tillage on Physical and Hydraulic Properties of the Soil
Overall, the lack of alternative soil tillage at Andrano, Gagliano del Capo, and Soleto sites prevented to establish comparisons. Consequently, only the sites of Neviano and Sternatia were considered for further investigation.
Results of Neviano are reported in Figure 6. Although the no-tilled sites of Neviano were very close (about 50 m), so we assumed they could have the same soil texture, discrepancies were detected in the particle size distribution, as they were sandy-loam and clay soils. Based on available information, we established that only one of the two sites (i.e., clay soil) was undisturbed for a long time, while the second one (sandy-loam soil) was periodically tilled until sometime before (4-5 years or more). This finding is plausible because, considering that the bedrock was quite superficial in the investigated soils, continuous soil tillage (with or without inversion of the soil layers) may have brought back relatively coarser soil particles to the surface, or determined the downward migration of finer ones. This would explain the highest sand (about 69%) and the lowest clay (about 10%) contents observed at the NE site (Table 1). However, NE sites were comparable in terms of BD, that is, they were not statistically different according to Manici et al. [4], and the probability that the soil was undisturbed for at least five years therefore appears to be founded. As a consequence, observed discrepancies in particle size distribution at Neviano prevented a homogeneous comparison in terms of physical and hydraulic properties of the soil. Specifically, the coarser soil showed BD values higher by a factor 1.1 than the finer one (1.26 and 1.15 g cm −3 , respectively for sandy-loam and clay); this resulted in a higher saturated soil water content for the former soil ( Figure 6). As a consequence, coarser soil was more conductive than the finer one by a factor 1.3 (Ks equal to 0.028 and 0.021 mm s −1 ). The shapes of the water retention curves were quite typical, as WRCs of coarser soils are generally less flat as compared to those of finer soils (in other words, the changes of water content as a function of pressure head are relatively larger; therefore, the curve is steeper). Conversely, relatively negligible differences were detected for HCF, with differences in unsaturated hydraulic conductivity, K, by a factor of 2.0 (NE-C/NE-SL) at field capacity (a factor 4.5 when an agronomic field capacity value, i.e., h = 300 cm, was considered).
Sternatia sites showed relatively lower discrepancies (Figure 7). Although ST-MT was sampled a few days after a minimum tillage, comparable BD values were detected (1.18 and 1.20 g cm −3 , respectively, under MT and NT). Therefore, a possible soil compaction during sampling may have occurred. However, the WRC modelled by BEST also takes into account soil texture and infiltrometric measurements. For this site, in fact, differences in Ks values by a factor 2.3 were detected (0.051 and 0.117 mm s −1 , respectively, for NT and MT). Consequently, discrepancies between soil management in WRC and HCF increased from saturated to unsaturated soil conditions (Figure 7). In agreement with the results discussed so far, HCF showed even smaller K discrepancies, equal to a factor 1.3 (ST-NT/ST-MT) at field capacity (a factor 1.6 at h = 300 cm). detected (0.051 and 0.117 mm s , respectively, for NT and MT). Consequently, discrepancies between soil management in WRC and HCF increased from saturated to unsaturated soil conditions (Figure 7). In agreement with the results discussed so far, HCF showed even smaller K discrepancies, equal to a factor 1.3 (ST-NT/ST-MT) at field capacity (a factor 1.6 at h=300 cm).
In summary, the differences in water retained in the two soils compared were not negligible (up to 18%), both at field capacity and wilting point (Figure 8). However, results of the homogeneous comparison in terms of soil texture carried out at ST, although they have quantified the maximum discrepancies due to alternative soil tillage (CT and NT), also confirms that optimal soil aeration can be reached with the natural restoration of the soil structure.   In summary, the differences in water retained in the two soils compared were not negligible (up to 18%), both at field capacity and wilting point ( Figure 8). However, results of the homogeneous comparison in terms of soil texture carried out at ST, although they have quantified the maximum discrepancies due to alternative soil tillage (CT and NT), also confirms that optimal soil aeration can be reached with the natural restoration of the soil structure.
The effects of soil management on capacitive-based soil indicators are depicted in Figure 10.

Discussion
The optimization of water resources is an open topic in the agriculture of the Mediterranean basin [42], and rational soil management of rainfed orchards needs the knowledge of the hydraulic functions of the soil [43]. For instance, Rallo et al. [44] investigated the relationship between plant and soil water status for mature olive orchards in Sicily using the AQUACROP model to establish critical thresholds of soil water content, in accordance with the vegetative stages. Moreover, in order to apply the selected model on physically based bases, the physical and hydraulic properties of the soil were measured directly using standard methods, and the model parameters were estimated accordingly. Findings allowed to establish two critical thresholds of soil water content at specific soil pressure head values: -40 and -200 m, equal respectively to 16% and 11% of volumetric soil water content. This result has practical implications for crop production because θ values higher than 16% can assure the absence of water stress, while for lower values (θ < 16%), the crop water stress increased with decreasing soil water contents. In another investigation by Autovino et al. [45], the Figure 8. Differences in soil water retention (∆θ) and in hydraulic conductivity (∆K) as a function of soil pressure head (h) for Neviano (differences between clay and sandy-loam) and Sternatia (differences between no-tillage and minimum tillage). Vertical bars represent θ values at field capacity and wilting point.

Discussion
The optimization of water resources is an open topic in the agriculture of the Mediterranean basin [42], and rational soil management of rainfed orchards needs the knowledge of the hydraulic functions of the soil [43]. For instance, Rallo et al. [44] investigated the relationship between plant and soil water status for mature olive orchards in Sicily using the AQUACROP model to establish critical thresholds of soil water content, in accordance with the vegetative stages. Moreover, in order to apply the selected model on physically based bases, the physical and hydraulic properties of the soil were measured directly using standard methods, and the model parameters were estimated accordingly. Findings allowed to establish two critical thresholds of soil water content at specific soil pressure head values: -40 and -200 m, equal respectively to 16% and 11% of volumetric soil water content. This result has practical implications for crop production because θ values higher than 16% can assure the absence of water stress, while for lower values (θ < 16%), the crop water stress increased with decreasing soil water contents. In another investigation by Autovino et al. [45], the performance of the physically based hydrological model HYDRUS-2D to predict the soil water contents and transpiration fluxes in a Sicilian irrigated olive orchard was assessed. Also in this investigation, accurate soil characterization has been carried out, and the BEST-procedure, coupled with the tension infiltrometer method, was applied to estimate the hydraulic conductivity of the soil. Main results showed that HYDRUS-2D, in general, was able to reproduce the trends of measured soil water contents at the different distances and depths from the plant row, and, in particular, it was also suitable (i) to estimate actual transpiration, (ii) to identify crop water status during the different stages of crop growth, and therefore (iii) can be used to identify irrigation strategies aimed to cope with water scarcity [45]. As a consequence, the cited studies confirm that the knowledge of the physical and hydraulic properties of the soil is the prerequisite for an adequate modeling of the soil-plant-atmosphere (SPA) system.
In this study, the impact of the soil tillage on physical and hydraulic soil properties was investigated by applying the BEST-procedure to obtain complete soil hydraulic characterization. A preliminary comparative evaluation of BEST among investigated soils allowed establishing that the procedure is applicable for the purposes of comparisons in rainfed orchards.
Specific insights for the soils of Neviano and Sternatia investigated the impact of soil texture and soil tillage, respectively, on soil hydraulic functions. Notably, results confirmed that differences in soil texture determined greater differences as compared to soil tillage. However, findings of this study quantified how the voids of the soil, mainly determined by the different particle size distributions (meso-micropores) or by the induced changes in pore size distribution during tillage (macropores), can modify the balance of water and air into the soil. Consequently, since a lack of information on this topic exists for the soils of Salento, this case study can represent a contribution of knowledge to better understand the soil tillage impact on water saving in dry-farmed olive orchards. Figure 9a shows the comparison between the retention curves of the investigated soils obtained by BEST and the "optimal" one, as suggested by Reynolds et al. [35]. Among those considered in the study, two curves (i.e., NE-SL and ST-NT) differed more than the literature reference, and the differences in water retention between estimated and optimal values (Figure 9b) may be summarized, on average, by the sequence NE-SL > ST-NT > SO > NE-C > GC > AN > ST-MT. Therefore, based on this criterion, the main result of this analysis suggests that, compared with no tillage, minimum tillage improved the soil water retention, as to obtain a curve with slightly greater degrees of saturation than that of the optimal curve (i.e., for pore sizes up to 1 µm, and lower ones for bigger soil pores).
On the other hand, the comparison between NT and MT at Sternatia highlighted the possible RFC range between a long-term no-tilled soil (ST-NT) and a recently tilled one (ST-MT) ( Figure 10). For this second case study, our results showed that ST-NT soil had optimal average RFC values, and that surface soil tillage (MT) induced an obvious increase in soil porosity. However, this relatively negative soil condition (RFC < 0.6) should be considered as temporary because such optimal RFC values (0.6-0.7) are theoretically achievable again, in accordance with the time for a natural re-consolidation of the soil. Moreover, MT is traditionally suggested in dry-farming as useful practice to break the continuity of the soil capillaries and decrease water losses by evaporation or to control the spread of weeds (both in conventional and biological agriculture). Moreover, according with the current European legislation, this option should ensure to control disease vectors in southern Italy [4,46]. Among capacitive-based soil indicators considered in this investigation, references of literature suggest that multivariate analysis has identified the relative field capacity (RFC) as a key soil physical indicator, as it partially combines the air capacity (AC) and plant available water capacity (PAWC) of the soil, thus expressing the optimal air/water ratio into the porous medium [31]. As a consequence, RFC was selected to compare the investigated olive orchards. Figure 10 clearly shows that, compared to NE-C, NE-SL showed an imbalance of the soil properties towards excessive aeration (RFC < 0.6), with potential negative effects such as the reduction in microbial activity and nitrate production because of insufficient soil water [30]. This result fully summarizes the possible effect of the different soil textures on the optimal water-air availability in the soil. Overall, several investigations in the Mediterranean basin have evaluated the impact of surface tillage (i.e., minimum or conventional) on soil erosion [47,48], pointing out that the shift to organic farming (i.e., reduction of soil tillage; spread of agronomic management systems that incorporate a vegetative cover) could assure an increased protection of the soil and a corresponding lowered erosion risk. However, Salento soils are generally flat, and under the experimental conditions of the Sternatia site, ST-CT has shown a lower saturated hydraulic conductivity than that of ST-MT by a factor of 2.3. This should allow a better wetting of the soil profile even for low rainy events (lower water losses due to deep percolation) and, consequently, the improvement of the soil water storage.

Conclusions
This study represents a contribution of knowledge for water resource optimization of the olive orchards of Salento, and the impact of minimum and no tillage on soil water retention was quantified through well-established guidelines of literature.
The BEST-procedure was able to model the hydraulic functions of the soil through a simple and relatively quick procedure, with an acceptable degree of accuracy, and it was able to provide indicators to quantify the impact of soil use on water saving. Therefore, its reliability was established for agro-environments different from those so far tested.
The positive effects of the minimum tillage, as compared to a long-term no-tilled soil, were clearly highlighted, and comparisons with optimal values in the literature strengthened the experimental information reached. Consequently, although periodic surface tillage seems an advisable solution to improve the water reserve of the soils, further investigations are necessary to verify such findings in order to improve the sustainability of the arid Mediterranean environments. On the other hand, the comparison between NT and MT at Sternatia highlighted the possible RFC range between a long-term no-tilled soil (ST-NT) and a recently tilled one (ST-MT) ( Figure 10). For this second case study, our results showed that ST-NT soil had optimal average RFC values, and that surface soil tillage (MT) induced an obvious increase in soil porosity. However, this relatively negative soil condition (RFC < 0.6) should be considered as temporary because such optimal RFC values (0.6-0.7) are theoretically achievable again, in accordance with the time for a natural re-consolidation of the soil. Moreover, MT is traditionally suggested in dry-farming as useful practice to break the continuity of the soil capillaries and decrease water losses by evaporation or to control the spread of weeds (both in conventional and biological agriculture). Moreover, according with the current European legislation, this option should ensure to control disease vectors in southern Italy [4,46].
Overall, several investigations in the Mediterranean basin have evaluated the impact of surface tillage (i.e., minimum or conventional) on soil erosion [47,48], pointing out that the shift to organic farming (i.e., reduction of soil tillage; spread of agronomic management systems that incorporate a vegetative cover) could assure an increased protection of the soil and a corresponding lowered erosion risk. However, Salento soils are generally flat, and under the experimental conditions of the Sternatia site, ST-CT has shown a lower saturated hydraulic conductivity than that of ST-MT by a factor of 2.3. This should allow a better wetting of the soil profile even for low rainy events (lower water losses due to deep percolation) and, consequently, the improvement of the soil water storage.

Conclusions
This study represents a contribution of knowledge for water resource optimization of the olive orchards of Salento, and the impact of minimum and no tillage on soil water retention was quantified through well-established guidelines of literature.
The BEST-procedure was able to model the hydraulic functions of the soil through a simple and relatively quick procedure, with an acceptable degree of accuracy, and it was able to provide indicators to quantify the impact of soil use on water saving. Therefore, its reliability was established for agro-environments different from those so far tested.
The positive effects of the minimum tillage, as compared to a long-term no-tilled soil, were clearly highlighted, and comparisons with optimal values in the literature strengthened the experimental information reached. Consequently, although periodic surface tillage seems an advisable solution to improve the water reserve of the soils, further investigations are necessary to verify such findings in order to improve the sustainability of the arid Mediterranean environments.