New Research In
Physical Sciences
Social Sciences
Featured Portals
Articles by Topic
Biological Sciences
Featured Portals
Articles by Topic
 Agricultural Sciences
 Anthropology
 Applied Biological Sciences
 Biochemistry
 Biophysics and Computational Biology
 Cell Biology
 Developmental Biology
 Ecology
 Environmental Sciences
 Evolution
 Genetics
 Immunology and Inflammation
 Medical Sciences
 Microbiology
 Neuroscience
 Pharmacology
 Physiology
 Plant Biology
 Population Biology
 Psychological and Cognitive Sciences
 Sustainability Science
 Systems Biology
Seasonality, density dependence, and population cycles in Hokkaido voles

Edited by Thomas W. Schoener, University of California, Davis, CA, and approved July 11, 2003 (received for review September 2, 2002)
Abstract
Voles and lemmings show extensive variation in population dynamics regulated across and within species. In an attempt to develop and test generic hypotheses explaining these differences, we studied 84 populations of the graysided vole (Clethrionomys rufocanus) in Hokkaido, Japan. We show that these populations are limited by a combination of densityindependent factors (such as climate) and densitydependent processes (such as specialist predators). We show that densitydependent regulation primarily occurs in winter months, so that populations experiencing longer winters tend to have a stronger delayed densitydependence and, as a result, exhibit regular density cycles. Altogether, we demonstrate that seasonality plays a key role in determining whether a vole population is cyclic or not.
Along controversy over the issue of densitydependent versus independent population regulation has led to the conclusion that both factors are important for understanding population fluctuations (1–8). It is, however, less clear how such densitydependent and independent factors interact with each other in shaping the dynamic pattern of populations across a larger part of the species range. In an attempt to disentangle these issues, we analyze a set of 30year seasonal (spring and fall) time series from 84 populations of the graysided vole [Clethrionomys rufocanus (Sundevall, 1846)] from Hokkaido, Japan (Fig. 1A) (9, 10). To investigate the role of seasonality in the generation of population cycles, we decompose the annual (falltofall) density dependence, as well as the densityindependent stochasticity into their seasonal components. The added detail provided by pinpointing the seasonal arena of population regulation (see supporting information on the PNAS web site, www.pnas.org) provides us with a better basis for suggesting and evaluating hypotheses about the biological mechanisms that cause density dependence, stochasticity, and population fluctuations.
A perennial problem in the study of population dynamics has been the relative lack of extensive and accurate data. In this study we attempt to reach more accurate conclusions than are usually possible by two means. First, we use comparative timeseries data from a large number of very similar populations. Second, we address the very substantial problem of biased and imprecise measures of population density through use of a statespace modeling approach (11–13), where timeseries observations are related to unobserved “states” of the real population through a probabilistic observation model accounting for sampling variation. Our study confirms and extends an earlier study of ours (14); whereas the earlier study used only fall data (and as a result could cover the entire island of Hokkaido), the present study used both fall and spring data. The greater detail of the data used in this paper makes a much more detailed analysis of the seasonal structure possible; a picture we assume, based on earlier studies (14), applies to the entire island of Hokkaido.
The Study System, Its Seasonal Structure, the Sample Sites, and the Data
Hokkaido is the northernmost island (41° 24′–45° 31′ N, 139° 46′–145° 49′ E) of Japan and covers 78,073 km^{2}. It neighbors the Asian continent and is surrounded by the Sea of Okhotsk, the Pacific Ocean, and the Sea of Japan. A southern warm current prevails in the Sea of Japan along the western shore of Hokkaido, whereas northern cold currents prevail in the Pacific Ocean along the northern and eastern shores. Mountain ridges run north–south through the middle of the island; in the southwestern part there is another mountain ridge. Except for the temperate deciduous forests of the southwestern peninsula, most of the natural forests in Hokkaido are regarded as the transition between the temperate and the subarctic zones (15). The dominant tree genera are Abies, Acer, Betula, Picea, and Quercus. Hokkaido represents the easternmost part of the graysided voles' distribution, which extends as far as Fennoscandia in the west (16, 17).
The lengths of the seasons vary from one year to the next. Information on this variability is not available. A related measure, an index of vegetation growth called the warmth index (WI) (see supporting information) is, however, available for a large number of meteorological sites, but not for all trapping sites used in our study. This index is used as a basis for developing a measure of relative length of the winter (τ_{w}), being 1 – (relative length of the summer) (for details, see supporting information).
The graysided vole (C. rufocanus) represents a pest on plantations of larch (Larix leptolepis) and todofir (Abies sachalinensis). Since 1954, the Forestry Agency of the Japanese Government has carried out censuses of vole populations for management purposes in forests all over Hokkaido. The forests managed by the Forestry Agency cover 28,400 km^{2} [21,500 km^{2} natural forests and 6,900 km^{2} planted forests (1992 figures)]. These forests were in 1992 managed by 76 district offices, which were further divided into several ranger offices, giving a total of 433 ranger offices. The individual ranger office, which also represented our basic unit of analysis, carried out the censuses; personnel at the ranger offices were regularly trained, including in species identification, for performing the censuses (16).
Rodent censuses were carried out twice a year [spring (May or June) and fall (September or October)] on a 0.5hectare (50 × 100 m) plot [the latter being used for the annual analysis (14, 17, 19)]. At each plot, 50 snap traps were set in 10 × 10m grids for five or three nights. For each site, one to three separate grids were monitored in fixed preselected natural forests. The census grids were occasionally relocated within the local area. Altogether, 84 time series for both spring and fall covering 30 years (1963–1992) in the central and northernmost part of Hokkaido (the Asahikawa Regional Office, Forestry Agency of the Japanese Government) are used; this region is the one for which cycles occur (10) and for which we have sufficient data. We have grouped these populations according to topographic characteristics in our previous studies (19). Trapping efforts were represented by multiplying the number of traps by the number of nights and the number of census plots (grids). For instance, when threenight censuses were carried out at three grids in a ranger office, the trapping efforts by the ranger office at the census was 450 (= 50 × 3 × 3). The springtofall season (summer) corresponds to the period of population growth (because of reproduction), whereas the falltospring season (winter) corresponds to the period of population decline [winter reproduction is negligible (refs. 17 and 20–25; see supporting information)].
Reliable data are available on food habits of voles (16, 20). C. rufocanus of Hokkaido prefer green plants as food (20). Leaves of bamboo grasses (Sasa senanesis, Sasa kurilensis, and Sasa nipponica) constitute a large part of the winter diet, and their shoots occupy a considerable part of the summer diet (although summer diet is much more varied than winter diet). Thus, bamboo leaves influence survival during the winter and bamboo shoots influence reproduction during the summer (most likely primarily during the spring because other plants do not grow as well at this time of the year). Bamboo grasses may be particularly important during winter when other plants wither, even though their leaves do seem to be nutritionally poor during the winter. Even though there are many predators such as snakes, birds, and mammals in Hokkaido (21), their activity under the snow cover is often restricted. An exception is the least weasel Mustela nivalis, which is a highly efficient predator under the snow cover.
The graysided vole is the most common rodent species in Hokkaido (20, 25). Three other microtine and murine rodent species are recorded as an integral part of the census: Clethrionomys rutilus (Pallas, 1779), Apodemus speciosus (Temminck, 1844), and Apodemus argenteus (Temminck, 1844). In addition, Clethrionomys rex Imaizumi, 1971, Apodemus peninsulae (Thomas, 1907), and shrews (Sorex spp.) are occasionally caught. Personnel of the Forestry Agency identified the specimens. Clethrionomys is easily distinguished from other genera. Distinguishing C. rufocanus from its congeners may, however, be difficult, although the abundance of C. rufocanus and C. rutilus was separately reported. The numbers for the other Clethrionomys species is, however, very low in Hokkaido (20, 26), reducing the impact of any misclassification.
The Ecological Model
Let spring and fall abundances in year t (logtransformed and centered around their means) be denoted x_{t} and y_{t}, respectively. The seasonal models for the net winter population growth rate, R_{w}_{t} (= x_{t} – y_{t}_{–1}), and the net summer population growth rate, R_{s}_{t} (= y_{t} – x_{t}), are given as (27) [1a] [1b] where ε_{w}_{t} is the timeindependent process noise during the winter and ε_{s}_{t} is the corresponding summer component . The parameters a_{w} and a_{s} define the seasonal densitydependent structure for the winter and summer growth rates, respectively. The annual net growth rate is obtained by rewriting Eqs. 1 and 2 as a model in y only (27): [2] where is the annual noise term (see supporting information). The two last autoregressive terms in Eq. 2 are typically negligible (as is also the case in this study), rendering a secondorder process an appropriate approximation (6, 8, 25): [3] where ε_{t} is the annual noise term [ε_{t} ∼ N(0, ρ^{2})] and the parameters α_{1} and α_{2} define the direct and delayed annual density dependencies; for this system to correspond to a stationary process, the two autoregressive parameters (1 + α_{1}) and α_{2} must fall within the triangle defined by (–2, –1), (0, 1), and (2, –1) (28).
The link between the annual and the seasonal parameters are given as α_{1} = a_{w1} + a_{s1} + a_{w2} + a_{s2} + a_{w1}a_{s1} and α_{2} = a_{w3} + a_{s3} + a_{w4} + a_{a4} + a_{w1}a_{s3} + a_{s1}a_{w3} – a_{w2}a_{s2}. Assuming a whitenoise term (i.e., no delayed noise in the noise process of Eq. 3), ρ^{2}, the annual noise variance may be expressed by given as a function (see supporting information).
Linking the Data to the Ecological Model
Sampling variance is known to bias the estimates of, and invalidate tests for, density dependence (5, 29); the degree of direct density dependence is typically overestimated when ignoring sampling variance. Despite this, sampling variance is generally ignored (25) in empirical studies, partly because of availability of data and methodological shortcomings. Sampling variance may be accounted for by adopting a statespace modeling approach where timeseries observations are related to unobserved “states” of the real population through a stochastic observation model (11–13). Let the number of voles caught during the spring and fall session at a given site in year t be given by z_{s,}_{t} and z_{f,}_{t}, respectively. Similarly, let the trapping effort at a given site in the spring and the fall of year t be given by T_{s,}_{t} and T_{f,}_{t}, respectively. Our core data consist of the total number of graysided voles caught (z_{s,}_{t} and z_{f,}_{t}) and the corresponding trapping effort (T_{s,}_{t} and T_{f,}_{t}). The expected number of voles caught during a trapping session will be related to the abundance of voles at the site (y_{t} and x_{t}, logabundances scaled around their respective mean values, μ_{y} and μ_{x}) and on the trapping effort. As a first approximation, the ecological process model as described by Eqs. 1a and 1b or Eq. 3 is assumed to be linked to the observed number of voles caught through a Poisson model. For the fall samples this will then be given as: P(z_{f,}_{t} = zy_{t}) = exp(–λ_{f,}_{t})λ^{z}_{f,}_{t}/z!, where the Poisson mean (assumed to be proportional to the trapping effort and the population abundance) is defined as λ_{f,}t = qT_{f,}_{t} exp(μ_{y} + y_{f,}_{t}). The trapping effort T_{f,}_{t} is a known parameter. The parameter μ represents the mean of the logtransformed “true” population abundances. Because of the nature of the data, the average level of the true population (μ) is confounded by the “trappability” (i.e., the average probability for any present individual to be trapped) in the given environment q. However, although μ cannot be directly estimated, we can estimate the product of the average abundance level and the “trappability,” defined as exp(γ) = q exp(μ). This yields the following reparameterization (for the fall; and equivalently for the spring): λ_{f,}_{t} = T_{f,}_{t} exp(γ_{f} + y_{f,}_{t}).
By using the software package bugs (Bayesian inference using Gibbs sampling; ref. 30), a Bayesian framework was adopted to estimate parameters in the statespace model. The resulting posterior distributions of parameters from the seasonal model are shown for two typical populations (Fig. 2). Only vague (i.e., essentially flat) prior distributions were used. For further details, see supporting information. Posterior means are used as parameter estimates and used in further analyses. Seasonal parameter estimates being judged as statistically appropriate were obtained from only 74 sites (see supporting information). All 84 time series are, however, used for obtaining the main results summarized in Table 1, as well as those depicted in Figs. 1 B and C. Notice that the annual density dependencies given by Eq. 3 are both estimated to be stronger when sampling error is not accounted for (Fig. 3A).
Results
The obtained parameter estimates (Table 1 and Figs. 1B and 2) suggest that a fourthorder process (with the existence of some long lags) underlies the dynamics in at least one of the seasons (see Eqs. 1 and 2). The direct falltospring effect (a_{w1}) is clearly stronger than the direct springtofall effect (a_{s1}), suggesting that winter processes are more important than summer processes in determining the annual direct density dependence. These, as well as falltofall effects (a_{s2}) and springtospring effects (a_{w2}) all contribute to the annual direct density dependence (see Eq. 2 and supporting information), a result consistent with earlier findings both for Hokkaido (see refs. 25, 31, and 32) and the same species in northern Finland (refs. 27 and 33; also see refs. 34–38). As seen from Fig. 3B, we are able to predict both the direct and delayed annual density dependencies from the seasonal components. Assuming the model given by Eq. 2, the seasonal components explain 70% and 79% of the (deterministic) amongpopulation variability in the direct (α_{1}) and delayed (α_{2}) annual density dependencies, respectively. A major part of the amongpopulations variation in the level of densityindependent stochasticity (i.e., amongpopulation variation in ρ) is, however, accounted for by variation in the level of summer stochasticity (i.e., amongpopulation variation in ρ_{s}; Fig. 4). Variation in degree of summer stochasticity accounts for 32.0% of the total variation in annual stochasticity, whereas the winter stochasticity accounts for 11.7%. If summer stochasticity already is accounted for, the winter stochasticity only accounts for 1.2% of the remaining variation. In the multiple regression of annual stochasticity on the summer and winter stochasticities, the effect of the winter stochasticity is not statistically significantly different from zero.
It should finally be noticed that the falltofall annual densitydependent structure corresponds to the springtospring densitydependent structure (see supporting information). This close correspondence clearly strengthens the validity of the obtained results.
Discussion
The Seasonal Structure of the Population Dynamics. The summer is the main reproductive season for the graysided vole; essentially no winter reproduction is observed (see supporting information and refs. 16 and 20). Hence, strong, direct density dependence during the winter must necessarily involve winter survival. Mechanistically, such density dependence may be generated through a functional response of predators or through the limited amount of food (such as bamboo; Sasa spp.) produced during the preceding summer. On the other hand, the strong influence of summer stochasticity (Fig. 4) may stem from a variety of summer events (such as plant production and generalist predation), all of which may be sensitive to erratic temporal and spatial variation in climatic conditions. Notice, however, that when interpreting the statistical parameter a_{w1}, given a mechanistic process model, this will typically involve direct density dependence of the vole during the winter, in addition to other processes from other seasons. Hence, biological interpretation of these parameters is difficult in the absence of a detailed mechanistic process model. The same applies to the other seasonal statistical parameters.
If we are to explain the observed phenomenological lagstructure, involving four lags (Table 1), in terms of a process model that only allow interactions within a season, we would have to postulate four underlying variables. One of these must be the density of voles. The character of the remaining three variables is unclear. We suggest three possibilities: (i) voles, winter food (bamboo grass), and a specialist predator (the least weasel, M. nivalis) with two age classes (the two age classes of the predator are reasonable on this time scale; weasels require one year to mature); (ii) voles, bamboo grass, specialist predator, and additional members of the smallmammal community that either compete directly with the graysided vole or induce apparent competition (39) through some common predator; and, finally, (iii) an interacting web consisting of voles, food plants (e.g., bamboo grass), weasels, and generalist predators. A range of other hypotheses (involving, e.g., intrinsic properties of the vole populations) are certainly possible; however, the three suggested hypotheses are more plausible based on general naturalhistory insights (see ref. 40 for details). We know, for instance, that the vole specialist (41–43) least weasel is found all across the island (41, 44) and that it may have a particularly strong effect during the winter, because it is able to hunt in the subnivean space where the voles are protected from other predators. Bamboo grass is common all across Hokkaido (45); their new shoots in the spring are known to affect reproduction during the summer, and leaves of bamboograss (laid down by the snow) are known to affect survival during the winter. It is difficult to explore these alternatives within a timeseries setting, because the number of possible interactions in a fourvariable model leads to a large number of parameters. Thus, to further test such process hypotheses, one would need to search for more direct evidence of interacting variables (or perform appropriate experiments; see ref. 40 for further discussion).
The striking difference between population regulation during summer and winter suggests that seasonality may hold a clue to understanding geographic differences in population dynamics. Notice, for instance, that the western populations (Group 1; see Fig. 1 A) tend to fall in the upper right part of the triangle of Fig. 1B (corresponding to stable populations), whereas eastern and the mountainous populations tend to fall further down toward the central part of the triangle (corresponding to periodically fluctuating populations). Because of warm currents along the western coast of Hokkaido, the Group 1 populations are exposed to the longer summers (17). On the other hand, colder currents along the eastern coast expose the Group 2 populations to longer winters (17). The more mountainous Group 5 populations are also exposed to long winters (17). On this basis, it is reasonable to hypothesize that the observed systematic differences in the population dynamics of these groups are influenced by differences in season length, a view that is substantiated by observing the relation between the annual density dependencies and the length of the winter, particularly so when using the pooled estimates (Fig. 1). What we in effect might see is a detailed example of more cyclic dynamics being observed toward the north (18), where the summer season is fairly short.
Incorporating the Length of the Seasons into the PopulationDynamics Model. To further investigate the hypothesis that season length is essential for determining the dynamics, we reparameterized the seasonal model to reveal the relative length of the winter as a parameter τ_{w}.†† This allows separation of the effects of season length per se from the effect of density dependence per time within the season. A seasonal model will then be given as [4a][4b] where seasonal density dependencies a_{wi} and a_{si} in Eqs. 1a and 1b now are given as a*_{wi}·τ_{w} and a*_{si}·(1 – τ_{w}), respectively. The annual direct (α_{1}) and delayed (α_{2}) density dependencies expressed in terms of seasonal components are and The linear effects of winter length on α_{1} and α_{2} are then and Given that density dependencies are stronger (i.e., more negative) during the winter than the summer, and that the winter is the longer season, ∂α_{1}/∂τ_{w} and ∂α_{2}/∂τ_{w} will quite often be negative.
Using the pooled estimates reported in Table 1, the annual direct and delayed density dependencies typically become smaller (i.e., the strength increases) as winter becomes longer (i.e., dα_{1}/dτ_{w} and dα_{2}/dτ_{w} both become negative) (Fig. 1C). On this basis, we predict that noncyclic dynamics are more likely in areas with short winters and that periodic cycles are more likely in areas with long winters, just as seen empirically (ref. 9; also see Fig. 1B). Although the densitydependent structure is best explained by whether the population belongs to Group 1, 2, or 5 (Fig. 1B), the length of winter (τ_{w}) seems to be a good predictor of the densitydependent structure. This dependency of the relative winter length can, however, only be substantiated when using the (better) pooled estimates reported in Table 1, not when using the estimates for the individual populations.
Conclusion
Altogether we have shown that winter length has a clear effect on the population dynamics of Hokkaido voles. Essentially, increasing winter length may be seen as a bifurcation parameter shifting the dynamics from an intrinsically stable regime with irregular fluctuations (generated by densityindependent mechanisms) to largeramplitude, periodic cycles influenced by densitydependent mechanisms. As the required delayed densitydependent structure is widespread in voles (8), such a seasonbased bifurcation process may apply quite generally. The underlying mechanisms generating this densitydependent structure are the subject of contention, and may differ from one system to another; however, the relative length of the seasons will be a generic bifurcation parameter irrespective of the underlying mechanism. Hence, we suggest that seasonality represent a key for unlocking the mystery of population cycles.
Acknowledgments
We thank KungSik Chan for many very helpful discussions during the formulation of the model and the linking of the seasonal and annual processes. We also thank Torbjørn Ergon, Tero Klemola and Ole Chr. Lingjærde. Comments by Tim Coulson and two additional anonymous reviewers were greatly appreciated. We are indebted to the Forestry Agency of Japanese Government for providing the data analyzed in this paper. We acknowledge the support we, over the years, have received from Japan Society for the Promotion of Science, the Forestry and Forest Products Research Institute of Japan, the Norwegian Research Council (NFR/NT), the Center for Advanced Studies of the Norwegian Academy of Science and Letters, and the University of Oslo. This study was also supported in part by GrantinAid 14340240 from the Ministry of Education, Science, and Culture of Japan (to T.S.).
Footnotes

↵‡ To whom correspondence should be addressed. Email: n.c.stenseth{at}bio.uio.no.

This paper was submitted directly (Track II) to the PNAS office.

↵†† Our estimates for the relative length of the winter (τ_{w}) provide the best prediction [(WI – ŴI)^{2}] for Groups 1 and 2. To observe that Group 5 (the more mountainous group) is more difficult to estimate seems reasonable because topographic effects in mountainous regions (north/south sloping) are likely to be more pronounced.
 Received September 2, 2002.
 Copyright © 2003, The National Academy of Sciences
References
 ↵
Nicholson, A. J. (1933) J. Anim. Ecol., Suppl. 2, 132–178.

Andrewartha, H. G. & Birch, L. C. (1954) The Distribution and Abundance of Animals (Univ. of Chicago Press, Chicago).
 ↵
Rothery, P. (1998) in Insect Populations in Theory and in Practice, eds. Dempster, J. P. & McLean, I. F. G. (Kluwer Academic, Dordrecht, The Netherlands), pp. 97–133.
 ↵
Turchin, P. (1995) in Population Dynamics, eds. Cappuccino, N. & Price, P. (Academic, New York), pp. 19–40.
 ↵
 ↵
Bjørnstad, O. N., Champley, O. N., Stenseth, N. C. & Saitoh, T. (1996) Philos. Trans. R. Soc. London B 351, 867–875.
 ↵
 ↵
Fahrmeir, L. & Tutz, G. (1994) Multivariate Statistical Modelling Based on Generalized Linear Models (Springer, New York).
 ↵
 ↵
Stenseth, N. C., Kittilsen, M. O., Hjermann, D., Viljugrein, H. & Saitoh, T. (2002) Proc. R. Soc. London Ser. B 269, 1853–1863.
 ↵
Tatewaki, M. (1958) J. Fac. Agric. Hokkaido Univ. 50, 371–472.
 ↵
 ↵
 ↵
 ↵
 ↵
Ota, K., ed. (1984) Study on Wild Murid Rodents in Hokkaido (Hokkaido Univ. Press, Sapporo, Japan) (Japanese).
 ↵
Henttonen, H., Hansson, L. & Saitoh, T. (1992) Ann. Zool. Fenn. 29, 1–6.

Abe, H. (1976) J. Mammalog. Soc. Japan 7, 17–30.

Fujimaki, Y. (1972) Bull. Hokkaido For. Exp. Stn. 10, 59–67 (Japanese with English summary).

Fujimaki, Y. (1975) Bull. Hokkaido For. Exp. Stn. 13, 38–45 (Japanese with English summary).
 ↵
Stenseth, N. C., Bjørnstad, O. N. & Saitoh, T. (1996) Proc. R. Soc. London Ser. B 263, 1117–1126.
 ↵
 ↵
 ↵
Royama, T. (1992) Analytical Population Dynamics (Chapman & Hall, London).
 ↵
 ↵
Spiegelhalter, D. J., Thomas, A., Best, N. & Gilks, W. R. (1999) winbugs version 1.2 User Manual (MRC Biostatistics Unit, Inst. of Public Health, Cambridge, U.K.).
 ↵
 ↵
 ↵
Hansen, T. F., Stenseth, N. C., Henttonen, H. & Tast, J. (1999) Proc. Natl. Acad. Sci. USA 96, 986–991.pmid:9927680
 ↵

Korpimäki, E., Norrdahl, K., Klemola, T., Pettersen, T. & Stenseth, N. C. (2002) Proc. R. Soc. London Ser. B 269, 991–997.
 ↵
 ↵
 ↵
Saitoh, T., Stenseth, N. C., Viljugrein, H. & Kittilsen, M. O. (2003) Popul. Ecol., in press.
 ↵
Abe, H., ed. (1994) A Pictorial Guide to the Mammals of Japan (Tokai Univ. Press, Tokyo) (Japanese).

Turchin, P. & Hanski, I. (1997) Am. Nat. 149, 267–276.
 ↵
Turchin, P. (2003) Complex Population Dynamics: A Theoretical/Empirical Synthesis (Princeton Univ. Press, Princeton).
 ↵
The Mammalogical Society of Japan (1997) Red List of Japanese Mammals (Buninchi Sogo Press, Tokyo) (Japanese).
 ↵
Toyooka, K., Sato, A. & Ishizuka, M. (1981) North. For. 33, 3–6 (Japanese).