Next Article in Journal
Enriched Co-Treatment of Pharmaceutical and Acidic Metal-Containing Wastewater with Nano Zero-Valent Iron
Previous Article in Journal
Mineralogical and Chemical Specificity of Dusts Originating from Iron and Non-Ferrous Metallurgy in the Light of Their Magnetic Susceptibility
Previous Article in Special Issue
Evaluation of a Multi-Isotope Approach as a Complement to Concentration Data within Environmental Forensics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Soil Biochemical Indicators and Biological Fertility in Agricultural Soils: A Case Study from Northern Italy

by
Livia Vittori Antisari
1,2,
Chiara Ferronato
2,
Mauro De Feudis
1,2,*,
Claudio Natali
3,
Gianluca Bianchini
4 and
Gloria Falsone
1,2
1
Department of Agricultural and Food Sciences, Alma Mater Studiorum-University of Bologna, Via Fanin 40, 40127 Bologna, Italy
2
Centro Sperimentale per lo Studio e l’Analisi del Suolo (CSSAS), Alma Mater Studiorum-University of Bologna, Via Fanin 40, 40127 Bologna, Italy
3
Department of Earth Science, University of Florence, Via La Pira 4, 50121 Firenze, Italy
4
Department of Physics and Earth Science, University of Ferrara, Via Saragat 1, 44122 Ferrara, Italy
*
Author to whom correspondence should be addressed.
Minerals 2021, 11(2), 219; https://doi.org/10.3390/min11020219
Submission received: 26 January 2021 / Revised: 16 February 2021 / Accepted: 18 February 2021 / Published: 20 February 2021
(This article belongs to the Special Issue Elemental and Isotope Geochemistry of the Earth’s Critical Zone)

Abstract

:
Industrial farming without considering soil biological features could lead to soil degradation. We aimed to evaluate the biochemical properties (BPs) and biological fertility (BF) of different soils under processing tomato cultivation; estimate the BF through the calculation of a simplified BF index (BFIs); determine if the crop was affected by BP and BF. Three farms were individuated in Modena (MO), Ferrara (MEZ) and Ravenna (RA) provinces, Italy. Soil analysis included total and labile organic C, microbial biomass-C (Cmic) and microbial respiration measurements. The metabolic (qCO2), mineralization (qM) and microbial (qMIC) quotients, and BFIs were calculated. Furthermore, plant nutrient contents were determined. The low Cmic content and qMIC, and high qCO2 found in MEZ soils indicate the occurrence of stressful conditions. The high qMIC and qM, and the low qCO2 demonstrated an efficient organic carbon incorporation as Cmic in MO soils. In RA soils, the low total and labile organic C contents limited the Cmic and microbial respiration. Therefore, as confirmed by the BFIs, while MO showed the healthiest soils, RA soils had an inefficient ecophysiological energy state. However, no effects on plant nutrient contents were observed, likely because of masked by fertigation. Finally, BP monitoring is needed in order to avoid soil degradation and, in turn, crop production decline.

1. Introduction

Anthropogenic actions deeply influence processes occurring in the Earth’s critical zone (CZ), and a typical example is what is observed in agroecosystems where soils are processed for crop production. In fact, soil features reflect parent rock and climatic characteristics but also the effect of biological activities, including those triggered by man-made actions. These features are important to understand the soil fertility, i.e., the soil attitude to sustain crop productivity, a characteristic that has to be maintained to carry sustainable agricultural practices. In this view, studies must be carried out with a multidisciplinary approach and there is a need to develop analytical methods that measure biogeochemical indexes for a detailed CZ understanding.
In industrialized countries, a large part of the soil suitable for cultivation is used for intensive crop production [1]. Agricultural production must, however, consider its sustainability in order to reduce soil degradation processes, such as soil erosion, loss of fertility and biodiversity, accumulation of agrochemicals and salinization [2,3,4,5]. To date, conventional and industrial agriculture are considered unsustainable because they are eroding natural resources (i.e., soil and water) faster than the capacity of ecosystem to regenerate them.
Soil microbial biomass and its activity are driving forces for the cycling of nutrients in agricultural and natural ecosystems, soil processes and plant performance [4,6], and soil microbial biomass, together with soil organic carbon, is largely related to different soil ecological functions (e.g., sustain crops, accumulate and filter water, act as carbon sink).
Soil variability and specific environmental conditions may, however, control changes in soil processes and, consequently, soil microbial variations [6]. In fact, the considerable effects of pedological and climatic conditions on soil microbial biomass and activity have been documented [7,8].
The prompt response of soil microbial populations to changes in soil properties and land use (e.g., [8,9,10,11]) pressed the scientific community to focus on microbial biomass and activity as accurate and sensitive indicators of the fitness of soil to perform ecological functions [12]. In this context, several papers used soil biochemical properties to monitor soil fertility, (e.g., [12,13]), and degradation, (e.g., [14,15,16]), in agricultural lands. Widdig et al. [17], in an experiment conducted on six grassland sites from South Africa, the USA and the UK, found a reduction in soil microbial respiration and growth in soils added with nitrogen. Nunes et al. [18], in an investigation conducted in lands with different degradation degree in Northeast Brazil, showed that soil degradation, through the organic C loss, reduces both soil microbial biomass and activity. However, although preservation and improvement of soil biochemical parameters are signs of better soil ecological functions, in agricultural ecosystems, these ecological indicators have to face up with crop yields. From this point of view, previous studies highlighted the pivotal role of soil microbial community to sustain plant growth and yield [19,20,21,22].
Thus, the strong connection between soil edaphic properties and crop production through the intermediary role of the soil microbial community would suggest the use of biochemical properties as a useful tool to evaluate if a soil can be suitable for the cultivation of a given crop without having a negative impact on soil ecological functions. In this sense, we hypothesized that (i) the biochemical properties and related index of biological fertility are prompt and valuable indicators of soil stress conditions caused by industrial agriculture; (ii) soils with low biological fertility could negatively impact crop yield and quality.
Therefore, the aims of this study were (i) to evaluate the biochemical properties (i.e., microbial biomass and activity) and biological fertility of soils with contrasting characteristics under processing tomato cultivation; (ii) to estimate the soil biological fertility through the calculation of a simplified biological fertility index; (iii) to determine if the quality of the crop production (i.e., nutrients content in plants) was affected by soil biochemical properties and biological fertility.

2. Materials and Methods

2.1. Study Areas and Experimental Layout

The selected study areas were located in the eastern Padania Plain, and specifically in Modena (MO sites), Ferrara (MEZ sites) and Ravenna (RA sites) provinces of Emilia Romagna region, Italy (Figure 1). These sites were selected for their high representativeness of the Padania Plain of Emilia Romagna region, which is characterized by lands addressed for industrial crops.
The MO soils are located in the northernmost area and developed on silt and clay sediments of the Po River in the area between river course and the Apennines chain; the soils are classified as Irragric Vertisols [23]. The MEZ soils are located in a lowland reclaimed for agricultural purpose in the 1950s and are characterized by the common presence of peaty outcrops; the soils are classified as Thionic Sapric Histosols [23]. The RA soils are in the southernmost area consisting of silt and sand fluvial deposits flowing from the Apennines chain; the soils are classified as Endogleyic Fluvic Cambisols [23].
Temperature and precipitation were monitored during 2015 cultivation period through three climatic stations located in Guagnino (MEZ), San Felice sul Panaro (MO) and San Pietro in Vincoli (RA). For the considered year, the total annual precipitation was 621, 630 and 909 mm for MO, MEZ and RA, respectively. In the three sites, the minimum average annual temperature was recorded on February and ranged between 1.9 and 2.6 °C. The maximum average annual temperature was recorded in July and ranged between 22.8 and 24.0 °C. Finally, the mean annual temperatures were 13.3, 13.3 and 12.9 °C for MO, MEZ and RA, respectively. The climate data are in agreement with those recorded in the fifteen-year period (2001–2015) by the regional service in the climatic station located in Guagnino and Finale Emilia for MEZ and MO sites, respectively, and San Pietro in Trento and San Pancrazio for RA site [24].
For each study area, farms with similar land management were selected. In particular, 3 farms for MO, 3 for MEZ and 2 for RA were individuated. In each farm, 3 processing tomato fields (2.5 ha each) were chosen. In the fields, the tomato (Solanum lycopersicum L.) crop, which is a widely distributed crop of the three study areas, has been cultivated as monoculture since 2014 and was previously used for maize cultivation.
In the investigated year, the fields were thoroughly rototilled at a depth of 20 cm and then the tomato seedlings (hybrid Heinz 1015 F1) were transplanted within the first week of May (with a few days of variation among the fields). The planting density was 45,000 plants per hectare with row spacing of 0.5 m and plant spacing of 0.4 m. During the growing season, no more tillages were performed, while weed and pest controls were carried out according to the integrated crop management specifications of Emilia Romagna Region [25]. Plants were watered through drip irrigation lines and nutrients were provided through fertigation. Roughly, the amounts of nitrogen, phosphorus and potassium provided to the soils were 236, 167 and 62 kg ha–1, respectively, for MO, 260, 160 and 150 kg ha–1, respectively, for RA and 94, 117 and 174 kg ha–1, respectively, for MEZ (Table S1 of the Supplementary Materials). The higher amounts of nitrogen provided in MO and RA compared to MEZ, while the higher amounts of potassium provided in RA and MEZ comparing to MO were to counterbalance the lower amounts of those nutrients in the selected soils.

2.2. Soil and Tomato Plant Sampling

Two weeks after the tomato transplanting, in each field, soil sampling was carried out by opening two minipits down to 20 cm (i.e., soil samples at 0–20 cm) and by drilling the bottom part of the minipit down to 40 cm by auger (i.e., soil samples at 20–40 cm). For each field, the samples collected at each soil depth were mixed to obtain composite samples of 0–20 and 20–40 cm.
At harvesting, five randomly selected tomato plants and their rhizospheric soil (the soil still adherent to the roots after gently shaking and removing the loosely attached soil) were collected for each field. Once in the laboratory, the rhizospheric soil was separated from the plant roots and the organs of tomato plants (i.e., young and old leaves, stems, roots and tomato fruits) were separated from each other.
The composed 0–20 and 20–40 cm soil samples and the rhizosphere soils were air dried, sieved to 2 mm and an aliquot was finely ground. The tomato fruits were freeze- dried, while the other plant parts were dried at 60 °C for 48 h. All plant samples were then finely ground.

2.3. Soil and Tomato Plant Analyses

2.3.1. Soil Physicochemical Properties

The pH was determined potentiometrically in a 1:2.5 soil:water suspension (w/v). The electrical conductivity (EC) was measured in a 1:5 soil:water suspension (w/v). The cation exchange capacity (CEC) and the exchangeable cations were determined after exchange with hexamminecobalt(III) chloride and measuring the excess of Co3+, and Ca2+, Mg2+, K+ and Na+ in the extracts by inductively coupled plasma–optical emission spectrometry (ICP-OES, Ametek, Germany) [26]. The carbonate content was quantified by volumetric method through the reaction of soil sample with 6 M HCl [27]. The particle-size distribution was measured by pipette method [28]. The amount of total organic carbon (OC) and nitrogen (TN) was determined by dry combustion (EA 1110, Thermo Fisher, Waltham, MA, USA) after removing inorganic carbon by HCl. The amount of stable carbon isotope (13C) was determined by an isotope ratio mass spectrometer (IRMS Delta C or DELTA+XL, Thermo Finnigan MAT, Bremen, Germany), and it was then expressed as δ13C (‰) with respect to the V-PDB universal reference standard for 13C. Total P (TP) was determined after treating samples with suprapure HCl and HNO3 mixture (3:1 v/v) in microwave oven (Millestone 1200, Milestone Inc., Sorisole, Bergamo, Italy) [29]. The available P (POlsen) was extracted using 0.5 M NaHCO3 at pH 8.5 [30]. The amount of TP and POlsen was then measured by ICP-OES (Ametek, Germany).

2.3.2. Soil Biochemical Parameters

The C and N of microbial biomass (Cmic and Nmic, respectively) were determined as the difference between C and N extracted with 0.5 M K2SO4 (1:4 soil:solution ratio) from chloroform-fumigated and not-fumigated soil samples [31]. For Cmic and Nmic, a correction factor of 0.45 and 0.54, respectively, has been then applied, as proposed by Vance et al. [31]. The C and N extracted from not-fumigated soil samples represented the labile pool of C and N (K2SO4-C and K2SO4-N).
The microbial respiration was determined placing the soil samples in closed glass jars, incubating them at 25 °C in the dark at field capacity and trapping the evolved CO2 by 0.5 N NaOH. The C-CO2 evolved after 1, 2, 4, 7, 10, 14, 17, 21, 28 days was determined by 0.05 N HCl titration of the excess of NaOH [32]. From data respiration, the following indicators have been calculated: basal respiration (BR), as the hourly respiration rate; CO2 cumulative (Ccum), as the cumulated CO2 evolved during the incubation period; metabolic quotient (qCO2), as the hourly CO2 evolved per unit of microbial biomass, and used to evaluate the effects of external disturbances [33]; the mineralization quotient (qM), as the ratio between the Ccum and the OC, indicating the efficiency of micro-flora in metabolizing soil OC [34]. The microbial quotient (qMIC) was also calculated, as the Cmic:OC ratio and it represents the percentage of microbial C with respect to the total organic C [34].
Additionally, a simplified biological fertility index (BFIs), was calculated. This type of index was used to discriminate soil biological fertility status under different agrosystems [35,36,37].

2.3.3. Macronutrient Content in Rhizosphere and Tomato

The concentration of P, K, S, Na, Ca, and Mg in rhizospheric soil was detected by ICP-OES (Ametek, Germany) after treating samples with suprapure HCl and HNO3 mixture (3:1 v/v) in microwave oven (Millestone, 1200) [29].
The macronutrients concentration was also analyzed in plant tissues and fruits by ICP-OES, after digesting of finely ground sample in microwave oven with suprapure HNO3 and 30% H2O2 mixture (4:1.5 v/v).

2.4. Data Treatment

Differences in the physiochemical and biochemical soil properties between sites have been verified by Kruskal–Wallis test using R software 4.0.2 (R Foundation for Statistical Computing, Vienna, Austria). We chose to use a non-parametric test because in our dataset some deviations from normality occurred, as well as the residues of the different groups had different ranges of variation (heteroskedastic), as verified by aov R implementation. Therefore, the key assumptions for ANOVA test calculation (i.e., normality and homoskedasticity of the subgroups) were not respected. The p-value greater than 0.05 indicates that there was no difference in ranks between the different groups and the significant codes.
In order to obtain the BFIs, the correlation coefficients were calculated among the soil parameters included in the biological fertility index (BFI) proposed by Benedetti and Mocali [38] to identify and remove the redundant variables. Furthermore, the BFIs calculation was based on the sum of the scores given by non-correlated parameters (Table S2 of the Supplementary Materials).

3. Results

3.1. Soil Physicochemical Properties

The investigated agricultural sites showed a different particle-size distribution; in particular, we found a silty loam texture in MEZ, clay/silty clay in MO and sand/sandy loam in RA (data not shown).
Table 1 and Table 2 report the values of soil physicochemical characteristics in the three different sites. Generally, no differences were detected in the physicochemical soil properties between top- and subsoil (0–20 and 20–40 cm, respectively), while significant differences were evident among the pedoclimatic areas.
The pH ranged from neutral (MEZ) to subalkaline (RA, MO), while the EC varied from moderate saline (MEZ) to no saline soils (RA and MO) (Table 1). MEZ and RA soils were poor in carbonates (<10 g kg–1), while a slight increase in MO soils was evident, where the CaCO3 content reached 22 and 21 g kg–1, respectively (Table 1). Both for top- and subsoils, the CEC values trend was MEZ > MO >RA (Table 1). Calcium was, in general, the most represented cation on the exchange complex and its concentrations showed the highest and the lowest values in MEZ and RA, respectively. The content of exchangeable Mg and K was similar in MEZ and MO and higher than in RA. Besides calcium, MEZ soils showed the highest concentrations of exchangeable Na, while the lowest one was detected in MO.
As expected, due to ancient land reclamation, peaty MEZ soils were enriched in OC (99.1 and 96.5 g kg–1, respectively, in top- and subsoil) (Table 2). In MO site, soil showed significantly higher OC content (14.4 g kg–1) than RA (7.5 g kg–1). Similarly, total N showed the highest concentrations in MEZ while the lowest ones in RA. The OC:TN ratio showed the highest values in MEZ, while no differences occurred between RA and MO (Table 2). The isotopic C signature showed high values in RA, and low ones in MEZ.
Due to the high variability, few differences occurred for total P and the POlsen contents. In fact, TP showed some differences only in topsoil where MO had a higher TP content than MEZ and RA. For POlsen, the topsoil showed higher values in MO than in MEZ, while the subsoil showed the lowest values in MEZ (Table 2).

3.2. Soil Biochemical Properties

Soil biochemical properties are shown in Table 2. Similar to physicochemical data, in general no significant differences in the microbial parameters between the soil depths were found, while several differences among the three pedoclimatic areas were observed. As expected, MEZ showed the highest values of C and N labile pools (K2SO4-C and K2SO4-N). Furthermore, higher concentrations of K2SO4-C were found in MO than in RA. The MEZ soils had the highest C-CO2 emissions, both as basal and cumulative soil respiration, while no differences occurred between MO and RA. The microbial C content trend was MO > MEZ = RA in topsoils, while the microbial N values were similar among sites, both in top- and subsoils.
The qMIC trend was MO > RA > MEZ in the topsoil, while in the subsoil MEZ showed the lowest qMIC values and no differences occurred between MO and RA (Figure 2). For both soil layers, the qCO2 showed the highest and lowest values in MEZ and MO, respectively (Figure 2). The qM of the topsoil showed the highest values in MO and the lowest ones in RA. In subsoil, instead, MEZ had a higher qM compared to MO and RA (Figure 2).
Checking the correlations among the values of indicators used for BFI calculation (Table 3), we found a strong positive correlation between OC amount and the CO2 evolved during the 28-day incubation time (r = 0.62 and 0.63 for BR and Ccum, respectively) as well as between BR and Ccum (r = 0.99).
Therefore, the BFIs took into account the following indicators: OC and microbial C amount, mineralization quotient and metabolic quotient. In the BFI formulation, the range of soil fertility classes is (i) BFIs = 4 stressed soils with very low fertility; (ii) 4 < BFIs ≤ 8 pre-stress soils; (iii) 8 < BFIs ≤ 12 soils with intermediate fertility; (iv) 12 < BFIs ≤ 16 good fertility soils; (v) 16 < BIFs ≤ 20 soils with very high fertility. Thus, the values of BFIs, calculated for investigated soils, in both 0–20 and 20–40 cm soil depth were MO = 15 corresponding to IV fertility class (good), MEZ = 11 in III fertility class (intermediate) and RA = 6 in II fertility class (pre-stress) (Table 4).

3.3. Nutrient Content into Rhizosphere and Tomato Organs

Figure 3 shows the nutrient content in rhizosphere and in different organs of tomato plants. With the exception of Mg, significant differences in nutrient content in the rhizospheric soil were detected among sites (Table S3 of the Supplementary Materials); in particular, MO soils had the highest content in Ca, K, Na and P while MEZ soils had the highest content in S. The lowest concentrations of Ca and Na were found in the rhizosphere of MEZ and RA, respectively. In the tomato organs, instead, very few differences were found and only S seemed to maintain the detected difference among rhizospheric soils, being the highest in the young leaves and stem of plants grown in MEZ soils (Figure 3 and Table S3 of the Supplementary Materials).

4. Discussion

4.1. Processing Tomato Cultivation Causes Higher Stress Conditions in Soils of MEZ and RA Study Sites Than in Those of MO Study Site

The highest OC content and the most negative δ13C value (ca. −27‰) in MEZ soils was in agreement with the presence of ancient reclaimed peats. Indeed, in peatlands, accumulation of OC occurs and the δ13C values of the accumulated plant debris are preserved due to anoxic conditions that limit the organic matter degradation in the peat and, therefore, no isotopic C fractionation occurs [39]. The highest OC content resulted in a higher CEC, which is an important soil feature in cropland ecosystems because it allows retention of the nutrients added through fertilization.
It is noteworthy to mention that, despite MEZ soils showing the highest concentrations of soluble C and N, which are the most easily available forms for microorganisms [40] promoting soil microbial biomass growth [41], the Cmic content was low. Moreover, MEZ soils had the lowest qMIC and the highest qCO2. While the highest labile organic matter concentration can be assigned to the highest OC content [42], the low Cmic concentrations and the high qCO2 should be attributed to the high OC:TN ratio of the soil organic matter. In fact, the high OC:TN ratio observed in MEZ soils, suggesting a low organic matter quality [43], likely prevented the microbial biomass growth [44]. However, in MEZ, we cannot exclude the possible depressing effect of salinity on microbial community [45]. The more stressful conditions occurring in MEZ study sites were, however, confirmed by the high value of qCO2. In fact, the qCO2 depicts the energetic efficiency of the microbial community [34]; therefore, a high qCO2 indicates a high energy demand for microbial biomass maintenance [46]. The lowest qMIC in MEZ soils would further suggest a scarce efficiency of microbes to convert the soil organic carbon into microbial biomass [47]. Additionally, in MEZ soils, the qMIC value was greatly lower than 2, which is considered a critical stress threshold for soils with neutral pH [34]. These findings supported the stressful condition for microbial biomass growth in MEZ soils.
In contrast to MEZ soils, RA and MO had different C use efficiency and the microbial respiration. Soils of RA study sites, having the coarser texture, showed the lowest OC and K2SO4-C contents which might explain the low values of microbial respiration and Cmic. However, RA soils had the highest δ13C values which would suggest how the small microbial population harboring in RA soils possess a high metabolism. This was also showed by the quite high qMIC, which resulted in an accelerated soil organic matter decomposition. MO soils are characterized by a scarce OC content and high incorporation of soil organic C in the microbial biomass, as highlighted by the generally higher values of qMIC and qM. Furthermore, MO soils showed the lowest qCO2 values suggesting relatively low stress levels [48], lower also than in RA.
The low soil organic matter quality of the MEZ soils and the imbalance between the current use and properties of these soils should contribute to reduction in the C use efficiency and growth of microbial population with respect to RA and MO soils. We cannot in fact exclude the effect of diverse fertilization strategies on the greater C use efficiency and the low microbial respiration occurring in MO and RA soils compared to MEZ ones. In fact, both in MO and RA, a higher amount of N fertilizers was added compared to MEZ which likely had a positive effect on soil microbial respiration [49,50]. Additionally, the contrasting soil texture between MO and RA soils might, at least partially, explain the lower microbial stress condition in MO than in RA—it being well known that finer soil texture promotes the microbial biomass growth [51] in contrast to coarser ones [52].
Overall, our findings illustrated that, in MEZ soils, despite the high concentration of organic C also in labile form, the microbial community does not immobilize that C into the microbial biomass, but it is lost through CO2 emission (i.e., high basal and cumulate respiration values) and/or leached. On the contrary, MO soils, showed a less stressed microbial community characterized by a more efficient use of C resources. RA soils, instead, although showing similar microbial features as the MO ones, cannot be considered a healthy system. Indeed, taking into account the qCO2:OC ratio, a qualitative index proposed by Dilly [53], both RA and MEZ soils had values greater than 400, which testifies an inefficient ecophysiological energy state of the system [53].
The better soil biological and biochemical quality of MO study sites is further confirmed by the BFI value which classifies it as a soil with a “good” biological fertility. Conversely, RA soils converged in a class characterized by “pre-stress” conditions, while MEZ soils were featured by an “intermediate” biological fertility. It is interesting to observe that the differentiation occurring for BFIs between MEZ and MO is missing for BFI. In particular, according to BFI, both MEZ and MO soils converged in “good” class, suggesting an overestimation of the index likely due to the presence of redundant and strongly correlated variables in its calculation. The proposed BFIs well describe the interaction between microbial communities and their activity even in soils where the high organic matter can mask some processes. Furthermore, the indications recorded by BFIs were coherent to the suggestions given by the qCO2:OC ratio proposed by Dilly [53].

4.2. Crop Quality

Despite the different microbial replies to the diverse soil types, very few differences of nutrient concentrations in tomato plant organs were detected. The general lack of differences of nutrient concentrations in plant organs can be attributed to the balanced use of chemical fertilizers applied by the farmers. Specifically, the amount of nutrients added to the soils was adjusted on the base of nutrient contents already present in soil with the aim to satisfy the plant nutritional requirements. However, noteworthy were the higher yield and brix grades of processing tomatoes grown on MO soils compared to MEZ and RA (Table S1 of the Supplementary Materials). Due to the role of soil nutrients on tomato yield [54,55], the highest yield of processing tomatoes grown on MO could be attributed to the highest concentrations of Ca, K and P in the rhizosphere of tomatoes located in MO. Instead, the lower yield of processing tomatoes cultivated in MEZ could be due to the higher salinity [56]. However, the differences of yield and brix grades might also be assigned to the diversity of the soil biological and biochemical fertility, which is a key component for crop production and quality [57,58]. Hence, it appears that the critical soil conditions detected by the BFIs relating to inefficient C use by soil microbial communities can lead to a worsening trend of crop production. Overall, in our opinion the data indicated a good suitability of MO area to the studied farming system, but an environmental risk related to soil degradation exists in the MEZ and RA sites.

5. Conclusions

Soil is an essential part of the CZ and is deeply influenced by anthropogenic actions, such as agricultural activities. Therefore, there is a need for measurable proxies that indicate the health status of agricultural soils. In this light, our study provided evidence about the usefulness of the biochemical properties to evaluate the stressful conditions of soils. Specifically, low values of qMIC and Cmic and high values of qCO2 indicate soils where the microbial community has a low fitness to incorporate carbon in their own structures and, therefore, a scarce growth of microbial biomass. In this sense, MEZ and RA sites showed less biologically active soils. Furthermore, our study highlighted how the BFIs, where the highly correlated parameters were removed, could be an effective indicator of soil stress conditions.
The better yield and brix grades of processing tomatoes observed in those soils with higher values of BFIs, qMIC and Cmic and lower values of qCO2 (MO sites) confirmed our second hypothesis about the influence of soil biological fertility on crop yield and quality.
Finally, the present paper shows that agroecosystems, i.e., important CZ components, have to be investigated with suitable biogeochemical indexes that help to understand the soil functioning, giving indication for a sustainable agricultural management. Indeed, a soil characterized by scarce biological fertility could indicate the occurrence of degradation processes.

Supplementary Materials

The following are available online at https://www.mdpi.com/2075-163X/11/2/219/s1, Table S1: Amount of nitrogen (N), phosphorus (P2O5) and potassium (K2O) applied by fertigation to soil, the mean yield and brix degree of tomatoes of processing tomato fields located in Modena (MO), Ferrara (MEZ) and Ravenna (RA) provinces, Emilia Romagna Region, Italy, Table S2: The scores of indicators for BFI calculation, Table S3: Results of Kruskal–Wallis test on the nutrients content in rhizospere (Rhizo) and tomato organs (Roots, Old and Young leaves, Stem, Fruit) in processing tomato fields located in Modena (MO), Ferrara (MEZ) and Ravenna (RA) provinces, Emilia Romagna Region, Italy.

Author Contributions

Conceptualization, L.V.A. and G.F.; methodology, L.V.A. and G.F.; formal analysis, C.F. and G.F.; investigation, C.N., C.F. and G.B.; writing—M.D.F. and G.F.; writing—review and editing, L.V.A., M.D.F., G.B. and G.F.; visualization, M.D.F.; supervision, L.V.A. and G.F.; project administration, L.V.A. and G.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank Carla Scotti for her useful contribution in the identification and involvement of farmers.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. FAO. ITPS Status of the World’s Soil Resources (SWSR)—Main report. In Food and Agriculture Organization of the United Nations and Intergovernmental Technical Panel on Soils; FAO: Rome, Italy, 2015; p. 608. [Google Scholar]
  2. Foley, J.A.; DeFries, R.; Asner, G.P.; Barford, C.; Bonan, G.; Carpenter, S.R.; Chapin, F.S.; Coe, M.T.; Daily, G.C.; Gibbs, H.K.; et al. Global consequences of land use. Science 2005, 309, 570–574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Tsiafouli, M.A.; Thébault, E.; Sgardelis, S.P.; de Ruiter, P.C.; van der Putten, W.H.; Birkhofer, K.; Hemerik, L.; de Vries, F.T.; Bardgett, R.D.; Brady, M.V.; et al. Intensive agriculture reduces soil biodiversity across Europe. Glob. Chang. Biol. 2015, 21, 973–985. [Google Scholar] [CrossRef] [PubMed]
  4. Bender, S.F.; Wagg, C.; van der Heijden, M.G.A. An underground revolution: Biodiversity and soil ecological engineering for agricultural sustainability. Trends Ecol. Evol. 2016, 31, 440–452. [Google Scholar] [CrossRef] [PubMed]
  5. Vittori Antisari, L.; Speranza, M.; Ferronato, C.; De Feudis, M.; Vianello, G.; Falsone, G. Assessment of water quality and soil salinity in the agricultural coastal plain Ravenna, North Italy. Minerals 2020, 10, 369. [Google Scholar] [CrossRef] [Green Version]
  6. Dilly, O.; Pompili, L.; Benedetti, A. Soil micro-biological indicators separated land use practices in contrast to abiotic soil properties at the 50 km scale under summer warm Mediterranean climate in northern Italy. Ecol. Indic. 2018, 84, 298–303. [Google Scholar] [CrossRef]
  7. Lori, M.; Symnaczik, S.; Mäder, P.; De Deyn, G.; Gattinger, A. Organic farming enhances soil microbial abundance and activity—A meta-analysis and meta-Regression. PLoS ONE 2017, 12, e0180442. [Google Scholar] [CrossRef] [PubMed]
  8. Richter, A.; Huallacháin, D.; Doyle, E.; Clipson, N.; Van Leeuwen, J.P.; Heuvelink, G.B.; Creamer, R.E. Linking diagnostic features to soil microbial biomass and respiration in agricultural grassland soil: A large-scale study in Ireland. Eur. J. Soil Sci. 2018, 69, 414–428. [Google Scholar] [CrossRef] [Green Version]
  9. Urra, J.; Alkorta, I.; Lanzén, A.; Mijangos, I.; Garbisu, C. The application of fresh and composted horse and chicken manure affects soil quality, microbial composition and antibiotic resistance. Appl. Soil Ecol. 2019, 135, 73–84. [Google Scholar] [CrossRef]
  10. Ghosh, A.; Singh, A.B.; Kumar, R.V.; Manna, M.C.; Bhattacharyya, R.; Rahman, M.M.; Sharma, P.; Rajput, P.S.; Misra, S. Soil enzymes and microbial elemental stoichiometry as bio-indicators of soil quality in diverse cropping systems and nutrient management practices of Indian Vertisols. Appl. Soil Ecol. 2020, 145, 103304. [Google Scholar] [CrossRef]
  11. Kabiri, V.; Raiesi, F.; Ghazavi, M.A. Tillage effects on soil microbial biomass, SOM mineralization and enzyme activity in a semi-arid Calcixerepts. Agric. Ecosyst. Environ. 2016, 232, 73–84. [Google Scholar] [CrossRef]
  12. Schloter, M.; Nannipieri, P.; Sørensen, S.J.; van Elsas, J.D. Microbial indicators for soil quality. Biol. Fertil. Soils 2018, 54, 1–10. [Google Scholar] [CrossRef] [Green Version]
  13. Trivedi, P.; Delgado-Baquerizo, M.; Anderson, I.C.; Singh, B.K. Response of soil properties and microbial communities to agriculture: Implications for primary productivity and soil health indicators. Front. Plant Sci. 2016, 7, 990. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Guillaume, T.; Maranguit, D.; Murtilaksono, K.; Kuzyakov, Y. Sensitivity and resistance of soil fertility indicators to land-use changes: New concept and examples from conversion of Indonesian rainforest to plantations. Ecol. Indic. 2016, 67, 49–57. [Google Scholar] [CrossRef]
  15. Arif, M.S.; Riaz, M.; Shahzad, S.M.; Yasmeen, T.; Ashraf, M.; Siddique, M.; Mubarik, M.S.; Bragazza, L.; Buttler, A. Fresh and composted industrial sludge restore soil functions in surface soil of degraded agricultural land. Sci. Total Environ. 2018, 619, 517–527. [Google Scholar] [CrossRef] [PubMed]
  16. De Santiago, A.; Recena, R.; Perea-Torres, F.; Moreno, M.T.; Carmona, E.; Delgado, A. Relationship of soil fertility to biochemical properties under agricultural practices aimed at controlling land degradation. L. Degrad. Dev. 2019, 30, 1121–1129. [Google Scholar] [CrossRef]
  17. Widdig, M.; Schleuss, P.M.; Biederman, L.A.; Borer, E.T.; Crawley, M.J.; Kirkman, K.P.; Seabloom, E.W.; Wragg, P.D.; Spohn, M. Microbial carbon use efficiency in grassland soils subjected to nitrogen and phosphorus additions. Soil Biol. Biochem. 2020, 146, 107815. [Google Scholar] [CrossRef]
  18. Nunes, J.S.; Araujo, A.S.F.; Nunes, L.A.P.L.; Lima, L.M.; Carneiro, R.F.V.; Salviano, A.A.C.; Tsai, S.M. Impact of land degradation on soil microbial biomass and activity in Northeast Brazil. Pedosphere 2012, 22, 88–95. [Google Scholar] [CrossRef]
  19. Kumar, A.; Maurya, B.R.; Raghuwanshi, R. Isolation and characterization of PGPR and their effect on growth, yield and nutrient content in wheat (Triticum aestivum L.). Biocatal. Agric. Biotechnol. 2014, 3, 121–128. [Google Scholar] [CrossRef]
  20. Plenchette, C.; Clermont-Dauphin, C.; Meynard, J.M.; Fortin, J.A. Managing arbuscular mycorrhizal fungi in cropping systems. Can. J. Plant Sci. 2005, 85, 31–40. [Google Scholar] [CrossRef]
  21. Singh, J.S.; Kumar, A.; Rai, A.N.; Singh, D.P. Cyanobacteria: A precious bio-resource in agriculture, ecosystem, and environmental sustainability. Front. Microbiol. 2016, 7, 529. [Google Scholar] [CrossRef] [Green Version]
  22. Vimal, S.R.; Singh, J.S.; Arora, N.K.; Singh, S. Soil-plant-microbe interactions in stressed agriculture management: A Review. Pedosphere 2017, 27, 177–192. [Google Scholar] [CrossRef]
  23. Regione Emilia Romagna. Cartografia dei Suoli della Regione Emilia Romagna 1:50,000 (1994 updates 2000). Available online: https://geo.regione.emilia-romagna.it/cartografia_sgss/user/viewer.jsp?service=pedologia&bookmark=1%22 (accessed on 10 October 2019).
  24. ARPAE, Regione Emilia Romagna. Tabelle Climatologiche. Servizio Idro-Meteo-Clima. Available online: https://www.arpae.it/sim/?osservazioni_e_dati/climatologia (accessed on 10 October 2019).
  25. Regione Emilia Romagna. Disciplinari di Produzione Integrata e Deroghe. Available online: https://agricoltura.regione.emilia-romagna.it/fitosanitario/temi/difesa-sostenibile-delle-produzioni/disciplinari-di-produzione-integrata/disciplinari-di-produzione-integrata-1 (accessed on 3 October 2020).
  26. Orsini, L.; Rémy, J. Utilisation du chlorure de cobaltihexamine pour la détermination simultanée de la capacité d’échange et des bases échangeables des sols. Bulletin de l’Association Française d’Etude du Sol. Bull. Association Française Etude Sol 1976, 4, 269–279. [Google Scholar]
  27. Loeppert, R.H.; Suarez, D.L. Carbonate and Gypsun; USDA-ARS/UNL Faculty: Lincoln, NE, USA, 1996. [Google Scholar]
  28. Gee, G.W.; Bauder, J.W. Methods of Soil Analysis: Part 1—Physical and Mineralogical Methods; SSSA Book Series; Soil Science Society of America, American Society of Agronomy: Madison, WI, USA, 1986; ISBN 978-0-89118-864-3. [Google Scholar]
  29. Vittori Antisari, L.; Carbone, S.; Gatti, A.; Vianello, G.; Nannipieri, P. Uptake and translocation of metals and nutrients in tomato grown in soil polluted with metal oxide (CeO2, Fe3O4, SnO2, TiO2) or metallic (Ag, Co, Ni) engineered nanoparticles. Environ. Sci. Pollut. Res. 2015, 22, 1841–1853. [Google Scholar] [CrossRef]
  30. Olsen, S.R.; Cole, C.V.; Watandbe, F.; Dean, L. Estimation of available phosphorus in soil by extraction with sodium bicarbonate. J. Chem. Inf. Model. 1954, 53, 1689–1699. [Google Scholar]
  31. Vance, E.D.; Brookes, P.C.; Jenkinson, D.S. An extraction method for measuring soil microbial biomass C. Soil Biol. Biochem. 1987, 19, 703–707. [Google Scholar] [CrossRef]
  32. De Feudis, M.; Falsone, G.; Vittori Antisari, L. Mid-term (30 years) changes of soil properties under chestnut stands due to organic residues management: An integrated study. Catena 2020, 198, 105021. [Google Scholar] [CrossRef]
  33. Anderson, T.H.; Domsch, K.H. The metabolic quotient for CO2 (qCO2) as a specific activity parameter to assess the effects of environmental conditions, such as ph, on the microbial biomass of forest soils. Soil Biol. Biochem. 1993, 25, 393–395. [Google Scholar] [CrossRef]
  34. Anderson, T.H. Microbial eco-physiological indicators to asses soil quality. Agric. Ecosyst. Environ. 2003, 98, 285–293. [Google Scholar] [CrossRef]
  35. Pompili, L.; Mellina, A.S.; Benedetti, A.; Bloem, J. Microbial indicators in three agricultural soils with different management. Fresenius Environ. Bull. 2008, 17, 1128–1136. [Google Scholar]
  36. Renzi, G.; Canfora, L.; Salvati, L.; Benedetti, A. Validation of the soil Biological Fertility Index (BFI) using a multidimensional statistical approach: A country-scale exercise. Catena 2017, 149, 294–299. [Google Scholar] [CrossRef]
  37. Francaviglia, R.; Renzi, G.; Ledda, L.; Benedetti, A. Organic carbon pools and soil biological fertility are affected by land use intensity in Mediterranean ecosystems of Sardinia, Italy. Sci. Total Environ. 2017, 599–600, 789–796. [Google Scholar] [CrossRef] [PubMed]
  38. Benedetti, A.; Mocali, S. Analisi a livello di suolo. In Indicatori di Biodiversità per la Sostenibilità in Agricoltura. Linee Guida, Strumenti e Metodi per la Valutazione della Qualità degli Agroecosistemi; Caporali, F., Vazzana, C., Benedetti, A., Calabrese, J., Mancinelli, R., Lazzerini, G., Mocali, S., Campiglia, E., et al., Eds.; ISPRA: Rome, Italy, 2008; pp. 159–208. [Google Scholar]
  39. Krüger, J.P.; Leifeld, J.; Alewell, C. Degradation changes stable carbon isotope depth profiles in palsa peatlands. Biogeosciences 2014, 11, 3369–3380. [Google Scholar] [CrossRef] [Green Version]
  40. Blagodatsky, S.; Blagodatskaya, E.; Yuyukina, T.; Kuzyakov, Y. Model of apparent and real priming effects: Linking microbial activity with soil organic matter decomposition. Soil Biol. Biochem. 2010, 42, 1275–1283. [Google Scholar] [CrossRef]
  41. Cookson, W.R.; Abaye, D.A.; Marschner, P.; Murphy, D.V.; Stockdale, E.A.; Goulding, K.W.T. The contribution of soil organic matter fractions to carbon and nitrogen mineralization and microbial community size and structure. Soil Biol. Biochem. 2005, 37, 1726–1737. [Google Scholar] [CrossRef]
  42. Rizinjirabake, F.; Tenenbaum, D.E.; Pilesjö, P. Sources of soil dissolved organic carbon in a mixed agricultural and forested watershed in Rwanda. Catena 2019, 181, 104085. [Google Scholar] [CrossRef]
  43. Breulmann, M.; Schulz, E.; Weißhuhn, K.; Buscot, F. Impact of the plant community composition on labile soil organic carbon, soil microbial activity and community structure in semi-natural grassland ecosystems of different productivity. Plant Soil 2012, 352, 253–265. [Google Scholar] [CrossRef]
  44. Eiland, F.; Klamer, M.; Lind, A.M.; Leth, M.; Bååth, E. Influence of initial C/N ratio on chemical and microbial composition during long term composting of straw. Microb. Ecol. 2001, 41, 272–280. [Google Scholar] [CrossRef]
  45. Rietz, D.N.; Haynes, R.J. Effects of irrigation-induced salinity and sodicity on soil microbial activity. Soil Biol. Biochem. 2003, 35, 845–854. [Google Scholar] [CrossRef]
  46. Okolo, C.C.; Dippold, M.A.; Gebresamuel, G.; Zenebe, A.; Haile, M.; Bore, E. Assessing the sustainability of land use management of northern Ethiopian drylands by various indicators for soil health. Ecol. Indic. 2020, 112, 106092. [Google Scholar] [CrossRef]
  47. Anderson, T.H.; Domsch, K.H. Ratios of microbial biomass carbon to total organic carbon in arable soils. Soil Biol. Biochem. 1989, 21, 471–479. [Google Scholar] [CrossRef]
  48. Pabst, H.; Gerschlauer, F.; Kiese, R.; Kuzyakov, Y. Land use and precipitation affect organic and microbial carbon stocks and the specific metabolic quotient in soils of eleven ecosystems of Mt. Kilimanjaro, Tanzania. L. Degrad. Dev. 2016, 27, 592–602. [Google Scholar] [CrossRef]
  49. Spohn, M.; Pötsch, E.M.; Eichorst, S.A.; Woebken, D.; Wanek, W.; Richter, A. Soil microbial carbon use efficiency and biomass turnover in a long-term fertilization experiment in a temperate grassland. Soil Biol. Biochem. 2016, 97, 168–175. [Google Scholar] [CrossRef]
  50. Ward, D.; Kirkman, K.; Hagenah, N.; Tsvuura, Z. Soil respiration declines with increasing nitrogen fertilization and is not related to productivity in long-term grassland experiments. Soil Biol. Biochem. 2017, 115, 415–422. [Google Scholar] [CrossRef]
  51. Bauhus, J.; Paré, D.; Côté, L. Effects of tree species, stand age and soil type on soil microbial biomass and its activity in a southern boreal forest. Soil Biol. Biochem. 1998, 30, 1077–1089. [Google Scholar] [CrossRef]
  52. Blanco-Moure, N.; Gracia, R.; Bielsa, A.C.; López, M.V. Soil organic matter fractions as affected by tillage and soil texture under semiarid Mediterranean conditions. Soil Tillage Res. 2016, 155, 381–389. [Google Scholar] [CrossRef] [Green Version]
  53. Dilly, O. Microbial energetics in soils. In Microorganisms in Soils: Roles in Genesis and Functions; Springer: Berlin/Heidelberg, Germany, 2005; pp. 123–138. [Google Scholar]
  54. Heeb, A.; Lundegårdh, B.; Savage, G.; Ericsson, T. Impact of organic and inorganic fertilizers on yield, taste, and nutritional quality of tomatoes. J. Plant Nutr. Soil Sci. 2006, 169, 535–541. [Google Scholar] [CrossRef]
  55. Fandi, M.; Muhtaseb, J.; Hussein, M. Effect of N, P, K concentrations on yield and fruit quality of (Solanum lycopersicum L.) in tuff culture. J. Cent. Eur. Agric. 2010, 11, 179–184. [Google Scholar]
  56. Wold, A.B.; Rosenfeld, H.J.; Baugerød, H.; Blomhoff, R. The effect of fertilization on antioxidant activity and chemical composition of tomato cultivars (Lycopersicon esculentum Mill.). Eur. J. Hortic. Sci. 2004, 69, 167–174. [Google Scholar]
  57. Jannoura, R.; Joergensen, R.G.; Bruns, C. Organic fertilizer effects on growth, crop yield, and soil microbial biomass indices in sole and intercropped peas and oats under organic farming conditions. Eur. J. Agron. 2014, 52, 259–270. [Google Scholar] [CrossRef]
  58. Tejada, M.; Gonzalez, J.L.; García-Martínez, A.M.; Parrado, J. Effects of different green manures on soil biological properties and maize yield. Bioresour. Technol. 2008, 99, 1758–1767. [Google Scholar] [CrossRef]
Figure 1. Map of Italy with magnification of the Emilia Romagna region and indication of the study sites.
Figure 1. Map of Italy with magnification of the Emilia Romagna region and indication of the study sites.
Minerals 11 00219 g001
Figure 2. Microbial (qMIC), mineralization (qM) and metabolic (qCO2) quotients and in 0–20 and 20–40 cm soil layers of processing tomato fields located in Modena (MO), Ferrara (MEZ) and Ravenna (RA) provinces, Emilia Romagna Region, Italy. Error bars represent standard error of the mean. Different letters mean significant differences (Kruskal–Wallis test, p < 0.05) among sites within 0–20 and 20–40 cm soil layers.
Figure 2. Microbial (qMIC), mineralization (qM) and metabolic (qCO2) quotients and in 0–20 and 20–40 cm soil layers of processing tomato fields located in Modena (MO), Ferrara (MEZ) and Ravenna (RA) provinces, Emilia Romagna Region, Italy. Error bars represent standard error of the mean. Different letters mean significant differences (Kruskal–Wallis test, p < 0.05) among sites within 0–20 and 20–40 cm soil layers.
Minerals 11 00219 g002
Figure 3. Ca, K, Mg, Na, P and S concentrations in rhizosphere soil, root, old and young leaves, stem and fruit of processing tomato plants cultivated in Modena (MO), Ferrara (MEZ) and Ravenna (RA) provinces, Emilia Romagna Region, Italy. Error bars represent standard error of the mean.
Figure 3. Ca, K, Mg, Na, P and S concentrations in rhizosphere soil, root, old and young leaves, stem and fruit of processing tomato plants cultivated in Modena (MO), Ferrara (MEZ) and Ravenna (RA) provinces, Emilia Romagna Region, Italy. Error bars represent standard error of the mean.
Minerals 11 00219 g003
Table 1. Mean and standard error (SE) values of soil pH, electrical conductivity (EC), CaCO3 content, cation exchange capacity (CEC), contents of exchangeable Ca, Mg, K and Na in 0–20 and 20–40 cm soil layers of processing tomato fields located in Modena (MO), Ferrara (MEZ) and Ravenna (RA) provinces, Emilia Romagna Region, Italy. Different letters mean significant differences (Kruskal–Wallis test, p < 0.05) among sites within 0–20 and 20–40 cm soil layers.
Table 1. Mean and standard error (SE) values of soil pH, electrical conductivity (EC), CaCO3 content, cation exchange capacity (CEC), contents of exchangeable Ca, Mg, K and Na in 0–20 and 20–40 cm soil layers of processing tomato fields located in Modena (MO), Ferrara (MEZ) and Ravenna (RA) provinces, Emilia Romagna Region, Italy. Different letters mean significant differences (Kruskal–Wallis test, p < 0.05) among sites within 0–20 and 20–40 cm soil layers.
Soil Depth MOMEZRA
MeanSEMeanSEMeanSE
0–20 cmpH 8.3 a0.17.1 c0.27.7 b0.1
ECμS cm−1690 b1581465 a221975 b72
CaCO3 g kg−121.1 a11.4 c0.27.1 b0.8
CECcmol(+) kg−121.5 b1.940.2 a5.35.5 c2.6
Cacmol(+) kg−117.1 b0.532.9 a0.93.5 c1.4
Mgcmol(+) kg−12.3 a0.13.5 a0.10.2 b<0.1
Kcmol(+) kg−11.4 a0.10.4 a0.10.1 b<0.1
Nacmol(+) kg−10.3 c0.052.5 a0.31.4 b0.2
20–40 cmpH 8.4 a0.16.8 c0.37.8 b0.1
ECμS cm−1847 b2162774 a598931 b99
CaCO3 g kg−122.1 a0.71.5 c0.27.7 b0.9
CECcmol(+) kg−120.8 b2.139.8 a4.78.5 c3.5
Cacmol(+) kg−117.4 b0.532.4 a2.93.1 c1.4
Mgcmol(+) kg−12.4 a0.11.6 a0.70.2 b0.1
Kcmol(+) kg−10.4 a0.10.5 a0.10.1 b<0.1
Nacmol(+) kg−10.3 c0.14.5 a0.81.2 b0.2
Table 2. Mean and standard error (SE) values of soil organic C (OC) and total N contents, C:N ratio, δ 13C, total and available P contents (TP and POlsen, respectively), labile organic C and N contents (K2SO4-C and K2SO4-N, respectively), basal and cumulative respiration (BR and Ccum, respectively), microbial biomass C and N contents (Cmic and Nmic, respectively) in 0–20 and 20–40 cm soil layers of processing tomato fields located in Modena (MO), Ferrara (MEZ) and Ravenna (RA) provinces, Emilia Romagna Region, Italy. Different letters mean significant differences (Kruskal–Wallis test, p < 0.05) among sites within 0–20 and 20–40 cm soil layers.
Table 2. Mean and standard error (SE) values of soil organic C (OC) and total N contents, C:N ratio, δ 13C, total and available P contents (TP and POlsen, respectively), labile organic C and N contents (K2SO4-C and K2SO4-N, respectively), basal and cumulative respiration (BR and Ccum, respectively), microbial biomass C and N contents (Cmic and Nmic, respectively) in 0–20 and 20–40 cm soil layers of processing tomato fields located in Modena (MO), Ferrara (MEZ) and Ravenna (RA) provinces, Emilia Romagna Region, Italy. Different letters mean significant differences (Kruskal–Wallis test, p < 0.05) among sites within 0–20 and 20–40 cm soil layers.
Soil Depth MOMEZRA
MeanSEMeanSEMeanSE
0–20 cmOC g kg−114.4 b0.599.1 a10.27.5 c0.6
TN g kg−12.1 b0.16.1 a0.61.1 c0.1
OC:TN 6.8 b0.516.2 a0.26.8 b0.1
δ 13C −26.2 b0.1−26.9 c0.1−25.7 a0.2
Ptotg kg−11.2 a0.20.8 b0.10.7 b0.1
P Olsenmg kg−14.6 a1.62.7 b0.22.9 ab0.9
K2SO4-Cmg kg−192 b3378 a3070 c4
K2SO4-Nmg kg−128 b567 a718 b1
BRmg C-CO2 kg−115 b2120 a513 b2
Ccummg C-CO2 kg−1688 b572567 a125125 b49
Cmicmg kg−1438 a18143 b35121 b19
Nmicmg kg−1212255153
20–40 cmOC g kg−114.1 b0.795.6 a15.87.9 c0.6
NT g kg−12.1 b0.16.5 a0.61.1 c0.1
OC:TN 6.7 b0.514.7 a0.27.2 b0.1
δ 13C −26.5 ab0.6−26.8 b0.3−25.7 a0.2
Ptotg kg−11.10.10.800.70.1
P Olsenmg kg−13.8 a1.31.5 b0.34.3 a0.1
K2SO4-Cmg kg−196 b5405 a3580 c6
K2SO4-Nmg kg−123 b260 a531 b6
BRmg C-CO2 kg−112 b358 a711 b3
Ccummg C-CO2 kg−1443 b752843 a188143 b88
Cmicmg kg−1208171463113220
Nmicmg kg−1172256153
Table 3. Correlation matrix among the indicators used for biological fertility index (BFI) calculation in the three investigated sites. In bold, significant (p < 0.05) correlations are reported.
Table 3. Correlation matrix among the indicators used for biological fertility index (BFI) calculation in the three investigated sites. In bold, significant (p < 0.05) correlations are reported.
OCBRCumCmicqCO2qM
OC1
BR0.621
Ccum0.630.991
Cmic0.190.060.061
qCO2−0.29−0.09−0.090.631
qM0.69−0.32−0.32−0.320.681
OC: organic C; BR: basal respiration; Ccum: cumulative respiration; Cmic: microbial C; qCO2: metabolic quotient; qM: mineralization quotient.
Table 4. The simplified biological fertility index (BFIs) scores and corresponding fertility classes calculated in 0–20 and 20–40 cm soil layers of processing tomato fields located in Modena (MO), Ferrara (MEZ) and Ravenna (RA) provinces, Emilia Romagna Region, Italy.
Table 4. The simplified biological fertility index (BFIs) scores and corresponding fertility classes calculated in 0–20 and 20–40 cm soil layers of processing tomato fields located in Modena (MO), Ferrara (MEZ) and Ravenna (RA) provinces, Emilia Romagna Region, Italy.
SchemeSiteBFIs ScoreBFIs Class
MO15IV (good)
0–20 cmRA6II (pre-stress)
MEZ11III (intermediate)
MO15IV (good)
20–40 cmRA6II (pre-stress)
MEZ11III (intermediate)
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Vittori Antisari, L.; Ferronato, C.; De Feudis, M.; Natali, C.; Bianchini, G.; Falsone, G. Soil Biochemical Indicators and Biological Fertility in Agricultural Soils: A Case Study from Northern Italy. Minerals 2021, 11, 219. https://doi.org/10.3390/min11020219

AMA Style

Vittori Antisari L, Ferronato C, De Feudis M, Natali C, Bianchini G, Falsone G. Soil Biochemical Indicators and Biological Fertility in Agricultural Soils: A Case Study from Northern Italy. Minerals. 2021; 11(2):219. https://doi.org/10.3390/min11020219

Chicago/Turabian Style

Vittori Antisari, Livia, Chiara Ferronato, Mauro De Feudis, Claudio Natali, Gianluca Bianchini, and Gloria Falsone. 2021. "Soil Biochemical Indicators and Biological Fertility in Agricultural Soils: A Case Study from Northern Italy" Minerals 11, no. 2: 219. https://doi.org/10.3390/min11020219

APA Style

Vittori Antisari, L., Ferronato, C., De Feudis, M., Natali, C., Bianchini, G., & Falsone, G. (2021). Soil Biochemical Indicators and Biological Fertility in Agricultural Soils: A Case Study from Northern Italy. Minerals, 11(2), 219. https://doi.org/10.3390/min11020219

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop