Skip to main content
  • Research article
  • Open access
  • Published:

Narrow environmental niches predict land-use responses and vulnerability of land snail assemblages

Abstract

Background

How land use shapes biodiversity and functional trait composition of animal communities is an important question and frequently addressed. Land-use intensification is associated with changes in abiotic and biotic conditions including environmental homogenization and may act as an environmental filter to shape the composition of species communities. Here, we investigated the responses of land snail assemblages to land-use intensity and abiotic soil conditions (pH, soil moisture), and analyzed their trait composition (shell size, number of offspring, light preference, humidity preference, inundation tolerance, and drought resistance). We characterized the species’ responses to land use to identify ‘winners’ (species that were more common on sites with high land-use intensity than expected) or ‘losers’ of land-use intensity (more common on plots with low land-use intensity) and their niche breadth. As a proxy for the environmental ‘niche breadth’ of each snail species, based on the conditions of the sites in which it occurred, we defined a 5-dimensional niche hypervolume. We then tested whether land-use responses and niches contribute to the species’ potential vulnerability suggested by the Red List status.

Results

Our results confirmed that the trait composition of snail communities was significantly altered by land-use intensity and abiotic conditions in both forests and grasslands. While only 4% of the species that occurred in forests were significant losers of intensive forest management, the proportion of losers in grasslands was much higher (21%). However, the species’ response to land-use intensity and soil conditions was largely independent of specific traits and the species’ Red List status (vulnerability). Instead, vulnerability was only mirrored in the species’ rarity and its niche hypervolume: threatened species were characterized by low occurrence in forests and low occurrence and abundance in grasslands and by a narrow niche quantified by land-use components and abiotic factors.

Conclusion

Land use and environmental responses of land snails were poorly predicted by specific traits or the species’ vulnerability, suggesting that it is important to consider complementary risks and multiple niche dimensions.

Background

Land use disturbs natural environments, changes local geographical landscape structure and alters local biotic and abiotic conditions, e.g. microclimate [1,2,3,4,5,6]. Reduction of habitat and microhabitat heterogeneity may lead to a homogenization of plant and animal communities, trigger a reduction in functional diversity and thus lower the capacity of an ecosystem to buffer disturbances [7, 8]. Homogenization of animal communities by increasing land-use intensity has been shown for several taxa; e.g., in managed grasslands, 34% of plant- and leafhoppers species were significant losers (i.e. species that were significantly less abundant under conditions of high land-use intensity) of land-use intensification, particularly increases in mowing frequency had a negative effect [9].

Land snails are an important macroinvertebrate group that is directly and indirectly involved in ecosystem processes such as litter decomposition or nutrient cycling [10, 11]. There is a natural north–south and west–east gradient of snail species distributions and abundances within Europe; species richness increases from north to south and to a lesser extent from west to east which is linked to regional and ecological differences and the land-use history [12]. Snail species also differ in their tolerance to abiotic factors (pH, soil moisture), and vary greatly in life-history parameters (e.g., lifespan, development, number of offspring, food requirement, shell size) and general behavior [13] which also affect their distribution. Variation in body size and diet seems to be especially important for structuring snail communities [14] as well as the species-specific tolerance to a variety of environmental factors which can result in nested communities at a specific site [15, 16].

Studies on trait composition of snail communities in Sweden pointed to the importance of the species’ niche-width and the importance of local environmental conditions over spatial variables [17]. While tolerance-related traits such as humidity preference or inundation tolerance were positively associated with abiotic soil moisture, a large amount of variation remained unexplained [17], which may be related to land use. The impact of land use and its intensity on land snail communities is less intensively investigated although most land snail species are characterized by a limited mobility and therefore are vulnerable to human introduced habitat changes [15, 18,19,20]. Changes in abiotic factors such as soil pH, soil moisture, soil calcium content, leaf litter depth, soil surface structure or the type of vegetation have been shown to alter snail communities [15, 21,22,23,24,25]. Also land-use factors such as the proportion of wood harvested in forests or the amount of grazing livestock in grasslands can influence snail communities directly and/or indirectly [20, 26, 27]. In addition, disturbances by different land-use types and intensities may alter the trait composition of snail communities on the regional level; i.e. the presence of coniferous timber may favor snail communities with differing traits than communities in natural deciduous stands.

In the present study, we investigated land snail communities at forest and grassland sites in different regions of Germany, which were characterized by different land-use types and intensities. We aimed to test whether the trait composition of the snail community is influenced by land-use intensity (and soil conditions). We then tested the responses of each snail species to land-use intensity; ‘winners’ significantly increase in abundance and occurrence with land-use intensity, whereas ‘losers’ significantly decrease compared to the null model [9, 39]. We than compared these responses with the snail species’ habitat association; i.e. we asked whether species that only occasionally occur in forests are more affected by forest management than species that are specialized to forest habitats. On the other hand, do species that are grassland specialist suffer less from grassland management than those only occasionally occurring in grasslands? Finally, we compared our findings of the land-use effects and the ‘winner/loser’ status of a species with its putative vulnerability (Red List status), to test if losers of land-use intensifications in forests and grasslands are those species that are classified as vulnerable.

Results

Response to land use

The trait composition of land snail communities differed strongly between forests and grasslands within regions, indicated by a strong differentiation of community-weighted mean trait values (CWMs). Assemblages of forest species consisted of larger species, consistently showed lower light and higher humidity preference, lower drought resistance and mostly lower inundation tolerance than grassland assemblages; differences in the number of offspring were inconsistent among forest and grassland habitats (Fig. 1).

Fig. 1
figure 1

Trait distribution (a shell size, b number of offspring, c light preference, d humidity preference, e drought resistance, f inundation tolerance) of snail communities among forest (grey) and grassland (white) habitats in the Swabian Alb, the Hainich-Dün and the Schorfheide-Chorin. Traits are given as community weighted mean (CWM). Difference among habitats per region are tested using an ANOVA (asterisks), differences between regions are tested by a posthoc Tukey test (letters). Significances: ns not significant, *p < 0.05, **p < 0.01, ***p < 0.001

In forests, land-use intensity and abiotic conditions significantly influenced the CWMs of all traits investigated, although often in a different way across regions (Table 1, Additional file 1: Appendix 1; see interaction terms with region). Similarly, in grasslands the trait composition of snail communities was significantly influenced by most land-use components and abiotic conditions (Table 2, Additional file 1: Appendix 1).

Table 1 Influence of land-use parameter and abiotic factors on the trait composition of snail communities in forest habitats
Table 2 Influence of land-use parameter and abiotic factors on the trait composition of snail communities in grassland habitats

In forest habitats, some 4% of all species were ‘losers’ of the combined forest management index (i.e. they were significantly less common in intensively used forests), whereas 12% were ‘winners’ and thus increased with forest management intensity (Table 3). The proportions of non-native trees (4% losers vs. 8% winners) and the proportion of dead wood with saw cuts (6% losers vs. 8% winners) revealed a similar pattern, but for the proportion of wood harvested the percentage of losers (12%) exceeded that of winners (8%).

Table 3 Red list status, occurrence and total abundance of snail species in the Swabian Alb (A), the Hainich-Dün (H) and the Schorfheide-Chorin (S) in forest habitats

In grasslands, many species were predominantly found at low land-use intensities (LUI); 21% of all species were significant losers and only Monacha cartusiana profited from high LUI (Table 4). However, single land-use components in grasslands had only weak effects. Grazing intensity positively affected Cecilioides acicula and Cepaea hortensis, but showed no negative impact. Similarly, mowing (2% losers and 2% winners) and fertilization (4% losers and 4% winners) had a very little impact compared to the combined LUI.

Table 4 Red list status, occurrence and total abundance of snail species in the Swabian Alb (A), the Hainich-Dün (H) and the Schorfheide-Chorin (S) in grassland habitats

However, in both forests and grasslands, species’ land-use responses (i.e. their ‘winner/loser’ status) were independent of their traits; i.e. losers in forests or grasslands were neither characterized by a smaller or larger shell size nor by lower or higher numbers of offspring nor by lower or higher light preference etc. (Additional files 215: Appendix 2–15).

Response to abiotic factors

Although niches of common land snail species for soil pH and soil moisture were generally broad, some differentiation was found in the communities of both habitats. In forests, Aegopinella pura, the genus Carychium, Cochlicopa lubrica, Ena montana and Vitrea contracta were significantly associated with higher pH values (Table 3) and Cepaea hortensis, Euconulus fulvus, Nesovitrea hammonis, Vallonia pulchella and Vitrinobrachium breve were found at sites with low pH (Table 3). Furthermore, A. pura and Carychium tridentatum were associated with high soil moisture in forests and Ceciliodes acicula, E. fulvus, N. hammonis, Punctum pygmaeum, Trochulus striolatus and V. pulchella were found at low soil moisture values (Table 3).

Grassland sites had a higher mean pH (6.7) as compared to forest soils, and many snail species (e.g., Candidula unifasciata, the genus Carychium, Granaria frumentum, Pupilla muscorum, Vertigo antivertigo) were associated with higher pH values (Table 4). Only N. hammonis was significantly more common on sites with low pH. Soil moisture niches of grassland species were even broader than those of pH. The genus Carychium, Trochulus hispidus and Vallonia pulchella were found at high moisture values, while C. unifasciata, Discus rotundatus, Truncatellina cylindrica, V. excentrica were associated with low soil moisture (Table 4).

Habitat association

Snail species differed in their habitat association and their distribution among regions (Fig. 2). However, effects of land-use management components and abiotic factors in forests were independent of the species’ habitat association, i.e. species that occurred in forests at low frequencies (e.g., 25% of the individuals in Cochlicopa lubrica; Fig. 2) were equally affected by land-use intensification as species that are exclusively found in forests (e.g., Cepaea hortensis) (F1,49 = 0.14, p = 0.71, Fig. 2, Additional file 14: Appendix 14). In contrast, grassland species that predominately prefer grassland habitats were less tolerant to fertilization than species that also occur in forests (F1,50 = 5.84, p = 0.019, Fig. 3a, Additional file 15: Appendix 15). Furthermore, grassland “specialists” were significantly associated with higher pH values (F1,49 = 9.21, p = 0.004, Fig. 3b).

Fig. 2
figure 2

Relation between the responses (abundance-weighted mean) of each snail species to fertilization (a) and soil pH (b) and their proportional occurrence in forests. Indicated species above the line are significant “winners” for fertilization respective soil pH, indicated species below the line (in italic) are significant “losers”

Fig. 3
figure 3

Proportional distribution of land snail species in the Schorfheide-Chorin, the Hainich-Dün and the Swabian Alb. Grasslands are given in light grey, forests in dark grey. The three most abundant species are symbolized by big circles, less abundant species by small circles. Species that are underlined are specific for the respective region. Percentages in brackets indicate the proportional occurrence of species of the same genus

Species’ vulnerability

Across forests and grasslands, 75% of the 61 snail species found are currently not threatened or endangered according to their Red List status (Tables 3, 4). Nevertheless, Nesovitrea petronella, Candidula unifasciata and Granaria frumentum are regarded as ‘endangered’ while Vallonia enniensis is ‘highly endangered’ and V. angustior is listed on the FFH directive.

There was no statistical support that a negative response to land-use intensity of a certain species (“loser”) is associated with a high vulnerability of the species, neither in forests nor in grasslands (Table 5). A better predictor for the species’ vulnerability in forests was a relatively low number of sites in which the species occurred, and in grasslands both a low occurrence and a low total abundance corresponded to a higher vulnerability (Table 5). Furthermore, the 5-dimensional niche hypervolume based on the species’ tolerance to land-use components and abiotic conditions was significantly correlated with the species’ vulnerability, hence species with a small niche hypervolume are more vulnerable in both forests (Spearman rank test: S = 20,091, p = 0.0004; Fig. 4a) and grasslands (Spearman rank test: S = 15,547, p = 0.003, Fig. 4b).

Table 5 Statistical p values of a general linearized model with Poisson distribution testing the influence of land-use parameters and abiotic factors on species vulnerability
Fig. 4
figure 4

Species vulnerability in relation to the five-dimensional niche hypervolume in forest (a) and grassland (b). The hypervolume was the product of the abundance-weighted standard deviations (AWSDs) of all single land-use components as well as pH and soil moisture in forests or grasslands, respectively

Discussion

Response to land use and abiotic factors

Land snail species are slow-dispersing organisms, and historical influences are of general importance for their distribution [28]. Their diversity and heterogeneity is modified by predation, parasitism, competition, abiotic environmental gradients, natural barriers and disturbances [16]. While abiotic and vegetation parameters can be used to predict snail communities, disturbances by human land use are less frequently discussed. Our previous study [27] focused on land snail density, diversity and species composition and emphasized that direct impacts of land use on snail communities were on average lower than the impact of abiotic drivers and biotic substrates. However, unlike several studies on insects, few direct effects have been shown for wood harvesting in forests and mowing in grasslands on snail diversity [27]. How these direct land-use effects influence populations of single species and whether these effects are related to species-specific traits remains largely unclear.

Our study showed that snail assemblages varied consistently in their trait composition (shell size, number of offspring, light and humidity preference, drought resistance and inundation tolerance) across regions and among the two habitats, forests and grasslands. The variation between regions is consistent with a biogeographic gradient of increasing land snail diversity from the north to south caused by historical and ecological factors (temperature, moisture) [12, 22] and snail species responded differently to variable physical environments [13]. Local environmental conditions have been shown to explain about 19% of the trait variability of a snail metacommunity in Sweden [17], where the authors suggested that the unexplained variation may mirror land use. Our results confirmed that land-use intensity significantly influenced the trait distribution of snail communities, a pattern that was more pronounced in forest habitats than in grasslands. Since snail species in forest communities seem to be more specialized than those of grassland communities [12, 28], they may suffer more from habitat changes. For example, as the activity level of snails is temperature-dependent, thinning the canopy by wood harvesting or a high amount of non-native trees can enhance solar irradiance and the enhanced snail locomotion allows the exploitation of ambient heterogeneity [29] and may favor species with higher light preferences. This hypothesis is consistent with results from snail assemblages in our study, since the community-weighted mean (CWM) of light preference increased with the amount of non-native, mainly coniferous trees that may not have a closed canopy. Furthermore, changes of the community trait composition are not only directly caused by land-use parameters, but also by indirectly changing abiotic factors such as soil pH and soil moisture although most snail species exhibit broad niches for these abiotic factors.

In our study, 4% of all forest and 21% of all grassland snail species were significant losers concerning the compound indices of land-use intensity, including three land-use components in the forests or in the grasslands, respectively. The proportion of losers among grassland snail species was lower than the level found for grasshoppers (about 52%) [30] and plant- and leafhoppers (about 34%) [9], but similar to that for moths (28%) [31], confirming that snails are a suitable indicator for habitat quality and land-use intensity [17, 22, 32, 33]. The low proportion of loser species may be explained by their ground-living behavior (intangible for combine harvesters), the presence of a shell (protection against exposure and predation) and a larger diet breath compared to insect taxa (omnivory for flexibly changing food resources). However, we may have underestimated the amount of loser species since we did not distinguish between living individuals and empty shells. Empty shells decay at different rates under different ecological conditions [44]. Therefore, in some cases we may have evaluated shells of species which can no longer be found alive in the respective places. Keeping this in mind, our methodological approach may have ramifications on the conclusions drawn.

While increasing land-use intensity in open habitats is known to trigger a decline of pollinator species, and such losses were associated with species-specific trait attributes such as a narrow diet breadth, climate specialization, a large body-size and low fecundity [33,34,35,36,37,38,39], we did not find traits for snail species to correspond with their land-use response at species level. This is surprising, given that particularly those traits that are associated with soil moisture (drought resistance, inundation tolerance), body size or reproductive outcome are likely to respond to human-mediated disturbances. Furthermore, land-use effects in forests were independent of the species habitat association (i.e. forests specialists were equally affected as non-forests specialists), but grassland specialists suffered more from land use (i.e. fertilization) and were more dependent on high soil pH.

Note that single land-use parameters and abiotic conditions are often confounded in real landscapes as in our study, and thus responses of some snail species may not always correspond to single environmental dimensions as known from their global distribution or other sources. For example, Cochlicopa lubricella is a xerophilic land snail [42] whereas our data showed a neutral response to soil moisture.

Species’ vulnerability

The range of resources and the ecological conditions generally define the niche breadth and determine the geographical area of a species at the small or large scale [40]. Specialists are expected to be more vulnerable to habitat loss and climate change due to synergistic effects of a narrow niche and a small range size.

Only a few snails in our study across managed forests and grasslands are considered threatened or endangered according to the national Red List. Consistent with the expectation based on their environmental niche breadth, the species’ vulnerability status was significantly predicted by a particularly narrow niche hypervolume—an index that includes single land-use components as well as pH and soil moisture in each habitat. The smaller the hypervolume of a species, the higher its vulnerability according to the Red List. In addition, rarity was important: in forests, the most important predictor for their vulnerable status was a low number of sites in which they occurred. In grasslands, both their restricted occurrence and low total abundance predicted the species’ vulnerability.

Conclusion

In summary, our results indicate that the trait composition of snail communities was significantly altered by land-use intensities and abiotic conditions, and several species especially in grasslands were losers of intensive land use. These land-use and environmental responses were largely independent of specific traits and the species’ Red List status—this suggests that complementary risks may be important for predicting a species’ vulnerability. Instead, species vulnerability was mirrored in the species’ rarity and its overall niche hypervolume including single land-use components and abiotic factors.

Methods

Data origin

Data for this study were already part of a previous analysis of biodiversity and community composition, i.e. Wehner et al. [27] and are available at https://www.bexis.uni-jena.de/PublicData/PublicDataSet.aspx?DatasetId=24986. Wehner et al. [27] collected 15,607 snail individuals belonging to 71 taxa in three regions in Germany in the framework of the Biodiversity Exploratories Project (http://www.biodiversity-exploratories.de) [2]. The collaborative research unit addresses effects of land-use on biodiversity and biodiversity-related ecosystem processes in three regions: the Swabian Alb (ALB), a low-mountain range in South-West Germany (460–860 m a.s.l., 09° 10′ 49″–09° 35′ 54″ E/48° 20′ 28″–48° 32′ 02″ N), the Hainich-Dün (HAI), a hilly region in Central Germany (285–550 m a.s.l., 10° 10′24″–10° 46′ 45″ E/50° 56′ 14″–51° 22′ 43″ N) and the Schorfheide-chorin (SCH), a glacial formed landscape in North-East Germany (3–140 m a.s.l., 13° 23′ 27″–14° 08′ 53″ E/52° 47′ 25″–53° 13′ 26″ N). SCH is characterized by the lowest annual precipitation (520–580 mm), with a mean annual temperature of 6–7 °C. It is followed by HAI (630–800 mm, 6.5–8 °C) and ALB (800–930 mm, 8–8.5 °C).

In each region, 100 experimental plots (50 in forests and 50 in grasslands) were setup in 2008 along a land-use gradient covering different management types and intensities including mowing, grazing and fertilization in grasslands and the proportion of non-native trees, the proportion of dead-wood with saw cuts and the proportion of wood harvested in forests (Table 6). Forest plots have a size of 1 ha and grassland plots are 0.5 ha in size.

Table 6 Description and origin of land-use parameter and abiotic factors

In June 2017, Wehner et al. [27] took five replicated surface samples from all 50 forest and 50 grassland experimental plots (EPs) in the Swabian Alb and the Hainich, and from 49 forest and 34 grassland plots in the Schorf-heide due to constrained accessibility (1415 samples in total). Shelled snails were subsequently determined to the species, genus or family level using [41,42,43]. Although suggested elsewhere [e.g., 44], [27] did not distinguish between empty shells and living snail individuals.

As our current study focuses on species-level responses, only those individuals that could be assigned to the species level were used (ALB grasslands: 36, ALB forests: 37, HAI grassland: 31, HAI forest: 35, SCH grassland: 24, SCH forest: 21, 61 different land snail species in total). Grassland plots (although not permanently flooded) in one region (Schorfheide) harbored large numbers of aquatic and semi-aquatic snails. In contrast to our previous analysis that covered all snails recorded [27], we excluded aquatic snails from the analyses since their role and responses to terrestrial environmental variables such as land-use in grasslands remain unclear,

Statistical analyses

All statistical analyses were performed in R 3.5.2 [45] using the main packages “car” [46], “dplyr” [47], “lme4” [48] and “SMDTools” [49].

Trait composition of snail communities

Morphological and life-history trait values for all snail species were obtained from an established trait database by Falkner et al. [50] and compared to findings of [51] whenever possible; see Astor et al. [17] for a similar approach based on [50]. Traits for the set of species in our study are summarized in Table 7. Note that these traits are either continuous variables (size), integers (offspring) or ranks (all others); ranks can been treated as integers or continuous variables for an analysis based on community weighted mean (CWM, see below); the resulting distribution of the CWM in species-rich communities and across a large number of plots typically approach a Gaussian distribution. Moreover, to explore the response to potential environmental filtering, traits with different meaning are treated independently for the following analysis (a common practice, although some traits, e.g. shell size and number of offspring, may be correlated, see [17]).

Table 7 Characterization of snail traits according to Falkner et al. 2001 [50]

For comparing snail communities among habitats and regions, the community weighted mean (CWM) of each trait was calculated as CWM per plot p

$${CWM}_{p}= \sum_{i=1}^{I}{T}_{i}\bullet \frac{{a}_{i,p}}{{A}_{p}}$$

where Ti is the trait value of species i, ai,p is the abundance of species i in plot p and Ap the total abundance of all snails in plot p (total I species).

Environmental niches

We characterized the environmental conditions of each forest or grassland plot by its land-use intensity and two abiotic soil parameters (pH and soil moisture; Table 6) [52, 53]. Data were obtained from the BExIS database (Table 6).

We tested the response of the CWM of each trait to variation in environmental conditions using linear regressions. Values for grazing and fertilization were square root transformed before statistical analyses.

In order to characterize the snail species’ responses to environmental conditions (land-use gradient, soil conditions), we calculated each species’ “environmental niche”. The method has been established in the context of the Biodiversity Exploratories and was applied to several taxa such as grasshoppers [30], cicadas, moths [31], bumblebees [54] or plants [55]. The “niche optimum” was calculated as the abundance weighted mean (AWM) for species i as

$${AWM}_{i}= \sum_{p=1}^{{n}_{p}}{L}_{p}\bullet \frac{{a}_{i,p}}{{A}_{i}}$$

where np is the number of plots investigated, Lp is the land-use gradient value of plot p, ai,p the abundance of species i in plot p and Ai the total abundance of species i across all 149 forest or 134 grasslands sites, respectively. Hence, the CWM characterizes the plots by the trait distribution of snails, and the AWM characterizes snail species by the environmental conditions of the plot, and the snail abundance ai,p is used to weight either species or plot, respectively.

In addition to the AWM as a niche optimum, we also characterized the “niche breadth” of each species to a single environmental variable using the abundance-weighted standard deviation (AWSD) [30]. To test whether AWMs and AWSDs statistically deviate from an expected random distribution, we compared the calculated values against the expected values obtained from a null model that distributes each species across Ni sites with the same probability, with Ni being the number of sites in which species i was found. The null model thus chooses values of the focal land-use parameter (LUI, Formi, single components, pH, soil moisture) of Ni sites and calculates a distribution of predicted AWMs and AWSDs values for each species based on 10,000 iterations. The null model was restricted to the one, two, or three regions in which the species was recorded to consider potential distribution boundaries of each species in Germany that may not be related to plot conditions [30].

As in any randomization model, the proportion of AWMs or AWSDs from 10,000 null models with greater or smaller expected values respectively than the observed value, provides the p value for the significance of the deviation between observed and expected values. A ‘winner’ is defined as a species with an observed AWM larger than the upper 5% of the distribution of AWMs obtained by the null models (i.e. adapted on higher-than average land-use intensity), a ‘loser’ shows an observed AWM smaller than the lower 5% (low land-use intensity specialist). For species which could be classified neither as ‘losers’ nor as ‘winners’, we tested whether they are specialized on intermediate land-use or abiotic levels, that is, whether they have an intermediate AWM with a narrower niche than expected. We standardized the niche breadth as weighted coefficient of variation (CV = AWSD/AWM) to account for the increase in SD with increasing mean, and compared observed CV and expected CV from the null models. This comparison allows us to distinguish ‘opportunists’ (observed CV ≥ expected CV) from species that are ‘specialized’ on intermediate land-use intensities (observed CV < expected CV and species not only occurring on one site, i.e., CV ≠ 0) [30]. The environmental niche (AWM, AWSD) and the assignment of low- and high-gradient specialists were also calculated for soil pH and soil moisture, although we did not adopt the ‘loser’/’winner’ terminology here unlike for land-use intensity.

Species vulnerability

Vulnerability (classified as a rank variable comparable to IUCN categories: least concern, endangered to unknown extent, very rare, near threatened, critically endangered, endangered, vulnerable) of land snail species was obtained from the Red List 2011 (according to [56]; see Table 3). We tested the relation of vulnerability with the species’ habitat association by calculating the proportional occurrence in either forest or grassland habitats of a certain species’ presence; a ‘specialist’ was defined if more than 90% of all individuals found were present in one habitat (forest or grassland). The relation between vulnerability and species’ habitat association was tested by a linear model using the land-use management components and the abiotic conditions as fixed factors and the proportional occurrence as explanatory factor.

To further test if a species’ vulnerability can be predicted by its land-use response (‘winner’ or ‘loser’ status) and its relation to abiotic soil conditions, we used a general linearized model with Poisson distribution including vulnerability as response factor and the respective land-use parameter or abiotic factor, the number of plots where the species occurred and its total abundance as explanatory factors. Values for grazing and fertilization were square-root transformed prior to statistical analyses and data on abundances and occurrence were log transformed because of data structure.

Finally, we calculated a five-dimensional niche hypervolume (consistent with Hutchinson's n‐dimensional niche concept) as a proxy for the total ‘niche breadth’ of each snail species by multiplying the abundance-weighted standard deviations (AWSD) of all three single land-use components as well as of pH and soil moisture, respectively. The hypervolume was defined for forests and grasslands separately.

Whether the total niche breadth can predict vulnerability was tested using a Spearman rank correlation between the vulnerability and the five-dimensional niche hypervolume.

Availability of data and materials

Snail data obtained by [27] and used in this study are available online under https://www.bexis.uni-jena.de/PublicData/PublicDataSet.aspx?DatasetId=24986. Data on snail vulnerability were obtained from the Red List 2011 according to [43] and snail traits were extracted from [38]. Environmental data and those for land-use intensity in grasslands and forests were obtained from the BExIS database (see Table 6).

References

  1. Poschlod P, Bakker JP, Kahmen S. Changing land use and its impact on biodiversity. Basic Appl Ecol. 2005;6:93–8.

    Article  Google Scholar 

  2. Fischer M, Bossdorf O, Gockel S, Hänsel F, Hemp A, Hessenmöllerd D, Weisser WW, et al. Implementing large-scale and long-term functional biodiversity research: the biodiversity exploratories. Basic Appl Ecol. 2010;11:473–85.

    Article  Google Scholar 

  3. Steinhäußer R, Siebert R, Steinführer A, Hellmich M. National and regional land-use conflicts in Germany from the perspective of stakeholders. Land Use Policy. 2015;49:183–94.

    Article  Google Scholar 

  4. Axelsson R, Angelstam P, Svensson J. Natural forest and cultural woodland with continuous tree cover in Sweden: how much remains and how is it managed? Scand J Forest Res. 2007;22:545–58.

    Article  Google Scholar 

  5. Socher AS, Prati D, Boch S, Müller J, Klaus VH, Hölzel N, Fischer M. Direct and productivity-mediated indirect effects of fertilization, mowing and grazing on grassland species richness. J Ecol. 2012;100:1391–9.

    Article  Google Scholar 

  6. Simons NK, Gossner MM, Lewinsohn TM, Boch S, Lange M, Müller J, Weisser WW, et al. Resource-mediated indirect effects of grassland management on arthropod diversity. PLoS ONE. 2014;9:e107033.R.

    Article  CAS  Google Scholar 

  7. Hooper DU, Chapin FS, Ewel JJ, Hector A. Effects of biodiversity on ecosystem functioning: a consensus of current knowledge. Ecol Monogr. 2005;75:3–35.

    Article  Google Scholar 

  8. Dormann CF, Schweiger O, Augenstein I, Bailey D, Billeter R, De Blust G, Zobel M, et al. Effects of landscape structure and land-use intensity on similarity of plant and animal communities. Global Ecol Biogeogr. 2007;16:774–87.

    Article  Google Scholar 

  9. Chisté MN, Mody K, Kunz G, Gunczy J, Blüthgen N. Intensive land use drives small-scale homogenization of plant- and leafhopper communities and promotes generalists. Oecologia. 2018;186:529–40.

    Article  PubMed  Google Scholar 

  10. Astor T, Lenoir L, Berg MP. Measuring feeding traits of a range of litter-consuming terrestrial snails: leaf litter consumption, faeces production and scaling with body size. Oecologia. 2015;178:833–45.

    Article  PubMed  Google Scholar 

  11. Cameron R. Slugs and Snails. Collins New Naturalist Library, Book 133; HarperCollins Publishers, ePub edition; 2006.

  12. Limondin-Lozouet N, Preece RC. Quaternary perspectives on the diversity of land snail assemblages from northwestern Europe. J Mollus Stud. 2014;80:224–37.

    Article  Google Scholar 

  13. Randolph PA. Influence of environmental variability on land snail population properties. Ecology. 1973;54:933–55.

    Article  Google Scholar 

  14. Schamp B, Horsák M, Hájek M. Deterministic assembly of land snail communities according to species size and diet. J Anim Ecol. 2010;79:803–10.

    PubMed  Google Scholar 

  15. Hylander K, Nilsson C, Jonsson BG, Göther T. Differences in habitat quality explain nestedness in a land snail meta-community. Oikos. 2005;108:351–61.

    Article  Google Scholar 

  16. Hovermann JT, Davis CJ, Werner EE, Skelly DK, Relyea RA, Yurewicz KL. Environmental gradients and the structure of freshwater snail communities. Ecography. 2011;34:1049–58.

    Article  Google Scholar 

  17. Astor T, von Proschwitz T, Strengbom J, Berg MP, Bengtsson J. Importance of environmental and spatial components for species and trait composition in terrestrial snail communities. J Biogeogr. 2017;44:1362–72.

    Article  Google Scholar 

  18. Goodfried GA. Variation in land-snail shell form and size and its causes: a review. System Zool. 1986;35:204–23.

    Article  Google Scholar 

  19. Baur A, Baur B. Individual movement patterns of the minute land snail Punctum pygmaeum (Draparnaud) (Pulmonata: Endodontidae). Veliger. 1988;30:372–6.

    Google Scholar 

  20. Kappes H, Jordaens K, Hendrickx F, Maelfait J-P, Lens L, Backeljau T. Response of snails and slugs to fragmentation of lowland forests in NW Germany. Landsc Ecol. 2009;24:685–97.

    Article  Google Scholar 

  21. Wäreborn I. Changes in the land mollusc fauna and soil chemistry in an inland district in southern Sweden. Ecography. 1992;15:62–9.

    Article  Google Scholar 

  22. Nekola JC. Large-scale terrestrial gastropod community composition patterns in the Great Lakes region of North America. Divers Distrib. 2003;9:55–71.

    Article  Google Scholar 

  23. Martin K, Sommer M. Relationships between land snail assemblage patterns and soil properties in temperate humid ecosystems. J Biogeogr. 2004a;31:531–45.

    Article  Google Scholar 

  24. Martin K, Sommer M. Effects of soil properties and land management on the structure of grassland snail assemblages in SW Germany. Pedobiologia. 2004b;48:193–203.

    Article  Google Scholar 

  25. Horsák M. Mollusc community patterns and species response curves along a mineral richness gradient: a case study in fens. J Biogeogr. 2006;33:98–107.

    Article  Google Scholar 

  26. Denmead LH, Barker GM, Standish RJ, Didham RK. Experimental evidence that even minor livestock trampling has severe effects on land snail communities in forest remnants. J Appl Ecol. 2013;52:161–70.

    Article  Google Scholar 

  27. Wehner K, Renker C, Brückner A, Simons NK, Weisser WW, Blüthgen N. Land-use affects land snail assemblages directly and indirectly bymodulating abiotic and biotic drivers. Ecosphere. 2019;10(5):e02726.

    Article  Google Scholar 

  28. Cameron RAD, Down K, Pannett DJ. Historical and environmental influences on hedgerow snail faunas. Biol J Linn Soc. 1980;13:75–87.

    Article  Google Scholar 

  29. Chapperon C, Seuront L. Space-time variability in environmental thermal properties and snail thermoregulatory behavior. Funct Ecol. 2011;25:1040–50.

    Article  Google Scholar 

  30. Chisté M, Mody K, Gossner MM, Simons NK, Köhler G, Weisser WW, Blüthgen N. Losers, winners, and opportunists: How grassland land-use intensity affects orthopteran communities. Ecosphere. 2016;7(11):e01545.

    Article  Google Scholar 

  31. Mangels J, Fiedler K, Schneider FD, Blüthgen N. Diversity and trait composition of moths respond to land-use intensification in grasslands: generalists replace specialists. Biodivers Conserv. 2017;26:3385–405.

    Article  Google Scholar 

  32. Čejka T, Hamerlík L. Land snails as indicator of soil humidity in Danubian woodland (SW Slovakia). Pol J Ecol. 2009;57:741–7.

    Google Scholar 

  33. Banaszak-Cibicka W, Żmihorski M. Wild bees along an urban gradient: winners and losers. J Insect Conserv. 2011;16:331–43.

    Article  Google Scholar 

  34. Douglas DD, Brown DR, Pederson N. Land snail diversity can reflect degrees of anthropogenic disturbance. Ecosphere. 2013;4:1–14.

    Article  CAS  Google Scholar 

  35. McKinney ML, Lockwood JL. Biotic homogenization: a few winners replacing many losers in the next mass extinction. TREE. 1999;11:450–3.

    Google Scholar 

  36. Williams P, Colla S, Xie Z. Bumblebee vulnerability: common correlates of winners and losers across three continents. Conserv Biol. 2008;23:931–40.

    Article  Google Scholar 

  37. Rader R, Bartomeus I, Tylianakis JM, Lalibert E. The winners and losers of land use intensification: pollinator community disassembly is non-random and alters functional diversity. Divers Distrib. 2014;20:908–17.

    Article  Google Scholar 

  38. Weiner CN, Werner M, Linsenmair KE, Blüthgen N. Land use impacts on mutualistic networks: disproportional declines in specialized pollinators via changes in flower composition. Ecology. 2014;95:466–74.

    Article  PubMed  Google Scholar 

  39. Kühsel S, Blüthgen N. High diversity stabilizes the thermal resilience of pollinator communities in intensively managed grasslands. Nat Commun. 2015;6:7989.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Slatyer RA, Hirst M, Sexton JP. Niche breadth predicts geographical range size: a general ecological pattern. Ecol Lett. 2013;16:1104–14.

    Article  PubMed  Google Scholar 

  41. Welter-Schultes F. European non-marine molluscs. A guide for species identification. Göttingen: Planet Poster Editions; 2012.

    Google Scholar 

  42. Wiese V. Die Landschnecken Deutschlands. 2nd ed. Wiebelsheim: Quelle & Meyer; 2016.

    Google Scholar 

  43. Süßwassermollusken GP. Ein Bestimmungsschlüssel für die Muscheln und Schnecken im Süßwasser der Bundesrepublik Deutschland. Göttingen: Deutscher Jugendbund für Naturbeobachtungen; 2017.

    Google Scholar 

  44. Pearce TA. When a snail dies in the forest, how long will the shell persist? Effect of dissolution and micro-bioerosion. Am Malacol Bull. 2008;26:111–7.

    Article  Google Scholar 

  45. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2010. http://www.R-project.org/. ISBN 3-900051-07-0.

  46. Fox J, Weisberg S. An R companion to applied regression. 2nd ed. Thousand Oaks: Sage; 2011.

    Google Scholar 

  47. Wickham H, Franҫois R, Henry L, Müller K. “dplyr”: A grammar of data manipulation. 2019. http://dplyr.tidyverse.org, https://github.com/tidyverse/dplyr.

  48. Bates D, Maechler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1):1–48.

    Article  Google Scholar 

  49. VanDerWal J, Falconi L, Januchowski S, Shoo L, Storlie C. `SMDTools`: species distribution modelling tools: tools for processing data associated with species distribution modelling exercises. 2014. http://tidyr.tidyverse.org, https://github.com/tidyverse/tidyr.

  50. Falkner G, Obrdlík P, Castella E, Speight MCD. Shelled gastropoda of Western Europe. Munich: Friedrich-Held-Gesellschaft; 2001.

    Google Scholar 

  51. Frömming E. Biologie der mitteleuropäischen Landgastropoden. Berlin: Duncker & Humblot; 1954.

    Google Scholar 

  52. Kahl T, Bauhus J. An index of forest management intensity based on assessment of harvested tree volume, tree species composition and dead wood origin. Nat Conserv. 2014;7:15–27.

    Article  Google Scholar 

  53. Blüthgen N, et al. A quantitative index of land-use intensity in grasslands: integrating mowing, grazing and fertilization. Basic Appl Ecol. 2012;13:207–20.

    Article  Google Scholar 

  54. Kämper W, Weiner C, Kühsel S, Storm C, Thomas ELTZ, Blüthgen N. Evaluating the effects of floral resource specialisation and of nitrogen regulation on the vulnerability of social bees in agricultural landscapes. Apidologie. 2017;48(3):371–83.

    Article  CAS  Google Scholar 

  55. Busch V, Klaus VH, Schäfer D, Prati D, Boch S, Müller J, Hölzel N, et al. Will i stay or will i go? Plant species-specific response and tolerance to high land-use intensity in temperate grassland ecosystems. J Veg Sci. 2019;30(4):674–86.

    Article  Google Scholar 

  56. Jungbluth J, von Knorre D, Bößneck U, Groh K, Hackenberg E, Kobialka, Zettler M, et al. Rote Liste und Gesamtartenliste der Binnenmollusken (Schnecken und Muscheln; Gastropoda et Bivalvia) Deutschlands. 6. überarbeitete Fassung, Stand Februar 2019. Naturschutz Biolog Vielfalt. 2011;70(3):647–708.

Download references

Acknowledgements

We thank the managers of the three Exploratories, Kirsten Reichel-Jung, Iris Steitz, and Sandra Weithmann, Juliane Vogt, Miriam Teuscher and all former managers for their work in maintaining the plot and project infrastructure; Christiane Fischer for giving support through the central office, Andreas Ostrowski for managing the central data base, and Markus Fischer, Eduard Linsenmair, Dominik Hessenmöller, Daniel Prati, Ingo Schöning, François Buscot, Ernst-Detlef Schulze, and the late Elisabeth Kalko for their role in setting up the Biodiversity Exploratories project. Many thanks to all research assistants: Kevin Frank, Wiebke Kämper, Jessica Schneider, Andrea Hilpert, Matteo Trevisan, Matthias Brandt, Sebastian Schmelzle, Tewannakit Mermagen, Kathrin Ziegler, Annika Keil, Andreas Kerner, Katja Gruschwitz, and Kimberly Adam.

Funding

Open Access funding enabled and organized by Projekt DEAL. The work has partly been funded by the DFG Priority Program 1374 “Infrastructure-Biodiversity-Exploratories” (DFG BL860/8-3).

Author information

Authors and Affiliations

Authors

Contributions

KW did the fieldwork, collected and determined snail species, performed the statistical analyses and wrote the manuscript. CR assisted in the species determination and commented on the manuscript. NKS assisted in the statistical analyses and commented on the manuscript. WWW and NB designed the study, NB also assisted in the statistical analyses and the paper writing. All authors have approved to the final version.

Corresponding author

Correspondence to Katja Wehner.

Ethics declarations

Ethics approval and consent to participate

The study complied the fundamental principles of the Basel declaration for research in animals. The investigated species are not at risk of extinction. Fieldwork permits were issued by the responsible state environmental offices of Baden-Württemberg, Thüringen, and Brandenburg.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Appendix 1.

Summary of significant effects of land-use parameters and abiotic factors in forests (forest management index Formi, proportion of non-native tress, proportion of dead wood with saw cuts, proportion of wood harvested, pH and soil moisture) and grasslands (land-use index LUI, mowing, grazing, fertilization, pH and soil moisture) on the community weighted mean of the maximum shell size, the number of offspring, light preference, humidity preference, drought resistance and inundation tolerance. * p < 0.05, ** p < 0.01, *** p < 0.001. ↓ negative effect, ↑ positive effect.

Additional file 2: Appendix 2.

Influence of the abundance-weighted mean (AWM) of the forest management index on the maximum shell size, number of offspring, light preference, humidity preference, drought resistance and inundation tolerance in forests. Species in italics are land-use “winners”, species in bold are land-use “losers”.

Additional file 3: Appendix 3.

Influence of the abundance-weighted mean (AWM) of the proportion of non-native trees on the maximum shell size, number of offspring, light preference, humidity preference, drought resistance and inundation tolerance in forests. Species in italics are land-use “winners”, species in bold are land-use “losers”.

Additional file 4: Appendix 4.

Influence of the abundance-weighted mean (AWM) of the proportion of deadwood with saw cuts on the maximum shell size, number of offspring, light preference, humidity preference, drought resistance and inundation tolerance in forests. Species in italics are land-use “winners”, species in bold are land-use “losers”.

Additional file 5: Appendix 5.

Influence of the abundance-weighted mean (AWM) of the proportion of wood harvested on the maximum shell size, number of offspring, light preference, humidity preference, drought resistance and inundation tolerance in forests. Species in italics are land-use “winners”, species in bold are land-use “losers”.

Additional file 6: Appendix 6.

Influence of the abundance-weighted mean (AWM) of soil pH on the maximum shell size, number of offspring, light preference, humidity preference, drought resistance and inundation tolerance in forests. Species in italics are land-use “winners”, species in bold are land-use “losers”.

Additional file 7: Appendix 7.

Influence of the abundance-weighted mean (AWM) of soil moisture on the maximum shell, size number of offspring, light preference, humidity preference, drought resistance and inundation tolerance in forests. Species in italics are land-use “winners”, species in bold are land-use “losers”.

Additional file 8: Appendix 8.

Influence of the abundance-weighted mean (AWM) of land-use intensity (LUI) on the maximum shell size, number of offspring, light preference, humidity preference, drought resistance and inundation tolerance in grasslands. Species in italics are land-use “winners”, species in bold are land-use “losers”.

Additional file 9: Appendix 9.

Influence of the abundance-weighted mean (AWM) of mowing on the maximum shell size, number of offspring, light preference, humidity preference, drought resistance and inundation tolerance in grasslands. Species in italics are land-use “winners”, species in bold are land-use “losers”.

Additional file 10: Appendix 10.

Influence of the abundance-weighted mean (AWM) of grazing on the maximum shell size, number of offspring, light preference, humidity preference, drought resistance and inundation tolerance in grasslands. Species in italics are land-use “winners”, species in bold are land-use “losers”.

Additional file 11: Appendix 11.

Influence of the abundance-weighted mean (AWM) of fertilization on the maximum shell size, number of offspring, light preference, humidity preference, drought resistance and inundation tolerance in grasslands. Species in italics are land-use “winners”, species in bold are land-use “losers”.

Additional file 12: Appendix 12.

Influence of the abundance-weighted mean (AWM) of soil pH on the maximum shell size, number of offspring, light preference, humidity preference, drought resistance and inundation tolerance in grasslands. Species in italics are land-use “winners”, species in bold are land-use “losers”.

Additional file 13: Appendix 13.

Influence of the abundance-weighted mean (AWM) of soil moisture on the maximum shell size, number of offspring, light preference, humidity preference, drought resistance and inundation tolerance in grasslands. Species in italics are land-use “winners”, species in bold are land-use “losers”.

Additional file 14: Appendix 14.

Relation of the abundance-weighted means (AWM) of the forest management index, proportion of non-native trees, proportion of dead wood with saw cuts, proportion of wood harvested, pH and soil moisture and the proportional occurrence of a certain species in forests.

Additional file 15: Appendix 15.

Relation of the abundance-weighted means (AWM) of the land-use intensity, mowing, grazing, fertilization, pH and soil moisture and the proportional occurrence of a certain species in forests.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Wehner, K., Renker, C., Simons, N.K. et al. Narrow environmental niches predict land-use responses and vulnerability of land snail assemblages. BMC Ecol Evo 21, 15 (2021). https://doi.org/10.1186/s12862-020-01741-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12862-020-01741-1

Keywords