On the basis of questionnaires to moose hunters the colonisation pattern of reintroduced beaver in the province of Värmland was mapped. The beaver spread rapidly within water drainage systems resulting initially in scattered colonisation and very low population densities in relation to range. Spread was slowed by watershed divides. Illustrating the shape of the expansion wave front by plotting the relative increase of the population and range against time, shows a very flat front characterised by long dispersal distances.
Reintroducing exterminated species to areas within their former ranges is a conservational tool of increasing importance (Stanley Price 1989, Stuart 1991). Reintroduction programmes have not been numerous and most of them are fairly recent. To be able to predict the future development of reintroductions, it is important to study earlier cases and events with similar processes, such as deliberate or accidental introductions of non-native species and spontaneous invasions.
Organisms spreading into uninhabited areas is a phenomenon that has received much attention in spite of the relative scarcity of good data sets. The general picture of a spreading population is that of a rolling wave where the distribution increases linearly with time. There have been many attempts to model the process of spread of organisms (e.g. Skellam 1951, Okubo 1980, Hengeveld 1989, Kareiva 1990, van den Bosch et al. 1990, Holmes et al. 1994). Most models are based on the assumption that movements of individuals are indistinguishable from diffusion on a given geographical scale. Thus, the adequacy of a particular model can only be evaluated relevant to the purposes at hand, and to particular scales of investigation (Andow et al. 1990, Gardner et al. 1991, Hengeveld 1989, Johnson et al. 1992, Kareiva 1990).
Several models take special factors affecting the process into consideration, such as topographic irregularities (Andow et al. 1990, Gardner et al. 1991), finite velocity and non-random movement (Holmes 1993), Allee effects (Lewis & Kareiva 1993), and density-dependent movements (Namba 1980). The main goal has almost always been to predict the wave velocity. This has been criticised, however, because the trend of the radial equivalent on time does not reveal much about the processes by which populations spread. It may allow rejection of extreme hypotheses, but it does not differentiate between a number of very different processes (Caughley 1977). Several authors have emphasised the importance of the shape of the wave front, which is much determined by the frequency of long distance dispersers (e.g. Mollison 1986, Hengeveld 1989, Holmes et al. 1994, Gardner et al. 1991). Some have also warned about the danger in neglecting the tip of the tail of the dispersal distribution, since even a small frequency of long distance dispersers may be of importance (Andow et al. 1990, van den Bosch et al. 1992, Porter & Dooley 1993, Goldwasser et al. 1994). Thus, there is a need for more information on shapes of expansion wave fronts and the underlying factors determining these shapes.
The last well-documented observations of European beaver Castor fiber from the original Swedish population date back to the 1870s (Ekman 1910). By the time authorities banned beaver hunting in 1873, it was already too late. The cause of extermination was apparently overex-ploitation because of the high prices of pelts and castoreum. By the beginning of this century, only a small remnant of the Scandinavian beaver population remained in southern Norway. In 1922, a pair of beavers were purchased from Norway and released in the province of Jämtland in central Sweden. About 80 Norwegian beavers were released at 19 different sites in Sweden between 1922 and 1939. Reproduction was observed at 11 of these sites (Fries 1940, Fig.1). Since then beavers have spread, and are still spreading, over large parts of the country, aided partially by translocations. In 1992, the population size was estimated at about 100,000 individuals (Hartman 1995).
The aim of this study was to describe the pattern of spread of a reintroduced beaver population in the province of Värmland in southwestern Sweden, and to discuss the results in a general context with special reference to the shape of expansion fronts and its importance for introductions.
The province of Värmland comprises an area of about 17,600 km2 (see Fig. 3). Boreal forest covers about 75% of the area. Altitudes vary from 40 to 690 m, but only 23% of the area is higher than 200 m above sea level. The area can be divided into nine major water drainage systems which mainly flow in a north/southerly direction (see Fig. 3). The last beavers of the original population in Värmland were probably killed during the 1830s (Ekman 1910). In 1925, two beaver pairs were imported from southern Norway and reintroduced to the eastern parts of the province. This first introduction was followed by another (2 pairs) at the same site in 1927. Other introductions at a site 40 km southeast of the first site took place in 1928 (1 pair) and 1930 (1 pair), but these introductions probably were unsuccessful (Fries 1940). In 1960, 10 animals were translocated from the eastern parts of the province to a site approximately 45 km west of the expansion front (see Fig. 3).
The administration of moose hunting in Sweden is based on a system of license areas that consist of the land of one or several landowners or part of a larger property where hunting rights are rented by a hunting team. These hunting teams usually rent the same land year after year and often include the landowners. In 1987, I performed a beaver survey based on questionnaires sent to hunters in 426 moose license areas. One of the questions asked in the questionnaire was when beavers first settled in the hunting area. The accuracy of the answer depends mainly on the memory of the observer. As public interest in the species has been high since the first introduction and as the occurrence of beaver in new areas is often discussed and has frequently been reported in newspapers, this is an event likely to be remembered by local landowners and hunters. 1 received answers to the questionnaires from 322 areas covering 9,127 km2 (x = 28.3 km2). I studied the pattern of spread by marking colonisation for each year on computorised maps run in sequence. To determine the shape of an expansion wave front, when determined more by life history traits than by topography, one has to look within water drainage systems. I did this by plotting the relative increase in the population and in the range against time (referred to as p.r-plots). In the plots I used a population size index (the number of areas colonised in a certain year as a proportion of the number colonised by 1987), and a relative measure of range (the maximum distance between the two colonised areas furthest apart in a certain year as a proportion of the same measure in 1987). As I wanted to document the slope of the front, range expansion was simplified and viewed as if the spread was primarily one-dimensional. I did this by drawing a line through the middle of a water drainage system in its longitudinal direction, projecting the positions of colonised areas perpendicular to the line, and measuring distances between projection points (Fig. 2). If all dispersers settle in the first empty area they encounter (steep wave front), the plotted lines of population size and range against time will track each other closely. If, however, the wave front is flat, the range index will initially increase faster than the population index, and the longer the tail of the front, the larger a gap between the lines should be expected.
In 1945, 20 years after the first introduction, beavers were found 122 km northwest and 38 km southeast of the release point (Fig. 3a). This was also the first year, when beavers were observed west of the Klarälven drainage system (area A).
By 1950 the range had expanded even further north-west, and was now 306 km. Expansions in other directions were more limited, and the population was distributed in a very patchy way (Fig. 3b).
In 1960, the range covered most of the eastern drainage systems and the southern front had moved to 155 km south of the release point (Fig. 3c). At this time a considerable filling in of gaps in the distribution had occurred. The western front, however, had hardly moved at all and remained approximately 100 km west of the release point. This means that the yearly spread had been ca 3 km since 1925. In 1960, 10 beavers were translocated to area B (see star in Fig. 3c). The most striking change in the following decade was the resulting rapid spread of beavers in area B, probably due to the above-mentioned translocation.
In 1967, a beaver colony was found in the northwestern corner of area B; thus the yearly spread in area B was 19.7 km. This means that beavers in area B dispersed much faster than beavers in area A, where maximum dispersal rates were 12–13 km/year.
In 1970, the distribution pattern in area B was very patchy, similar to that observed earlier in the eastern parts of the province. By 1970, colonised areas could be found in almost the entire province, except for the southern parts along the shores of Lake Vänern, and area C in the south-western corner (Fig. 3d).
During 1970–1980 distribution became more continuous (Fig. 3e,f). A sharp front along the watershed divide formed the border of area C. This last hindrance to dispersal was passed in the mid-1980s, and by 1987, when the survey was performed, the beaver population range covered the whole province. There were, however, great differences in population density reflecting differences in time passed since colonisation (Hartman 1994).
Drawing p.r.-plots of the two water drainage systems with the largest number of reporting areas (A and B) showed that the relative increase in range was considerably faster than the relative increase in number of colonised areas (Fig. 4), thus indicating a flat expansion front. The whole province has been surveyed on five occasions after the reintroduction (Faxén 1939, Lundberg et al. 1965, Weinberg 1971, Lavsund 1979, Hartman 1994). A comparison of the relative population increase in the whole province with the relative increase in number of occupied areas, showed that the p.r.-plots using number of occupied areas as an index for population size under-estimated the difference between population and range increase (Fig. 5).
In spite of the well-documented ability of beavers Castor fiber and C. canadensis to travel over land (Berghofer 1961, Wilsson 1971. Hodgdon 1978, Smith & Peterson 1988) the velocity of dispersal within and between drainage systems differs. Watershed divides were by no means absolute dispersal barriers, but they slowed the spread in certain directions and thereby had a considerable influence on the pattern of spread. The colonisation maps and p.r.-plots revealed a pattern with a flat wave front formed by long distance dispersers obviously passing large areas of suitable habitat before settling down. German and Polish reports of beavers moving long distances in the early colonisation phase after reintroduction also exist (Zuppke 1989, Zurowski & Kaspercyk 1988). The prevalence of long distance dispersers in expansion fronts has likewise been recorded in other species, for example, nine-banded armadillo Dasypus novemcinctus (Humphrey 1974), raccoon dog Nyctereutes procyonoides (Helle & Kauhala 1991) and shiny cowbird Molothrus bonariensis (Post et al. 1993).
It seems reasonable to assume that in most species, dispersal behaviour is adapted to finding vacant territories and mates within a relatively stable species range, rather than invading vast areas of suitable habitat, unoccupied by conspecifics even though so-called super-tramp species (Diamond 1975) exploiting temporary habitats, form an exception. This means that dispersal patterns may be different between more normal situations and an expansion front. If dispersers at the wave front search not only for an unoccupied site but also for a mate, they might keep moving further away from the core population. This might explain the long dispersal distances of beavers in the early stages of the colonisation of Värmland. Many areas that were passed in these early stages of colonisation later reached high population densities. This indicates that dispersing beavers passed them, not because they lacked suitable habitats, but for some other reason, e.g. because they lacked mates.
The risk of undercrowding effects, sometimes referred to as Allee effects (Allee 1938), has sometimes been mentioned in connection with introductions (Williamsson 1989, May 1991) and has been suggested as the reason why many attempts have failed (Hopper & Roush 1993). High dispersal rates from the reintroduction site were considered to be a major cause behind the problems with re-introducing beaver in Bavaria (Reicholf 1977). Even if an introduction is successful, population development might be limited by Allee effects in the expansion front (Hopper & Roush 1993). Caughley (1970) found no evidence of an Allee effect in the expansion front of Himalayan thar Hemitragus jemlahicus in New Zealand. His explanation was that, as dispersal in Himalayan thar is sex-biased, females dispersed into areas with a sufficient density of males. This situation is not likely to occur in a beaver population because dispersal distances are only weakly biased towards males (Berghofer 1961, Leege 1968).
Long-distance dispersal and flat wave fronts are by nomeans a general pattern in introductions. A p.r.-plot of the development of the reintroduced sea otter Enhydra lutris on the North American west coast (Lubina & Levin 1988) shows a radically different pattern with lines lying close to each other. In the case of the sea otter, the relative population increase was even faster than the relative increase in range, thus indicating a very steep wave front. Another illustrating example of how the pattern of spread of two different species may differ, is the invasions of American mink Mustela vison and coypu Myocastor coypus in England (Usher 1986). The coypu has a closed wave front leaving almost no suitable habitat behind it, whereas the spread of the mink is characterised by long-distance dispersers making the front very flat and broken up.
There is no way of knowing if the beaver population in Värmland initially suffered Allee effects. The rapid spatial dilution of the population observed within the drainage systems in the early colonisation phase, suggests that it did. The plot of the population increase for the whole province (Fig. 5) showed the expected long initial lag phase, which also supports the hypothesis, albeit weakly. If there was an initial risk that the reintroduction could have failed due to Allee effects, then the occurrence of watershed divides, which slowed down the spread and thereby increased local density, could have been a crucial factor.
I thank the hunters in Värmland for supplying me with the necessary information on beaver colonisation. I also want to thank Erik Lindström, Jon Swenson, Allan Carlsson and Lennart Hansson for valuable comments on the manuscript and Christina Hultgren for locating the watershed divides. The Swedish Environmental Protection Board funded the project.
- Allee, W.C. 1938: The social life of animals. - W.W. Norton and Co., New York, 293 pp. Google Scholar
- Andow, D.A., Kareiva, P.M., Levin, S.A. & Okubo, A. 1990: Spread of invading organisms. - Landscape Ecology 4: 177–188. Google Scholar
- Berghofer, C.B. 1961: Movement of beaver. - 40th. Annual Conference of Western Association of State Game and Fish Commissioners, pp. 177–181. Google Scholar
- Caughley, G. 1970: Liberation, dispersal and distribution of himalayan thar (Hemitragus jemlahicus) in New Zealand. - New Zealand Journal of Science 13: 200–239. Google Scholar
- Caughley, G. 1977: Analysis of vertebrate populations. - John Wiley & Sons, 234 pp. Google Scholar
- Diamond, J.M. 1975: Assembly of species communities. - In: Cody, M.L. & Diamond, J.D. (eds.); Ecology and evolution of communities. Harvard University Press, Cambridge, U.S.A, pp. 342–444. Google Scholar
- Ekman, S. 1910: Norrlands jakt och fiske. - Norrländskt handbibliotek, No. 4. Uppsala, Sweden, 473 pp. (In Swedish). Google Scholar
- Faxén, L. 1939: Data ur bäverns spridningshistoria. - Jämtlands läns jaktvårdsförenings årsbok 1939: 157–163. (In Swedish). Google Scholar
- Fries, C. 1940: Bäverland. - Nordisk Rotogravyr, Stockholm, Sweden, 246 pp. (In Swedish). Google Scholar
- Gardner, R.H., Turner, M.G., O'Neill, R.V. & Lavorel, S. 1991: Simulation of scale-dependent effects of landscape boundaries on species persistence and dispersal. - In: Holland, M.M., Risser P.G. & Naiman R.J. (eds.); Ecotones; The role of landscape boundaries in the management and restoration of changing environments. Chapham and Hall, New York, pp. 76–89. Google Scholar
- Goldwasser, L., Cook, J. & Silverman, E.D. 1994: The effects of variability on metapopulation dynamics and rates of invasion. - Ecology 75: 40–47. Google Scholar
- Hartman, G. 1994: Long-term population development of a reintroduced beaver (Castor fiber) population in Sweden. - Conservation Biology 8: 713–717. Google Scholar
- Hartman, G. 1995: Population development of European beaver in Sweden. - Third Nordic Beaver Symposium 1992, Ilomantsi, Finland, 54 pp. Google Scholar
- Helle, E. & Kauhala, K. 1991: Distribution history and present status of the raccoon dog in Finland. - Holarctic Ecology 14: 278– 286. Google Scholar
- Hengeveld, R. 1989: Dynamics of Biological Invasions. - Chapham and Hall, London, 160 pp. Google Scholar
- Hodgdon, H.E. 1978: Social dynamics and behavior within an unexploited beaver (Castor canadensis) population. - Ph.D.thesis, University of Massachusetts, 292 pp. Google Scholar
- Holmes, E.H. 1993: Are diffusion models too simple? A comparison with telegraph models of invasion. - American Naturalist 142: 779–795. Google Scholar
- Holmes, E.E., Lewis, M.A., Banks, J.E. & Veit, R.R. 1994: Partial differentiation equations in ecology: Spatial interactions and population dynamics. - Ecology 75: 7–29. Google Scholar
- Hopper, K.R. & Roush, R.T. 1993: Mate finding, dispersal, number, released, and the success of biological control introductions. - Ecological Entomology 18: 321–331. Google Scholar
- Humphrey, S.R. 1974: Zoogeography of the nine-banded armadillo (Dasypus novemcinctus) in the United States. - BioScience 24: 457–462. Google Scholar
- Johnson, A.R., Milne, T.B. & Wiens, J.A. 1992: Diffusion in fractal landscapes: Simulations and experimental studies of Tenebroid beetle movements. - Ecology 73: 1968–1983. Google Scholar
- Kareiva, P. 1990: Population dynamics in spatially complex environments: theory and data. - In: Hassell, M.P. & May, R.M. (eds.); Population, regulation and dynamics. Philosophical Transactions of the Royal Society of London B 330: 175–190. Google Scholar
- Lavsund, S. 1979: Bäverinventering i Värmlands län. - Department of Wildlife Ecology, Swedish University of Agricultural Sciences, Uppsala, Sweden, 28 pp. (In Swedish). Google Scholar
- Leege, T.A. 1968: Natural movements of beavers in southeastern Idaho. - Journal of Wildlife Management 32: 973–976. Google Scholar
- Lewis, M.A. & Kareiva, P. 1993: Allee dynamics and the spread of invading organisms. - Theoretical Population Biology 43: 24–37. Google Scholar
- Lubina, J.A. & Levin S.A. 1988: The spread of a reinvading species: range expansion in the Californian sea otter. - American Naturalist 131:526–543. Google Scholar
- Lundberg, S., Petterson, J.O. & Weinberg, U. 1965: Bäver. - Sveriges Natur 56: 134–142. (In Swedish). Google Scholar
- May, R. 1991: The role of ecological theory in planning re-introduction of endangered species. - In: Gipps, J.W.H. (ed.); Beyond captive breeding. Reintroducing endangered mammals to the wild. - Zoological Society of London Symposia 62: 145–163. Google Scholar
- Mollison, D. 1986: Modelling biological invasions: chance, explanation, prediction. - In: Kornberg, H. & Williamson, M.H. (eds.); Quantitative aspects of the Ecology of biological invasions. Philosophical Transactions of the Royal Society of London B 314: 675–693. Google Scholar
- Namba, T. 1980: Density-dependent dispersal and spatial distribution of a population. - Journal of Theoretical Biology 86: 351– 363. Google Scholar
- Okubo, A. 1980: Diffusion and Ecological problems: Mathematical models. - Springer Verlag, Berlin, 254 pp. Google Scholar
- Porter, J.H. & Dooley J.L. Jr . 1993: Animal dispersal patterns: A re-assessment of simple mathematical models. - Ecology 74: 2436– 2443. Google Scholar
- Post, W., Cruz, A. & McNair, D.B. 1993: The North American invasion pattern of the shiny cowbird. - Journal of Field Ornithology 64: 32–41. Google Scholar
- Reicholf, J.H. 1977: The reintroduction of the beaver (Castor fiber) in Bavaria: some preliminary results. - Proceedings from Seminar on reintroductions: techniques and ethics, Rome, WWF, pp. 49–53. Google Scholar
- Skellam, J.G. 1951: Random dispersal in theoretical populations. - Biometrika 38: 196–218. Google Scholar
- Smith, D.W. & Peterson, R.O. 1988: The effects of regulated lake levels on beaver in Voyageurs national park, Minnesota. - U.S. Department of the Interior, National Park Service, Research/Resources Management Report MWR-1. Midwest Regional Office, Omaha, Nebraska 68102, 84 pp. Google Scholar
- Stanley Price, M. 1989: Animal reintroductions: The Arabian oryx in Oman. - Cambridge University Press, 291 pp. Google Scholar
- Stuart, S. N. 1991: Reintroductions: to what extent are they needed? - In: Gipps, J.W.H. (ed.); Beyond captive breeding. Reintroducing endangered mammals to the wild. Zoological Society of London Symposia 62: 27–37. Google Scholar
- Usher, M.B. 1986: Invasability and wildlife conservation: invasive species on nature reserves. - In: Kornberg, H. & Williamson, M.H. (eds.); Philosofical Transactions of the Royal Society of London B 314: 695–710. Google Scholar
- van den Bosch, F., Metz J.A. & Diekmann, O. 1990: The velocity of spatial population expansion. - Journal of Mathematical Biology 28: 529–565. Google Scholar
- van den Bosch, F., Hengeveld, R. & Metz, J.A. 1992: Analysing the velocity of animal range expansion. - Journal of Biogeography 19: 135–150. Google Scholar
- Weinberg, U. 1971: Bäverinventering. - Swedish Environmental Protection Agency, Solna, Sweden, 69 pp. (In Swedish). Google Scholar
- Williamsson, M. 1989: Mathematical models of invasion. - In: Drake, J.A., Mooney, H.A., di Castri, F., Groves, R.H., Kruger, F.J., Rejmánek, M. & Williamson, M. (eds); Biological Invasions: a global perspective, SCOPE 37, John Wiley & Sons, pp. 329–350. Google Scholar
- Wilsson, L. 1971: Observations and experiments on the ethology of the European beaver (Castor fiber L.). - Viltrevy 8, 266 pp. Google Scholar
- Zuppke, U. 1989: Besiedlungstendenzen des Elbebibers, (Castor fiber albicus), im Kreis Wittenberg. - Hercynia, N.F. Leipzig. 26: 351–361. (In German). Google Scholar
- Zurowski, W. & Kaspercyk, B. 1988: Effects of reintroduction of European beaver in lowlands of the Vistula basin. - Acta Theriologica 33: 325–338. Google Scholar