Introduction

Dying and standing dead trees (snags) have been recognized as a crucial resource for primary cavity excavators, particularly woodpeckers, providing foraging, and roost and nest sites (Short 1979; Winkler and Christie 2002). The relationship between dead-wood supply and the abundance of some woodpecker species has been well documented through studies demonstrating low woodpecker abundance when snags occur in very low densities, or are experimentally removed or created (Zarnowitz and Manuwal 1985; McPeek et al. 1987; Hutto and Gallo 2006; Kilgo and Vukovich 2014). However, other studies have found no such relationship across a range of measured snag densities, snag sizes, and snag conditions (e.g. Sedgewick and Knopf 1986; Gunn and Hagan 2000; Spiering and Knight 2005; Smith 2007).

The presence of standing dead wood may be especially beneficial for species with weaker excavating abilities, i.e. ‘snag specialists’ (Hågvar et al. 1990; Smith 2007; Lõhmus et al. 2010). However, past and present human activities, e.g. deforestation, impoverishment of tree species diversity, short-rotation schemes, active fire suppression, and removal of dead wood, have led to a decline in the numbers of dead and decaying trees in forest ecosystems (Angelstam and Mikusiński 1994; Mikusiński and Angelstam 1997; Mikusiński 2006). Some of the specialist, dead-wood dependent woodpecker species have declined dramatically in their geographic ranges (e.g. Bütler et al. 2004; Czeszczewik and Walankiewicz 2006; Mikusiński 2006; Roberge et al. 2008). Since the availability of snags suitable for cavity excavation in managed forests is low as a direct result of forestry practice, snags are often considered a limiting resource for woodpeckers (Newton 1994).

The great spotted woodpecker Dendrocopos major and the middle spotted woodpecker Leiopicus medius coexist in old, European deciduous forests (Winkler and Christie 2002; Kosiński 2006). Both species are primary excavators. The middle spotted woodpecker feeds predominantly on bark- and foliage-living arthropods by gleaning and probing their prey. The great spotted woodpecker exploits food resources at deeper levels of wood by excavating and probing on trees. During the nestling period, both species feed heavily on defoliating caterpillars. The beginning of the breeding season in both species is highly synchronized and does not vary between the species (Michalek and Miettinen 2003; Pasinelli 2003; Kosiński and Ksit 2006).

The mechanisms that allow for the coexistence of ecologically similar species include habitat partitioning, use of resources during different times, and/or differential use of food resources (Gause 1934; Schoener 1974). Earlier studies show that separation of the nesting niches of great spotted and middle spotted woodpeckers in riverine forests is expressed by the preference for different tree species, tree condition, and nest placement in relation to the forest edge (Kosiński and Winiecki 2004). Middle spotted woodpeckers show an affinity to nest near forest edges and favor nesting in snags, but great spotted woodpeckers more frequently nest in closed forests, and prefer nesting in live oaks Quercus sp. and alders Alnus glutinosa. However, other studies suggest that the type of forest stand, proportion of different tree species, age of tree stand, and share of dead trees may modify intra- and interspecific nest-site characteristics and niche partitioning of both woodpecker species (Wesołowski and Tomiałojć 1986; Kosiński et al. 2006; Pasinelli 2007; Domokos and Cristea 2014). To our knowledge, only one long-term study has examined the population response of three woodpecker species to increased snag density in Europe (Smith 2007). Furthermore, the ecological consequences of mass mortality of trees, such as those resulting from insect and disease outbreaks, for woodpecker nest-site selection in the forests of continental Europe are unknown.

In this study, we investigated the abundance and nest-site use of the great spotted and middle spotted woodpeckers in riverine forests in Central Europe as a response to a pulse in snag availability. This increased snag density, mainly oaks Quercus robur and ashes Fraxinus excelsior, was driven primarily by drought-mediated tree mortality (Schwartz 2008; Kamiński et al. 2011; Miler 2016). For oaks, mortality of stressed trees was affected by lethal attacks from oak tree borers, i.e. Agrilus biguttatus, in 2005 and 2006 (Wargo 1996; Thomas et al. 2002; Schwartz 2008). In an area of the oldest stands (59 ha), the density of dead oaks with a diameter at breast height (Dbh) ≥14 cm increased by 38%, from 5.4 snags 1 ha−1 in 2006 to 7.5 snags 1 ha−1 in 2007, and remained stable in the years 2009–2012 (Jarocin Forest District, unpubl. data). The causal agent of ash dieback is fungal disease caused by the ascomycete Hymenoscyphus fraxineus (Pautasso et al. 2013; Kowalski et al. 2016). Large-scale ash dieback, which first appeared in Poland in the 1990s (Przybył 2002), spread rapidly. The density of dead ashes (Dbh ≥14 cm) increased by 45%, from 1.6 snags 1 ha−1 in 2006 to 2.3 snags 1 ha−1 in 2007. However, in 2009 the density of dead ashes reached 9.8 snags 1 ha−1, and then decreased to 6.1 snags 1 ha−1 in 2012. Tree dieback was observed in variable intensity across the study area, probably as an effect of varying tree density and species composition, as well as the position of stands above the groundwater table.

Great spotted and middle spotted woodpeckers differ in their adaptations for drilling and drumming, i.e. the middle spotted woodpecker’s neck muscles are less suitable for excavation and drumming. These differences have been well documented not only through knowledge of species skeletomuscular systems (Jenni 1981), but also through studies demonstrating that both species differ in terms of nest-site characteristics; the latter species clearly prefers dead trees or excavates its cavities only in soft or decayed parts of living trees (Kosiński and Winiecki 2004; Kosiński et al. 2006; Pasinelli 2007). Although the role of dead or decaying wood as a resource for cavity excavation is well-documented, its impact on the abundance of the middle spotted woodpecker is controversial. Hertel (2003) suggests that lowland beech forests with a lot of standing dead wood can be colonized by middle spotted woodpeckers, even when the density of large oaks is low. Weggler et al. (2013) assume that increased amounts of standing dead wood are one of the possible explanations for the striking population growth of the middle spotted woodpecker in Zurich canton. However, other studies suggest that snags can be used for nesting if available, but their presence is not a necessary condition for great spotted and middle spotted woodpeckers’ occurrence and abundance (Lovaty 2002; Kosiński and Kempa 2007; Smith 2007; Stachura-Skierczyńska and Kosiński 2016; but see Spühler et al. 2014). As the middle spotted woodpecker seems to be largely confined to dead wood, we hypothesized that a pulse in snag availability should positively affect their abundance. Moreover, we expected that dead trees should be used with higher frequency rates in the outbreak than in pre-outbreak period.

The great spotted and middle spotted woodpeckers are the two most common woodpecker species in old-growth deciduous forests in Central Europe (Kosiński and Kempa 2007; Domokos and Cristea 2014; Wesołowski et al. 2015). As major creators of cavities, both species provide a critical resource for secondary cavity users (Mazgajski 2003; Smith 2005; Wesołowski 2007; Lučan et al. 2009). Moreover, the middle spotted woodpecker is listed in Annex I of the European Union Birds Directive, and is the subject of special conservation measures concerning their habitat to ensure their survival and reproduction in their area of distribution. Thus, identifying the role of habitat elements that are potentially critical requirements for woodpeckers species, e.g. dead wood, is an important step in developing promising conservation strategies, and this is essential information for policy planning bodies and decision-makers (Zawadzka et al. 2013).

In this article, we examine whether: (1) the abundance of the great spotted and middle spotted woodpeckers is associated with dead tree availability; (2) increases in dead wood availability affect intraspecific variation in nest-site use; and (3) changes in dead wood availability modify niche partitioning between two ecologically related species.

Materials and methods

Study area

We established a study site in the riverine forest of the Warta river valley, western Poland, near Czeszewo (17°31′E, 52°09′N), 50 km south east of Poznań. This woodland is a remnant of flood-plain forests, which in the past covered narrow alluvial strips along major European rivers. In some parts of Europe, this forest has completely disappeared. Currently, the total area of riverine forests in Poland is estimated at 27 km2 (Matuszkiewicz, 2001). Field work was carried out on a 222.6 ha plot, which encompasses 186 ha of forest. The vegetation consists of Quercus-Fraxinus-Ulmus (Fraxino-Ulmetum) woodland in the flooded parts and Quercus-Carpinus (Stellario-carpinetum) forest on the higher ground. The rest of the study area (39 ha, 17%) is formed by oxbow lakes and extensively used meadows. The study forest differs in terms of age, species composition and spatial structures because of previous human impact. The oldest, near-natural stands are distinguished from younger and human-altered stands with regards to the height of trees (up to 33 m in ashes, and 32 m in oaks), the multi-story profile of stands (4–5 layers), and the high diversity of tree species (up to 7 species). In 2000, about 40% of tree stands on the study plot (74 ha) were ca. 155–165 years old, 24% (43 ha)—81–120 years old, 28% (51 ha)—41–80 years old and 4% (10 ha)—≤40 years old; ca. 4% (8 ha) was covered by coniferous stands. By 2013, the proportion of the oldest forest stands (>120 years old) had increased to 62%, and the maximum age of stands reached 184 years. Since 1959, the oldest near-natural stands have been protected within two reserves (“Czeszewo” and “Lutynia”), and have been almost unmanaged. In 2004, the whole study area was declared a nature reserve: “Czeszewski Las”. Our study area is a part of both the NATURA 2000 Special Protection Area Middle Warta River Valley PLB300002, with the middle spotted woodpecker as one of the target species, and the Ostoja Nadwarciańska Special Area of Conservation PLH300009. Further characteristics of the forest stands and their management can be found elsewhere (Kosiński et al. 2006; Kamiński et al. 2011).

Woodpecker occurrence data

We counted middle spotted woodpeckers from 2000 to 2015 and great spotted woodpeckers in 2000–2009 and 2013–2015. However, to examine the potential effect of dead wood availability on nest-site characteristics we used data collected in 2000–2005 for middle spotted and in 2001–2005 for great spotted woodpeckers, and in 2013–2015 for both species. To estimate the number and distribution of middle spotted woodpecker territories, the study plot was surveyed in each of the pre-breeding seasons (from the second half of March till the end of April), when aggressive interactions among individuals are most common (Pasinelli et al. 2001). We performed four to five complete censuses during 2000–2005, and two to three censuses in other years. Censuses were conducted with a minimum interval of 1 week. All counts were based on the same call-playback method, which has been recommended for middle spotted woodpecker censuses as a way of overcoming difficulties in detection and mapping (Kosiński and Winiecki 2003; Kosiński et al. 2004; Stachura-Skierczyńska and Kosiński 2016). The survey was conducted by listening for and recording vocal responses of pairs/individuals at selected points. The number of playback locations amounted to 58 during a single census. To provoke responses from territorial birds, we used taped calls (rattle- and advertising-calls of males). To reduce the probability of some individuals being attracted away from their territories by the playback technique, we applied the following parameters: a minimum distance of 150–200 m between points, a maximum time of 40 s for the stimulation in one bout, and a minimum interval of 1–2 min for listening. After the first bird’s response, the taped calls were stopped to determine the initial location of the bird. We recorded and mapped the woodpecker’s initial position, behavior, type of call, direction of movements, and sex, if possible. We paid special attention to register simultaneously active birds to separate neighboring pairs/territories. To locate the positions of birds precisely, we walked along forest roads, division lines, and shore lines of oxbow lakes. Therefore, points of stimulation were distributed unevenly. The observations began in the morning at 7–8 AM and continued until late afternoon. Our experience suggests that time of day does not affect the detection of middle spotted woodpeckers in the pre-breeding season. Analysis of the data enabled an estimation of the number of birds holding territory on the study plot. We use the term “territory”, because middle spotted woodpeckers defend almost exclusive territories during March and April, with both sexes participating to similar degrees in territorial defense (Pasinelli et al. 2001). On the basis of our long-term experience, we assumed that two registrations in the pre-breeding period was sufficient to accept a territory (Lõhmus et al. 2016; Stachura-Skierczyńska and Kosiński 2016), and to exclude transient/floater birds (Pasinelli et al. 2001; Robles and Ciudad 2012). Moreover, two visits allow an experienced observer to find almost 100% of territories in the pre-breeding period (Z. Kosiński unpubl. data). All middle spotted woodpecker territories were checked to find a nest-hole, except for the years 2010–2012. A search for nest-holes began in the pre-breeding period (2000–2005 and 2008–2010); however, in other years we focused our field procedure mainly on the detection of calling nestlings. We conducted extensive nest searches by systematically walking each territory in the second part of May and at the beginning of June. In total, in both periods we found nest-holes in 60 and 77% territories, respectively. These values are most likely close to the number of breeding attempts that survive to fledgling period (Kosiński and Ksit 2006). However, the lower value in the first period may be affected by the lower level of experience of observers in nest-searching. The distribution of cavities reflected the territory distribution (Kosiński and Winiecki 2003).

No playback censuses were conducted for great spotted woodpeckers, because of its inappropriateness for this species owing to its high density in riverine forests and difficulties in separation of neighboring pairs (Tomiałojć 1980). In 2001–2005, the nest-holes of great spotted woodpeckers were found at the beginning of the breeding season either by the sound of excavation or signs of wood-boring, as well as by detection of adults and calling nestlings during the latter period. In 2013–2015, we detected nest-holes by listening for begging calls of the young. Since in this period nest-holes were sought on a smaller forest area of 154 ha, data from the pre-outbreak period were limited to the same area. We aborted attempts to check part of the study area due to the high number of fallen logs and leaning trees, which rendered the field work time-consuming and risky. Moreover, this area was rarely used as a nesting habitat by the middle spotted woodpecker. Based on the high survival rates of great spotted woodpecker nests (e.g. Mazgajski 2002; Kosiński and Ksit 2006), and the high efficiency of nest searching methods, we assumed that probably all the nesting attempts were recorded in each year. Almost every year we found a few pairs with newly fledged offspring, with poor flying ability outside the nest-holes. To assess the total number of breeding pairs, these families were added to the number of active cavities. Adults with fledglings that flew well were not included in the total number of pairs. Such a treatment reduces the risk of an erroneous population estimate.

Description of nest-sites and changes in dead tree numbers

In this analysis, we used nine parameters for nest-sites: tree species, tree condition (live vs. dead), part of tree (trunk vs. limb/branch), condition of tree fragment (live vs. dead), diameter at breast height, nest height above the ground, presence or absence of wounds/cracks at nest- site, presence or absence of fruiting bodies of polypores within 1 m of the cavity entrance, and location of nesting tree in relation to distance from the forest edge (forest edge ≤25 m vs. closed forest >25 m from an edge).

To describe the changes in the availability of standing dead trees in pre-outbreak and outbreak periods, we used data obtained from randomly selected circular plots with a radius of 10 m (0.033 ha) established in 46 territories of great spotted woodpeckers in 2001–2002 and 91 territories of middle spotted woodpeckers in 2000–2004, as well as in 137 and 72 territories of both species in 2013–2015, respectively. The maximum distance between nesting trees and selected points for each woodpecker species was defined as the radius of a circle with an area equal to the male core area in the pre-breeding season, i.e. 7.3 ha for the middle spotted woodpecker (Pasinelli et al. 2001), and 3.0 ha for the great spotted woodpecker (Bachmann and Pasinelli 2002). The detailed procedure for randomization was described thoroughly by Kosiński and Winiecki (2004). In each plot, we measured all trees with a diameter at breast height (DBH) sufficiently large for hole excavation, i.e. ≥14 cm for middle spotted woodpeckers and ≥19 cm for great spotted woodpeckers (Kosiński and Ksit 2007). In total, in the ‘pre-outbreak’ and ‘outbreak’ periods we measured 395 and 1066 available trees in the great spotted woodpecker territories, as well as 951 and 670 trees in the middle spotted woodpecker territories, respectively.

Statistical analysis

We used generalized linear models (GLMs), with binomial error distribution, and logit link function to test for differences in nest-sites of the great spotted and middle spotted woodpeckers before the outbreak and during the pulse in dead wood supply, with period (‘pre-outbreak’ and ‘outbreak’) as a response variable, and characteristics of nest-sites as the explanatory variables. Interspecific differences in the nest-sites of the two woodpecker species were compared using species as the categorical dependent variable (Pasinelli 2007; Vierling et al. 2009). To test the effect of pre-outbreak and outbreak conditions on nest-site use, we included all the original variables of interest and their interactions with the study period.

We used generalized linear mixed effects models (GLMMs) to test the effect of tree condition (‘alive’ and ‘dead’) as a fixed effect on the utilization of trees (‘available’ and ‘used’) as a response variable. Separate models were run for the great spotted and middle spotted woodpeckers. We treated the study period (‘pre-outbreak’ and ‘outbreak’) as a random effect. This approach was used to control for the overall increase in available dead trees when comparing selection between pre-outbreak and outbreak periods. To test the significance of the random effect, we performed the likelihood ratio test comparing the log likelihood functions of the full model and the model without the random effect.

Results of multiple regression analyses may be affected by strong correlations among the explanatory variables (Dormann et al. 2013). Therefore, the correlation structure of ordinal and continuous variables was investigated (separately for each woodpecker species) before GLM analysis. Spearman rank correlations indicated one pair of moderately/strongly correlated explanatory variables (p < 0.05), these being tree condition with the condition of the tree fragment (great spotted woodpecker—|r S | = 0.62; middle spotted woodpecker—|r S | = 0.59). Since in fact the condition of the tree fragment is largely affected by tree condition, we analyzed reduced models by including only one member of the most related variables along with all the other explanatory variables. We developed a set of 5 candidate models (Table 1). Candidate models consisted of: the intercept-only model (model 1); two models including tree characteristics with only one member of the two most related variables, i.e. tree condition and condition of tree fragment (models 2–3); a model with a combination of nest-tree characteristics (except for condition of the tree fragment) and distance to forest edge (model 4); and the full model with all explanatory variables (model 5).

Table 1 Candidate model set for the analysis of changes in nest-site selection of great spotted and middle spotted woodpecker in response to dead wood supply

Following the information-theoretic approach (Burnham and Anderson 2002), multimodel inference and model averaging were used to test for differences between the nest sites of the two woodpecker species in pre-outbreak and outbreak periods of dead wood supply, and to test for interspecific differences in nest sites. The relative support for each model was assessed using Akaike’s Information Criterion corrected for small sample size (AICc) (Burnham and Anderson 2002). The subset of models with substantial support (delta AICc <2 compared to the best model) was selected, with their parameter estimates and standard errors averaged (Burnham and Anderson 2002; Bartoń 2009). For all models, Akaike weights were calculated to provide a measure of the probability that a given model was the most strongly supported in the model set, with higher weights indicating better explanatory power. We calculated odds ratios and their 95% confidence intervals (CI) to characterize the effect of explanatory variables (nest-site characteristics) on binary responses (the two study periods or woodpecker species) for model parameters based on model-averaging (Rita and Komonen 2008). In subsequent analyses, we used only complete records for nest-sites. All analyses were carried out in ‘MuMin’ (Bartoń 2009), and ‘lme4’ (Bates et al. 2014) packages using R v.3.3.2 (R Development Core Team 2016).

Results

Changes in woodpecker abundance

The number of great spotted woodpecker pairs varied between 22 in 2001 (1.4 pairs 10 ha−1 of forest area) and 51 in 2015 (3.3 pairs 10 ha−1). The number of middle spotted woodpecker pairs varied between 29 in 2003 and 2009 (1.6 pairs 10 ha−1 of forest area), and 38 in 2000 (2.0 pairs 10 ha−1) (Fig. 1). A simple linear regression analysis suggests that temporal trends in the abundance of both woodpeckers from 2000 to 2015 were not significant (great spotted woodpecker: slope ± SE = 1.009 ± 0.484, t = 2.09, df = 11, p = 0.061; middle spotted woodpecker: slope ± SE = −0.057 ± 0.136, t = −0.42, df = 14, p = 0.680). Quadratic analysis of woodpecker abundance across time did not improve results (both trends were not significant).

Fig. 1
figure 1

Changes in abundance of great spotted and middle spotted woodpeckers in Nature Reserve “Czeszewski Las” in 2000–2015. In 2010–2012, the abundance of the great spotted woodpecker was not assessed. Arrow indicates the start of lethal attacks from oak tree borers and divide study period into two parts—before and after a mass mortality of trees

Changes in standing dead tree availability in pre-outbreak and outbreak periods

The percentage of plots with standing dead trees available for cavity excavation increased from 2.2% (N = 46) in the pre-outbreak period to 54.0% (N = 137) in the outbreak period for the great spotted woodpecker, and from 17.6% (N = 91) to 55.6% (N = 72) for the middle spotted woodpecker, respectively. The density of dead trees increased significantly between the pre-outbreak and outbreak periods for both the great spotted woodpecker (Mann–Whitney test, Z = 6.03, p < 0.0001), and middle spotted woodpecker (Z = 5.09, p < 0.0001). Mean snag density available for the great spotted woodpecker increased by 37-fold, and for middle spotted woodpecker by 4-fold (Fig. 2). In the pre-outbreak period, we found only one dead ash sufficiently large for cavity excavation by the great spotted woodpecker; however, in the outbreak period oaks comprised 56.9%, ashes—25.7%, and other tree species—17.4% of standing dead trees (N = 109). For the middle spotted woodpecker, oaks in the pre-outbreak period represented 33.3%, ashes—37.5%, and other tree species—29.2% of standing dead trees available for cavity excavation (N = 24). In the outbreak period, these species represented 35.5, 40.8, and 23.7% (N = 76) of trees, respectively. The breast height diameter (Dbh) of standing dead trees in pre-outbreak and outbreak periods was 38.2 cm (N = 1) and 42.0 ± 18.6 cm (range 19–108; N = 109) for great spotted woodpeckers. The mean Dbh of dead trees potentially available for middle spotted woodpeckers in successive periods was 30.9 ± 11.2 cm (range 14–54; N = 24) and 36.9 ± 21.3 cm (range 14–95; N = 76), respectively.

Fig. 2
figure 2

Mean density of snags (±SE) sufficiently large for hole excavation for great spotted and middle spotted woodpeckers before and after the mass mortality of trees in Nature Reserve “Czeszewski Las”

GLMMs show that the tree condition affected the utilization of trees by the middle spotted woodpecker irrespective of the study period (Table 2). It was found that the study period (‘pre-outbreak’ and ‘outbreak’), when treated as a random effect, was inert (likelihood ratio tests, p > 0.05). Dead trees were used for nesting significantly more frequently according to their availability.

Table 2 Model coefficients from generalized linear mixed models explaining the utilization of dead trees according to their availability for great spotted and middle spotted woodpeckers

Characteristics of nest-sites

Univariate comparisons showed that some nest site characteristics differed significantly between the pre-outbreak and dead tree outbreak, for both great spotted and middle spotted woodpeckers (Table 3). Both woodpecker species used dead trees twice as often in outbreak as in pre-outbreak periods. Cavities made by great spotted woodpeckers were found in wounds more frequently in the outbreak period. Middle spotted woodpeckers nested significantly more frequently in closed forests in the dead tree outbreak. Other nest-site characteristics did not differ between the pre-outbreak and outbreak of dead trees.

Table 3 Characteristics of cavity placement of great spotted and middle spotted woodpeckers before the outbreak and during the pulse in dead wood supply

Intraspecific variation in nest-site use

For the great spotted woodpecker, the three best models (models 5, 2, and 4; delta AICc <2 compared to the best model) accounted for most of the Akaike model weight (approximately 0.88), indicating they were better supported by the data than the remaining models (Table 4). The first-ranked model, with an Akaike weight of 0.413, was 2.2 times more likely than the third-ranked model, with an Akaike weight of 0.186. The parameters based on model-averaging show that three variables had effect sizes larger than their standard errors (SE), and their 95% confidence intervals (CI) of estimates did not contain the value 0 (Table 5). This suggests tree species other than oak and ash, dead trees, and places with wounds were more important for differentiating between the nest sites of the great spotted woodpecker in both study periods. The odds ratios indicated that in the outbreak period other tree species, dead trees, and tree fragments with wounds were selected 2.89, 4.27, and 2.71 times more often relative to the reference categories (oak, living trees, places without wounds respectively) than in the pre-outbreak period (Table 5).

Table 4 Model-selection statistic for the analysis of the changes in nest-site selection of great spotted and middle spotted woodpeckers in response to dead wood supply
Table 5 Model coefficients and odds ratios of the environmental variables included in the best ranked models explaining changes in nest-site selection of great spotted and middle spotted woodpeckers in response to dead wood supply

For the middle spotted woodpecker, the three best models (models 4, 2, and 5; delta AICc <2 compared to the best model) accounted for most of the Akaike model weight (approximately 0.95) (Table 4). The first-ranked model, with an Akaike weight of 0.418, was 2.31 times more likely than the third-ranked model, with an Akaike weight of 0.181. The parameters based on model-averaging suggest tree species other than oak and ash, and dead trees were more important for differentiating between the nest sites of the middle spotted woodpecker in both study periods. The odds ratio of a nest-site being excavated in other tree species in the outbreak period was 4.17 times less than in the pre-outbreak period. In the outbreak period, dead trees were selected 3.78 times often than living trees than was the case in the pre-outbreak period (Table 5).

Interspecific variation in nest-site use

Of the five models examined, only one (model 4) was well supported by the data, as indicated by both delta AICc value and model weight (Table 6). The model parameters suggest that the probability of a nest-hole having been excavated by the middle spotted woodpecker increased with the nest-hole being in ash, in a dead tree, in a limb/branch, and decreased with a nest-hole being in a closed forest (Table 7). The odds ratios indicated that for the middle spotted woodpecker ashes, dead trees, and limbs/branches were selected 2.63, 3.25, and 2.73 times more likely relative to the reference categories (oak, living tree, limb/branch respectively) than was the case for the great spotted woodpecker. Moreover, the odds ratio for a nest site being in a forest edge was 2.17 times less for the great spotted woodpecker than the odds of the event occurring for the middle spotted woodpecker (Table 7). The interaction terms including characteristics of nest-sites and study periods were not significant.

Table 6 Model-selection statistic for the analysis of niche differentiation between great spotted and middle spotted woodpeckers in response to dead wood supply
Table 7 Model coefficients and odds ratios of the environmental variables included in the best ranked model explaining niche differentiation between great spotted and middle spotted woodpeckers in response to dead wood supply

Discussion

Our results suggest that in mature riverine forests woodpecker abundance does not appear to be limited by the presence of snags. We did not find any evidence that the great spotted and middle spotted woodpeckers tended to increase in abundance over the 16-year period. The density of standing dead trees found in our study area varied greatly between the pre-outbreak and outbreak periods, from a value typical for managed forests to a value close to that reported from primeval lowland forests of temperate Europe, respectively (Walankiewicz et al. 2002). These results are in line with previous studies which examined the role of snags for the great spotted and middle spotted woodpeckers. In two mature oak woodlands in England the great spotted woodpecker showed strong selection for nesting in dead snags (Smith 2007). However, the availability of dead snags did not seem to affect population increase, i.e. the number of breeding pairs increased (Smith 2005), whereas the number of cavities in dead snags remained constant and no trend was discernible over time (Smith 2007). It was found that standing dead trees in the studied forests varied between 6 and 12 and 2 and 6 dead stems ha−1, respectively, with a trend for increasing over a 20-year period (Smith 2007). However, the density of the great spotted woodpecker (individuals km−1) in coniferous stands in Białowieża Forest was dependent on the basal area of dead trees, it being the highest in primeval stands where the basal area of dead snags was 20 times larger than that in the nature reserves and managed stands (Czeszczewik et al. 2013). This suggests that retention of dead trees in coniferous forests may substantially increase habitat quality, i.e. this may be a driver of patch-occupancy dynamics and species density (van Horne 1983; Johnson 2007; Robles and Ciudad 2012), which is generally unsuitable even for a generalist woodpecker species, because of the low number of native deciduous trees suitable for excavation (Barrientos 2010). It has been suggested that with increasing numbers, the great spotted woodpecker is able to exploit living trees as well as dead snags that are available (Smith 2007). This is in line with our results. We found that the proportion of cavities excavated in snags by this species increased two-fold between the pre-outbreak and outbreak periods, despite finding no preference for this type of trees.

For the middle spotted woodpecker, Pasinelli (2000) did not find a significant relationship between the amount of dead wood and the territory size. Lovaty (2002) found high densities of this species in lowland mature oak forest, more than 1.5 territory 10 ha−1 on average, despite all standing dead wood being removed from the study area. Moreover, Stachura-Skierczyńska and Kosiński (2016) did not find a positive relationship between the density of snags and the occurrence and abundance of the middle spotted woodpecker in large-scale studies. Both great spotted and middle spotted woodpeckers show no preferences for foraging on dead trees (Pasinelli and Hegelbach 1997; Robles et al. 2007; Smith 2007), except for in old, near-natural beech forests, where dead wood plays an essential role as a foraging substrate for both species (Hertel 2003). The preferences for dead trees have also been reported during autumn–winter when arthropod availability decreases (Pettersson 1983; Domínguez et al. 2016). Taken together, these studies suggest that the amount of standing dead wood does not affect the abundance of great spotted and middle spotted woodpeckers, irrespective of the spatial scale. On a habitat and landscape scale, middle spotted woodpecker abundance (and density) increases along with the growing number of large living trees and the share of old growth stands (Kosiński and Winiecki 2005; Stachura-Skierczyńska and Kosiński 2016). The most likely mechanism that permits this woodpecker species to increase in density is the reduction of its home range size owing to high levels of availability of large rough-barked trees, mainly oaks, related to potential food abundance (Pasinelli 2000). Earlier generalizations suggest that densities of the great spotted woodpecker increase with increasing proportions of old broad-leaved trees, and with increasing amounts of dead wood (Michalek and Miettinen 2003). However, to our knowledge, there are no quantitative models describing relationships between densities and habitat parameters for this species. In our study area, both woodpecker species reached very high densities in both pre-outbreak and outbreak periods compared to other forest stands (Michalek and Miettinen 2003; Pasinelli 2003; Kosiński and Winiecki 2005; Wesołowski et al. 2015). This suggests that the availability of suitable nest sites does not limit the abundance of either woodpecker species in mature riverine forests. Moreover, previous studies did not find unequivocal evidence of the influence of food supply on the abundance of either species. However, changes in the abundance of middle spotted woodpecker have been positively correlated with winter temperatures (Wesołowski and Tomiałojć 1997). Our results are in line with earlier suggestions that the response of cavity nesters to snag density in hardwood stands seems to be weaker than in coniferous forests (Bunnell 2013). The lack of a positive effect of experimentally increased hardwood snag density, from 15 to 33 snags 10 ha−1, on the abundance of cavity nesters, including four woodpecker species, has been found in a mixed mesophytic forest in North America (McPeek et al. 1987). On the other hand, it has been found that snag creation in managed pine forests, from 3 to 40 snags ha−1 on control and treatment plots respectively, led to an almost 6-fold increase in red-headed woodpecker Melanerpes erythrocephalus populations than on control plots (Kilgo and Vukovich 2014).

Although there were no positive numerical responses from great spotted and middle spotted woodpeckers on the supply of standing dead trees, our data show that standing dead trees were more likely to be used for cavity excavation in the outbreak than pre-outbreak period. Earlier research indicates that woodpeckers appear to select nest sites in dead, soft or decayed parts of trees whose heartwood has been infected by fungi (Jackson and Jackson 2004; Zahner et al. 2012). The preference for dead trees with soft interiors is probably adaptive for woodpeckers, by offering savings in the time and energy required for cavity excavation (Wiebe et al. 2007; Zahner et al. 2012). Presumably, all woodpecker species can identify variations in wood hardness, i.e. they can distinguish between harder sapwood and softer heartwood, both between trees and within the stems (Matsuoka 2008, 2010; Lorenz et al. 2015). Moreover, recent study shows that wood hardness is the most important predictor of nest-site selection at the nest tree and territory scale for six North American woodpecker species inhabiting coniferous forests (Lorenz et al. 2015). Earlier data show that the majority of the nest-holes of great spotted and middle spotted woodpeckers are excavated in limb-holes, open wounds/scars, close to the fruiting body of tree fungus, and in the base of dead or broken limbs (Kosiński et al. 2006; Kosiński and Kempa 2007; Pasinelli 2007; Ćiković et al. 2014). These structures are the primary means by which the fungus gets inside trees causing heart rot, i.e. decay of wood at the center of the trunk and branches. However, trees with characteristics suited for cavity excavation seem to be scarce, even in old growth deciduous forests (Kosiński and Winiecki 2004; Pasinelli 2007; Stachura-Skierczyńska and Kosiński 2016) and coniferous forests (Lorenz et al. 2015). If such structures are rare, then the woodpeckers could benefit from an increasing number of snags, among which at least some should have soft spots suitable for cavity excavation. We did not measure tree hardness; however, since wood hardness decreases from live trees to dead trees (Schepps et al. 1999), we can suppose that dead trees in our study area are more suitable for cavity excavation, and should be preferred over living trees (Kosiński and Winiecki 2004). However, the ability of woodpeckers to excavate in trees of different hardness is related to species’ morphological adaptations to drilling. The neck muscles of the middle spotted woodpecker are less suitable for drilling and drumming, but are more capable of precise lateral probing movements compared with the great spotted woodpecker (Jenni 1981). Thus, an increase in dead tree density may confer a greater benefit to middle spotted woodpeckers than to great spotted woodpeckers. This supposition is confirmed by our data. The proportion of cavities excavated in snags increased from 18 to 40% for the middle spotted woodpecker, reaching the value found in primeval forest (Hebda et al. 2016), and from 5 to 13% for the great spotted woodpecker. Lorenz et al. (2015) suggest that ease of excavation, i.e. preference for sites with soft heartwood, is a major factor driving nest-site selection by woodpeckers in the coniferous forests of northwestern USA. They did not find differences in internal wood hardness of nesting trees; however, they found differences in external (sill) wood hardness (see also Matsuoka 2008). It is likely that the soft external wood of snags is more appropriate for cavity excavation for the middle spotted woodpecker than wounded, decayed or dead places in living trees. This supposition can also be confirmed by the fact that the model including the condition of the tree fragment instead of tree condition received little support. In addition, we found that dead trees were preferred by the species in both pre-outbreak and outbreak periods. Moreover, the presence of wounds and conks of fungi does not significantly affect the nest-site use of the middle spotted woodpecker between study periods. External indicators of wood decay, i.e. presence of the fruiting bodies of polypores, can be fairly unreliable indicators of nest-site availability (Gunn and Hagan 2000; Jusino et al. 2014; Lorenz et al. 2015; but see Pasinelli 2007).

In the outbreak period, tree species other than oak and ash were less likely to be used by middle spotted woodpeckers, and more likely to be used by great spotted woodpeckers compared with the pre-outbreak period. For middle spotted woodpeckers, this was affected by the shift of nest sites into dead trees more suitable for cavity excavation. Among trees other than oaks and ashes, only one tree was classified as dead. Nest niche shifts of great spotted woodpeckers toward tree species other than oak and ash are difficult to explain. Among 21 trees used by this woodpecker species only two were dead. However, most of the cavities (67%) were excavated in living softwood trees, e.g. birch Betula sp., poplar Populus sp., willow Salix sp., and lime Tilia sp. Since great spotted woodpeckers are able to detect variation in wood hardness among trees (Matsuoka 2010), it is likely that these tree species, if available, play an important role in nest site selection (Wesołowski and Tomiałojć 1986; Stenberg 1996; Smith 2007).

Our data show that great spotted woodpeckers were more likely to excavate cavities in wounds in the outbreak than pre-outbreak period. This could be caused by the increasing number of wounds between study periods. Earlier studies have indicated that time is required for the creation of wounds/cracks and holes, and the development of decay processes in living trees (Ranius et al. 2009; Regnery et al. 2013). The sensitivity of trees to injury seems to increase with the age of the tree, due to the higher probability of microhabitat-creating events caused by any number of mechanical factors (Vuidot et al. 2011). The increasing number of fallen dead trees and limbs observed recently probably plays an important role in the creation of wounds. Injuries created in places on a branch that has been broken off are frequently used for cavity excavation by great spotted woodpeckers (Kosiński et al. 2006; Kosiński and Kempa 2007; Pasinelli 2007; Ćiković et al. 2014).

Our study based on long-term data provides information about patterns of niche partitioning between the great spotted and middle spotted woodpeckers in mature riverine forests. We found that the middle spotted woodpecker was more likely to nest in ashes, dead trees, limbs/branches, and less likely to do so in closed forests compared to the great spotted woodpecker. These results expand our earlier knowledge based on univariate analyses (Kosiński and Winiecki 2004; Kosiński et al. 2006). We did not find any significant effects of interaction between nest-site characteristics and the study period which express differences in dead wood availability. Therefore, we suggest that strong increases in standing dead trees did not modify the general pattern of niche partitioning between the great spotted and middle spotted woodpeckers. In contrast to the pre-outbreak period, the frequency of cavities excavated by middle spotted woodpeckers in closed forests increased significantly. This was probably affected by small-scale gap disturbances caused by the death of individual or multiple trees with subsequent falls.

In contrast to our study, Pasinelli (2007) revealed that cavities created by the middle spotted woodpecker were more often freshly excavated, situated close to sporophores, and in dead trees rather than cavities, as is the case for the great spotted woodpecker. However, this study was based on a slightly different set of nest-site attributes, i.e. tree species other than oaks were excluded due to their scarcity, but the position of substrates containing the cavity, and cavity age were included in the analysis. Moreover, Pasinelli (2007) did not include nest location in relation to the forest edge. We did not use cavity age, or the position of the substrate with the cavity in our analysis, since it was not possible to assess this reliably, especially for cavities placed in upper parts of tree crowns. For example, even if cavities seemed to be newly excavated there were no fresh excavation chips on the ground, probably as an effect of dispersion by wind. Moreover, some individuals may begin a cavity that will be completed in following years (Matsuoka 2010; Zahner et al. 2012). Since we did not observe reuse of the same cavities as nest sites in subsequent years by middle spotted woodpeckers, in contrast to a few cases with great spotted woodpeckers (Z. Kosiński unpubl. data; see also Hebda et al. 2016), we can suggest that this species excavates new cavities more frequently than the latter species. This tendency can be enhanced by reuse of the middle spotted woodpecker cavities by the great spotted woodpecker, and nest usurpation by common starlings Sturnus vulgaris (Mazgajski 2003; Kosiński and Ksit 2006; Smith 2006). Thus, the mechanism forcing middle spotted woodpeckers to excavate new holes may be motivated by interspecific competition. Moreover, it has also been found that the lifespan of cavities excavated by woodpeckers in dead trees or dead part of trees is much shorter than in trees with living sapwood (Wesołowski 2011); cavities excavated by middle spotted woodpeckers persisted for a shorter time than those of great spotted woodpeckers, even if excavated in the same type of substrate. All of these results suggest that the middle spotted woodpecker preference for new cavities might be a by-product of interspecific competition, and the shorter lifespan of their cavities.

In our study area, the presence of conks of tree fungi near the cavity did not affect niche differentiation of great spotted and middle spotted woodpeckers. As stated earlier, the presence of conks can be an unreliable indicator of tree hardness (Jusino et al. 2014; Lorenz et al. 2015). Fungi that enter the tree center through wound sites are less destructive for the wood body than heart rot fungi that enter the tree center through old broken limbs which are most frequently seen in standing dead trees (Butin 1983; Zahner et al. 2012). Moreover, primary excavators probably do not select specifically for rot or sporophores, but rather they select trees with considerable wood hardness variation, and soft interiors (Matsuoka 2010; Lorenz et al. 2015).

Standing dead trees have been recognized as an optimal habitat for primary cavity excavators, particularly woodpeckers (Winkler and Christie 2002; Cockle et al. 2011). Our findings show that although an increase in the standing dead wood availability does not affect woodpecker population growth in mature riverine forests, it seems to provide habitats that are more favorable, especially for weak excavators, i.e. the middle spotted woodpecker. Snags and cavities are only suitable for nesting for a few years after their creation (Smith 2007; Wesołowski 2011). Therefore, responsible forest management should consider the necessity to leave standing dead wood, dying trees, and trees with low commercial value, and dead wood dynamics, e.g. the rate between the creation and removal/collapse of snags (Smith 2007). Earlier data suggest that it is possible to simultaneously preserve suitable habitat conditions for middle spotted woodpecker populations and maintain the commercial exploitation of a forest with a high economic return when dead wood availability is relatively low, i.e. 5 dead trees ha−1 (Stachura-Skierczyńska and Kosiński 2016; see also Lõhmus et al. 2010). Although large-scale dieback of common ashes may have a negative impact on forest ecosystems throughout temperate Europe (Pautasso et al. 2013), our data suggest that this process may locally provide habitats more favorable for middle spotted woodpecker, e.g. in riverine forests where the common ash occurs as an important admixed tree species. However, the decreasing number of ashes and large oaks, and the lack of natural regeneration of the latter species, may negatively affect the middle spotted woodpecker in the future.

Conclusions

The abundance of great spotted and middle spotted woodpeckers did not increase between 2000 and 2015. This result suggests that the abundances of both species are not affected by an increase in the density of standing dead trees. It is likely that living trees with characteristics suited for cavity excavation preferred by woodpeckers exist in managed mature riverine forests well above the given critical threshold and, therefore, do not limit woodpecker abundance. However, our study demonstrates that differences in standing dead tree availability affect intraspecific nest-site use. The most important factor is that both species, i.e. the generalist forest dweller, the great spotted woodpecker, and the specialist forest dweller, the middle spotted woodpecker, were significantly more likely to use snags in outbreak periods than in pre-outbreak periods. This increase was most visible for the middle spotted woodpecker, i.e. the species with weaker excavating abilities, which clearly prefers nesting in dead trees. An increase in the number of dead trees used by the great spotted woodpeckers should be treated as a random process, as this species did not select snags more frequently according to their availability. This type of substrate probably reduces the energy expenditure of woodpeckers and may be adaptive by offering time and energy savings. However, no study into the fitness consequences of preferred dead trees has yet been conducted and it needs to be determined whether excavating in dead and soft wood is helpful for woodpeckers, e.g. in terms of laying dates and clutch size. It is likely that excavation decisions among the great spotted woodpecker involve the selection of trees with a harder sapwood and a softer heartwood to optimize the energy expenditure in the cavity excavation and the need for nest security.