Short- and Medium-Term Effects of On-Farm Compost Addition on the Physical and Hydraulic Properties of a Clay Soil

: Compost amendment is a widespread agronomic practice, but little information is available about the short- and medium-term effects on clay soils. In this investigation, we selected three soil compost rates (treatments, T), i.e., fertilizer (T2 = 1.5 kg m − 2 ), amending (T3 = 15 kg m − 2 ), and organic (T4 = 75 kg m − 2 ), that were compared to a control (T1). Our research accounts for the effects reached on representative large boxes (about 0.75 m 3 ), without crops and for about eighteen months. An overall assessment of the physical and hydraulic properties of the soil was made, including inﬁltration rate or saturated hydraulic conductivity ( K s ), water content and water retention ( θ ) of the soil, bulk density and structure, and several physical quality indicators obtained from the water retention curve, accounting for the optimal balance between water/air into the soil, pore volume distribution function, and soil features in the inﬂection point of the soil water retention curve. Additionally, the temporal changes of K s were evaluated. The main results showed that (i) after eighteen months, and regardless of T, θ signiﬁcantly improved by a factor of 1.2–1.3, but these improvements (up to 1%) were detected only close to water saturation (i.e., until 6 cm of soil pressure head) when little (T2) or no compost (T1) addition was considered, while a larger range was detected (until 60 cm) when higher rates (T3–T4) were used; (ii) K s determination allowed to establish that compost effects vanished after about eleven months, but it was not possible to verify that composting increased the permeability of the investigated clay soil within that time frame; (iii) some signiﬁcant correlations between K s and some soil physical indicators estimated from both the inﬂection point of the water retention curve or bulk density suggested possible improvements in soil permeability. Because some factors (water content of the soil above all) could have affected the comparison of K s measurements, further research on this topic is needed.


Introduction
Intensive agriculture can lead to soil degradation, and the increase in world food demand is expected to accentuate this issue [1]. One of the main causes of soil degradation is loss of soil fertility due to the progressive organic matter depletion of soil medium [1,2]. Therefore, with a decrease in structural porosity, loss of soil structure is expected, which can lead to a worsening of hydrodynamic soil properties [3,4]. Consequently, farm cropping systems should take account of such risk and manage, or stabilize, soil organic matter levels over time.
Horticultural crops typically produce a relatively large amount of organic waste due to the difference between marketable yields and plant residues. The on-farm composting of organic waste has been found to be environmentally sustainable in many investigations, both in terms of energy consumption and carbon emissions [5][6][7]. Therefore, using zerokilometer compost may represent an agronomic practice to be suggested, especially when combined with other agroecological techniques [7,8]. and north-west. Because of heavy rainfall events, frequent flooding can occur with a rapid increase in groundwater levels.
The soil texture of the upper layer of the experimental farm (0-0.5 m) was a clay with 58.8 and 35.8% of clay and silt, respectively [19]. Further information on the study site can be found in Persiani et al. [7].
At the beginning of autumn 2019, four plastic boxes of about 0.67 m 3 (1.1 m × 1.1 m × 0.55 m) were filled with air-dried soil collected on the farm, and four different soil and compost mixtures (i.e., soil treatments, Ts) were considered. To do this, about four cubic meters of soil were collected from the MITIORG field experiment on the experimental farm. This is an organic long-term experiment (started in 2014), where different combined agroecological practices are tested in horticultural rotations to study the adaptation of the vegetable systems to extreme rainfall and flooding events [8].
To assess the impact of compost addition on the physical and hydraulic properties of the soil, an on-farm compost produced in a small-scale composting plant setup at the experimental farm was used. Crop residues used for the composting process were: vicia faba minor L., plants and fruits of Cucurbita pepo L., Festuca arundinacea Schreb/Trifolium repens L. (ratio 50/50%), Festuca arundinacea Schreb/Trifolium repens L. (ratio 70/30%), and compost as a starter. Three compost-soil mixtures, i.e., three treatments (T2 to T4), were compared with a control without compost addition (T1). In particular, a fertilizer (T2, equal to 1.5 kg m −2 ), amending (T3, equal to 15 kg m −2 ), and organic (T4, 75 kg m −2 ) dose were considered. Practically, after filling the four boxes with untreated (unamended) soil, the same soil of the three boxes (i.e., T2 to T4) corresponding to a depth of 30 cm deep was moved to the work surface, where it was manually mixed to different air-dried compost amounts, namely 2 kg (T2), 18 kg (T3), and 91 kg (T4). Boxes were then filled again with the obtained soil mixtures, and no particular layer-by-layer compaction procedure was adopted, in order to reproduce relatively high soil porosity conditions. During box filling, holes were created on one side of each box with a drill for the (horizontal) insertion of 10 cm probes for measuring the soil water content (θ p ) at two depths, i.e., 10 and 20 cm, from the soil surface. Therefore, an unmonitored buffer soil layer (i.e., between 20 and 30 cm) was created over the untreated soil.
Measurements were acquired and recorded in a data logger, and the on-farm weather station provided the hourly rainfall. An image of the experimental device is reported as Supplementary Material ( Figure S1). Monitoring of soil water content and rainfall began at the end of November 2019 and ended at the end of May 2021.
The main information on the compost used in this investigation is summarized in Table 1. Further details on the composting process can be found in Persiani et al. [7]. Table 1. Days of composting process, and mean values (and standard deviations) of dry matter, pH, electrical conductivity (EC), total nitrogen (N) content, total organic carbon (C) content, and C/N ratio of the compost used in the investigation.

Soil Measurements
Soil physical properties, i.e., water retention and bulk density, were determined both at the beginning and at the end of the investigation. Specifically, to characterize the soil water retention curve, i.e., the relationship between volumetric water content (θ) and soil pressure head (h), and soil dry bulk density (BD), two soil cores (5 cm in height by 8 cm in diameter, equal to 251 cm 3 ) were collected in the surface soil layer of each box at the beginning of the investigation (i.e., in January 2020, two months after box preparation). This relatively small sample size was considered adequate given the lack of structure of Agronomy 2022, 12, 1446 4 of 16 the disturbed soil considered. After core sampling, more soil was poured into the holes to restore the flat soil surface. A similar soil sampling was carried out at the end of the investigation (after 20 months), and four soil cores were collected for each soil treatment.
Soil water retention was determined by the tension hanging water column apparatus [20] for pressure head h (cm) values ranging from −5 to −100 cm and the pressure plate extractors [21] for h values ranging from −1030 to −15,300 cm. Following a standard methodology, the soil cores were first saturated within the buchners (funnels) by capillary rising. After saturation, a drainage process was imposed, and seven h values (−5, −10, −20, −40, −60, −80, and −100 cm) were set to determine the corresponding volumetric soil water content (θ) values [22]. In addition, three other θ values were determined on the repacked soil samples at low h values (−1030, −3060, and −15,300 cm). BD was determined at the end of the tension experiments (i.e., at h = −100 cm) assuming a particle density of 2.65 g cm -3 [10,22]. The pairs of θ-h values were fitted by the van Genuchten model (vG) using the SWRC Fit (web interface) by Seki [23], and the soil water retention curve, i.e., the θ(h) function, was obtained.
Saturated soil hydraulic conductivity (K s ) was determined from early October 2020 (i.e., eleven months after box filling) to mid-May 2021 by the SSBI method [24] over six dates (hereafter referred to as sampling dates (SDs)). For each box, six infiltrations were carried out on the soil surface. Consequently, 144 (4 treatments × 6 sampling dates × 6 replicates) infiltrations were carried out in this investigation.
The SSBI method involves carrying out infiltration tests of the Beerkan type. In detail, Beerkan tests were performed using an 8 cm inner-diameter steel cylinder, inserted to a depth of about 1 cm to avoid lateral loss of ponded water. A known volume of water (55 mL) was repeatedly poured into the cylinder, establishing a height of water of about 1 cm, and the time needed for complete infiltration was logged. The procedure was repeated for a total of 18 water volumes, and experimental cumulative infiltration, I(t), and infiltration rate, i r , were thus deduced. In accordance with Bagarello et al. [24], saturated soil hydraulic conductivity was deduced by the equation: where i s (L T −1 ) is the slope of the straight line defined by the last part of the cumulative infiltration describing the steady-state conditions, r (L) is the radius of the ring, 1.3635 is a coefficient accounting for the shape of the wetting-drainage front and the geometrical correction of the infiltration front shape [24], and α* (L −1 ) is a parameter expressing the relative importance of gravity and capillary forces during the infiltration process [25]. A value of α* equal to 4 m -1 , suggested for soils with strong capillarity (i.e., porous materials that are both fine textured and massive, including unstructured clayey and silty soils, as well as very fine to fine structureless sandy materials), was used in this investigation. For further details, please refer to the original manuscript by Bagarello et al. [24], while for a brief overview of current applications, see the review by Castellini et al. [26]. Total organic carbon (TOC), determined only at the end of the investigation, was quantified through the dry combustion method using a TOC Vario Select analyzer (Elementar, Germany) [27].

Soil Physical Quality Determination
For each treatment, averaged retention curves (and bulk density values) were used to estimate several soil physical quality indicators accounting for (i) soil compaction or balance between the liquid/aeriform phase of the soil (Group 1); (ii) the pore volume distribution function (Group 2); and (iii) the inflection point of the "S-shaped" water retention function, i.e., the point that discriminates the non-capillary (macropores) from capillary porosity (Group 3). Specifically, Group 1 includes BD and four widely used soil indicators such as macroporosity (P MAC ), air capacity (AC), relative field capacity (RFC), and plant available water capacity (PAWC) [13]; Group 2 includes the "location parameters", i.e., modal (d mod ), median (d med ) and mean (d m ) pore diameter, and the "shape parameters", i.e., standard deviation (Sd), skewness (Sk), and kurtosis (Ku) of the pore volume distribution function [13]; Group 3 includes the soil pressure head at the inflection point of the water retention curve (h inf ), the corresponding volumetric soil water content (θ inf ), soil porosity at h inf (POR inf ), S-index by Dexter [28], namely the slope of the gravimetric water retention curve at the inflection, and the pore distribution index, λ inf by Han et al. [29] (i.e., the ratio between S-index and θ inf ). According to the procedure indicated by Reynolds et al. [13], the pore volume distribution functions of the considered treatments were calculated. These waterretention-curve-derived indicators are widely used and suggested in the literature [30,31], also to verify plausible correlations with respect to independent variables (for example, K s ) and check the consistency of the results [32,33].

Data Analysis
For each main variable considered in this investigation, i.e., BD, θ-h, θ i , P MAC , AC, RFC, PAWC, d mod , d med , h inf , θ inf , POR inf , S-index, and λ inf , arithmetic means and corresponding coefficients of variation were always calculated, with the exception of K s , d m , Sd, Sk, and Ku because of log-normal distributions, and, therefore, geometric means and associated coefficients of variation are generally suggested for these soil properties [34].
A Tukey HSD (honestly significant difference) test was used to evaluate temporal changes by comparing mean values and to establish multiple comparisons among treatments [35], whereas a two-tailed paired t-test was used to check the statistical significance of the soil water retention data. A probability level of p = 0.05 was assumed.
Regression analysis was carried out between saturated hydraulic conductivity and independent soil variables, namely soil physical quality indicators obtained from the soil water retention curve or bulk density. With this aim, the full set of data (i.e., soil physical quality indicators obtained after early and final soil sampling) or partial set of data (only after early or final soil sampling) were correlated with corresponding K s values, and only the correlations that were both significant and expected from a physical point of view were considered for comparison between treatments [32,33].

Impact on Water Content, Water Retention, and Bulk Density of the Soil
The temporal changes in mean soil water content (SWC) corresponding to the four treatments (T1 to T4) and the hourly rainfall events are reported in Figure 1. Except for the first months (September-October 2020), the soil water content was generally always higher than the field capacity and lower than the full saturation. Overall, an inadequate soil-probes contact probably caused SWC underestimation in Treatments 1 and 4 (T1 and T4) at the beginning of monitoring, and a probe malfunction did not allow to obtain reliable measurements from March onwards in T2 and T4. Therefore, a full dataset for specific analysis (correlations) was obtained for only four out of six sampling dates (i.e., SD1 to SD4). Both measured field capacity and estimated saturated soil water contents account for the soil condition at the end of the investigation, i.e., after a plausible complete soil structure restoration, because they were determined at the end of the research. In particular, mean values of field capacity of about 0.30 or 0.32 cm 3 cm −3 were detected for T1 and T2, or T3 and T4, while estimated saturated soil water content increased with increasing compost concentration, ranging from 0.53 cm 3 cm −3 (T1) to 0.60 cm 3 cm −3 (T4) (Figure 1). The comparison between treatments in terms of soil water retention in the early (2nd month) and final (20th month, after SD6) stages of the investigation is reported in Figure  2. Overall, a quite clear induced effect of compost addition (i.e., addition of increasing doses of compost) on soil water retention capacity improvement was always detected (i.e., in both stages). However, with reference to the final phase, a similarity between pairs of curves was identified (T1-T2 and T3-T4), and some improvements in soil water retention (up to 1%) were detected only close to water saturation (i.e., range of soil pressure heads of about 0-6 cm) when little (T2) or no compost (T1) addition was considered, while a larger range was detected (0-60 cm) when higher rates (T3-T4) were used ( Figure 2). As expected, significantly higher soil water contents were always detected at the end of the experiment, according to the two-tailed paired t-test (p = 0.05). Regardless of the stages of the investigation (early or final), no particular correlation between van   The comparison between treatments in terms of soil water retention in the early (2nd month) and final (20th month, after SD6) stages of the investigation is reported in Figure 2. Overall, a quite clear induced effect of compost addition (i.e., addition of increasing doses of compost) on soil water retention capacity improvement was always detected (i.e., in both stages). However, with reference to the final phase, a similarity between pairs of curves was identified (T1-T2 and T3-T4), and some improvements in soil water retention (up to 1%) were detected only close to water saturation (i.e., range of soil pressure heads of about 0-6 cm) when little (T2) or no compost (T1) addition was considered, while a larger range was detected (0-60 cm) when higher rates (T3-T4) were used ( Figure 2). The comparison between treatments in terms of soil water retention in the early (2nd month) and final (20th month, after SD6) stages of the investigation is reported in Figure  2. Overall, a quite clear induced effect of compost addition (i.e., addition of increasing doses of compost) on soil water retention capacity improvement was always detected (i.e., in both stages). However, with reference to the final phase, a similarity between pairs of curves was identified (T1-T2 and T3-T4), and some improvements in soil water retention (up to 1%) were detected only close to water saturation (i.e., range of soil pressure heads of about 0-6 cm) when little (T2) or no compost (T1) addition was considered, while a larger range was detected (0-60 cm) when higher rates (T3-T4) were used ( Figure 2). As expected, significantly higher soil water contents were always detected at the end of the experiment, according to the two-tailed paired t-test (p = 0.05). Regardless of the stages of the investigation (early or final), no particular correlation between van   As expected, significantly higher soil water contents were always detected at the end of the experiment, according to the two-tailed paired t-test (p = 0.05). Regardless of the stages of the investigation (early or final), no particular correlation between van Genuchten model parameters (θ r , α, and n) of the soil water retention curve and compost rate was identified. Conversely, the expected increasing relationships between the selected main values of the soil water content (θ s , θ 10 , and θ 100 , respectively, equal to the water contents of the saturated soil, matrix, and at field capacity) and compost rate were always verified. Moreover, when similar correlations were checked using some capacitive indicators obtained from the soil retention curve (i.e., aeration capacity (AC), plant available water capacity (PAWC), relative field capacity (RFC)) or soil bulk density (BD), we verified that AC and PAWC increased with the rate, while RFC and BD instead decreased. These results were logical from a soil physics point of view, because it is expected that available water or soil air can increase with the compost rate; likewise, relative field capacity and soil density will decrease.
T4 showed the greatest impact as compared to the remaining treatments. For instance, when we considered BD changes over time, i.e., between the beginning and the end (shortly after SD6) of the investigation, we detected a higher increase in BD values under T1-T3 (by a factor of 1.31-1.36) than T4 (a factor of 1.11); BD values ranged from 0.81 to 1.11 g cm −3 under T1-T3 and 0.77 to 0.86 g cm −3 under T4. Therefore, treatment 4 experienced less compaction. A similar comparison was carried out for soil water retention close to water saturation (range: 5-100 cm), and the results showed that the differences (i.e., improvements) in water retention between the beginning and the end of the investigation ranged between 8-9% (T2-T3) and 11% (T4), as compared to T1 (7%); consequently, an increase between 2 and 4% could be attributed to the effect of compost addition.

Impact on Infiltration Rate and Saturated Hydraulic Conductivity of the Soil
The infiltration rate curves, plotted against cumulative infiltration in order to simplify comparisons between treatments, revealed a curve overlap from SD3 onwards ( Figure 3). However, detectable differences were obtained at the beginning of the trial with the reference soil (control), T1, which showed both the lowest (SD1) and the highest (SD2) infiltration rates, when compared with the composted soils ( Figure 3). Therefore, the impact of the amendment waned over time. Genuchten model parameters (θr, α, and n) of the soil water retention curve and compost rate was identified. Conversely, the expected increasing relationships between the selected main values of the soil water content (θs, θ10, and θ100, respectively, equal to the water contents of the saturated soil, matrix, and at field capacity) and compost rate were always verified. Moreover, when similar correlations were checked using some capacitive indicators obtained from the soil retention curve (i.e., aeration capacity (AC), plant available water capacity (PAWC), relative field capacity (RFC)) or soil bulk density (BD), we verified that AC and PAWC increased with the rate, while RFC and BD instead decreased. These results were logical from a soil physics point of view, because it is expected that available water or soil air can increase with the compost rate; likewise, relative field capacity and soil density will decrease. T4 showed the greatest impact as compared to the remaining treatments. For instance, when we considered BD changes over time, i.e., between the beginning and the end (shortly after SD6) of the investigation, we detected a higher increase in BD values under T1-T3 (by a factor of 1.31-1.36) than T4 (a factor of 1.11); BD values ranged from 0.81 to 1.11 g cm -3 under T1-T3 and 0.77 to 0.86 g cm -3 under T4. Therefore, treatment 4 experienced less compaction. A similar comparison was carried out for soil water retention close to water saturation (range: 5-100 cm), and the results showed that the differences (i.e., improvements) in water retention between the beginning and the end of the investigation ranged between 8-9% (T2-T3) and 11% (T4), as compared to T1 (7%); consequently, an increase between 2 and 4% could be attributed to the effect of compost addition.

Impact on Infiltration Rate and Saturated Hydraulic Conductivity of the Soil
The infiltration rate curves, plotted against cumulative infiltration in order to simplify comparisons between treatments, revealed a curve overlap from SD3 onwards ( Figure 3). However, detectable differences were obtained at the beginning of the trial with the reference soil (control), T1, which showed both the lowest (SD1) and the highest (SD2) infiltration rates, when compared with the composted soils ( Figure 3). Therefore, the impact of the amendment waned over time.    Figure 4 shows an illustrative example of the equilibration time (t s ) and infiltrated depth at the equilibration time (I(t s )) determination from cumulative infiltration, using the common criteria by Angulo-Jaramillo et al. [36] for analyzing cumulative infiltration data,  Table 2 summarizes the corresponding results obtained in this investigation. Similar to the information reported, t s and I(t s ) changed mainly in the first two sampling dates, then remained constant and comparable in the following ones. A certain inaccuracy of estimates was obtained for SD2-T2, with values of t s that, on average, were two or three times higher.
Agronomy 2022, 12, x FOR PEER REVIEW 8 of 16 Figure 4 shows an illustrative example of the equilibration time (ts) and infiltrated depth at the equilibration time (I(ts)) determination from cumulative infiltration, using the common criteria by Angulo-Jaramillo et al. [36] for analyzing cumulative infiltration data, while Table 2 summarizes the corresponding results obtained in this investigation. Similar to the information reported, ts and I(ts) changed mainly in the first two sampling dates, then remained constant and comparable in the following ones. A certain inaccuracy of estimates was obtained for SD2-T2, with values of ts that, on average, were two or three times higher.  Box plots of saturated hydraulic conductivity, Ks, reported in Figure 5, show some significant differences only in SD1 and SD2. Overall, a comparable variability among sampling dates was detected, as mean values of the coefficient of variation were within the range of 26-47%. However, an almost double variability was detected in SD3, as a probable effect of the previous rainfall events, which impacted differently on the different treatments. Regarding the first two dates, relatively lower Ks mean values were detected for T1 under SD1 and T2 under SD2 ( Figure 5). According to the Tukey HSD test (p = 0.05), saturated hydraulic conductivity values obtained from the fertilizer dose were  Box plots of saturated hydraulic conductivity, K s , reported in Figure 5, show some significant differences only in SD1 and SD2. Overall, a comparable variability among sampling dates was detected, as mean values of the coefficient of variation were within the range of 26-47%. However, an almost double variability was detected in SD3, as a probable effect of the previous rainfall events, which impacted differently on the different treatments. Regarding the first two dates, relatively lower K s mean values were detected for T1 under SD1 and T2 under SD2 ( Figure 5). According to the Tukey HSD test (p = 0.05), saturated hydraulic conductivity values obtained from the fertilizer dose were significantly higher than control in SD1 and significantly lower in SD2; after that, the small discrepancies observed vanished in the following sampling dates (SD3 to SD6) ( Figure 5). significantly higher than control in SD1 and significantly lower in SD2; after that, the small discrepancies observed vanished in the following sampling dates (SD3 to SD6) ( Figure 5). Although there was a lack of differences between treatments on sampling dates 3 to 6, the results showed some temporal changes, especially for T1, T2, and T4. As an example of the extreme treatments, i.e., control and organic dose (T1 and T4), the results of the Tukey HSD test highlighted significantly higher Ks values in SD3 and SD2, respectively, for T1 and T4 (Table 3).  Although there was a lack of differences between treatments on sampling dates 3 to 6, the results showed some temporal changes, especially for T1, T2, and T4. As an example of the extreme treatments, i.e., control and organic dose (T1 and T4), the results of the Tukey HSD test highlighted significantly higher K s values in SD3 and SD2, respectively, for T1 and T4 (Table 3). Adding compost to the soil increases its water content, and this can affect the hydrodynamic properties of the porous medium. Limited to the first four sampling dates, due to the lack of the remaining soil moisture values, the results showed a plausible positive relationship between saturated soil hydraulic conductivity and volumetric soil water content (i.e., K s increases with increasing soil moisture), with a coefficient of determination that increased with increasing compost rate (R 2 = 0.46, 0.62, 0.91 and 0.96). However, only the two highest compost concentrations were found to be significant (T3 and T4), thus suggesting a negligible effect for T1 ( Figure 6). 5 vs. 6 0.900 n.s. 0.900 n.s. 0.900 n.s. 0.522 n.s.
Adding compost to the soil increases its water content, and this can affect the hydrodynamic properties of the porous medium. Limited to the first four sampling dates, due to the lack of the remaining soil moisture values, the results showed a plausible positive relationship between saturated soil hydraulic conductivity and volumetric soil water content (i.e., Ks increases with increasing soil moisture), with a coefficient of determination that increased with increasing compost rate (R 2 = 0.46, 0.62, 0.91 and 0.96). However, only the two highest compost concentrations were found to be significant (T3 and T4), thus suggesting a negligible effect for T1 ( Figure 6).

Impact on Soil Physical Quality
Regardless of the soil treatment, the transition from an unstructured to a more structured soil resulted in an increase in bulk density, relative field capacity, plant available water capacity, soil pressure head at water retention curve inflection point, and corresponding porosity, as well as shape parameters of pore volume distribution function; detected discrepancies ranged by a factor of 1.0 (Ku) to 3.4 (PAWC) ( Table 4). Conversely, soil indicators accounting for the air capacity, macroporosity, or pore diameter, decreased accordingly by a factor ranging from 1.4 (AC) and 3.8 (dm). Moreover, a halving of Dexter's S-index value was detected, since it changed from 0.16 to 0.08 (Table 4). Table 4. Soil physical quality indicators obtained for initial (start) and final (end) stage of investigation, in the treatments T1 to T4. T1: control without compost; T2: fertilizer dose; T3: amending dose; T4: organic dose.

Impact on Soil Physical Quality
Regardless of the soil treatment, the transition from an unstructured to a more structured soil resulted in an increase in bulk density, relative field capacity, plant available water capacity, soil pressure head at water retention curve inflection point, and corresponding porosity, as well as shape parameters of pore volume distribution function; detected discrepancies ranged by a factor of 1.0 (Ku) to 3.4 (PAWC) ( Table 4). Conversely, soil indicators accounting for the air capacity, macroporosity, or pore diameter, decreased accordingly by a factor ranging from 1.4 (AC) and 3.8 (d m ). Moreover, a halving of Dexter's S-index value was detected, since it changed from 0.16 to 0.08 (Table 4).  The comparison between normalized pore volume distributions of considered treatments and the differences between start and end of the investigation is reported as Supplementary Material (Figures S2 and S3, respectively).
The mean values of the seventeen soil physical quality indicators estimated from the water retention curve and bulk density data were correlated with the corresponding K s mean values, and only the significant correlations were evaluated to establish their consistency with theoretical expectations. The overall data set, i.e., characterized by data corresponding to both the beginning and the end of the investigation, never provided significant correlations. However, when the two sub-datasets were considered separately, plausible and significant correlations between variables were detected only for the structured soil (Figure 7). Specifically, when the soil at the start of the investigation (i.e., poorly structured soil) was considered, K s was never significantly correlated with independent soil variables. Conversely, increasing significant correlations were identified for the betterstructured soil, as K s increased with increasing soil bulk density (BD) and modal diameter (d mod ) of the pore volume distribution function; porosity at the inflection point of the water retention curve (POR inf ) increased significantly as well, while consequently, θ inf decreased as K s decreased ( Figure 7). Finally, according to the observed sequence of Tn values around the regression line, although the expected sequence (i.e., T1 towards T4 or vice versa) was not always detected, T4 was always associated with the lowest K s values and corresponding correlated variables, then generally followed by T3 ( Figure 7); an inversion between T1 and T2 also suggests that the lowest compost rate had a negligible effect on the changes in the soil physical and hydraulic properties.
The comparison between normalized pore volume distributions of considered treatments and the differences between start and end of the investigation is reported as supplementary material (Figures S2 and S3, respectively).
The mean values of the seventeen soil physical quality indicators estimated from the water retention curve and bulk density data were correlated with the corresponding Ks mean values, and only the significant correlations were evaluated to establish their consistency with theoretical expectations. The overall data set, i.e., characterized by data corresponding to both the beginning and the end of the investigation, never provided significant correlations. However, when the two sub-datasets were considered separately, plausible and significant correlations between variables were detected only for the structured soil ( Figure 7). Specifically, when the soil at the start of the investigation (i.e., poorly structured soil) was considered, Ks was never significantly correlated with independent soil variables. Conversely, increasing significant correlations were identified for the betterstructured soil, as Ks increased with increasing soil bulk density (BD) and modal diameter (dmod) of the pore volume distribution function; porosity at the inflection point of the water retention curve (PORinf) increased significantly as well, while consequently, θinf decreased as Ks decreased (Figure 7). Finally, according to the observed sequence of Tn values around the regression line, although the expected sequence (i.e., T1 towards T4 or vice versa) was not always detected, T4 was always associated with the lowest Ks values and corresponding correlated variables, then generally followed by T3 ( Figure 7); an inversion between T1 and T2 also suggests that the lowest compost rate had a negligible effect on the changes in the soil physical and hydraulic properties.

Discussion
Soil amendment with compost is a widespread agronomic practice with a view to restore the chemical-physical fertility of soil and, therefore, to counterbalance the risks of soil degradation due to intensive agriculture [9,37,38]. Furthermore, this practice takes on greater environmental significance when compost is produced directly on the farm, starting from the use of waste produced [7]. However, the impacts of this practice on the hydrodynamic properties of clay soils are significant and require experimental checks under real (or representative) conditions.
We investigated the short-and medium-term effects of on-farm compost addition on the physical and hydraulic properties of a clay soil with a relatively low organic fertility.
By making use of some widely applied literature indicators, the investigated reference soil (control without amendment) was structurally degraded, as the structure stability index (SSI) by Pieri [39] was lower than the minimum threshold suggested in the literature (SSI < 5%). Following the suggestions by Prout et al. [2], which relate observed organic matter with clay content, we point out that the investigated soil had about a quarter of the minimum optimal organic level to be defined as "not degraded". TOC levels a little higher than 1%, however, are quite common in Mediterranean agro-environments [31], and for this reason, farmers should be encouraged to progressively improve such soil quality parameters over time.
Compost addition affected the water retention and bulk density of the investigated soil. By way of example, Figure 8 shows a summary of the effects in the final stage of the investigation: the addition of increasing compost doses to the investigated soil obviously increased the total organic carbon content that, consequently, positively affected the air capacity of the soil. The highly significant linear regression obtained provided further evidence of the positive effect of soil amendment on water retention. Overall, our results showed that only the highest amendment (T4) was effective in improving the structural stability of the soil (SSI = 5.4%). Since SSI assessment refers to the end of the investigation, it is possible to state that the highest compost rate (75 kg m −2 ) was adequate to guarantee minimum structural stability for eighteen months at least.

Discussion
Soil amendment with compost is a widespread agronomic practice with a view to restore the chemical-physical fertility of soil and, therefore, to counterbalance the risks of soil degradation due to intensive agriculture [9,37,38]. Furthermore, this practice takes on greater environmental significance when compost is produced directly on the farm, starting from the use of waste produced [7]. However, the impacts of this practice on the hydrodynamic properties of clay soils are significant and require experimental checks under real (or representative) conditions.
We investigated the short-and medium-term effects of on-farm compost addition on the physical and hydraulic properties of a clay soil with a relatively low organic fertility.
By making use of some widely applied literature indicators, the investigated reference soil (control without amendment) was structurally degraded, as the structure stability index (SSI) by Pieri [39] was lower than the minimum threshold suggested in the literature (SSI < 5%). Following the suggestions by Prout et al. [2], which relate observed organic matter with clay content, we point out that the investigated soil had about a quarter of the minimum optimal organic level to be defined as "not degraded". TOC levels a little higher than 1%, however, are quite common in Mediterranean agro-environments [31], and for this reason, farmers should be encouraged to progressively improve such soil quality parameters over time.
Compost addition affected the water retention and bulk density of the investigated soil. By way of example, Figure 8 shows a summary of the effects in the final stage of the investigation: the addition of increasing compost doses to the investigated soil obviously increased the total organic carbon content that, consequently, positively affected the air capacity of the soil. The highly significant linear regression obtained provided further evidence of the positive effect of soil amendment on water retention. Overall, our results showed that only the highest amendment (T4) was effective in improving the structural stability of the soil (SSI = 5.4%). Since SSI assessment refers to the end of the investigation, it is possible to state that the highest compost rate (75 kg m −2 ) was adequate to guarantee minimum structural stability for eighteen months at least. Regardless of the statistical significance detected, the discrepancies in terms of Ks among treatments, at the time of sampling, were generally low and probably negligible from a practical point of view, as the maximum discrepancies among treatments were lower than a factor of 6 in sampling dates 1 and 2; our results also suggest that differences further reduced in the following, as they were never greater than a factor of 2. We hypothesized that (1) the different soil water content at the time of sampling was a key factor Regardless of the statistical significance detected, the discrepancies in terms of K s among treatments, at the time of sampling, were generally low and probably negligible from a practical point of view, as the maximum discrepancies among treatments were lower than a factor of 6 in sampling dates 1 and 2; our results also suggest that differences further reduced in the following, as they were never greater than a factor of 2. We hypothesized that (1) the different soil water content at the time of sampling was a key factor capable of affecting the comparison of K s measurements, and (2) some effects may not have been detected as the first soil measurement (SD1) was carried out approximately eleven months after box filling.
Soil water content at the time of measurements is a main factor affecting the saturated hydraulic conductivity measurements [36]. Our investigation accounted for soil water content that was generally always higher than the field capacity and lower than the water saturation. In this respect, our investigation is placed in a medium-high range of soil water conditions at the time of measurements.
Although it is quite expected that compost may improve (increase) soil water retention (e.g., [37]), temporal changes in hydraulic conductivity (or in infiltration rate) are less obvious and/or easy to interpret overall. It is worth noting that only the highest compost rates, namely amending (T3) and organic (T4) doses, were significantly correlated with saturated hydraulic conductivity, as K s increased with increasing soil moisture values ( Figure 6). However, the results do not allow us to draw general conclusions, as some signals, although not statistically significant, would suggest a relatively greater permeability of unamended soil (i.e., for T1, regression of K s vs. VSWC had the same slope as compared to T3 and T4 but was placed to the left and higher up than the regressions shown in Figure 6).
When K s was correlated with independently obtained soil variables, i.e., BD or with indicators obtained from the water retention curve (i.e., POR inf , θ inf or d mode ), we were also able to obtain information from the relative position of the treatments considered. In this regard, it is worth noting that while the expected positive relationship between compost rate and TOC or AC returns an increasing sequence of treatments (i.e., T1 < T2 < T3 < T4), a different positioning around the regression line is obtained when the effects on conductive features are considered (Figures 7 and 8, respectively), probably due to a negligible effect of the fertilizer dose (T2). The rationale used to interpret our findings, successfully applied in other investigations both in natural [33,40] and agricultural [41] environments, however, was not applicable when the full dataset (initial + ending soil sampling) or when only initial soil condition (initial soil sampling) was considered (Table 4). This result seems obvious enough and may be attributed to the lack of soil structure after box setup. In order to not alter the soil surface of boxes during the investigation, it was not possible to determine water retention at intermediate dates; this prevented establishing after how long a more stable soil structure was reached. However, hydraulic conductivity of saturated soil is, in itself, an indicator that accounts for the soil structure, and checking the short-or medium-term effects of compost addition on soil structure deserves further study. Our results, in fact, highlighted some convincing results, especially those related to the extreme treatments (e.g., the control showed significantly lower K s values in SD1, as compared to subsequent dates), while no temporal changes were detected for the amending dose (T3). Heavy rainfall events may cause compaction of the upper layers of agricultural soils [42], and some studies have recently investigated such effects, highlighting negative effects on both the K s [43] and the physical quality of soil [44]. In our investigation, some rainfall events occurred just before the third sampling date, but their effects could not be quantified. Since implementing organic matter into the soil, or creating a mulching layer, is suggested as a viable solution to mitigate such risk factors, more articulated research is desirable to increase our knowledge on these topics.

Conclusions
The water retention and bulk density of a clay soil were improved by compost implementation at a rate of 75 kg m −2 (organic dose, T4) and 15 kg m −2 (amending dose, T3), while a lower dosage, equal to 1.5 kg m −2 (fertilizer dose, T2), was found to be comparable with the control (soil not amended, T1). The results allowed to verify that some soil indicators obtained from the water retention curve showed the expected trends, because air capacity (AC) and plant available water capacity (PAWC) increased with the compost rate, while relative field capacity (RFC) and soil bulk density (BD) instead decreased. However, only T4 adequately increased soil organic levels, as suggested, for example, by the structural stability index.
The saturated hydraulic conductivity (K s ) determination allowed to unequivocally establish that compost effects vanished after about eleven months. Conversely, it was not possible to state that composting significantly improved (increased) the permeability of the investigated soil. In fact, although some indirect information would suggest that higher compost rates (T3 or T4) may have improved the soil hydrodynamic properties (e.g., K s monotonically increased with increasing soil water contents), the results of the first two sampling dates were uncertain, and, moreover, the main effects probably were lost, as K s measurement began eleven months after box filling. However, significant correlations between K s and the soil variables independently obtained (soil physical quality indicators carried out from the water retention curve or BD) were detected only after a relatively long time (after 20 months, rather than 2 months). Consequently, if the purpose is to evaluate the impact of compost implementation on K s improvement, we could conclude that the time frame is limited to about six months after box filling or, by analogy, after conventional soil tillage.
In conclusion, the results obtained were more convincing for the effects on hydrostatic properties than hydrodynamic ones. The experimental factor that probably contributed to this was the different water content of the soil, as compost is known to help improve soil moisture. Further studies could help improve our knowledge of the effects of compost incorporation on soil structure, hydraulic properties, and physical quality.
Supplementary Materials: The following supporting information can be downloaded at: https:// www.mdpi.com/article/10.3390/agronomy12061446/s1, Figure S1. View of experimental boxes with indication of the four treatments (T) considered. Note that boxes arrangement in figure corresponds to schematization on the side. Figure S2. Normalized pore volume distributions of the four treatments obtained at the start and end of investigation. Figure S3. Comparison of normalized pore volume distributions between start and end of investigation for each considered treatment (T1 to T4). The reference curve by Reynolds et al. [13] was also reported.