Spatially Explicit Analysis of Biodiversity Loss Due to Global

Feb 25, 2016 - Influence of spatial differentiation in impact assessment for LCA-based decision support: Implementation of biochar technology in Indon...
0 downloads 0 Views 1MB Size
Subscriber access provided by EMORY UNIV

Article

Spatially Explicit Analysis of Biodiversity Loss due to Global Agriculture, Pasture and Forest Land Use from a Producer and Consumer Perspective Abhishek Chaudhary, Stephan Pfister, and Stefanie Hellweg Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.5b06153 • Publication Date (Web): 25 Feb 2016 Downloaded from http://pubs.acs.org on February 28, 2016

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 27

Environmental Science & Technology

3

Spatially Explicit Analysis of Biodiversity Loss due to Global Agriculture, Pasture and Forest Land Use from a Producer and Consumer Perspective

4

Abhishek Chaudhary1*, Stephan Pfister1, and Stefanie Hellweg1

1 2

5 6 7

1

8

ABSTRACT

9

Anthropogenic land use to produce commodities for human consumption is the major driver of

10

global biodiversity loss. Synergistic collaboration between producers and consumers in needed to

11

halt this trend. In this study, we calculate species loss on 5 min × 5 min grid level and per

12

country due to global agriculture, pasture and forestry by combining high-resolution land use

13

data with countryside species area relationship for mammals, birds, amphibians, and reptiles.

14

Results show that pasture was the primary driver of biodiversity loss in Madagascar, China and

15

Brazil, while forest land use contributed the most to species loss in DR Congo and Indonesia.

16

Combined with the yield data, we quantified the biodiversity impacts of 1 m3 of roundwood

17

produced in 139 countries, concluding that tropical countries with low timber yield and a large

18

presence of vulnerable species suffer the highest impact. We also calculated impacts per kg for

19

160 crops grown in different countries and linked it with FAO food trade data to assess the

20

biodiversity impacts embodied in Swiss food imports. We found that more than 95% of Swiss

21

consumption impacts rest abroad with cocoa, coffee and palm oil imports being responsible for

22

majority of damage.

Institute of Environmental Engineering, ETH Zurich, 8093 Zurich, Switzerland, (*Corresponding author phone: 41-44-6330254; fax: 41-44-6331061; e-mail: [email protected])

23 24 25

ACS Paragon Plus Environment

Environmental Science & Technology

26

INTRODUCTION

27

Global biodiversity faces unprecedented extinction crisis1 mainly owing to habitat loss caused by

28

conversion of natural forests to anthropogenic land use.2 Growing population, increasing per

29

capita consumption, a shift to meat-based diets and biofuel production are the leading factors

30

responsible for this trend.3 More than 75% of the Earth’s ice-free land is already being affected

31

by human activities.4 The FAO forecasts a ∼70% increase in global food demand from 2000 to

32

2050, perhaps leading to further habitat loss to make way for additional cropland.5 In addition to

33

deforestation for agriculture and pasture land use, forest exploitation for commercial purposes is

34

another threat to biodiversity especially in the tropics.6

35

International efforts aimed at halting rate of global biodiversity loss have failed to meet their

36

targets.1,7 Apart from traditional measures such as setting aside areas rich in biodiversity for

37

conservation purposes,8 novel policies aimed at directly addressing the human drivers of

38

biodiversity loss (e.g. reducing the intensity of damaging land use types or changing

39

consumption patterns) need to be implemented in parallel.9 Identification of hotspots and the

40

land use types causing high biodiversity damage can help producer nations develop local

41

strategies and control further damage. On the other hand, as the international trade of food and

42

timber products continues to increase,10 informing the consumers regarding the environmental

43

impacts ‘hidden’ behind the imported products and identifying commodities causing high

44

damage in producing nations may induce sustainable consumption patterns and help design

45

demand-side mitigation measures.

46

Life cycle assessment (LCA) is increasingly used to evaluate the cradle to grave environmental

47

impacts of products.11 The main advantage of LCA studies over traditional impact assessments is

48

that it establishes a link between the final commodity (e.g. one kg of tea) and the associated

ACS Paragon Plus Environment

Page 2 of 27

Page 3 of 27

Environmental Science & Technology

49

biodiversity loss.12 However, within LCA, biodiversity impacts due to land use are often poorly

50

quantified.13-15

51

Land use impact assessment within LCA

52

Traditional methods for land use impact assessment within LCA are often based on limited

53

biodiversity datasets from specific world regions (such as Europe) or taxa (mainly vascular

54

plants).16 Only recently, globally applicable and spatially differentiated methods have been made

55

available. de Baan et al.17 were the first to provide local characterization factors (CFs, i.e. the

56

factors indicating species loss caused by unit area of a particular land) for different taxa in

57

different terrestrial biomes. However, local CFs do not inform regarding the contribution of land

58

use towards potential irreversible, global extinction of rare and threatened species due to habitat

59

loss/degradation. Avoiding global species extinctions is important to preserve the evolutionary

60

and genetic diversity of life on Earth. For predicting regional and global species extinctions

61

resulting from habitat loss, species area-relationships (SARs) have commonly been used in

62

LCA.18,19

63

Countryside SAR has recently been shown to perform better than other SAR models in

64

predicting species extinction from habitat loss in a heterogeneous, human modified landscape

65

(i.e. the countryside) and recognizes the fact that species adapted to human-modified habitats

66

also survive in the absence of natural habitat.20 Recently Chaudhary et al.21 provided ecoregion

67

specific regional CFs for six land use types using countryside SAR (e.g. regional species lost per

68

m2 of annual crop) for different taxa. The global CFs were then calculated by weighting the

69

regional CFs with a ‘vulnerability score’ (VS) of that ecoregion.21

70

Global scale land use impact assessment

ACS Paragon Plus Environment

Environmental Science & Technology

71

Many studies have assessed the land use impacts of individual crops or food items from

72

particular countries using different methods. For example, Mattsson et al.22 assessed Swedish

73

rapeseed, Brazilian soybean and Malaysian palm oil, de Baan et al.12 assessed impacts of tea,

74

coffee and tobacco in East Africa, Chaudhary et al.21 assessed sugarcane, wheat, sugarbeet and

75

maize from Brazil, France and USA. However, an explicit assessment of species loss due to all

76

crops from all countries is lacking. Similarly, for impacts due to forest or pasture land use,

77

studies have focused on individual countries, e.g. for Norway23 or Ghana24, while a global

78

analysis is lacking. In order to assess the impacts of different land use types on biodiversity with

79

global coverage, spatially differentiated land use inventory data need to be combined with

80

regionalized characterization factors.25

81

On the inventory side, global maps of agricultural crops and pasture land use have been compiled

82

by Monfreda et al.26 at 5 × 5 arc minute grid level. For forest land use, global maps also at 5 arc

83

minute grid level are available.27,28 The advantage of these maps is that they combine ‘land

84

cover’ data from remote sensing with statistical data on human activities, and thereby provide

85

information on ‘land use’ (e.g. forests for timber extraction) and not only land cover (e.g.

86

coniferous forests - as given by Hansen et al.29). On the impact assessment side, the

87

characterization factors (CFs) provided by Chaudhary et al.21 can be used to assess impacts of

88

different land use classes such as arable or permanent crops on a global scale.

89

Objectives & Scope

90

The overall goal of manuscript is to provide a spatially-explicit analysis of land-use driven

91

biodiversity loss from worldwide agriculture, pasture and forestry, using the recently published

92

impact assessment method of Chaudhary et al.21in combination with high resolution land use

93

maps. From a producer’s perspective, we identify the most damaging land use types causing high

ACS Paragon Plus Environment

Page 4 of 27

Page 5 of 27

Environmental Science & Technology

94

species loss for mammals, birds, amphibians and reptiles globally and also on a regional scale at

95

5 arc minute resolution. We then aggregate the results to country level and in combination with

96

crop and timber yield data, quantify species loss per kg of 160 global crops and per m3 of

97

roundwood production. These results can be applied in combination with trade data to assess

98

biodiversity impacts embodied in traded food, fibre and wood items and to determine

99

biodiversity impacts from a consumer perspective. To illustrate how the newly calculated

100

country-crop impact factors can be used to quantify the impacts imported and exported by

101

different countries, we carry out a case-study on Swiss food consumption.

102

MATERIALS AND METHODS

103

Land use biodiversity impacts at 5 arc minute resolution

104

We imported the global land occupation characterization factors (CFs) per ecoregion for the land

105

use types - annual crops, permanent crops, pasture, extensive forest and intensive forest from

106

Chaudhary et al.21. The global CFs were derived using countryside SAR and weighted with

107

vulnerability score for each of the four vertebrate taxa (mammals, birds, amphibians and reptiles)

108

in 804 terrestrial ecoregions. The global CFs are higher for ecoregions hosting more endemic and

109

threatened species and provide an estimate of permanent (irreversible) species extinction caused

110

by an unit area of land use.21 See supporting information-1 (SI-1) for more details and equations

111

used for calculation of CFs.

112

We combine these CFs with high-resolution (5 x 5 arc minute) land use maps of global

113

agriculture26, pasture26 and forest27,28 land use to calculate biodiversity impacts per grid cell

114

(pixel). It was assumed that value of characterization factors in each pixel, CFp, is the same for

115

all pixels p situated within an ecoregion j (CFj). The biodiversity damage for each taxon g due to

116

land use type i per grid cell (,, ) is given by:

ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 27

,, =   ,,, × , 117

Equation 1

118

Where , is the area occupied by land use type i in the pixel p (in m2·years). The global CFs

119

(   ,,, ) provide biodiversity damage in the units- global species equivalents lost per m2 of

120

land use type i in pixel p (hereafter species eq. lost/m2, see Chaudhary et al.21 for full details).

121

Therefore, equation-1 provides the biodiversity damage (,, ) caused by land use type i in

122

units of species eq. lost·years for each taxonomic group g.21 Additionally, we also assessed the

123

impacts using the taxa-aggregated CFs for each land use type per ecoregion provided by

124

Chaudhary et al.21 that give the biodiversity impacts in the units- potentially disappeared fraction

125

per m2 of land use (PDF/m2), which can be used in LCA studies (see SI-1 for details).21

126

The pasture and agricultural land area inside each pixel was imported from Monfreda et al.26

127

while the area occupied by forest land use per pixel was taken as the average of managed forest

128

area provided by global land use maps ANTHROME28 and LADA27, both at 5 arc minute grid

129

level.

130

For global agricultural land use, the harvested area and yield maps are available for each of 160

131

crops on a 5 min grid level from Monfreda et al.26 Pfister et al.30 adjusted this crop area per pixel

132

for multiple cropping, using length of growing season estimates of each crop in different agro-

133

ecological zones. We used these adjusted values for area occupied by a crop c in each pixel p

134

(denoted as , hereafter). The biodiversity damage (,, ) per taxa g due to crop c is

135

calculated by multiplying the characterization factor of agriculture land with the area occupied

136

by the crop within each pixel p: ,, =   ,,, × ,

137

Equation 2

ACS Paragon Plus Environment

Page 7 of 27

Environmental Science & Technology

138

Land use biodiversity impacts per country

139

The total biodiversity damage caused by each land use type in a country is calculated as the sum

140

of biodiversity impacts per pixel ,, calculated in equation 1 above: 







,, =  ,, =    ,,, ∙ , 141

Equation 3

142

Here n is the total number of pixels within the country k. Analogously, the total biodiversity

143

damage caused by each of the 160 crops in each country (,, ) is obtained by summing the

144

impacts per pixel (,, ) calculated in equation 2 above.

145

Country-specific characterization factors

146

As life cycle inventory and global trade databases often report the country but not the ecoregion

147

of production, we also calculated the CFs per country which will be more convenient than those

148

at ecoregion level previously provided by Chaudhary et al.21We calculated the country-specific

149

CFs for the land use types – extensive forestry, intensive forestry and pasture for all countries k

150

and taxa g by dividing the total impact due to these land use types by the total area occupied by

151

the them in country k.   ,,,

∑   ,,, ∙ , = ∑ ,

152

Equation 4

153

Here n denotes the number of pixels within the country k,   ,,, are CFs on pixel level

154

derived by Chaudhary et al.21 Similarly, we also calculated country-specific CFs for each of 160

155

global crops by dividing their total impact (from equation-2 above) by the total area occupied by

156

them in country k.

157

Biodiversity impacts per m3 of roundwood production per country

ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 27

158

We calculate the biodiversity impacts per m3 of roundwood by dividing the country-specific

159

extensive and intensive forestry CFs (species eq. lost per m2, from equation 4) with the

160

corresponding harvesting intensity (m3/ha/year) in the units species eq. lost·years per m3: ,,  =

  ,

!,,

× 10000 × (1 − % , )   , ,, × 10000 × % , + &' !, &' ,

161

Equation 5

162

Where index ext denotes extensive forestry, pla for planted (intensively managed) forests. The

163

forest plantations often consisting of fast growing timber species were assumed to fall in

164

intensive forestry category as they simplify the forest structure leading to higher biodiversity

165

damage.6 All other types of natural managed forests were considered in extensive forestry

166

category. Wood harvesting intensity per country (&' ) in m3/ha/year was derived for both

167

managed natural forests and planted forests by dividing the total annual roundwood production

168

by their respective area per country using values from FAO forest resource assessment report.31

169

The fraction of a country’s roundwood coming from planted forests (% , ) was obtained from

170

Jürgensen et al.32

171

Biodiversity impacts per kg of crop

172

The biodiversity damage per kg of crop c of country k (species eq. lost·years per kg) is calculated

173

as the total impact due to crop land use in that country (equation 3) divided by the total crop

174

production (in kg):  ,,,  =

,,  ∑ ), ∙ ,

175

Equation 6

ACS Paragon Plus Environment

Page 9 of 27

Environmental Science & Technology

176

Here n is the total pixels within the country k, , is the area of crop c in pixel p (ha) and ), is

177

the yield (kg/ha) for crop c in pixel p, also obtained from Monfreda et al.26

178

Case-study: Swiss food consumption impacts

179

The weight of each food item imported by Switzerland in 2011 was obtained from FAOSTAT

180

trade database.33 For processed food products, we used conversion factors provided by FAO34 to

181

estimate the weight of crop required to produce the item (e.g. 4 kg of oranges required per kg of

182

orange juice).

183

However, the FAO data only reports the last country from which the food item is traded but not

184

the actual country where the item was produced. An example is Belgium as an intermediate

185

trader, which exports tea to a lot of EU countries but does not produce it. We assumed that if a

186

country produces the exported crop then the land use occurred there. However if it does not

187

produce the product (e.g. tea from Belgium), the imported quantity was allocated to the biggest

188

producers of this product worldwide in the same proportion as their global export share (data

189

from FAOSTAT).33

190

More sophisticated approaches exist to trace the actual country of origin such as MRIO databases

191

that utilize monetary transaction and economic data (e.g. Eora9, EXIOBASE35, or Kastner et

192

al.36). However, one of the shortcomings of existing MRIO tables is the grouping of multiple

193

crops (e.g. staple crops) or regions (e.g. South-east Asia) into a common category. As the crops

194

of one category might differ significantly in terms of their environmental impacts or the

195

countries within one big region might differ in terms of impacts per unit land use, such grouping

196

might lead to under/over estimation of their impacts. We used an alternative approach and utilize

197

the physical supply chain data available from FAOSTAT33 to allocate the imported food item to

198

correct for intermediate trading countries. We only considered the food items directly derived

ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 27

199

from crops and did not consider the imported livestock products as this was beyond the scope of

200

this study.

201

Once the food item was allocated to countries of origin, the biodiversity damage associated with

202

its import was calculated by multiplying its mass (kg) with newly calculated per kg impacts for

203

that combination of crop and country (equation 6). The imported impacts were also compared

204

with impacts occurring due to net agriculture land use within Switzerland for Swiss consumption

205

(i.e. total agricultural land use minus land use for exports).

206

RESULTS

207

Biodiversity damage due to crop, pasture and forest land use

208

Figure 1(a) shows the total global mammal species loss per pixel due to agriculture, pasture and

209

forestry land use combined. Geographic hotspots of land-use driven mammal species loss are

210

located in Indonesia, Madagascar, Philippines, Brazil, Papua New Guinea, China, India, DR

211

Congo and Mexico. Figures 1(b) and 1(c) show the zoomed in maps for Brazil and Indonesia

212

respectively to help identify hotspots of mammal species loss within these countries at a 5 arc

213

minute resolution. Managed forest is the main driver of mammal species loss in Indonesia,

214

whereas pasture is most damaging land use type in Brazil. Agriculture land use was identified as

215

main driver in Philippines India, and Sri Lanka. See supporting information-1 for all maps

216

showing impacts due to each land use type per taxa in each pixel.

217

ACS Paragon Plus Environment

Page 11 of 27

Environmental Science & Technology

218

219 220

Figure 1. a). Total mammal impacts (global species eq. lost·years) at 5 arc minute resolution due

221

to current anthropogenic land occupation. Impacts per pixel due to agriculture, pasture and forest

222

land use were calculated separately using equation 1 and then summed, b). Zoomed in map of

223

Brazil showing high damage areas lie mostly along the east coast and in the Andes, c). Zoomed

224

in map of Indonesia showing global hotspots of mammal species loss (dark brown pixels). See

225

supporting information-1 for all maps showing impacts due to each land use type per taxa.

226

The impact due to total agriculture land use per country for all taxa was calculated using

227

equation 3. Table S1 of supporting information-2 (SI-2) shows the detailed results of biodiversity

228

loss due to each land use type per country and taxa. Here again, although countries with large

ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 27

229

agriculture area were expected to incur high impacts, smaller countries such as Sri Lanka,

230

Malaysia, Philippines, all featured in top 10 countries suffering high biodiversity loss due to high

231

species richness and vulnerability scores for ecoregions within these countries (Table S1, SI-2).

232

Table S1 further shows that for pasture land use, highest mammal species loss was found to

233

occur in Madagascar, China, Brazil, Australia and Colombia causing ~45% of the total mammal

234

species loss due to global pasture land use. Together these five countries account for 28% of

235

global pasture area. For birds, in addition to the above countries, New Zealand also showed high

236

species loss. For amphibians, Colombia, Brazil, Ecuador, Peru and Venezuela were identified as

237

hotspots of species loss due to pasture land use.

238

Regarding the managed forests, Russia, Brazil, Canada, USA and China are the top five

239

countries in terms of area. However they do not feature in the top 5 countries with highest

240

species loss. For instance, despite the fact that Russia has highest amount of forest land use, it

241

ranks 21, 87, 101 and 73rd on the mammal, birds, amphibians and reptile species lost

242

respectively. This is due to the fact that these countries have low species richness to start with

243

and the low vulnerability score of taxa hosted by them, thereby having low characterization

244

factors (equation 1).21Indonesia, Papua New Guinea, Madagascar, DR Congo, Brazil and

245

Malaysia suffer the highest species loss due to anthropogenic forest land use (Table S1, SI-2).

246

Table S2 in SI-2 shows the total biodiversity loss associated with each of the 160 crops. Rice,

247

maize and wheat were expected to contribute most to biodiversity loss because together these

248

three crops occupy ~40% of global agricultural land. However, crops such as coffee, rubber, tea,

249

palm oil and soybean have a disproportionally high biodiversity footprint considering the fact

250

that they only occupy less than 10% of global agricultural land. This is because these crops

ACS Paragon Plus Environment

Page 13 of 27

Environmental Science & Technology

251

occupy biodiversity-rich regions hosting high number of endemic and threatened species (i.e.

252

with high CFs, equation 1).

253

Table S3 in SI-2 throws further insight into the biodiversity burden of global agricultural land

254

use by listing the impacts for each of the 7,666 crop × country combinations. While wheat from

255

the USA, Canada, and Russia occupy large agricultural areas on earth, their contribution to

256

global land-use related biodiversity loss is meagre. Land use for rice, coconut, rubber and palm

257

oil production in South-east Asian countries Indonesia, Malaysia and Philippines were found to

258

contribute the most to biodiversity loss among all combinations of crop and country. For

259

amphibians and reptiles, maize production in Mexico and tea in Sri Lanka also contributed

260

significantly to species loss. For each taxa, just the top 30 combinations of crop × country are

261

responsible for ~50% of global agriculture impact.

262

Country-specific characterization factors (CFs)

263

The calculated CFs for pasture, extensive and intensively managed forest land use for each

264

country (from equation 4) are listed in SI-2, Table S4. Pasture land use had in general higher

265

CFs, reflecting the relatively low affinity of species to them as compared to natural managed

266

forests but were close to intensively managed forest. All three CFs were within one order of

267

magnitude for most countries. The CFs were in general higher for tropical countries and small,

268

island countries (such as in Caribbean); lower for countries in temperate and boreal regions and

269

varied over 4 orders of magnitude across 250 countries. For agricultural land use, CFs were

270

calculated for 7,666 crop and country combinations and are shown in SI-2, Table S5.

271

Impacts per m3 of roundwood production

272

Figure 2 and Table S6 in SI-2 show that for an equal amount of roundwood production, the

273

biodiversity impacts range over 4 orders of magnitude across different countries. Island and

ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 27

274

tropical countries such as Madagascar, Comoros, Sao Tome and Dominican Republic suffer

275

highest species lost per m3 of roundwood produced. However most of these countries have little

276

forest area and thus low wood production (Table S6, SI-2). Countries producing majority of

277

global roundwood such as the US, Canada, Russia, Brazil -rank very low in terms of impact per

278

m3. The impacts differ because of interplay of factors such as harvesting intensity, CFs and forest

279

management regime prevalent in different countries. For example, the forestry CFs for mammal

280

species loss (equation 4) for Indonesia is around 500 times higher than that for Germany (both

281

for plantations and natural managed forests, Table S6). Further, the average wood yield from

282

forests in Germany is more than four times higher than in Indonesia. Consequently, mammal

283

species lost per m3 of roundwood is ~2000 times higher in Indonesia (4.7 ×10-7) than in Germany

284

(2.5 ×10-10). In general, countries with low yields per hectare along with high biodiversity and

285

presence of endemic and threatened species (high global CFs) performed the worst in terms of

286

impacts per m3.

287 288

Figure 2. Total mammal impacts per m3 of roundwood produced in 139 countries (unit- global

289

species eq. lost·years/m3) calculated using equation 5 above. See Table S6 in supporting

ACS Paragon Plus Environment

Page 15 of 27

Environmental Science & Technology

290

information-2 for impacts on other taxa per country. NA denotes the countries with negligible

291

roundwood production.

292

Biodiversity impacts per kg of crop

293

Table S5 in SI-2 shows the biodiversity impacts per kg of crop in each country. For all four taxa,

294

the impacts per kg for a particular crop ranged over five orders of magnitude (~10-11 to 10-16

295

species eq. lost·years/kg) depending on the country. For example, 1 kg of wheat from Brazil

296

results in 1.05 x 10-11 mammal species eq. lost·years as compared to 5.94 x 10-13 in Canada

297

(Table S5). The impacts per kg were generally higher for amphibians and reptiles followed by

298

mammals and least for birds but were mostly within an order of magnitude.

299

Case-study: Swiss food import impacts

300

Table 1 shows the top 10 crop × country combinations causing highest mammal impacts along

301

with their net mass and area imported (i.e. after subtracting the exports) into Switzerland for the

302

year 2011. Cocoa, coffee and soybean imports resulted in highest mammal species loss in

303

countries of origin.

304

Table 1. Top 10 crop items imported into Switzerland for the year 2011 with highest embodied

305

mammal impacts (in units – species eq. lost·years) along with their country of origin. Mass of

306

crop items (in tons) and corresponding area embodied (in ha) is also shown. Taxa-aggregated

307

impacts are in the units - potentially disappeared fraction·years (PDF·years). See Table S7 in

308

supporting information-2 for full list of imported crop items and their embodied impacts. Crop

cocoa

From

Ecuador

Mass

Land

Mammal

Birds

Amphibians

Reptiles

Aggregated

imported

imported

(Species eq.

(Species eq.

(Species eq.

(Species eq.

(PDF·years)

(Tons)

(ha)

lost·years)

lost·years)

lost·years)

lost·years)

4

1.1×10

5.1×10

4

9.3×10

-3

5.7×10

-3

ACS Paragon Plus Environment

-3

74×10

4.6×10

-3

7.2×10

-6

Environmental Science & Technology

4

1.2×10

5

6.8E×10

3

6.9×10

3

2.7×10

4

4.3×10

4

2.50×10

3

6.0×10

3

2.50×10

4

2.4×10

4

2.30×10

3

2.4×10

3

1.90×10

3

9.2×10

3

1.70×10

4

2.2×10

4

1.40×10

5

3.2×10

4

1.30×10

cocoa

Ghana

3.5×10

cocoa

Cameroon

2.4×10

sunflower

Tanzania

2.6×10

coffee

India

5.3×10

coffee

Brazil

1.9×10

coconut

Philippines

9.6×10

coffee

Colombia

8.4×10

cocoa

Ivory coast

1.5×10

soybean

Brazil

2.5×10

-3

-3

1.1×10

-3

Page 16 of 27

5.6×10

-3

-3

2×10

-4

4.3×10

-4

1.40×10

-4

7.30×10

-4

1.50×10

-4

1.50×10

-4

1.30×10

-4

2.70×10

0.33×10

-3

7.90×10

-3

4.10×10

-3

7.60×10

-3

5.90×10

-3

8.40×10

-3

2.20×10

-3

4.80×10

-3

-3

3.2×10

-3

1.7×10

-6

0.2×10

-3

4.8×10

-7

-4

7.1×10

-7

-2

3.0×10

-6

-4

8.4×10

-7

-3

7.3×10

-7

-4

1.4×10

-6

-4

3.5×10

-7

-4

3.8×10

-7

5.30×10

-2

1.10×10

-3

2.40×10

-3

2.40×10

-2

8.90×10

-3

6.10×10

-3

1.80×10

309 310

Total Swiss crop imports resulted in loss of 0.0556 mammals, 0.019 birds, 0.189 amphibians and

311

0.046 reptile species eq. lost·years (Table S7, SI-2). The total imported impacts are 20-300 times

312

higher than those due to domestic crop land use for Swiss consumption in Switzerland for

313

different taxa considered (Table S7). This is much higher in contrast to the ratio of total land

314

embodied in imported products and net domestic agricultural land used for consumption (= 3.5).

315

This suggests embodied land is not a good proxy for embodied biodiversity impacts. Overall

316

>95% of biodiversity impacts of Swiss food consumption occur outside its border.

317

DISCUSSION

318

The study is first to present the spatially-explicit biodiversity impacts due to global agriculture,

319

pasture and forest land use at 5 arc minute grid cell level (Figure 1, Figure S1 to S3 in supporting

320

information-1) and aggregated to country level (Table S1, SI-2). As land use decisions are

321

mostly made at national and sub-national level, country-specific impacts and their distribution

ACS Paragon Plus Environment

Page 17 of 27

Environmental Science & Technology

322

within the nation for pasture, forestry land and specific crops can be used by producer nations to

323

identify geographical hotspots and most damaging land use types. This in turn can induce the

324

changes in production methods and other measures to control further damage (e.g. by shifting

325

from high to low intensity agriculture or forestry or by protecting ecologically valuable habitats).

326

On the consumption side, the study is first to provide ready to use factors quantifying

327

biodiversity impacts per m3 of roundwood (Table S5, SI-2) and per kg of crop (Table S6, SI-2)

328

from each country. The Swiss case study demonstrated how these factors can be linked with

329

trade data to identify commodities causing high biodiversity loss (e.g. imports of cocoa from

330

Ecuador, coffee from India or soybean from Brazil). Reducing the volume of imported trade

331

commodities that cause high species loss and raising consumers’ awareness about the

332

biodiversity damage caused by the products they buy can go a long way in reducing the existing

333

rate of biodiversity loss.

334

Overall the study marks a significant improvement over previous land use impact assessment

335

methods within LCA that provide characterization factors (CFs) at an ecoregion level for broad

336

land use categories.19,21 Combining the high-resolution yield and area maps of global crops26

337

with the ecoregion specific CFs for annual and permanent crops from Chaudhary et al.21, enabled

338

us to calculate the impacts and CFs for each of the 160 individual crops from 250 different

339

countries (Table S2, S3 in SI-2). Within LCA inventory databases, the land use flow information

340

is often available at country level rather than ecoregion, therefore the country-specific CFs are

341

more useful and can be directly used for impact assessment (Table S4, SI-2).

342

High biodiversity impacts associated with certain land use types in particular regions can be

343

traced back to high CFs for these regions (see all maps in SI-1). To calculate biodiversity

344

damage, we used the CFs from Chaudhary et al.21which were calculated by combining

ACS Paragon Plus Environment

Environmental Science & Technology

Page 18 of 27

345

countryside species area relationship (SAR) and taxa-specific vulnerability score (VS). The SAR

346

model includes aspects of ecosystem vulnerability (i.e. how much an ecosystem is already

347

affected by land use pressures), and the VS addresses the vulnerability of species inhabiting a

348

particular region to future land use pressures. The VS give particular weight to regions hosting

349

range-restricted and threatened species that are near extinction and whose loss can result in

350

permanent loss of unique evolutionary history associated with them.21 The CFs are higher for

351

ecoregions hosting biodiversity that is unique and endemic and is found nowhere else, and they

352

are lower for the ecoregions that contain only tiny fractions of species’ range (mostly range

353

edges). As the VS also incorporates IUCN assigned species threat level, it means that for two

354

regions hosting equal endemic richness, the CFs will be higher for regions containing more

355

threatened species than those containing non-threatened species.21 The CFs therefore help flag

356

regions with high land use impacts on species requiring immediate conservation attention.

357

Tropical and island countries were found to suffer the highest biodiversity damage due to land

358

use. For example, Ecuador ranks 61st globally in terms of total forest land use area. However, it

359

ranks 9th on the country suffering highest amphibian species loss due to forest land use (Table

360

S1, SI-2). This is because it is home to 336 amphibian endemic species that are not found

361

anywhere else in the world and thus any land use poses high extinction risk to them. Our results

362

show that for the biodiversity impacts, the region where the land use takes place is more

363

important than the total area occupied.

364

Regarding the calculated impacts per m3 of roundwood, the results showed that countries with

365

high forestry CFs and low wood harvesting yields fare the worst (Table S6, SI-2). Interestingly,

366

we found that some sub-tropical, species rich countries hosting many endemic and threatened

367

species (i.e. with high CF) performed very well in terms of biodiversity loss per m3 of

ACS Paragon Plus Environment

Page 19 of 27

Environmental Science & Technology

368

roundwood production. For example, in India, the majority of roundwood is sourced from

369

planted forests (94%) which have 50 times higher yield (3.8 m3/ha/year) than corresponding

370

natural managed forests (0.07 m3/ha/year), leading to low impacts per m3. The case-study results

371

showed that food imports to Switzerland are responsible for ~20-300 times the biodiversity loss

372

occurring domestically due to its total agricultural land use (Table S7, SI-2). This ratio actually

373

might be even higher as we did not include finished livestock products imported from abroad and

374

just considered crop products in our analysis. Overall the results corroborate with unequal-

375

exchange theory37 and previous studies who also found that food consumption in industrialized

376

nations drives biodiversity loss in tropical developing countries through international trade.9

377

Limitations and data gaps

378

The input data used to calculate the characterization factors through SAR model as well as the

379

global crop, pasture and forestry land use maps used come with uncertainties and limitations that

380

should be considered when interpreting the results. For example, the species affinity estimates

381

(ℎ,,+ , equation S1 of SI-1) fed into the SAR model were derived from empirical data from

382

literature review (see supporting information-1).21 As more plot-scale local biodiversity

383

monitoring data comes along (e.g. PREDICTS database)38, these estimates can be updated and

384

accuracy of results can be improved. Further, the yield and harvested area maps used to derive

385

impacts per kg (equation 6) were based on the year 2000.26 Yields and area of some crops might

386

have increased (or decreased) over last decade in different world regions.39 Thus we might have

387

over or under estimated some of the impacts but currently these maps are the best available

388

source for high resolution inventory of global crop land use.

389

We used species richness loss as an indicator of land use driven biodiversity loss. However, this

390

indicator masks the biodiversity damage due to changes in species composition that can take

ACS Paragon Plus Environment

Environmental Science & Technology

Page 20 of 27

391

place following disturbances. Future studies should explore the alternative measures of

392

biodiversity loss such as beta diversity index 40 or mean species abundance41 which might reveal

393

alternate hot-spots of land use impacts. Another limitation of the analysis is that we only

394

considered four taxa and owing to lack of data, impacts on other species groups such as

395

arthropods, bacteria, fungi that make up majority of global terrestrial species richness could not

396

be calculated. These species groups perform several important ecological functions and the data

397

gaps should be filled through future research efforts.

398

For roundwood production, we took average yield of plantation forests and managed natural

399

forest per country. However, biodiversity damage and yields differ greatly depending upon the

400

harvesting techniques used in the region. In this study all managed natural forests were grouped

401

into broad category of ‘extensive forestry.’ In reality, low yields observed in some of the

402

countries might be because of low intensity harvesting techniques rather than technological

403

limitations and therefore our calculated impacts might be overestimated for these countries.

404

Many studies have shown that forests managed using low intensity harvesting techniques such as

405

reduced impact logging42 or retention harvesting43 result in negligible species loss as compared

406

to conventional selective logging or clear-cut regimes. However, spatially-explicit maps

407

depicting forest management regimes and their harvesting intensity are currently unavailable on

408

a global scale.

409

Finally, any application of the methods and results presented should be accompanied by an

410

analysis of other impact categories. In this paper, we calculate biodiversity impacts due to land

411

use which is the main global driver of biodiversity loss. However, methods compatible with our

412

approach are still missing and more effort is needed to harmonize the methods and metrics so

413

that the biodiversity impacts from other stressors can be jointly assessed. We used global

ACS Paragon Plus Environment

Page 21 of 27

Environmental Science & Technology

414

(permanent) species extinctions as a metric for biodiversity loss, but other impact categories

415

within LCA typically quantify impacts as local or regional species loss.13 We assess the impacts

416

for four vertebrate taxa while the LCA methods for other stressors such as acidification44

417

consider impacts on vascular plants, while eutrophication impacts are usually assessed on fish

418

species45 due to data availability restrictions. To allow for a comparison of biodiversity impacts

419

from different impact categories, consistent approaches for harmonizing the metrics and

420

aggregating the impacts across different taxa are needed.13,46 Efforts are underway in this

421

direction as Verones et al.46 proposed approaches for harmonized assessment of biodiversity loss

422

from land and water use, but future studies should aim to include additional stressors. Similarly,

423

other impacts of land use such as soil erosion47 or appropriation of net primary productivity48

424

must also be taken into account in environmental decision making. For the Swiss case-study,

425

including these impacts will further increase the imported biodiversity impacts.

426

Outlook

427

The overall results are useful for producer nations and can be a first step for further in-depth

428

investigations aimed at sustainable land use management on a global and regional level, as they

429

reveal current geographical hotspots and the drivers of biodiversity loss. Results are also relevant

430

to consumers, global retailers and food processing companies who are increasingly interested in

431

the environmental product information. By quantifying biodiversity impacts per kg of crop from

432

different locations, the results can help improve the life-cycle based product information, which

433

currently often only address carbon emission impacts (e.g. UK carbon reduction label).9 The

434

calculated country, crop and taxa–specific impacts can also be used as a basis for compensatory

435

mechanisms or offsetting programs. For example, food products with high biodiversity impacts

436

could be made more expensive and the premium can go towards financing ecosystem service or

ACS Paragon Plus Environment

Environmental Science & Technology

Page 22 of 27

437

biodiversity conservation programs. Such biodiversity compensation programs could

438

complement already existent efforts to reduce and compensate for greenhouse gas emissions.9 To

439

ensure progress toward reducing the rate of global biodiversity loss, policies aimed at all the

440

actors in the supply chain (be it exporters, traders or consumers) would have to be implemented

441

in parallel. 9

442

The calculated impacts per m3 of roundwood production demonstrated how species rich regions

443

with currently high forest land use impacts can potentially benefit by producing wood from high-

444

yield planted forests. Planted forests are less biodiversity-benign than extensively managed

445

natural forests6, but, if grown on previously degraded lands and inducing no land-use change

446

from natural ecosystems, they can serve dual purpose of timber production and alleviating the

447

pressure on the remaining natural forests which in turn can spared for biodiversity conservation.

448

While several trade databases9,33,35,49 documenting the flow of commodities between countries

449

exist, the factors giving biodiversity impacts per unit commodity have not been available. Our

450

estimates of biodiversity impacts per kg of crop or per m3 of roundwood fill these gaps. The

451

Swiss case study showed how the results can be combined with trade data to identify the location

452

and severity of environmental impacts caused by imported goods in a country. Future studies

453

should expand this analysis to quantify biodiversity impacts embodied in global trade of food

454

and forestry products.

455

Supporting Information

456

The supporting information-1 contains additional methods and maps. Supporting information-2

457

excel file contains Tables S1 to S7. This information is available free of charge via the Internet at

458

http://pubs.acs.org.

459

Corresponding Author

ACS Paragon Plus Environment

Page 23 of 27

Environmental Science & Technology

460

* E-mail: [email protected]; phone: +41-44-633-02-54; fax: +41-44-633-10-61

461

ACKNOWLEDGEMENTS

462

This research was funded within the National Research Programme “Resource Wood” (NRP 66)

463

by the Swiss National Science Foundation (project no. 136612).

464

REFERENCES

465

1

Ceballos, G.; Ehrlich, P. R.; Barnosky, A. D.; García, A.; Pringle, R. M.; Palmer, T. M.

466

Accelerated

467

extinction. Science advances, 2015, 1(5), e1400253.

468

2

modern

human–induced

species

losses:

Entering

the

sixth

mass

Millennnium Ecosystem Assessment Ecosystems and human wellbeing: biodiversity synthesis; Washington, DC, 2005.

469 470

3

Foley, J.A., et al. Global consequences of land use. Science 2005, 309,570–574.

471

4

Sanderson, E.W., Jaiteh, M., Levy, M.A., Redford, K.H., Wannebo, A.V., & Woolmer, G. The human footprint and the last of the wild. Bioscience 2002, 52, 891–904.

472 473

5

Food and Agriculture Organization, Rome, 2012.

474 475

Alexandratos, N.; Bruinsma, J. World Agriculture Towards 2030/2050: The 2012 Revision,

6

Gibson, L.; Lee, T. M.; Koh, L. P.; Brook, B. W.; Gardner, T. A.; Barlow, J.; ... Sodhi, N.

476

S. Primary forests are irreplaceable for sustaining tropical biodiversity. Nature 2011,

477

478(7369), 378-381.

478

7

Tittensor, D. P.; Walpole, M.; Hill, S. L.; Boyce, D. G.; Britten, G. L.; Burgess, N. D.; ...

479

Visconti, P. A mid-term analysis of progress toward international biodiversity targets.

480

Science 2014, 346(6206), 241-244.

481

8

hotspots for conservation priorities. Nature, 2000, 403(6772), 853-858.

482 483

9

Lenzen, M.; Moran, D.; Kanemoto, K.; Foran, B.; Lobefaro, L.; Geschke, A. International trade drives biodiversity threats in developing nations. Nature 2012 486(7401), 109-112.

484 485

Myers, N.; Mittermeier, R. A.; Mittermeier, C. G.; Da Fonseca, G. A.; Kent, J. Biodiversity

10

Erb, K. H.; Krausmann, F.; Lucht, W.; Haberl, H. Embodied HANPP: Mapping the spatial

486

disconnect between global biomass production and consumption. Ecol. Econ. 2009, 69(2),

487

328-334.

ACS Paragon Plus Environment

Environmental Science & Technology

488

11

cycle assessment. Science 2014, 344(6188), 1109-1113.

489 490

Hellweg, S.; i Canals, L. M. Emerging approaches, challenges and opportunities in life

12

de Baan, L.; Curran, M.; Rondinini, C.; Visconti, P.; Hellweg, S.; Koellner, T. High-

491

resolution assessment of land use impacts on biodiversity in life cycle assessment using

492

species habitat suitability models. Environ. Sci. Technol., 2015, 49(4), 2237-2244.

493

13

Curran, M.; de Baan, L.; de Schryver, A.M.; van Zelm, R.; Hellweg, S.; Koellner, T.;

494

Sonnemann, G.; Huijbregts, M.A.J. Toward meaningful endpoints of biodiversity in Life

495

Cycle Assessment. Environ. Sci. Technol., 2011, 45 (1), 70–79.

496

14

Koellner, T., de Baan, L., Beck, T., Brandão, M., Civit, B., Margni, M., ... & Müller-Wenk,

497

R.. UNEP-SETAC guideline on global land use impact assessment on biodiversity and

498

ecosystem services in LCA. Int. J. Life Cycle Assess. 2013, 18 (6), 1188−1202.

499

15

with Life Cycle Assessment: are we there yet?. Glob. Chang. Biol. 2015, 21(1), 32-47.

500 501

Souza, D. M.; Teixeira, R. F.; Ostermann, O. P. Assessing biodiversity loss due to land use

16

Goedkoop, M.; Heijungs, R.; Huijbregts, M.; De Schryver, A.; Struijs, J.; van Zelm, R.

502

ReCiPe Main Report, Part 1: Characterization, 1st ed.; Ministerie van Volkshuisvesting,

503

Ruimtelijke Ordening en Milieubeheer (VROM): The Hague, 2009; http://www.pre-

504

sustainability.com/download/misc/ReCiPe_main_report_final_27-02 2009_web.pdf.

505

17

approach. Int. J. Life Cycle Assess. 2013a, 18 (6), 1216−1230.

506 507

de Baan, L.; Alkemade, R.; Koellner, T. Land use impacts on biodiversity in LCA: a global

18

Olson, D.; Dinerstein, E.; Wikramanayake, E.; Burgess, N.; Powell, G.; Underwood, E.;

508

D’Amico, J.; Itoua, I.; Strand, H.; Morrison, J.; Loucks, C.; Allnutt, T.; Ricketts, T.; Kura,

509

Y.; Lamoreux, J.; Wettengel, W.; Hedao, P.; Kassem, K. Terrestrial ecoregions of the

510

worlds: A new map of life on Earth. BioScience 2001, 51 (11), 933− 938.

511

19

de Baan, L.; Mutel, C. L.; Curran, M.; Hellweg, S.; Koellner, T. Land use in life cycle

512

assessment: global characterization factors based on regional and global potential species

513

extinction. Environ. Sci. Technol. 2013b, 47 (16), 9281-9290.

514

20

Pereira, H. M.; Ziv, G.; Miranda, M. Countryside Species–Area Relationship as a Valid

515

Alternative to the Matrix‐Calibrated Species–Area Model. Conserv Biol. 2014, 28 (3), 874-

516

876.

ACS Paragon Plus Environment

Page 24 of 27

Page 25 of 27

517

Environmental Science & Technology

21

Chaudhary, A.; Verones, F.; de Baan, L.; Hellweg, S. Quantifying Land Use Impacts on

518

Biodiversity: Combining Species–Area Models and Vulnerability Indicators. Environ. Sci.

519

Technol., 2015, 49(16), 9987-9995.

520

22

case studies of three vegetable oil crops. J. Cleaner Prod. 2000, 8(4), 283-292.

521 522

23

24

Eshun, J. F., Potting, J., Leemans, R. LCA of the timber sector in Ghana: preliminary life cycle impact assessment (LCIA). Int. J. Life Cycle Assess. 2011, 16(7), 625-638.

525 526

Michelsen, O. Assessment of land use impact on biodiversity. Int. J. Life Cycle Assess. 2008, 13(1), 22-31.

523 524

Mattsson, B.; Cederberg, C.; Blix, L. Agricultural land use in life cycle assessment (LCA):

25

Mutel, C. L.; Pfister, S.; Hellweg, S. GIS-based regionalized life cycle assessment: how

527

big is small enough? Methodology and case study of electricity generation. Environ. Sci.

528

Technol. 2011, 46(2), 1096-1103.

529

26

Monfreda, C., Ramankutty, N., Foley, J. A. Farming the planet: 2. Geographic distribution

530

of crop areas, yields, physiological types, and net primary production in the year 2000.

531

Glob. Biogeochem. Cycles, 2008, 22(1), GB1022.

532

27

LADA, Mapping Land use Systems at global and regional scales for Land Degradataion

533

Assessment Analysis. Nachtergaele, F., Petri, M. In LADA Technical Report n.8, version

534

1.1; UNEP/GEF; 2008.

535

28

Frontiers Ecol. Environ. 2008, 6 (8), 439−447.

536 537

Ellis, E.; Ramankutty, N. Putting people in the map: anthropogenic biomes of the world.

29

Hansen, M. C., Potapov, P. V., Moore, R., Hancher, M., Turubanova, S. A., Tyukavina, A.,

538

... & Townshend, J. R. G. High-resolution global maps of 21st-century forest cover

539

change. Science 2013, 342(6160), 850-853.

540

30

Pfister, S., Bayer, P., Koehler, A., & Hellweg, S. Environmental impacts of water use in

541

global crop production: hotspots and trade-offs with land use. Environ. Sci. Technol. 2011,

542

45(13), 5761-5768.

543

31

assessment 2015: Main report. Food and Agriculture Organization of the United Nations.

544 545

Food and Agriculture Organization of the United Nations. Global forest resources

32

Jürgensen, C., Kollert, W. and Lebedys, A. Assessment of industrial roundwood

546

production from planted forests. FAO Planted Forests and Trees Working Paper FP/48/E.

547

Rome. 2014.

ACS Paragon Plus Environment

Environmental Science & Technology

548

33

Page 26 of 27

FAOSTAT, 2015. FAO Statistical Databases: Agriculture, Fisheries, Forestry, Nutrition.

549

Statistics Division, Food and Agriculture Organization of the United Nations, Rome

550

(Available at http://faostat.fao.org/).

551

34

Organization of the United Nations, Rome. 2003.

552 553

35

36

Kastner, T., Kastner, M., & Nonhebel, S. Tracing distant environmental impacts of agricultural products from a consumer perspective. Ecol. Econ. 2011, 70(6), 1032-1040.

556 557

Wood, R.., et al. Global sustainability accounting-developing EXIOBASE for multiregional footprint analysis. Sustainability 2015, 7 (1), pp. 138-163.

554 555

FAO. Technical Conversion Factors for Agricultural Commodities. Food and Agriculture

37

Shandra, J. M., Leckband, C., McKinney, L. A., & London, B. Ecologically Unequal

558

Exchange, World Polity, and Biodiversity Loss A Cross-National Analysis of Threatened

559

Mammals. Int. J. Comp. Sociol. 2009, 50(3-4), 285–310.

560

38

Hudson, L. N., Newbold, T., Contu, S., Hill, S. L., Lysenko, I., De Palma, A., ... &

561

Choimes, A. The PREDICTS database: a global database of how local terrestrial

562

biodiversity responds to human impacts. Ecol. Evol. 2014, 4 (24), 4701-4735.

563

39

Mueller, N. D.; Gerber, J. S.; Johnston, M.; Ray, D. K.; Ramankutty, N.; Foley, J. A.

564

Closing yield gaps through nutrient and water management. Nature, 2012, 490(7419), 254-

565

257.

566

40

Fattorini, S.; R.L.H. Dennis.; Cook, L.M. Conserving organisms over large regions

567

requires multi-taxa indicators: one taxon’s diversity-vacant area is another taxon’s

568

diversity zone. Biological Conservation, 2011, 144, 1690-1701.

569

41

Alkemade, R.; van Oorschot, M.; Miles, L.; Nellemann, C.; Bakkenes, M.; ten Brink, B.

570

GLOBIO3: A framework to investigate options for reducing global terrestrial biodiversity

571

loss. Ecosystems 2009, 12 (3), 374−390.

572

42

opportunities. Forest Ecol. Manag. 2008, 256, 1427–1433.

573 574

Putz, F. E.; Sist, P.; Fredericksen, T.; Dykstra, D. Reduced-impact logging: Challenges and

43

Fedrowitz, K., Koricheva, J., Baker, S. C., Lindenmayer, D. B., Palik, B., Rosenvald, R., ...

575

& Gustafsson, L. Can retention forestry help conserve biodiversity? A meta‐analysis. J.

576

Appl. Ecol. 2014, 51(6), 1669-1679.

ACS Paragon Plus Environment

Page 27 of 27

577

Environmental Science & Technology

44

Azevedo, L. B.; van Zelm, R.; Hendriks, A. J.; Bobbink, R.; Huijbregts, M. A. J. Global

578

assessment of the effects of terrestrial acidification on plant species richness. Environ.

579

Pollut. 2013a, 174, 10-15.

580

45

Azevedo, L. B.; Henderson, A. D. ; van Zelm, R. ; Jolliet, O.; Huijbregts, M. A. J.

581

"Assessing the Importance of Spatial Variability versus Model Choices in Life Cycle

582

Impact Assessment: The Case of Freshwater Eutrophication in Europe." Environ. Sci.

583

Technol. 2013b, 47(23), 13565-13570.

584

46

Verones, F. ; Huijbregts, M. A. ; Chaudhary, A. ; de Baan, L. ; Koellner, T.; Hellweg, S.

585

Harmonizing the assessment of biodiversity effects from land and water use within

586

LCA. Environ. Sci. Technol. 2015, 49(6), 3584-3592.

587

47

Núñez, M. ; Antón, A. ; Muñoz, P.; Rieradevall, J. Inclusion of soil erosion impacts in life

588

cycle assessment on a global scale: application to energy crops in Spain. Int. J. Life Cycle

589

Assess. 2013, 18(4), 755-767.

590

48

Haberl, H. ; Erb, K. H. ; Krausmann, F. ; Gaube, V., Bondeau, A., Plutzar, C., ... &

591

Fischer-Kowalski, M. Quantifying and mapping the human appropriation of net primary

592

production in earth's terrestrial ecosystems. Proc. Natl. Acad. Sci. U.S.A. 2007, 104(31),

593

12942-12947.

594 595

49

Moran, D.; Petersone, M.; Verones, F. On the suitability of input–output analysis for calculating product-specific biodiversity footprints. Ecol. Indic. 2016, 60, 192-201.

596 597

TOC art:

598 599

ACS Paragon Plus Environment