Drivers of nest survival rate in a southern Tunisian population of Laughing Doves(Spilopelia senegalensis)

Jihen Boukhriss and Slaheddine Selmi

Abstract

Keywords: Laughing Dove, Laying date, Nest age, Nest microhabitat, Nest survival, Oasis, Tunisia

Background

Nest survival is a major component of breeding success in birds, with direct consequences on their population dynamics (Cresswell 2011). The ornithological literature often emphasizes that nest predation is the major cause of bird nest failure (Ricklefs 1969; Martin 1992). It also points out that the risk of nest predation is particularly affected by habitat parameters, through their effects on nest detectability, accessibility, and predator abundance (Martin and Roper 1988; Filliater et al. 1994; Latif et al. 2012; Guan et al. 2018). Nest predation risk is also known to vary with nesting date due to seasonal variation in predator abundance and activity (Schaub et al. 1992;Roos 2002; Cox et al. 2012; Ma et al. 2019). In addition,nest age (i.e. days after clutch initiation) is often mentioned as one of the key factors determining the risk of predation, through its possible effects on the probability that the nest is detected by predators (Schaub et al.1992; Sockman 1997; Roos 2002). Investigating the relative importance of these factors in shaping nest survival is of great importance for assessing the breeding success of bird populations and for understanding their spatiotemporal dynamics.

Tunisian oases are semi-natural wood lots surrounded by an inhospitable desert environment, and host an original diversity of migratory and resident birds (Selmi 2000). Although the Palearctic character of this avifauna is clear (Isenmann et al. 2005), an Afro-tropical species,namely the Laughing Dove (Spilopelia senegalensis), has successfully colonized these oases and proliferated (Selmi 2000; Isenmann et al. 2005). The breeding ecology of this expanding Afro-tropical bird in the Tunisian oasis habitat has previously been described by Boukhriss and Selmi(2009) who reported a modal clutch size of two eggs(more than 95% of studied clutches) and an estimated chick productivity of one chick per breeding attempt.However, the key factors affecting nest survival have not been investigated.

The aim of this study was thus to investigate key factors affecting nest survival in a southern Tunisian population of Laughing Doves. More specifically, we assessed the relevance of nest age, laying date, nest height and vegetation structure around the nest tree as possible predictors of daily nest survival rate. Previous studies have shown that bird nest predators in the oasis habitat included snakes, Semi-domestic Cats (Felis domestica),Gennet (Genetta genetta), but the most notable impact seemed to be that of the Black Rat (Rattus rattus), which is a non-native rodent species in southern Tunisia Ruffino and Vidal (2010). This arboreal rodent is abundant in Tunisian oases and seems to more frequently prey on bird eggs than on nestlings (Selmi 2004; Boukhriss et al.2009; Boukhriss and Selmi 2018). Because the risk of nest predation is often expected to decrease with increasing concealment and decreasing accessibility (Colombelli-Négrel and Kleindorfer 2009), daily nest survival rate was hypothesized to increase with increasing vegetation density around the nest tree, and also with increasing nest height. In addition, since parental activity around the nest is generally assumed to be higher during the chick rearing phase than during egg incubation (Skutch 1949;Roper and Goldstein 1997; Martin et al. 2000), the risk of the nest being detected and attacked by predators was expected to increase as the nesting attempt progressed.A negative relationship between age and nest survival was therefore expected. Nest survival rate could also be expected to vary across the breeding season, mainly because of possible seasonal changes in nesting microhabitat, and hence nest concealment (Winter et al. 2005;Ludwig et al. 2012), and in the activity of nest predators(Sperry et al. 2008; Husby and Hoset 2018).

Methods

Study area and field procedure

This work was carried out in the oasis of Kettana, in the south-east of Tunisia (33° 45′N, 10° 13′E). This palm grove covers an area of 464 ha and is composed of a mosaic of small private fields where farmers practice a rural agricultural activity based on irrigation and crop diversification. Within this oasis the plant cover is organized into three main layers: herbaceous plants (cultivated and spontaneous), fruit trees, mainly olive trees (Olea europaea) and pomegranate trees (Punica granatum),and date palm trees (Phoenix dactylifera).

From early March to mid-August 2005, nests were searched for by checking potential nest sites available in the area, while paying attention to the behaviour of the birds that can sometimes inform about the location of the nest. When an active nest was found, we immediately checked its content and tried to determine nest age and laying date. Nest age was precisely recorded for nests whose laying date was known, i.e. nests discovered before the laying of the first egg. However, for nests discovered during the incubation or chick rearing stages, nest age at the moment of discovery could not be determined with precision. We thus relied on the aspect of eggs or chicks and we used descriptive accounts from known-aged nests in the same population to determine approximately nest age and laying date. Then the nest was marked by placing one coloured plastic band at 1 m from the nest tree and was regularly checked until failure or fledging. Nests were monitored every 5 days during incubation and daily after eggs hatched until the end of the breeding attempt.To avoid disturbance, nests were usually observed from a distance to ensure that they were empty before approaching them. Nests were considered successful if at least one young fledged. Following Manolis et al. (2000), nest exposure was determined as the number of days since nest discovery until the halfway point between the last day that the nest was known to be active and the day that the breeding attempt was observed to have finished (fledging or loss).

After the nest was inactive (i.e. fledging or failure),we measured nest height (m), and we assessed visually the covers (%) of palm trees, fruit trees and herbaceous plants in a 5 m-radius circle. These vegetation variables provided descriptors of habitat structure in the close nest tree environment. We are aware that the possible change in foliage density between the day of nest discovery and the end of the nesting attempt may have introduced some bias in our data, but we did so to avoid disturbing birds and attracting predators to the nest. All measurements were usually conducted by the same observer (J. Boukhriss) to avoid a possible observer effect.

Data analyses

We first conducted a principal component analysis (PCA)on the three assessed vegetation variables (covers of palm trees, fruit trees and herbaceous plants) to summarize them into a composite index of vegetation structure at the close nest environment. This PCA was conducted using the FACTOR procedure in SAS software (1998).

Nest monitoring data were first used to estimate daily nest survival rate using the Mayfield method, which accounts for heterogeneity in exposure period among nests (Mayfield 1975). The rate of nesting success was then estimated by raising the obtained estimate to the power of the number of days of nesting (Mayfield 1975),which is habitually 28 days in the studied Laughing Dove population (14 days of egg incubation and 14 days of nestling rearing; personal data).

We also investigated the relevance of nest height, vegetation structure in the close nest tree environment (i.e.a 5 m-radius circle around the nest tree), nest age (days after clutch initiation) and laying date (Julian date of the first egg laying) as predictors of nest daily survival rate, by using the modelling approach of Shaffer (2004).We used the GENMOD procedure in SAS (1998) to fit logistic exposure models with a binomial error distribution and a logit-link function. We constructed 16 candidate models, including a constant model (i.e.,without covariates) and models considering each of the four covariates by itself or in association with one,two or three other covariates. The candidate models were then ranked according to their AICc scores (Akaike’s Information Criterion corrected for small sample size), with lowest AICc indicating the best approximating model for the data (Burnham and Anderson 2002).Models with ΔAICc ＞ 4 were considered to have little to no empirical support (Burnham and Anderson 2002).We also took into account model selection uncertainty by using the model-averaging technique (Burnham and Anderson 2002) for the assessment of the effects of covariates on daily nest survival rate. The relative importance of each variable was obtained by summing up the weights (W) of all models including that variable(Burnham and Anderson 2002), and its averaged effect(and associated 95% confidence interval) was estimated following the formula: βaverage= ΣWiβi, where Wiis the weight of model i and βiis the estimated effect of the considered covariate according to model i (Burnham and Anderson 2002). The effect of a given covariate on daily nest survival rate was considered as significant at the 5% level if the 95% confidence interval of its βaveragedid not overlap zero.

Finally, as we were interested in exploring more the relationship between nest age and daily survival rate, we sought to investigate how daily nest survival rate varied between the pre-hatching (incubation) and post-hatching(chick rearing) stages. To do so, we subdivided our nest sample into two sub-samples: (1) nests discovered during the pre-hatching stage and (2) nests discovered during the post-hatching stage. For the pre-hatching sub-sample, the clutches that succeeded to hatch eggs were considered successful, while the success corresponded to the fledging of chicks for the post-hatching sub-sample. In other words, we estimated the rate of clutch hatching and the rate of chick fledging in hatched clutches. The estimations were conducted using the Mayfield method (1975),and the comparison between the obtained estimates was conducted by means of a z-test following the procedure described in Johnson and Shaffer (1990).

Results

In total, 120 Laughing Dove nests were found and monitored. The exposure period varied among nests, from 2 to 35 days, with an average (± SE) of 13 ± 2 days. The earliest egg laying occurred on March 16, while the latest one was on August 8. The age of the nest at its discovery varied between 4 and 28 days, with an average (± SE) of 22 ± 0.78 days.

Most nests were placed on pomegranate trees (55%)and olive trees (40%). The remaining 5% of nests were placed on four species of fruit trees (apricot, pear, mulberry and grape). Nest height ranged from 0.80 to 5.98 m,with an average (± SE) of 2.59 ± 0.09 m. With regard to vegetation structure around the nest tree, the PCA conducted on the three assessed vegetation variables (covers of palm trees, fruit trees and herbaceous plants) summarized them into one factor (first component) providing a composite index of vegetation structure at the close nest environment. This factor had an eigenvalue of 1.41 and accounted for 47% of the variance in the original data,while the second component derived from the PCA had an eigenvalue less than 1 (0.90) and was not retained.The index of vegetation structure derived from the PCA was positively correlated with the covers of palm trees(r = 0.58, p ＜ 0.0001) and herbaceous plants (r = 0.69,p ＜ 0.0001) and negatively correlated with the cover of fruit trees (r = - 0.78, p ＜ 0.0001). Low scores characterized dense plantations of fruit trees where palm trees and herbaceous plants were rare. However, high scores characterized plots where fruit trees were replaced by palm trees, giving a more open appearance and allowing the development of an important herbaceous cover.

Among the 120 nests monitored, 40 were depredated,while 80 were successful to raise at least one fledgling.This gave an apparent nest success rate (i.e., proportion of successful nests) of 67% (95% CI 58-75%). However,when the heterogeneity in exposure period among monitored nests was accounted for, by applying the Mayfield method, the daily nest survival rate was estimated at 0.975 (95% CI 0.968-0.983), corresponding to an estimated rate of nesting success of 49% (95% CI 40-62%).

Results of logistic exposure regressions showed that among the 16 competing models, four models were supported by the data and provided parsimonious models in explaining variation in daily survival rates among the studied nests (Table 1). Together, these highly-ranked models summed up a weight exceeding 99%, while the remaining models had negligible weights and could not be considered.

All four selected models included the effects of vegetation structure around the nest tree and nest age (Table 1),suggesting that these two parameters were the best predictors of daily nest survival rate in our studied population. Moreover, these two parameters showed weights exceeding 99%, and were the only investigated variables to have estimated effects whose confidence intervals did not include zero (Table 2). Nest daily survival rate was negatively related to the index of vegetation structure(Table 2; Fig. 1), indicating that the presence of palm trees in the close nest tree environment (i.e. high score of vegetation structure index) was associated with reducedchance of nesting success and increased risk of predation.It was however positively related to nest age (Table 2),suggesting that the risk of nest failure decreased as the nesting attempt progressed, throughout the egg incubation and nestling stages. The latter suggestion was also supported by the comparison of the estimated daily survival rates of nests in the two breeding stages (Fig. 2;z = 5.376, p ＜ 0.0001).

Table 1 Results of model selection for daily nest survival rate as a function of nest age, nesting date, nest height and vegetation structure in the close nest environment

Table 2 Model-averaged parameter estimates and associated 95% confidence intervals derived from the 16 candidate logistic-exposure models for daily nest survival rate as a function of nest age, laying date, nest height and vegetation structure in the close nest environment

Fig. 1 Plot of predicted daily nest survival rate as a function of vegetation structure in the close nest tree environment. Dashed lines represent the limits of 95% confidence interval for the estimated daily nest survival rates

Discussion

Fig. 2 Comparison of the estimated daily nest survival rates between the pre-hatching (egg incubation) and post-hatching (chick rearing)stages. Numbers on the bars indicate sample sizes

In this study we investigated factors shaping the survival of Laughing Dove nests in a southern Tunisian oasis habitat. By doing so, we mainly aimed to provide information helping to understand the breeding success and the expansion of this afro-tropical bird species in a North African area. Our results showed that the estimated daily survival rate of Laughing Dove nests in the studied oasis was similar to those previously estimated for some Palearctic passerine species nesting in the same oasis habitat, i.e. the Rufous Bush Robin (0.970) and Woodschat Shrike (0.980) (Boukhriss and Selmi 2018).It was also equal to that reported for the same dove species by Brahmia et al. (2015) in an agricultural landscape in northern Algeria (0.975). We also found that vegetation structure in the close nest tree environment and nest age provided important predictors of daily nest survival rate in the studied Laughing Dove population, whereas neither nest height nor laying date showed significant effects.

Our results showed that the presence of palm trees near the nest tree was associated with decreased daily nest survival rate. Dense plantations of fruit trees offered safer nesting sites compared to plots with a mixture of fruit trees and palm trees. This result could be explained by increased abundance and/or activity of nest predators around palm trees compared to fruit tree plantations. Palm trees seemed to provide refuges for small nest predators, increasing the probability of the nests to be detected and depredated. Indeed, with their multiple suckers, palm trees often form dense tufts serving as refuges for wild mammals, notably the Black Rat which seemed to be the main nest predator in the oasis habitat(Selmi 2004; Boukhriss et al. 2009; Boukhriss and Selmi 2018). Overall, this explanation is consistent with the general trend of decreased bird nest survival in preferred habitats for nest predators (Johnson and Temple 1990;Dion et al. 2000; Whittingham and Evans 2004; Hammond 2016). Moreover, given the arboreal habits of this rodent and the low height at which dove nests were built,the lack of relationship between nest survival and height is no longer surprising.

With regard to nest age, our results showed that daily nest survival rate increased as the nesting attempt progressed, a pattern that has also been reported in passerines breeding in the same habitat (Boukhriss and Selmi 2018) and in other different habitats (e.g., Martin 1992;Roper and Goldstein 1997). Nest survival was higher during the nestling period compared to the egg incubation period, which is contrary to the widespread opinion of increased nest mortality during the post-hatching period(Skutch 1949; Roper and Goldstein 1997; Martin et al.2000).

Several non-exclusive hypotheses could be proposed to explain the positive relationship of daily nest survival with nest age. First, the Black Rat seemed to more frequently consume bird eggs than nestlings in the oasis habitat (Boukhriss et al. 2009; Boukhriss and Selmi 2018). Its impact on bird nests would thus be higher during egg incubation compared to the nestling stage. An alternative explanation could be the possible increase in parental defence behaviour with the progress in the nesting attempt. Indeed, birds are known to take higher risk against predators and to defend more intensely their nests as the breeding attempt progresses and its reproductive value increases (Montgomerie and Weatherhead 1988; Brunton 1990; Boukhriss and Selmi 2010). We thus hypothesise that the observed increase in daily nest survival rate with nest age could be at least partly due to the increase in the intensity of nest defence behaviour exhibited by parents with the progress in the nesting attempt.Lastly, the positive relationship between nest age and nest survival rate may be a simple artifact. Indeed, nests easily detectable by predators (i.e. less concealed and more accessible nests) are more likely to be eliminated early (i.e. during egg laying or early incubation), resulting in a low estimate of the daily nest survival rate. On the other hand, nests arriving at an advanced stage would be nests that had escaped predation because they had a low probability of being detected. The last nests would necessarily have a high probability of survival.

Conclusions

In conclusion, our results suggested that the presence of palm trees in the close environment of the nest tree was associated with increased predation risk. Clumps of palm tree suckers seemed to provide refuges for small nest predators, notably the Black Rat, which has been reported to be the main nest predator in the oasis habitat. The predatory activity of this rodent seemed more directed against eggs than nestlings, which may explain the observed increase in daily survival rate with nest age.Investigations of the ecology and feeding behavior of this nest predator in the oasis habitat are however needed for firmer conclusions about these issues.

Acknowledgements

We thank H. Belwafi and A. Belwafi who helped a lot in nest finding.

Authors’ contributions

JB and SS conceived and designed the study. JB performed the field work. JB and SS processed and analyzed the data. JB wrote the manuscript. SS made comments to the manuscript. Both authors read and approved the final manuscript.

Funding

The work was not funded by any specific body.

Availability of data and materials

The datasets used in the present study are available from the corresponding author on reasonable request.

Ethics approval and consent to participate

This work complies with the current laws of Tunisia. It was based on simple field observation without any experimental manipulation or prejudice to the studied birds.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Received: 8 July 2019 Accepted: 4 November 2019