Long-term patterns in European brown hare population dynamics in Denmark: effects of agriculture, predation and climate
© Schmidt et al; licensee BioMed Central Ltd. 2004
Received: 24 June 2004
Accepted: 12 October 2004
Published: 12 October 2004
In Denmark and many other European countries, harvest records suggest a marked decline in European brown hare numbers, a decline often attributed to the agricultural practice. In the present study, we analyse the association between agricultural land-use, predator abundance and winter severity on the number of European brown hares harvested in Denmark in the years 1955 through 2000.
Winter cereals had a significant negative association with European brown hare numbers. In contrast to this, root crop area was positively related to their numbers. Remaining crop categories were not significantly associated with the European brown hare numbers, though grass out of rotation tended to be positively related. The areas of root crop production and of grass out of rotation have been reduced by approximately 80% and 50%, respectively, while the area of winter cereals has increased markedly (>70%). However, European brown hare numbers were primarily negatively associated with the number of red fox. Finally, we also found a positive association between mild winters and European brown hare numbers.
The decline of Danish European brown hare populations can mainly be attributed to predation by red fox, but the development in agricultural land-use during the last 45 years have also affected the European brown hare numbers negatively. Additionally, though mild winters were beneficial to European brown hares, the increasing frequency of mild winters during the study period was insufficient to reverse the negative population trend.
In most countries in Western Europe, the landscape has undergone dramatic changes during the last century due to changes in the agricultural practices. In Denmark, both the uncultivated land and the semi-cultivated land, such as permanent grass areas, have been reduced dramatically, reflecting the general intensification of agriculture . Additionally, fields have become larger, which has resulted in widespread fragmentation of remaining habitats, and today the landscape appears as a mosaic of natural habitats surrounded by cultivated land [1, 2]. These changes in agriculture have affected a number of wildlife species living in this man-made landscape. For instance, the shift in agricultural practice has severely influenced the diversity and abundance of insects with concomitant consequences for the dynamics of a wide range of farmland birds [3, 4], including the grey partridge (Perdix perdix) . Among mammals, the European brown hare (Lepus europaeus) in particular has experienced a dramatic decline in many European countries [reviewed by ], including Denmark [7, 8].
Despite its currently declining numbers, the European brown hare is still common and one of the most important game species throughout the country . The dynamics of European brown hares seem resilient to even heavy hunting pressure , though local population dynamic data may be needed to ensure sustainable harvest . In Denmark hunting of European brown hares is generally assumed to be without regulating effect . The European brown hare is a typical grass steppe herbivore, and inhabits primarily open landscapes, including cultivated farmland , which is the predominant landform in Denmark . The species is rather sedentary, and has generally small home ranges [e.g. ]. This site fidelity makes European brown hares highly susceptible to changes in their surrounding habitats, and the general decline in the European brown hare populations in Europe is mainly being attributed to changes in agriculture practice and land-use [reviewed by [6, 12]].
European brown hares are important prey primarily for mammalian predators. In Northern Europe, the red fox (Vulpes vulpes) is the main predator on European brown hares, and foxes have been reported to influence the dynamics of several European brown hare populations substantially [13–17]. Also, infectious diseases such as the European brown hare syndrome virus, pseudotuberculosis, pasteurellosis and coccidiosis are present in many European brown hare populations [18–20]. Haerer et al. , however, concluded that diseases were not responsible for the decline of brown hare populations in Switzerland. Similarly, Frölich et al.  found that compared to red foxes, infectious diseases seemed to play a minor role in the dynamics of European brown hare populations in Germany.
An increasing number of papers have documented the importance of climate for a number of life history traits and abundance of many terrestrial vertebrates , and European brown hare populations are affected negatively by cold winters and cold springs [12, 22–24]. Factors regulating vertebrate populations may, however, exhibit large spatial variation, and even among proximate populations spatial variation and gradients in vertebrate population dynamics may exist [e.g. ].
In this study, we analyse and contrast the simultaneous associations between agricultural land-use, the number of red fox harvested, winter severity and the population dynamics of European brown hares across 11 Danish districts during 45 years.
We found statistical significant direct density-dependence (X t-1 ) in the European brown hare time series (G = 281.4, df = 2, P < 0.0001). Additionally, the effect of district was also significant (G = 769.0, df = 1, P < 0.0001).
Summarised results from the analysis of the impact of agricultural land-use, red fox and winter climate on European brown hare harvest records from 11 Danish counties, 1955–2000. Also given is the change in model deviance (Δ deviance) explained by the variable when fitted last. Type 3 sums of squares.
Winter cereals (t)
Spring cereals (t)
Grass in rotation (t)
Grass out of rotation (t)
Winter rape (t)
Spring rape (t)
Root crops (t)
The number of red foxes harvested in year t-1 had a marked negative effect on the number of European brown hares harvested the following year (Table 2), whereas red fox number in year t seemed unimportant for the European brown hare numbers.
Despite its high reproductive potential [e.g. ], the Danish European brown hare has, according to annual harvest records, declined dramatically since 1955. Studies of Danish European brown hare populations indicate that its reproductive success is low compared to that of con-specifics in other countries , and has, in turn, declined from the 1950s to the 1990s . Hansen  suggested that the low reproductive success observed in Danish European brown hares might be attributed to the agricultural practice. Using data covering almost half a century, our analyses suggest that the dramatic decline in the Danish European brown hares can be attributed mainly to the negative effect of red foxes, but also to the agricultural land-use. The area of winter cereal production has increased dramatically during the last century, and became the dominating crop in the early 1990s. We also found a significant positive association between root crops, a crop type that has declined dramatically in the second half of the last century, and European brown hare numbers.
European brown hares mainly forage on grasses and herbs , and cereals such as wheat are preferred during winter [29, 30], which seems to contradict the results of our analyses. However, as European brown hares choose to feed on specific crops according to plant phenology [12, 31], cereals, although important in winter, still occupy large areas when no longer of nutritional value to European brown hares, which may result in low availability of food especially during summer. Similarly, rape is avoided in the diet , but European brown hares may spend a substantial fraction of their time in rape fields during winter  prior to the development of glucosinulates [see ]. Apart from winter cereals and root crops, the crop categories did not affect the European brown hare numbers significantly. The lacking effect of grass and green fodder areas, especially those out of rotation, was unexpected as other studies have shown that hares prefer grass areas year-round [12, 30]. This, however, may be attributed to the fact that we were unable to separate grass areas into those e.g. with and without cattle, factors that might affect European brown hare use of grass areas .
In a recent study, Fox  showed that farmland birds seemed to benefit from the reduced application of pesticides and inorganic fertilisers seen in Denmark since the early 1980s. The continuing decline in European brown hare numbers in that period therefore suggests that hares are not directly affected by the use of such substances, but rather respond to the loss of suitable habitat and space. European brown hares mainly move along field margins [12, 33], and the increasing field size  combined with the general loss of suitable habitats possibly force hares to aggregate in the remaining patches of non-agricultural and non-urbanized land. Here, density-dependence may perpetuate the negative population development as European brown hares aggregate in the few, remaining pockets of suitable habitat. In line with this, Frylestam  reported an inverse relationship between fertility and density in European brown hares, which he attributed to shortage of food at least in some parts of the year, which also has been suggested in a number of other studies [12, 24, 27]. Hence, agricultural land-use, especially the increasing use of winter cereals and the marked reduction in root crops, but probably also grass areas out of rotation, seem to have contributed to the European brown hare decline in Denmark. However, the single-most important parameter for the number of European brown hares was the number of red foxes. Hence, our results are consistent with other studies reporting that red foxes may affect European brown hare populations negatively through predation [13–17]. This relationship is also particular evident from the positive development in the European brown hare population on the island Bornholm following the outbreak of sarcoptic mange there . To elaborate the fox-hare interaction further, we reran the analyses including all variables but the red fox variables, and added the delayed AR term (i.e. X t-2 ). In seven of the 11 districts the delayed density dependence was significant (P < 0.05), suggesting important inter-specific interactions [36, 37]. The significant association with red fox numbers in year t-1 (and not year t ) most likely reflects that compared to adults, juvenile European brown hares suffer more from predation [e.g. ]. Hence, high red fox numbers in year t result in low harvest of European brown hares in year t , which affects the reproductive potential of the populations, and, hence, the number of European brown hares the following year (year t+1 ). Also, the significant effect of district may point to differences in habitat quality, but also differences in the history of the sarcoptic mange, i.e. time since the eruption of the mange, among districts [see ].
Both European brown hare over-winter survival  and reproductive rate  increases with temperature, which may be attributed to improved forage availability following mild winters [see [38, 39]]. Our analyses revealed similar results as mild winters affected the European brown hares positively. There may, however, also be a negative effect of mild winters, namely through transmission of diseases and parasites, which may be enhanced under mild climatic conditions . Nevertheless, the overall effect of mild winters seemed positive. The upward trend in the NAO since the 1960s , however, was not sufficient to reverse the European brown hare decline, but may have decelerated it.
Our analyses have provided important insight into the structure of the European brown hare dynamics in Denmark, and documented important patterns within the mechanisms regulating European brown hares. Hence, we have documented that the decline of European brown hares in Denmark mainly can be attributed to predation by red fox, but also to changes in agricultural land-use. Additionally, though mild winters were beneficial to European brown hares, the increasing frequency of mild winters during the study period was insufficient to reverse the negative population trend.
European brown hare and red fox density indices
As indices of European brown hare and red fox density, we used the harvest records from 14 Danish counties from 1955–2000. In 1970, the geographical borders of some of the counties were changed, which in two cases lead to substantial changes in area . This, together with inconsistency in crop registration, forced us to lump together some counties, and we therefore present analyses of European brown hare harvest records from 11 districts (Fig. 1). European brown hare and red fox harvest records were expressed as the number of animals shot per hectare per year. Harvest records may seem rather crude indicators of density compared to e.g. spotlight surveys and line transects [see ]. Unfortunately, no alternative indices of European brown hare density in Denmark are available. Harvest records, however, may still be useful indicators of the long-term trends in European brown hare numbers .
Sarcoptic mange was first encountered in Denmark in the early 1980s, and is now present in seven of the 11 districts . In one district (Bornholm; Fig. 1), the disease has practically eliminated the entire red fox population, and red fox hunting has been prohibited here since 1991. Consequently, we restrict our analyses of the Bornholm district to 1955 through 1990. In all other districts, red fox hunting is open from 1 September to 31 January, and for European brown hares from 1 October to 31 December. Neither red fox nor European brown hare hunting is quota-based.
Description of model variables used in the analyses of European brown hare harvest records from Denmark, 1955–2000.
Winter cereals, including winter wheat, winter barley, and winter rye
Spring cereals, including spring wheat, spring barley, spring rye, and oat
Grass in rotation
Grass and green fodder in rotation, including grass, clover, maize, and Lucerne
Grass out of rotation
Grass areas permanently out of rotation
Sugar beets for sugar production and fodder, mangolds, turnips, and carrots for feed
The number of red fox harvested
The winter state of the North Atlantic Oscillation
No data on the fraction of winter and spring cereals, respectively, were available for 1955–61, and these data points were omitted in the analyses. Similarly, prior to 1982 no data on the fraction of winter barley and spring barley were available. However up until 1968, winter barley was only sown on a few thousand hectares in Denmark, when it was forbidden due to its function as reservoir for mildew attacking spring barley . The use of mildew resistant winter barley was permitted in 1982. Hence, from 1955–81 barley was regarded as spring barley. Finally, no data on the fraction of winter and spring rye are available from 1955–61 and 1979–98, but in Denmark spring rye is generally sown on very few hectares only , and we regarded rye as winter rye unless otherwise specified [see also ]. Oat is sown as spring crop only. No data on the areas covered with grass and green fodder in or out of rotation were available for 1994, and the data point was interpolated, i.e. 1994 equalled the average of 1993 and 1995.
As an index of winter severity, we used the winter state of the North Atlantic Oscillation [NAO; ]. The climatic phenomenon NAO is defined as the difference in sea-level pressure between Portugal and Iceland, affecting the temperature, precipitation and wind across the Northern Hemisphere [26, 43, 44]. The NAO, thus, integrates the effects of a number of abiotic factors, and therefore seems particularly useful when analysing the dynamics of larger, endothermic animals with relatively large buffer-capacity against climatic variability.
The relationship between the winter NAO (December t-1 through March t ) and temperature, precipitation and wind is particularly strong in Northern Europe . When NAO is in its high state, the winter climate in Northern Europe is warm, wet, and windy, while when in its low state winters are cold and dry . The NAO winter index has shown an upward trend since the 1960s, which accounts for a substantial fraction of the general warming in the Northern Hemisphere . The NAO index data are available from the Climate Analysis Section, National Center for Atmospheric Research, USA http://www.cgd.ucar.edu/~jhurrell/nao.html.
To remove heteroscedasticity in the system, data were loge-transformed, i.e. X (t) = ln(N (t) + 1), where N (t) is the number of European brown hares harvested in year (t) . To allow the direct comparison of regression coefficients, variables were standardized prior to the analyses (i.e., [X (t) -µ(1955–2000)]/s(1955–2000)). Regression coefficients therefore express the rate of change in standard deviation units of the independent variable per one standard deviation unit of the dependent variable . To obtain stationarity in the time series [see  for details], data were detrended following  by including year as covariate in all models.
For each district, we tested for the presence of multi-collinearity among parameters prior to the analyses by means of condition indices and variance proportions , but multi-collinearity was not observed in any districts (condition index < 12.47; [see  for details]).
The European brown hare time series were tested for non-linearity in X (t) vs. X (t-1) by applying the likelihood ratio test . In all districts, linearity in X (t) vs. X (t-1) was not rejected λ< 7.76, P > 0.05). We analysed the data using a first-order autoregressive (AR(1)) mixed model approach with district as fixed factor:
X (t) ~(µ,V),
where LAND_USE represents the seven crop categories (see Table 1). Because of the inclusion of red fox numbers, we did not include a delayed AR term (i.e. X t-2 ), as inclusion of both delayed density dependence and predator abundance can be seen as redundant . All analyses were performed in SAS 8.2, using the PROC MIXED procedure with restricted maximum-likelihood estimation of regression coefficients . Model reduction was conducted using log-likelihood ratio tests .
We extend our sincere thanks to Eric Post and two anonymous referees for valuable comments and improvements to the manuscript.
- Bach H, Christensen N, Kristensen P: 2002, Copenhagen, Denmark, National Environmental Research Institute, 409: 1-368.The state of the environment in Denmark,Technical report,Google Scholar
- Caspersen OH: Landbrug og landskab fra 1800 til 2000. Grænser i landskabet. Edited by: HelsT, FrandsenJN, FritzbøgerB and OlesenCR. 2001, Denmark, Odense University Press, 111-134.Google Scholar
- Fuller RJ, Gregory RD, Gibbons DW, Marchant JH, Wilson JD, Carter N: Population declines and range concentrations among lowland farmland birds in Britian. Conserv Biol. 1995, 9: 1425-1441. 10.1046/j.1523-1739.1995.09061425.x.View ArticleGoogle Scholar
- Chamberlain DE, Fuller RJ, Bunce RGH, Duckworth JC, Shrubb M: Changes in abundance of farmland birds in relation to the timing of agricultural intensification in England and Wales. J Appl Ecol. 2000, 37: 771-788. 10.1046/j.1365-2664.2000.00548.x.View ArticleGoogle Scholar
- Benton TG, Bryant DM, Cole L, Crick HQP: Linking agricultural practice to insect and bird populations: a historical study over three decades. J Appl Ecol. 2002, 39: 673-687. 10.1046/j.1365-2664.2002.00745.x.View ArticleGoogle Scholar
- Edwards PJ, Fletcher MR, Berny P: Review of the factors affecting the decline of the European brown hare, Lepus europaeus (Pallas, 1778) and the use of wildlife incident data to evaluate the significance of paraguat. Agr Ecosyst Environ. 2000, 79: 95-103. 10.1016/S0167-8809(99)00153-X.View ArticleGoogle Scholar
- Strandgaard H, Asferg T: The Danish bag record II. Danish Review of Game Biology. 1980, 11: 1-112.Google Scholar
- Madsen J, Asferg T, Clausager I, Noer H: 1996, Copenhagen, Denmark, National Environmental Research Institute, 6: 1-112.Status og jagttider for danske vildtarter,Technical report,Google Scholar
- Pepin D: Variation in survival of brown hare (Lepus europaeus) leverets from different farmland areas in the Paris Basin (France). J Appl Ecol. 1989, 26: 13-24.View ArticleGoogle Scholar
- Marboutin E, bray Y, Péroux R, Mauvy B, Lartiges A: Population dynamics in European hare: breeding parameters and sustainable harvest rates. J Appl Ecol. 2004, 40: 580-591. 10.1046/j.1365-2664.2003.00813.x.View ArticleGoogle Scholar
- Frylestam B: Utilisation of farmland habitats by European hares (Lepus europaeus, Pallas) in southern Sweden. Swedish Wildlife Research. 1980, 11: 271-284.Google Scholar
- Tapper SC, Barnes RFW: Influence of farming practice on the ecology of the Brown hare (Lepus europaeus). J Appl Ecol. 1986, 23: 39-52.View ArticleGoogle Scholar
- Pielowski Z: The role of foxes in the reduction of the European hare population. Ecology and management of European hare populations. Edited by: PielowskiZ and PucekZ. 1976, Warszawa, Polish Hunting Association, 135-148.Google Scholar
- Lindström E, Andrèn H, Angelstam P, Widèn P: Influence of predators on Hare populations in Sweden: a critical review. Mammal Rev. 1986, 16: 151-156.View ArticleGoogle Scholar
- Goszczynski J, Wasilewski M: Predation of foxes on a hare population in central Poland. Acta Theriol. 1992, 4: 329-338.View ArticleGoogle Scholar
- Reynolds JC, Tapper SC: Predation by foxes Vulpes vulpes on brown hares Lepus europaeus in central southern England, and its potential impact on annual population growth. Wildlife Biol. 1995, 1: 145-158.Google Scholar
- Reynolds JC, Tapper SC: The ecology of the red fox Vulpes vulpes in relation to small game in rural southern England. Wildlife Biol. 1995, 1: 105-117.Google Scholar
- Rattenborg EE: Diseases in the hare (Lepus europaeus) population in Denmark - epidemiological description and analyses of post-mortem data. 1994, The Royal Veterinary & Agricultural University, Copenhagen, Denmark, 1-114.Google Scholar
- Haerer G, Nicolet J, Bacciarini L, Gottstein B, Giacometti M: Causes of mortality, zoonoses and reproductive performance in European brown hare in Switzerland. Schweizer Archiv Für Tierheilkunde. 2001, 143: 193-201.PubMedGoogle Scholar
- Frolich K, Wisser J, Schmüser H, Fehlberg U, Neubauer H, Grunow R, Nikolaou K, Priemer J, Thiede S, Streich WJ, Speck S: Epizootiologic and ecologic investigations of European brown hares (Lepus europaeus) in selected populations from Schleswig-Holstein, Germany. J Wild Dis. 2003, 39: 751-761.View ArticleGoogle Scholar
- Mysterud A, Stenseth NC, Yoccoz NG, Ottersen G, Langvatn R: The response of terrestrial ecosystems to climate variability associated with the North Atlantic Oscillation. The North Atlantic Oscillation - climatic significance and environmental impact. Edited by: Hurrell James W, Kushnir Yochanan, Ottersen Geir and Visbeck Martin. 2003, Washington, DC, American Geophysical Union, 235-262.View ArticleGoogle Scholar
- Andersen J: Studies in Danish Hare populations. I. Population fluctuations. Danish Review of Game Biology. 1957, 3: 85-131.Google Scholar
- Frylestam B: Structure, size, and dynamics of three European Hare populations in Southern Sweden. Acta Theriol. 1979, 24: 449-464.View ArticleGoogle Scholar
- Marboutin E, Hansen K: Survival rates in a nonharvested Brown hare population. J Wildl Manage. 1998, 62: 772-779.View ArticleGoogle Scholar
- Stenseth NC, Bjørnstad ON, Saitoh T: Seasonal forcing on the dynamics of Clethrionomys rufocanos: Modelling geographic gradients in population dynamics. Res Popul Ecol. 1998, 40: 85-95.View ArticleGoogle Scholar
- Hurrell JW: Decadal trends in the North Atlantic Oscillation: Regional temperatures and precipitation. Science. 1995, 269: 676-679.View ArticlePubMedGoogle Scholar
- Hansen K: Reproduction in European hare in a Danish farmland. Acta Theriol. 1992, 37: 27-40.View ArticleGoogle Scholar
- Hansen K, Hartmann P: 1994, Copenhagen, Denmark, Ministry of Environment and Energy, 125: 1-20.Unge harer i jagtudbyttet - en indikator for bestandens tilstand,Technical report,Google Scholar
- Frylestam B: Agricultural land use effects on the winter diet of Brown Hares (Lepus europaeus Pallas) in southern Sweden. Mammal Rev. 1986, 16: 157-161.View ArticleGoogle Scholar
- Hansen K, Sørensen PL, Daugaard O: Bestandsstørrelse og afgrødevalg hos hare (Lepus europaeus). 1989, Vildtbiologisk Station, Kalø, Landbrugsministeriets Vildtforvaltning, 1-12.Google Scholar
- Chapuis JL: Comparison of the diets of two sympatric lagomorphs, Lepus europaeus Pallas and Oryctolagus cuniculus (L.) in an agrosystem of the Ile-de-France. Zeitschrift für Säugetierkunde. 1990, 55: 176-185.Google Scholar
- Fox AD: Has Danish agriculture maintained farmland bird populations?. J Appl Ecol. 2004, 41: 427-439. 10.1111/j.0021-8901.2004.00917.x.View ArticleGoogle Scholar
- Lewandowski K, Nowakowski JJ: Spatial distribution of Brown hare Lepus europaeus populations in habitats of various types of agriculture. Acta Theriol. 1993, 38: 435-442.View ArticleGoogle Scholar
- Brandt J, Holmes E, Larsen D, Madsen MM: Småbiotoper i det danske agerland 1991. Landscape Ecological Papers. 1996, 3: 1-153.Google Scholar
- Forchhammer MC, Asferg T: Invading parasites cause a structural shift in red fox dynamics. Proc R Soc Lond B. 2000, 267: 779-786. 10.1098/rspb.2000.1071.View ArticleGoogle Scholar
- Stenseth NC, Falck W, Chan KS, Bjørnstad ON, O' Donoghue M, Tong H, Boonstra R, Boutin S, Krebs CJ, Yoccoz NG: From patterns to processes: Phase and density dependencies in the Canadian lynx cycle. Proc Natl Acad Sci. 1998, 95: 15430-15435. 10.1073/pnas.95.26.15430.PubMed CentralView ArticlePubMedGoogle Scholar
- Post E, Forchhammer MC: Pervasive influence of large-scale climate in the dynamics of a terrestrial vertebrate community. BMC Ecol. 2001, 1:5:Google Scholar
- Forchhammer MC, Clutton-Brock TH, Lindstrom J, Albon SD: Climate and population density induce long-term cohort variation in a northern ungulate. J Anim Ecol. 2001, 70: 721-729. 10.1046/j.0021-8790.2001.00532.x.View ArticleGoogle Scholar
- Mysterud A, Stenseth NC, Yoccoz NG, Langvatn R, Steinheim G: Nonlinear effects of large-scale climatic variability on wild and domestic herbivores. Nature. 2001, 410: 1096-1099. 10.1038/35074099.View ArticlePubMedGoogle Scholar
- Hurrell JW, Kushnir Y, Visbeck M: The North Atlantic Oscillation. Science. 2001, 291: 603-605. 10.1126/science.1058761.View ArticlePubMedGoogle Scholar
- Langbein J, Hutchings MR, Harris S, Stoate C, Tapper SC, Wray S: Techniques for assessing the abundance of Brown Hares Lepus europaeus. Mammal Rev. 1999, 29: 93-116. 10.1046/j.1365-2907.1999.00040.x.View ArticleGoogle Scholar
- Jensen KM, Reenberg A: Landbrugsatlas Danmark. 1986, København, Det Kongelige Geografiske Selskab og C.A. Reitzels Forlag, 1-118.Google Scholar
- Lamb PJ, Peppler RA: North Atlantic Oscillation: Concept and an application. Bull Amer Meteorol Soc. 1987, 68: 1218-1225. 10.1175/1520-0477(1987)068<1218:NAOCAA>2.0.CO;2.View ArticleGoogle Scholar
- Hurrell JW, Van Loon H: Decadal variations in climate associated with the North Atlantic Oscillation. Clim Change. 1997, 36: 301-326. 10.1023/A:1005314315270.View ArticleGoogle Scholar
- Ottersen G, Planque B, Belgrano A, Post E, Reid PC, Stenseth NC: Ecological effects of the North Atlantic Oscillation. Oecologia. 2001, 128: 1-14. 10.1007/s004420100655.View ArticleGoogle Scholar
- Sokal RR, Rohlf FJ: Biometrics. 1995, New York, W.H. Freeman and Company, 1-887. ThirdGoogle Scholar
- Priestley MB: Spectral analysis and time series. 1984, London, Academic Press, 1-890. ThirdGoogle Scholar
- Forchhammer MC, Stenseth NC, Post E, Langvatn R: Population dynamics of Norwegian red deer: density-dependence and climatic variation. Proc R Soc Lond B. 1998, 265: 341-350. 10.1098/rspb.1998.0301.View ArticleGoogle Scholar
- Inc. SAS Institute: The SAS system for Windows, version 8e. 2000, Cary, NCGoogle Scholar
- Tong H: Non-linear time series. 1993, Oxford, Clarendon Press, 1-564.Google Scholar
- Littel RC, Milliken GA, Stroup WW, Wolfinger RD: SAS Systems for Mixed Models. 1996, Cary, NC, USA, SAS Institute Inc.Google Scholar
This article is published under license to BioMed Central Ltd. This is an open-access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.