Inﬂuence of Tillage and Crop Rotations in Organic and Conventional Farming Systems on Soil Organic Matter, Bulk Density and Enzymatic Activities in a Short-Term Field Experiment

: Intensive agricultural practices are leading to loss of soil fertility and overexploitation of natural resources which cause nutrients imbalance and further impair ecosystem services. Organic farming (OF), also coupled with minimum tillage and crop rotations, represents one of the strategies to limit this process and maintain soil functions. In a two-year ﬁeld trial, organic farming practices, including a set of fertilizations combined with crop rotations and association with nitrogen ﬁxing cover crops, were compared. The aim of this research was to assess in the short-term the effects on soil organic carbon, aggregate stability, and soil enzymes activities of using a combination of promising management practices in the delta region of the Po river. Results did not show improvements in organic carbon content and soil aggregate stability. Conversely, enzymatic activities were always signiﬁcantly higher in OF treatments than the conventional one. Crop rotation and associated legumes were effective in enhancing β -glucosidase and P ﬁxation through phosphatases activities. The present work suggests that an effective choice of crop species coupled with legumes can enhance biological activity re-starting main mechanisms of microbial development even without a contemporary increase of organic matter.


Introduction
The growing threats of climate change and unsustainable farming systems require a critical review of agronomic strategies needed to secure food production, while mitigating greenhouse gasses (GHGs) emissions, enhancing biodiversity, restoring quality of soil and water resources, and reducing soil erosion [1][2][3]. In seeking for the best agriculture practice, main issues to be addressed include: (i) depletion of soil organic matter (SOM), (ii) nutrient imbalance, (iii) accelerated soil erosion (with the oncoming risk of desertification), (iv) degradation of soil structure leading to, e.g., crusting and compaction, and (v) decline in soil water and nutrient retention capacities with reduction in resource use efficiency [4].
Organic farming (OF), which is legally defined in European Union (EU), has been proposed as a suitable alternative to conventional farming (CF) [5], with the ambition to drive the transition towards a sustainable food system that enhances the environment. Within the framework of the European Green Deal (European Commission, 2020) [6], the EU aims to achieve the target of 25% of agricultural land under OF by 2030. Organic farming relies on the SOM management to enhance soil physical, chemical, biological properties, in order to optimize crop production [7]. In OF systems, plant production temperature reaches minimum values in January (0.7 • C, on average) and maximum values in July and August (29.0 • C, on average). Air relative humidity ranges from 60 to 97%.
The characterization of the soil initial conditions for the experiment is reported in Table 1. The soil was a fairly homogeneous sandy loam classified as Calcaric Gleysol [25] characterized by a sub-neutral/alkaline pH, low SOC content, and low CEC.

Experimental Design
The experiment was conducted on plots, consisting each of 240 m 2 (8 m × 30 m), which were conducted in accordance to organic EC Regulation since 2000. Experimental design consisted of one conventional (Conv) and four organic treatments (Org) in a randomized block design with three replicates.
In the two years preceding the experiment (2006 and 2007), the area was managed as follows: the Conv plots were left uncultivated for 3 out of 4 semesters and for the remaining semester cultivated with radicchio (Chicorium intybus L., 1753); the Org plots were uncultivated for 2 semesters, then cultivated again with wheat (Triticum aestivum L.) and cabbage (Brassica oleracea var. cauliflower) and intercropped with clover (Trifolium repens L.).

Crops Yield of the Experimental Fields
The yield achieved in each Org management was recorded for the two-year experiment and total production is reported in Table 3  The data set shows the whole production of the entire experimental period (2008 and 2009). Clover biomass was not considered.

Soil Sampling
Disturbed soil samples were collected in the top 25 cm layer before starting the experiment (January 2008) for preliminary analyses reported in Table 1. Sampling on January 2008 was performed for initial soil characterization (main soil physical-chemical properties). During March, July, and November for both 2008 and 2009 years for a total of six sampling dates, samples for SOM and enzyme activities were collected. Field-moist bulk soil samples were split into two subsamples. One soil subsample was sieved at 2 mm and stored at 2 • C for enzyme activities analyses, while the other was used for chemical and physical analyses. Soil samples from each plot were composed of five subsamples that were bulked to obtain a sample of about 3 kg. Sampling involved the central area of the plots excluding a 2-m edge to avoid soil movement effects [26]. Additionally, a total of 30 undisturbed soil cores, 8 cm diameter, one per plot each year, in November 2008 and 2009, were taken from the topsoil (0-25 cm) for bulk density calculation.

Physical-Chemical Parameters
Undisturbed samples were weighed and a fraction was oven-dried at 105 • C for 24 h to calculate bulk density according to the core method [27]. The pH was measured potentio-metrically on pulverized samples by adding deionized water (1:2.5 w/v, dry weight basis). The electrical conductivity (EC) was potentiometrically determined after water extraction (2:5 w/v) and filtration through Whatman filter N. 2 paper. Particle size analysis was performed according to the hydrometer method, using sodium hexametaphosphate as a dispersant (Gee and Bauder 1986). Exchangeable bases, including calcium (Ca), magnesium (Mg), potassium (K), and sodium (Na), were extracted using barium chloride as reported in [28], and their concentrations were determined by inductively coupled plasma-optical emission spectroscopy (ICP-OES) on a SPECTRO CIROS (Spectro Analytical Instruments, Kleve, Germany). Cation exchange capacity CEC was calculated as reported in [29]. Plant available phosphorus (P) was calculated using the Olsen method [30]. Carbonate content (CaCO 3 ) was determined by the volumetric calcimeter method [31].
Soil organic C (SOC) was determined by a wet oxidation titration procedure using an acid dichromate system [32].
Total N was measured on air-dried 0.5 mm sieved samples using the flash combustion method (CNS Elemental Analyzer, Vario Max, Elementar Inc., Langenselbold, Germany).
For particle size distribution determination, the standard pipette method was used as described by USDA-NRCS (2004) [33]

Soil Organic Matter Parameters
Humic substances were extracted from 2 mm-sieved, air-dried samples with 0.1 M KOH (1:10 w/v) at room temperature for 16 h under an N 2 atmosphere and freed from the suspended material by centrifugation at 7.000× g for 20 min and filtration on Whatman (Whatman, Maidstone, England) N • 42 filter paper. Organic carbon content was assayed by dichromate oxidation [32]. Molecular-weight distribution and gel-permeation chromatography of each humic extract was carried out on a Sephadex G-100 gel packed in a 70 cm × 1.6 cm Pharmacia column (Pharmacia, Uppsala, Sweden) as reported in Cardinali et al. [34]. The gel packing solution and the eluant were 20 mM Na 2 B 4 O 7 . The apparent molecular weight of the fractions was >100, 10-100, and <10 kDa. The column calibration was based on standard kit for molecular weights (Sigma-Aldrich Gel Filtration Molecular Weight Markers MWGF200). All chemical determinations were performed in triplicate.

Enzyme Activities
Enzyme activities were determined on fresh moist sieved (<2 mm) soils within 10 days from the collection of the samples. Triplicates were performed for each activity assay. Control assays were performed in all cases, except for dehydrogenase, by adding the substrate after the reaction was stopped and before filtration of the soil suspension. All data were corrected to soil dry weight by means of soil moisture content data.
Urease activity (E.C. 3.5.1.5) was assayed as described by Kandeler and Gerber [35] using 0.08 M urea as substrate while the samples were incubated at 30 • C for 2 h. N-αbenzoyl-L-argininamide (BAA) hydrolyzing protease activity was determined in 0.1 M phosphate buffer at pH 7; 0.03 M BAA was used as substrate [36]. The sample was incubated at 40 • C for 2 h. The resulting suspensions were filtered, and filtrates analyzed for ammonia by the following colorimetric procedure. One milliliter of filtrate was diluted to 10 mL with distilled H 2 O and, successively, 5 mL Na salicylate solution (100 mL sodium nitroprussiate 1.2 g L −1 , 100 mL Na salicylate 17 g L −1 , and 100 mL distilled H 2 O mixed and the moment of use) and 2 mL 10% Na dichloroisochanurate were added. Solutions were incubated for 30 min at room temperature and solution absorbance measured spectrophotometrically at 600 nm [35]. Urease and protease activities were determined as µg NH 4 + g −1 dry soil h −1 released in the hydrolysis reaction.
Dehydrogenase activity was determined following a modification of Masciandaro et al. [37]. Soil was mixed with 0.2 mL of 0.4% (7.91 mM) 2-p-iodo-phenyl-3p-nitrophenyl-5-phenyltetrazolium chloride (INT) solution (in distilled water, w/v). The control was the soil treated with distilled water (0.2 mL), instead of INT. Soils were incubated for 20 h at 20 • C in darkness; pH values did not vary during the incubation time. The iodonitrotetrazolium formazan (INTF) formed by the reduction of INT was extracted by adding 10 mL of a mixture of 1:1.5 tetrachloroethylene and acetone; then, the soil mixture was vigorously shaken by hand for 1 min and finally filtered through a Whatman no. 41 filter paper. The INTF concentration was measured spectrophotometrically at 490 nm, and the results were expressed as µg INTF g -1 dry soil h -1 .
Phosphodiesterase activity was determined following Browman and Tabatabai [38]. Briefly, 4 mL 0.05 M TRIS-buffer pH 8 and 1 mL bis-p-nitrophenyl phosphate 0.005 M were added to 1 g of soil and incubated at 37 • C for 1 h. The reaction was stopped by adding 1 mL 0.5 M CaCl 2 and 4 mL 0.1 M TRIS-NaOH pH 12 buffer.
For the determination of β-glucosidase activity, four milliliters of modified universal buffer (MUB), pH 6, and 1 mL of substrate p-nitrophenyl-b-D-glucopyranoside (PNG) 0.025 M were added to 1 g of soil and incubated at 37 • C for 1 h The reaction was stopped by adding 1 mL 0.5 M CaCl 2 and 4 mL 0.1 M TRIS-NaOH pH 12 buffer [39].The p-nitrophenol (pNP) formed in alkaline phosphatase, phosphodiesterase, and β-glucosidase activities was determined by spectrophotometry at 400 nm and reported as µg p-nitrophenol (pNP) g −1 dry soil h −1 .
The fluorescein diacetate hydrolase (FDAH) activity was assessed as described by Adam and Duncan [40]. A 2-g soil sample was shaken for 15 min with 15 mL of 60 mM potassium phosphate pH 7.6 buffer and with 0.2 mL of 1000 µg mL −1 FDA solution for 20 min at 30 • C. 15 mL of 2:1 chloroform/methanol solution was added, and the contents shaken thoroughly by hand before filtration (Whatman, No 2). The absorbance was measured spectrophotometrically at 490 nm on the filtrate. The concentration was expressed in mg fluorescein released g −1 dry soil h −1 .

Soil Aggregate Stability
Samples were air-dried. During the drying operation, the bulk soil was broken along natural fissures into small pieces by gentle hand manipulation [40]. Once dried, the 1-2 mm aggregate fraction was separated and subjected to 200-µm wet sieving, applying vertical shaking for 30 min with a frequency of 30 oscillations per minute (900 oscillations in total) at a 3 cm amplitude [41]. The material was then carefully transferred to a second beaker by a washing flask and oven-dried at 105 • C to constant mass. The coarse sand (200-2000 µm) content was determined by adding 25 mL of a solution containing 40 g sodium hexametaphosphate and 10 g sodium carbonate in 1 L to each dry sample (about 7 g), allowing to stand for 24 h, washing (with demineralized water) on a sieve (200 µm), and oven-drying the material retained.
The wet aggregate stability index (WASI) was calculated as specified in Dal Ferro et al. [13]: where WSA is the dry mass of water stable aggregates >200 µm, W agg the total dry mass of the sieved aggregates, and W sand the dry mass of coarse sand particles measured after dispersion. Before wet sieving, soil aggregates were subjected to fast wetting obtained by immersing 10 g soil aggregates in deionized water for 10 min to simulate the effect of heavy rain storms [42]. To distinguish slaking from physical-chemical dispersion, before fast wetting soil aggregates were subjected to: (i) no pretreatment (WSA), (ii) pretreatment with benzene (BSA = benzene stable aggregates) to coat the aggregates and emphasize the hydrophobic nature of SOM [12], (iii) pretreatment with ethanol (ASA = alcohol stable aggregates) to test the wet mechanical cohesion by substituting air with a polar liquid before energy is applied [43]. Immersion of ethanol-pretreated soil aggregates in water minimizes slaking of dried soil aggregates and emphasizes the resistance of soil aggregates to dissolution and the dispersive action of water.

Statistical Analyses
One-way ANOVA with repeated measures were performed in order to compare treatments among the sampling periods; for those parameters that have available less than 3 sampling times, a two-way ANOVA was performed, including the sampling time as a fixed factor. Duncan test was performed for post-hoc comparisons. Paired T test was performed to compare enzymatic activities between the same period of different years (e.g., March 2008 vs. March 2009). Levene and Mauchly's tests were applied to check homoscedasticity and sphericity, respectively, to ensure that assumptions of the model were met as recommended in Field [44]. Pearson's correlation was performed to test correlations between organic carbon content and soil stability indexes. All these procedures were performed by IBM SPSS 26.0.
A canonical discriminant analysis (CDA) was performed by SAS software (SAS Institute, 2004) to determine which variables best predicts group membership and to visualize the data by grouping the variables in one or more axes [34]. The analysis selects a limited set of discriminant functions that best describe the separation of the studied plots. A stepwise procedure was also applied to choose the set of variables characterized by the highest correlation with the factor and assess the relative importance of each independent variable. The consistency and overall robustness of the separation was evaluated by the Mahalanobis quadratic distance, an index that measures the squared distances between the centroids of each cloud, thus reflecting the actual separation between the different managements. The Wilks Lambda, performed to evaluate the quality of the multivariate statistics, shows significant differences among population means.

Soil Physical Parameters
Bulk density (Table 4) showed significant differences only in 2008 between Org PBF and the other treatments. It was also noted a convergence to a narrow range of bulk density values in 2009 compared to the previous year. Although this behavior appears clear, from a statistically point of view, there is no evidence about a significant effect of the factor time. Wet aggregate stability indexes (WASI) (Table 5A) pretreated with ethanol (ASA) led to the highest aggregate stability in both years, showing a 48.6% biennial mean value, while benzene pretreatment (BSA) and non-pretreated aggregates (WSA) showed lower and similar values, respectively 8.1% and 12.8%. The aggregate structure was not degrading significantly in the different years, despite a tendency to lower WASI was observed in 2009 than 2008, regardless the tested treatments.  WASI indexes were generally significantly lower in Conv WR than Org managements, especially after ASA, resulting over the years consistently <30% compared to values >45% in other treatments.

Soil Organic Carbon and Humic Fraction Content
Soil organic carbon content (Table 5B) was 9.5 g kg −1 in March 2008, on average, afterwards it decreased to 8.5 g kg −1 (July 2008) which was roughly maintained until the end of the experiment (November 2009). The only significant difference was observed between Conv and Org treatments, the former showing lower values (5.2 ± 0.4 g kg −1 ) than Org managements (10 ± 0.2 g kg −1 ) across the years.
Humic matters were characterized by a decrease of high (F1) and low (F3) apparent molecular weight fractions and an increase of intermediate ones (F2), regardless the agricultural system management (Table 5C). Time seems to be the most important factor in humic matter evolution (p < 0.01) with no interaction with management. Differences among treatments were statistically significant only for F1 and F3, particularly emphasized by a lower F1, and higher F3, in Org PBF than other organic or conventional treatments.
While the humic matter evolution process was affected by agricultural management systems, the humic carbon content (Table 5D) did not differ significantly across the tested dates. It showed a common pattern among all the organic treatments, with slight increases of humic carbon (HC) content in July and November 2008 and a following decrease, while in Conv WR content did not show significant variations.

Interaction between Soil Organic Matter and Soil Stability
Soil organic carbon and HC showed significant positive correlations with all aggregate stability indexes. Ethanol-pretreated aggregates (ASA) were the most correlated with SOC (r = 0.80) and HC (r = 0.73), followed by BSA (r = 0.63 for both SOC and HC) and WSA (r = 0.59 for both SOC and HC), confirming the role of organic matter in giving stability to soil aggregates (Table 5A). The high molecular weight humic fraction (F1) was the most positively correlated with aggregate stability indexes, in particular with BSA (r = 0.70) and ASA (r = 0.73) that emphasized the humic matter hydrophobicity and ability to enhance chemical bonds with soil mineral particles (Table 5B).

Enzymatic Activity
Average values of conventional management were always significantly lower than those of organic managements (M average). Among organic managements, dehydrogenase, phosphodiesterase, protease, and urease activities were similar (Tables 6-9, respectively) while phosphomonoesterase, FDA, and β-glucosidase (Tables 10-12, respectively) activities were higher in thesis with legumes employed as cover/intercrops.
Among the seasons and years (T average), FDA and dehydrogenase showed highest values in summer of both years while β-glucosidase showed a regular increase from 2008 to 2009 with a peak in July 2009; the two phosphoesterases presented a peak in summer with an irregular trend across the years; protease and urease were constant with a peak in July 2008.
Within the seasons, except for FDA and β-glucosidase that presented significant differences between March 2008 and March 2009 for Org SGR and Org SGR + C, respectively (Tables 11 and 12), all treatments showed significant differences only between July 2008 and 2009, and/or between November 2008 and 2009. Dehydrogenase activity remained stable in all the situations with the unique difference between CF and OF.

Canonical Discriminant Analysis
CDA analysis was performed (1) to confirm the differences in soil fertility parameters between CF and OF systems and (2) to evaluate the effects of diverse crop rotations on soil biochemical properties.
The scatter plot (Figure 1) of the first two extracted canonical roots showed that CAN 1 (accounting for 80.09% of the model variability) well discriminated Conv WR and Org SGR from Org PBF + P and Org PBF, whereas CAN 2 (explaining 17.48%) discriminated Org SGR from the other managements. The canonical scores are reported in supplementary data. The overall Wilk's Lambda was 0.0074 (F = 10.01; p < 0.0001). The Mahalanobis quadratic distances of the 5 clouds and their significance (Tables A1 and A2) are reported in Appendix A. Values indicated that only Org PBF + P and Org PBF were not significantly different. Among the seasons and years (T average), FDA and dehydrogenase showed highest values in summer of both years while β-glucosidase showed a regular increase from 2008 to 2009 with a peak in July 2009; the two phosphoesterases presented a peak in summer with an irregular trend across the years; protease and urease were constant with a peak in July 2008.
Within the seasons, except for FDA and β-glucosidase that presented significant differences between March 2008 and March 2009 for Org SGR and Org SGR + C, respectively (Tables 11 and 12), all treatments showed significant differences only between July 2008 and 2009, and/or between November 2008 and 2009. Dehydrogenase activity remained stable in all the situations with the unique difference between CF and OF.

Canonical Discriminant Analysis
CDA analysis was performed (1) to confirm the differences in soil fertility parameters between CF and OF systems and (2) to evaluate the effects of diverse crop rotations on soil biochemical properties.
The scatter plot (Figure 1) of the first two extracted canonical roots showed that CAN 1 (accounting for 80.09% of the model variability) well discriminated Conv WR and Org SGR from Org PBF + P and Org PBF, whereas CAN 2 (explaining 17.48%) discriminated Org SGR from the other managements. The canonical scores are reported in supplementary data. The overall Wilk's Lambda was 0.0074 (F = 10.01; p < 0.0001). The Mahalanobis quadratic distances of the 5 clouds and their significance (Tables A1 and A2) are reported in Appendix A. Values indicated that only Org PBF + P and Org PBF were not significantly different.  Variables that significantly contributed to the CDA model were the enzymatic activities β-glucosidase, phosphomonoesterase, phosphodiesterase, the humic carbon content, and the soil physical properties of aggregate stability (ASA) and bulk density (Figure 1). These variables were not equally important in the construction of the canonical functions, as demonstrated by Pearson correlation coefficients and significance levels (Table A2). Although, contributing to the model, bulk density values did not correlate with the canonical axes, while all other factors significantly correlate with the two canonical functions.

Soil Organic Matter and Stability Indexes
Since the first sampling to the end of the experiment, Conv WR showed significantly lower topsoil content of both SOC and HC (Table 5B). These results were likely due to the combined effect of different fertilization (mineral in Conv vs. organic in Org) and tillage operations (intensive in Conv vs. minimum in Org), that were primarily involved in SOC change dynamics as also observed in other studies [45]. However, all treatments showed SOC and HC reduction, highlighting the poor ability of coarse texture soils to form stable organo-mineral complexes that would protect the organic matter from mineralization. These results corroborated previous findings in similar pedo-climatic conditions, where different soil improving cropping systems were compared for SOC stock in a 6-year experiment, and results revealed that they suffered from low soil reactivity to SOC accumulation dynamics [46]. These dynamics are in contrast with long-term studies where SOC accumulation dynamics occurred in different soil types, either clayey, loamy, or sandy [47]. In our experiment, the topsoil SOC depletion could be related with: (i) the short time of the experimental monitoring; (ii) the effects of soil disturbance and SOC mineralization with minimum tillage that counteracted the topsoil organic carbon accumulation with no tillage in previous years in the same plots. The latter seems in agreement with previous studies [48,49], where a topsoil SOC content was partially lost just after some greater intensity in soil tillage operations, suggesting that aggregate breakdown stimulation and reduced SOC protection occurred [50].
The results about WASIs corroborated previous findings on the strong relationship between SOC content and structure [51,52]: here, a trend towards a loss of aggregate structure stability was observed from early 2008 to 2009, which was consistent with the SOC and HC depletion. In addition, significant WASIs reductions were observed when more intensive tillage operations-under Conv-were compared with minimum tillageunder Org.
The wet mechanical cohesion between soil and organic matter particles assumed particular importance, being ASA stability indexes always greater in Org compared to Conv. In contrast, the BSA index had little differentiation between treatments, suggesting that the hydrophobic behavior of soil organic matter to protect the structure was likely a secondary mechanism.
Humic matter hydrophobicity may not be the most significant factor in improving aggregates stability [38], this is consistent with our data: although with the clear evolution of humic substances (Table 5B) towards the intermediate molecular fraction, there is no evidence of increased aggregate stability through WASI tests (Table 5A); those, conversely, showed a general decrease of stability in agreement with OC and HC decrease. On the other side, the clear correlation among WASIs indexes (especially ASA) and high molecular weight fraction (Table 13) can suggest that F1 can exert an influence on the chemical bonding that affects aggregates stability. This is corroborated by previous findings [9] that reported strong positive correlations (r > 0.90) between F1 and carboxyl groups, which in turn are responsible for strong soil-SOM chemisorptive bondings, especially in soils with low SOC content <1% [14]. We argue that the decrease of bulk density (Table 4), which could be a signal of enhanced porosity, was the consequence of conventional and light ploughing that broke down the superficial layer after the no cultivation period. A slow reaction to improve soil structure properties in similar silty soils with low organic carbon content was also observed in previous studies, [53], suggesting that long-term strategies to increase the chemical and physical soil fertility are required.

Enzymatic Activities
All enzymes activities considered in this work show higher values in organic treatments than conventional one but with noticeable differences among themselves.
Fluorescein diacetate (FDA) hydrolysis and β-glucosidase are indicators of total enzymatic activity due to the capability of many enzymes free or bound to cell membranes to hydrolyze fluorescein diacetate, and the presence of β-glucosidase in circulant solution or in fungi: agricultural practices involving minimum soil disturbance, maximizing soil cover, and the use of organic amendments can contribute to improve β-glucosidase activity through the increase of fungi presence and organic C but FDA also displayed modifications caused by climate [54][55][56]. Our data showed a peak of activity for FDA in correspondence with the warmer periods of the experiment (July 2008 and 2009) while β-glucosidase showed also peak in November 2009. Interestingly, our data do not match with works cited above about direct relationship between C increase and enzymatic activity increase as we recorded a general OC and HC decrease, in particular from March 2008 to July 2008 and a following steady state until the end of the trial. Although in this trend, enzymatic activity of FDA hydrolases and β-glucosidase increased in all thesis except Conv WR for β-glucosidase. An explanation could be one as described in Weerasekara [57]: under the legume hairy vetch and cereal rye employed as cover crops in silty loam or sand, β-glucosidase increased its activity while FDA hydrolases significantly decreased. This is consistent with thesis Org SGR + C that employs a cereal (sorghum) as one of main crops plus a legume as cover and associate crop (clover) in a sandy loam soil: not only average β-glucosidase values are the highest but it shows the higher increase between the same seasons of the two years (Table 13). Conversely, FDA activity increased and remained higher in thesis Org SGR + C and Org PBF + P characterized by cereals + legumes (Table 11).
Plant diversity can increase microbial activity and carbon storage as result mainly of new carbon integration by microbial community [58] but there seems the need of the "right combination" of crops, soil, management, and climate in order to achieve some specific results in term of biological activity.
In our work, dehydrogenase activity did not show any significant change within each thesis during the experimental period (Table 6) nor among organic managements. In the present work, this enzymatic activity seems to depend on the organic carbon content that remained roughly stable in the present experiment, after a first decrease. Previous long-term data confirmed the relation among dehydrogenase activity and OC content [59]. Lange at al. indicated that C storing in soil is mainly the effect of new C integration by the microbial community [58] and that the dehydrogenase activity is strictly related to microbial biomass and not only to total organic matter content [60]. It is reasonable to hypothesize the conditions of our experiment were not enough to cause an increase of dehydrogenase activity in two years, although data on FDA indicates that organic managements were able to increase the microbial community.
Phosphodiesterase and phosphomonoesterase activities are known to act sequentially in soil. It has been reported that the lower activity of phosphodiesterases than acid and alkaline phosphomonoesterase causes the production of P monoesters from P diesters that may stimulate the microbial synthesis of more phosphomonoesterases [61]. This conclusion seems in agreement with our data where phosphomonoesterase shows an activity up to four times higher than phosphodiesterase. Phosphodiesterase activity is significantly higher in organic managements characterized by light tillage than in conventional tillage which is consistent with previous trials that reported a higher phosphodiesterase activity in light tillage + mulching than in conventional and no tillage [59,60,62]. Phosphodiesterase activity was in average higher in 2009 than in 2008: due to the stable OC content, the main cause of this increase could be the rotation and in particular the combination cropassociated legume. The clover-garlic (CG) combination resulted in the highest values in summer among all managements and within the Org SGR + C thesis, this combination gave always the highest results. The presence of garlic seems to be correlated with higher phosphodiesterase activity in Org SCG also when not associated with clover. The positive action of garlic root exudates was reported also in short-term continuous cropping (5-10 years) characterized by a significant increase in bacteria, fungi, phosphatase, and urease activity in the rhizosphere [63]. A positive effect on soil enzyme activities of garlic was also reported by Zhou at al. [64] for soil urease activity; Wang et al. [65] for urease and alkaline phosphatase activities. Yiang et al. [66] reported that garlic increased soil enzymatic activities and fungi populations also when rotated with watermelon. In both summer 2008 and 2009, the phosphodiesterase activity (Table 7), although without statistic evidence, seems to be higher in theses managed with legumes than in those organic theses without clover and pea: this is consistent with [67] who found that nodulated legumes species growing in savanna soils stimulated the acid phosphatase activity of the rhizosphere soil.
Although, phosphodiesterase and phosphomonoesterase act sequentially [52], when comparing the same period within the two years of experiment, we found phosphodiesterase increased in 2009, while phosphomonoesterase decreased (Table 10). We hypothesize the soil was enriched in organic and inorganic P through the years, causing a decrease of phosphomonoesterase activity. This is in agreement with Kandeler et al. [68] who found an inverse correlation between the phosphomonoesterase activity and the P organic content of the rhizosphere soil sampled from Triticum aestivum and Trifolium alexandrium.
Protease activity evidenced a decrease from 2008 to 2009 except for the march datapoint (Table 8), with the higher values in July 2008 and July 2009 when samples were collected during fertilizers addition and higher temperatures: the sudden increase in C and N availability could have increased the microbial community growth rate with a consequent higher need of N. As reported by Geisseler and Horwath [69], the level of protease synthesis is determined by the need for C rather than N: this could be in agreement with OC decrease during the experiment.
Urease activity followed the same trend of protease (Table 8): most authors, see for example [70], found an increase of urease activity correlated with zero/light tillage and the consequent increase of organic matter as confirmed in the present work on comparison with conventionally managed soils. For the year-to-year trend among the organic samples, we hypothesize, like for protease, a disbalance in C and N content. However, it is interesting to note that, although not significant, urease activity is higher in bare soil when compared with same soil covered by crop ( Table 9).
The canonical discriminant analysis successfully discriminates conventionally managed soils and organic managed theses. This is in line with the known action of organic management in soil activities [71]. The analysis also distinguishes Org SGR and Org SGR + C from Org PBF treatments with or without pea intercropping and/or associated crop. This confirms the importance of the crops in determining the overall soil quality. Present results indicate that garlic plays an important role in these phenomena in association with the presence of leguminous crop as clover. Not all the factors were equally important in the CDA algorithm ( Figure 1, Table A1); this is in agreement with our analysis: indeed phosphodiesterase, phosphomonoesterase, and β-glucosidase gave the main contribution in separating the 4 clusters. In particular, the two phosphatases showed the main role in separating groups being the two enzymes apparently not dependent on C content decrease/stability, probably due to the priority to store P to build ATP.

Conclusions
In a sandy loamy soil, characterized by low OC content, organic management, coupled with minimum tillage, maintained a higher level of biological activity in comparison with conventional farming systems. Due to the short-term period of the experiment, the main physical parameters depending on organic matter content, like bulk density, did not significantly change. In contrast, the soil aggregate stability seemed affected by the quantity and quality of organic matter, suggesting that tested OF systems are promising for improvements of soil physical fertility in the long-term.
Crop rotation and the association with a legume gave the higher contribution in enhancing the biological activity in soil increasing "phosphatase mechanism" as a main biological source of energy. Garlic showed the best results in term of phosphatases activity.
Microbial phosphomonoesterases are probably the major contributing enzymes to the measured enzyme activity in soil [72,73]. This work suggests that the right choice of an effective rotation with the right crop species is fundamental to re-start biological activity in impoverished soil or to rapidly ameliorate its structure. Finally, from a methodological point of view, this work confirms the pivotal role played by the enzymatic indexes for monitoring the effects of management systems, being sensitive soil attributes that reflect the capacity of a soil to function and can be used as indicators of soil quality also in the short term.