Spatial Risk Assessment of Alien Invasive Plants in China

Smithsonian Tropical Research Institute, P.O. Box 0843-03092, Balboa, Ancón, Republic of Panamá. Environ. Sci. Technol. , 2013, 47 (14), pp 7624–7...
1 downloads 6 Views 880KB Size
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