Next Article in Journal
The Role of Education in Sustainable Dietary Patterns in Slovenia
Next Article in Special Issue
Estimating Productivity, Detecting Biotic Disturbances, and Assessing the Health State of Traditional Olive Groves, Using Nondestructive Phenotypic Techniques
Previous Article in Journal
Framework Proposal for Achieving Smart and Sustainable Societies (S3)
Previous Article in Special Issue
Microaerobic Digestion of Low-Biodegradable Sewage Sludge: Effect of Air Dosing in Batch Reactors
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Patterns and Scales of Collembola Taxonomic and Functional Diversity in Urban Parks

by
Alessandro Bellino
1,
Daniela Baldantoni
1,*,
Vittoria Milano
2,3,
Lucia Santorufo
2,4,
Jérôme Cortet
3,† and
Giulia Maisto
2,4,†
1
Department of Chemistry and Biology “Adolfo Zambelli”, University of Salerno, Via Giovanni Paolo II 132, 84084 Fisciano, Italy
2
Department of Biology, University of Naples Federico II, Via Cinthia, 80126 Napoli, Italy
3
Centre d’Ecologie Fonctionnelle et Evolutive (CEFE) UMR5175, University Paul Valéry Montpellier, University of Montpellier, EPHE, CNRS, IRD, Route de Mende, 34000 Montpellier, France
4
BAT Center-Interuniversity Center for Studies on Bioinspired Agro-Environmental Technology, University of Naples Federico II, 80100 Naples, Italy
*
Author to whom correspondence should be addressed.
These authors equally contributed to the research.
Sustainability 2021, 13(23), 13029; https://0-doi-org.brum.beds.ac.uk/10.3390/su132313029
Submission received: 22 October 2021 / Revised: 19 November 2021 / Accepted: 21 November 2021 / Published: 25 November 2021

Abstract

:
Anthropogenic pressures can affect the distribution of species and elicit the appearance of spatial patterns that provide insights into the species’ responses to environmental filtering, mediated by their functional traits. Due to the functional redundancy in ecological communities, the spatial variations of species and functional traits can occur at different scales, but little is known on this topic, especially for soil arthropods with limited dispersion capabilities and highly dependent on environmental characteristics. The present research aims at shedding light on the spatial ecology of both the taxonomic and functional biodiversity of collembolan communities colonizing urban parks, adopted as model taxa for their functional traits, diversity and sensitivity to environmental drivers. To this end, the spatial patterns and scales of collembolan communities from 8 parks in Naples (Italy) and 14 in Montpellier (France) were investigated through an approach based on Moran eigenvector maps, modified to allow for evaluating the community spatial connectivity and the scales underpinning the spatial variation of each species and functional trait. The obtained findings demonstrate a limited spatial connectivity of collembolan communities in terms of both taxonomic and functional diversity, with mostly species-specific micro-scale variations that may be shaped by environmental constraints.

1. Introduction

The sustainable development of urban ecosystems pivots on networks of green infrastructures that allow preserving key ecosystem processes in otherwise constructed environments [1]. The provision of suitable habitats for flora and fauna, the regulation of organic matter and biogeochemical cycles, the climate mitigation and the environmental remediation are just some of the showcase processes supported by semi-natural environments in urban ecosystems [2,3]. At landscape scale, the interaction between these areas and the surrounding environment can be analyzed within the classical framework of island biogeography [4,5,6], with the urban tissue representing a hostile matrix that surrounds patches of habitat suitable for most of the flora and fauna. Urban parks represent most of these patches, with characteristics shaped by their size, topology, management and interaction with the surrounding matrix. Even in large parks, with limited protrusion of the urban tissue within their borders, management can have pervasive effects on their capability to act as biodiversity reservoirs. For example, complex canopies are able to host richer communities and counteract biotic homogenization [7,8,9]. Conversely, canopy simplification or reduction and the realization of passageways or fences can constrain the movements and dispersal of organisms. Management-induced fragmentation within urban parks has dramatic effects, especially on organisms with small body size and limited movements, which are common traits within soil communities [10]. The perceived spatial scale is what exacerbates the effects of management on these organisms that may not only act at small temporal scale, but also develop into temporal legacies exerting their effects through decades or centuries. A clear example is the demonstrated effect of park history or age on several taxa, including enchytraeids, nematodes and springtails [11,12,13,14].
On a theoretical level, environmental drivers may exert different effects on community taxonomic and functional diversities, as well as on their spatial variations. However, in spite of the well-documented effects on the taxonomic diversity, those on the functional diversity are comparatively less studied, especially in the case of soil fauna [15]. In this context, even fewer studies tackle the issue from a spatial perspective, investigating the connectivity of communities and their scales of spatial variation.
The present research aims at shedding light on these topics through a comparative analysis of the spatial patterns and scales of the taxonomic and functional biodiversity of soil collembolan communities of urban parks, and the evaluation of the main drivers underpinning these patterns. The diversity in functional traits expressed by this group, with large variations in body size, feeding strategies, colonized environments, speed of movements or dispersal, as well as their sensitivity to environmental drivers, make Collembola exceptional models for such purposes. Under the community convergence hypothesis [16], the study focuses in parallel on two urban systems, in southern Italy and southern France, in order to strengthen the comparative analysis by stressing the differential responses of species vs. functional assemblages to the geographical distance.

2. Materials and Methods

2.1. Study Area and Sampling Design

The study was carried out in 8 urban parks of Naples (40°50′ N; 14°15′ E, Italy) and 14 of Montpellier (43°36’ N; 03°52′ E, France), two cities lying on the coast of the western Mediterranean Sea and characterized by similar climate, but ~900 Km apart from each other (great circle distance).
Naples is an historical city and current megalopolis mainly lying on Andosols [17], characterized by an uninterrupted urban tissue, whereas Montpellier is a recently expanded city mainly lying on Calcisols [17]. According to the Köppen climate classification, both sites are classified as Csa (hot climate with drought summer). In Naples, mean temperatures range from 9.7 °C in winter to 25.7 °C in summer, with a mean annual temperature of 17.4 °C, and the annual precipitation is ~865 mm (http://www.meteo.unina.it/clima-di-napoli; accessed on 19 November 2021). In Montpellier, mean temperatures range from 5.9 °C in winter to 22.3 °C in summer, with a mean annual temperature of 13.9 °C, and the annual precipitation is ~745 mm (https://fr.climate-data.org/location/1126/; accessed on 19 November 2021). Woody habitats within the parks are dominated by Quercus ilex L. and, to a lesser extent, by Quercus petrea (Matt.) Lielb., Platanus occidentalis L. and Pinus pinaster Aiton. Brushwood layer, when present, is mostly represented by Laurus nobilis L., Nerium oleander L. and Ruscus aculeatus L.
Soil sampling, microarthropod extraction, Collembola species identification and counting, as well as the determination of environmental parameters and their coding are extensively reported in Milano et al. [13,14] and briefly reported hereafter.
Soil samplings were collected [18] after litter removal (if present) from four different plots (one sample per plot), either under canopy or from open spaces. Samples were then stored at 4 °C until microarthropod extraction, using a Macfaiden apparatus [19] and the following identification and counting of Collembola specimens.
In order to describe the environment experienced by Collembola, several descriptors of either park or sample properties, summarized in Table S1, were adopted. In particular, the age since park opening to the public (very old, old, young), previous land use (forest, wilderness, tuff quarry, agricultural land), elevation and urban density, where employed at the park level, coupled to characteristics of the site where samples were collected: presence of canopy, brushwood and litter; soil chemical and physical properties. The latter encompassed dissolved organic carbon and C/N ratio, as proxies of quantity and quality of trophic resources, and pH and bulk density as general descriptors of soil characteristics. The determination of soil properties was performed on the soil cores after microarthropod extraction and oven drying, according to the protocols described in Milano et al. [13,14]. Species that were observed in more than 2% of the samples, reported in Table 1, were adopted for all the spatial analyses.

2.2. Species Traits

Collembola species traits were obtained from the Betsi online database (http://betsi.cesab.org/; accessed on 16 September 2021) for soil invertebrates [20], integrated with information from various dichotomous keys [21,22,23,24,25,26,27], coded according to the Betsi rules. Each species was described by morphological proxies of life-form sensu Gisin [28]: body length, shape, pigmentation, presence of visual organs, scales, postantennal organ (PAO) and length of the furca. The set of functional traits adopted aimed at including different aspects of the species’ ecological niches that may be directly or indirectly related to their spatial distribution. Indeed, whereas body length is an overall proxy of food resource use, traits such as body shape, the presence of pigmentation for camouflage or as a means of protection against excessive light and the presence of organs for sensing the environment such as ocelli or PAO are directly related to the habitat preferences [29]. In terms of movements, instead, the length of the furca is directly involved in dispersion and predation evasion, to which the presence of scales also contributes [30].
In order to limit the dimensionality of the functional traits dataset and compress the information into few informative variables, the traits were coded as in Table 2 and transformed into dummy variables (presence–absence) according to Podani [31]. In the case of multilevel variables (body and furca lengths), each level defined a new presence–absence variable. Traits were then employed in calculating the community weighted mean (CWM) as a measure of functional composition of communities in each soil sample [32].

2.3. Data Analysis

The spatial analysis relied on the algorithm described in Bellino and coworkers [33], allowing investigation of the spatial organization and scales of variation of each species or trait within communities. In particular, 88 spatial weighting matrices (SWMs) were defined according to 4 connectivity diagrams with varying degrees of connectedness (Delaunay triangulation, Gabriel’s graph, minimum spanning tree and relative neighbors graph) and 22 weighting schemes. The schemes comprised linear, concave and convex functions, the second and third with varying exponential parameters in the set η = {0, 0.1, 0.2, …, 1.0}. The selection of the SWMs describing the spatial organization of communities and CWMs in Naples and Montpellier relied on the definition of least absolute shrinkage and selection operator (LASSO) models relating each species or trait to the Moran eigenvector maps (MEM) derived through eigen analysis of each SWM. The models relative to the species and the CWMs will be hereafter referred to as S-models and T-models, respectively. The active set of the models was chosen based on the Akaike information criterion (AIC) profile. The AICs relative to the active sets of the LASSO models were then ranked and aggregated by calculating the mean rank, and the SWM identified by the lowest mean rank was chosen as the consensus SWM, describing the spatial patterns of either species or CWMs. The MEMs derived from the consensus SWM associated to non-zero LASSO coefficients were considered the spatial scales describing the spatial variations of each species and functional trait.
As a means to detect the presence of relationships between the number of MEMs and the spatial distribution of species, i.e., the hypothesis of common species being described by more spatial scales, the Spearman correlation between the number of parks in which species in Naples and Montpellier were observed and the number of their associated MEMs was adopted.
Finally, in order to evaluate the effects of environmental drivers on the spatial scales of species and traits, an approach based on fourth corner analysis was adopted. In particular, the analysis was carried out using either species or CWMs as the community elements, the MEMs associated to each species and CWM as functional traits, and the descriptors of park/sample characteristics as environmental variables. A permutational approach was adopted in deriving P-values, which were adjusted using a false discovery rate criterion.
All the analyses were carried out within the R 4.0.5 [34] programming environment, using the core functions employed in Bellino et al. [33] and the function “fourthcorner” from the ade4 package [35].

3. Results

The mean rank of each SWM based on either the species or the CWMs in Naples and Montpellier is shown in Figure 1. Overall, the mean ranks show a consistent trend of decreasing values associated to concave weighting schemes and high η, with a progressive reduction in their differentiation among topologies. The consensus SWM, always weighted according to a concave function, has different topologies and η: relative neighbor graph with η = 0.9 and Gabriel’s graph with η = 1.0 in the case of the S-models for Naples and Montpellier, respectively, and Gabriel’s graph with η = 0.2 and Delaunay triangulation with η = 1.0 in the case of the T-models. The MEMs derived from the consensus SWMs are shown in Figure 2. Approximately half of the species and all the CWMs in both Naples and Montpellier are associated to non-zero LASSO coefficients, i.e., have non-random spatial patterns that, in addition, are described by several MEMs spanning wide ranges of scales. The number of MEMs describing species’ spatial variations (Figure 2) does not depend on species distribution, since no relationship (for α = 0.05) between the number of parks in which species were observed and number of MEMs describing their variations was found. Species with random spatial patterns belong to a few genera: Brachystomella, Desoria, Lathriopyga, Lepidocyrtus, Sminthurinus, Sminthurus, Willemia, as well as Folsomia and Isotomurus. In the latter two cases, however, the genus comprises also species with spatial variations described by several MEMs. The few species that occur in both Naples and Montpellier tend to have spatial variations described by similar scales (albeit with the limitations in comparing scales derived from different SWMs) or to exhibit random spatial variations in both locations. Fourth corner analysis indicates that none of the environmental drivers was significant (for α = 0.05) in explaining the spatial variations of species and CWMs in Naples and Montpellier.

4. Discussion

Overall, the patterns observed in the ranking of SWMs and in the choice of MEMs from the consensus SWMs, as well as the lack of relationships between MEMs and potential environmental drivers, provide clear insights into the spatial ecology of Collembola and of their functional traits. In terms of spatial patterns, species abundances and CWMs exhibit mostly small-scale variations, with limited spatial connectivity. Indeed, on the one hand, the degree of connectedness has little effect on the choice of SWMs in both S- and T-models, as indicated by the reduced differentiation among the mean ranks of SWMs based on different topologies. On the other hand, the generally high η of the concave weighting functions indicates that the abundance of species or the values of the CWMs change abruptly in space, actually removing (weighting down to zero) edges even in highly connected topologies. Taken together, these results suggest that collembolan communities could exhibit either random spatial structures or spatial structures at scales lower than the ones investigated in the present study, with patterns that disappear at larger spatial scales.
These deductions are in line with the fourth corner analysis, indicating that none among elevation, land use, park age, urban density, presence of litter, brushwood, or canopy cover, soil origin, C/N ratio, pH or bulk density explain the spatial variations of species and CWMs. Indeed, such an occurrence would be a direct consequence of a mismatch between the spatial scales at which the environmental parameters were measured and the scales of their variations that are relevant in determining Collembola spatial patterns. In other words, the spatial variations in environmental characteristics may structure collembolan communities at scales coherent with the perceived spatial scales of most of the species, but not at larger ones [36,37]. Such a hypothesis is confirmed by several studies on the drivers of diversity in edaphic communities, indicating either the presence of stochastic spatial distribution in the absence of environmental gradients [38] or micro-scale spatial structures shaped by both environmental filtering and species interactions [39,40,41,42].
An additional explanation for the independence of spatial patterns from environmental drivers refers, instead, to the effects of adaptation and amplitude of species’ ecological niches, which may actually decouple their distribution from environmental constraints. In this context, with wide variations in environmental parameters within and between parks/cities [13,14], the distribution of species can be considered a proxy of their ecological niche amplitude. A notable example of widely distributed species among the investigated parks is Cryptopygus thermophilus, which, unsurprisingly, colonizes a wide variety of environments worldwide [43,44]. Other species observed in both Naples and Montpellier are Isotomiella minor, Ceratophysella gibbosa, Protaphorura aurantiaca, Parisotoma notabilis, Proisotoma minuta, Lepidocyrtus lanuginosus, Pseudosinella immaculata and Folsomides parvulus. Interestingly, in spite of their wide distribution, these species encompass taxa with different patterns of spatial scales, from species characterized by both large- and small-scale variations, such as C. thermophilus or P. minuta, to species without clear spatial structure, such as L. lanuginosus. On the one hand, the lack of correlation between the number of significant MEMs and the number of parks where species were observed guarantees MEMs’ independence from species spatial distribution, and thus their ability to unbiasedly model species’s spatial structure. On the other hand, it is remarkable that species tend to show similar spatial structure, or the lack thereof, in both Naples and Montpellier, with patterns that are congruent with their taxonomy. For example, and considering all the observed species, those belonging to the Isotomidae family, with few exceptions, are all spatially structured, whereas those belonging to other families, such as Sminthuridae and Entomobryidae, are not. The appearance of spatial structure and the scales at which it appears seem thus to be species-specific rather than environment-driven characteristics [45]. In other words, the spatial patterns of Collembola, in terms of the scales describing their spatial variations, appear to be primarily determined by intrinsic factors, such as functional traits affecting the perceived spatial scale [36] and the dispersion or adaptability to different environments [37]. Species’ adoption of different dispersion strategies, such as wind transport or active locomotion [46], also supports this hypothesis, since species morphology or behavior may directly affect the efficiency of different dispersion processes. As a corollary for this explanation, the environment plays a minor role, likely contributing to shape species distribution at micro-scale only [36,40,42]—a conclusion in agreement with the deductions derived from the fourth corner analysis.
Overall, Collembola spatial patterns appear thus at micro-scale, where they are likely affected by environmental constraints, with scales of variation defined by species characteristics. In this respect, we supposed that CWMs would exhibit spatial patterns with higher connectivity and variations at larger scales than species, as for the expected behavior of functional diversity vs. taxonomic diversity according to the community convergence hypothesis [16]. However, the deduced scenario for species’ spatial variations holds true for the spatial patterns and scales of the functional traits, with mostly abrupt changes in CWM values from site to site and scales independent from the analyzed environmental parameters. In terms of scales, the functional traits investigated vary, with few exceptions, not only at large, but also at local scales, an occurrence that may be explained by the peculiar environment where the collembolan communities were studied. Indeed, urban parks are intrinsically more spatially heterogeneous than natural areas, with dramatic changes in the environment even at very fine scales, which may reflect the changes of not only the species composition of communities, but also their functional traits. One notable exception corroborating such a hypothesis is the pigmentation, which varies at large to medium scales, as expected from the scales of variation in the canopy cover or other means of shading modifying the light environment to which species are subjected, that affect this functional trait [47].

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/su132313029/s1, Table S1: Characteristics of parks from Naples (N) and Montpellier (M). The coding for park names, soil origin (N: native, T: technosol), previous land use (A: agricultural land, Ab: abandoned area, F: forest, L: leisure estate, Tq: tuff quarry, U: urban) as well as the abbreviations for dissolved organic carbon (DOC), urban density (UD) and bulk density (BD) follow the ones adopted in Milano et al. [13,14]. The presence (P) or absence (A) of litter, canopy cover and brushwood is also reported, as well as the soil carbon to nitrogen ratio (C/N).

Author Contributions

Conceptualization, A.B., D.B., J.C. and G.M.; methodology, A.B.; formal analysis, A.B. and D.B.; investigation, V.M., J.C. and G.M.; resources, J.C. and G.M.; data curation, A.B. and V.M.; writing—original draft preparation, A.B. and D.B.; writing—review and editing, A.B., D.B., V.M., L.S., J.C., G.M.; funding acquisition, J.C. and G.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding, but was carried out thanks to the utilities provided by the University Paul Valéry in Montpellier (France) and the University of Naples Federico II (Italy).

Data Availability Statement

The data collected and analyzed in this research are available upon request from the authors.

Acknowledgments

The authors would like to thank Bruno Buatois and Raphaelle Leclerc of the PACE platform (Montpellier, France) for facilitating physico-chemical analyses, Michel Bonnet and Sthéphanie Grosset of Montpellier municipality for providing cadastral data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Vergnes, A.; Le Viol, I.; Clergeau, P. Green corridors in urban landscapes affect the arthropod communities of domestic gardens. Biol. Conserv. 2012, 145, 171–178. [Google Scholar] [CrossRef] [Green Version]
  2. Morel, J.L.; Chenu, C.; Lorenz, K. Ecosystem services provided by soils of urban, industrial, traffic, mining, and military areas (SUITMAs). J. Soils Sediments 2015, 15, 1659–1666. [Google Scholar] [CrossRef]
  3. Guilland, C.; Maron, P.A.; Damas, O.; Ranjard, L. Biodiversity of urban soils for sustainable cities. Environ. Chem. Lett. 2018, 16, 1267–1282. [Google Scholar] [CrossRef]
  4. Medeiros-Sousa, A.R.; Fernandes, A.; Ceretti-Junior, W.; Wilke, A.B.B.; Marrelli, M.T. Mosquitoes in urban green spaces: Using an island biogeographic approach to identify drivers of species richness and composition. Sci. Rep. 2017, 7, 17826. [Google Scholar] [CrossRef] [Green Version]
  5. Fattorini, S.; Mantoni, C.; De Simoni, L.; Galassi, D. Island biogeography of insect conservation in urban green spaces. Environ. Conserv. 2018, 45, 1–10. [Google Scholar] [CrossRef]
  6. Mardiastuti, A. Urban park design for bird diversity: Theory and application in landscape and site scales. IOP Conf. Ser. Earth Environ. Sci. 2020, 501, 012026. [Google Scholar] [CrossRef]
  7. Bräuniger, C.; Knapp, S.; Kühn, I.; Klotz, S. Testing taxonomic and landscape surrogates for biodiversity in an urban setting. Landsc. Urban Plan. 2010, 97, 283–295. [Google Scholar] [CrossRef]
  8. Huang, Y.; Zhao, Y.; Li, S.; von Gadow, K. The effects of habitat area, vegetation structure and insect richness on breeding bird populations in Beijing urban parks. Urban For. Urban Green. 2015, 14, 1027–1039. [Google Scholar] [CrossRef]
  9. Threlfall, C.G.; Williams, N.S.G.; Hahs, A.K.; Livesley, S.J. Approaches to urban vegetation management and the impacts on urban bird and bat assemblages. Landsc. Urban Plan. 2016, 153, 28–39. [Google Scholar] [CrossRef]
  10. Vergnes, A.; Pellissier, V.; Lemperiere, G.; Rollard, C.; Clergeau, P. Urban densification causes the decline of ground-dwelling arthropods. Biodivers. Conserv. 2014, 23, 1859–1877. [Google Scholar] [CrossRef]
  11. Amossé, J.; Dózsa-Farkas, K.; Boros, G.; Rochat, G.; Sandoz, G.; Fournier, B.; Mitchell, E.A.D.; Le Bayon, R.-C. Patterns of earthworm, enchytraeid and nematode diversity and community structure in urban soils of different ages. Eur. J. Soil Biol. 2016, 73, 46–58. [Google Scholar] [CrossRef] [Green Version]
  12. Francini, G.; Hui, N.; Jumpponen, A.; Kotze, D.J.; Romantschuk, M.; Allen, J.A.; Setälä, H. Soil biota in boreal urban greenspace: Responses to plant type and age. Soil Biol. Biochem. 2018, 118, 145–155. [Google Scholar] [CrossRef] [Green Version]
  13. Milano, V.; Cortet, J.; Baldantoni, D.; Bellino, A.; Dubs, F.; Nahmani, J.; Strumia, S.; Maisto, G. Collembolan biodiversity in Mediterranean urban parks: Impact of history, urbanization, management and soil characteristics. Appl. Soil Ecol. 2017, 119, 428–437. [Google Scholar] [CrossRef]
  14. Milano, V.; Maisto, G.; Baldantoni, D.; Bellino, A.; Bernard, C.; Croce, A.; Dubs, F.; Strumia, S.; Cortet, J. Impact of the neighboring cityscape on Collembola biodiversity: A Mediterranean case study. Landsc. Urban Plan. 2018, 180, 135–147. [Google Scholar] [CrossRef]
  15. Li, X.; Zhu, H.; Geisen, S.; Bellard, C.; Hu, F.; Li, H.; Chen, X.; Liu, M. Agriculture erases climate constraints on soil nematode communities across large spatial scales. Glob. Change Biol. 2020, 26, 919–930. [Google Scholar] [CrossRef]
  16. Ricklefs, R.E. Evolutionary diversification and the origin of the diversity-environment relationship. Ecology 2006, 87, S3–S13. [Google Scholar] [CrossRef]
  17. IUSS Working Group WRB. World Reference Base for Soil Resources 2014, Update 2015; International Soil Classification System for Naming Soils and Creating Legends for Soil Maps; World Soil Resources Reports No. 106; FAO: Rome, Italy, 2015. [Google Scholar]
  18. ISO 23611-2. Soil Quality. Sampling of Soil Invertebrates. Part 2. Sampling and Extraction of Micro-Arthropods (Collembola and Acarina); International Organization for Standardization: Geneva, Switzerland, 2006. [Google Scholar]
  19. Macfadyen, A. Improved funnel-type extractors for soil arthropods. J. Anim. Ecol. 1961, 30, 171–184. [Google Scholar] [CrossRef]
  20. Pey, B.; Nahmani, J.; Auclerc, A.; Capowiez, Y.; Cluzeau, D.; Cortet, J.; Decaëns, T.; Deharveng, L.; Dubs, F.; Joimel, S.; et al. Current use of and future needs for soil invertebrate functional traits in community ecology. Basic Appl. Ecol. 2014, 15, 194–206. [Google Scholar] [CrossRef] [Green Version]
  21. Bretfeld, G. Synopses on Palaeartic Collembola Pt 2: Symphpleona; Staatliches Museum für Naturkunde: Görlitz, Germany, 1999. [Google Scholar]
  22. Dunger, W.; Schlitt, B. Synopses on Palaeartic Collembola Pt 6/1: Tullbergiidae; Senckenberg Museum of Natural History: Görlitz, Germany, 2011. [Google Scholar]
  23. Hopkin, S.P. Biology of the Springtails (Insecta: Collembola); Oxford University Press: Oxford, UK, 1997. [Google Scholar]
  24. Jordana, R. Synopses on Palaeartic Collembola Pt 7/1: Capryinae & Entomobryini; Senckenberg Museum of Natural History: Görlitz, Germany, 2012. [Google Scholar]
  25. Potapov, M. Synopses on Palaearctic Colembola Pt 3: Isotomidae; Museum of Natural History: Görlitz, Germany, 2001. [Google Scholar]
  26. Thibaud, J.M.; Schulz, H.J.; da Gama Assalino, M.M. Synopses on Palaeartic Collembola Pt 4: Hypogastruridae; Museum of Natural History: Görlitz, Germany, 2004. [Google Scholar]
  27. Zimdars, B.; Dunger, W. Synopses on Palaeartic Collembola Pt 1: Tullbergiinae; Museum of Natural History: Görlitz, Germany, 1994. [Google Scholar]
  28. Gisin, H. Okologie und Lebengemeinschaften der Collembolen im Schweizerischen Exkursionsgebiet Basels. Rev. Suisse Zool. 1943, 50, 131–224. [Google Scholar]
  29. Joimel, S.; Schwartz, C.; Bonfanti, J.; Hedde, M.; Krogh, P.H.; Pérès, G.; Pernin, C.; Rakoto, A.; Salmon, S.; Santorufo, L.; et al. Functional and taxonomic diversity of Collembola as complementary tools to assess land use effects on soils biodiversity. Front. Ecol. Evol. 2021, 9, 450. [Google Scholar] [CrossRef]
  30. Hawes, T.C. The unusually colourful ‘clothing’ changes of Lepidophorella australis Carpenter, 1925 (Collembola: Tomoceridae). N. Z. Entomol. 2017, 40, 22–29. [Google Scholar] [CrossRef]
  31. Podani, J. Extending Gower’s general coefficient of similarity to ordinal characters. Taxon 1999, 48, 331. [Google Scholar] [CrossRef]
  32. Garnier, E.; Cortez, J.; Bill, S.G.; Navas, M.-L.; Roumet, C.; Debussche, M.; Laurent, G.R.; Blanchard, A.; Aubry, D.; Bellmann, A.; et al. Plant functional markers capture ecosystem properties during secondary succession. Ecology 2004, 85, 2630–2637. [Google Scholar] [CrossRef]
  33. Bellino, A.; Alfani, A.; De Riso, L.; Gregorio, R.; Pellegrino, T.; Baldantoni, D. Long-established and new active biomonitors jointly reveal potentially toxic element gradients across spatial scales in freshwater ecosystems. Ecol. Indic. 2020, 118, 106742. [Google Scholar] [CrossRef]
  34. R Development Core Team. A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2020. [Google Scholar]
  35. Dray, S.; Dufour, A.-B. The ade4 package: Implementing the duality diagram for ecologists. J. Stat. Softw. 2007, 22, 1–20. [Google Scholar] [CrossRef] [Green Version]
  36. Zinger, L.; Taberlet, P.; Schimann, H.; Bonin, A.; Boyer, F.; De Barba, M.; Gaucher, P.; Gielly, L.; Giguet-Covex, C.; Iribar, A.; et al. Body size determines soil community assembly in a tropical forest. Mol. Ecol. 2019, 283, 528–543. [Google Scholar] [CrossRef]
  37. Piano, E.; Souffreau, C.; Merckx, T.; Baardsen, L.F.; Backeljau, T.; Bonte, D.; Brans, K.I.; Cours, M.; Dahirel, M.; Debortoli, N.; et al. Urbanization drives cross-taxon declines in abundance and diversity at multiple spatial scales. Glob. Change Biol. 2020, 26, 1196–1211. [Google Scholar] [CrossRef]
  38. Bahram, M.; Kohout, P.; Anslan, S.; Harend, H.; Abarenkov, K.; Tedersoo, L. Stochastic distribution of small soil eukaryotes resulting from high dispersal and drift in a local environment. ISME J. 2016, 10, 885–896. [Google Scholar] [CrossRef] [Green Version]
  39. Coulson, S.J.; Hodkinson, I.D.; Webb, N.R. Microscale distribution patterns in high Arctic soil microarthropod communities: The influence of plant species within the vegetation mosaic. Ecography 2003, 26, 801–809. [Google Scholar] [CrossRef]
  40. Gao, M.; He, P.; Zhang, X.; Liu, D.; Wu, D. Relative roles of spatial factors, environmental filtering and biotic interactions in fine-scale structuring of a soil mite community. Soil Biol. Biochem. 2014, 79, 68–77. [Google Scholar] [CrossRef]
  41. Gao, M.; Liu, D.; Lin, L.; Wu, D. The small-scale structure of a soil mite metacommunity. Eur. J. Soil Biol. 2016, 74, 69–75. [Google Scholar] [CrossRef]
  42. Arribas, P.; Andújar, C.; Salces-Castellano, A.; Emerson, B.C.; Vogler, A.P. The limited spatial scale of dispersal in soil arthropods revealed with whole-community haplotype-level metabarcoding. Mol. Ecol. 2021, 30, 48–61. [Google Scholar] [CrossRef] [PubMed]
  43. Hopkin, S.P. A Key to the Collembola (Springtails) of Britain and Ireland; Field Studies Council: Telford, UK, 2007. [Google Scholar]
  44. Fjellberg, A. The Collembola of Fennoscandia and Denmark—Part II: Entomobryomorpha and Symphypleona; Fauna Entomologica Scandinavica; Brill: Leiden, The Netherlands, 2007; Volume 42. [Google Scholar]
  45. Collins, G.E.; Hogg, I.D.; Convey, P.; Barnes, A.D.; McDonald, I.R. Spatial and temporal scales matter when assessing the species and genetic diversity of springtails (Collembola) in Antarctica. Front. Ecol. Evol. 2019, 7. [Google Scholar] [CrossRef] [Green Version]
  46. Dunger, W.; Schulz, H.-J.; Zimdars, B. Colonization behaviour of Collembola under different conditions of dispersal. Pedobiologia 2002, 46, 316–327. [Google Scholar] [CrossRef] [Green Version]
  47. Salmon, S.; Ponge, J.F. Species traits and habitats in springtail communities: A regional scale study. Pedobiologia 2012, 55, 295–301. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Mean ranks ( ϑ ^ SWM) relative to each SWM for the S-models (upper panels) and the T-models (lower panels) in Naples (left panels) and Montpellier (right panels). Mean ranks are coded in different colors according to the connectivity topology (red: Delaunay triangulation, green: Gabriel’s graph, cyan: minimum spanning tree, violet: relative neighbour graph) and are ordered according to the weighting schemes, reported on the x-axis along with their exponential parameter (η). The SWM scoring the lowest mean rank for each model was adopted as the SWM describing the spatial pattern of Collembola communities (for the S-models) and CWMs (for the T-models).
Figure 1. Mean ranks ( ϑ ^ SWM) relative to each SWM for the S-models (upper panels) and the T-models (lower panels) in Naples (left panels) and Montpellier (right panels). Mean ranks are coded in different colors according to the connectivity topology (red: Delaunay triangulation, green: Gabriel’s graph, cyan: minimum spanning tree, violet: relative neighbour graph) and are ordered according to the weighting schemes, reported on the x-axis along with their exponential parameter (η). The SWM scoring the lowest mean rank for each model was adopted as the SWM describing the spatial pattern of Collembola communities (for the S-models) and CWMs (for the T-models).
Sustainability 13 13029 g001
Figure 2. Heatmaps of the LASSO coefficients associated to the MEMs deriving from the consensus SWMs for the S-models (upper panels) and the T-models (lower panels) in Naples (left panels) and Montpellier (right panels). MEMs are ordered from the largest to the finest scale. The number of colored squares for each column indicates the number of MEMs describing species or trait spatial variations. Labels on the x-axes are coded according to Table 1 and Table 2.
Figure 2. Heatmaps of the LASSO coefficients associated to the MEMs deriving from the consensus SWMs for the S-models (upper panels) and the T-models (lower panels) in Naples (left panels) and Montpellier (right panels). MEMs are ordered from the largest to the finest scale. The number of colored squares for each column indicates the number of MEMs describing species or trait spatial variations. Labels on the x-axes are coded according to Table 1 and Table 2.
Sustainability 13 13029 g002
Table 1. Species of Collembola observed in urban parks of Naples (N), Montpellier (M) or both the cities (B), grouped according to the family they belong to and with indication of the identifying label.
Table 1. Species of Collembola observed in urban parks of Naples (N), Montpellier (M) or both the cities (B), grouped according to the family they belong to and with indication of the identifying label.
Isotomidae Entomobryidae
Cryptopygus thermophilus AxelsonC_theBLepidocyrtus cyaneus TullbergL_cyaM
Desoria tigrina NicoletD_tigNLepidocyrtus lanuginosus GmelinL_lanB
Folsomia penicula BagnallF_penNLepidocyrtus lignorum FabriciusL_ligN
Folsomia quadrioculata TullbergF_quaNPseudosinella alba PackardP_albM
Folsomia sexoculata TullbergF_sexNPseudosinella immaculata Lie-PettersenP_immB
Folsomia similis BagnallF_simM
Folsomides parvulus StachF_parB
Isotomiella minor SchaefferI_minB
Isotomodes templetoni BagnallI_temN
Isotomurus maculatus SchaefferI_macN
Parisotoma notabilis SchaefferP_notB
Proctostephanus stuckeni BornerP_stuM
Proisotoma minuta TullbergP_minB
Hypogastruridae Tullbergiidae
Ceratophysella gibbosa BagnallC_gibBParatullbergia callipygos BornerP_calM
Schoettella ununguiculata TullbergS_unuMParatullbergia caroli Luciañez, Ruiz & SimonP_carN
Willemia denisi MillsW_denMParatullbergia macdougalli BagnallP_macM
Sminthuridae Onychiuridae
Sminthurinus elegans FitchS_eleNProtaphorura armata TullbergP_armN
Sminthurus viridis LinnaeusS_virNProtaphorura aurantiaca RidleyP_aurB
Neanuridae Brachystomellidae
Lathriopyga longiseta CaroliL_lonNBrachystomella parvula SchaefferB_parN
Table 2. Functional traits and the relative coding adopted in the present research, together with the rationale behind trait choice.
Table 2. Functional traits and the relative coding adopted in the present research, together with the rationale behind trait choice.
TraitCodingRationale
Body lengthShort (<1.5 mm)Food resource use
Medium (1.5–3.0 mm)
Long (>3.0 mm)
Body shapeSphericalHabitat preference
Cylindrical
FurcaAbsent/VestigialMovement/dispersion
Short
Long
PigmentationPresentHabitat preference
Absent
EyesPresentHabitat preference
Absent
PAOPresentHabitat preference
Absent
ScalesPresentPredation avoidance
Absent
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bellino, A.; Baldantoni, D.; Milano, V.; Santorufo, L.; Cortet, J.; Maisto, G. Spatial Patterns and Scales of Collembola Taxonomic and Functional Diversity in Urban Parks. Sustainability 2021, 13, 13029. https://0-doi-org.brum.beds.ac.uk/10.3390/su132313029

AMA Style

Bellino A, Baldantoni D, Milano V, Santorufo L, Cortet J, Maisto G. Spatial Patterns and Scales of Collembola Taxonomic and Functional Diversity in Urban Parks. Sustainability. 2021; 13(23):13029. https://0-doi-org.brum.beds.ac.uk/10.3390/su132313029

Chicago/Turabian Style

Bellino, Alessandro, Daniela Baldantoni, Vittoria Milano, Lucia Santorufo, Jérôme Cortet, and Giulia Maisto. 2021. "Spatial Patterns and Scales of Collembola Taxonomic and Functional Diversity in Urban Parks" Sustainability 13, no. 23: 13029. https://0-doi-org.brum.beds.ac.uk/10.3390/su132313029

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