CC BY
82
Rcqiilatorv Mechanisms I
mTnosystems
Regulatory Mechanisms
in Biosystems
ISSN 2519-8521 (Print) ISSN 2520-2588 (Online) Regul. Mech. Biosyst., 9(2), 198-202 doi: 10.15421/021829
Fitting competing models and evaluation of model parameters of the abundance distribution of the land snail Vallonia pulchella (Pulmonata, Valloniidae)
O. N. Kunakh*, S. S. Kramarenko**, A. V. Zhukov*, A. S. Kramarenko**, N. V. Yorkina***
*Oles Honchar Dnipro National University, Dnipro, Ukraine
**Mykolayiv National Agrarian University, Mykolayiv, Ukraine
***Bohdan Khmelnitsky Melitopol State Pedagogical University, Melitopol, Ukraine
Article info
Received 15.03.2018 Received in revised form
21.04.2018 Accepted 25.04.2018
Oles Honchar Dnipro National University, Gagarin av., 72, Dnipro, 49010, Ukraine. Tel.: +38-098-507-96-82. E-mail: [email protected]
Mykolayiv National Agrarian University, Georgiy Gongadze st., 9, Mykolayiv, 54020, Ukraine. E-mail: [email protected]
Bohdan Khmelnitsky Melitopol State Pedagogical University, Hetmanska st., 20, Melitopol, 72312, Ukraine. E-mail:
Introduction
Kunakh, O. N., Kramarenko, S. S., Zhukov, A V., Kramarenko, A S., & Yorkina, N. V. (2018). Fitting competing models and evaluation of model parameters of the abundance distribution of the land snail Vallonia pulchella (Pulmonata, Valloniidae). Regulatory Mechanisms in Biosystems, 9(2), 198-202. doi:10.15421/021829
This paper summarizes the mechanisms behind the patterning of the intra-population abundance distribution of the land snail Vallonia pulchella (Müller, 1774). The molluscs were collected in recultivated soil formed on red-brown clays (Pokrov, Ukraine). Data obtained in this study reveal that V. pulchella population abundance ranges from 1 to 13 individuals per 100 g of soil sample. To obtain estimates of the mean, three models were used: the model of the arithmetic mean, the Poisson model and a log-normal model. The arithmetic mean of the occurrence of this species during the study period was 1.84 individuals/sample. Estimation of the average number of molluscs in one sample calculated using the Poisson model is lower and equals 1.40 individuals/sample. The distribution of the number of individuals in a population was described by the graphics "rank - abundance". The individual sample plot sites with molluscs may be regarded as equivalents of individual species in the community. For the analysis, the following models were used: broken sticks model, niche preemption model, log-normal model, Zipf model, and Zipf-Mandelbrot model. Applying the log-normal distribution gives a lower estimate of the mean density at 1.28 individuals/sample. Median value and mode is estimated at 1.00 individuals/sample. The Zipf-Mandelbrot model was shown as the most adequate to describe distribution of the V. pulchella population within the study area. The Zipf-Mandelbrot model belongs to the family of so-called non-Gaussian distributions. This means that the sample statistics do not possess asymptotic properties and by increasing the sample size, they tend to infinity, and are not close to the values of the general population. Therefore, the average value of the random variable that describes the non-Gaussian distribution has no statistical meaning. From an environmental point of view, this means that within the study area the capacity of the habitat is large, and for some combination of environmental conditions the rapid growth of the abundance of a given species is possible.
Keywords: abundance; intra-population structure; land snails
There are plenty of publications devoted to habitat selection by land snails based on the investigation of mollusc communities from spatially different biotopes which differ in vegetation cover, soil type, moisture level (Millar & Waite, 1999; Martin & Sommer, 2004; Muller et al., 2005; Weaver et al., 2006). The number of invertebrates on a particular plot changes randomly. It is impossible to predict in advance for any moment of time at a given plot which of the factors will operate within the optimum range and which at an extreme range for any given species of invertebrate (Brygadyrenko, 2015a). The calcium content in the soil, pH, the soil texture (Ondina et al., 2004), and the content of exchangeable cations and aluminum (Ondina et al., 1998) are the most important among the edaphic factors which affect the molluscs. The soil moisture plays a significant role in forming ecological niches of molluscs (Nekola, 2003). However, P. Ondina et al. (2004) noted the restrictions of the importance of soil moisture data at any given moment because of the considerable changeability of this parameter. To help address this situation phytoindication data application is an appropriate method for assessment of land snails' ecological characteristics and community properties (Horsak et al., 2007; Dvorakova & Horsak, 2012). Ellenberg phytoindication scales have been effectively used to
explain the habitat preferences of the land snail Vertigo geyeri in Poland and Slovakia (Schenkova et al., 2012). Interaction of litter macrofauna with herbaceous plants has been investigated in natural forests of Ukraine's steppe zone (Brygadyrenko, 2015b) and in poplar plantations within city parks (Faly & Brygadyrenko, 2014). The importance of soil factors in the spatial distribution, abundance and diversity of mollusc communities has been revealed on a large-scale level of research (Nekola & Smith, 1999; Jufickova et al., 2008; Szybiak et al., 2009). The issue of the spatial scale and the hierarchy of factors acting on molluscs are of special interest (Nekola & Smith, 1999; Bohan et al., 2000; McClain & Nekola, 2008; Mysak et al., 2013).
Habitat is attributed by the presence of a certain range of resources and ecological conditions which make it possible for a given species to occupy, survive and reproduce in a given area (Hall et al., 1997). For identification of the environmental characteristics which make a territory appropriate for the existence of a species, it is important to investigate the selection of habitat types (Calenge & Basille, 2008).
The mollusc Vallonia pulchella (Muller, 1774) is Holarctic species distributed around the world at high latitudes. In Europe it prefers wetter habitats such as wet grasslands and swamps, as well as dry dunes and meadows (Kerney & Cameron, 1979). In Ukraine it is found in moderately dry and wet grasslands habitats (Gural-Sverlova & Gural, 2012).
Theie have been few publications on the results of aut- and demecologi-cal research on the land snail V. pulchella. In the 1930s, research was carried out on the reproduction of this mollusc (Whitney, 1937, 1941), which was continued in this century (Kuznik-Kowalska & Prockow, 2016). It was revealed that the density of V. pulchella is affected by many soil physico-chemical properties, especially calcium and magnesium concentration and pH of water extract (Hermida et al., 1993, 2000). Furthermore, interrelation was found between the abundance of V. pulchella and pH and electrical conductivity of the water in marshy lowland central Scandinavia (Schenkova et al., 2015). Local trends, as well as the mosaic nature of the organization of the soil body determine the structure of the vegetation cover. This explains the role of indicators in the organization of ecological niche of mollusc V. pulchella, the dimension of which exceeds the size of the ecological space of an individual of this species. Soil mechanical impedance, as well as phytoindication indicators of vegetation, determines the peculiarities of marginality and specialize-tion of the ecological niche of V. pulchella. It is the measurement of these indicators that makes it possible to draw a map of habitat preferences for V. pulchella. Thus, ecological niche optima may be presented by integral variables such as marginality and specialization axes and may be plotted in geographic space by means of habitat suitability index reproduction (Yorkina et al., 2018).
Thus, the main aim of our work is to select the optimal model and to evaluate model parameters of the abundance distribution of the land snail V. pulchella (Müller, 1774) (Gastropoda, Pulmonata, Valloniidae).
Material and methods
The studies were conducted at the research station of Dnipro Agrarian and Economic University in Pokrov in June 2012. Sampling was carried out on a type of technogenic soil formed on red-brown clays (geographic coordinates of the southwest corner of the polygon are 47°38'55" N, 34°08'33" E) (Zhukov et al., 2013; Zverkovskyi et al., 2017; Kharytonov et al., 2018; Yorkina et al., 2018). According to WRB 2007 (IUSS Working group WRB, 2007), the examined soil belongs to the RSG Technosols. A perennial legume-grass agrophytocenosis was cultivated at this research station between 1995 and 2003, following which the process of naturalization of vegetation began (Zhukov et al., 2017).
The area within which the samples were taken consisted of seven transects of 15 samples each. The sampling points formed a regular grid of a mesh size of 3 m. From the center of each sampling site a soil sample weighing 100 g was collected. Each sample was examined in the laboratory with a x10 binocular microscope MBS-9 and the number of live specimens of V. pulchella was noted.
To obtain estimates of the mean, three models were used: the model of the arithmetic mean, the Poisson model and a log-normal model. Estimation of the mean from assumptions about the distribution of the random variable according to Poisson was obtained by the formula (Shebanin et al., 2008):
d=- ln n», a)
n
where D is the mean estimation; n0 is number of samples without the presence of the species; n is total number of samples.
Statistical calculations were performed using Statistica 7.0 program (Khalafyan, 2007) and a 95% confidence interval for the mean estimation obtained using different models has been calculated on the basis of the bootstrap approach using bootES package (Kirby & Ger-lanc, 2013). The distribution of the number of individuals in a population may also be described by the graphics "rank - abundance", which are often used in community ecology (Whittaker, 1965). In this case, the individual sample plot sites with molluscs may be regarded as equivalents of individual species in the community. For the analysis in this case, the following model can be used:
- broken sticks model (MacArthur, 1957):
a - N h1 ; (2)
a- s hk
- Motomura model (the Whittaker niche preemption model) (Mo-tomura, 1932):
ar - Na{1 -a)'1 ; (3)
(5)
- log-normal model (Preston, 1948, 1962):
ar = exp[log(^) + log(o-)®];
- Zipf model (Zipf, 1949):
ar = Np1rr;
- Zipf-Mandelbrot model (Mandelbrot, 1983):
ar = Nc(r + P)r ; (6)
where a is the expected abundance of species of rank r; S is the number of species; N is the number of individuals; 0 is a standard normal distribution function; pi is the estimated proportion of the most abundant species; a, /, a, y, f and c are the parameters in each model.
The degree of adequacy of the model was evaluated using Akaike's information criterion (AIC) and the Bayesian information criterion (BIC). The best model has the lowest AIC and BIC.
Plant and soil data used for distance matrix calculation were discussed previously in our publication (Yorkina et al., 2018). Statistical calculations were performed with the help of the Statistica 7.0 program and the project for statistical computations R (www.r-project.org) using vegan (Oksanen et al., 2017).
Results
In examining 105 samples, 193 V. pulchella were detected. Thus, the arithmetic mean of the occurrence of this species in sod lithogenic soils on red-brown clays during the study period was 1.84 individuals per of 100 g soil sample (Table 1). Thus, the observed distribution of V. pulchella abundance in one sample is adequately described by the Poisson law (Kolmogorov-Smirnov test: dKS = 0.089, P > 0.05) (Fig. 1). This type of discrete distribution is described by a random flow (Puza-chenko, 2004; Shebanin et al., 2008). Estimation of the average number of molluscs in one sample calculated using the Poisson model is lower and equals 1.40 individuals/example (Table 1).
Table 1
Estimates of V. pulchella population density based on different models, individuals per of 100 g soil sample
Model D Dboot ± SD 95% CI
Arithmetic mean 1.84 1.79 ± 0.23 1.35 - 2.27
Poisson distribution 1.40 1.43 ± 0.20 1.08 - 1.85
Log-normal distribution 1.28 1.26 ± 0.14 1.03 - 1.52
Notes: D - density estimates of the population; D^u -estimates of the population; SD - standard deviation; CI -
bootstrap-based density ■ confidence interval.
The pronounced asymmetric character of the distribution makes the arithmetic mean an insignificant value of the general population mean. It is quite logical that applying the log-normal distribution gives a lower estimate of the mean density at 1.28 individuals/sample (Table 1). This value is closest to the median value and mode, which is estimated at 1.00 individuals/sample.
dKS = 0,089 (ns)
45
40
35
30
o C 25
u
20
u
s-
Hh 15
10
0 1 2 3 4 5 6 7 8 9 10 11 12 13 Individuals/example
Fig. 1. Histogram of the empirical distribution of V. pulchella per 100 g of soil sample (dashed line marks the theoretical Poisson distribution for X = 1.40)
The use of curves "ranked-abundance" shows that the Zipf-Man-delbrot model is the most adequate to describe the abundance of the
population of V. pulchella within the study area, as the Akaike information criterion for this model is the smallest among all least studied (Table 2).
Table 2
Statistics for the various models
of distribution "rank-abundance" experimental data
Model Deviation Information criterion
Akaike (AIC) Bayesian (BIC)
Broken sticks 30.30 232.29 232.29
Niche preemption 16.83 220.82 223.19
Log-normal 8.78 214.77 219.51
Zipf 7.23 213.23 217.96
Zipf-Mandelbrot 3.37 211.36 218.47
It should be noted that the Zipf-Mandelbrot model belongs to the family of so-called non-Gaussian distributions (Khaytun, 2005). This means that the sample statistics do not possess asymptotic properties
Table 3
Correlations of distance matrices (Mantel test)
and by increasing the sample size, they tend to infinity, and are not close to the values of the general population. Therefore, the average value of the random variable that describes the non-Gaussian distribution has no statistical meaning. From an environmental point of view, this means that within the study area the capacity of the habitat is large, and for some combination of environmental conditions the rapid growth of the abundance of a given species is possible.
The Mantel test allows us to compare the correlation matrices among themselves (Table 3). The Mantel test for matrices of distance measures based on the number of molluscs (Euclidean distance) and measured ecological characteristics (Euclidean distance) indicates a reliable correlation of these matrices (rm = 0.15, P = 0.02). Mantel tests indicate that in a complex of environmental indicators that affect molluscs, a key role is played by edaphic factors (positive correlation) and phytoindication indicators (negative correlation). The space (the matrix of geographical distances) does not play an important role in the variability of the molluscs' density.
Factor All environmental Edaphic Phytoindication scales Spatial
indicators indicators all climate edaphic coordinates
Mantel Test r=0,15, r = 0,15, r = -0,11, r=-0,05, r=-0,08, r = 0,04,
p = 0,014 p = 0,016 p = 0,99 p = 0,79 p = 0,93 p = 0,120
A partial Mantel test with control variables
Edaphic r = -0,12, r = -0,13, r=-0,06, r=-0,09, r = 0,02,
indicators p = 0,997 p = 0,994 p = 0,859 p = 0,96 p = 0,329
Phytoindication scales r=0,16, p = 0,009 r = 0,16, p = 0,009 - - - r = 0,07, p = 0,077
Climate r=0,15, r = 0,15, r = -0,10, r=-0,08, r = 0,04,
p = 0,016 p = 0,015 p = 0,978 p = 0,937 p = 0,153
Edaphic r=0,15, r = 0,15, r = -0,08, r=-0,05, r = 0,05,
p = 0,015 p = 0,009 p = 0,954 p = 0,784 p = 0,123
Spatial r=0,14, r = 0,14, r = -0,12, r=-0,05, r=-0,09,
coordinates p = 0,023 p = 0,028 p = 0,987 p = 0,806 p = 0,961
It should be noted that the information value for describing the variation in the number of molluscs has a complex of phytoindication indicators, and does not separate the matrices of climatic and edaphic indicator values. The use of edaphic indicators as a control variable in the partial Mantel test leads to a change in the correlation coefficient of the distance matrices in terms of the number of molluscs and the matrix of ecological indicators (rm = -0.12, P = 0.997). This confirms the result that edaphic factors and phytoindication indicators mark opposite tendencies of influence of environmental factors on the mollusc population, while the Mantel test between the matrices of edaphic and phytoindicating distances is positive (rm = 0.11, P = 0.013).
Phytoindication scales as drivers in the partial Mantel test do not affect the correlation of molluscs and both environmental indicators in general and edaphic indicators.
Thus, the analysis of the results of Mantel tests indicates a specific character of the influence of edaphic factors and vegetation on the population of the molluscs V. pulchella and the absence of influence of spatial factors on the chosen scale. It is likely that the selected set of indicators can fully describe the spatial variation of the population of V. pulchella within the studied test site.
Discussion
Distribution of the species in the loci can be described by a number of models. Abundance curves are widely used to describe dominance structures in communities of living organisms. This tool can also be used to describe the statistical distribution of the individuals in a population. The difference consists in the fact that the x-axis will not be delayed in the rank of species abundance and species abundance levels in the locus of a given population. Ideally, the extent to which the actual distribution of a particular model makes it possible to make a choice in favour of one of the theories that explain the structuring of the population (or community) (Pielou, 1975). Geometric series (Motomura, 1932) was proposed to describe the benthic communities of a lake. The parameter k of the distribution is a measure of the complexity of the species composition of the community. According to
Tokeshi (1993), Motomura applied a geometric series as the simplest mathematical form for the description of an ecological community, although later this model was interpreted as a reflection of the distribution of resources among species (Ferreira & Petrere-Jr., 2008). From this viewpoint, the distribution should be described by a situation in which the dominant species uses a proportion k of the whole initially available resource, leaving accessible proportion 1-k. The next dominant species uses the same fraction k of the remaining available resource. This goes for all species of the community. May (1975) believes that the geometric series is a situation in which all species of a community are energy equivalent and value of ecological interactions is proportional to the abundance of species - more abundant species require more energy from the system. Geometric series take into account only the number of species and neglect the influence of the size of the animals on the energy needs of a species (Ferreira, Petrere-Jr., 2008).
Logarithmic series distribution proposed by Fischer et al. (1943) slightly resembles a hyperbola decreasing with increasing number of species and can be predicted from the expression: ax, ax2/2, ax3/3, ..., ax"/n, where ax7n represents the number of species with n number of individuals. The number x can be obtained from the equation
S/N = [(1 -x)/x] *[-ln (1-x)] and in practice is in the range 0.9-1.0 (Magurran, 2004). The constant a is independent of the sample size and may be viewed as an index of diversity, while remaining a robust assessment even in the case where the data is not very well matched logarithmic series (F erreira & Petrere-Jr., 2008).
Truncated lognormal distribution was proposed by Preston (1948). With the species abundances in the histogram in logarithmic scale, the author obtained a curve that is a good description of a large number of community data. The R classes resulting logarithms are called octaves, wherein each class represented twice the number of species of the previous class (1, 2, 4, 8, 16 ...). The lognormal model has two parameters: ^ - mean value of the logarithmic score abundance and a -standard rating average logarithmic abundance deviation (Ferreira & Petrere-Jr., 2008). For real samples distribution is truncated on the left.
The area behind this point represents species that is not included in the sample, and it decreases with increasing sample size. MacArthur (1957) suggested that the ecological niches in a community can be compared with a stick of unit length, at which n-1 are randomly generated points form n segments of length proportional to the abundance of each species in the community. MacArthur's broken stick model does not contain parameters for evaluation.
Model Zipf (Zipf, 1949) has two parameters - abundance estimation dominant form y and the parameter p, which indicates the probability of occurrence of the species in the community. The Model Zipf-Mandelbrot (Mandelbrot, 1983) has a third parameter a, which may be regarded as ecological variety of potential environmental or ecological niche diversity.
The geometrical model of the distribution of species abundance is consistent with the hypothesis interception niches Whittaker (1972). It assumes that every species of community sites in decreasing order of abundance (increase in rank) uses a constant (K-th) of the remaining resources of the community. For example, if the strongest competitor (dominant) occupies 70% of the niche space using the appropriate community resources share, the second most important species of able to take a similar share of the remainder of the first niche space, the third type - the same piece of space left over from the first and the second kind, etc. According to this model, the proportion of community resources used by the most abundant species (dominant) is not a special case (the result of the biological and ecological characteristics or accidental circumstances), and reflects the general character of the distribution of niche space between species in different conditions. Hence the community with a higher dominance of the most competitive species must be characterized as having not only fewer resources available for related species, but also a more "rigid" means of distribution of resources among them (i.e. higher K values) that ceteris paribus conditions may affect their species richness.
There are other types (models) of species abundance patterns, a review of which is available in many publications (Whittaker, 1972; Megarran, 2004; Ferreira & Petrere-Jr., 2008). They assume a more or less uniform, compared to the geometric model, allocation of resources and thus more or less marked predominance of one species. The former include, for example, log-normal model and the "broken stick" MacArthur; to the second - hyperbolic model.
According to the log-normal model (for the analysis of biological communities, this was first used by Preston), K values in the species of the first rank higher than the species next few ranks, and in accordance with the hypothesis of random boundaries between niches MacArthur vice versa. At the same time as the lognormal model and MacArthur model typical for communities with a relatively low level of dominance. Therefore, the allocation of resources between the associated views in the three models (hyperbolic lognormal and MacArthur) is related to the level less than the first domination rank species than in geometric.
Particular attention is paid to the types of distribution in the community ecology and population ecology due to the need to open or create a picture of the mechanisms of the relationship of the community or population with the resources offered by the habitat. It is expected that the community or the population is subject to the theoretical concepts, on the basis of which are derived the corresponding types of distributions. The fact that the population size distribution of V. pulchella is not subject to a fully stochastic processes (Poisson distribution model of the broken sticks) indicates that the selected site is not entirely uniform for the species under study. It should be noted that an alternative statistical broken stick model is a little more realistic, which can only be the basis for the assumption of the presence of any processes that may underlie the alternative models.
The best of the studied models to describe the number of molluscs is the model of the Zipf-Mandelbrot. This type of distribution is often found in the description of complex organized systems (Khaytun, 2005).The non-Gaussian nature of this distribution makes it difficult to use traditional statistics to describe the corresponding processes. With respect to the type of population, Mandelbrot-Zipf distribution systems can arise as the result of significant increase in population numbers in the formation of optimal environmental conditions. In other words, the
number is not in this situation, a marker of environmental conditions in the sense that the relationship is not proportional. Then a great informational value is simply the fact of the presence (or absence) of the form.
Studies of terrestrial micromolluscs has always been accompanied by significant methodological difficulties related to their small size. For example, the diameter of the shell of many species of Vallonia Risso, 1926 does not exceed 2-3 mm (Gerber, 1996). However, the procedure for assessing the abundance of species, using a widespread approach based on test plots, characterized by significant errors and obtained on the basis of such an assessment approach, cannot always be trusted unconditionally. In addition, these estimates may differ even by two-four orders of magnitude and cannot always be clear, these differences are fundamental and are the consequences of violations of population estimation procedure. For example, the population density of V. pulchella in alder and oak forests of Belarus was 4-8 individuals/m2 (Zemoglyadchuk, 2005), on artificial Robinia stands in a recultivated area near the town of Zhovti Vodi it was 5.6 individuals/m2 (Kul'bachko & Unkov'ska, 2008), whereas in ash-elm forests in Poland it did not exceed the average of 0.13 individuals/m2 (Koralewska-Batura & Bloszyk, 2007), and in floodplain forests of Slovakia it was 0.07 individuals/m2 (Cejka & Hamerlik, 2009). On the other hand, Hermida et al. (1993) estimated the average density for three studied populations of V. pulchella in Spain (meadow and forest habitats, as well as near the river bank) to be 5.9-10.1 individuals/m2, reaching on separate test plots values of the order of 200 individuals/m2. In bush willow depressions in Kazakhstan Uvalieva (1990) a density of 224 individuals/m2 was estimated.
Therefore, in principle other methods have been proposed for abundance estimates for litter micromolluscs, for example, collecting a significant amount of topsoil, moss, lichen, leaf litter (volume 10-12 liters), followed by hand sorting of these samples in a laboratory. In addition, for complete extraction of the shells (if possible), a set of sieves with decreasing mesh sizes may be used (Horsak, 2003).
For typical soil species it has been proposed to conduct sampling of soil (specific weight), followed by its analysis under a microscope in the laboratory. This method originated in paleozoology (Evans, 1972). Using this approach, Davies et al. (1996) found the abundance of V. pulchella to be 0.4-40.4 individuals per 100 g sample of soil on chalky soils in UK, and for the Jurassic limestone in the center of Krakow (Poland), its abundance was 1-22 individuals per 100 g of soil sample (Golas-Siarzewska, 2013). The results thus obtained are entirely consistent with the data obtained in this study - 1-13 individuals per 100 g of soil sample, which may indicate the higher accuracy of this method.
References
Bohan, D. A., Glen, D. M., Wiltshire, C. W., & Hughes, L. (2000). Parametric intensity and the spatial arrangement of the terrestrial mollusc herbivores
Deroceras reticulatum and Arion intermedial. Journal of Animal Ecology, 69(6), 1031-1046.
Brygadyrenko, V. V. (2015a). Influence of moisture conditions and mineralization of soil solution on structure of litter macrofauna of deciduous forests of Ukraine's steppe zone. Visnyk of Dnipropetrovsk University. Biology, Ecology, 23(1), 50-65 (in Ukrainian). Brygadyrenko, V. V. (2015b). Influence of tree crown density and density of the herbaceous layer on the structure of litter macrofauna of deciduous forests of Ukraine's steppe zone. Visnyk of Dnipropetrovsk University. Biology, Ecology 23(2), 134-148 (in Ukrainian). Calenge, C., & Basille, M. (2008). A general framework for the statistical exploration of the ecological niche. Journal of Theoretical Biology, 252(4), 674-685. Cejka, T., & Hamerlik, L. (2009). Land snails as indicators of soil humidity in Da-nubian woodland (SW Slovakia). Polish Journal of Ecology, 57(4), 741-747. Davies, P., & Grimes, C. J. (1999). Small-scale spatial variation of pasture mollus-can faunas within a relic watermeadow system at Wylye, Wiltshire, UK. Journal of Biogeography, 26(5), 1057-1063. Dvorakova, J., & Horsak, M. (2012). Variation of snail assemblages in hay meadows: Disentangling the predictive power of abiotic environment and vegetation. Malacologia, 55(1), 151-162. Evans, J. G. (1972). Land snails in archaeology. Seminar Press, London and New York. Faly, L. I., & Brygadyrenko, V. V. (2014). Patterns in the horizontal structure of litter invertebrate communities in windbreak plantations in the steppe zone of the Ukraine. Journal of Plant Protection Research, 54(4), 414-420.
Ferreira, F. C., & Petrere Jr., M. (2008). Comments about some species abundance patterns: Classic, neutral, and niche partitioning models. Brazilian Journal of Biology, 68, 1003-1012.
Fisher, R. A., Corbert, A S., & Williams, C. B. (1943). The relation between the number of species and the number of individuals in a random sample of animal population. The Journal of Animal Ecology, 12, 42-58.
Gerber, J. (1996). Revision der Gattung Vallonia Risso 1826 (Mollusca: Gastropoda: Valloniidae). Hemmen, Wiesbaden (Schriften zur Malakozoologie aus dem Haus der Natur-Cismar, Haus der Natur Cismar. Issue 8).
Golas-Siarzewska, M. (2013). Malacofauna of the Wawel Hill in Cracow (Poland) - a quarter of a century after its first description. Folia Malacologica, 21(1), 19-23.
Gural-Sverlova, N. V., & Gural, R. I. (2012). Vyznachnyk nazemnyh moljuskiv Ukrai'ny [Guide to the terrestrial mollusks of Ukraine]. Lviv (in Ukrainian).
Hall, L. S., Krausman, P. R., & Morrison, M. L. (1997). The habitat concept and a plea for standard terminology. Wildlife Society Bulletin, 25(1), 173-182.
Hermida, J., Ondina, M. P., & Rodriguez, T. (2000). The relative importance of edaphic factors on the distribution of some terrestrial gastropod species: Aut-ecological and synecological approaches. Acta Zoologica Academiae Scien-tiarum Hungaricae, 46, 265-274.
Hermida, J., Outeiro, A., & Rodriguez, T. (1993). Aspectos ecologicos de Vallonia pulchella (Muller, 1774) (Gastropoda, Pulmonata). Nova Acta Cientifica Compostelana (Bioloxia), 4, 129-135.
Horsak, M. (2003). How to sample mollusc communities in mires easily. Malaco-logica Bohemoslovaca, 2, 11-14.
Horsak, M., Hajek, M., Tichy, L., & Jufickovä, L. (2007). Plant indicator values as a tool for land mollusc autecology assessment. Acta Oecologica, 32(2), 161-171.
IUSS Working Group WRB (2007). World reference base for soil resources 2006, first update 2007. World Soil Resources Reports No 103. FAO, Rome.
Jufickova, L., Horsak, M., Cameron, R., Hylander, K., Mikovcova, A., Hlavac, J. C., & Rohovec, J. (2008). Land snail distribution patterns within a site: The role of different calcium sources. European Journal of Soil Biology, 44(2), 172-179.
Kearney, M. P., & Cameron, R. A. D. (1979). A field guide to the land snails of Britain and North West Europe. William Collins Sons and Company, Ltd., London.
Khalafyan, A. A. (2007). Statistica 6: Statisticheskij analiz dannyh [Statistica 6: Statistical Analysis of the Data]. Bean-Pres, Moscow (in Russian).
Kharytonov, M., Babenko, M., Velychko, O., & Pardini, G. (2018). Prospects of medicinal herbs management in reclaimed minelands of Ukraine. Ukrainian Journal of Ecology, 8(1), 527-532.
Khaytun, S. D. (2005). Kolichestvennyy analiz sotsial'nykh yavleniy. Problemy i perspektivy [Quantitative analysis of social phenomena. Problems and prospect]. KomKniga Publisher, Moscow (in Russian).
Kirby, K. N., & Gerlanc, D. (2013). BootES: An R package for bootstrap confidence intervals on effect sizes. Behavior Research Methods, 45(4), 905-927.
Koralewska-Batura, E., & Bloszyk, J. (2007). Stability of malacocoenoses in an ash-elm forest based on two-year observations. Folia Malacologica, 15(4), 175-180.
Kul'bachko, Y. L., & Unkovs'ka, S. I (2008). Kharakteristika fauny nazemnykh bryu-khonogikh mollyuskov iskusstvennykh drevesnykh nasazhdeniy i rekul'tiviro-vannykh territoriy v g Zheltye Vody [Characteristic of terrestrial snails fauna of artificial stands and rehabilitated territories in the town of Zhovti Vody]. Visnyk of Dnipropetrovsk University. Biology, Ecology, 16(1), 128-132 (in Russian).
Kuznik-Kowalska, E., & Prockow, M. (2016). Reproductive biology and growth of two Vallonia species in laboratory conditions (Gastropoda: Eupulmonata: Valloniidae). Folia Malacologica, 24(4), 265-273.
MacArthur, R. H. (1957). On the relative abundance of bird species. Proceedings of the National Academy of Sciences of the United States of America, 43(3), 293-295.
Magurran, A. E. (2004). Measuring biological diversity. Blackwell Publishing, Oxford.
Mandelbrot, B. B. (1983). The fractal geometry of nature. W. H. Freeman and Company, New York.
Martin, K., & Sommer, M. (2004). Relationships between land snail assemblage patterns and soil properties in temperate-humid forest ecosystems. Journal of Biogeography, 31(4), 531-545.
May, R. M. (1975). Patterns of species abundance and diversity. In: Cody, M. L., & Diamond, J. M. (Eds.). Ecology and evolution of communities. Harvard University Press, Cambridge.
McClain, C. R., & Nekola, J. C. (2008). The role of local-scale processes on terrestrial and deep-sea gastropod body size distributions across multiple scales. Evolutionary Ecology Research, 10(1), 129-146.
Millar, A. J., & Waite, S. (1999). Mollusks in coppice woodland. Journal of Con-chology, 36, 25-48.
Motomura, I. (1932). On the statistical treatment of communities. Zoological Magazine, 44, 379-383.
Müller, J., Strätz, C., & Hothorn, T. (2005). Habitat factors for land snails in European beech forests with a special focus on coarse woody debris. European Journal of Forest Research, 124(3), 233-242.
Mysak, J., Horsak, M., Svobodova, E., & Cernohorsky, N. (2013). Small-scale distribution of terrestrial snails: Patterns of species richness and abundance related to area. Journal of Molluscan Studies, 79(2), 118-127.
Nekola, J. C. (2003). Large-scale terrestrial gastropod community composition patterns in the Great Lakes region of North America. Diversity and Distributions, 9(1), 55-71.
Nekola, J. C., & Smith, T. M. (1999). Terrestrial gastropod richness patterns in Wisconsin carbonate cliff communities. Malacologia, 41(1), 253-270.
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O'Hara, R B., Simpson, G. L., Solymos, P., Stevens, M. H. H., & Wagner, H. (2017). Community Ecology Package. R package version 2.4-5.
Ondina, P., Hermida, J., Outeiro, A., & Mato, S. (2004). Relationships between terrestrial gastropod distribution and soil properties in Galicia (NW Spain). Applied Soil Ecology, 26(1), 1-9.
Ondina, P., Mato, S., Hermida, J., & Outeiro, A. (1998). Importance of soil exchangeable cations and aluminium content on land snail distribution. Applied Soil Ecology, 9(1), 229-232.
Pielou, E. C. (1975). Ecological diversity. Wiley, New York.
Preston, F. W. (1948). The commonness, and rarity, of species. Ecology, 29(3), 254-283.
Preston, F. W. (1962). The canonical distribution of commonness and rarity. Ecology, 43, 185-215.
Puzachenko, Y. G. (2004). Matematicheskie metody v ekologicheskikh i geogra-ficheskikh issledovaniyakh [Mathematical Methods in Ecological and Geographical Studies]. Academy, Moscow (in Russian).
Schenkova, V., Horsak, M., Hajek, M., Hajkova, P., & Dite, D. (2015). Mollusc assemblages of Scandinavian fens: Species composition in relation to environmental gradients and vegetation. Annales Zoologici Fennici, 52, 1-16.
Schenkova, V., Horsak, M., Pleskova, Z., & Pawlikowski, P. (2012). Habitat preferences and conservation of Vertigo geyeri (Gastropoda: Pulmonata) in Slovakia and Poland. Journal of Molluscan Studies, 78(1), 105-111.
Shebanin, V. S., Melnik, S. I., Kramarenko, S. S., & Ganganov, V. M. (2008). Analiz strukturi populyatsiy [Analysis of population structure]. MSAU Publishing House, Mykolayiv (in Ukrainian).
Szybiak, K., Bloszyk, J., Koralewska-Batura, E., & Goldyn, B. (2009). Small-scale distribution of wintering terrestrial snails in forest site: relation to habitat conditions. Polish Journal of Ecology, 57(3), 525-535.
Tokeshi, M. (1993). Species abundance patterns and community structure. Advances in Ecological Research, 24, 111-186.
Uvalieva, K. K. (1990). Nazemnye mollyuski Kazakhstana i sopredel'nykh territoriy [Terrestrial mollusks of Kazakhstan and adjacent territories]. Science, Alma-Ata (in Russian).
Weaver, K. F., Anderson, T., & Guralnick, R (2006). Combining phylogenetic and ecological niche modeling approaches to determine distribution and historical biogeography of Black Hills mountain snails (Oreohelicidae). Diversity and Distributions, 12(6), 756-766.
Whitney, M. E. (1937). Some observations on the reproductive cycle of a common land snail, Vallonia pulchella Mull.; influence of environmental factors. Proceedings of the Indiana Academy of Science, 47, 299-307.
Whitney, M. E. (1941). The hermaphrodite gland and germ cells of Vallonia pulchella. Michigan Academy of Science, Arts and Letters, 26, 311-338.
Whittaker, R. H. (1965). Dominance and diversity in land plant communities. Science, 147, 250-260.
Whittaker, R. H. (1972). Evolution and measurement of species diversity. Taxon, 21(2/3), 213-251. '
Yorkina, N., Maslikova K., Kunah, O., & Zhukov, O. (2018). Analysis of the spatial organization of Vallonia pulchella (Muller, 1774) ecological niche in Technosols (Nikopol manganese ore basin, Ukraine). Ecologica Monteneg-rina, 17, 29-45.
Zemoglyadchuk, K. V. (2005). Vidovoy sostav nazemnykh mollyuskov Berezin-skogo gosudarstvennogo biosfernogo zapovednika [Species composition of terrestrial mollusks of the Berezinsky state biosphere reserve]. Proceedings of the National Academy of Sciences of Belarus, Biological Series, 5(1), 87-90 (in Russian).
Zhukov, A. V., Zadorozhnaya G. O., Bets, T. U., & Lyadskaya, I. V. (2013). Dy-namika usadky dernovo-litogennyh g'runtiv na chervono-buryh glynah za sharamy [Dynamics of the shrinkage of the sod-lithogenic soils on the red-brown clays over the layers]. Scientific Herald of Chernivtsi University. Biology (Biological Systems), 5(3), 425-430 (in Ukrainian).
Zhukov, O. V., Zadorozhnaya, G. O., Maslikova, K. P., Andrusevich, K. V., & Lyadskaya, I. V. (2017). Ekologija tehnozemiv [Ecology of Technosoils]. Zhurfond, Dnipro (in Ukrainian).
Zipf, G. K. (1949). Human behavior and the principle of least effort. Addison-Wesley, Cambridge.
Zverkovskyi, V. M., Sytnyk, S. A., Lovynska, V. M., Kharytonov, M. M., & Mykolenko, S. Y. (2017). Remediation potential of forest-forming species in the reclamation planting. Ukrainian Journal of Ecology, 7(3), 64-72.
