Next Article in Journal
Recovery of 15N Labeled Nitrogen Fertilizer by Fertigated and Drip Irrigated Greenhouse Vegetable Crops
Next Article in Special Issue
Management Intensification of Hay Meadows and Fruit Orchards Alters Soil Macro- Invertebrate Communities Differently
Previous Article in Journal
Seed Morphology in the Vitaceae Based on Geometric Models
Previous Article in Special Issue
Do Long-Term Continuous Cropping and Pesticides Affect Earthworm Communities?
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Response of the Arthropod Community to Soil Characteristics and Management in the Franciacorta Viticultural Area (Lombardy, Italy)

1
Department of Molecular and Translational Medicine, University of Brescia, 25123 Brescia, Italy
2
Sata Studio Agronomico, 25121 Brescia, Italy
3
Department of Agricultural and Environmental Sciences—Production, Landscape, Agroenergy, University of Study of Milan, 20133 Milan, Italy
*
Author to whom correspondence should be addressed.
Submission received: 12 April 2020 / Revised: 13 May 2020 / Accepted: 18 May 2020 / Published: 21 May 2020

Abstract

:
Soil represents an important pool of biodiversity, hosting about a quarter of the living species on our planet. This soil richness has led to increasing interest in the structural and functional characteristics of its biodiversity. Studies of arthropod responses, in terms of abundance and taxon richness, have increased in relation to their ecological value as bioindicators of environmental change. This research was carried out over the 2014–2018 period with the aim to better understand arthropod taxa responses in vineyard soils in Franciacorta (Lombardy, Italy). To determine the biological composition in terms of arthropod taxa presence, one hundred soil samples were analysed. Environmental characteristics, such as chemical composition, soil moisture and temperature, and soil management were characterized for each soil sample. A total of 19 taxa were identified; the NMDS model analysis and the cluster analysis divided them into five groups according to their co-occurrence patterns. Each group was related to certain abiotic conditions; of these, soil moisture, temperature and organic matter were shown to be significant. A decision tree analysis showed that a longer period since conversion from conventional to organic farming lead to a higher arthropod biodiversity defined as a higher number of taxa.

1. Introduction

Soil has recently been described as the most complex and diverse ecosystem in the world [1], and it represents an important pool of biodiversity. It is indeed one of the richest habitats of terrestrial ecosystems in terms of species diversity [2,3]. The European Commission [4] estimates that about a quarter of living species on our planet are found in the soil, and the importance of this biodiversity has already been described in relation to the functional roles that the soil biota plays in regulating ecosystem processes [5]. However, despite the increasing number of studies on soil biodiversity, many structural and functional aspects of this biodiversity remain largely unexplored [5,6]. In this context, the investigation of the relationship between soil arthropod communities, in terms of abundance and taxa richness, and environmental conditions played an important role. The sensitivity of soil arthropods to environmental conditions [7,8], soil properties [9] and soil management practices [10] allows them to be considered as bioindicators of environmental change [11,12].
The relative importance of the factors influencing soil arthropod diversity and abundance in agroecosystems is still far from being understood. Indeed, the influence of abiotic and biotic variables and their interactions [13,14] varies according to the climate, type of soil and agricultural practices. The influence of meteorological variables (i.e., precipitation and air temperature), soil moisture and temperatures on soil arthropods has already been evaluated in different habitats [15,16,17,18,19,20,21,22,23]. Soil moisture and soil temperature have emerged as important factors that determine arthropod distribution [16,19,21,24,25], but the response of soil arthropods to soil water availability and temperature has been shown to vary between taxa [26,27,28,29]. In general, the positive effect of soil moisture on the abundance of soil arthropod communities has been emphasised [16,30] and an optimum temperature range of between 5 °C and 10 °C was identified for species active in winter, and between 10 °C and 18 °C for those active in summer [31] (p. 6). Soil chemical and physical characteristics have been identified as important drivers in soil arthropod distribution and abundance [9,11,13,20,32,33,34,35,36,37,38]. Soil texture [39], soil organic matter content [40,41], pH [11,37,38] and heavy metal concentration [13,20] have been shown to have a major influence on soil biota. Of these, soil pH and soil organic matter represent the most significant drivers in relation to the influence of soil pH variation on soil arthropods presence [11,37] and to arthropods role in soil organic matter degradation [34,35]. Soil arthropods contribute, in fact, to nutrient cycling as secondary decomposers, conditioning litter through comminution and passage through the gut, for further breakdown by the microflora [34] and stimulating microbial mineralisation of nutrients through grazing activity [35].
The influence of management on soil arthropods has been investigated in different agricultural contexts [10,42,43,44,45,46]. In particular, different studies have focused on the effect of organic viticulture on soil arthropod communities; different authors [15,47,48] show the general positive effect of organic management on soil arthropod abundance and distribution. However, the results have varied for each taxon investigated [49,50]. Furthermore, only few studies have been carried out evaluating the role of time of organic practice application on soil arthropod biodiversity [51]. Further research is therefore needed to assess the medium and long-term effects of organic agriculture on soil biodiversity [52].
In this paper, we report on the results of a 5-year investigation into the responses of the arthropod community to soil characteristics and vineyard management in the Franciacorta viticultural area (Lombardy, Italy). The diversity and co-occurrence patterns of different taxa were analysed in relation to abiotic factors, such as soil temperature, soil moisture and soil chemical properties. Moreover, the influence of vineyard management (conventional vs. organic) and the time of conversion from conventional to organic on arthropod biodiversity was investigated.

2. Materials and Methods

2.1. Study Sites

This study was carried out in a major Italian winemaking area. Franciacorta is the most famous Italian wine region for the production of sparkling wine using the champenoise method and is located in the Lombardy Region (Figure 1). The zone covers a total area of 2615 ha (as of 2018) and hosts 117 wineries (as of 2019). This research collected a total of 100 soil samples from 100 different vineyards over the period 2014–2018. Eighty-five per cent of samples were collected in spring (May or June) and 15% of samples in autumn (September, October or November). All the samples were characterised by presence of arthropods and the chemical characteristics of soil.
In Figure 1, the location of each sampling site is shown.
Vineyard management systems were classified in two main groups: conventionally managed vineyards without any specific environmental certification (conventional) and organic vineyards managed in compliance with the European Regulation on organic farming (reg EC n. 2018/848 and subsequent amendments and additions) (organic). For organic vineyards, we refer to the presence of this certification that implies compliance with the provisions of the law. In addition to this, we have verified a minimal set of conditions occurred in each farm monitored in organic farming. These actions refer to: no use of synthetic chemicals for plant protection and for fertilizing the vineyard; the integration of organic matter into the soil through the supply of organic matrices; the total absence of use of herbicides and the management of the sub-row through mechanical intervention; the preservation of the herbaceous covering on the ground; the minimum tillage adoption. Organic vineyards were then further divided into three subgroups, on the basis of how long ago they had been converted from conventional to organic farming: 3 years or less (organic ≤ 3), between 4 and 9 years (4 ≤ organic ≤ 9), and 10 years or more (organic ≥ 10).

2.2. Environmental Variables

Soil moisture (SM) and soil temperature (ST) data for the Franciacorta area from 2014 to 2018 were obtained from the National Centers for Environmental Predictions [53]. These data were then re-analysed using the Weather Research and Forecasting (WRF) simulations [54]. The WRF model (version 4.02) was applied to a high spatial resolution grid (each cell of the grid representing a 2 × 2 km area) to generate hourly data. In particular, the Noah scheme [55] has been used as land surface model (LSM) scheme (i.e., Noah, Noah-MP, and CLM4) to assess detailed multi-layer soil moisture and soil temperature. We focused on a depth of 0–15 cm below land surface. Each vineyard was associated with the nearest grid node to allow extraction of the specific soil temperature and moisture values.
To assess the influence of environmental variables on the presence of soil arthropods, SM and ST were evaluated for each vineyard in a 30-days reference period prior to the sampling date (Table A1). Two thermal thresholds ( ST low and ST up ) were considered to define two intervals of temperature. These intervals characterise organisms that prefer lower temperature features (taxa occurring more frequently in the interval [ ST low , ST up ]) or higher temperature features (taxa occurring more frequently when soil temperature is higher than ST up ). The lower threshold ( ST low ) was set at 10 °C and the upper threshold ( ST up ) was set at 20 °C [31] (p. 6). Soil temperatures lower than 10 °C were not included since they can be considered to be below the lower development threshold for most taxa. Considering ST low and ST up , two variables related to soil temperature were calculated for the reference period (720 h):
TL = 1 24 i = 1 720 ( ST i ST low ) ,   ST low ST i < ST up
TH = 1 24 i = 1 720 ( ST i ST up ) ,   ST i ST up  
where ST i is hourly soil temperature. TL is the daily cumulative soil temperature degrees exceeding 10 °C when ST i is between 10 °C and 20 °C; TH is the daily cumulative soil temperature degrees exceeding 20 °C when ST i is greater than 20 °C.
A soil moisture threshold ( SM ) was defined to discriminate between organisms that prefer drier conditions, i.e., taxa occurring more frequently when soil moisture ranges in the interval [ 0 , SM ], or wetter conditions, i.e., taxa that more frequently occur when soil moisture is in the range [ SM ,1]. SM was set equal 0.35 (corresponding to 35%), which represents a reference value that can be associated, in different ecological contexts, to a status that satisfies the requirements in terms of humidity of soil arthropods [21,22]. Two variables associated with soil moisture were calculated for the reference period (720 h):
MD = 1 24 i = 1 720 | SM i SM | , 0 SM i SM
MH = 1 24 i = 1 720 ( SM i SM ) ,   SM i > SM  
where SM i is the hourly soil moisture. MD is the daily sum of absolute deviations in soil moisture values from the threshold value when SM i is lower than 0.35; MH is the daily cumulative soil moisture exceeding 0.35, when SM i is higher than 0.35.

2.3. Chemical Characterisation of Soils

Chemical analysis of soils was performed according to the Italian regulation (DM 13 September 1999). Soil samples were taken at a depth of 0–15 cm and mixed homogeneously. Leaf litter was excluded, as it is not part of the soil itself. The collected soil samples were air-dried, homogenized and passed through a 2 mm sieve for chemical analysis.
Characterisation of the soil chemistry involved measuring soil texture (TXT), pH, active limestone (expressed in g CaCO3/kg of soil) (AL), organic matter content (expressed in g/kg of soil) (SOM), available phosphorus (mg P2O5/kg of soil) (P), available potassium (mg K2O/kg of soil) (K), available magnesium (mg MgO/kg of soil) (Mg) and copper content (mg/kg) (Cu). Soil texture was classified following the USDA soil texture triangle classification [56] (p. 125).

2.4. Soil Arthropods Identification

A cubic sample of soil (with a dimension of about 30 cm3) was collected at the same depth described for chemical soil analysis, at each vineyard. Arthropods were extracted by placing the soil sample in a Berlese–Tüllgren funnel under a 60 W incandescence bulb, leading soil arthropods to migrate towards the damp part of the soil sample (away from the light). The soil arthropods fell through the cavity, into a preserving solution (2/3 alcohol and 1/3 glycerol). Determination of biological forms was carried out according to the QBS-ar (Soil Biological Quality-arthropod) method as proposed by [57], and the definition of the taxonomic entities and the biological stages is in agreement with the one reported in the same paper.

2.5. Data Analysis

2.5.1. Taxa Co-Occurrence Patterns

To measure soil arthropods biodiversity a taxa co-occurrence approach was used. For each vineyard, a taxa presence profile was defined, i.e., a vector indicating the presence or absence of the taxa in each vineyard. The presence profile did not consider population abundance. Based on the presence profiles, vineyards and taxa were described in a J dimensional space (J is the number of taxa considered), allowing taxa to be ordered by their vineyard presence profiles. Two taxa are close to each other if they share a similar pattern of co-occurrence in the vineyards, they are far from each other if one is present in the vineyards where the other is absent and vice versa.
To allow easy visualisation and interpretation of dissimilarity in soil biodiversity and taxa co-occurrence, it is useful to represent these profiles in a two-dimensional space, called an ordination plane. Non-metric multidimensional scaling (NMDS) can be used to summarise information and reduce the dimensionality of profiles [58]. By applying NMDS, vineyards and taxa can be ordered by the dissimilarity of the presence profiles. Bray–Curtis dissimilarity [59], used extensively in the ecological field, was adopted. NMDS analysis was performed using the metaMDS function of the vegan package in R [60]. Loss of information due to a reduction in dimensionality is assessed by the stress value, which refers to the disagreement between 2-D representation and original positions of taxa in multidimensional space.
To test which environmental drivers (Cu, pH, AL, SOM, P, K, Mg, TL, TH, MD and MH) are significantly correlated to the first two axes of the NMDS ordination plane, we applied the envfit function of the vegan R package [60]. Each variable was correlated independently and plotted on the plane as a vector. The direction of the vector represents the gradient direction of the environmental driver, while the length of the vector is proportional to the correlation of the ordination system and the environmental driver.
Taxa were grouped into clusters as homogeneous as possible in terms of co-occurrence patterns, based on taxa ordination results [59]. To perform hierarchical cluster analysis, the hclust function of R software [61] was applied.

2.5.2. Vineyard Management Impact

To assess the impact of vineyard management on the biodiversity of soil biota, decision tree analysis was performed. The number of taxa present in each soil sample was considered as a measurement of edaphic biodiversity, and three categories of soil biodiversity were defined: ‘low’ when the number of taxa was lower or equal to 4, ‘medium’ when the number of taxa in the soil sample was between 5 and 8, and ‘high’ when the number of taxa was greater than 8. A classification decision tree allowed to split the soil samples into homogeneous groups according to edaphic biodiversity based on the different vineyard management classes. Recursive partitioning and regression tree (RPART) analysis were performed by applying the rpart package of R software [62]. The fitting of the model was investigated using the accuracy index that corresponds to the percentage of cases correctly classified.

3. Results

3.1. Descriptive Analysis

3.1.1. Environmental and Vineyard Management Variables

The descriptive statistics for environmental variables included in the full model are shown in Table 1.
Seven types of soil texture were considered: clay, clay loam, silty clay loam, sandy clay loam, loam, silt loam and sandy loam.
Vineyard management was categorised into four classes: conventional management (7% of the sample), vineyards converted to organic farming in the last three years (45% of the sample), vineyards converted between 4 and 9 years ago (31% of the sample), and vineyards converted at least 10 years ago (17% of the sample).

3.1.2. Taxa Identification

A total of 19 taxa were identified in the soil samples. In case of Diptera and Coleoptera, the biological stage of larvae were also detected (Table 2).
Collembola, Acari and Hymenoptera recorded the highest frequency of presence in the soil samples analysed. Collembola and Acari were reported in 89 of the 100 vineyards, Hymenoptera in 80 vineyards. The lowest frequency of occurrence was recorded for Psocoptera, Thysanoptera and Isopoda (8/100, 7/100, 6/100 respectively).

3.2. Co-Occurrence Pattern Identification

Taxa dispersion in the non-metric multidimensional scaling plane is shown in Figure 2. Taxa were ordered according to their co-occurrence profiles. Neighbouring taxa in the plane were characterised by the presence in the same vineyards (e.g., Collembola and Coleoptera larvae, Psocoptera and Pseudoscorpionida); the more distant are two taxa, greater is the difference in terms of their presence in the vineyards (e.g., Diptera and Psocoptera, Acari and Pauropoda). The stress value estimated for the model was equal to 0.2, indicating the model has good ability to predict data in the reduced space.
The results of analysis of the correlation between environmental drivers and the NMDS plane are shown in Figure 2.
p-Values of the correlation coefficients were used to discriminate the intensity of the relationship between environmental drivers and the taxa ordering system (Table 3): strong correlation for SOM, TL, TH and MH (p-value < 0.05); medium intensity correlation for pH (p-value < 0.1); low intensity correlation for MD (p-value < 0.15). The other environmental drivers were not significantly correlated with the first two axes of the NMDS system.
The results obtained from NMDS and cluster analysis (Figure 3) allowed the taxa to be divided into five groups according to their co-occurrence pattern. The five clusters shown in the cluster dendrogram correspond to the clusters identified by the green circles in the NMDS plane (Figure 2).
Group A included the largest number of taxa and specifically the Pseudoscorpionida, Psocoptera, Protura, Diplura Chilopoda, Symphyla and Pauropoda. Group B was made up of Diptera, Hemiptera and Isopoda taxa, while the larval form of Diptera was located in group D, together with Coleoptera, both as larvae and other biologic forms, and Collembola. The Acari, Hymenoptera, Thysanoptera and Diplopoda taxa made up group C. Group E is only represented by the taxa defined as ‘Other_holometabolous’.

3.3. Vineyard Management

The results of the classification tree showed that variable vineyard management could be useful for discriminating different categories of soil biodiversity. In particular, the analysis of the tree shown in Figure 4 showed that conventionally managed vineyards were associated with a low level of biodiversity, vineyards that had adopted organic management for a maximum of three years were associated with a medium level of biodiversity, and vineyards that had adopted organic management for at least four years were associated with a high level of biodiversity.
The accuracy index showed a good fit of the model as 57% of cases was correctly classified.

4. Discussion and Conclusions

The results obtained in this study allowed to identify the co-occurrence pattern for 19 taxa of soil arthropods on the basis of a 5-year investigation carried out in the Franciacorta viticultural area (Lombardy, Italy). The NMDS showed significant relationships between investigated soil arthropod taxa and soil moisture (MD and MH), soil temperature (TL, TH), soil organic matter (SOM) and pH. The decision tree showed an increased taxa diversity in relation to organic vineyard management and to the increase of time period of conversion from conventional to organic management.
In line with the expectations, Collembola and Acari were the most frequent of the 19 taxa identified, confirming that they are the most present groups of arthropods in soil [36,47,51]. The high level of presence of Hymenoptera recorded in our analysis is in agreement with other studies carried out in different agricultural contexts and reporting a significant presence of this taxon, mostly represented by Formicidae, in vineyard soils [63].
Based on taxa co-occurrence patterns, identified through NMDS analysis, five groups were found. Moreover, NMDS analysis made it possible to explore the relationship between soil abiotic variables and the aggregation of arthropod taxa in groups. In particular, the results obtained from our study pointed out that presence patterns characterising group A showed only one significant correlation (p-value < 0.15) with low soil moisture (MD). This result is compatible with the hypothesis that the taxa included in group A were relatively less dependent on high humidity values. The taxa in group B and the Coleoptera and Diptera larvae taxa (group D) were associated with higher pH (p-value < 0.1) and higher soil temperatures (TH) (p-value < 0.05), in line with the possible thermophilic habit of some representatives of these taxa [22,31,64]. The relationship with a higher pH level is more evident for the Isopoda and this is in agreement with van Straalen [11], who underlined weakly alkaliphilous or sub-neutral behaviour for some species of Isopoda. The detected ubiquitous presence of Collembola (group D) could be partially explained by the significant variability of responses to soil temperature, moisture and chemical properties of the different species of this taxon. In particular, the effect of soil moisture on Collembola has been documented by different authors [25,65], while species-specific responses have been reported [21]. Furthermore, Heiniger et al. [66] highlight that the role of microclimate for Collembola could be less important for their distribution than the role of trophic resources and competition. The presence of taxa in groups C and E is mostly determined by soil organic matter (SOM), soil moisture value higher than threshold level of 0.35 (MH) and lower temperature (TL) (p-values < 0.05). The relationship with SOM can be related to the involvement of these taxa in the soil food webs that starts from decomposition of dead organic matter generated by the activity of bacteria and fungi [32]. Diplopoda (Millipedes) are involved in SOM degradation, as their feeding activity is focused on dead organic matter [35,67]. A significant influence of soil nitrogen on species richness and biodiversity has been observed for this taxon [68], while Hymenoptera are involved in the decomposition of organic substances [35]. In relation to the positive response of group C to soil moisture increase, some authors have underlined that soil water availability is an important factor controlling presence of mites (Acari) [69]. Other authors have showed that Oribatid mites (Acari: Oribatida) are positively influenced by soil temperature [27] and that their distribution is dependent on soil moisture [35]. The relationship observed between group E and soil moisture can be related for Mecoptera (included in Other holometabolous taxon) with data reported for pre-imaginal stages of this order which develop in the soil and showed preference for high soil moisture [70].
The co-occurrence pattern of the taxa identified in our study is in line with similar pattern reported in the literature. Taxa co-occurrence in group B agrees with the results in [71] that confirmed Diptera and Isopoda co-existence in some specific habitats. Acari and Hymenoptera (group C) have also been grouped together by other authors [41]. The composition of groups A and C suggests that the co-occurrence pattern can also be influenced by biotic relationships among taxa. According to Eisenbeis and Wichard [31] (p. 192), the trophic niche of Diplura includes Symphyla, while Weygoldt [72] noted that Pseudoscorpionida feed on different orders of small soil arthropods, including Psocoptera. All these associations support the taxa co-occurrence in group A. Similarly, Coleoptera contain taxa (e.g., Carabid beetle) that have been described as predators of Collembola [73]. This association is in line with the co-occurrence of these two taxa in group D.
The results obtained analysing the role of vineyard management on soil arthropods diversity allows to identify an increase of taxa diversity in relation to organic vineyard management. This is consistent with previous studies, which reported a general increase of arthropod biodiversity [48] and arthropod abundance [47] associated to organic vineyard management. This effect was evident even before a 3-year period after conversion. The effect on arthropod biodiversity markedly increases with the length of the period since organic farming adoption.
The results obtained in this study provide additional knowledge supporting the interpretation of diversity and co-occurrence patterns in soil Arthropoda in vineyard. The importance of abiotic variables together with the interpretation of the possible role of biotic relationship among taxa have been explored in the specific geographic context of the Franciacorta viticultural area. Furthermore, our study confirmed the effect of organic vineyard management in increasing arthropod taxa diversity and, most importantly, it showed the critical role of the time of conversion from conventional to organic farming in increasing arthropod biodiversity. Further experiments are needed to extend these results to other viticultural contexts.

Author Contributions

Conceptualization, I.G., P.D., M.T., L.V. and G.G.; Data curation, I.G., A.S., F.O., M.T. and G.G.; Formal analysis, A.S.; Investigation, I.G., P.D., M.T. and L.V.; Methodology, I.G., A.S. and P.D.; Supervision, G.G.; Visualization, I.G. and A.S.; Writing—original draft, I.G. and A.S.; Writing—review & editing, G.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to thank all the wine companies assisting with data collection. We thank the ‘Consorzio per la tutela del Franciacorta’, which played a key role in the development of new operational protocols for sustainable viticulture by promoting and inspiring collaboration between universities and productive excellence. We thank the Chemical Laboratory and Organic Agriculture Unit of the Technology Transfer Center (CTT) of the Fondazione Edmund Mach (San Michele all’Adige, Trento, Italy), which carried out chemical analysis of soils and arthropods groups identification. We also thank TerrAria s.r.l. for providing high-resolution weather and soil data.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Table A1. Average values of soil temperature (ST °C) and soil moisture (SM measured in the interval 0–1) of the 30-days reference period prior to the sampling date for each site.
Table A1. Average values of soil temperature (ST °C) and soil moisture (SM measured in the interval 0–1) of the 30-days reference period prior to the sampling date for each site.
ID FieldDateSTSMID FieldDateSTSM
114 May20.200.0875115 June24.900.299
214 May20.210.0905215 June24.610.221
314 May19.650.0855315 June24.850.227
414 May19.720.0915415 June24.540.244
514 May20.270.1065515 June26.140.225
614 May20.310.0825615 June25.210.226
714 May20.280.0825715 June25.210.225
814 May20.190.0885815 June25.970.223
914 May20.180.0865915 June25.970.222
1014 May20.280.0876016 June21.710.340
1114 May20.290.0846116 June21.710.341
1214 May20.700.1006216 June21.710.340
1314 May20.710.1036316 June21.710.339
1414 May20.520.1016416 June25.750.370
1514 May20.520.1026516 June25.760.370
1614 May19.260.1286616 June20.290.344
1714 May19.260.1276716 June26.050.369
1814 May19.220.1456816 June25.750.370
1914 May19.530.1216916 June26.050.369
2014 May19.440.1207016 June26.050.369
2114 May19.580.1227116 June25.750.370
2214 May20.520.1007216 June26.020.369
2314 May20.520.1007316 June20.290.344
2414 May20.520.0997416 June26.050.369
2514 May19.630.0877516 June26.560.370
2614 May20.170.0877616 June20.290.344
2714 May20.170.0887716 June26.040.369
2814 September23.470.1907816 June21.540.337
2914 September23.580.3257916 June25.140.371
3014 September24.560.1888016 June21.380.326
3114 September24.560.1858116 June18.980.352
3215 May18.700.3538216 June18.980.350
3315 May19.970.3368316 June18.980.352
3415 June22.980.3618416 June18.980.354
3515 June22.980.3458516 June18.980.352
3615 June22.920.3618616 June18.980.351
3715 June24.740.2308717 September27.970.349
3815 June24.740.2288817 September27.840.349
3915 June25.020.3008917 September22.950.353
4015 June23.890.2759017 September21.670.164
4115 June24.740.1769117 September21.670.162
4215 June24.700.1789217 September21.670.159
4315 June24.700.1779317 September25.840.339
4415 June24.750.1659417 September25.840.340
4515 June25.080.2999517 September21.670.162
4615 June24.900.3029617 September25.840.340
4715 June24.140.3369717 September25.840.340
4815 June24.130.3389818 June24.370.323
4915 June24.140.3449918 June24.370.322
5015 June24.110.33410018 June24.370.323

References

  1. Kopittke, P.M.; Menzies, N.W.; Wang, P.; McKenna, B.A.; Lombi, E. Soil and the intensification of agriculture for global food security. Environ. Int. 2019, 132, 105078. [Google Scholar] [CrossRef] [PubMed]
  2. Wolters, V. Biodiversity of soil fauna and its function. Eur. J. Soil Biol. 2001, 37, 221–227. [Google Scholar] [CrossRef] [Green Version]
  3. Decaëns, T.; Jiménez, J.J.; Gioia, C.; Measey, G.J.; Lavelle, P. The values of soil animals for conservation biology. Eur. J. Soil Biol. 2006, 42, S23–S38. [Google Scholar] [CrossRef]
  4. European Commission. The Factory of Life. Why Soil Biodiversity Is so Important; Office for Official Publications of the European Union: Luxembourg, 2010; p. 22.
  5. Geisen, S.; Briones, M.J.I.; Gan, H.; Behan-Pelletier, V.M.; Friman, V.-P.; de Groot, G.A.; Hannula, S.E.; Lindo, Z.; Philippot, L.; Tiunov, A.V.; et al. A methodological framework to embrace soil biodiversity. Soil Biol. Biochem. 2019, 136, 107536. [Google Scholar] [CrossRef]
  6. Cameron, E.K.; Martins, I.S.; Lavelle, P.; Mathieu, J.; Tedersoo, L.; Gottschall, F.; Guerra, C.A.; Hines, J.; Patoine, G.; Siebert, J.; et al. Global gaps in soil biodiversity data. Nat. Ecol. Evol. 2018, 2, 1042–1043. [Google Scholar] [CrossRef]
  7. Juan-Ovejero, R.; Benito, E.; Barreal, M.E.; Rodeiro, J.; Briones, M.J.I. Tolerance to fluctuating water regimes drives changes in mesofauna community structure and vertical stratification in peatlands. Pedobiologia 2019, 76, 150571. [Google Scholar] [CrossRef]
  8. Wallwork, J.A. Distribution Patterns and Population Dynamics of the Micro-Arthropods of a Desert Soil in Southern California. J. Anim. Ecol. 1972, 41, 291–310. [Google Scholar] [CrossRef]
  9. Ruf, A.; Beck, L.; Dreher, P.; Hund-Rinke, K.; Römbke, J.; Spelda, J. A biological classification concept for the assessment of soil quality: “biological soil classification scheme” (BBSK). Agric. Ecosyst. Environ. 2003, 98, 263–271. [Google Scholar] [CrossRef]
  10. Marasas, M.E.; Sarandón, S.J.; Cicchino, A.C. Changes in soil arthropod functional group in a wheat crop under conventional and no tillage systems in Argentina. Appl. Soil Ecol. 2001, 18, 61–68. [Google Scholar] [CrossRef]
  11. van Straalen, N.M. Evaluation of bioindicator systems derived from soil arthropod communities. Appl. Soil Ecol. 1998, 9, 429–437. [Google Scholar] [CrossRef]
  12. Büchs, W.; Harenberg, A.; Zimmermann, J.; Weiß, B. Biodiversity, the ultimate agri-environmental indicator?: Potential and limits for the application of faunistic elements as gradual indicators in agroecosystems. Agric. Ecosyst. Environ. 2003, 98, 99–123. [Google Scholar] [CrossRef]
  13. Migliorini, M.; Pigino, G.; Bianchi, N.; Bernini, F.; Leonzio, C. The effects of heavy metal contamination on the soil arthropod community of a shooting range. Environ. Pollut. 2004, 129, 331–340. [Google Scholar] [CrossRef] [PubMed]
  14. Xu, G.-L.; Kuster, T.M.; Günthardt-Goerg, M.S.; Dobbertin, M.; Li, M.-H. Seasonal exposure to drought and air warming affects soil Collembola and mites. PLoS ONE 2012, 7, 1–9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Ghiglieno, I.; Simonetto, A.; Donna, P.; Tonni, M.; Valenti, L.; Bedussi, F.; Gilioli, G. Soil Biological Quality Assessment to Improve Decision Support in the Wine Sector. Agronomy 2019, 9, 593. [Google Scholar] [CrossRef] [Green Version]
  16. Frampton, G.; Van den Brink, P.; Gould, P. Effect of spring drought and irrigation on farmland arthropods in southern Britain. J. Appl. Ecol. 2001, 37, 865–883. [Google Scholar] [CrossRef]
  17. Uvarov, A.V. Effects of diurnal temperature fluctuations on population responses of forest floor mites. Pedobiologia 2003, 47, 331–339. [Google Scholar] [CrossRef]
  18. Jucevica, E.; Melecis, V. Global warming affect Collembola community: A long-term study. Pedobiologia 2006, 50, 177–184. [Google Scholar] [CrossRef]
  19. Clapperton, M.J.; Kanashiro, D.A.; Behan-Pelletier, V.M. Changes in abundance and diversity of microarthropods associated with Fescue Prairie grazing regimes. Pedobiologia 2002, 46, 496–511. [Google Scholar] [CrossRef]
  20. Holmstrup, M.; Maraldo, K.; Krogh, P.H. Combined effect of copper and prolonged summer drought on soil Microarthropods in the field. Environ. Pollut. 2007, 146, 525–533. [Google Scholar] [CrossRef]
  21. Tsiafouli, M.A.; Kallimanis, A.S.; Katana, E.; Stamou, G.P.; Sgardelis, S.P. Responses of soil microarthropods to experimental short-term manipulations of soil moisture. Appl. Soil Ecol. 2005, 29, 17–26. [Google Scholar] [CrossRef]
  22. Reddy, M.V.; Venkataiah, B. Seasonal abundance of soil-surface arthropods in relation to some meteorological and edaphic variables of the grassland and tree-planted areas in a tropical semi-arid savanna. Int. J. Biometeorol. 1990, 34, 49–59. [Google Scholar] [CrossRef]
  23. Grear, J.S.; Schmitz, O.J. Effects of Grouping Behavior and Predators on the Spatial Distribution of a Forest Floor Arthropod. Ecology 2005, 86, 960–971. [Google Scholar] [CrossRef] [Green Version]
  24. Choi, W.I.; Ryoo, M.I.; Kim, J.-G. Biology of Paronychiurus kimi(Collembola: Onychiuridae) under the influence of temperature, humidity and nutrition. Pedobiologia 2002, 46, 548–557. [Google Scholar] [CrossRef]
  25. Choi, W.I.; Moorhead, D.L.; Neher, D.A.; Ryoo, M.I. A modeling study of soil temperature and moisture effects on population dynamics of Paronychiurus kimi (Collembola: Onychiuridae). Biol. Fertil. Soils 2006, 43, 69–75. [Google Scholar] [CrossRef]
  26. O’Lear, H.A.; Blair, J.M. Responses of soil microarthropods to changes in soil water availability in tallgrass prairie. Biol. Fertil. Soils 1999, 29, 207–217. [Google Scholar] [CrossRef]
  27. Liu, J.-L.; Li, F.-R.; Liu, L.-L.; Yang, K. Responses of different Collembola and mite taxa to experimental rain pulses in an arid ecosystem. Catena 2017, 155, 53–61. [Google Scholar] [CrossRef]
  28. Ikemoto, T. Intrinsic Optimum Temperature for Development of Insects and Mites. Environ. Entomol. 2005, 34, 1377–1387. [Google Scholar] [CrossRef]
  29. Ikemoto, T. Possible existence of a common temperature and a common duration of development among members of a taxonomic group of arthropods that underwent speciational adaptation to temperature. Appl. Entomol. Zool. 2003, 38, 487–492. [Google Scholar] [CrossRef] [Green Version]
  30. Rapoport, E.H.; Tschapek, M. Soil water and soil fauna. Rev. Ecol. Biol. Sol. 1967, 4, 1–58. [Google Scholar]
  31. Eisenbeis, G.; Wichard, W. Atlas on the Biology of Soil Arthropods; Springer: Berlin/Heidelberg, Germany, 1987. [Google Scholar]
  32. Petersen, H.; Luxton, M. A Comparative Analysis of Soil Fauna Populations and Their Role in Decomposition Processes. Oikos 1982, 39, 288–388. [Google Scholar] [CrossRef]
  33. Wu, P.; Liu, X.; Liu, S.; Wang, J.; Wang, Y. Composition and spatio-temporal variation of soil microarthropods in the biodiversity hotspot of northern Hengduan Mountains, China. Eur. J. Soil Biol. 2014, 62, 30–38. [Google Scholar] [CrossRef]
  34. Culliney, T. Role of arthropods in maintaining soil fertility. Agriculture 2013, 3, 629–659. [Google Scholar] [CrossRef] [Green Version]
  35. Bagyaraj, D.J.; Nethravathi, C.J.; Nitin, K.S. Soil Biodiversity and Arthropods: Role in Soil Fertility. In Economic and Ecological Significance of Arthropods in Diversified Ecosystems; Chakravarthy, A.K., Sridhara, S., Eds.; Springer: Singapore, 2016; pp. 17–51. [Google Scholar]
  36. André, H.M.; Noti, M.-I.; Lebrun, P. The soil fauna: The other last biotic frontier. Biodivers. Conserv. 1994, 3, 45–56. [Google Scholar] [CrossRef]
  37. van Straalen, N.M.; Verhoef, H.A. The Development of a Bioindicator System for Soil Acidity Based on Arthropod pH Preferences. J. Appl. Ecol. 1997, 34, 217–232. [Google Scholar] [CrossRef]
  38. Lavelle, P.; Chauvel, A.; Fragoso, C. Faunal activity in acid soils. In Plant-Soil Interactions at Low pH: Principles and Management. Developments in Plant and Soil Sciences; Date, R.A., Grundon, N.J., Rayment, G.E., Probert, M.E., Eds.; Springer: Dordrecht, The Netherlands, 1995; Volume 64, pp. 201–211. [Google Scholar]
  39. Andrés, P.; Moore, J.C.; Simpson, R.T.; Selby, G.; Cotrufo, F.; Denef, K.; Haddix, M.L.; Shaw, E.A.; de Tomasel, C.M.; Molowny-Horas, R.; et al. Soil food web stability in response to grazing in a semi-arid prairie: The importance of soil textural heterogeneity. Soil Biol. Biochem. 2016, 97, 131–143. [Google Scholar] [CrossRef]
  40. Potapov, A.; Goncharov, A.; Semenina, E.; Korotkevich, A.; Tsurikov, S.; Rozanova, O.; Anichkin, A.; Zuev, A.; Samoylova, E.; Semenyuk, I.; et al. Arthropods in the subsoil: Abundance and vertical distribution as related to soil organic matter, microbial biomass and plant roots. Eur. J. Soil Biol. 2017, 82, 88–97. [Google Scholar] [CrossRef]
  41. Shakir, M.M.; Ahmed, S. Seasonal abundance of soil arthropods in relation to meteorological and edaphic factors in the agroecosystems of Faisalabad, Punjab, Pakistan. Int. J. Biometeorol. 2015, 59, 605–616. [Google Scholar] [CrossRef]
  42. Sapkota, T.B.; Mazzoncini, M.; Bàrberi, P.; Antichi, D.; Silvestri, N. Fifteen years of no till increase soil organic matter, microbial biomass and arthropod diversity in cover crop-based arable cropping systems. Agron. Sustain. Dev. 2012, 32, 853–863. [Google Scholar] [CrossRef]
  43. Bedano, J.C.; Cantú, M.P.; Doucet, M.E. Soil springtails (Hexapoda: Collembola), symphylans and pauropods (Arthropoda: Myriapoda) under different management systems in agroecosystems of the subhumid Pampa (Argentina). Eur. J. Soil Biol. 2006, 42, 107–119. [Google Scholar] [CrossRef]
  44. Chauvat, M.; Wolters, V.; Dauber, J. Response of collembolan communities to land-use change and grassland succession. Ecography 2007, 30, 183–192. [Google Scholar] [CrossRef]
  45. Gkisakis, V.; Volakakis, N.; Kollaros, D.; Bàrberi, P.; Kabourakis, E.M. Soil arthropod community in the olive agroecosystem: Determined by environment and farming practices in different management systems and agroecological zones. Agric. Ecosyst. Environ. 2016, 218, 178–189. [Google Scholar] [CrossRef]
  46. Hadjicharalampous, E.; Kalburtji, K.L.; Mamolos, A.P. Soil Arthropods (Coleoptera, Isopoda) in Organic and Conventional Agroecosystems. Environ. Manag. 2002, 29, 683–690. [Google Scholar] [CrossRef] [PubMed]
  47. Gagnarli, E.; Goggioli, D.; Tarchi, F.; Guidi, S.; Nannelli, R.; Vignozzi, N.; Valboa, G.; Lottero, M.R.; Corino, L.; Simoni, S. Case study of microarthropod communities to assess soil quality in different managed vineyards. Soil 2015, 1, 527–536. [Google Scholar] [CrossRef] [Green Version]
  48. Caprio, E.; Nervo, B.; Isaia, M.; Allegro, G.; Rolando, A. Organic versus conventional systems in viticulture: Comparative effects on spiders and carabids in vineyards and adjacent forests. Agric. Syst. 2015, 136, 61–69. [Google Scholar] [CrossRef]
  49. Puig-Montserrat, X.; Stefanescu, C.; Torre, I.; Palet, J.; Fàbregas, E.; Dantart, J.; Arrizabalaga, A.; Flaquer, C. Effects of organic and conventional crop management on vineyard biodiversity. Agric. Ecosyst. Environ. 2017, 243, 19–26. [Google Scholar] [CrossRef]
  50. Seniczak, A.; Seniczak, S.; García-Parra, I.; Ferragut, F.; Xamaní, P.; Graczyk, R.; Messeguer, E.; Laborda, R.; Rodrigo, E. Oribatid mites of conventional and organic vineyards in the Valencian Community, Spain. Acarologia 2018, 58, 119–133. [Google Scholar]
  51. Costantini, E.A.C.; Agnelli, A.E.; Fabiani, A.; Gagnarli, E.; Mocali, S.; Priori, S.; Simoni, S.; Valboa, G. Short-term recovery of soil physical, chemical, micro- and mesobiological functions in a new vineyard under organic farming. Soil 2015, 1, 443–457. [Google Scholar] [CrossRef] [Green Version]
  52. Döring, J.; Collins, C.; Frisch, M.; Kauer, R. Organic and Biodynamic Viticulture Affect Biodiversity and Vine and Wine Properties: A Systematic Quantitative Review. Am. J. Enol. Vitic. 2019, 70, 221–242. [Google Scholar] [CrossRef]
  53. National Centers for Environmental Prediction; National Weather Service; NOAA; U.S. Department of Commerce. NCEP GDAS/FNL 0.25 Degree Global Tropospheric Analyses and Forecast Grids 2015. Research Data Archive at the National Center for Atmospheric Research, Computational and Information Systems Laboratory. Available online: https://0-doi-org.brum.beds.ac.uk/10.5065/D65Q4T4Z (accessed on 5 November 2019).
  54. Powers, J.G.; Klemp, J.B.; Skamarock, W.C.; Davis, C.A.; Dudhia, J.; Gill, D.O.; Coen, J.L.; Gochis, D.J.; Ahmadov, R.; Peckham, S.E.; et al. The Weather Research and Forecasting Model: Overview, System Efforts, and Future Directions. Bull. Am. Meteorol. Soc. 2017, 98, 1717–1737. [Google Scholar] [CrossRef]
  55. Ek, M.B.; Mitchell, K.E.; Lin, Y.; Rogers, E.; Grunmann, P.; Koren, V.; Gayno, G.; Tarpley, J.D. Implementation of Noah land surface model advances in the National Centers for Environmental Prediction operational mesoscale Eta model. J. Geophys. Res. Atmos. 2003, 108, 8851. [Google Scholar] [CrossRef]
  56. Soil Science Division Staff. Soil Survey Manual; Ditzler, C., Scheffe, K., Monger, H.C., Eds.; USDA Handbook 18; Government Printing Office: Washington, DC, USA, 2017; p. 125.
  57. Menta, C.; Conti, F.D.; Pinto, S.; Bodini, A. Soil Biological Quality index (QBS-ar): 15 years of application at global scale. Ecol. Indic. 2018, 85, 773–780. [Google Scholar] [CrossRef]
  58. Kenkel, N.C.; Orloci, L. Applying Metric and Nonmetric Multidimensional Scaling to Ecological Studies: Some New Results. Ecology 1986, 67, 919–928. [Google Scholar] [CrossRef]
  59. Clarke, K.R. Non-parametric multivariate analyses of changes in community structure. Aust. J. Ecol. 1993, 18, 117–143. [Google Scholar] [CrossRef]
  60. Oksanen, J.; Blanchet, F.G.; Friendly, M.; Kindt, R.; Legendre, P.; McGlinn, D.; Minchin, P.R.; O’Hara, R.B.; Simpson, G.L.; Solymos, P.; et al. Vegan: Community Ecology Package, R Package Version 2.5; 2018. Available online: http://CRAN.Rproject.org/package=vegan (accessed on 19 May 2020).
  61. R Core Team. R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing; R Core Team: Vienna, Austria, 2019; Available online: https://www.R-project.org/ (accessed on 19 May 2020).
  62. Therneau, T.; Atkinson, B. rpart: Recursive Partitioning and Regression Trees, R Package Version 4.1-13; 2018. Available online: http://CRAN.Rproject.org/package=rpart (accessed on 19 May 2020).
  63. Pérez Bote, J.; Romero, A. Epigeic soil arthropod abundance under different agricultural land uses. Span. J. Agric. Res. 2012, 10, 55–61. [Google Scholar] [CrossRef] [Green Version]
  64. Zhu, G.; Luo, Y.; Xue, M.; Zhao, H.; Xia, N.; Wang, X. Effects of high-temperature stress and heat shock on two root maggots, Bradysia odoriphaga and Bradysia difformis (Diptera: Sciaridae). J. Asia-Pac. Entomol. 2018, 21, 106–114. [Google Scholar] [CrossRef]
  65. Sjursen, H.; Holmstrup, M. Cold and drought stress in combination with pyrene exposure: Studies with Protaphorura armata (Collembola: Onychiuridae). Ecotoxicol. Environ. Saf. 2004, 57, 145–152. [Google Scholar] [CrossRef]
  66. Heiniger, C.; Barot, S.; Ponge, J.-F.; Salmon, S.; Meriguet, J.; Carmignac, D.; Suillerot, M.; Dubs, F. Collembolan preferences for soil and microclimate in forest and pasture communities. Soil Biol. Biochem. 2015, 86, 181–192. [Google Scholar] [CrossRef]
  67. Bertrand, M.; Lumaret, J.P. The role of Diplopoda litter grazing activity on recycling processes in a Mediterranean climate. Vegetatio 1992, 99, 289–297. [Google Scholar] [CrossRef]
  68. Stašiovstasiov, S.; Stašiová, A.; Svitok, M.; Michalková, E.; Slobodník, B.; Lukáčikarbbh, I. Millipede (Diplopoda) communities in an arboretum: Influence of tree species and soil properties. Biologia 2012, 67, 945–952. [Google Scholar]
  69. Badejo, M.A. Seasonal Abundance of Soil Mites (Acarina) in Two Contrasting Environments. Biotropica 1990, 22, 382–390. [Google Scholar] [CrossRef]
  70. Byers, G.W.; Thornhill, R. Biology of the Mecoptera. Annu. Rev. Entomol. 1983, 28, 203–228. [Google Scholar] [CrossRef]
  71. Varga, I. Structure and changes of macroinvertebrate community colonizing decomposing rhizome litter of common reed at Lake Fertő/Neusiedler See (Hungary). Hydrobiologia 2003, 506, 413–420. [Google Scholar] [CrossRef]
  72. Weygoldt, P. The Biology of Pseudoscorpions; Harvard University Press: Cambridge, Massachusetts, 1969. [Google Scholar]
  73. Bauer, T. Predation by a carabid beetle specialized for catching collembola. Pedobiologia 1982, 24, 169–179. [Google Scholar]
Figure 1. Map of the Franciacorta DOCG (Designation of Controlled and Guarantee Origin) winegrowing area. The locations of vineyards where samples were collected are indicated with red dots.
Figure 1. Map of the Franciacorta DOCG (Designation of Controlled and Guarantee Origin) winegrowing area. The locations of vineyards where samples were collected are indicated with red dots.
Agronomy 10 00740 g001
Figure 2. Results of non-metric multidimensional (NMDS) analysis: dispersion of taxa (red points) according to their co-occurrence profiles (NMDS1 and NMDS2 are the two axes of the ordination plane). Blue arrows refer to the correlation of environmental drivers and soil characteristics with NMDS ordination pattern (solid line Agronomy 10 00740 i001 p-value < 0.5, dashed line Agronomy 10 00740 i002 p-value < 0.1, dotted line Agronomy 10 00740 i003 p-value < 0.15). The five clusters of taxa according their presence pattern are highlighted with the green circles.
Figure 2. Results of non-metric multidimensional (NMDS) analysis: dispersion of taxa (red points) according to their co-occurrence profiles (NMDS1 and NMDS2 are the two axes of the ordination plane). Blue arrows refer to the correlation of environmental drivers and soil characteristics with NMDS ordination pattern (solid line Agronomy 10 00740 i001 p-value < 0.5, dashed line Agronomy 10 00740 i002 p-value < 0.1, dotted line Agronomy 10 00740 i003 p-value < 0.15). The five clusters of taxa according their presence pattern are highlighted with the green circles.
Agronomy 10 00740 g002
Figure 3. Dendrogram of hierarchical cluster analysis of taxa based on NMDS results. The five clusters are highlighted in green.
Figure 3. Dendrogram of hierarchical cluster analysis of taxa based on NMDS results. The five clusters are highlighted in green.
Agronomy 10 00740 g003
Figure 4. Classification decision tree of soil samples to predict soil biodiversity according to vineyard management. The predicted level of soil biodiversity (low, medium, high) is reported in the squared box, together with the percentage of soil samples included in that node. The paths from the initial box (with 100% of cases) to the final boxes represent the classification rules.
Figure 4. Classification decision tree of soil samples to predict soil biodiversity according to vineyard management. The predicted level of soil biodiversity (low, medium, high) is reported in the squared box, together with the percentage of soil samples included in that node. The paths from the initial box (with 100% of cases) to the final boxes represent the classification rules.
Agronomy 10 00740 g004
Table 1. Descriptive statistics of continuous variables (soil characteristics and environmental drivers) included in the analysis.
Table 1. Descriptive statistics of continuous variables (soil characteristics and environmental drivers) included in the analysis.
Unit of MeasureMean ± Standard DeviationMedian (Q25–Q75)MinMax
Cu(mg/kg)58.68 ± 32.8155.40 (36.9–72.2)4.20170.00
pH 7.10 ± 0.877.30 (6.35–7.9)5.308.20
AL(g CaCO3/kg)10.16 ± 22.620.00 (0.00–13.50)0.00130.00
SOM 1(g/kg)21.94 ± 9.0823.00 (15.00–25.00)5.0042.00
P(mg P2O5/kg)54.47 ± 40.2051.00 (26.00–64.00)9.00222.00
K(mg K2O/kg)148.52 ± 67.53145.00 (94.00–178.00)60.00354.00
Mg(mg MgO/kg)165.75 ± 75.37138.00 (117.00–210.00)66.00383.00
TL 2°C68.76 ± 44.9069.75 (33.67–104.39)0.00161.46
TH 3°C106.05 ± 59.06121.76 (55.73–153.38)14.20241.77
MD 4Pure number10.39 ± 26.035.16 (1.51–7.88)0.00135.53
MH 5Pure number0.11 ± 0.200.00 (0.00–0.156)0.000.63
1 SOM: soil organic matter. 2 TL: daily cumulative soil temperature exceeding 10 °C when soil temperature is between 10 and 20 °C. 3 TH: daily cumulative soil temperature exceeding 20 °C when soil temperature is higher than 20 °C. 4 MD: daily sum of absolute deviations in soil moisture when soil moisture is lower than 0.35. 5 MH: daily cumulative soil moisture exceeding 0.35 when soil moisture is higher than 0.35.
Table 2. Distribution of taxa according to stages considered in the analysis and presence (i.e., number of soil samples in which the taxon has been identified).
Table 2. Distribution of taxa according to stages considered in the analysis and presence (i.e., number of soil samples in which the taxon has been identified).
TaxaLarvaeN° of PresencesOther Stages 1N° of Presences
Acari x89
Myriapoda—Diplopoda x12
Myriapoda—Chilopoda x17
Myriapoda—Symphyla x56
Myriapoda—Pauropoda x32
Hymenoptera x80
Thysanoptera x7
Pseudoscorpionida x11
Psocoptera x8
Coleoptera x31
Coleoptera larvaex39
Collembola x89
Diptera x31
Diptera larvaex32
Protura x25
Diplura x27
Hemiptera x12
Isopoda x6
Other_holometabolous2 x20
1 Other stages include all forms that produce active participation in soil cycles (e.g., pupae are excluded). In the case of the ‘Other_holometabolous’ taxon, the pupal stage is also included. 2 Other_holometabolous taxa include Mecoptera, Neuroptera and Raphidioptera orders in agreement with QBS-ar (Soil Biological Quality-arthropod) method [57].
Table 3. Correlation analysis of environmental drivers and soil characteristics with NMDS ordination pattern.
Table 3. Correlation analysis of environmental drivers and soil characteristics with NMDS ordination pattern.
VariableSquared Correlation Coefficientp-Value 6 of Correlation Coefficient
Cu0.050.17
pH0.060.09**
AL0.020.48
SOM 10.080.05***
P0.020.43
K0.010.66
Mg0.010.64
TL 20.080.04***
TH 30.150.01***
MD 40.050.15*
MH 50.120.01***
TXT0.040.83
1 SOM: soil organic matter. 2 TL: daily cumulative soil temperature exceeding 10 °C when soil temperature is between 10 and 20 °C. 3 TH: daily cumulative soil temperature exceeding 20 °C when soil temperature is higher than 20 °C. 4 MD: daily sum of absolute deviations in soil moisture when soil moisture is lower than 0.35. 5 MH: daily cumulative soil moisture exceeding 0.35, when soil moisture is higher than 0.35. 6 * p-value < 0.15, ** p-value < 0.1, *** p-value < 0.05.

Share and Cite

MDPI and ACS Style

Ghiglieno, I.; Simonetto, A.; Orlando, F.; Donna, P.; Tonni, M.; Valenti, L.; Gilioli, G. Response of the Arthropod Community to Soil Characteristics and Management in the Franciacorta Viticultural Area (Lombardy, Italy). Agronomy 2020, 10, 740. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy10050740

AMA Style

Ghiglieno I, Simonetto A, Orlando F, Donna P, Tonni M, Valenti L, Gilioli G. Response of the Arthropod Community to Soil Characteristics and Management in the Franciacorta Viticultural Area (Lombardy, Italy). Agronomy. 2020; 10(5):740. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy10050740

Chicago/Turabian Style

Ghiglieno, Isabella, Anna Simonetto, Francesca Orlando, Pierluigi Donna, Marco Tonni, Leonardo Valenti, and Gianni Gilioli. 2020. "Response of the Arthropod Community to Soil Characteristics and Management in the Franciacorta Viticultural Area (Lombardy, Italy)" Agronomy 10, no. 5: 740. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy10050740

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