Tracking Changes on Soil Structure and Organic Carbon Sequestration after 30 Years of Different Tillage and Management Practices

: Long-term ﬁeld trials are essential for monitoring the effects of sustainable land management strategies for adaptation and mitigation to climate change. The inﬂuence of more than thirty years of different management is analyzed on extensive crops under three tillage systems, conventional tillage (CT), minimum tillage (MT), and no-tillage (NT), and with two crop rotations, monoculture winter-wheat ( Triticum aestivum L.) and wheat-vetch ( Triticum aestivum L.- Vicia sativa L.), widely present in the center of Spain. The soil under NT experienced the largest change in organic carbon (SOC) sequestration, macroaggregate stability, and bulk density. In the MT and NT treatments, SOC content was still increasing after 32 years, being 26.5 and 32.2 Mg ha − 1 , respectively, compared to 20.8 Mg ha − 1 in CT. The SOC stratiﬁcation (ratio of SOC at the topsoil/SOC at the layer underneath), an indicator of soil conservation, increased with decreasing tillage intensity (2.32, 1.36, and 1.01 for NT, MT, and CT respectively). Tillage intensity affected the majority of soil parameters, except the water stable aggregates, inﬁltration, and porosity. The NT treatment increased available water, but only in monocropping. More water was retained at the permanent wilting point in NT treatments, which can be a disadvantage in dry periods of these edaphoclimatic conditions.


Introduction
Herbaceous crops occupy most of the agricultural surface of the planet, and the choice of soil management practice can significantly modify the quality and physical properties of soil, one of which is the soil organic carbon (SOC) content [1][2][3][4], particularly in the context of adaptation to climate change and considering that the accumulation of SOC is a positive indicator of soil conservation. When exploring the maximum potential for SOC sequestration, site-specific conditions must be considered as SOC dynamics vary depending on the edaphic [5] and geographical context. Temperate soils are especially sensitive to the SOC losses associated with some agricultural management systems. Tillage practices can influence SOC storage and dynamics, which is important for both agricultural productivity and carbon sequestration [6,7].
Agricultural soils are usually managed using the 'conventional tillage' (CT) method, based on practices that use the moldboard as the main agricultural implement. However, CT is unsustainable, as it causes reductions in SOC content, destroys soil structure [8,9], Table 1. Analytical data from the trial soil pit in 1985. (SOC = soil organic carbon; BD = bulk density).

Experimental Design
Random blocks were organized on an area of 3.6 hectares, with 18 trial plots of 720 m 2 (18 m × 40 m), resulting in the experimental design of three repetitions of 2 rotations × 3 tillage treatments. The study was conducted between 1985 and 2017. Two groups of trials were established (T1 and T2 trials). In each trial, three tillage treatments were tested: CT, MT, and NT with three repetitions (blocks). The T1 trial used winter wheat (Triticum aestivum L.) monoculture, while the T2 trial used a wheat-vetch (Triticum aestivum L.-Vicia sativa L.) rotation. There was no fallow period at any time in any of the plots. Annual fertilization for T1 trial was 80 units (kg) of nitrogen per hectare in the T1 trial and 40 units per hectare in the trial T2. Before drilling, weeds were controlled with a non-selective herbicide (glyphosate) in both trials by adding 1.6 to 2.0 L per hectare depending on the growing of weeds. In the T2 trial, other specific granule post-emergence herbicides based on tribenuron methyl or on bromoxynil were used to control broadleaved weeds.
In the CT treatment, the soil was managed with one pass of moldboard plowing at an average depth of 30 cm, followed by one or two passes with a spring tine cultivator (15-20 cm depth) for sowing preparation. Whenever needed, disk harrowing was applied before moldboard plowing for weed control after harvesting the previous crop. Sowing was performed in early November, after tillage with the chisel. In the MT treatment, the soil was managed with a chisel plow to an average depth of 15-20 cm. Weeds were controlled before tillage, and seedbed preparation and sowing were carried out according to that described for the CT treatment. The NT treatment consisted of sowing with a no-till drill. In all treatments, crops were harvested in early July.
This experiment started in 1985 with the aim to evaluate changes in SOC; originally CT was compared with MT and NT, considering three different rotations systems, wheat monocropping, vetch-wheat, and pea (Pisum sativum)-wheat, the latter was changed to vetch in 2002 [55]. This research follows the same experimental design albeit includes only plots having monocropping and wheat-vetch rotations over time.

Soil Sampling and Laboratory Methods
In 2017, soil samples (composite samples of three subsamples) were randomly taken in each of the plots at two depths: between 0 and 10 cm and between 10 and 20 cm. Soil stratification (SR), calculated as the ratio between the soil parameter values in the topsoil (0-10 cm) and those at 10 to 20 cm depth, was determined for SOC as an index of soil quality [56]. For bulk density (BD) determination, soil cores were taken (with a 100 cm 3 cylinder) and subsequently oven-dried until they reached a constant weight.
Samples were air-dried and sieved according to different procedures, to obtain specificsized fractions to carry out the soil tests described below.
The thirty-six sieved soils (<2 mm) were used for texture analysis using the pipette method [57]; SOC (%) was analyzed using the Walkley and Black method [58]. Following the FAO recommendations [59], we calculated the stock of C (Mg ha −1 ) based on BD and 30 cm thickness, following the Intergovernmental Panel on Climate Change [60], Equation (1).
where Stock is the stock of C (Mg ha −1 ), conc. = C concentration (%), BD = bulk density (Mg m −3 ), d = soil thickness (m), and δ2mm = proportion of gravel (>2 mm), which was negligible in the soils of the study. Parts per thousand SOC variation over time were calculated using the 4/1000 initiative recommendations [61] (Equation (2)). The 4/1000" Initiative is part of the Global Climate Agenda for Action, promoting scientific cooperation and actions towards reducing greenhouse gas emissions through protecting and increasing SOC stocks yearly at a rate of 4/1000.
where SOC f = final SOC concentration (%), SOC i = initial SOC concentration (%), and yr = years considered for the calculation. Total nitrogen was determined using the Kjeldahl method [62], available phosphorus content using the Olsen method [63], and exchangeable potassium [64] was measured with atomic emission spectrometry. The water-stable microaggregates (WSA) test was calculated as the percentage of aggregates <2 mm diameter resistant to wet sieving following the Kemper and Rosenau method [65]. Microaggregate samples were submerged and emerged over a 0.25 mm sieve at a rate of 30 oscillations per minute for 3 min. These calculations were corrected for sand content. The WSA test was conducted on three subsamples of 5 g from each of the soil samples.
Sieved soils between 4 and 4.75 mm were analyzed for macroaggregate stability using the CND method (number of drop impacts needed to break aggregates) [66,67]. Thirty macro-aggregates per soil sample were used.
To determine the macro, meso, and microporosity; total porosity; available water capacity (AWC); and BD, 36 undisturbed soil cores were taken (1 sample per treatment × 3 tillage treatments × 3 blocks × 2 rotations × 2 depths). Cores of 5 cm in diameter and 100 cm 3 were used. These cores were water saturated via capillarity in a sandbox, and a progressive series of increasing suctions was applied afterward. The sandbox was used for suctions ranging from 0 to 2 pF (0-10 kPa). The Richards' pressure plate system was used for pressures ranging from 2.54 to 4.2 pF (34 to 1500 kPa). Finally, samples were completely oven-dried (24 h at 105 • C) to calculate dry weight and BD. Different pore size groups were established following the classification suggested by Taboada et al. [68,69] and Bienes et al. [70]. These authors separated macropores (pores with diameter > 60 µm; corresponding to suctions between 0.1 and 6.3 kPa), mesopores (pores between 60 and 10 µm; corresponding to suctions between 6.3 and 34 kPa), and micropores (pores ≤ 10 µm; ≥34 kPa). AWC was calculated as the difference between the moisture held in soil microp- Infiltration rates were obtained using a simple-ring infiltrometer of 12.7 cm (5 in) in diameter and recording the time necessary to infiltrate 25 mm of water [71] by performing 10 repetitions at the same position. Overall, 18 infiltration tests were carried out (one test per treatment × 3 tillage treatments × 3 blocks × 2 rotations).

Statistical Analysis
The number of repetitions was different for different variables due to the diverse number of tests performed for each measurement. Twelve repetitions were performed for physical-chemical variables and laboratory analysis was in duplicates; six repetitions were performed for variables related to soil porosity due to their little variability, and 3 infiltrations tests were performed per treatment due to operational reasons.
Kruskal-Wallis and Kolmogorov-Smirnov tests, which are based on ranks rather than means, were used as non-parametric alternatives to a one-way variance analysis [72]. The analyses were developed in different stages using Statistica 10 software [73] and PAST 3.2 [74]. Figure 1 shows the evolution of carbon sequestration at depths of 0-10 and 10-20 cm, under different tillage systems since the beginning of the experiment over the elapsed period of 30 years; the figure does not include the rotation effects. Results are shown at 10-year intervals, as this period is considered appropriate for observing soil changes [75]. Data for the first two decades (historic data) were obtained from Hernanz et al. [55]. The SOC in 2016 was higher in the topsoil (0-10 cm), with the highest values obtained under the NT treatment (22.2 Mg ha −1 ). SOC was significantly lower (p < 0.05) under the MT (14.4 Mg ha −1 ) and the CT (10.3 Mg ha −1 ) treatments. In the first 10 cm, there was a gradual and significant increase in SOC in both the MT and NT treatments, which was faster in NT as can be observed for the last period of study in Figure 1. However, there were no major changes in the 10-20 cm depth, and there were no significant differences between tillage treatments. SOC managed using CT showed no change over time.  7.5 ± 0.7 c 10.6 ± 1.9 b 15.4 ± 2.2 a 10.8 ± 3.5 a 11.5 ± 3.9 a SOC 10-20 cm (g kg −1 ) 7.5 ± 0.9 ab 7.8 ± 0.7 a 6.6 ± 0.9 b 7.2 ± 0.9 a 7.4 ± 1.0 a Historical data were taken from Hernanz et al. [55], while the 2016-2017 data are new. Under the NT treatment, the SOC stock grew rapidly in the first years. However, in subsequent years, although SOC continued to increase, the rate of growth slowed. Under the MT treatment, the response was slower, and over the last 20 years, the increase in C Agronomy 2021, 11, 291 6 of 15 stock was almost stable. Considering the 0 to 20 soil thickness, there was a significant increase in SOC stock as the sum of the two layers (0 to 10 and 10 to 20 cm) ranged from 20.8 Mg ha −1 in CT, compared to 26.5 in MT and 32.2 in NT, however, increases in SOC were only attributed to the topsoil (0-10 cm depth) as variations in the 10-20 cm depth were negligible. Consequently, over the 32 years, soil management (tillage and crop rotation factors) led to a stratification (SR) in the topsoil. The SR was more evident under the NT treatment, and the SOC content was particularly affected ( Table 2). The SR revealed the effect of tillage intensity on SOC accumulation (Table 2). This soil parameter was higher in the 0-10 cm soil depth compared to the 10-20 cm depth, and the accumulation in the topsoil gradually increased as the tillage intensity decreased (CT < MT < NT). Furthermore, regardless of the type of tillage, under the wheat-vetch rotation, there was an accumulation of SOC in the upper 10 cm of soil. Table 2. Stratification relationships (SR: contents at 0-10 cm/contents at 10-20 cm) of SOC, total nitrogen (N), and available phosphorus (P) and potassium (K) under different tillage types and crop rotations. Data obtained from the last two years of the experiment.

Tillage (n = 12)
Crop rotation (n = 18) .5 a 11.5 ± 3.9 a SOC 10-20 cm (g kg −1 ) 7.5 ± 0.9 ab 7.8 ± 0.7 a 6.6 ± 0.9 b 7.2 ± 0.9 a 7.4 ± 1.0 a SR SOC Similarly, an increase of N was observed in NT treatments at 0-10 cm depth, in this case, the SR indicated that N content at the topsoil nearly doubled the value found at 10 to 20 cm depth (SR = 1.96). The effect of rotation also showed an increasing trend although variability impeded statistical significance.
The content of P and K followed the same pattern: higher concentration at the surface (0-10 cm) and higher concentration as tillage intensity decreased. This pattern yielded progressively higher SR for P and K due to tillage treatments, but it was not detected for rotation treatments.

Changes in Physical and Chemical Properties
In order to understand the combined effects of rotation and tillage, changes in soil parameters are collated in Table 3, in which soil depths between 0 and 20 cm are considered together. Soils managed using monocropping methods are more sensitive to tillage treatments. The CT treatment significantly decreased aggregate stability, SOC, TN, and P, and especially noticeable was the decreased AWC. On the contrary, CT increased macroporosity and mesoporosity and, therefore, decreased BD compared to the other two tillage systems. On average, there were no significant differences in total porosity between treatments. Table 3. Results of Kolmogorov-Smirnov test, showing the influence of tillage and rotation on plots with wheat-wheat monocropping and a wheat-vetch rotation in 2015 and 2016, after 32 years of treatments in 0 to 20 cm of soil thickness. Significant differences are marked at p < 0.05. Average values and standard deviation are stated. Different letters in each row indicate differences between different tillage types for the same cultivation system (monocrop or rotation); any differences found between rotation systems under the same tillage system are marked by an asterisk (*).  The wheat-vetch rotation seemed to mitigate the effects of tillage treatments documented above. There were then no differences in aggregate stability or BD between tillage methods when rotation wheat-vetch was performed.

Factors (Rotation and Tillage
When comparing the tillage treatment under different rotation systems (different capital letters in Table 3), we found that rotation produced higher mesoporosity, which ranged from 12.5 to 15.5 under monocropping and 18.8 to 20.5 under rotation. No differences were found for nutrients due to rotation.
The intensity of tillage progressively increases from NT to MT and CT; this rank is reflected in the results of the variables studied, and the MT methods usually showed intermediate values between NT and CT.
From the methodological point of view, the changes in microaggregate stability were not detected using the WSA method.

Changes in Water Infiltration and Water Availability
The average infiltration measured using a single ring infiltrometer was 96 mm h −1 . No significant differences were noticed between the CT and MT treatments, or between crop rotations (wheat-wheat, wheat-vetch). However, under the NT treatment, there was less water infiltration in soils.
With regard to AWC, the CT treatment reduced significantly this parameter compared to NT in monocropping (14.6% compared to 16.6% respectively). Rotation with vetch reduced even more this AWC although only for NT treatment (16.6% in W-W compared to 13% in W-V). As mentioned before, total porosity was similar between treatments, but when considering the ratio AWC/TP, there was a significant decrease in soils managed using CT and monocropping; this ratio was less than 30% for CT treatment, compared to 33.3 and 37.8% in MT and NT respectively. Therefore, this management practice reduced the ability of soil to hold water as the AWC/TP ratio provides information about the percentage of TP that is occupied by available water for crops. There were then no differences in AWC between CT and MT methods when rotation wheat-vetch was performed. The water retention curves of soils under different treatments (tillage and rotation) can be observed in Figure 2. On the left, soils under monocropping experienced significant differences between porosity when soils were not wet; these differences were only noticed between NT and CT, as is also indicated in Table 3. Particularly, the little variations in porosity for the NT treatments were noted. On the right, soils under rotation showed similar values of porosity, and therefore, similar AWC, with no significant differences between tillage treatments. vetch reduced even more this AWC although only for NT treatment (16.6% in W-W compared to 13% in W-V). As mentioned before, total porosity was similar between treatments, but when considering the ratio AWC/TP, there was a significant decrease in soils managed using CT and monocropping; this ratio was less than 30% for CT treatment, compared to 33.3 and 37.8% in MT and NT respectively. Therefore, this management practice reduced the ability of soil to hold water as the AWC/TP ratio provides information about the percentage of TP that is occupied by available water for crops. There were then no differences in AWC between CT and MT methods when rotation wheat-vetch was performed. The water retention curves of soils under different treatments (tillage and rotation) can be observed in Figure 2. On the left, soils under monocropping experienced significant differences between porosity when soils were not wet; these differences were only noticed between NT and CT, as is also indicated in Table 3. Particularly, the little variations in porosity for the NT treatments were noted. On the right, soils under rotation showed similar values of porosity, and therefore, similar AWC, with no significant differences between tillage treatments.

Discussion
There were no significant differences in crop yield between treatments (average 2.5 ± 1.3 Mg ha −1 , Hernanz et al. [76]). However, several changes were observed in soil parameters after 32 years of different management practices.

Discussion
There were no significant differences in crop yield between treatments (average 2.5 ± 1.3 Mg ha −1 , Hernanz et al. [76]). However, several changes were observed in soil parameters after 32 years of different management practices.

Changes in SOC
In 1985, the soils had 10.6 Mg ha −1 of SOC in the upper 10 cm layer and 9.4 Mg ha −1 in the 10 to 20 cm layer (Figure 1). Considering the upper 20 cm layer of soil, the SOC stock in 2017 had reached an average of 20.8 Mg C ha −1 under CT management, 26.5 Mg C ha −1 under MT management, and 32.2 Mg C ha −1 under NT management. This implies an annual average increase of 10 per thousand under MT management and 14 per thousand under NT management over the 32 years considered. The progression of carbon accumulation suggests that these soils have not yet reached their maximum accumulation potential or their stabilization curve. However, the accumulation rate decreased after 11-12 years [55].
In the upper 20 cm, SOC under the CT treatment was approximately 7 g kg −1 , which is below the limit of 10 g kg −1 considered typical of poor soils that are vulnerable to desertification [77]. Under the NT system that preserves soils, the current carbon stock is 32.2 Mg C ha −1 , which is higher than the value reported by Rodríguez Martín et al. [78] for central Spain, but much lower than the average value for Spain (56 Mg C ha −1 , ranging from 20 Mg C ha −1 in the center and south of the country, and 150 Mg ha −1 in the cooler and wetter northern zones). The low initial rates of SOC in the topsoil (0-10 cm depth), may explain the high carbon sequestration figures found in this updated study, which shows that the initial value has doubled over the last 32 years of NT management. As expected, the greatest increases in SOC were observed in the topsoil, under the NT and MT treatments, although in the latter case the changes were not as pronounced. This may be a consequence of the accumulation of plant litter in the topsoil. Continuous carbon sequestration for more than 20 years may be common in certain conditions such as those described in this study and others in the literature. For example, Minasny et al. [79] reported cases where there was effective carbon sequestration over 30 or 40 years after different stages of accumulation and equilibrium [80,81]. Other research, based on modeling the impacts of land use change on SOC, had demonstrated potential carbon sequestration for even longer periods, for example, 50 or 100 years [82]. Furthermore, these authors highlighted the important role that the conversion of conventional agriculture to non-tillage could achieve. It is recognized that reducing emissions coming from agriculture is imperative [81] and sustainable agricultural practices have a significant influence on SOC sequestration [59], hence, changes to NT systems could make up for the total carbon emissions of agriculture globally, especially if conserving SOC is a "least-cost opportunity" for climate change adaptation and mitigation.
Natural ecosystems and managed grasslands and forests show a marked horizontal distribution of SOC and other variables. This is also frequently observed in previously degraded croplands which have been restored using conservation agriculture, in particular, NT management practices [75]. For this reason, the degree of stratification is considered an index of soil improvement and soil quality [55]; this is due to the positive effect of surface organic matter for erosion control, water infiltration, and the conservation of nutrients. With this in mind, the NT treatment produced the highest carbon stratification rate (SRSOC), which was 2.32 times higher than under CT, and in turn, SR for MT treatment was 1.36 times higher than in CT (Table 2). Consequently, in this specific context, NT can be regarded as a management practice that improves soil quality and provides a lower susceptibility to erosion.
According to Peregrina et al. [83], enrichment rates greater than two for SOC indicate improvements in soil quality. This limit was only exceeded under the NT treatment. Barbera et al. [47] reported that the annual sequestration of SOC with wheat monocropping was 2.75 times higher than with a wheat-faba bean rotation. However, in this study, considering the data available from the last years (2016-2017), and after 32 years, we did not find significant differences in SOC between wheat monocropping and wheat-vetch rotation.

Changes in Physical and Chemical Properties
The accumulation of vegetation residues on the topsoil is greater for the NT treatment; such organic material is a source of SOC, N, P, and K for the soil. This accumulation leads to the stratification found between topsoil (0-10 cm) and the layer underneath (10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20), which was remarkable for different tillage treatments. As expected, CT impeded stratification, and SR was not observed in this treatment, as SR values were 1.01, 1.16, 0.92, and 1.21 for SOC, N, P, and K respectively. On the contrary, NT treatment favored stratification, and SR values were 2.32, 1.96, 1.80 and 1.78 for these variables. The effects of MT were in between.
The pattern of concentrations found for SOC and total N were similar due to the high correlation between these variables, as usually organic forms of N represent more than 90% of total N in soils [84]. In addition, plant activity includes translocation of nutrients from deeper soil layers to the topsoil, these nutrients are not totally used by the crop near the soil surface [85], in this case from 0 to 10 cm.
Changes were also evident for P; the MT and NT treatments showed higher concentration of P and higher stratification compared to CT. It is well known that the availability of this nutrient is limited in calcareous soil suffering high fixation and little penetration [86] and hence high stratification.
Concerning K, it was significantly higher in the topsoil of NT treatment; this available form is bound to the surface of clays and humic substances and may be exchanged with other cations. In this semi-arid environment, the movement of K is limited by low moisture [87], and consequently, lixiviation is not favored, which leads again to stratification in less intensive tillage practices.
Soil structure plays a critical role in soil functions and can be measured using aggregate stability tests [88,89], which in turn can be used as soil quality indicators [90,91]. The CND test is sensitive to soil management practices [92]. In this study, the CND test values gradually decreased with increasing intensity of tillage systems following the order of macroaggregate resistance, NT > MT > CT, which coincided with gradual decreases in SOC content.
The CND test values changed with different soil tillage practices [27], but tests considering other structural properties, such as WSA (which defines the microaggregate stability), were not sensitive to these tillage practices. The aggregate stability also influenced, either directly or indirectly, other physical and chemical properties of the soils, so it can also be used as an indicator of soil degradation [90]. The macroaggregate stability was related to SOC and its protection from the activity of microorganisms [88]. Since roots and hyphae stabilize macroaggregates, macroaggregation is also affected by soil management. On the contrary, the water stability of microaggregates is strongly conditioned by persistent organic binding agents and seems to be independent of management [93]. Microaggregate stability was also related to clay content and probably to cationic links [94] that may be important in these argillo-calcareous soils.
Both TP and BD were highly influenced by tillage treatment and SOC content [23,95]. Usually, CT reduces BD in the topsoil [4,14]. In this study, BD increased under NT treatment. Similar observations have been made by other studies, for example, Costa et al. (2015) [96] . Increased BD was related to a decreased percentage of macro and mesoporosity. The lack of rotation exacerbates this influence, and BD was lower under the CT wheat monoculture treatment (1.36 Mg m −3 ) and higher under the NT monoculture (1.51 Mg m −3 ) treatment. Soil compaction has been frequently mentioned as one negative effect in soils under semiarid climates [97]. These differences were not as strong in plots under rotation, where there were no significant differences in BD between tillage treatments. Other variables followed the same pattern, which leads us to conclude that the effect of rotation (vetch-wheat in this study) may cushion, and even cancel out, the effects of different tillage systems in these particular soils.

Changes in Water Infiltration and Water Availability
Increased macroporosity as a result of CT is expected to influence the movement of water through soil [98,99]. However, several authors have found that NT systems yield higher water infiltration than CT systems [100,101]; other meta-analysis found that no-till did not consistently improve infiltration rates [102]. In this study, under the NT treatment, there was less infiltration than under the CT or MT treatments, despite the few differences in porosity. Infiltration tests were conducted in April 2017, a particularly dry year, during the period from January to April, when only 77 mm of rain was recorded in the area (60% lower than average according to the National Meteorological Agency records). Owing to the drought, the soil had multiple and large cracks which may overestimate infiltration rates. Under the NT system, soils had fewer cracks and therefore less infiltration than observed under the CT and MT systems. Vetch root systems may be less efficient than wheat root systems in creating effective connectivity between pores. Therefore, wheat-vetch rotations can hamper infiltration, when compared to monoculture in the local conditions studied here. The dry meteorological conditions during the experimental tests could have influenced these results; further research is needed to ascertain this hypothesis.
Under monocropping, AWC was higher under the NT treatment and lower under the CT treatment (p < 0.05; Table 3) owing mainly to the increased water content at field capacity. This suggests that in rainy years, soils managed under NT can hold more water than soils under CT.
Soils under NT showed PWP values higher than under other tillage systems. This can be a disadvantage in dry years, therefore, it may be a problem for adaptation under future scenarios of climate change. This was not observed in plots with rotations, therefore, we can infer then that crops managed by rotation may be able to adapt more efficiently under future climate change.
When considering the ratio of AWC and TP (AWC/TP), we found that values were higher (although not significantly different) under the MT and NT treatments, compared to those under the CT treatment. This may indicate that under the CT treatment soils have a lower percentage of porosity able to hold water, which again, may be related to SOC content under moderate or low-intensity tillage systems. Water availability can be improved by NT systems in spite of the increase of bulk density, which may be due to the development of a new pore-size distribution under NT compared to CT [103]. We observed that CT reduced AWC if soils were managed using monocropping, but this negative effect was not observed in treatments with rotations. In this later case, rotation and NT dry soils increased the water retained at the permanent wilting point and this reduced AWC.
The NT may not be always beneficial under all soil types and climate scenarios [104,105] and should be studied under local conditions, particularly in semi-arid conditions [106].

Conclusions
The agricultural soils managed using no-tillage (NT) in central Spain increased their soil organic carbon (SOC) content for 32 years. These changes in SOC are significant for the upper 10 cm of topsoil but are less evident from the 10 to 20 cm depth; such a difference defined as SOC stratification is a positive indicator related to soil conservation. Due to this change, the NT treatment improved soil properties related to structure and porosity and showed higher available water capacity (AWC).
However, it has been argued that NT may reduce yield in dry years (in which soil moisture does not reach the field capacity) because these NT soils held more water at wilting point compared to those managed under minimum tillage (MT) or conventional tillage (CT). To improve water availability, it is essential that NT systems are combined with crop rotations, as under this combined management, no significant differences were found in AWC. Another negative effect of NT was the decrease in infiltration, which can be important in this semi-arid climatic condition.
The effects of tillage intensity were more pronounced under monocropping treatments, e.g., bulk density, SOC, field capacity, permanent wilting point, and AWC. This suggests that wheat-vetch rotations may mitigate the negative effects of tillage on soil structure, stability, and water holding capacity. Under the MT system, most of the soil parameters had intermediate values compared to those recorded under the other two tillage systems.
Further long-term studies are needed in other different areas to confirm the effects of management practices on soil quality and water availability, to persuade growers to adopt more sustainable techniques as an alternative to CT in semi-arid areas and to determine the advantages of combining NT with rotations.