Long-Term Effects of the Use of Organic Amendments and Crop Rotation on Soil Properties in Southeast Spain

: The evolution of soil chemical properties over 20 years was monitored to assess the effects of the change in soil management from a rainfed to an irrigated model and the use of organic amendments and crop rotation. Intensive agriculture has been the activity that has caused most degradation and contamination of this soil. Long-term monitoring of the soil proﬁle made it possible to assess its response to the application of sustainable agricultural techniques intended to offset these effects. Three proﬁles of the same soil were studied—P1 (1998), P2 (2003), P3 (2017)—to show the evolution in time and space. An incipient degradation process was detected in the ﬁrst ﬁve years, veriﬁed by increases in salinity (2.3 dS m − 1 ), exchangeable Na (0.5 g kg − 1 ), and TN (1.3 g kg − 1 ) in P2 in comparison with P1 (1.0, 0.2, and 1.1, respectively). There was also leaching towards the deep horizons for TN (0.4, 0.9, and 0.7 g kg − 1 for P1, P2, and P3, respectively), and for assimilable elements such as P (1.1, 6.4, and 3.8), Fe (2.0, 2.1, and 5.6), Mn (0.3, 6.5, and 1.9), Zn (0.3, 0.5, and 0.9), and Cu (0.5, 0.6, and 1.3) (all mg kg − 1 , for P1, P2, and P3, respectively). Between 2004 and 2017, organic amendments (sheep manure) were reduced by 50%, crop rotation was intensiﬁed, and green fertilization and forage maize cultivation were included. As a result, P3 showed an improvement in comparison with P2, with decreases in EC (1.4 dS m − 1 ), exchangeable Na (0.2 g kg − 1 ), and TN (0.8 g kg − 1 ). The change in soil management enhanced some soil functions (carbon sink and chemical fertility) and attenuated degradation.


Introduction
Since the Second World War, in developed countries and progressively in developing countries, the agricultural sector has undergone major changes aimed at increasing the production and profitability of farmland and, therefore, the ability to compete in an increasingly globalized market. Intensive agriculture has probably been the activity that has caused the most degradation and pollution in the soil and associated ecosystems [1]. In arid and semiarid areas, 30% of irrigated agricultural soils have been rapidly degraded due to the trend to excessively apply inorganic fertilizers and manures [2], to utilize poor quality water for irrigation [3][4][5][6][7], and to perform inadequate tillage [8,9].
As a result of such problems, society in general and, on its behalf, its leaders, are promoting other models of land management that are more respectful of the environment. In this scenario, the use of organic substrates, as substitutes for chemical fertilizers, and crop rotation are promoted by the proponents of sustainable agriculture and form part of the recommendations of the European Union in its Common Agricultural Policy (CAP). A significant proportion of the cultivated soils have been subjected to the shift from rainfed, extensive agriculture to an irrigated model. Concurrently, traditional rainfed crops (cereals and almonds) have been replaced with horticultural crops (broccoli, tomato, celery, potato, etc.), whose management implies changing from an agrosystem that is only slightly aggressive with respect to the environment to another in which processes such as soil degradation, pollution, and overexploitation of aquifers are frequent. To investigate the effects of management on the soil properties, the horizons of three soil profiles were sampled. Profile 1 (P1) was sampled in 1998 and corresponded to a soil used for rainfed cereal cultivation. Profiles 2 and 3 (P2 and P3) were sampled in the same location, but in 2003 and 2017, respectively. The use of irrigation water on the farm began at the end of 1998, after the sampling of P1.
The hypotheses that justified this study were: -The change in the use to which the studied soil has been subjected, as a result of the introduction of irrigation and horticultural crops, entails a potential risk of accelerated soil degradation, as has happened in other, similar soil and climate zones. This must be manifested in the long term in some morphological properties, such as the color and thickness of the surface horizons, chemical properties (organic carbon (OC), electrical conductivity (EC), total nitrogen (TN), assimilable elements, etc.), and mineralogical properties of the soil. However, if the management involves sustainable agricultural techniques, such as the use of organic amendments and crop rotation, and if the evolution of soil horizons is monitored in the medium and long term, the degradation of this resource can be minimized.
Although agricultural activity takes place mainly in the upper horizons of the soil, changes in soil use can also affect deeper horizons, which are rarely studied. Therefore, it is necessary to monitor the soil profile as a whole.
The objectives of this study were as follows: The hypotheses that justified this study were: -The change in the use to which the studied soil has been subjected, as a result of the introduction of irrigation and horticultural crops, entails a potential risk of accelerated soil degradation, as has happened in other, similar soil and climate zones. This must be manifested in the long term in some morphological properties, such as the color and thickness of the surface horizons, chemical properties (organic carbon (OC), electrical conductivity (EC), total nitrogen (TN), assimilable elements, etc.), and mineralogical properties of the soil. However, if the management involves sustainable agricultural techniques, such as the use of organic amendments and crop rotation, and if the evolution of soil horizons is monitored in the medium and long term, the degradation of this resource can be minimized.
Although agricultural activity takes place mainly in the upper horizons of the soil, changes in soil use can also affect deeper horizons, which are rarely studied. Therefore, it is necessary to monitor the soil profile as a whole.
The objectives of this study were as follows: -To establish the influence of organic amendments and crop rotation on some soil properties over a period of 20 years and to determine the intensity of degradation processes in the different horizons of the profile. - To promote a management model that enhances the soil ecological functions, such as the role of soil as a carbon sink, and minimizes soil degradation.
Regarding the innovative aspects of this work that were necessary to achieve the aforementioned objectives, the first concerns the sampling and monitoring unit (the pedon instead of arable layer samples). The soil was considered as a natural three-dimensional body formed by horizons, and its properties, including mineralogical and morphological ones, were monitored over two decades, in all horizons. This is vital because, although anthropic activity is carried out in the upper horizons, the consequences can be manifested in both the surface and deep horizons, and study of the latter is essential to reveal the processes of salinization and contamination that can affect the soil and aquifers. The second innovative aspect is the long-term nature of the study, which will provide reliable data on the soil response to the treatments and the management model applied (irrigation, organic amendment, subsoiling, crop rotation, green fertilizer, etc.), permitting the definition of a set of sensitive indicators for soil monitoring. This will allow managers to take measures to correct those variables that deviate from the desired value, and thus avoid soil degradation.

Materials and Methods
The experimental plot is located in the southeast of Spain (Figure 1), with an average annual precipitation of around 303 mm and an average annual temperature of 12.7 • C, and is classified within the supramediterranean bioclimatic stage. The geographic coordinates are P1: 37 • 21.20" W, the altitude is 948 m above sea level, and the area covers approximately 1 hectare with a natural slope of less than 2% (leveled to 0%) and has developed on alluvial sediments. The weather of the study area has a highly seasonal nature, with minimum temperatures of around −10 • C and maximum temperatures of around 40 • C. The soil moisture regime is xeric and the temperature mesic [23]. This plot is representative of the type of soil on which the region's irrigation plan was carried out in the 1990s, so the results obtained in this study can be extrapolated.
In the sampling design, the evolution of the soil properties was monitored, in addition to changes in them due to irrigation. In this way, it was possible to consider a change in the management model to avoid soil degradation, if necessary. The pedon was adopted as the sampling and monitoring unit, as defined by [24]. During the study period, three trial pits, approximately 200 cm deep, 80 cm wide, and 300 cm long, were opened by means of a backhoe and all the horizons were sampled on the three faces of each trial pit, obtaining three samples from each horizon, one frontal and two laterals, thus ensuring the representativeness of the samples.
Three complete soil profiles (pedons) were sampled, following the FAO Guidelines for Soil Description [25]. Profile 1 (P1) was sampled in 1998 and corresponded to a soil used for rainfed cereal cultivation. It represents the initial starting point before the change in land use. Profiles 2 and 3 (P2 and P3) were sampled in the same location, but in 2003 and 2017, respectively ( Figure 1). Both P2 and P3 corresponded to soils dedicated to summer vegetable cultivation under irrigated conditions and so had properties that, in part, were a consequence of the agronomic management that had been carried out on the farm for 5 and 19 years after the change of use, respectively. The distance among sampling points was close to 20 m, in order to avoid perturbations derived from previously made trial pits.
From 1998 to 2017, the plot was used, basically, to cultivate vegetables throughout the summer season: broccoli (3 years), tomato (5 years), pepper (2 years), lettuce (7 years), and corn (2 years). In addition, during the autumn and winter seasons, arable crops were cultivated with legumes (Vicia sativa L.) and oats (Avena sativa L.) as green manure (50 kg of oats and 125 kg of legume seeds were sown per hectare or even the soil was left bare (fallow) to meet CAP requirements). The water supply varied from 2500 to 6500 m 3 ha −1 , depending on the summer crop in the plot and on the cycles of repetition of the crop, because the autumn-winter crops did not receive any irrigation ( Figure 2). The irrigation water used was hard, with a pH of 7.6 and an EC of 0.8 dS m −1 , so it can be considered of good quality according to the combined Wilcox and Green index, based on the ratio of % Na to dissolved salts in the soil water. The irrigation water was of subterranean origin and its composition is presented in Table A1. The irrigation and fertilization were carried out simultaneously (fertigation) using drip irrigation. The cultivation system had a 2000 l tank for the nutrient solution with programmed stirring; electrovalves, filters, and polyethylene tubes (25 mm in diameter) were used to distribute the water or nutrient solution, and 14 mm diameter tubing was used to supply the self-compensating emitters (2 l h −1 ). Irrigation control was achieved with a set of Watermark electrical resistance Soil samples were air dried, sieved with a 2 mm sieve, and processed in or analyze their properties and features [25]. The description of the consistency as dry ness), moist (friability), or wet (adherence and plasticity) was carried out according indications of [25]. Soil color was described using a Munsell Soil Color Chart, soil t was determined by the Robinson pipette method [26], and TOC and TN were deter in an elemental analyzer (Leco, model CHNS-932). The estimation of total carbon c Soil samples were air dried, sieved with a 2 mm sieve, and processed in order to analyze their properties and features [25]. The description of the consistency as dry (hardness), moist (friability), or wet (adherence and plasticity) was carried out according to the indications of [25]. Soil color was described using a Munsell Soil Color Chart, soil texture was determined by the Robinson pipette method [26], and TOC and TN were determined in an elemental analyzer (Leco, model CHNS-932). The estimation of total carbon capture was estimated from the OC and CaCO 3 contents, respectively, according to the equations: where OCc is organic carbon capture in Mg ha −1 ; OC (g kg −1 ); BD is bulk density (1.25 g cm −3 ); depth is the depth of each horizon (m); and where ICc is the inorganic carbon capture in Mg ha −1 ; CaCO 3 (g kg −1 ); AWC is atomic weight of carbon; BD is bulk density (1.25 g cm −3 ); depth is the depth of each horizon (m); MMCaCO 3 is the molecular mass of CaCO 3 . Moreover, pHw was determined in a 1:1 suspension of soil in water and pHKCl in a 1:1 suspension of soil in 1 M KCl [27], and EC in a saturated paste extract [28]; total carbonates (CaCO 3 ) were estimated by volumetric analysis using a Bernard calcimeter [28]; the cation exchange capacity (CEC) and the exchangeable Na, K, and Mg were determined by means of the method described in [28]; the available contents of the cations Na, K, and Mg were measured by atomic absorption; and the available P was determined by Watanabe and Olsen's method [28]. Atomic absorption was used to determine the available contents of the elements Fe, Cu, Mn, and Zn-that is, the contents extracted with a solution of 0.05 M DTPA, 0.01 M CaCl 2 , and 0.1 M triethanolamine at pH 7.3 [28].
Regarding the mineralogical analysis, the clay fraction was subjected to the following treatments prior to its analysis in oriented aggregates: removal of organic matter and carbonates, saturation with Mg 2+ and ethylene glycol, and heating to 350 and 500 • C [29]. All samples were analyzed using an X-ray diffractometer and the minerals were identified using the database [30]. The working conditions were radiation Kα Cu; intensity 24 mA; 40 kV; Ni filter; window slit of 1; meter gap of 0.1; scanning speed, 1 • C 2Θ min −1 ; and sensitivity, 5 × 10 3 .
The data were analyzed using the General Linear Model of the SPSS Version 22.0 statistical package (IBM SPSS,). A one-way ANOVA and Tukey's test were carried out in order to ascertain significant differences among the P1, P2, and P3 horizons.

Results
Profiles P1, P2, and P3 were relatively well developed, containing an Ap-Bw-C sequence of horizons. According to the results of the macromorphological study of the complete profiles (Tables A2-A4), land use change caused modifications of differing importance in some soil properties. There were changes in the depth of the Ap horizons. In fact, in P2 and P3 the Ap horizon was divided into Ap1 and Ap2, and was more than 50 cm deep, whereas the Ap horizon of P1 was around 30 cm deep. The depths at which the rest of the horizons appear in profiles P2 and P3 have also been affected by surface leveling.
The soil texture in P1, P2, and P3 was silty clay loam for Ap and Bw, clay loam for C1, and loam for C2, except in P1, for which it was clay loam. The soil structure was blocky subangular for Ap and Bw in the three profiles.
Regarding soil color, the Ap1 horizon darkened very quickly from P1 to P3, changing from pale brown (10YR 6/3-dry) and dark yellowish-brown (10YR 4/3.5-moist) to brown (10YR 5/3-dry) and dark brown (10YR 3/3-moist). A whitish mottling was observed in the Ap1, Ap2, and C1 horizons of profiles P2 and P3, and analysis showed that this comprised soluble salts ( Figure A1). Furthermore, the Bw horizon had reddish tones, and a lower CaCO 3 content (Table 1) than the underlying horizons, so it was denoted as a cambic horizon. This feature was relevant for the soil taxonomic classification. Table 1. Total organic carbon (TOC), total nitrogen (TN), C/N ratio, electrical conductivity (EC), CaCO 3 content, and pH (water and KCl). Profiles P1, P2, and P3 had low to moderate TOC contents that decreased progressively with depth ( Table 1). The values were highest in the upper horizons (especially in P2 and P3, ≈10 g kg −1 ) and lowest in the C horizons (<3 g kg −1 ). The statistical analysis revealed significant differences in the TOC contents among P1, P2, and P3, specifically in the Ap1, Ap2, and C1 horizons.
The TN contents (Table 1), in general, did not vary significantly with the change in use of the cropland. As happened with TOC, the highest values of TN were found in the upper horizons of P1 and P2 (1.1 and 1.3 g kg −1 , respectively). However, N accumulation was detected in the C1 and C2 horizons of P2 and P3, with a maximum of 0.9 g kg −1 . This value even exceeded that measured in the surface horizon of P3 (0.8 g kg −1 ).
The C/N ratio differed significantly among Ap1, Ap2, and C2; it had values between 7 and 9 in horizons Ap and Bw of P1 and P2, whereas the maximum value (16.5) was found in horizon Ap2 of P3 (Table 1). The lowest C/N ratio found was in the C horizons of P2 and P3, a result of the increase in N mentioned in the previous section. The C/N ratios in the upper horizons (Ap1 and Ap2) of P3 were increased in comparison with P2 and were the closest to the values of the surrounding natural soils (C/N≈10), due to the presence of high concentrations of TOC.
The EC values varied significantly in all horizons with the change in the use of the cropland (Table 1). In P1, the EC was lower than 2 dS m −1 in all horizons of the profile. However, in P2, the EC was significantly (p < 0.05) increased in comparison with P1 in all horizons except Bw (1.2 dS m −1 ); therefore, the first 5 years of vegetable cultivation led to statistically significant increases in the EC. This increase in the EC was accentuated in the deep horizons (C1 and C2) of P2 that, in comparison to those of P1, reached EC values around four times higher. This aspect confirms the leaching of salts to the deep horizons during this first stage and an incipient salinization of the soil.
The dynamic of the total CaCO 3 in the soil profile suggests its mobilization from the upper horizons, where its concentration was lower in P3 than in P1 (by 6% in Ap1 and by 7% in Bw), to the underlying C2 horizon, where its concentration rose by 20.8%. However, as can be seen, the decrease in surface CaCO 3 was smaller than the increase at depth, indicating that inorganic carbon sequestration was occurring in the deep horizons ( Table 2). The incorporation of carbonates and bicarbonates with the irrigation water, and the bioformation of CaCO 3 due to the root respiration of the crops, may explain these results. By contrast, organic carbon sequestration occurred mainly in the upper soil horizons, turning this soil into a carbon sink. This is supported by the fact that during the study period an increase in OCc of 35.0 Mg OC ha −1 was observed, most of which occurred in the upper horizon (31.5 Mg OC ha −1 ).
The pHw values of P1 were basic, ranging between 8.3, for the Ap1 horizon, and 8.6, for the C1 horizon ( Table 1). The pHw values decreased slightly with the change in land use. Such decreases in P2 and P3, compared to the P1 values, were found in all horizons, both upper and deep.
The CEC values were relatively high in the A horizons, due to the presence of 2:1 clays ( Figure 3) and humus, and decreased with increasing depth in the profile. Regarding the concentrations of basic cations, the change in soil use caused a short-term enrichment in Na, K, and Mg in P2 and a decrease in P3 in comparison with P2 (Table 3), with this trend being statistically significant in the Ap, Bw, and Ck1 horizons throughout the study period.
Agronomy 2021, 11, x FOR PEER REVIEW 9 By contrast, organic carbon sequestration occurred mainly in the upper soil horiz turning this soil into a carbon sink. This is supported by the fact that during the st period an increase in OCc of 35.0 Mg OC ha −1 was observed, most of which occurre the upper horizon (31.5 Mg OC ha −1 ).
The pHw values of P1 were basic, ranging between 8.3, for the Ap1 horizon, and for the C1 horizon ( Table 1). The pHw values decreased slightly with the change in l use. Such decreases in P2 and P3, compared to the P1 values, were found in all horiz both upper and deep.
The CEC values were relatively high in the A horizons, due to the presence of clays ( Figure 3) and humus, and decreased with increasing depth in the profile. Regard the concentrations of basic cations, the change in soil use caused a short-term enrichm in Na, K, and Mg in P2 and a decrease in P3 in comparison with P2 (Table 3), with trend being statistically significant in the Ap, Bw, and Ck1 horizons throughout the st period.    The vertical distributions of the available P concentration and the available micronutrient concentrations (Fe, Mn, Cu, and Zn) were similar in P1, P2, and P3; they were highest in the upper horizons, with significant increases throughout the study period (Tables 3 and 4). In the case of P, this implies an increase in chemical fertility, due to its functions as a macronutrient. In the case of Fe, the initial values tripled with the change of use, both in the arable layer and in the deep horizons, whereas the concentration of Mn in the arable layer doubled over the years of the study, with extremely high values in the deep horizons of P2, which subsequently stabilized in P3. The Cu values were doubled in all soil horizons and the Zn values were increased up to six-fold, when comparing P3 with P1. Finally, a statistically significant correlation was found between OC and available P and available micronutrients, so that the maximum concentrations of available micronutrients and P occurred in the horizons with the highest organic matter content.
The clay fraction of P1 and P3 was characterized using the method described in Section 2. As can be seen in Figure 3, there were no appreciable differences between the corresponding horizons of P1 and P3. The Ap horizons of both profiles featured the presence of illite (reflection at 1.0 nm) and kaolinite (0.7 nm). In the Bw and C1 horizons of P1 and P3, together with illite and kaolinite, there were trioctahedral smectites, characterized by their reflection at 1.4 nm in the Mg diagram, which expanded at 1.7 nm upon solvation in EG and disappeared in the heat treatment. Chlorites were observed in these horizons, characterized by the presence of a reflection at 1.4 nm, which persisted in the heat treatment. Of these minerals with basal spacing at 1.4 nm, only traces of smectites were seen in the Ap horizon ( Figure 3). All of these can be considered as minerals derived from illite.
Profiles P1 and P2 were classified, according to [31], as Calcaric Vertic Cambisols (Loamic) because they were calcareous and possessed one cambic and one protovertic horizon. In contrast, P3 contained the above-mentioned horizon plus an A mollic horizon and a calcic horizon; consequently, it was classified as a Vertic Calcic Kastanozem (Cambic, Loamic).

Discussion
The change of use and management to which the soil was subjected has had consequences for its macromorphology, with the thickness of the A horizon and the depths of the different B and C horizons being affected. These were due to the subsoiling carried out at a depth of 50-60 cm every two or three years to eliminate the hardened tillage layer and improve the physical properties (in P2 and P3). This caused part of the Bw horizon in P1 to become part of the Ap horizon in P2 and P3, while the leveling of the soil affected both the thickness of the superficial horizon and the depth at which the B and C horizons appeared.
The surface soil horizon darkened extremely quickly from P1 to P3. This was due to the increase in the content of humic substances, following the addition of manure and crop residues, as reported by other authors [32]. The mineralization of these organic amendments, together with the mineral fertilizers provided by fertigation (nitrate, phosphates, sulfates, etc.) caused mottling in the Ap1, Ap2, and C1 horizons of the P2 and P3 profiles ( Figure A1A) and contributed to the vertical mobilization of salts in the profile. However, the spots in the C2 horizons of the profiles ( Figure A1B,C) have their origin in the genesis of the soil and correspond to CaCO 3 accumulations. These spots, together with the oxidoreduction spots ( Figure A1B) and other macromorphological and mineralogical characteristics, indicate specific conditions of hydromorphism, typical of soils with vertic or protovertic horizons and similar to those described by [33] in Vertisols in the province of Granada (Spain).
The TOC contents were low to moderate and decreased progressively with depth, as evidenced by [34] for other soils of the region. They increased as a consequence of the change in soil management, due to the incorporation of sheep manure, and also due to the application of nitrogen fertilizers, which caused an increase in the root biomass of the crops in the Ap horizons. These results are comparable to those obtained by [13]. Hence, there was a significant increase in TOC in the Ap1 horizon of the P2 and P3 profiles in comparison with P1 (33.3 and 21.7%, respectively). The authors of [19] found that organic fertilization can enhance the formation of organo-Ca complexes, which in turn improve soil aggregation and thus contribute to long-term TOC sequestration. The subsoiling to which the soil was subjected led to an accumulation of TOC in Ap2, the TOC content in this horizon being higher in P3 (significantly so) than in P2. Finally, the differences observed in the TOC contents of the B and C horizons among the three profiles must be related to the fluvial origin of this soil, rather than to a possible mobilization of soluble organic compounds.
The increased presence of C in the form of CaCO 3 in the deep horizons of the P2 and P3 profiles must have arisen from the incorporation of carbonates and bicarbonates with the irrigation water and the bioformation of CaCO 3 driven by the root respiration of the crops [35,36]. This CaCO 3 washing process was probably favored by the use of acidic fertilizer solutions, a very frequent practice in fertigation for cleaning irrigation equipment. Together with the application of manure, this would have led to a decrease in soil pHw and pHKCl, as reported by [37,38], among others.
These results confirm that the change in management practices favored C sequestration in this soil. Similar results were reported by [36,39,40] in the case of OCc, and by authors such as [41] and [35] in the case of ICc. All this reinforces the role of soils in arid and semiarid climates as C sinks and validates the compatibility of the management system used throughout the study with trends in sustainable agriculture [42,43], which the European Union is promoting in terms of soil protection [44]. The evolution over time of the C/N ratio, which is an indicator of the capacity of ecosystems to store and recycle energy and nutrients [45], towards values similar to those of the soil under the climax vegetation, also appears to validate the management practices.
The inputs of N, whether from manure or mineral fertilizers, contributed to an increase in biomass in the long term [46], in the form of the harvested crop or the crop residues, but at the same time also caused vertical leaching of N to deep horizons. This process was also observed by other authors [9,[47][48][49][50], who attributed it to tillage practices, in addition to excess irrigation and fertilization. If this is not remedied, it can lead to the contamination of the surrounding ecosystems.
In relation to salts and their dynamics in time and space, it should be pointed out that the increases in EC observed between 1998 and 2003 were due to the incorporation of soluble mineral fertilizers-such as sulfate, chloride, and nitrate-as well as organic amendments [9]. Similarly, the attenuation of EC and even its decrease in P3 may be explained by the decrease in sheep manure inputs, the introduction of other winter crops such as cereals and legumes (or even maize during the summer), and the adjustment of the fertigation programs [17,18,21]. The same factors would partly explain the rapid increases in K, Mg, and Na (from P1 to P2) and their subsequent declines (from P2 to P3), to which crop uptake and salt flushing would also have contributed.
The availability of other essential elements (P, Fe, Cu, Mn, and Zn) is also strongly affected by soil management and certain soil properties, such as pH and the amounts of CaCO 3 and organic matter [51,52]. Thus, the strong increase in available P over the study period is justified not only by the P content of the manure added, but also by the chelating effect of the organic matter present in the manure, which favors the mobilization of the non-available fraction of P present in carbonate rich soils [53]. In our work, the availability of Fe was also strongly influenced by the microenvironment, with results similar to those obtained by other authors [37,51,52]. These authors concluded that the concentration of available Fe in the soil is the result of the balance between antagonistic factors. A basic or slightly basic reaction of the soil and high concentrations of CaCO 3 (as in our work) cause immobilization of Fe in the form of insoluble salts and thus low Fe availability [54,55]. Conversely, organic matter favors the mobilization of Fe and similar metals through the formation of organometallic complexes [37,56,57]. A predominant role of the latter process would explain the increase in the availability of these essential elements over time in this work and their incorporation in available forms into the deep horizons. Because of this, it can be stated that the chemical fertility of the soil was significantly improved by the addition of organic amendments. This is very important for the good development of the crops, because in these soils with basic pH values and high limestone contents most of the divalent cations are present in precipitated forms and, therefore, very often the crops are deficient in these elements.
Crop rotation, especially since 2004, has increased the agronomic biodiversity [58], with the introduction of new crops such as cereals, broccoli, maize, or pepper and the use of oat-legumes mixtures as green manure, an objective promoted in the EU within its CAP. Finally, at the taxonomic level, the soil management implemented has led to a progressive evolution of the soil (Calcisol → Kastanozem) that is the opposite of that described by [59] for similar soils (Kastanozem → Calcisol). This progressive evolution of the soil, together with the other results already mentioned, validates the sustainable agronomic management model that has been implemented.
In summary, we consider that irrigation and the introduction of horticultural crops in soils that had traditionally been dedicated to rainfed cereal crops caused changes in the color and thickness of the upper horizon and, as a consequence, a taxonomic change. Most of the chemical properties analyzed underwent variations, while mineralogical properties were not affected. The results obtained confirm that soil degradation studies should be carried out in the medium and long term and should consider the entire soil profile, because, although anthropogenic activity takes place at the surface, the consequences are manifested both at the surface and in the deep horizons, the analysis of the latter being fundamental to reveal the processes of salinization and contamination. Based on the above, it can be affirmed that the hypotheses that gave rise to this research work were adequate.

Conclusions
Because of the management model implemented after the commencement of irrigation, between 1998 and 2003 an incipient salinization process was observed throughout the soil profile, which was even more intense in the Ck1 and C2 horizons than at the surface itself. Leaching of N was also detected, shown by the increase in the concentration of this constituent in the deep horizons. Macromorphologically, these processes were reflected in the soil profile by the appearance of white spots in the most affected horizons. However, from 2004 onwards, corrective measures implemented in the sustainable soil management model-with the introduction of winter crops, the use of green manures, and the adjustment of fertilization and irrigation-slowed these degradation and contamination processes. As a result, a positive evolution of the soil was observed, as manifested in its morphological characteristics and in some chemical properties. There was a darkening and an increase in the thickness of the A horizon, to the detriment of the underlying Bw horizon, and an increase in its fertility of the soil, and in the functional capacity of the soil as a carbon sink. Together, these have led to a change at the taxonomic level and to the attenuation of the incipient degradation detected in the first five years.
This study shows that the agronomic soil management carried out since 2004 constitutes a model of sustainable agriculture that can be extrapolated to regions with similar cultivation problems, and that most of the chemical and physical-chemical properties analyzed can be considered as sensitive indicators of the management model implemented. These properties should be monitored in all soil horizons to effectively monitor soil degradation and associated ecosystems. Clear and smooth limit.