Biodiversity Data Journal :
Research Article
|
Corresponding author: Lucas Lamelas-Lopez (lucaslamelaslopez@gmail.com)
Academic editor: Krizler Tanalgo
Received: 17 Sep 2023 | Accepted: 27 Nov 2023 | Published: 22 Dec 2023
© 2023 Lucas Lamelas-Lopez, Marco Ferrante, Paulo Borges, Isabel Amorim do Rosário, Veronica Neves
This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Citation:
Lamelas-Lopez L, Ferrante M, Borges PAV, Amorim do Rosário I, Neves V (2023) Nest predation of Cory's shearwater Calonectris borealis (Aves, Procellariiformes) by introduced mammals on Terceira Island, Azores. Biodiversity Data Journal 11: e112871. https://doi.org/10.3897/BDJ.11.e112871
|
|
The Azores holds the largest population of Cory's shearwater Calonectris borealis (Cory, 1881) (Aves, Procellariiformes) in the world. Apart from a few mammal-free islets, the bulk of the population breeds in coastal areas on the main human-inhabited islands, where several non-native predators have been introduced. Throughout the entire year of the 2019 breeding season, we used motion-triggered cameras and regularly visited three colonies of Cory's shearwater to identify nest predators and the factors affecting nest predation. A total of 292,624 photos were obtained, of which 97.7% were of Cory’s shearwaters, 1.7% of non-target species (e.g. other birds, rabbits) and 0.52% of potential predators. Of the monitored nests, 25.7% were predated (n = 9), mainly by black rats (n = 8), but also by cats (n = 1). The relative abundance of black rats in the nests was the main factor explaining nest mortality. This variable was significantly and negatively related with the daily survival rate of Cory’s shearwater nestlings. Identification of the main nest predators is crucial for the management and conservation of native bird populations, particularly on oceanic islands, which harbour an important number of threatened and endemic species.
biological invasions, camera traps, invasive species, nest characteristics, oceanic island, seabirds
Biological invasions together with habitat degradation and climatic changes are currently the most important biodiversity erosion drivers and a significant component of human-caused global environmental change (
The Azores Archipelago, located in the North Atlantic Ocean, is one of the most isolated archipelagos worldwide and a high priority area for seabird conservation. Regular breeding species are the Monteiro's storm-petrel (Hydrobates monteiroi (Bolton et al., 2018)), an endemic bird considered as Vulnerable (
Although only one mammal species was present in the Azores before the Portuguese arrival in the 15th century, the Azores noctule Nyctalus azoreum (Thomas, 1901) (
In the present study, we monitored nests of Cory’s shearwater on three colonies on Terceira Island (Azores) using motion-triggered cameras. Over the last years, this technology has become an important tool for understanding the impact of invasive mammals on seabird breeding success (
Our objectives were: (1) to identify Cory’s shearwater nest predators and (2) to determine the abiotic and biotic factors that may directly or indirectly affect nest predation. In particular, we measured variables related with nest characteristics, as the maximum width and height of the cavity entrance, the nest entrance area, the nest depth and the amount of nest entrance coverage by vegetation; and surrounding habitat, as cover of tree and shrubs, herbaceous plants, rocky areas and barren areas within a 25 m radius of the nest, as well as distance to the nearest urban areas, according to another studies about nest predation in the Azores (
The Azores Archipelago is a group of nine volcanic islands and 26 islets, located in the North Atlantic Ocean about 1500 km from Europe and 1900 km from North America (37° to 40° N latitude and 25° to 31° W longitude). The Archipelago is divided into three groups, the Western group, formed by Corvo and Flores Islands, the Central group, formed by Faial, Pico, São Jorge, Graciosa and Terceira Islands and the Eastern group, formed by São Miguel and Santa Maria Islands. The Azorean climate is characterised by stable temperatures, a high annual rainfall and relative air humidity and persistent winds (e.g. http://www.climaat.angra.uac.pt/). The Azorean landscape has suffered severe modifications as a consequence of human settlement, mainly associated with land-use changes by the replacement of native forests by agricultural fields, forestry plantations, namely of Japanese cedar Cryptomeria japonica D.Don and Australian cheesewood Pittosporum undulatum Vent. and urban areas (
The Azores Archipelago harbour regular breeding species such as the Monteiro's storm-petrel, the band-rumped storm-petrel, the Bulwer's petrel, the Cory’s shearwater, the Barolo's shearwater, the Manx shearwater, the yellow-legged gull, the sooty tern, the roseate tern and the common tern. The colonies of these species are generally located on small islets close to the main islands and coastal cliffs dominated by low forests of native and exotic vegetation and rocky/uncovered areas (
The study was performed on three colonies of Cory’s shearwater, on Terceira Island (27º10' W, 38º40' N) (Fig.
Chanoca is located in the southern coast of Terceira Island. This area is formed by cliffs enclosing small rocky bays and it is covered by rocky and barren areas and a low density of herbaceous plants. Sour fig Carpobrotus edulis (L.) N.E.Br. is the dominant plant in the area, being relatively abundant and widespread. Chanoca is located next to an urban area (< 200 m) and it is frequently used by local fishermen to gain access to the sea.
Raminho is located in the north-western part of Terceira Island. The area is formed by cliffs, dominated by a low forest of native vegetation, namely Erica azorica and small and scarce individuals of Morella faya (Aiton) Wilbur. The nests are located in the rocky cliffs and within the forest. Although the colony is relatively inaccessible, there is a narrow dirt road that goes through the forest to the lowest part of the colony, near to the coast. Raminho is located approximately 700-800 m from the nearest urban area.
Agualva is located in the northern part of Terceira Island. The area consists of a rocky area scarcely covered by native vegetation, mainly Erica azorica Hochst. ex Seub. Although it is relatively isolated (> 1400 m from the nearest urban area), it is crossed by a walking trail (PR2TER - Baías da Agualva).
Fieldwork was performed between 4 April and 25 October 2019, covering the entire breeding season of Cory’s shearwater in the Azores (
From camera-trap records, we obtained the Relative Abundance Index (RAI) (
For each nest, we measured the maximum width and height of the cavity entrance, the nest entrance area (calculated by multiplying nest width and height) and the nest depth (distance from the entrance to the back of the cavity) (see
In addition, we calculated the cumulative proportion of covered (trees, shrubs and herbs cover) and uncovered areas (cover of rocky and barren areas) within a radius of 25 m from the nest. Nest UTM coordinates were recorded using a GPS navigator (Garmin eTrex). We also measured the distance from the nest to the nearest urban area using the ArcGIS software (
We calculated the daily mortality rate (DME) as (P/D), where P stands for “Predation” and is the status of the nest at the end of the breeding season (i.e. predated or not predated) and D stands for “Days” and is the number of days when the camera was recording. Consequently, the daily mortality rate is a continuous value between 0 (i.e. the nest survived) and 1 (i.e. the nest was predated after the first day of recording). The DME was used as a response variable in two zero-inflated beta-regressions; one including only abiotic factors and the other only biotic factors. The beta-regression distribution is suitable for continuous proportional data with a high frequency of zeroes, as in our dataset. Log-transformation was applied to all variables with outliers (i.e. nest area; nest depth; distance from the nearest urban centre and the RAIs of Rattus sp.). Before creating the models, we calculated the Pearson Correlation Coefficient (r) to test for potential correlation between the abiotic factors and, independently, between the biotic factors.
Since we found high correlation between nest area and nest depth (r = 0.45, Suppl. material
No concerning correlations were detected between biotic factors (Suppl. material
Moreover, we tested which factors affected the activity of rats around Cory’s shearwater nests formulating a linear model were the log-transformed RAI of Rattus sp. was the response and the presence of vegetation at the entrance of the nest (binary), the percentage of vegetation near the nest, the log-transformed nest area, the log-transformed nest depth, the log-transformed distance from the nearest urban centre and the RAIs of potential predators of rats (cat and least weasel) were the fixed factors and the site was the random factor. This model was simplified using stepwise selection with the R function step. The final model included the presence of vegetation at the entrance of the nest, the percentage of vegetation near the nest and the log-transformed distance from the nearest urban centre and did not show any concerning pattern in the residuals (Suppl. material
An overview of the model can be found in Table
Overview of the statistical models with the daily mortality estimates (DME) and the log-transformed RAI of Rattus sp. (logRAIrr) as responses. Fixed factors are coded as follows: log-transformed nest area (logArea), log-transformed distance from the nearest urban centre (logDistUrb), log-transformed RAI of Rattus sp. (logRAIrr), presence of vegetation at the entrance of the nest (VegEntr), the percentage of vegetation near the nest (Cover).
Model |
Response |
Distribution |
Fixed effects |
Random effects |
Abiotic factor |
DME |
Zero-inflated beta distribution |
logArea + logDistUrb |
Site |
Biotic factor |
DME |
Zero-inflated beta distribution |
logRAIrr |
Site |
RAI Rattus sp. |
logRAIrr |
Gaussian |
VegEntr + Cover + logDistUrb |
Site |
A total of 292,624 photos were obtained, of which 286,061 (97.7%) portrayed Cory’s shearwaters, 5045 (1.7%) non-target species (e.g. other birds, rabbits) and 1518 (0.52%) potential predators (Table 1). Overall, nine of the thirty-five successful monitored nests (25.7%) across the three monitored colonies were predated throughout the breeding season. Of these, two were predated by black rats (22.2%) and one (11.1%) by feral cats. Predators were not identified in six nests (66.7%), but the high activity of rats in these nests suggests that they were the most likely predators (Table
Detected potential predator species, prey age-class, total number of photos recorded, number of predated nests and mean Relative Abundance Index (RAI) of potential predators. *Black rats were only recorded predating two nests, but the high RAI detected nearby the other predated nests suggest that they were likely responsible for the predation of eight nests; N/A – Non applicable.
Predator | Prey age-class | Number of records | Number of predated nests | Mean RAI |
Domestic cat - Felis catus | Chick | 176 | 1 | 7.26 |
Black rat - Rattus rattus | Egg | 973 | 8* | 42.62 |
House mouse - Mus musculus | N/A | 363 | 0 | 10.10 |
Least weasel - Mustela nivalis | N/A | 6 | 0 | 0.20 |
The highest nest predation was recorded in Agualva (33.3% of the nests), followed by Raminho (30.8%) and Chanoca (15.4%). The mean RAI of Cory’s shearwater in nests that were predated (12,798.2) was higher than in nests that were not predated (9806.0). The same was true for the RAI of recorded predators, such as black rats (81.8 vs. 29.1 for predated and not predated nests, respectively) and cats (8.4 vs. 6.9), but not for potential predators of eggs or chicks such as the least weasel (0.0 vs. 0.3) and the house mouse (8.1 vs. 10.8).
No abiotic factor significantly affected the DME, while the log-transformed RAI of Rattus sp. was significantly and negatively related to the nest survival (p < 0.001, GLMM; Fig.
Outputs of the three statistical models with daily mortality estimates (DME) or the log-transformed RAI of Rattus sp. (logRAIrr) as response variables.
Response |
Conditional model |
Variance |
SD |
||
DME |
Site |
6.42E-11 |
8.01E-06 |
||
Estimate |
SE |
z-value |
Significance |
||
Intercept |
-3.56791 |
1.80982 |
-1.971 |
0.0487 * |
|
logDistUrb |
-0.18838 |
0.16936 |
-1.112 |
0.266 |
|
logArea |
0.08966 |
0.17658 |
0.508 |
0.6116 |
|
Zero-inflation model |
|||||
Intercept |
1.0609 |
0.3867 |
2.743 |
0.00609 ** |
|
DME |
Site |
6.70E-01 |
8.19E-01 |
||
Estimate |
SE |
z-value |
Significance |
||
Intercept |
-6.2332 |
0.55015 |
-11.33 |
< 2e-16 *** |
|
logRAIrr |
0.60395 |
0.07612 |
7.934 |
2.12e-15 *** |
|
Zero-inflation model |
|||||
Intercept |
1.0609 |
0.3867 |
2.743 |
0.00609 ** |
|
logRAIrr |
Site |
4.29E-09 |
6.55E-05 |
||
Estimate |
SE |
z-value |
Significance |
||
Intercept |
-0.35059 |
1.2469 |
-0.281 |
0.779 |
|
Cover |
0.01662 |
0.01149 |
1.446 |
0.148 |
|
VegEntr |
-1.15858 |
0.49786 |
-2.327 |
0.020 * |
|
logDistUrb |
0.39758 |
0.21946 |
1.812 |
0.070 . |
The best model explaining the RAI of Rattus sp. included the percentage of vegetation nearby the nest (despite this factor being non-significant), the log-transformed distance from the nests to the urban areas and the presence of vegetation at the entrance of the nest. The RAI of Rattus sp. was marginally significantly (p = 0.07, LM; Fig.
The impacts of rodents, particularly rats and cats on insular biota have been widely documented (
The abundance of introduced mammal predators (
Other potential nest predators of Cory’s shearwater eggs and chicks include house mice and mustelids. In our study, we detected both species in the nests, but predation was never recorded. The impact of the house mouse as bird nest predator has rarely been reported for the Azores (
The identification of the main nest predators is crucial for the management and conservation of native bird populations (
We gratefully acknowledge the Natural Park of Terceira Island for logistic support and to Jose Ortolá, Natalia Fierro, Ana Sánchez, Giulia Spadoni and Clara Polaino for field assistance We also thank to Rui Medeiros for help with Figure 1.
LLL was supported by a grant from the Fundação para a Ciência e Tecnologia - FCT (SFRH/BD/115022/2016); VN was supported by a contract Investigadores MarAz – PO Azores 2020; IRA was funded by Portuguese national funds through FCT – Fundação para a Ciência e a Tecnologia, I.P., under the Norma Transitória – DL57/2016/CP1375/CT0003. Open access was funded by the project FCT-UIDB/00329/2020-2024 (Thematic Line 1 – integrated ecological assessment of environmental change on biodiversity).
Our methods met all ethical guidelines stipulated by the standards and policies of Portuguese law and special care was taken in order to minimise disturbance of biota. Fieldwork was conducted under permit from the Azores Government, Regional Environmental and Climate Changes Directorate.
LL-L, MF and VN contributed to the conception/design of the study. LL-L performed the fieldwork. LL-L and MF performed the data analysis. PAVB and IRA contributed with logistical support. All authors contributed to manuscript writing.
Pearson Correlation Coefficient between the abiotic factors: the log-transformed distance from the nearest urban centre (logDistUrb), the log-transformed nest area (logArea), the log-transformed nest depth (logDepth), the percentage of vegetation near the nest (Covered) and the presence of vegetation at the nest entrance as a binary variable (Veg_Entr2).
Model validation plot created using the simulateResiduals function in DHARMa of the daily mortality rate model including abiotic factors.
Pearson Correlation Coefficient between the biotic factors: the Relative Abundance Index (RAI) of the Cory’s shearwater (RAI_Cd), the log-transformed RAI of Rattus sp. (logRAIrr), the log-transformed RAI of the least weasel (logRAImn), the RAI of the domestic cat (RAI_Fc) and the RAI of the house mouse (RAI_Mm).
Model validation plot created using the simulateResiduals function in DHARMa of the daily mortality rate model including biotic factors (i.e. the log-transformed RAI of Rattus sp.).
Model validation plot created using the simulateResiduals function in DHARMa of the log-transformed RAI of Rattus sp. model including the presence of vegetation at the entrance of the nest, the percentage of vegetation near the nest and the log-transformed distance from the nearest urban centre.