| Literature DB >> 32489622 |
Pavel Pipek1,2, Tim M Blackburn3,4, Steven Delean5, Phillip Cassey5, Çağan H Şekercioğlu6,7,8, Petr Pyšek1,2.
Abstract
Invasive alien species are a major threat to biodiversity and human activities, providing a strong incentive to understand the processes by which alien invasion occurs. While it is important to understand the determinants of success at each of several invasion stages-transport, introduction, establishment, and spread-few studies have explored the first of these stages. Here, we quantify and analyze variation in the success of individual animals in surviving the transport stage, based on shipping records of European passerines destined for New Zealand. We mined the original documents of Acclimatisation Societies, established in New Zealand for the purpose of introducing supposedly beneficial alien species, in combination with recently digitized newspaper archives, to produce a unique dataset of 122 ships that carried passerines from Europe to New Zealand between 1850 and 1885. For 37 of these shipments, data on the survival of individual species were available. Using generalized linear mixed models, we explored how survival was related to characteristics of the shipments and the species. We show that species differed greatly in their survival, but none of the tested traits accounted for these differences. Yet, survival increased over time, which mirrors the switch from early haphazard shipments to larger organized shipments. Our results imply that it was the quality of care received by the birds that most affected success at this stage of the invasion process.Entities:
Keywords: 19th century; New Zealand; alien birds; invasions; shipping; survival; transport
Year: 2020 PMID: 32489622 PMCID: PMC7244811 DOI: 10.1002/ece3.6143
Source DB: PubMed Journal: Ecol Evol ISSN: 2045-7758 Impact factor: 2.912
Figure 1Between 1850 and 1884, passerines were shipped at least 122 times from Europe, mostly from port cities in England, to New Zealand. One ship was wrecked not far from the start. The thickness of the lines corresponds to the number of ships that travelled between the ports which are connected by them, and the colors serve only to distinguish the lines from different ports from each other. The brown and the purple correspond to England and Germany, respectively, and the port is unknown. The destinations are abbreviated as follows (NZ regions are in brackets): AUK—Auckland (Auckland), BLU—Bluff (Southland), GLA—Glasgow, LIV—Liverpool, LON—London, LYT—Lyttelton (Canterbury), NAP—Napier (Hawke's Bay), NEL—Nelson (Nelson), NPL—New Plymouth (Taranaki), PCH—Port Chalmers (Otago), PLY—Plymouth, SOU—Southampton, WEL—Wellington (Wellington)
Summary statistics for different NZ ports
| Port | Region | Number of ships | Number of species | Median survival |
|---|---|---|---|---|
| Auckland | Auckland | 42 | 23 | 0.18 (18) |
| Bluff | Southland | 3 | 7 | |
| Lyttelton | Canterbury | 32 | 29 | 0.29 (7) |
| Napier | Hawke's Bay | 4 | 15 | |
| Nelson | Nelson | 6 | 12 | 0.33 (5) |
| Port Chalmers | Otago | 18 | 28 | 0.60 (3) |
| New Plymouth | Taranaki | 2 | 8 | |
| Wellington | Wellington | 14 | 16 | 0.00 (3) |
Median survival was counted only for a subset of shipments (n = 37) and for taxa as specified in methods (number of shipments shown in brackets).
Analysis of deviance table for nested model comparisons using likelihood ratio tests
| Model comparison | Number of parameters | Log‐likelihood | Deviance |
|
|
|---|---|---|---|---|---|
| Duration + year + season | 9 | −448.64 | 4.94 | 6 | .55 |
| ns(Duration, | 15 | −446.17 | |||
| Duration | 5 | −452.29 | 7.28 | 4 | .12 |
| Duration + year + season | 9 | −448.64 | |||
| Year | 5 | −450.00 | 2.72 | 4 | .61 |
| Duration + year + season | 9 | −448.64 | |||
| Season | 7 | −450.09 | 2.89 | 2 | .24 |
| Duration + year + season | 9 | −448.64 | |||
| Intercept only | 4 | −452.43 | 0.29 | 1 | .59 |
| Duration | 5 | −452.29 | |||
| Intercept only | 4 | −452.43 | 4.86 | 1 | .03 |
| Year | 5 | −450.00 | |||
| Intercept only | 4 | −452.43 | 4.69 | 3 | .20 |
| Season | 7 | −450.09 |
All models, that is, generalized linear mixed models (with beta‐binomial logit function) with proportion of surviving bird individuals out of the total loaded as a response variable, contain an intercept term; “ns” indicates natural spline. Difference in model deviance is assumed to be distributed as chi‐squared with specified degrees of freedom (df).
Figure 2Conditional estimates of probability of survival (log odds scale) for ship voyages to five different regions in New Zealand; error bars are standard errors. Positive and negative values with standard error not overlapping the zero line indicate significantly higher survival or mortality, respectively. Zero value indicates no deviation from the average survival conditional on the other modeled effects. Shipments highlighted in blue were organized and accompanied by Richard Bills or his son Henry
Figure 3Conditional estimates of probability of survival (log odds scale) for bird species; error bars are standard errors. Positive and negative values with standard error not overlapping the zero line indicate significantly higher survival or mortality, respectively. Zero value indicates no deviation from the average survival conditional on the other modeled effects
Analysis of deviance tables for nested model comparisons for species trait variables, one continuous—body mass, three binaries—sexual dimorphism, sedentary, and short‐distance migrant (the latter together describe migratory behavior—two zeroes equal to long‐distance migration), and one categorial—the diet (possible values insectivorous, omnivorous, plant/seed eater)
| Model comparison | Number of parameters | Log‐likelihood | Deviance |
|
|
|---|---|---|---|---|---|
| Year | 5 | −450.00 | 1.04 | 1 | .31 |
| Year + body mass | 6 | −449.49 | |||
| Year | 5 | −450.00 | 1.84 | 1 | .17 |
| Year + sexual dimorphism | 6 | −449.08 | |||
| Year | 5 | −450.00 | 0.20 | 1 | .65 |
| Year + sedentary | 6 | −449.90 | |||
| Year | 5 | −450.00 | 2.96 | 1 | .09 |
| Year + short‐distance migrant | 6 | −448.53 | |||
| Year | 5 | −450.00 | 1.72 | 2 | .42 |
| Year + diet | 7 | −449.15 |
All models, that is, generalized linear mixed models (with beta‐binomial logit function) with proportion of surviving bird individuals out of the total loaded as a response variable, contain an intercept term. Difference in model deviance is assumed to be distributed as chi‐squared with specified degrees of freedom (df).