1. Introduction
The soil seed bank represents a natural storage of plant seeds in soil, in which seeds may persist and remain viable for many years [
1]. Depending on the above-ground vegetation dynamics, soil seed banks contain different quantities of seeds. These quantities are largest in arable land, where numbers may exceed 10
5 seeds m
−2, and in grasslands (up to 10
4 seeds m
−2) [
2,
3]. Soil seed banks thus provide a “back-up” for future situations when environmental conditions become more suitable for germination (e.g., soil and canopy disturbance) and may prevent local extinction.
Plant species largely differ in the persistence of seed banks. Some species form only a transient seed bank, as seeds that do not germinate quickly do not usually survive more than a year [
2,
4]. Seeds of other species remain alive for a long time, even hundreds of years [
5], and form so-called persistent seed banks [
2]. Longevity correlates with seed morphology [
6,
7,
8,
9], taxonomic affiliation [
10,
11] and plant ecology [
12]; for example, arable species tend to have more persistent seeds compared to forest species [
4]. However, within evolutionary constraints of a species, seed persistence is not a fixed trait, and may vary with soil disturbance, type and pH, nutrient content in soil and (micro-) climatic conditions such as soil moisture and temperature and depth at which the seeds are buried [
13,
14,
15,
16,
17,
18,
19,
20,
21,
22]. Thus, the realized seed persistence is a combination of seed characteristics and the immediate environment [
12].
A variety of approaches have been applied to determine seed persistence in soil [
2]. Studies of changes in seed bank composition over time often rely on analyses of series of samples from sites with known time from the last disturbance [
5], which is used as a proxy of the minimum seed age. The limitations of these studies for understanding the dynamics of seed persistence in soil lie in the fact that the age determination of the retrieved seeds is very rough and that the initial seed cohort size is unknown. Another common approach includes the creation of artificial seed banks by either sowing the seeds in soil followed by repeatedly taking soil cores and viability testing the retrieved seeds [
23,
24] or by burying seed bags followed by periodic exhumation and viability testing [
2,
14,
20,
25]. With these approaches, the deterioration of the seed population can be followed, parametrized and used for making predictions. Such studies do not usually cover more than 5–6 years [
14,
21], and those that would span over longer periods of time usually do not collect data on yearly basis [
22,
26].
In addition to changes in viability, seeds in soil also undergo changes in external morphology [
27] caused by soil chemistry and activity of soil microorganisms. There is a lack of literature records that would provide a qualitative description of these changes, which may include changes in colour and loss of trichomes and other extremities or of the seed coat entirely. The little attention paid to these changes may result from the fact that the research has focused mainly on seed persistence per se and that morphological changes have been perceived as too descriptive and thus unimportant.
However, there are several ecological as well as practical reasons to have these changes described for an array of seed species. In the first place, many seeds undergo physical dormancy, i.e., they can germinate only if the seed coat is scarified, i.e., damaged, such as by soil microorganisms or abrasion during movement through soil. Knowing over the course of morphological changes then may become useful for predicting the germinability of the seeds. In the second place, knowledge of the morphological changes can be used in seed identification from soil samples, particularly in species in which these changes are substantial. Knowledge on changes in seed morphology may become relevant in connection with seed predation of exhumed seeds [
27,
28,
29] because predators are known to select seeds based on seed morphology [
30,
31]. The aim of this paper is therefore to describe and analyse (i) seed persistence and (ii) changes in morphological seed traits in 26 species of herb seeds in the course of 6–8 years of burial in artificial soil seed banks. Seed persistence was determined over up to 8 years spent in soil, 50 and 5% persistence times were predicted for each seed species based on logistic regression modelling of the proportion of persistent seeds, and changes in morphological traits that can be perceived by seed predators, such as mass, volume, density and form, were examined.
2. Materials and Methods
Seed material. Seeds of 26 species of herbs, mostly weeds, were introduced into artificial seed banks in 2005 and 2006 (
Table 1). The species of seeds were selected so they differed in morphology, taxonomic affiliation and presumed persistence in soil, and based on availability. Seed materials were sampled in Prague-Ruzyně (western Czech Republic) in a c. 10 × 10 km area surrounding Crop Research Institute (CRI) and centred at 50.09 N and 14.30 E. The seeds were harvested from mother plants by hand at full ripeness in July-October 2005 and October 2006 and stored dry at 5 °C until the burial. Furthermore, a subset of the seeds was stored at −20 °C as a control cohort (0 years). Seeds of 9 species were buried on 8 November 2005, and the remaining 17 species were buried on 24 October 2006 (
Table 1). Seed batches destined for burial for the time period of a particular length (seed cohorts) were prepared in a standard way. Approximately 10
4 seeds were mixed with soil dug from a 0.6 m depth and sieved through 0.05 mm mesh, placed in bags of nylon fabric and buried at a 20 cm depth under grassland on a ground of CRI. Mixing seeds with soil is important to prevent excessive degradation [
32], and fine sieving of the soil facilitates separation of the seeds after exhumation. Seed material was divided into 8 (2005) or 6 (2006) batches per species. Each seed batch packed in a bag with soil was buried separately and connected by a nylon cord with a label on the ground surface. Every year, one batch per seed species was exhumed, and recognizable seeds were separated from soil, dried at 25 °C and 40% r.h., and then stored at −20 °C for experimental use (a cohort). The last batches of seeds buried in 2005 were thus excavated in 2013 after 8 years of burial, and those buried in 2006 were excavated in 2012 after 6 years of burial, so there were 8 or 6 burial cohorts per species, respectively. For
Crepis biennis 2011 excavation and for
Plantago lanceolata 2012 excavation, no seeds were available due to complete deterioration of the respective batches.
Changes in seed morphology with duration of burial. Prior to measuring changes in seed morphology with duration of burial, subsamples of seeds from each cohort were cleaned from fine soil particles in an ultrasound cleaner (Sonorex RK 31, Bandelin electronic, Berlin, Germany), submerged in water for 2 min, and dried in an oven for 24 hr at 75 °C. The following measurements were made on seeds from each cohort:
Seed mass (M)—For each species and cohort, the average seed mass was determined based on five batches of 20 seeds using analytical balances (CP225D-0CE, Sartorius AG, Göttingen, Germany) with a precision of 0.00001 g.
Seed dimensions—following [
19], five seeds per cohort were measured by digital scales NTD12P-15CX (Mitutoyo Corp., Kawasaki, Japan) with a precision of 0.01 mm. The following dimensions were measured: length (
L; the longest dimension of the seed), width (
W; the longest dimension perpendicular to
L within the same plane), and height (
H; the longest dimension perpendicular to the plane of
L and
W). Based on these measurements, the following metrics of the seed form were calculated for each of the five seeds per combination of species and cohort, except for density, which could only be estimated as a mean value:
Volume—calculated as
[
34]
Density—calculated as
[
34]; only mean values of mass and volume were used, as these were measured on different seeds.
Shape—according to [
19], the seed shape can be expressed as a dimensionless measure:
, where x represents a division of either
L,
W, and
H through
L and
as their mean, and
n is 3.
Vs ranges from 0 for perfectly spherical seeds to 0.2 shaped as a thin disc or spindles.
Flatness Index—calculated as
[
34]. It ranges from 1 for a complete sphere to greater values for plane or spindle like shaped seeds.
Eccentricity Index—calculated as
[
34]. It ranges from 1 for round seeds to values greater than 2 for spindle like seeds.
Seed persistence in soil. Seed persistence was measured on another subsample of seeds. We used the so-called imbibed seed crush test (ICT), which was found to be reliable for estimation of the true viability of weed seeds [
35,
36,
37]. This test was performed in order to determine the proportion of fresh [
38], i.e., persistent and viable seeds. In these seeds the seed coat does not collapse when crushed with the tips of a pair of forceps, or have apparent and intact cotyledons or embryos when the seed coat is broken. In other cases, the seeds were considered dead. In this test, 20 seeds per cohort tested, and were left imbibe for 3 days in laboratory conditions (20 °C, 12 h light: 12 h dark) prior testing.
Statistical analysis. Prior to analysis, the duration of burial was expressed as a continuous numerical vector, year after burial (
YAB), and used as the explanatory variable in most of the analyses. Frozen (control) seeds were given age 0, those exhumed after one year were given age 1, etc., so the maximum age was 8 for seeds buried in 2005 and exhumed in 2013 and age 6 for those buried in 2006 and exhumed in 2012.
YAB 1 then denotes a cohort of seeds that were buried for 1 year and so on. The analyses were performed in R version 3.3.1 [
39]. The variation in morphological seed traits with
YAB was tested within seed species by regression methods [
40]. For comparative reasons across species of seeds, the values of
M,
V and
D were converted to relative scale against the mean values of the respective trait for control seeds and denoted
rM,
rV and
rD (i.e., relative mass, relative volume and relative density). The three seed form indexes were not relativized because they are inherently dimensionless. For each species of seeds, the
rM,
rV and
rD were regressed against
YAB by fitting three different models were used to fit:
- (a)
Linear: , where y is the respective trait on relative scale, a is the intercept, and b is the slope of the linear regression line;
- (b)
3-parameter logistic: , where y is the respective trait on relative scale, A is the upper asymptote, C is the scale parameter on the x-axis, and D is an inflexion point of the curve;
- (c)
4-parameter logistic: , where y is the respective trait on relative scale, A is the upper asymptote, B is the lower asymptote, C is the scale parameter on the x-axis, and D is an inflexion point of the curve.
The non-linear curves were fitted as the self-starting functions SSlogis and SSfpl for 3- and 4-parameter logistic functions (nls package of R). The explanatory power of the three models was compared based on the Akaike Information Criterion (AIC) [
40], and the one with the lowest AIC value (∆AIC = 2 as threshold) was chosen as best and presented. Only the linear model was fitted to
Shape,
FI and
EI based on a priori visual data inspection. In case the diagnostic graphs [
40] suggested that outlying
YAB values could influence the parameter estimates, these were removed from model fitting. To better interpret the general pattern in seed form variation with
YAB, we performed principal component analysis by implementing the
prcomp function from stats package of R and including
rM,
rV,
rD,
Shape,
FI and
EI as variates, and the resulting principal scores for the first and second axes (
PC1 and
PC2) were used for interpretation of changes in the proportion of persistent seeds [
9].
The variation in the proportion of persistent seeds was assessed based on the logistic regression:
[
41], and linearized by implementing generalized linear modelling, with
a being the intercept and
b the slope of the regression, using binomial distribution of errors and logit link function (GLM-b) [
17]. This approach takes the nature of the data (proportions of viable seeds in case of germination test and binary response variable: 0—dead, 1—live into account without a need for transformation [
40] and allows for making predictions [
17]. The parameter
x is either
YAB,
PC1,
PC2 and
(i.e., mean relative seed mass per
YAB cohorts) and the models were compared based on the AIC as above. The seed persistence in the soil was estimated for each species of seed based on the parameters of the model with
YAB, predicting the 50 and 5% persistence time in
YAB (
PT50 and
PT05, respectively), i.e., the time after burial when 50 and 5% of the initial seed population were still viable. To identify which morphological traits might have affected the predicted seed persistence across species, we regressed
PT50 and
PT05 separately against all traits measured on control seeds against
rM,
rV and
rD for each separate
YAB cohort and against slopes of change in
Shape,
FI and
EI across
YAB.
Seed longevity literature survey. Seed longevity data for some of the study species can be excerpted from the literature. The monograph [
2] was used as a starting point for the survey. We used only seed persistence data originating from burial experiments (method coded 1–3 after work [
2]), which were available for only 14 out of 26 species of seeds included in this study. From the literature data, we calculated two indexes, the longevity index (
LI) and seed persistence index (
SPI).
LI was calculated as:
, where
Rt,
Rsp and
Rlp are the proportion of records classifying the species as transient (persisting less than 1 year), the proportion of records classifying the species as short persistent (persistent for more than 1 year but less than 5 years), and the proportion of records classifying the species as long persistent (persistent for more than 5 years), respectively [
4].
SPI was calculated as:
, where
T,
SP and
LP represent the number of records reporting transient, short persistence and long persistence, respectively [
8]. The values of
LI and
SPI with the number of literature records used for calculations are listed in
Table 1.
Seed mass (as
rM) decreased with
YAB in 25 out of 26 species of seeds available. The only exception was
H. niger, the seeds of which did not change their
rM with
YAB (
Figure 1a). Logistic curves appeared to fit the decline in
rM better than the simple linear fit in 17 cases, suggesting that in the majority of species, the rate of mass reduction is slow during the first years after burial, then accelerates and finally slows down again. The 4-parameter logistic curve that included a lower asymptote described the data best in 10 cases, and the 3-parameter curve without the lower asymptote curve did so in seven cases. Examples are in
Supplementary Materials. Seed volume (as
rV) declined with
YAB in 14 species, of which the linear decrease was found in 10 and logistic in four species, respectively (
Figure 1b). The surprising observation that relative seed mass or volume increased after burial in some species of seeds (
Figure 1a,b) has to be viewed as an artefact of natural variability among the seeds. Seed density (as
rD) was even less variable with
YAB in terms of the number of significant changes (in 11 species only), but the changes were both positive (1 case) and negative (10 cases), suggesting that the mechanisms that lie behind the variation in morphological traits over time spent in soil are diverse across species. Three species of seeds showed a significant linear change, and eight species showed a significant logistic course of change in
rD with
YAB (
Figure 1c).
3. Results
The interspecific interrelationship of the three morphological traits and temporal dynamics of their change across
YAB is exemplified in
Figure 2 using fitted values from the above models for
YAB 1,
5 and
8. After one year in soil, there was little variation in
rM,
rV and hence
rD (
Figure 2a). With the course of time spent in soil, it becomes evident that the rate of change in
rM and
rV and hence in
rD becomes more variable across species (
Figure 2b,c)—species distributed along the line representing slope = 1 through the origin are those in which both mass and volume decreased in similar rate, species below the line are those in which the volume decreased more rapidly than mass and species above the line are those in which seed mass declined more rapidly than volume.
The shape of seeds remained similar overall, as the three indexes of the seed form varied with
YAB in only a few species (
Figure 3a–c). Principal component analysis, however, revealed patterns in seed form changes across seed species and
YAB (
Figure 4). The PC1 axis explained 48.4% of the variance and presented the change in external seed form as
Shape,
FI and
EI correlated positively and
rV negatively with this axis. This suggests that seeds that change their shape also lose their volume as they either shrink and become flatter or lose their extremities. The PC2 axis explained 27.3% of the variance and was well correlated positively with
rM and
rD, i.e., represented by seeds that did not change their volume but changed their mass and hence density. The individual variation in all morphological traits with
YAB for each of the seed species can be found in the
Appendix S1: panels A–F.
Alternative models for the proportion of persistent seeds provided mixed results. While the two models based on
PC1 and
PC2 together improved the explanatory power of the variation compared to the ICT model in two cases only (
PC1:
L. cardiaca, AIC = 185.28;
PC2:
C. album, AIC = 63.86), models based on
rM explained the change in the proportion of persistent seeds better than did the initial models in 10 cases (
Table 2) suggesting that the decline in the proportion of persistent seeds may be associated with the decline in seed mass more tightly than solely to
YAB, especially when large or inconsistent scatter in
rM across
YAB appeared. The change in the proportion of persistent seeds with
YAB for each of the individual species can be found in the
Appendix S1: panels G, and with relative seed mass in the
Appendix S1: panels H.
The persistence times of particular seed species in soil predicted based on the ICT are shown in
Table 2. In some species, such as
A. sagittata, the
PT50 was predicted to be as low as 1 year, while in
G. urbanum, the
PT50 reached 10 years. The
PT05 ranged from 2 (
C. glaucum) to more than 20 years (
U. dioica and
P. major). As all seeds were found to be viable in
T. arvense, the persistence in soil could not be estimated for this seed based on our data.
Across the species, predicted persistence was not related to
M,
V,
D nor any of the three indexes for seed forms, using data for control seeds (results not shown). On the other hand, the
PT50 was negatively related to the
rM of buried seeds (
Figure 5a). The relationship of the
rM with
PT50 for
YAB 1 and
YAB 2 cohorts was much steeper compared to the remaining
YAB cohorts, suggesting that the seed species losing their mass in the first two years of burial survive in soil relatively shorter. The results were also similar for
PT05, with the exception of the insignificant relationship with
YAB 1 (not shown). The relative volume of neither cohort significantly explained the
PT50 nor
PT05 (not shown), and neither did the relative density of
YAB 1 and
YAB 2. In contrast, seed persistence was negatively affected by the decrease in relative density in cohorts of
YAB 3–8 (
Figure 5b), suggesting that seeds that lose relatively more mass than volume suffer from shorter lifespan in soil. Change in seed form indexes seemed to affect the predicted persistence (
Figure 5c), but the effect was due to the outlying
G. urbanum seed, which greatly changed its shape as a result losing seed coat and its appendages, so the trends were no longer significant when this species was removed from the analyses.
4. Discussion
The seed characteristics investigated in this study involve traits that directly affect the fitness of plant species as well as traits that are most likely relevant only with respect to the perception and attractiveness of seeds to seed predators. The seeds were generally losing viability gradually over time when buried in the soil, and this loss of viability was accompanied by more or less noticeable morphological changes. Of these changes, the loss of mass and volume were the most prominent, although changes in other traits (seed density and shape) appeared to be significant in some species of seeds as well.
Seed persistence. Seed viability is known to non-linearly decline with time spent in soil as seeds decay or age [
14,
17,
25,
42], and this was true for the vast majority of species included in this study. Exponential decay rates were previously used for fitting decline in seed viability and for estimating the annual decay rate [
14,
43]. However, the visual inspection of the plots in the abovementioned studies clearly shows that the exponential curve was not capable of catching the uneven decline rate over time (see, e.g.,
Figure 1 in Conn et al. [
43], plots for
C. bursa-pastoris or
T. inodorum, among others), resulting in the prediction of unrealistic annual decay rates. In this paper, we took advantage of using GLM with a binomial distribution of errors, which is a linearized form of a logistic regression, for describing the course of decline in seed viability, approximated by using the change in the proportion of persistent seeds in artificial seed banks over time. The method for analysing and predicting the persistence of seeds from soil seed banks at the population level was previously used, e.g., in work of Pakeman et al. [
17], and we encourage others to do so as well.
We found this method to be robust, well grounded in mathematical theory, and easy to perform and interpret, and more importantly, it can be used for modelling seed persistence in soil. The greatest advantage of this approach is that one can describe the dynamics of seed persistence in time (or with any other variable that may have influenced the persistence, such as mass also used in this study) in a realistic way and estimate the persistence times. Using this method more frequently, the estimation of seed persistence in soil would approximate the reality closer, and would be more comparable among studies. The prerequisite for using this method is that a proportion of a seed population at each time point must be known, as estimated by a subset of the retrieved seeds, as in this study. Another advantage is that by using the appropriate distribution and link function, no data transformation is needed. Other researchers used various logistic regression models for analysing the decay/survival of seeds in soil [
42,
44] or complicated probit-based models [
25]. Depending on the nature of their data (either Poisson or binomial distribution of errors), various data transformations were applied to approach counts (Poisson) or proportions (binomial) to Gaussian distribution, allowing for fitting the chosen model. Although fitting some of these models may provide seed persistence curves and hence help with making predictions, we argue that GLM with binomial error structure shall be the preferred method in such studies due to its simplicity and accuracy, and this should be taken into consideration already when seed persistence data are collected.
We consider that wider application of logistic regression models in seed bank persistence studies would greatly increase our understanding of the dynamics of seed banks and that this approach represents a major step forward from classifying species into groups according to the detected persistence of their seeds, e.g., transient (< 1 year), short persistent (1–5 years) and long persistent (> 5 years) [
2]. Such classification might erroneously suggest that seed persistence in soil is a discrete process, although the logistic approach has demonstrated that losing viability is continuous process and that a minor proportion of seed population may persist for a long time. For example, it can be predicted for
C. trachelium, one of the least persistent species in our study, that 2.1 seeds out of a hypothetical 10
5 initial seeds will survive for five years in the soil, using the parameters describing the course of change in the proportion of persistent seeds estimated here, and their detection (and classification in the respective persistence class) depends on the sampling method, sample sorting precision and, of course, coincidence of finding viable seeds in soil samples. Therefore, literature data on seed persistence are extremely variable for each species (see, e.g., database of Thompson et al. [
2]), and it is not rare that a species is classified as transient in one study and long persistent in another (e.g.,
A. retroflexus used also in this study according to the database of Thompson et al. [
2]). There are obvious ecological reasons that may have caused such differences, including interpopulation or geographical variation, soil conditions and microflora [
9,
12,
15,
17,
45], but such variation can also be largely attributed to the methods used for data collection (burial vs. emergence, soil depth) and evaluation. It is a well-established fact that seeds tend to germinate earlier and thus remain in a seed bank for less time when buried at shallow depths compared to deeper burial [
25]. It also seems that seed emergence data are less reliable than seed burial experiments basically because species with larger seed production are more likely to be detected in soil for a longer period of time [
46]. To address the variability among studies, several authors developed indexes [
4,
8,
19] that consider the relative frequency of studies reporting a class of seed persistence and unite the known persistence data in one number. However, these indexes are inherently biased when different methods of seed persistence evaluation are combined for calculation [
46] and therefore ecologically meaningless, despite their frequent use in ecological papers. For the sake of demonstration, we calculated the longevity index [
4,
19] and seed persistence index [
8] for the 14 species for which data were available using only the burial data from the database of Thompson et al. [
2]. The variation in the indexes values were in fact minor among species (
Table 1), and all of them could be classified as long persistent if one would wish to do so, despite the variation in persistence times found in our study. In fact, papers that would provide reliably calculated predictions of seed persistence times are scarce, making difficult any comparison with persistence times
PT50 and
PT05 predicted in this study. For the species included in this study, we found only one reference [
17], and interestingly, by using binomial GLM, the estimated
PT50 for
P. lanceolata (2.35 ± 0.43) was nearly identical to that from our study (2.2 ± 0.18). Others provided persistence times based on an exponential decay formula, but the resulting values are not comparable due to lack of fit (see above). We attempted to find patterns in the variation in persistence times and ecological explanation of these patterns by using morphological data of the control seeds as well as data describing the morphology change of seeds over time spent in soil. We found that predicted persistence times were unrelated to any of the morphological traits of the control seeds available, which contradicts the previously published literature data (e.g. [
6,
7,
19,
45,
47]). Nevertheless, we found that a relative change in seed mass affected the predicted persistence. Seed species that lose mass rapidly already in the first two years of burial survive in soil relatively shorter. Seeds that lose relatively more mass than volume, i.e., seeds that have a rather strong seed coat, would suffer from a shorter life span in the soil as the embryo deteriorates. Additionally, relative seed mass explained the decline in the proportion of persistent seeds better than
YAB in a significant portion of the species studied, mainly in those in which an increase in the proportion of persistent seeds occurred after a period of decline. Such reversal periodically appears in the literature [
14,
21,
22,
26]) and is sometimes explained as a result of microsite variation in soil conditions favourable or disfavouring seed persistence. Additionally, the persistent times estimated in this study might be overestimating those that would occur in the field since these were based on identifiable seeds only and neglecting those which deteriorated completely. As seed mass generally declines with
YAB non-linearly, using seed mass as an explanatory variable for the decline in the proportion of persistent seeds partly removed this scatter from the data. Congeneric species included in this study were expected to have similar persistence times because of common evolution and thus presumed taxonomic constraints on seed biology and similar life histories [
12]. Contrary to this expectation, the variation among congeners was quite large in most cases; only the species of
Amaranthus had rather similarly short persistence times. It is difficult to provide an explanation for such a scatter based on a limited number of species per genus (or family), and more likely, functional traits of plants and their habitat requirements might help to disentangle the ecological significance, if it exists, of such a striking scatter.
Seed morphology. The limited availability of descriptions of changes in external morphology and seed mass through time spent in soil is evident because these traits have not seemed to be important for species’ ecology or economic importance. They may, however, have vital importance in determining the consumption rate of seed predators [
48] and for identifying seeds from seed banks. We provide evidence that over the study period, the seed mass, volume and density largely become altered. There are several possible causes that result in changes in seed mass, volume, density and form through seed aging and decay. These relate to deterioration of the external structures, decay of the embryo and shrinking of the seed. As changes in specific descriptors of seed morphology were often correlated for a particular species of seed, the nature of the correlations provides insight into the manner of morphological changes the seeds underwent during soil burial. Significant change in seed shape indicators with
YAB was caused by the reduction in volume. This can be a consequence of shrinking soft seeds due to death and deterioration of the embryo (e.g.,
C. trachelium), so the density remained unchanged. In a special case of
G. urbanum, the seed coat of the control seeds is soft and equipped with a hook. Both structures deteriorate over time, which changes the overall seed shape and even increases seed density. Our results suggest that the seed coat does not have protective function in this species since
G. urbanum was also one of the longest persistent species in our study. Many seeds did not change their volume over time in soil but mass (e.g.,
A. powellii), which indicates that these seeds are hard-coated with embryos dying and deteriorating. The fact that in some species of seeds the relative mass and volume was observed to “increase” after burial has to be viewed as a result of natural variability in seed size; alternatively it may partly be an artefact that the smallest seeds deteriorated first, so these were not available for taking morphological measurements.
The morphological changes and changes in seed mass over time spent in soil can also be well described by logistic regression in many species, and it seems plausible that this approach will find wide application in seed ecology studies. In addition to seed persistence [
17,
42,
49] and morphology (this study), logistic models were recently applied for describing germination response to environmental conditions [
50]. Thus, we propose using these models in situations in which the response of this type of seed to a focal variable can be expected.