Introduction
Seabirds are often abundant on islands (Mulder et al., Reference Mulder, Jones, Kameda, Palmborg, Schmidt, Ellis, Mulder, Anderson, Towns and Bellingham2011). However, colonial nesting habits, low reproductive output and extended periods of parental care at breeding grounds make seabirds vulnerable to predators (here we use the term predator to describe only non-native mammalian predator species; Furness & Camphuysen, Reference Furness and Camphuysen1997; Baillie et al., Reference Baillie, Hilton-Taylor and Stuart2004; Wolf et al., Reference Wolf, Keitt, Aguirre-Muñoz, Tershy, Palacios and Croll2006; Croxall et al., Reference Croxall, Butchart, Lascelles, Stattersfield, Sullivan, Symes and Taylor2012). Introduced predators are among the most acute and widespread threat to seabirds; for example, rats Rattus spp., which depredate seabirds at every life stage, have invaded at least 80% of island groups (Towns et al., Reference Towns, Atkinson and Daugherty2006; Jones et al., Reference Jones, Tershy, Zavaleta, Croll, Keitt, Finkelstein and Howald2008; Varnham, Reference Varnham2010; Croxall et al., Reference Croxall, Butchart, Lascelles, Stattersfield, Sullivan, Symes and Taylor2012; Spatz et al., Reference Spatz, Newton, Heinz, Tershy, Holmes, Butchart and Croll2014). Consequently, seabirds are one of the most globally threatened groups of animals. Of c. 365 species worldwide, 29% have been categorized as globally threatened, and 10% as Near Threatened on the IUCN Red List; 5% are categorized as Critically Endangered (IUCN, 2015; BirdLife International, 2015). Removal of predators is considered to be one of the most effective strategies for seabird conservation (Jones, Reference Jones2010; Towns et al., Reference Towns, West and Broome2013; Jones et al., Reference Jones, Holmes, Butchart, Tershy, Kappes and Corkery2016), and can have positive feedbacks for many biota and ecosystem processes affected by predator presence (Towns et al., Reference Towns, Atkinson and Daugherty2006; Towns, Reference Towns2009; Lavers et al., Reference Lavers, Wilcox and Donlan2010; Jones & Kress, Reference Jones and Kress2012; Le Corre et al., Reference Le Corre, Danckwerts, Ringler, Bastien, Orlowski and Rubio2015; Jones et al., Reference Jones, Holmes, Butchart, Tershy, Kappes and Corkery2016).
Despite considerable investment in the removal of introduced predators from islands, the way seabirds respond subsequently remains poorly understood (but see Lavers et al., Reference Lavers, Wilcox and Donlan2010; Buxton et al., Reference Buxton, Jones, Moller and Towns2014, Reference Buxton, Taylor, Jones, Lyver, Moller, Cree and Towns2016). Available research has described variable responses among species and islands (Gaze, Reference Gaze2000; Lavers et al., Reference Lavers, Wilcox and Donlan2010; Ismar et al., Reference Ismar, Baird, Gaskin, Taylor, Tennyson and Rayner2014). The same traits that make seabirds vulnerable to predation can also inhibit their natural recovery following the removal of predators (Jones et al., Reference Jones, Towns, Bodey, Miskelly, Ellis, Rauzon, Mulder, Anderson, Towns and Bellingham2011; Buxton et al., Reference Buxton, Jones, Moller and Towns2014). Species-specific differences in reproductive output, philopatry and behavioural characteristics influence seabirds’ recovery and recolonization of newly available habitat (Danchin et al., Reference Danchin, Boulinier and Massot1998; Jones et al., Reference Jones, Towns, Bodey, Miskelly, Ellis, Rauzon, Mulder, Anderson, Towns and Bellingham2011; Buxton et al., Reference Buxton, Jones, Moller and Towns2014); for example, there was no observed increase in populations of flesh-footed shearwaters Ardenna carneipes or sooty shearwaters Puffinus griseus following the removal of Norway rats Rattus norvegicus from Titi Island, Marlborough Sounds, in the 1970s (Gaze, Reference Gaze2000). Similarly, on Raoul Island in the Kermadec Islands there was no observed response of white-naped petrels Pterodroma cervicalis following the removal of Norway rats, kiore Rattus exulans, and cats Felis catus during 2002–2004 (Veitch et al., Reference Veitch, Gaskin, Baird, Ismar, Veitch, Clout and Towns2011). Conversely, on the same island there was a notable increase in nesting by black-winged petrels Pterodroma nigripennis and wedge-tailed shearwaters Ardenna pacifica, and recolonization by Kermadec petrels Pterodroma neglecta, Kermadec little shearwaters Puffinus assimilis kermadecensis and Kermadec storm petrels Pelagodroma albiclunis (Gaskin, Reference Gaskin2012; C. Gaskin, unpubl. data).
Twenty-five species of invasive predators of seabirds have been eradicated successfully from 1,016 islands worldwide (Island Conservation and Invasive Species Specialist Group, 2015; Jones et al., Reference Jones, Holmes, Butchart, Tershy, Kappes and Corkery2016). Responses of seabirds and island ecosystems to these eradications are now being evaluated to identify how seabird colonies recover following predator removals, and whether they can recover to a fully restored state (Jones, Reference Jones2010). Understanding the recovery of seabirds through space and time following predator eradication is crucial for informing conservation management of seabirds, particularly for species experiencing population declines exacerbated by threats such as fisheries bycatch, marine plastic pollution and climate change (Rolland et al., Reference Rolland, Nevoux, Barbraud and Weimerskirch2009; Croxall et al., Reference Croxall, Butchart, Lascelles, Stattersfield, Sullivan, Symes and Taylor2012).
The equilibrium theory of island biogeography states that species diversity maintains a dynamic equilibrium over time, influenced by area and ecological diversity (Preston, Reference Preston1962; MacArthur & Wilson, Reference MacArthur and Wilson1967; Simberloff, Reference Simberloff1974). We would therefore expect seabird diversity on islands to follow this pattern. By comparing islands cleared of introduced predators to islands that have not been affected by the presence of predators or habitat modification we may be able to discern the state of recovery following predator eradication. Here, with a basis of island biogeography theory, we synthesize key influences on seabird recolonization and recovery reported in the literature. We present a case study to test bigeographical influences by examining whether there has been measurable change in seabird assemblages on a sample of islands in New Zealand that have been cleared of invasive predators. We discuss how these biogeographical influences and additional factors of behaviour, demography and within-island habitat availability may be driving seabird recovery. Ideally, pre-eradication population census data would be used for comparative evaluations of species’ responses to the removal of threats, but in practice such data are rarely available (Duffy, Reference Duffy, Nettleship, Burger and Gochfeld1994). We therefore evaluate changes in seabird richness on islands that have been cleared of invasive predators, and compare them to a subset of islands of comparable area that remained free of predators and for which historical seabird species composition data are available. We focus on the Hauraki Gulf, in northern New Zealand, which has 27 species of breeding seabirds (Gaskin & Rayner, Reference Gaskin and Rayner2013). We used survey data of seabird richness on islands where species were confirmed to be breeding as of the Austral winter of 2015.
Study area
The islands of the Hauraki Gulf are bounded by the Poor Knights to the north, the Great Barrier Island group to the east and the Aldermen (Ruamaahua) Islands to the south (Fig. 1). The islands included in our study represent a range of ecological states, from extensively modified by invasive mammals, either through previous or current influences (e.g. Rakitu Island), to near pristine (e.g. most islands in the Poor Knights group; Supplementary Material 1; Supplementary Table S1). We excluded islands of < 1 ha because of the small island effect (Burns et al., Reference Burns, McHardy and Pledger2009), and because the survey data on predator presence and seabird populations may be unreliable because of survey constraints (Supplementary Material 2). We also excluded islands that have a permanent residential population because of the potential effects of human activities (e.g. land use, presence of domestic animals, light pollution; Le Corre et al., Reference Le Corre, Ollivier, Ribes and Jouventin2002) on the establishment of seabird colonies. We grouped islands (n = 98) into three categories based on predator status: uninvaded (n = 30), invaded (n = 37) and cleared (n = 31; Fig. 1; Supplementary Table S4). The climate of the Hauraki Gulf is temperate–humid, with a mean annual temperature of 16°C, relative humidity of 87% and mean rainfall of 1,202 mm (CliFlo, 2016).
Methods
Data
Long-term survey data on species confirmed to be breeding on the islands as of the Austral winter of 2015 are from C. Gaskin (unpubl. data) and D.R. Towns (unpubl. data), supplemented by peer-reviewed literature (Fig. 1; Supplementary Table S4). The data are presence only. Absence data are rarely available because of the more labour intensive sampling methods and greater economic costs associated with identifying the absence of a species or community (Phillips et al., Reference Phillips, Anderson and Schapire2006). In most cases the data are based on observations from up to 30 years ago, with recent presence (i.e. 2004 to present) confirmed through acoustic surveys, occasional island visits and, in some cases, targeted species surveys (e.g. flesh-footed shearwater, Buller's shearwater Ardenna bulleri, grey-faced petrel Pterodroma gouldi, black petrel Procellaria parkinsoni, Cook's petrel Pterodroma cookii, New Zealand storm petrel Fregetta maoriana; C. Gaskin, unpubl. data) where additional data on other species have been collected. These reports are based on observations; impressions can vary in quality, depending on the methodology, taxon, and skill of the observer (Supplementary Material 2). These biases in sampling effort may influence our results. However, the majority of data were collected by the same group of researchers and any potential biases are likely to be consistent across our study area.
We followed BirdLife International (2015) for taxonomy and nomenclature (Supplementary Material 3). We included 16 seabird species in the study: 14 Procellariiformes, one Sphenisciform and one Pelecaniform (Supplementary Table S2). All exhibit natal site philopatry, and they have a range of population-level responses to introduced predators. Terns (Sternidae), gulls (Laridae) and shags and cormorants (Phalacrocoracidae) were discarded from the analysis because they exhibit ephemeral breeding site selection behaviour (Monaghan, Reference Monaghan1996; New Zealand Birds Online, 2015). We have generalized assumptions about behavioural traits, such as sex biases and habitat selection preferences.
We included rats (Rattus rattus, R. novegicus, R. exulans), mice Mus musculus, cats, and pigs Sus scrofa (Harris, Reference Harris1970; Medway, Reference Medway2001) as known predators of seabirds. Predator eradication data were obtained from the Department of Conservation, D.R. Towns (unpubl. data), Gaskin et al. (Reference Gaskin, Fitzgerald, Cameron and Heiss-Dunlop2011) and the Database of Island Invasive Species Eradications (Island Conservation and Invasive Species Specialist Group, 2015). Eradications were carried out over 77 years, from 1936 (pigs from Aorangi) to 2011 (R. exulans from Taranga). We used the date of the last successful eradication, where predators had not been identified in biosecurity follow-ups (Supplementary Table S4).
Given the paucity of pre-eradication seabird species composition data for cleared islands, temporal analysis of seabird recovery is challenging. We therefore used the available census data for uninvaded islands as a basis for comparisons with seabird species richness of cleared islands (i.e. this is essentially a natural experiment; Oksanen, Reference Oksanen2001). This approach simplified potential ecological complications, as we assumed that species richness for a given area would be comparable between islands, given the similar geology and habitat of the islands and the proximity of seabird populations. To assess the stability of seabird assemblages we used historical survey data for six uninvaded ‘reference’ islands (Table 1), with the earliest data collected in 1928. Although census surveys on the reference islands were not systematic, these data provide information on the observed stability of seabird species composition over the survey period. Species were included on the list only if they were recorded as breeding at that location. Disturbance on these islands varied from historical land clearance by burning, and agricultural use by Māori, to minimal use (Table 1). The reference islands were of pristine or outstanding quality (class I and II; Taylor, Reference Taylor1989; Supplementary Table S1), and therefore we assumed that there was no limitation of available breeding habitat for seabirds.
*McCallum (Reference McCallum1981) noted that seabirds that were recorded on the island previously were not present on his visit in 1981 because it did not coincide with the breeding times of these birds. (1981 not confirmed breeding).
Analysis
All statistical analyses were carried out in R v. 3.1 (R Development Core Team, 2013). We used generalized linear models with Poisson-distributed errors to estimate the species–area relationship by predator status for all islands (Supplementary Material 4). The models of species richness and island area assume there is a linear increase in the number of species as island area increases but in reality this is limited to the number of seabird species in the region. We assume that the maximum species richness is constrained by the number of species included in the analysis (n = 16). Model assumptions of a linear mean variance relationship and of spatially uncorrelated errors were checked using the dispersion test of the AER package in R (Kleiber & Zeileis, Reference Kleiber and Zeileis2008), and spatial correlation tests of the sp, ape and gstat packages (Supplementary Material 4; Paradis et al., Reference Paradis, Claude and Strimmer2004; Pebesma, Reference Pebesma2004; Bivand et al., Reference Bivand, Pebesma and Gomez-Rubio2013).
We considered three models for the species–area relationship by predator status, all being Poisson regressions of a semilog model (i.e. species count, S, as response, log10(A) as predictor). The three models are all Poisson generalized linear models: S ~ Poisson(μ) with log link, such that log μ(S) = βX i , where X and β are the design matrix and the parameter vector, respectively.
Effect of log(area) only:
Additive effects of log(area) and predator status (i.e. separate intercept per predator status but same slope):
Interaction of log(area) and predator status (i.e. separate slope and intercept per predator status):
Multi-model inference based on Poisson generalized linear models was used to explore possible predictors of species richness for the cleared islands (Burnham & Anderson, Reference Burnham and Anderson2002). The predictors explored were size, time since eradication, distance to the mainland (i.e. the nearest point on the coastline of the North Island) and distance to Auckland city centre (36.85°S, 174.76°E), the latter two being potential proxies for human disturbance and/or distance to offshore feeding grounds. Distances were calculated using the rgeos package (Bivand & Rundel, Reference Bivand and Rundel2016). We used the MuMIn package (Bartoń, Reference Bartoń2015) to generate a complete set of candidate models based on the above predictors, and ranked the resulting model fits by Akaike information criterion corrected for a small sample size (AICc).
We used χ2 tests (α = 0.05) of frequency of occurrence to evaluate the differences in species present on cleared, uninvaded and invaded islands (R Development Core Team, 2013). Given the poor representation or absence of some species within island categories, analyses used a subset of six widely distributed species (Fig. 3).
We used a Kruskal–Wallis rank sum test to test for differences in the mean age at first reproduction, a possible proxy for colonization ability, for each island assemblage. Species-specific values (Supplementary Table S2) were averaged over all breeding species on a given island.
Results
The mean area of uninvaded islands was 13.88 ± SD 31.43 ha (range 1.04–155.92 ha), biased towards < 10 ha (23 of 30 islands). The mean area of invaded islands was 24.79 ± SD 62.34 ha (range 1.14–312.33 ha). The mean area of cleared islands was 286.15 ± SD 675.7 ha (range 1–2,817 ha), and represented the most evenly spread area range of all three categories.
The final model for the species–area relationship included separate slopes and intercepts by predator status. No significant overdispersion was detected (overdispersion test: dispersion parameter = 0.99; P = 0.53). Moran's I (I = 0.1, P = 0.011) indicated weak but statistically significant spatial autocorrelation in the residuals, and visual inspection of spatial residuals indicated that this was because the model overpredicted species richness for the islands of the inner Hauraki Gulf (Supplementary Fig. S1). Generalized linear model parameter estimates are in Table 2, and model predictions are illustrated in Fig. 2.
For uninvaded islands the intercept (i.e. expected species count on a 1 ha island) is 2.43 (95% CI 1.77–3.28; $\hat \beta $ = 0.89, P < 0.001), and for a 10-fold increase in area the species count is expected to increase 1.81-fold (95% CI 1.36–2.39; $\hat \beta $ = 0.26, P < 0.001). For cleared islands the intercept is 4.31 (95% CI 3.12–5.84; $\hat \beta $ = 1.46, P < 0.001), and no significant relationship was found between species richness and area ( $\hat \beta $ = 0.04, P = 0.295). For invaded islands the intercept is 0.62 (95% CI 0.34–1.05; $\hat \beta $ = −0.48, P = 0.092), and for a 10-fold increase in area the species count is expected to increase 1.77-fold (95% CI 1.14–2.71; $\hat \beta $ = 0.25, P = 0.009). Regardless of size, invaded islands did not have more than four species (Fig. 2).
The largest suite of species included in the analysis was on cleared islands, with the smallest on invaded islands (Fig. 3). Standard residuals from the χ2 test showed a higher than expected proportion of invaded islands with grey-faced petrels (P = 0.001, df = 12, standard residual = 4.66). A higher than expected proportion of cleared islands had sooty shearwaters (P = 0.023, standard residual = 2.39), and a higher than expected proportion of uninvaded islands had fluttering shearwaters Puffinus gavia (P = 0.033, standard residual = 1.52). On invaded islands there were fewer than expected populations of common diving petrels Pelecanoides urinatrix (P = 0.027, standard residual = −1.92). Seven species of seabirds present on cleared and/or uninvaded islands were absent from invaded islands (Fig. 3).
Multi-model inference did not provide strong evidence for a link between the time since eradication and species richness. Distance to Auckland was retained as a predictor in all six models that were within 3 AICc units of the optimal model, with higher species richness on islands further from the city. Distance to mainland and island size were each retained in three of the top six models, but effect sizes were close to zero when both predictors were retained in the same model. Island size had a positive effect on species richness, whereas distance to mainland had a negative effect (Table 3; Supplementary Fig. S4).
*Distance to Auckland was retained as a predictor in all six models that were within 3 AICc units of the optimal model.
On the reference islands, species composition was stable for up to 85 years (Table 1). No species identified in early reports had disappeared from any island. In one group (Ruamaahua) up to four species not recorded previously were present in 2013 (Table 1).
Discussion
Eradicating invasive predators facilitates the recovery of resident seabird populations, and recolonization by species that have been extirpated (Kappes & Jones, Reference Kappes and Jones2014). According to the equilibrium theory of island biogeography, area should be the best predictor of species richness on islands (Preston, Reference Preston1962; MacArthur & Wilson, Reference MacArthur and Wilson1967). In our analysis we assumed that species richness of seabirds on cleared islands should resemble, or begin to resemble, that of uninvaded islands of a comparable area within a similar biogeographical region. Despite our data undoubtedly being influenced by variable survey effort, and temporal and spatial biases such as island size and a lack of pre-invasion census data, seabirds are recolonizing cleared islands in the Hauraki Gulf (Fig. 2). This recolonization is promising, especially as all of the species absent from invaded islands are now breeding on cleared islands (Fig. 3). Furthermore, the uninvaded islands Hongiora and Ruamahuanui appear to have recruited species that had not been observed there previously (Table 1), perhaps as spillover from populations increasing elsewhere as a part of natural assemblage equilibrium processes (C. Gaskin, unpubl. data; Simberloff, Reference Simberloff1974; Buxton et al., Reference Buxton, Taylor, Jones, Lyver, Moller, Cree and Towns2016). Although varying rates of recolonization and local extinction may result in a species richness rebound effect following the removal of predators (Cirtwill & Stouffer, Reference Cirtwill and Stouffer2016), ecological diversity is likely to exert a strong influence on the return of seabird assemblages to equilibrium on cleared islands (Preston, Reference Preston1962; MacArthur & Wilson, Reference MacArthur and Wilson1967; Simberloff, Reference Simberloff1974). We posit that within our case study system the seabird assemblage structure on cleared islands probably reflects the effects of five components: (1) spatial distribution of islands influencing proximity to source populations, foraging areas and/or human activities; (2) time since eradication; (3) behavioural influences, including interspecific interactions; (4) the availability of suitable habitats; and (5) the colonizing ability of individual species (Fig. 4). Our statistics support the influence of space (1) and, to a lesser extent, time (2). However, our model was not able to explain all of the variance. The literature provides some guidance for the behavioural (3), within habitat (4), and demographic (5) influences on seabird recolonization to islands cleared of invasive predators in the Hauraki Gulf.
Space
The rate and pattern of seabird recovery appear to be most influenced by the proximity of cleared islands to source populations, foraging areas and/or human activities. Distance to local populations was the strongest influence on seabird recruitment in the study of Buxton et al. (Reference Buxton, Jones, Moller and Towns2014), where the natural recolonization response fell below 50% when the source population was ≥ 25 km from the recruitment site. Borrelle et al. (Reference Borrelle, Buxton, Jones and Towns2015) found lower than expected species richness on Hauraki Gulf islands beyond the 25 km radius, thus supporting distance to source populations as a key driver of passive seabird recolonization (Buxton et al., Reference Buxton, Jones, Moller and Towns2014). Regarding proximity to foraging grounds, the ideal situation for species is the juxtaposition of resources that reduces costs to the individual (Estades, Reference Estades2001). Seabirds represent the extreme end of these habitat–resource spatial patterns, with foraging grounds often at considerable distances from breeding sites, and distributed discontinuously across space and time (Estades, Reference Estades2001). The disparity between breeding sites and the dynamic nature of foraging areas means that these relationships are not well understood (Fernández-Chacón et al., Reference Fernández-Chacón, Genovart, Pradel, Tavecchia, Bertolero and Piccardo2013). Although we found no compelling evidence of spatial correlation, we note that islands located in the Inner Gulf all had negative residuals (Supplementary Fig. S1). This may be explained by the results from multi-model inference, which suggested that proximity to Auckland negatively influenced the spatial distribution of seabird recovery (Table 3; Supplementary Fig. S4). This could be attributed to a combination of proximity to source populations, optimal foraging theory, and greater habitat modification and marine activity, such as commercial and recreational fisheries, and shipping, in these areas (Borrelle, Reference Borrelle2013; Gaskin & Rayner, Reference Gaskin and Rayner2013). More detailed evaluation of effect modifiers in relation to proximity to human activities and foraging patterns (e.g. using detailed telemetry data) may improve the predictability of the model.
Time
A positive but statistically non-significant increase in species richness was observed with time since eradication. This result may be influenced by a strong temporal artefact because our analysis of seabird recovery on some islands reflects predator eradication carried out c. 20 years ago (Supplementary Fig. S2). Our comparisons of species richness with island area indicate that not all cleared islands have a similar level of seabird richness to uninvaded islands of comparable size. This may indicate that there has been insufficient time for the recolonization of seabirds to those islands, particularly on larger islands that were cleared more recently following the resolution of logistical and operational challenges (Towns & Broome, Reference Towns and Broome2003). Variations in the speed and/or capacity of individual species to recolonize may also be a contributing factor. However, given the generalized linear model predictions of species richness on cleared islands compared to invaded islands it appears that on many cleared islands species richness has increased following eradication of predators (Fig. 2). On some islands these responses have been rapid and unpredicted; for example, Burgess Island has seen the recovery of seven species of Procellariiformes (Table 4; Ismar et al., Reference Ismar, Baird, Gaskin, Taylor, Tennyson and Rayner2014). Furthermore, since the eradications of cats (1970s) and kiore (2004) on Te-Hauturu-o-Toi (Little Barrier Island), the New Zealand storm petrel has reached sufficient numbers for identification of a breeding population of this species, which was hitherto presumed to be extinct (Rayner et al., Reference Rayner, Gaskin, Fitzgerald, Baird, Berg and Boyle2015). Cirtwill & Stouffer (Reference Cirtwill and Stouffer2016) found that species richness on disturbed mangrove islands increased temporarily once the disturbance (predation) was removed. Species assemblages then stabilized, with equilibrium between immigration and extirpation (Cirtwill & Stouffer, Reference Cirtwill and Stouffer2016). A similar process may be operating on some of our study islands, following predator eradication. However, little is known about the temporal dynamics of seabird communities during recolonization of islands. In the case of the volcanic island of Surtsey, in Iceland, nine seabird species established breeding sites successively over the course of 55 years. Only one of these, the Arctic tern Sterna paradisaea, failed to establish a permanent presence and went locally extinct after three breeding attempts across 4 years (Petersen, Reference Petersen2009).
Behaviour
Examples of intrinsic effects on seabird recovery can include situations where intra- and interspecific interactions exert strong influences on recovery rate and recruitment (Danchin et al., Reference Danchin, Boulinier and Massot1998; Parejo et al., Reference Parejo, Danchin and Avilés2005; Buxton et al., Reference Buxton, Jones, Moller and Towns2014). Recruitment to new breeding sites may be limited by the number of immature individuals in local metapopulations; because of the life history traits of seabirds, colony growth can be slow (Parejo et al., Reference Parejo, Danchin and Avilés2005). Also, the recruitment of immature individuals is influenced by the availability of suitable habitat and social cues (Danchin et al., Reference Danchin, Boulinier and Massot1998; Parejo et al., Reference Parejo, Danchin and Avilés2005). These interactions can influence breeding success or immigration rates as a result of limited access to nest sites. Passive recolonization of seabirds on 92 offshore islands in New Zealand following the removal of introduced predators was greater where more than two seabird species were present (Buxton et al., Reference Buxton, Jones, Moller and Towns2014), probably because of a conspecific attraction to preferential habitat (Danchin et al., Reference Danchin, Boulinier and Massot1998; Parejo et al., Reference Parejo, Danchin and Avilés2005; Buxton et al., Reference Buxton, Jones, Moller and Towns2014). In contrast, numbers of grey-faced petrels appeared to decline on the Poor Knights Islands as Buller's shearwaters increased in abundance, even though grey-faced petrels are abundant on Taranga Island in the neighbouring Hen and Chickens group (Harper, Reference Harper1983). Interspecific competition may play a greater role than predator presence in the distribution of grey-faced petrels on the Poor Knight's group, given the higher than expected proportion of invaded islands where this species is present (Fig. 3).
Size & habitat
We found that five of the 16 seabird species included in our analyses were absent from uninvaded islands (Fig. 3), whereas we found the full suite of species across cleared islands, which indicates rapid recolonization by seabirds following a release from predation pressure (Cirtwill & Stouffer, Reference Cirtwill and Stouffer2016), and/or area-related limitations on ecological diversity influencing the availability of suitable habitat (Simberloff, Reference Simberloff1974). Larger invaded islands may also serve as refugia for some species, despite the presence of predators, because of the availability of preferred habitat. For example, Cook's petrels and black petrels were not represented on uninvaded islands; however, both species were present on invaded and cleared islands. Cook's petrel was widely distributed on the mainland of New Zealand before the introduction of mammals (Imber et al., Reference Imber, West and Cooper2003b) but is now restricted to only three breeding sites: Codfish Island (southern New Zealand), Aotea (Great Barrier) and Te-Hauturu-o-Toi (Towns, Reference Towns2009). Nest burrows are found in tall forests, generally > 250 m elevation on the latter two islands, suggesting climate could also play a role in this species’ habitat preferences, with the largest population on predator-free Te-Hauturu-o-Toi (estimated 286,000 breeding pairs in 2007; Rayner et al., Reference Rayner, Clout, Stamp, Imber, Brunton and Hauber2007). Similarly, black petrels were once more widely distributed, breeding at five or more sites on the North Island before 1900 (Imber et al., Reference Imber, McFadden, Bell and Scofield2003a). Black petrels are now restricted to Te-Hauturu-o-Toi and Aotea (estimated 2,000 breeding pairs; E. Bell, Wildlife Management International Ltd, pers. comm., 2015), also in colonies largely restricted to high-altitude (> 400 m) tall forest (Francis & Bell, Reference Francis and Bell2010; Bell et al., Reference Bell, Sim, Scofield and Francis2012). For species such as Cook's and black petrels, habitat suitability and a lack of islands of sufficient size or elevation within their preferential climate envelope may be limiting factors in recruitment to newly predator-free sites. Conversely, more generalist species, such as grey-faced petrels, common diving petrels and fluttering shearwaters, may be less constrained by habitat suitability (Buxton et al., Reference Buxton, Taylor, Jones, Lyver, Moller, Cree and Towns2016).
Recruitment
We found that the proportion of islands with fluttering shearwaters was greatest for uninvaded islands (Fig. 3). Given that fluttering shearwaters are abundant in the region (estimated population of 20,000 pairs; C. Gaskin, unpubl. data), we might expect this species to be more widely distributed on cleared islands, as was found for sooty shearwaters, which have an estimated population of < 800 pairs (Fig. 3; C. Gaskin, unpubl. data). Their contrary distribution suggests that fluttering shearwaters may have a more limited capacity for recruitment to new breeding sites than species such as the sooty shearwater. Behavioural or sex differences between species may also influence recruitment rates (Dittmann et al., Reference Dittmann, Ezard and Becker2007) but there is limited understanding about how these influences may affect recruitment. Additional demographic traits, such as age at first breeding, may influence the species-specific rate of recovery (Nur & Sydeman, Reference Nur and Sydeman1999), where species that have shorter pre-breeding states may recolonize more quickly. On the Mokohinau Islands, for example, populations of common diving petrels and white-faced storm petrels Pelagodroma maoriana increased dramatically in the 2 decades following the removal of R. exulans (Ismar et al., Reference Ismar, Baird, Gaskin, Taylor, Tennyson and Rayner2014). For these two species the age at first reproduction is 2 or 3 years, compared to 5–10 years for black petrels and Pycroft's petrels Pterodroma pycrofti (New Zealand Birds Online, 2015). Mean assemblage age at first reproduction was lowest on invaded islands and highest on uninvaded islands, but these differences were not statistically significant (Kruskal–Wallis rank sum test, χ2 = 5.72, df = 2, P = 0.057; Supplementary Fig. S3a). Furthermore, across cleared islands the mean assemblage age at first reproduction increased with time since eradication, although again this result was not statistically significant (ordinary least squares regression, $\hat \beta $ = 0.013, P = 0.20; Supplementary Fig. S3b). Although we found no statistical evidence of age at first reproduction influencing recolonization potential in our study, the effect of age at first reproduction could have been masked by the small sample size, the relative recency of most eradications, and/or the lack of time-series observations of the recovery process on individual islands. Further research may provide insight into relationships between recovery rate and age at first reproduction.
Conclusion
Predator eradication is an effective conservation tool for protecting and enhancing seabird colonies (Jones, Reference Jones2010; Towns et al., Reference Towns, West and Broome2013; Jones et al., Reference Jones, Holmes, Butchart, Tershy, Kappes and Corkery2016). Such actions can result in fundamental changes to the structure and species richness of island seabird faunas; however, we also found that once assemblages form they undergo little change over many decades despite pervasive extrinsic influences, such as fisheries bycatch and marine pollution. Nonetheless there are numerous constraints on natural recolonization; for example, despite their mobility over water, seabirds rarely establish new breeding sites, because of K-selected traits (Schreiber & Burger, Reference Schreiber and Burger2001; Hamer et al., Reference Hamer, Schreiber, Burger, Schreiber and Burger2002), although some seabird colonies are influenced by emigration and immigration and thus are not closed (Igual et al., Reference Igual, Forero, Gomez and Oro2007; Lawrence et al., Reference Lawrence, Lyver and Gleeson2014). Demographic, intrinsic, extrinsic, habitat and spatial factors also complicate seabird recovery following predator eradications, with additional variation caused by strong species- and site-specific effects (e.g. Buxton et al., Reference Buxton, Jones, Moller and Towns2014). Given the complex influences on recolonization, seabird conservation efforts that involve eradication of invasive species should aim to improve understanding of the species-specific and external effects on demographic traits and population dynamics, as these may elicit responses peculiar to regional faunas. To achieve this, effective and consistent monitoring plans should be an integral component of any seabird island restoration project.
Acknowledgements
We thank all the dedicated bird counters who contributed to this dataset, and those staff that contributed to the successful eradication of predators from the islands in the Hauraki Gulf and New Zealand. SBB was partially supported by the Todd Foundation. PHBS was supported by the U.S. National Science Foundation (Grant PLR-1341649). We also thank H.P. Jones and two anonymous reviewers for their thorough and useful comments.
Author contributions
SBB and DRT conceived and designed the study. DRT and CPG provided the data and SBB and PHBS analysed data. All authors wrote and edited the manuscript.
Biographical sketches
Stephanie Borrelle’s research includes island restoration and recovery processes, and how demography and at-sea threats may be influencing seabird recovery following land-based conservation actions. Philipp Boersch-Supan is a quantitative ecologist focused on marine ecosystems. He is particularly interested in the ecology of oceanic predators and their prey. Chris Gaskin conducts surveys and research on seabirds in northern New Zealand. He coordinated the Important Bird Area programme for New Zealand's seabirds. David Towns’s research has been based on understanding the effects of invasive predators on island ecosystems, investigating methods for restoring modified island ecosystems, and the development of biological and social metrics for measuring progress of restoration projects.