Subscriber access provided by RYERSON UNIV
Article
Spatial Risk Assessment of Alien Invasive Plants in China Fan Bai, Ryan Chisholm, Weiguo Sang, and Ming Dong Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/es400382c • Publication Date (Web): 05 Jun 2013 Downloaded from http://pubs.acs.org on June 14, 2013
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Environmental Science & Technology is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 32
Environmental Science & Technology
1
Spatial Risk Assessment of Alien Invasive Plants in China
2
Fan Bai1, Ryan Chisholm2,3, Weiguo Sang1*, Ming Dong1*
3 4
1
5
Chinese Academy of Science, Beijing 100093, China
6
2
7
14 Science Drive 4, Singapore 117543
8
3
9
of Panamá
State Key Laboratory of Vegetation and Environmental Change, Institute of Botany, Department of Biological Sciences, Faculty of Science, National University of Singapore, Smithsonian Tropical Research Institute, P.O. Box 0843-03092, Balboa, Ancón, Republic
10 11
*Corresponding author:
12
Weiguo Sang, State Key Laboratory of Vegetation and Environmental Change, Institute of
13
Botany, Chinese Academy of Science, Beijing 100093, China, E-mail:
[email protected];
14
Phone: +86 10 62836278; FAX: +86 10 82599519.
15
Ming Dong, State Key Laboratory of Vegetation and Environmental Change, Institute of
16
Botany,
17
[email protected]; Phone: +86 10 82594676 .
Chinese
Academy
of
Science,
Beijing
100093,
18 19
Table: 2
20
Figure:6 (figure 6 is large)
21
Word Number: 6580
22
1
ACS Paragon Plus Environment
China,
E-mail:
Environmental Science & Technology
23
Abstract:
24
The large-scale distribution patterns of alien invasive plants (AIP) can provide key
25
information and a theoretical basis for management strategies, including the prevention of
26
invasions, the control and eradication of established AIPs, and the identification of areas at
27
high risk of invasion. This study aims to quantify distribution patterns of AIP in China, to
28
develop approaches that measure the social, economic and ecological impacts, and to
29
identify areas that are at higher risk of plant invasion. Based on published literature, there
30
were 384 AIPs in China, representing 233 genera from 66 families. Climatic factors were
31
among the primary factors determining AIPs’ overall distribution patterns. The majority of
32
AIPs were tropically distributed in China, meaning that they were mainly restricted to
33
southern China. Temperate-distributed AIPs, those distributed only or predominantly in
34
northern China, were fewer but had higher average rates of spread than tropically
35
distributed AIPs. Average ecological and economic impact per AIP was negatively
36
correlated with AIP richness, meaning that areas with relatively few AIPs nevertheless have
37
some of the most detrimental ones. Our comparative evaluation showed that the risk of
38
invasion differed among regions of China, with high-risk areas in southern China (Yunnan,
39
Guangxi and Guangdong) and central coastal areas of eastern China (Shandong, Hebei,
40
Jiangsu). In the context of climate change, areas around latitudes of 33°N, including Hebei,
41
Shandong, Henan and Jiangsu, should be given more attention for the control and
42
prevention of plant invasions. Predictions of high-risk areas for future invasions differed
43
depending on the scale of aggregation and the evaluation index, indicating that invasive
44
risk assessments should be based on multiple factors.
45
2
ACS Paragon Plus Environment
Page 2 of 32
Page 3 of 32
Environmental Science & Technology
46 47
INTRODUCTION
48
Alien Invasive Plants (AIPs) are defined as non-native plant species that have established
49
naturalized populations and subsequently spread over large geographic areas and have
50
negative economic and ecological impacts1, 2. Global environmental change and economic
51
development have accelerated the colonisation and spread of AIPs, many of which pose
52
serious threats to crop production and native biodiversity3, 4. In China, AIPs have caused
53
ecological disasters and economic losses5, 6. Many studies across the world have focused on
54
the mechanisms that facilitate invasion and on the prediction of future invasions for the
55
purposes of risk assessment7-9. The distribution patterns of AIPs, based on large-scale data
56
analysis, can provide key information and a theoretical basis for their management and
57
control10-12. In China, large-scale and comprehensive analyses of AIP distributions have
58
seldom been possible, because of a shortage of information13.
59 60
AIP distribution patterns depend on species’ traits and reflect the invasiveness of individual
61
plant species14, 15. Many factors have been associated with ecosystem invasibility and the
62
invasiveness of individual plant species. These include ecological factors (e.g., dispersal
63
opportunity, initial population sizes, residence times, numbers of introduction attempts,
64
biodiversity, geography and climate), economic activities (e.g., land use, agricultural intensity,
65
the presence of forestry and fishery industries, trade intensity and road density), and social
66
factors (e.g., demography and history)16, 17.
67 68
Despite decades of research, universal determinants of invasiveness and invasibility have yet 3
ACS Paragon Plus Environment
Environmental Science & Technology
69
to be discovered17. Therefore, the risk assessment of invasive species and invasibility of
70
invaded ecosystems must be considered jointly in multi-dimensional analyses that include
71
traits of AIPs, distribution range, richness, abundance, invasion rate, the degree of invasive
72
impact, and so on9, 18-21. However, most large-scale analyses of AIP distributions in China to
73
date have focussed solely on AIP richness as a metric of ecosystem invasion risk 5, 22.
74 75
Knowing which invaded areas are most at risk and which traits of invasive species
76
determine invasiveness is potentially of great value to invasion ecology and management23.
77
The pattern of AIP distribution in China is related to both the natural and economic
78
development conditions24. Southern China has been considered the area with highest
79
invasion risk in China because it has the highest AIP richness, has suffered large economic
80
losses caused by AIP, has suitable climatic conditions for many invasive species and is
81
undergoing an economic boom5,
82
invaded by several high-impact AIPs that are continuing to expand, it has been the subject
83
of relatively fewer studies26. Accordingly, investments in AIP prevention have been
84
directed mainly to southern China rather than northern China. Furthermore, these
85
investments have been directed towards AIPs with the largest economic impact, neglecting
86
those that may have potentially large ecological impacts.
10, 25
. Meanwhile, although Northern China has been
87 88
In view of the global acceleration of plant invasions and the limited knowledge about AIP
89
in China, the identification of AIP distribution hotspots in China is a key priority. By
90
collecting information on AIPs in China from published floras, internet databases and other 4
ACS Paragon Plus Environment
Page 4 of 32
Page 5 of 32
Environmental Science & Technology
91
literature, this paper builds a set of invasive indices to describe the AIP distribution patterns,
92
to shed light on the determining factors and to identify high-risk areas for invasion in China.
93
We hypothesised that (1) Southern China would have more AIPs and higher invasion risk
94
than Northern China, consistent with previous studies5,10; (2) Climatic factors would be the
95
primary determinants of AIP distributions in China, as they are across the world27and as
96
was found in one previous study in China24; (3) The invasion risk assessment would be
97
sensitive to the indices of invasion used, because different metrics convey different
98
information about the invasion process9,21.
99
METHODS
100
Study area
101
The study area covers the whole of China (18.05–54.55°N, 74.55–135.05°E), and employs
102
four scales of regional subdivision: 28 provinces, 6 geo-districts, and North/South (Figure
103
1). The species distribution data by provinces were the basic units used to estimate plant
104
invasions at the province scale. The four centrally administered municipalities (Beijing,
105
Tianjin, Shanghai and Chongqing), Hong Kong and Macao were merged into neighbouring
106
provinces (Hebei, Jiangsu, Sichuan and Guangdong), respectively. Based on similarities in
107
environmental and economic conditions, six aggregated geo-districts were defined as
108
Northeast (Liaoning, Jilin, and Heilongjiang), North center (Hebei, Shanxi, and
109
Neimenggu), Northwest (Shaanxi, Gansu, Qinghai, Ningxia, and Xinjiang), Southeast
110
(Jiangsu, Zhejiang, Anhui, Fujian, Taiwan, Jiangxi, and Shandong), South center (Henan,
111
Hubei, Hunan, Guangdong, Guangxi, and Hainan), and Southwest (Sichuan, Guizhou, 5
ACS Paragon Plus Environment
Environmental Science & Technology
112
Yunnan, and Xizang). A further scale of aggregation was attained by combining Northeast,
113
North center and Northwest into a North region (11 provinces) and Southeast, South center
114
and Southwest into a South region (17 provinces).
115
Classification of distributions of AIPs
116
Based on recorded presences and absences of invasive plants in provincial administrative
117
units, species were classified into five distribution types (SI 1): Cosmopolitan species that
118
occur in all 28 provinces (CO); Restricted southern species, that occur only in the southern
119
provinces (OS); Species that occur in a higher proportion of southern than northern
120
provinces (PS); Species that occur in a higher proportion of northern than southern
121
provinces (PN); Restricted northern species, that occur only in the northern provinces (ON).
122
The OS and PS were further classified as tropically distributed species; ON and PN were
123
classified as temperate-distributed species.
124
Data sources
125
A list of alien invasive higher plants in China (SI 1) was compiled from five sources (SI 2).
126
The distribution record of invasive plant presence and absence in provincial administrative
127
units and each AIP’s life form, origin, and first recorded year of presence were taken from
128
Flora of China28, Chinese Virtual Herbarium (CVH) (http://www.cvh.org.cn), Catalogue of
129
Life: Higher Plants in China (CNPC) (http://www.cnpc.ac.cn) and other literature published
130
before April 20115, 10, 29-34. The number of native higher plant species per province was also
131
obtained from CNPC as an indicator of total provincial species richness.
132 6
ACS Paragon Plus Environment
Page 6 of 32
Page 7 of 32
Environmental Science & Technology
133
Environmental variables were collected from the official Chinese government website (http:
134
//www.cma.gov.cn): mean annual temperature (MAT) (℃), mean annual precipitation
135
(MAP) (mm), mean temperature of the warmest month (MTWM) (℃), growing degree
136
days on a 0º basis (GDD0) ( ℃ ), photosynthetically active radiation (PAR), mean
137
precipitation in June, July and August (PJJA) (mm), climatic moisture index (MI), aridity
138
index (α) and forest cover (FC). Indices of economic and social activity, were obtained
139
from the official website (http://www.stats.gov.cn) of the National Bureau of Statistics of
140
China: gross domestic product (GDP), population density, total value of imports (log
141
transformed), inbound tourist arrivals (log transformed), total length of railway mileage
142
(RM), total length of inland waterway mileage, total length of substandard highway
143
mileage, and cargo moved per year.
144
Invasive indices
145
Four indices of spatial distribution patterns were used. Both invasive species richness and
146
relative invasive species richness (i.e., invasive species richness as a proportion of total
147
species richness) are metrics of levels of biological invasion (invasion level)9. Absolute
148
invasive species richness represents the number of invasive species and invasive species
149
richness relative to total species richness represents the proportion of invasive species.
150 151
To assess the economic and ecological impacts of individual AIPs, we constructed two
152
indices: an “Impact” index that encompasses social, economic and ecological impacts of an
153
invasive species, and a “Spread” index that reflects the rate of spread of a species across
7
ACS Paragon Plus Environment
Environmental Science & Technology
154
China. The impact index for a species was calculated by summing the number of categories
155
of source material in which it was mentioned (1 to 4; see SI 3 for details). The spread index
156
was measured as the rate of expansion of a species across geo-districts over time
157
(converted to a score from 1 to 4; see SI 3 for details).
158 159
All the provincial invasive indices were mapped using ArcGIS9.3 (ESRI, Redlands, CA,
160
USA). Each index was discretised onto a rank scale from 1 to 4, such that each rank
161
represented an equal interval of the index.
162
Data analysis
163
We simplified the classification of growth forms of AIPs into herbaceous and woody and
164
compared the prevalence of these two different growth forms across distribution types (i.e.,
165
CO, OS, PS, ON, PN) using a chi-squared test (R version 2.13.2). We compared the
166
prevalence of species with different impact and spread indices across geo-districts (i.e., NE,
167
NC, NW, SE, SC, SW), again using chi-squared tests.
168
From the full list of environmental and economic explanatory variables under “Data sources”
169
we selected two environmental variables (MAT and MAP) and two economic variables (GDP
170
and RM) that explained the largest proportion of the variance (using a Principal Components
171
Analysis (PCA); see SI 4 for details). These four explanatory variables were used in all
172
subsequent analyses. We investigated the relationship between our two invasion indices and
173
the environmental and economic variables across provinces with stepwise linear regression
174
analysis (SPSS 13.0). We also ran a Canonical Correspondence Analysis (CCA; Canoco for
175
Windows 4.5)35 to identify relationships between AIP distribution types and the
176
environmental and economic variables. For the CCA, we calculated mean values of the four 8
ACS Paragon Plus Environment
Page 8 of 32
Page 9 of 32
Environmental Science & Technology
177
explanatory variables for each species by taking averages across the provinces in which the
178
species occurred.
179 180
RESULTS
181
Patterns of AIPs
182
There were 384 higher plant species recorded as alien invasive in China. These represented
183
233 genera and 66 families (SI 1). The invasive plants mainly consisted of members of the
184
families Asteraceae, Poaceae, Fabaceae and Brassicaceae, and genera Amaranthus in
185
Amaranthaceae, Euphorbia in Euphorbiaceae, and Solanum in Solanaceae (SI 1). Both herb
186
(327 species) and woody (52 species) AIPs were more likely to have southern than northern
187
distributions, but herbaceous species were more likely to be cosmopolitan species
188
(chi-square test, df = 4, χ2 = 20.50, p < 0.001) (Figure 2).
189 190
The percentage composition of AIP impact indices was very similar across the six
191
geo-districts (chi-square test, χ2 =3.1605, df = 15, p = 0.9994) (Figure 3A). In every district,
192
fewer than ten species had the highest impact index (IM=4). In contrast, the spread indices
193
were significantly skewed towards higher values in northern China (chi-square test, χ2
194
=53.8344, df = 15, p < 0.001) (Figure 3B), indicating that AIP in the north (especially the
195
NE and NW geo-districts) spread faster and wider than those in the south.
196
Environmental parameters and patterns of AIPs
197
At a national level, invasive species richness was significantly positively related to total 9
ACS Paragon Plus Environment
Environmental Science & Technology
198
(invasive plus native) plant richness across provinces (R2=0.542, p=0.004) (Figure 4). Similar
199
qualitative patterns were observed within most regions.
200 201
The CCA of all the AIPs by environmental and economic factors (Figure 5) and stepwise
202
linear regression analyses (Table 1) both suggested that climatic factors are primary
203
determinants of the overall distribution patterns of AIPs. In the CCA ordination diagram
204
(Figure 5), the first axis (Axis 1) explained 60.3% of the variation in AIP distribution types.
205
This axis chiefly represented a climatic gradient (factor loadings on MAT and MAP were
206
-0.87 and -0.92). The second axis explained a further 17.2% of the plants’ distribution types,
207
and mainly represented an economic development gradient (factor loadings on GDP and
208
railway mileage were -0.85 and -0.40). The average score value (mean) and variation
209
(standard deviation, σ) of temperate-distributed species (mean=0.5422, σ=0.6332) on Axis 2,
210
especially PN (mean=0.9840, σ=0.8585) were both greater than those of tropically distributed
211
species (mean=0.2904, σ=0.3409).
212 213
Stepwise linear regression analyses (Table 1) showed that the various invasion indices for
214
each province were better predicted by climatic variables (MAT and MAP) than by economic
215
variables, although GDP was a significant predictor of absolute and relative AIP richness, and
216
railway mileage was a significant predictor of both the total and average impact scores.
10
ACS Paragon Plus Environment
Page 10 of 32
Page 11 of 32
Environmental Science & Technology
217
218
The invasive indices distribution pattern in China
219
AIP richness was highest in southern China (Yunnan, Guangxi and Guangdong)
220
6A), while relative richness tended to be higher in eastern China, especially in Hebei,
221
Shandong and Jiangsu (Figure 6B). The spatial patterns of total impact and total spread
222
closely tracked the spatial pattern of AIP richness (Figure 6C-D).
223
varied little across the country (Figure 6E). Average spread index was higher in the North
224
region than in the South region, and there was a hotspot on the eastern boundary between
225
North and South, including the provinces Liaoning, Shandong, Henan, Shanxi, Shaanxi and
226
Ningxia (Figure 6F).
(Figure
Average impact index
227 228
DISCUSSION
229
Distribution of alien invasive plants in China
230
Overall, the highest richness of invasive species in China was found in the south and in the
231
Hebei Province area, while the highest relative richness was found in central coastal areas of
232
eastern China (Figure 6A and B). Invasive species richness was positively correlated with
233
local species richness (Figure 4). The majority of AIP in China have tropical distributions
234
(OS and PS) and originate from the Americas (Figure 2, SI 1). Other regional characteristics
235
that predicted invasive species richness were mean annual temperature and GDP (Figure 5,
236
Table 1). These overall distribution patterns and ecosystem invasibility may be explained by a
237
combination of biogeographic, environmental and economic factors36.
238 11
ACS Paragon Plus Environment
Environmental Science & Technology
239
Globally, alien species richness tends to be positively related to native species richness10, 33, 37,
240
38
241
are both larger at lower latitudes) and environmental factors41-44 (the conditions that allow
242
large numbers of native species to persist also allow large numbers of alien species to persist).
243
In an exception to this general trend, alien species richness is usually seen to drop sharply in
244
tropical regions of the world, and this has been attributed to the complexity of biotic
245
interactions in high-diversity systems acting as a barrier to invasion. Our plant data from
246
China are consistent with the general global patterns, insofar as there is an overall positive
247
relationship between AIP richness and native richness (Figure 4) but with low relative AIP
248
richness in the tropics (Figure 6B). However, our result that absolute AIP richness is actually
249
highest in the tropics in China stands in contrast to data from most of the rest of the world37.
250
Further studies are needed to investigate this pattern, which, as an exception to the general
251
rule of lower alien richness in the tropics, may shed light on the mechanisms underlying
252
global biogeographical patterns of alien invasions.
and this is usually attributed to biogeographical factors39, 40 (native and alien species pools
253 254
Economic factors associated with human activity can also drive distribution patterns of alien
255
invasive plants6, 45. While AIP distributions between regions in China were differentiated
256
mainly by climatic variables, distributions within regions were differentiated more by
257
economic variables. In addition, economic variables appeared to have a stronger influence in
258
temperate than in tropical regions, suggesting that human activities have had a greater impact
259
in temperate regions (Figure 5). A caveat here is that there is a possible observation bias
260
resulting from greater research intensity, and hence higher detection rates of alien plant 12
ACS Paragon Plus Environment
Page 12 of 32
Page 13 of 32
Environmental Science & Technology
261
species, in economically developed areas (e.g., in the vicinity of Hebei province). Although
262
we cannot rule the possibility of an observation bias, our results are consistent with
263
global-scale data showing that invasive species richness is strongly positively associated with
264
international trade intensity46. Our observations from China reinforce the notion that invasive
265
species should receive particular attention in regions with booming economies33.
266 267
Climatic factors were found to be the primary factors impacting invasive plants’ distribution
268
patterns. In coming decades, global climate change will interact with biological invasions by
269
allowing new species to invade and by modifying the ecological and economic impacts of
270
existing invasive species47. Recent climate-vegetation models at a national scale in China
271
have indicated that forest ecosystems are moving northwards while alpine ecosystems are
272
shrinking48. These trends are consistent with observations in high northern latitudes across
273
the world49. In China, the responses to climate change of ecosystems at latitudes higher than
274
33°N are particular uncertain50. Because these latitudes coincide with the relative richness
275
and average spread index hotspot (Figure 6B, F) that includes the provinces Hebei, Liaoning,
276
Shandong, Jiangsu, Henan, Shanxi, Shaanxi and Ningxia, this region should receive
277
particular attention in invasive plant control and risk assessment.
278 279
The identification of areas at high risk of biological invasion can provide the basis for the
280
development of management strategies and policies. Based on our analysis of provincial
281
species richness, relative richness, spread and impact index, we suggest that the highest
282
invasion-risk areas are southern China (Yunnan, Guangxi and Guangdong) and the central 13
ACS Paragon Plus Environment
Environmental Science & Technology
283
Page 14 of 32
coastal areas of eastern China (Shandong, Hebei, Jiangsu).
284 285
Links between AIP impacts and rates of spread
286
The term “invasive” commonly refers to two different but related aspects of AIPs19: the first
287
is the environmental or socio-economic impact that is the chief concern of policy-makers and
288
managers51; the second is the rate at which a species establishes itself and spreads in a new
289
environment. The latter measure of invasiveness is generally used by conservationists and
290
scientists1. The impact index used in this paper is based on both the socio-economic and
291
environmental impacts of a species. Similar indices or species classifications based on the
292
same information, have been used widely in AIP prevention and management at large scales
293
by Chinese scientists and government agencies10, 52, 53. By contrast, indices similar to the
294
spread index used in this paper, based only on natural factors and measuring a species’
295
potential to rapidly colonize in a large area, have been applied in studies involving biological
296
mechanisms of invasions or invasion risk assessments54, but have seldom been applied to
297
large-scale invasive ecology studies in China.
298 299
Generally, the rate of spread of an alien plant should be correlated with its impact18. Most
300
species do not spread widely, nor do they cause substantial environmental or economic
301
damage to a local region. But many studies similar to this one show that individual species’
302
impacts are not strongly correlated with rate of spread (SI 1)
303
invasive species (impact index = 3 or 4) represented only a minority of the total number of
304
AIP and there were no obvious differences in the prevalence of such species across 14
ACS Paragon Plus Environment
19
. The most notorious
Page 15 of 32
Environmental Science & Technology
305
geo-districts (Figure 3A). The seven species with the maximum value of impact index were
306
not widely distributed and did not all spread rapidly (Table 2). For example, Mikania
307
micrantha is a serious weed worldwide: it not only disturbs the growth and development of
308
trees and crops, but also reduces the density of wild native herbaceous plants through
309
allelopathic effects55. Since it first invaded China, however, its distribution has remained
310
confined to three provinces: Guangdong, Guangxi and Yunan. Thus, taking account of
311
potential invasive risk at the national scale, the overall harm caused by Mikania micrantha
312
may be less than that of widely distribution weeds with the same degree of environmental
313
and economic impact, such as Ambrosia artemisiifolia, Ambrosia trifida and Lolium
314
temulentum. Therefore, AIP risk assessments should consider both the socio-economic
315
effects and the biological attributes of species on scales relevant to research and
316
management56.
317 318
Links between species invasiveness and the overall level of invasion
319
Invasiveness is determined by AIPs’ traits and adaptability to different environments.
320
Species with broader tolerance of environmental conditions have the potential to spread
321
further and cause more extensive environmental and economic harm57. We hypothesize that
322
drought tolerance may be a key trait of AIP that have a large-scale distribution in China. In
323
the north of China, where mean annual precipitation is relatively low (Figure 5, Table 1),
324
invasive species are mostly drought-tolerant and the majority of them could potentially
325
expand into southern China (Figure 2, Figure 6 F). In contrast, the tropically distributed
326
AIP tend to be drought-intolerant and are restricted to the south of China. Range expansion 15
ACS Paragon Plus Environment
Environmental Science & Technology
Page 16 of 32
327
seems to be a typical trait of species in response to climate change and resource availability
328
in low-resource environments58, 59. Following the recent trend of increasing drought-prone
329
areas across China60, the potential threat from drought-tolerant invasive species that are
330
currently restricted to temperate China should be given particular attention.
331 332
Community invasibility refers to characteristics of invaded communities that facilitate the
333
invasion of alien species61. Invasibility, in contrast to invasiveness, has rarely been quantified
334
because doing so requires that complex factors like invasion level, species invasiveness and
335
propagule pressure be controlled. Among these factors, invasion level reflects the degree to
336
which local environmental conditions and biotic interactions can accommodate new species9,
337
23
338
of invasion level were used in this study but were correlated differently with environmental
339
and economic factors and therefore convey independent information about the invasion
340
process9,
341
richness was correlated with economic factors (Table 1). The high AIP richness in southern
342
China means that these regions require management and prevention strategies that are more
343
diverse and comprehensive than in other regions (Figure 6A); the high relative AIP richness
344
in the central coastal areas of eastern China suggests that these ecosystems are more
345
vulnerable to plant invasion and that greater management resources need to be mustered here,
346
albeit against fewer species (Figure 6B).
and is best measured by invasive species richness or relative richness6.
62
These two metrics
: richness was correlated more with environmental factors whereas relative
347 348
The negative relationship of average alien plant spread rate to alien plant richness across 16
ACS Paragon Plus Environment
Page 17 of 32
Environmental Science & Technology
349
regions in China (Figure 6A) has two potential explanations. The first is that in a diverse
350
community it is harder for any single species to become abundant because of the complexity
351
of biotic interactions64, 65. On this view, invasive species in southern China are restricted in
352
their distributions by competitors, pathogens and herbivores, whereas invasive species in
353
northern China face fewer such constraints. The second explanation is simply that if one
354
considers all of the propagules of potential invasive species into a given region the species
355
richness of the pool is necessarily inversely related to the average number of propagules per
356
species63, 64. In any case, information on both the total alien richness and individual alien
357
invasive species’ impacts should be included in invasion risk assessments.
358 359
ACKNOWLEDGEMENTS
360
This study was supported by the grants the ‘Public Environment Benefits Programme’ of
361
Ministry of Environment Protection (201209027-2), the ‘111 Program’ of the Bureau of
362
China Foreign Experts and Ministry of Education (2008-B08044), National Basic Research
363
Program of China the Chinese (973 Program 2009CB119200) and the Academy of Science
364
and their Fellowship for International Scientists Programme (2011T2S18). Many critical
365
comments on the manuscript were received in the 2012 Smithsonian CTFS/SIGEO and
366
CForBio Workshop in Seattle (US NSF grant DEB-1046113). We greatly appreciate the many
367
useful criticisms and comments that the reviewers have provided.
368
Supporting Information
369
The List of alien invasive plants in China (SI 1), Literature sources for the list of alien
370
invasive plants in China (SI 2), The evaluated methods of AIP impact and spread indices (SI 17
ACS Paragon Plus Environment
Environmental Science & Technology
371
3) and PCA ordination diagram of environment variables and economic variables (SI 4) are
372
available free of charge via the Internet at http://pubs.acs.org/ in the online version of the
373
paper. The authors are solely responsible for the content and functionality of these materials.
374
Queries should be directed to the corresponding author.
375
REFERENCES
376 377 378 379 380 381 382 383 384 385 386 387 388 389 390 391 392 393 394 395 396 397 398 399 400 401 402 403 404 405 406 407 408 409
(1) Richardson, D.; Pysek, P.; Rejmanek, M.; Barbour, M.; Panetta, F.; West, C., Naturalization and invasion of alien plants: concepts and definitions. Divers. Distrib. 2000, 6, (2), 93-107. (2) Xu Z. H.; Wang Y. P., Disastrous mechanisms and control strategies of alien invasive plants. Chinese J. of Ecology 2004, 23, (3), 124-127. (3)Glowka, L.; Burhenne-Guilmin, F.; Synge, H., A guide to the Convention on Biological Diversity. World Conservation Union, 1994. (4) Balmford, A.; Bennun, L.; Brink, B.; Cooper, D.; Côté, I. M.; Crane, P.; Dobson, A.; Dudley, N.; Dutton, I.; Green, R., ECOLOGY: The Convention on Biological Diversity's 2010 Target. Science. 2005, 307 (5707), 212. (5) Weber, E.; Sun, S. G.; Li, B., Invasive alien plants in China: diversity and ecological insights. Biol Invasions 2008, 10(8), 1411-1429. (6) Xu, H.; Chen, K.; Ouyang, Z.; Pan, X.; Zhu, S., Threats of Invasive Species for China Caused by Expanding International Trade. Environ. Sci. Technol. 2012, 46(13), 7063-7064. (7) Pysek, P.; Richardson, D. M., The biogeography of naturalization in alien plants. J. Biogeogr. 2006, 33(12), 2040-2050. (8) Pysek, P.; Richardson, D. M.; Jarosik, V., Who cites who in the invasion zoo: insights from an analysis of the most highly cited papers in invasion ecology. Preslia 2006, 78(4), 437-468. (9) Catford, J. A.; Vesk, P. A.; Richardson, D. M.; Pysek, P., Quantifying levels of biological invasion: towards the objective classification of invaded and invasible ecosystems. Global Change Biol. 2012, 18(1), 44-62. (10) Liu, J.; Liang, S.; Liu, F.; Wang, R.; Dong, M., Invasive alien plant species in China: regional distribution patterns. Divers. Distrib. 2005, 11(4), 341-347. (11) Acosta, A.; Izzi, C. F.; Stanisci, A., Comparison of native and alien plant traits in Mediterranean coastal dunes. Community Ecol. 2006, 7(1), 35-41. (12) van Kleunen, M.; Weber, E.; Fischer, M., A meta-analysis of trait differences between invasive and non-invasive plant species. Ecol. Lett. 2010, 13(2), 235-245. (13) Shi, F. C.; Yuan, Z. N.; He, K. S.; Li, J. J.; Guo, S. Y., The Ecological Profit in Phenology and Pollination of Alien Species Spartina Alterniflora (Loisel) Introduced to Coastal Wetlands in Northern China. Pol. J. Ecol. 2009, 57(2), 283-289. (14) Lonsdale, W. M., Global patterns of plant invasions and the concept of invasibility. Ecology 1999, 80(5), 1522-1536. (15) Alpert, P.; Bone, E.; Holzapfel, C., Invasiveness, invasibility and the role of environmental stress in the spread of non-native plants. Perspect. Plant Ecol. 2000, 3(1), 52-66. 18
ACS Paragon Plus Environment
Page 18 of 32
Page 19 of 32
410 411 412 413 414 415 416 417 418 419 420 421 422 423 424 425 426 427 428 429 430 431 432 433 434 435 436 437 438 439 440 441 442 443 444 445 446 447 448 449 450 451 452 453
Environmental Science & Technology
(16) Richardson, D. M.; Pyšek, P., Plant invasions: merging the concepts of species invasiveness and community invasibility. Prog. Phys. Geog. 2006, 30(3), 409-431. (17) Catford, J. A.; Jansson, R.; Nilsson, C., Reducing redundancy in invasion ecology by integrating hypotheses into a single theoretical framework. Divers. Distrib. 2009, 15(1), 22-40. (18) McIntyre, S.; Martin, T. G.; Heard, K. M.; Kinloch, J., Plant traits predict impact of invading species: an analysis of herbaceous vegetation in the subtropics. Aust. J. Bot. 2005, 53(8), 757-770. (19) Ricciardi, A.; Cohen, J., The invasiveness of an introduced species does not predict its impact. Biol. Invasions 2007, 9(3), 309-315. (20) Culley, T. M.; Hardiman, N. A., The role of intraspecific hybridization in the evolution of invasiveness: a case study of the ornamental pear tree Pyrus calleryana. Biol. Invasions 2009, 11(5), 1107-1119. (21) Magee, T. K.; Ringold, P. L.; Bollman, M. A.; Ernst, T. L., Index of Alien Impact: A Method for Evaluating Potential Ecological Impact of Alien Plant Species. Environ. Manage. 2010, 45(4), 759-778. (22) Hou, Y. P.; Hu, Q. B.; Wan, F. H.; You, M. S., A conspectus of basic research on biological invasions in China supported by the National Natural Science Foundation. J. Biosafety 2011, 20(1), 20-28. (23) Crooks, J., Characterizing ecosystem-level consequences of biological invasions: the role of ecosystem engineers. Oikos 2002, 97(2), 153-166. (24) Feng, M. J.; Dong, X. D.; Xu, C. D., Spatial patterns of floristic composition of invasive alien plants in large scale and their climatic interpretation. J. Wuhan Botanical Res 2009, 27(2), 159-164. (25) Sang, W.; Zhu, L.; Axmacher, J., Invasion pattern of Eupatorium adenophorum Spreng in southern China. Biol. Invasions 2010, 12(6), 1721-1730. (26) Feng, J.; Zhu, Y., Alien invasive plants in China: risk assessment and spatial patterns. Biodivers .Conserv. 2010, 19(12), 3489-3497. (27) Thuiller, W.; Richardson, D. M.; PYŠEK, P.; Midgley, G. F.; Hughes, G. O.; Rouget, M., Niche‐based modelling as a tool for predicting the risk of alien plant invasions at a global scale. Global Change Biol 2005, 11, (12), 2234-2250. (28) Wu, Z. Y.; Raven, P. H.; Hong, D. Y., Flora of China. Science Press & Missouri Botanical Garden Beijing & St. Louis, 1994-2007. (29) Ding, J. Q.; Wang, R., Invasive alien plant species and their impact on biodiversity in China. In China's Biodiversity: A Country Study, Compilation Group of China' s Biodiversity, Ed. China Environmental Science Press: Beijing, 1998, 58-63. (30) Li, Z. Y.; Xie, Y., Invasive alien species in China. China Forestry Publishing House: Beijing, 2002. (31) Liu, J.; Dong, M.; Miao, S. L.; Li, Z. Y.; Song, M. H.; Wang, R. Q., Invasive alien plants in China: role of clonality and geographical origin. Biol. Invasions 2006, 8(7), 1461-1470. (32) Huang, Q. Q.; Wu, J. M.; Bai, Y. Y.; Zhou, L.; Wang, G. X., Identifying the most noxious invasive plants in China: role of geographical origin, life form and means of introduction. Biodivers. Conserv. 2009, 18(2), 305-316. (33) Wu, S. H.; Sun, H. T.; Teng, Y. C.; Rejmánek, M.; Chaw, S. M.; Yang, T. Y. A.; Hsieh, C. F., Patterns of plant invasions in China: Taxonomic, biogeographic, climatic approaches and anthropogenic effects. Biol. Invasions 2010, 12(7), 2179-2206. (34) Jiang, H.; Fan, Q.; Li, J. T.; Shi, S.; Li, S. P.; Liao, W. B.; Shu, W. S., Naturalization of alien plants in China. Biodivers. Conserv. 2011, 1-12. (35) ter Braak, C., The analysis of vegetation environment relationships by canonical correspondence analysis. Vegetatio 1987, 69, 69-77. (36) Naeem, S.; Knops, J.; Tilman, D.; Howe, K.; Kennedy, T.; Gale, S., Plant diversity increases resistance to 19
ACS Paragon Plus Environment
Environmental Science & Technology
454 455 456 457 458 459 460 461 462 463 464 465 466 467 468 469 470 471 472 473 474 475 476 477 478 479 480 481 482 483 484 485 486 487 488 489 490 491 492 493 494 495 496 497
invasion in the absence of covarying extrinsic factors. Oikos 2000, 91(1), 97-108. (37) Sax, D., Latitudinal gradients and geographic ranges of exotic species: implications for biogeography. J. Biogeogr. 2001, 28, (1), 139-150. (38) Kennedy, T. A.; Naeem, S.; Howe, K. M.; Knops, J. M. H.; Tilman, D.; Reich, P., Biodiversity as a barrier to ecological invasion. Nature 2002, 417(6889), 636-638. (39) Stohlgren, T.; Otsuki, Y.; Villa, C.; Lee, M.; Belnap, J., Patterns of plant invasions: a case example in native species hotspots and rare habitats. Biol. Invasions 2001, 3(1), 37-50. (40) Feng, J.; Zhang, Z.; Nan, R., The roles of climatic factors in spatial patterns of alien invasive plants from America into China. Biodivers. Conserv. 2011, 1-7. (41) Huang, J. H.; Chen, B.; Liu, C.; Lai, J.; Zhang, J.; Ma, K. P., Identifying hotspots of endemic woody seed plant diversity in China. Divers. Distrib. 2011. (42) Smith, M. D.; Knapp, A. K., Exotic plant species in a C-4-dominated grassland: invasibility, disturbance, and community structure. Oecologia 1999, 120(4), 605-612. (43) Dukes, J. S., Biodiversity and invasibility in grassland microcosms. Oecologia 2001, 126(4), 563-568. (44) Meng, T. T.; Ni, J.; Harrison, S. P., Plant morphometric traits and climate gradients in northern China: a meta-analysis using quadrat and flora data. Ann. Bot-London 2009, 104(6), 1217-1229. (45) Rodgers, J. C.; Parker, K. C., Distribution of alien plant species in relation to human disturbance on the Georgia Sea Islands. Divers. Distrib. 2003, 9(5), 385-398. (46) Westphal, M. I.; Browne, M.; MacKinnon, K.; Noble, I., The link between international trade and the global distribution of invasive alien species. Biol. Invasions 2008, 10, (4), 391-398. (47) Hellmann, J. J.; Byers, J. E.; Bierwagen, B. G.; Dukes, J. S., Five potential consequences of climate change for invasive species. Conserv. Biol. 2008, 22(3), 534-543. (48) Ni, J.; Sykes, M.; Prentice, I.; Cramer, W., Modelling the vegetation of China using the process-based equilibrium terrestrial biosphere model BIOME3. Global Ecol. Biogeogr. 2000, 463-479. (49) Walther, G.-R.; Post, E.; Convey, P.; Menzel, A.; Parmesan, C.; Beebee, T. J. C.; Fromentin, J.-M.; Hoegh-Guldberg, O.; Bairlein, F., Ecological responses to recent climate change. Nature 2002, 416(6879), 389-395. (50) Yu, M.; Gao, Q.; Xu, H. C.; Liu, Y. H., Responses of vegetation distribution and primary production of the terrestial ecosystems of China to climatic change. Quaternary Sci. 2001, 21, 281-293. (51) Kohli, R. K.; Batish, D. R.; Singh, H. P.; Dogra, K. S., Status, invasiveness and environmental threats of three tropical American invasive weeds (Parthenium hysterophorus L., Ageratum conyzoides L., Lantana camara L.) in India. Biol. Invasions 2006, 8(7), 1501-1510. (52) OTA, Harmful non-indigenous species in the United States. U.S. Government Printing Office: Washington, D.C., USA., 1993. (53) Xu, H. G.; Qiang, S.; Han, Z. M.; Guo, J. Y.; Huang, Z. G.; Sun, H. Y.; He, S. P.; Ding, H.; Wu, H. H.; Wan, F. H., The distribution and introduction pathway of alien invasive species in China. Biodiversity Sci. 2004, 12, (6), 626-638. (54) Lachmuth, S.; Durka, W.; Schurr, F. M., The making of a rapid plant invader: genetic diversity and differentiation in the native and invaded range of Senecio inaequidens. Mol. Ecol. 2010, 19, (18), 3952-3967. (55) Hou, Y. P.; Peng, S. L.; Chen, B. M.; Ni, G. Y., Inhibition of an invasive plant (Mikania micrantha H.B.K.) by soils of three different forests in lower subtropical China. Biol. Invasions 2011, 13, (2), 381-391. (56) Robertson, M. P.; Villet, M. H.; Fairbanks, D. H. K.; Henderson, L.; Higgins, S. I.; Hoffmann, J. H.; Le Maitre, D. C.; Palmer, A. R.; Riggs, I.; Shackleton, C. M.; Zimmermann, H. G., A proposed prioritization system for the management of invasive alien plants in South Africa. S. Afr. J. Sci. 2003, 99, (1-2), 37-43. 20
ACS Paragon Plus Environment
Page 20 of 32
Page 21 of 32
498 499 500 501 502 503 504 505 506 507 508 509 510 511 512 513 514 515 516
Environmental Science & Technology
(57) Anderson, B.; Akcakaya, H. R.; Araujo, M. B.; Fordham, D. A.; Martinez-Meyer, E.; Thuiller, W.; Brook, B. W., Dynamics of range margins for metapopulations under climate change. P. Roy. Soc. B-Biol. Sci. 2009, 276, (1661), 1415-1420. (58) Engelkes, T.; Morriën, E.; Verhoeven, K. J. F.; Bezemer, T. M.; Biere, A.; Harvey, J. A.; McIntyre, L. M.; Tamis, W. L. M.; Van der Putten, W. H., Successful range-expanding plants experience less above-ground and below-ground enemy impact. Nature 2008, 456(7224), 946-948. (59) Funk, J. L.; Vitousek, P. M., Resource-use efficiency and plant invasion in low-resource systems. Nature 2007, 446(7139), 1079-1081. (60) Zou, X. K.; Ren, G. Y.; ZHANG, Q., Droughts variations in China based on a compound index of meteorological drought. Clim. Environ. Res. 2010, 15(4), 371-378. (61) Belote, R. T.; Jones, R. H.; Hood, S. M.; Wender, B. W., Diversity-invasibility across an experimental disturbance gradient in Appalachian forests. Ecology 2008, 89(1), 183-192. (62) Thompson, K.; Petchey, O. L.; Askew, A. P.; Dunnett, N. P.; Beckerman, A. P.; Willis, A. J., Little evidence for limiting similarity in a long-term study of a roadside plant community. J. Ecol. 2010, 98(2), 480-487. (63) Colautti, R. I.; Grigorovich, I. A.; MacIsaac, H. J., Propagule pressure: A null model for biological invasions. Biol Invasions 2006, 8, (5), 1023-1037. (64) Trueman, M.; Atkinson, R.; Guezou, A.; Wurm, P., Residence time and human-mediated propagule pressure at work in the alien flora of Galapagos. Biol. Invasions 2010, 12(12), 3949-3960.
517
21
ACS Paragon Plus Environment
Environmental Science & Technology
518
Page 22 of 32
TABLES
519 520
Table 1 Results of stepwise linear regressions (adjusted standardized coefficients; β) of
521
Chinese provinces’ invasion indices on economic and environmental variables.
EV
Richness
R_rich
T_IM
T_IS
A_IM
A_IS
MAT
0.675**
-
0.573**
0.769**
0.621**
-
MAP
-
-
-
-
-
-0.782**
GDP
0.263*
0.579**
-
-
-
-
RM
-
-
0.396**
-
0.369*
-
85.272 + Final equation
4.587MAT + 0.001GDP R2=0.609**
522 523 524 525 526 527 528
0.028 + -7
8.806×10 GDP 2
R =0.309**
157.226 + 8.555MAT + 0.003 RM R2=0.639**
257.660 + 9.088MAT 2
R =0.591**
1.671 + 0.004MAT + 0.96 RM R2=0.336**
** Correlation is significant at the 0.01 level (2-tailed). * Correlation is significant at the 0.05 level (2-tailed). EV: environment variables; Richness: Invasive species richness; R_rich: Relative invasive species richness; T_IM: Provincial total value of impact index; T_IS: Provincial total value of spread index; A_IM: Provincial weighted average of impact index; A_IS: Provincial weighted average of spread index; MAT: mean annual temperature; MAP: mean annual precipitation; GDP: gross national product; RM: Railway mileage.
529
22
ACS Paragon Plus Environment
2.978 0.0003MAP R2=0.597**
Page 23 of 32
Environmental Science & Technology
530
Table 2 Alien invasive plants with highest scores on the impact index†.
531
532
†
†
†
Scientific name
Impact index
Spread index
Ambrosia artemisiifolia L. Ambrosia trifida L. Cenchrus echinatus Linnaeus Eichhornia crassipes (Martius) Solms Eupatorium odoratum L. Lolium temulentum Linnaeus Mikania micrantha Kunth
4 4 4 4 4 4 4
3 4 2 3 2 4 1
For definitions of the impact and spread indices, see text.
23
ACS Paragon Plus Environment
Environmental Science & Technology
533
FIGURE LEGENDS
534
Figure 1: Administrative areas and geographic zones in China. NE: Northeast (1:
535
Heilongjiang, 2: Jilin, 3: Liaoning); NC: North center (4: Hebei, 5: Neimenggu, 6:
536
Shanxi); NW: Northwest (7: Shaanxi, 8: Ningxia, 9: Gansu, 10: Qinghai, 11:
537
Xinjiang); SE: Southeast (12: Shandong, 13: Anhui, 14: Jiangsu, 15: Jiangxi, 16:
538
Zhejiang, 17: Fujian, 18: Taiwan); SC: South center (19: Henan, 20: Hubei, 21:
539
Hunan, 22: Guangdong, 23: Guangxi, 24: Hainan); SW: Southwest (25: Guizhou, 26:
540
Sichuan, 27: Yunnan, 28: Xizang).
541 542
Figure 2: Distribution types of alien invasive plant (AIP) in China (CO:
543
Cosmopolitan species; OS: Species distributed only in southern provinces; PS:
544
Species distributed predominantly in southern provinces; PN: Species distributed
545
predominantly in northern provinces; ON: Species distributed only in northern
546
provinces). Woody AIP tend to be more southerly distributed than herbs (p < 0.001;
547
chi-square (χ2) test).
548 549
Figure 3: Proportions of alien invasive plants with each impact score ( A) and spread
550
score ( B) across geo-districts in China. Differences in the distribution of impact or
551
spread scores across six geo-districts’ were tested with chi-square (χ2) tests.
552 553
NE:Northeast; NC: North center; NW: Northwest; SE: Southeast; SC: South center;
554
SW: Southwest.
24
ACS Paragon Plus Environment
Page 24 of 32
Page 25 of 32
Environmental Science & Technology
555 556
Figure 4: Alien invasive plant richness versus total species richness across Chinese
557
provinces.
558 559
SE: Southeast; SC: Southcenter; SW: Southwest; NE: Northeast; NC: Northcenter;
560
NW: Northwest.
561 562
Figure 5: Canonical Correspondence Analysis ordination diagram of alien invasive
563
plant species presence-absence data in China showing correlations with
564
environmental and economic variables. Each point represents a single species (but
565
note that all cosmopolitan (CO) species have the same provincial distribution and
566
thus map to a single point).
567
MAT: mean annual temperature; MAP: mean annual precipitation; GDP: gross
568
national product; RM: railway mileage.
569
OS: Species distributed only in southern provinces; PS: Species distributed
570
predominantly in southern provinces; PN: Species distributed predominantly in
571
northern provinces; ON: Species distributed only in northern provinces
572 573
Figure 6: Distributions of aggregate invasion indices by provinces across China
574
A: Invasive species richness; B: Relative invasive species richness; C: Provincial
575
total value of impact index (T_IM); D: Provincial total value of spread index (T_IS);
576
E: Provincial weighted average of impact index (A_IM); F: Provincial weighted
25
ACS Paragon Plus Environment
Environmental Science & Technology
577
average of spread index (A_IS); 1: Heilongjiang; 2: Jilin; 3: Liaoning; 4: Hebei; 5:
578
Neimenggu; 6: Shanxi; 7: Shaanxi; 8: Ningxia; 9: Gansu; 10: Qinghai; 11: Xinjiang;
579
12: Shandong; 13: Anhui; 14: Jiangsu; 15: Jiangxi; 16: Zhejiang; 17: Fujian; 18:
580
Taiwan; 19: Henan; 20: Hubei; 21: Hunan; 22: Guangdong; 23: Guangxi; 24: Hainan;
581
25: Guizhou; 26: Sichuan; 27: Yunnan; 28: Xizang.
582 583
FIGURE 1
584
26
ACS Paragon Plus Environment
Page 26 of 32
Page 27 of 32
Environmental Science & Technology
585 586
FIGURE 2
587 588 589
27
ACS Paragon Plus Environment
Environmental Science & Technology
590 591
FIGURE 3
592
28
ACS Paragon Plus Environment
Page 28 of 32
Page 29 of 32
Environmental Science & Technology
593 594 595
FIGURE 4
596
29
ACS Paragon Plus Environment
Environmental Science & Technology
597 598
FIGURE 5
599
30
ACS Paragon Plus Environment
Page 30 of 32
Page 31 of 32
Environmental Science & Technology
600 601
FIGURE 6
602
31
ACS Paragon Plus Environment
Environmental Science & Technology
Does alien invasive richness correlate with areas of the highest invasive risk in China? 450x306mm (100 x 100 DPI)
ACS Paragon Plus Environment
Page 32 of 32