Coupling in goshawk and grouse population dynamics in Finland

Different prey species can vary in their significance to a particular predator. In the simplest case, the total available density or biomass of a guild of several prey species might be most relevant to the predator, but behavioural and ecological traits of different prey species can alter the picture. We studied the population dynamics of a predator–prey setting in Finland by fitting first-order log-linear vector autoregressive models to long-term count data from active breeding sites of the northern goshawk (Accipiter gentilis; 1986–2009), and to three of its main prey species (1983–2010): hazel grouse (Bonasa bonasia), black grouse (Tetrao tetrix) and capercaillie (T. urogallus), which belong to the same forest grouse guild and show synchronous fluctuations. Our focus was on modelling the relative significance of prey species and estimating the tightness of predator–prey coupling in order to explain the observed population dynamics, simultaneously accounting for effects of density dependence, winter severity and spatial correlation. We established nine competing candidate models, where different combinations of grouse species affect goshawk dynamics with lags of 1–3 years. Effects of goshawk on grouse were investigated using one model for each grouse species. The most parsimonious model for goshawk indicated separate density effects of hazel grouse and black grouse, and different effects with lags of 1 and 3 years. Capercaillie showed no effects on goshawk populations, while the effect of goshawk on grouse was clearly negative only in capercaillie. Winter severity had significant adverse effects on goshawk and hazel grouse populations. In combination, large-scale goshawk–grouse population dynamics are coupled, but there are no clear mutual effects for any of the individual guild members. In a broader context, our study suggests that pooling data on closely related, synchronously fluctuating prey species can result in the loss of relevant information, rather than increased model parsimony.     

How reliable are harvesting data for analyses of spatio‐temporal population dynamics?

Harvest data are commonly used as proxy for count data, especially in studies of long‐term temporal and spatial patterns of population fluctuations. However, usually the concurrence of the conclusions based on different types of data is impossible to verify due to the lack of count data. Here, we use annual (1964–2004) harvest and population census data for capercaillie, black grouse and hazel grouse from 14 game management districts covering Finland, and demonstrate some mismatch in the information that these data sets provide. Overall, linear regressions of annual harvest against population count give a reasonable fit, but the slopes are less than 1 in every species. Harvest bags have been proportionally larger in north and eastern Finland than in southwestern Finland, with marked species‐specific differences. Considering population variation, the CV% in the census data (30–50%) is consistently smaller than it is in the harvest data (60–70%). Most importantly, conclusions on the spatio‐temporal patterns of the population dynamics are different if based on harvest rather than count data. In capercaillie, synchrony decreases faster with distance according to the harvest data, while in black grouse and hazel grouse the census data show the steeper decline. In addition, the autocorrelation coefficients in the census time series are higher in capercaillie and black grouse than in harvest data, but in hazel grouse the opposite is true. Finally, the parameter estimates for a second order autoregressive model using different data sets differ, and these differences are species‐specific. Despite the fact that annual harvest is a positive and linear function of annual grouse population density, the pattern of population dynamics derived from the bag data is different from that shown by the census data. This result urges caution in using wildlife bag data as reliable indices of population dynamics. deceased August 2008.     

Responses by the capercaillie Tetrao urogallus, and the willow grouse Lagopus lagopus, to the green matter available in early spring

Food matter eaten in the first snowless spots early in the spring by capercaillie Tetrao urogallus and willow grouse Lagopus lagopus was studied in Finnish Forest Lapland. When the snow disappears, both species change to feeding mainly on the plants available in the snowless spots, although male capercaillie does not exploit this nutritious diet to the same extent as female capercaillie and willow grouse. Selection in favour of nitrogen and phosphorus-rich food items, i.e. Betula pubescens , male catkins Eriophorum vaginatum , flower buds, and Equisetum spp., stems and tips, seems to be characteristic especially of the female capercaillie, whose winter diet is poor in these elements. Female capercaillie also feeds on more plant species or different parts of plants at one time than willow grouse, while no differences were recorded between male and female willow grouse in the composition of their spring food. The crowberry Empetrum hermaphroditum , almost regularly produces a bountiful berry crop in northern Finland, and since it overwinters well, it represents an easy source of energy and water available almost every spring. The quality and/or the quantity of the green food matter available early in the spring may fluctuate considerably and be of importance for short-term fluctuations in the populations of these tetraonid species.     

Sex and age influence intestinal parasite burden in three boreal grouse species

Parasite infections are often more common in male vertebrates than in females. Sexual selection leading to dimorphism can cause sexual differences in immune defence, behaviour or body size. Possible proximate explanations for male bias in parasitism are the immunosuppressive effect of male sex hormones and the large body size of males which increases the likelihood of being parasitized. To evaluate these hypotheses, we studied the prevalence and abundance of intestinal helminth parasites in three boreal grouse species, the capercaillie Tetrao urogallus , the black grouse Tetrao tetrix , and the hazel grouse Bonasa bonasia . The first two are sexually dimorphic polygynous species while the latter species is sexually monomorphic and monogamous. We found a male-bias in the prevalence and abundance of the nematode Ascaridia compar . The bias was most pronounced in the polygynous black grouse and capercaillie. In the monogamous hazel grouse, there was a slight male-bias in occurrence of ascarids, but no bias in abundance. In juvenile grouse, the male-bias was larger than in adult grouse. No sexual bias was found in regard to the cestodes ( Skrjabinia cesticillus , Paroniella urogalli and Hymenolepis spp.). However, age was a factor in cestode prevalence: juvenile grouse were more commonly infected than adults. Differences in growth rates and body size are potential factors that may lead to male-biased parasitism in these grouse species, and their impact requires further studies.     

Temporal dynamics of grouse populations at the southern edge of their distribution

Studies of grouse conducted at northern latitudes have shown that tetraonids frequently exhibit cyclic fluctuations in abundance but little is known about the dynamics of grouse species at the southerly edge of their range. Hunting statistics from four species of grouse based on 30 yr of data collected from 210 hunting areas were examined from the Dolomitic Alps in the province of Trentino. These data were summed to represent 18 time series from discrete mountain groups. Analyses identified cycles of ca 5 yr in the minority of rock ptarmigan Lagopus mutus and hazel grouse Bonasa bonasia populations. These cycles only showed significant negative autocorrelation at half the cycle period and were classified as phase-forgetting quasi-cycles. Cycles were not found in time series of black grouse Tetrao terix or capereaillie Tetrao urogallus. Correcting time series for hunting effort or hunting restrictions tended to increase the proportion of populations that exhibited cycles but no difference in the strength of second order density dependence, A linear first order density-dependent autoregressive model described the dynamics of most of the populations with the exception of a proportion of rock ptarmigan and black grouse populations where a non linear first order model provided the best tit. We compare the findings with studies conducted in Finland and suggest possible reasons for the reduced tendency to cycle in the populations of southern Europe.     

Diet of the golden eagle Aquila chrysaetos during the breeding season in Sweden

During a five-year period, 1975–1979, a total of 2881 prey individuals of 65 prey species were collected at 162 golden eagle nests from northern Sweden and from the island of Gotland. In northern Sweden birds are taken in higher numbers than mammals but calculated as weight the two categories are of equal importance, The main prey during the breeding season are capercaillie, black grouse, willow grouse, ptarmigan, mountain hare and reindeer fawns which together form 91% of the total food biomass. The capercaillie and the black grouse are taken more in the southern part of the coniferous region than in the northern. In contrast, in northern areas, reindeer fawns are more preyed upon than in the South, Ptarmigan and willow grouse are the most commonly captured prey species in mountain areas. The total number of reindeer fawns taken (dead and/or alive) by the Swedish golden eagle population during one summer is estimated at 600 individuals. On Gotland the golden eagles take mammals more often than in its northern distribution area. Rabbit and hedgehog arc the most important species.     

Searching for mechanisms of synchrony in spatially structured gamebird populations

1.  Time series data on five species of gamebird from the Dolomitic Alps were used to examine the relative importance of dispersal and common stochastic events in causing synchrony between spatially structured populations.
2.  Cross-correlation analysis of detrended time series was used to describe the spatial pattern of fluctuations in abundance, while standardized time series were used to describe both fluctuations and the trend in abundance. There were large variations in synchrony both within and between species and only weak negative relationships with distance.
3.  Species in neighbouring habitats were more likely to be in synchrony than species separated by several habitats. Species with similar density-dependent structure were more likely to be in synchrony.
4.  In order to estimate the relative importance of dispersal and environmental stochasticity, we modelled the spatial dynamics of each species using two different approaches. First, we used estimating functions and bootstrapping of time series data to calculate the relative importance of dispersal and stochastic effects for each species. Second, we estimated the intensity of environmental stochasticity from climatic records during the breeding season and then modelled the dispersal rate and dispersal distance for each species. The two models exhibited similar results for rock ptarmigan, black grouse, hazel grouse and rock partridge, while contrasting patterns were observed for capercaillie.
5.  The results suggest that environmental stochasticity plays the dominant role in synchronizing the fluctuations of these galliform species, although there will also be some dispersal between populations.     

Cycles in voles and small game in relation to variations in plant production indices in Northern Sweden

Population dynamics for voles (Cricetidae), Tengmalm's owl (Aegolius funereus (L.)), red fox (Vulpes vulpes (L.)) willow grouse (Lagopus lagopus (L.)), black grouse (Lyrurus tetrix (L.)), capercaillie (Tetrao urogallus L.), hazel hen (Tetrastes bonasia (L.)), mountain hare (Lepus timidus L.) and tularemia (Francisella tularensis (McCoy & Chapin)) and game bird recruitment were studied by index methods in northern Sweden. In addition contemporary temperature records and spruce (Picea abies (L.) Karst.) and pine (Pinus silvestris L.) cone crops (as indices for plant production) and the occurrence of forest damage, caused by voles eating bark, were studied.During 1970–80 two synchronous 4-year cycles were observed for voles, predators (Tengmalm's owl and red fox) and their alternative prey species (grouse and mountain hare). In grouse the change of numbers was correlated with that of recruitment. Autumn vole numbers peaked about a year before the other species and extensive forest damage occurred at winter peak densities of voles. These population fluctuations are consistent with a predator-prey model for their regulation. In short the model suggests that vole-food plant interactions trigger the cycle of voles, that voles generate the cycle of predators and that these in turn synchronize alternative prey populations to the others at vole declines.For voles, grouse and red fox the amplitude was higher in the first cycle compared to the second one whilst the opposite was true for the mountain hare. Although temperature and cone crops showed large interannual variations they still implied that herbivore food conditions were better during the former cycle. Hence, the reduction of the amplitude of the vole cycle may be explained by inter-cyclic differences in plant food conditions, implying food shortage (as indicated by bark-eating) at different population levels. The similar decrease of grouse and red fox populations may also be explained by deteriorated food conditions and/or for the fox by an outbreak of sarcoptic mange (Sarcoptes scabiae var. vulpes). The increased amplitude of the mountain hare cycle was part of a long-term rise in numbers after a tularemia epidemic in 1967. This is interpreted as a recovery, probably towards the generally higher pre-epidemic population level.     

Goshawk predation during winter, spring and summer in a boreal forest area of central Sweden

Predation by goshawks was studied in a central Swedish boreal forest area. Data were collected in winter (January–February) 1977-81 by tracking radio-tagged goshawks, and in the breeding season (April–July) by collecting prey remains at the nest. In the breeding season birds dominated the prey, amounting to 86% of prey number and 91% of prey biomass. Wood pigeon Columba palumbus , black grouse Tetrao tetrix , hooded crow Corvus corone cornix and jay Garrulus glandarius accounted for more than 50% of the prey animals, whereas capercaillie Tetrao urogallus and black grouse accounted for more than 50% of prey biomass. There was no functional response to black grouse density fluctuations. Every year goshawks killed significantly more females than males of both capercaillie and black grouse, due to high vulnerability of the grouse hens while laying and incubating. It was estimated that during spring and early summer goshawk predation removed 25% of the female, and 14% of the male black grouse population. In winter squirrel was the dominating prey, both in terms of number (79%) and weight (56%). The proportion of squirrel in the diet was equally high both in winters of low and high squirrel density. The high proportion of squirrel in the winterdiet, as compared to the breeding season, is believed to be due to squirrels having to accept an increased predation risk in winter, in order to feed efficiently enough.     

Looking into the past – the reaction of three grouse species to climate change over the last million years using whole genome sequences

Tracking past population fluctuations can give insight into current levels of genetic variation present within species. Analysing population dynamics over larger timescales can be aligned to known climatic changes to determine the response of species to varying environments. Here, we applied the Pairwise Sequentially Markovian Coalescent (psmc ) model to infer past population dynamics of three widespread grouse species; black grouse, willow grouse and rock ptarmigan. This allowed the tracking of the effective population size (N_e) of all three species beyond 1 Mya, revealing that (i) early Pleistocene cooling (~2.5 Mya) caused an increase in the willow grouse and rock ptarmigan populations, (ii) the mid‐Brunhes event (~430 kya) and following climatic oscillations decreased the N_e of willow grouse and rock ptarmigan, but increased the N_e of black grouse and (iii) all three species reacted differently to the last glacial maximum (LGM) – black grouse increased prior to it, rock ptarmigan experienced a severe bottleneck and willow grouse was maintained at large population size. We postulate that the varying psmc signal throughout the LGM depicts only the local history of the species. Nevertheless, the large population fluctuations in willow grouse and rock ptarmigan indicate that both species are opportunistic breeders while black grouse tracks the climatic changes more slowly and is maintained at lower N_e. Our results highlight the usefulness of the psmc approach in investigating species’ reaction to climate change in the deep past, but also that caution should be taken in drawing general conclusions about the recent past.     

Short-term storage and cryopreservation of black grouse Tetrao tetrix and capercaillie T. urogallus semen

Assisted reproductive technologies can be an important part of programs directed for maintenance and protection of genetic variability. The objective of this study was to develop methods for liquid storage and cryopreservation of semen capercaillie and black grouse semen. Our results provide for the first time evidence for successful short-term storage and cryopreservation of capercaillie and black grouse semen using criteria of sperm motility characteristics as quality indices. Sperm motility could be protected up to 48 h liquid storage; however, 24-h storage should be preferable. Cryopreservation secured 40–60% post-thaw motility (as compared with control), both for freshly collected and 24-h-stored semen. In conclusion, a procedure for short-term and cryopreservation of capercaillie and black grouse semen are now available. These assisted reproductive technologies can be implemented into captive breeding programs for these species.     

The importance of alder to hazel grouse in Fennoscandian boreal forest: evidence from four levels of scale

Buds and staminate catkins of alder (Alnus spp) form an important winter food for hazel grouse Bonasa bonasia in the Fennoscandian boreal forest Alder was found to be highly preferred over other deciduous trees, particularly alders ≥ 10 m tall and ≤15 m from spruce forest Winter territories were probably feeding territories, as size was correlated negatively with alder density and almost significantly correlated negatively with competitor density All winter territories were found to contain ample winter food resources for hazel grouse However, the distribution of territories was associated significantly with the distribution of alders at two levels of scale, the territory level and the landscape level Moreover, relationships between the abundance of alders and hazel grouse were found at two additional levels of scale the local patch level and the biogeographic region level This agreement of the results from four levels of scale strongly suggested that the abundance and distribution of alder was a major factor limiting hazel grouse winter territories within dense Norway spruce Picea abies forests in the boreal zone of Fennoscandia Alder was relatively uncommon and exhibited a clumped dispersion pattern at the local and landscape scales, being associated with wet and rich soils The close relationship to alder implies that hazel grouse winter habitats, even in natural forests, also should be distributed patchily Hazel grouse may select the catkins and buds of alder because it is a very nutntous food source, and small species, such as the hazel grouse, require more nutritious food than larger species     

Flowering of allergenically important plant species in relation to the North Atlantic Oscillation system and thermal time in the Czech Republic

This paper analyses long-term (1960–2015) onset of flowering in 16 native terrestrial plants (11 of them produce important allergens) recorded in different parts of the Czech Republic (southern, central and northern part) in relation to the North Atlantic Oscillation (NAO) index of the preceding winter and thermal data—growing degree-days (GDD) and soil temperature. Flowering occurred significantly earlier following positive winter NAO phases (causing spring to be warmer than normal in Central Europe) in nearly all early-flowering (February, March, April) species; high Pearson correlation values were recorded in, e.g. wood anemone, common snowdrop, goat willow, common hazel and common alder. There was found a difference between the southern and northern part of the country, e.g. in silver birch and pedunculate oak. Out of the later-flowering (May–July) plant species, black elder and meadow foxtail also significantly correlated with the winter NAO index, lime tree correlated less markedly. The best results of a threshold for calculation of GDD to onset of beginning of flowering were found in lime tree—it was 5 °C at all three stations. Results of other taxa were more variable (e.g. 4–7 °C in goat willow; 6–10 °C in silver birch). Pearson correlation coefficients between NAO index and GDD were negative in lime tree at all thresholds (5, 6, 7, 8, 9, 10 °C), while goat willow and silver birch were not so uniform (both positive and negative values). Correlation coefficients between phenophase onset and soil temperature (10 cm depth) had the highest values in silver birch, European larch and wood anemone. Stations situated at higher elevation showed negative correlation coefficient with soil temperature in common snowdrop, pedunculate oak, meadow foxtail and lime tree; other values were positive.     

Nest loss in capercaillie and black grouse in relation to the small rodent cycle in southeast Norway

Summary The relationship between nest loss in boreal forest grouse and the fluctuations in small rodents was studied at Varaldskogen in southeast Norway during 1979–1986, covering two complete rodent cycles. Nest loss in capercaillie (N=174) and black grouse (N=81) was calculated according to Mayfield (1975) based on nests from radio-equipped hens (N=77) and nests found by other methods (N=178). Small rodent density was measured by snap trapping during spring and autumn. Losses varied as predicted by the classical alternative prey hypothesis (Hagen 1952 and Lack 1954, as elaborated by Angelstam et al. 1984): high losses during rodent crash years (85.5% capercaillie, 51% black grouse), and smaller losses during peak years (54.5% capercaillie, 32.5% black grouse). Losses were inversely related to autumn abundance of rodents in capercaillie (P<0.05), but the correlation was not significant for black grouse (0.10<P<0.20). In capercaillie, the only species with an adequate sample for analysis, no relationship was detected between spring density of rodents and nest loss. Losses during the prepeak years were nearly as high as during crash years, a result inconsistent with the model. We conclude that the numerical response of predators to their cyclic main prey (i.e. small rodents) probably play a main role during the low phase and prepeak year, whereas the dietary shift is most important during the peak and crash year of the cycle.     

Guardian or threat: does golden eagle predation risk have cascading effects on forest grouse?

Previous studies on intraguild predation have mainly focused on within-class assemblages, even though avian top predators may also influence mammalian mesopredator prey. By using nation-wide long-term data from Finland, northern Europe, we examined the impacts of golden eagles (Aquila chrysaetos) together with red foxes (Vulpes vulpes) and pine martens (Martes martes) on forest-dwelling herbivores, black grouse (Tetrao tetrix) and hazel grouse (Tetrastes bonasia). We hypothesized that eagles may alleviate the overall predation pressure on grouse by imposing intraguild predation risk on mesopredators. The predation impact of eagle was modelled using eagle density estimates and distance to eagle nest. Wildlife triangle counts were used as predation impact proxies of mammalian mesopredators and as measures of response in grouse. Our results show that eagle density correlated negatively with black grouse abundance indices while being positively associated with the proportion of juveniles in both grouse species, irrespective of the abundance of mesopredators. Yet, foxes and martens alone had a negative effect on the abundance indices and the proportion of young in the two grouse species. This suggests that the possible cascading effects of eagles are not mediated by decreased mesopredator numbers, but instead by fear effects. Alternatively, they may be mediated by other species than fox or marten studied here. In conclusion, we found support for the hypothesis that eagles provide protection for juvenile black and hazel grouse, whereas they are a threat for adult grouse. This important information helps us to better understand the role of avian top predators in terrestrial ecosystems.     

Past and potential future population dynamics of three grouse species using ecological and whole genome coalescent modeling

Studying demographic history of species provides insight into how the past has shaped the current levels of overall biodiversity and genetic composition of species, but also how these species may react to future perturbations. Here we investigated the demographic history of the willow grouse (Lagopus lagopus), rock ptarmigan (Lagopus muta), and black grouse (Tetrao tetrix) through the Late Pleistocene using two complementary methods and whole genome data. Species distribution modeling (SDM) allowed us to estimate the total range size during the Last Interglacial (LIG) and Last Glacial Maximum (LGM) as well as to indicate potential population subdivisions. Pairwise Sequentially Markovian Coalescent (PSMC) allowed us to assess fluctuations in effective population size across the same period. Additionally, we used SDM to forecast the effect of future climate change on the three species over the next 50 years. We found that SDM predicts the largest range size for the cold‐adapted willow grouse and rock ptarmigan during the LGM. PSMC captured intraspecific population dynamics within the last glacial period, such that the willow grouse and rock ptarmigan showed multiple bottlenecks signifying recolonization events following the termination of the LGM. We also see signals of population subdivision during the last glacial period in the black grouse, but more data are needed to strengthen this hypothesis. All three species are likely to experience range contractions under future warming, with the strongest effect on willow grouse and rock ptarmigan due to their limited potential for northward expansion. Overall, by combining these two modeling approaches, we have provided a multifaceted examination of the biogeography of these species and how they have responded to climate change in the past. These results help us understand how cold‐adapted species may respond to future climate changes.     

Cartographic representation of the distribution of black grouse (<Emphasis Type="Italic">Lyrurus tetrix</Emphasis> L.) and hazel grouse (<Emphasis Type="Italic">Tetrastes bonasia</Emphasis> L.) in the West Siberian plain

Distribution of black grouse and hazel grouse was evaluated on the basis of the results of bird counts from July 16 to August 31, 1959–2005 over all natural geographical zones and subzones of West Siberia (with fractioning down to the type of landscape tract). Data averaging over groups of vegetation map units (dynamic series) and other specific averaging procedures revealed that black grouse is common in northern-, middle-, and southern-taiga bogs; in subtaiga pine forests; in grassy bogs combined with halophytic meadows within the boundaries of subtaiga forests and the northern forest-steppe; and in agricultural lands in the northern forest-steppe. Hazel grouse is numerous in southern-taiga and dark coniferous middle-taiga forests and their derivatives. The number of hazel grouse individuals in these habitats is larger than in other lands. In general, over the subzones, both black grouse and hazel grouse are most numerous in the middle and southern taiga and in subtaiga forests. To the north and to the south of these subzones, the abundance of these species decreases. They were not encountered in tundras and some open forest-steppe and steppe habitats.     

Use of hunters in wildlife surveys: does hunter and forest grouse habitat selection coincide?

Wildlife monitoring performed by volunteer personnel may suffer from bias with regard to their habitat use. Such errors can lead to erroneous population estimates, evidently influencing both management programmes and research that are based on the monitoring. We used a dataset on hunters’ habitat use in forest while searching for black grouse (Tetrao tetrix) and capercaillie (Tetrao urogallus) to test whether hunters’ utilisation of the habitat was independent of that of grouse or if it corresponded to the grouse habitat preference. Twenty volunteer hunters with dogs registered their tracks and all observations of capercaillie and black grouse in Østfold County, Norway, during August 2003 and 2004. We performed an ecological niche factor analysis (ENFA) and a K-select analysis with respect to hunters’ selection for habitat, described with ecogeographical variables related to forest stand characteristics, as well as the habitat preference of the observed grouse, conditional on the habitat utilisation of the hunters. Individual ENFA on the hunter's tracks revealed large variation in the habitat preferences of the hunters. The K-select indicated few overall patterns in the habitat characteristics of grouse observations, conditional on the hunters selected habitat. Accordingly, the results indicate that hunters’ observation of grouse prior to the hunt may give indicators of changes in grouse density unbiased by habitat preference due to the large between-hunter variation in habitat preference, given that a sufficient number of hunters is used. This suggests that such monitoring programmes can provide information about fluctuations in grouse population sizes valuable for both the management and research of forest grouse species.     

Seasonal changes in territorial behaviour of hazel grouse ( Bonasa bonasia) in a temperate forest of South Korea

We studied the seasonal changes in territorial behaviour of hazel grouse (Bonasa bonasia) from December 1999 to July 2001 in a temperate forest of South Korea. The study was based on the use of radio telemetry. Territory size and degree of overlap increased after the breeding period, and index of advertisement and border disputes decreased. In spring, the territory was defended by both males and females. During autumn and winter, this behaviour was relaxed to such a degree that hazel grouse cannot be considered territorial, but as forming a home range.     

Inference of hazel grouse population structure using multilocus data: a landscape genetic approach

In conservation and management of species it is important to make inferences about gene flow, dispersal and population structure. In this study, we used 613 georeferenced tissue samples from hazel grouse (Bonasa bonasia) where each individual was genotyped at 12 microsatellite loci to make inference on population genetic structure, gene flow and dispersal in northern Sweden. Observed levels of genetic diversity suggest that Swedish hazel grouse do not suffer loss of genetic diversity compared with other grouse species. We found significant F(IS) (deviation from Hardy-Weinberg expectations) over the entire sample using jack-knifed estimators over loci, which is most likely explained by a Wahlund effect. With the use of spatial autocorrelation methods, we detected significant isolation by distance among individuals. Neighbourhood size was estimated in the order of 62-158 individuals corresponding to a dispersal distance of 950-1500 m. Using a spatial statistical model for landscape genetics to infer the number of populations and the spatial location of genetic discontinuities between these populations we found indications that Swedish hazel grouse are divided into a northern and a southern population. We could not find a sharp border between these two populations and none of the observed borders appeared to coincide with any potential geographical barriers.These results imply that gene flow appears somewhat unrestricted in the boreal taiga forests of northern Sweden and that the two populations of hazel grouse in Sweden may be explained by the post-glacial reinvasion history of the Scandinavian Peninsula.