Land Use Intensity-specific Global ... - ACS Publications

types without differentiating between their management intensities, and have large. 13 uncertainties that limit .... Next, we utilize the land use int...
3 downloads 11 Views 780KB Size
Subscriber access provided by University of Groningen

Characterization of Natural and Affected Environments

Land Use Intensity-specific Global Characterization Factors to Assess Product Biodiversity Footprints Abhishek Chaudhary, and Thomas M. Brooks Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b05570 • Publication Date (Web): 12 Apr 2018 Downloaded from http://pubs.acs.org on April 13, 2018

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 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 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.

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

Land Use Intensity-specific Global Characterization

2

Factors to Assess Product Biodiversity Footprints

3

Abhishek Chaudhary1*, Thomas M. Brooks2

4 5 6

1

Institute of Food, Nutrition and Health, ETH Zurich, 8092 Zurich, Switzerland, (*Corresponding author phone: +41 76 757 7400; fax: +41 44 632 1155; e-mail: [email protected])

2

International Union for Conservation of Nature, 28 Rue Mauverney, 1196 Gland, Switzerland

7 8

ABSTRACT

9

The UNEP-SETAC life cycle initiative recently recommended use of the countryside species-

10

area relationship (SAR) model to calculate the characterization factors (CFs; potential species

11

loss per m2) for projecting the biodiversity impact of land use associated with a products’ life

12

cycle. However, CFs based on this approach are to date available for only six broad land use

13

types without differentiating between their management intensities, and have large

14

uncertainties that limit their practical applicability. Here we derive updated CFs for

15

projecting potential species losses of five taxa resulting from five broad land use types

16

(managed forests, plantations, pasture, cropland, urban) under three intensity levels (minimal,

17

light and intense use) in each of the 804 terrestrial ecoregions. We utilize recent global land

18

use intensity maps and International Union for Conservation of Nature (IUCN) habitat

19

classification scheme to parametrize the SAR model. As a case study, we compare the

20

biodiversity impacts of 1m3 of wood produced under four different forest management

21

regimes in India and demonstrate that the new land use intensity-specific CFs have smaller

22

uncertainty intervals and are able to discern the impacts of intensively managed land uses

23

from the low intensity regimes, which has not been possible through previous CFs.

24 25 26

1 ACS Paragon Plus Environment

Environmental Science & Technology

27

Table of content (TOC) art:

28 29 30

INTRODUCTION

31

Environmentally conscious purchasing behavior by individual consumers, businesses and

32

nations can contribute substantially towards mitigating the ongoing biodiversity crisis.1

33

Informing consumers about the biodiversity impacts ‘hidden’ in the life cycle of products

34

they consume is an important step towards transforming current consumption patterns.2-8 Life

35

cycle assessment (LCA) is increasingly used for estimating and comparing the ‘cradle to

36

grave’ environmental impacts of products and commodities.9 In LCA, biodiversity impacts of

37

land use are calculated by multiplying inventory (i.e. land occupation in m2∙year and land

38

transformation in m2; ) with characterization factors (CFs, i.e. the factors indicating

39

biodiversity damage caused by unit area of a particular land use).4

40

A suite of CFs have been proposed to assess the species loss within a product’s life cycle,

41

using different ecological models, and varying taxonomic coverage and spatial resolution.5, 9,

42

10

43

employed. A lack of consensus on which CFs best reflect the damage inflicted by a given

44

product system has, thus, hampered the application of these methods to inform decision

45

makers and stakeholders.11

46

To address this challenge, the land use task force of the United Nations Environment

47

Programme/Society of Environmental Toxicology and Chemistry (UNEP-SETAC) Life

48

Cycle Initiative was established to provide guidance and build scientific consensus on CFs

Different results are obtained when comparing the products depending upon which CFs are

2 ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32

Environmental Science & Technology

49

for assessing land use impacts on biodiversity.12 Over 2015-2017, the initiative organized

50

discussions between experts and stakeholders13 as well as conducted a critical review of the

51

existing framework and biodiversity models.5 It concluded that an acceptable model should

52

incorporate both the intensity of human land use as well as the vulnerability of the species

53

group in the region where the land use is taking place (e.g. if the region in question holds

54

many endemic or threatened species). The initiative culminated in a Pellston WorkshopTM,

55

which provisionally recommended the CFs calculated by Chaudhary et al.4 to assess land use

56

driven biodiversity impacts within LCA (see UNEP report).14

57

The CFs calculated by Chaudhary et al.4 represent biodiversity damage resulting from global

58

land use (in potential species lost per m2) for four taxa: mammals, birds, amphibians, and

59

reptiles. The CFs assess the potential damage caused by a unit area of six different land use

60

types (annual crops, permanent crops, pasture, intensive forestry, extensive forestry and

61

urban land) in each of the 804 terrestrial ecoregions.15 The CFs were calculated by combining

62

countryside species-area relationship (SAR) model16 with a vulnerability score of the species

63

group.4

64

However, the UNEP-SETAC also recognized that a number of aspects of these CFs need

65

improvement.14 As such, it recommended that until these improvements are made, the

66

method should only be used for identification of potential impact hotspots but not for product

67

comparison and labelling purposes. The recommended improvements were:

68

1. First and most importantly, expand land use classes and include different management

69

regimes. In other words, the CFs should be able to differentiate between biodiversity

70

damage caused by different land use classes (e.g. pasture vs. cropland) and also between

71

products coming from same land use class but varying in intensity of use (e.g.

72

extensive/intensive pasture).

3 ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 32

73

2. CFs for plants should be calculated and included in the impact assessment.

74

3. The uncertainty in the CFs, which is currently very high, should be reduced.

75

4. The model predicted species loss should be compared against the global rate of extinction

76 77 78

and/or the observed species loss within the studied systems. 5. Case studies should be conducted to examine the robustness and ability of the model to rank alternative product systems in terms of their biodiversity impacts.

79

In this study, we present first attempts towards addressing the above five issues, and calculate

80

updated CFs for five broad land use types and their three intensity levels.

81

MATERIALS AND METHODS

82

Countryside SAR model

83

We calculate the projected species loss for each species group (g = 1:5; mammals, birds,

84

amphibians, reptiles and plants) in each of the 804 terrestrial ecoregions (j = 1:804) due to

85

current land use through countryside SAR model:16

,,



 , + ∑  ℎ,, ∙ , =  ,, ∙ 1 −    ,

86 87 88 89 90 91 92 93

Equation 1

Here  ,, is the total number of species occurring in each ecoregion’s area ( , ) before any human intervention, , is the natural habitat area in the ecoregion currently (in m2),

, is the current area of land use type  ( =1:16; see supplementary Table S1) in m2,  is the SAR exponent for the ecoregion and ℎ,, is the affinity of the taxon g to the land use type  in ecoregion j. As in Chaudhary et al.4 we obtained the species richness of four

vertebrate taxa per ecoregion ( ,, ) from WWF Wildfinder database17; species richness of

plants from Kier et al.18 and the z-values per ecoregion ( ) from Drakare et al.19 The

4 ACS Paragon Plus Environment

Page 5 of 32

Environmental Science & Technology

94

procedure to calculate the land use areas per ecoregion and the taxon affinity to human land

95

uses differing in terms of intensity of use is described in detail below.

96

Land use areas per ecoregion

97

We follow a two-step approach to calculate the areas of different land use types and their

98

intensity levels in each ecoregion.

99

Areas of broad land use types per ecoregion

100

We first derive the total area ( , ), area of remaining natural habitat ( , ) and the areas

101

of four broad human land use types (secondary vegetation, pasture, cropland and urban) per

102

ecoregion (, !" ) by overlaying the ecoregion boundaries with recently available global

103

‘land use’ map of Hoskins et al.20 for the year 2005. We chose the Hoskins et al.20 map for

104

several reasons. First, unlike LADA21 and ANTHROME22 maps used by Chaudhary et al.4, it

105

is of approximately 10 times higher resolution (30’’). Next, it has been validated using

106

independent global datasets23,24 and is also compatible with global land use harmonization

107

dataset.25 Due to above reasons, it is expected to be more accurate than the former. Hoskins et

108

al.20 has also been the first choice in several recent global biodiversity impact assessment

109

studies.26-29Although more recent and high resolution maps of individual land uses such as

110

cropland30 or forests31 are available in isolation, to date, Hoskins et al.20 remains the only

111

global high resolution, harmonized and validated map providing areas of primary habitat and

112

four broad human land use types as a one consistent dataset.

113

Areas of broad land use types under three intensity levels

114

Next, we utilize the land use intensity information from Global Land Systems dataset (year

115

2000)32 and Arets et al.33 in order to derive the proportion (#,

116

(secondary habitat, pasture, cropland and urban land use) under three different intensity levels

117

(minimal, light and intense use) per ecoregion (Table S1). Here we assume that land use

$

%$5 ACS Paragon Plus Environment

) of broad land use types

Environmental Science & Technology

Page 6 of 32

118

intensity did not change significantly between the year 2000 and 2005 (year of Hoskins et

119

al.20 map). The area of a particular broad land use type under a particular intensity level in an

120

ecoregion j ( , ) is then simply calculated as: $

%$ , = , !" × #, 121

Equation 2

122

As an illustration, consider the ecoregion Biak-Numfoor Rain Forests [WWF ecoregion code:

123

AA0103] on the island of New Guinea.17 According to Hoskins et al.20 map, the total area of

124

this ecoregion ( , ) is 2793 km2, the remaining natural habitat area (, ) is 2047 km2

125

!" and the area of cropland use (' (!", ) is 168 km2. Now according to Global Land

126

Systems map32, the proportion of cropland (#' (!", ) under minimal, light and intense use

127

in this ecoregion is 0.14, 0.58 and 0.28 respectively. Therefore, by Eq. 2 above, the area of

128

cropland under the minimal use is 168 × 0.14 = 23 km2. Similarly, area of light and intense

129

use cropland is equal to 168 × 0.58 = 98 km2 and 168 × 0.28 = 47 km2 respectively.

130

We consider ‘minimal use secondary vegetation’ to have been abandoned and therefore to be

131

currently regenerating and not used for production (following definition of Newbold et al.26).

132

The ‘light’ and ‘intense’ use levels of secondary vegetation in the forest biomes are

133

considered as ‘managed forests’ class, which is further divided into four intensity classes:

134

clear-cut, selective logging, plantations and reduced impact logging based on their fraction in

135

different world regions as reported by Arets et al.33 that in turn are based on Integrated Model

136

to Assess the Global Environment (IMAGE) output.34 Similar to Newbold et al.26, we

137

assumed that the plantation area per ecoregion was composed of equal proportions under

138

minimal, light and intense use. This is because none of the Global Land Systems32 classes

139

could inform about the intensity of plantation forestry. We divided the secondary forest

140

vegetation further into four different well-known classes (plantations, clear-cut, selective

$

%$6 ACS Paragon Plus Environment

Page 7 of 32

Environmental Science & Technology

141

logging and reduced impact logging) because these are the major wood production systems

142

across the globe.35 Separate characterization factors for these systems are therefore needed to

143

differentiate the impacts caused by them.

144

Finally, this yielded the areas of 16 land classes, i.e. ‘regenerating secondary vegetation’ plus

145

areas of five broad land use types (managed forests, plantations, pasture, cropland and urban)

146

under three intensity levels each (minimal, light and intense use) (see Table S1 in

147

Supplementary Information for more details on land use classification).

148

Taxon affinity calculation

149

As a proxy for taxon affinity (0 < ℎ,, < 1) to different human land uses, past studies

150

(including Chaudhary et al.4) have relied on published literature documenting the ratio of

151

observed species richness in a particular land use type to that in a natural undisturbed area,

152

obtained from few plot-scale field monitoring within the region.4,16,36 However, there are

153

large gaps in our knowledge of how damaging a particular land use type is to different taxa in

154

different parts of the world. The uncertainty in ℎ estimates is therefore substantial and

155

although global field monitoring efforts are underway, these gaps are unlikely to be filled in

156

near future.4 For example, Chaudhary et al.4 derived the affinity of five taxa to six land use

157

types in 14 biomes15 through a meta-analysis of 861 field monitoring data points comparing

158

species richness in human land use and nearby natural (primary) vegetation. For most taxa,

159

land use and ecoregion combinations, there were no field study available and therefore a

160

biome-average value was taken as a proxy. Furthermore, the data represented only a fraction

161

of the total number of species of a taxa present globally (e.g. the affinity values for mammals,

162

birds and amphibians were just based on 66, 211 and 11 data points globally).4 This rendered

163

the affinity values with large uncertainty intervals.

7 ACS Paragon Plus Environment

Environmental Science & Technology

164

Here, we therefore apply a novel hybrid two-step approach to calculate the affinity of taxa

165

(for which the data is available, i.e. mammals, birds and amphibians) to different land use

166

types and their intensity levels.

167

Taxon affinity to broad land use types

168

We first calculate the affinity of taxa to five broad land use types (managed forests,

169

plantations, pasture, cropland and urban) in each ecoregion through the recently proposed

170

approach by Chaudhary & Brooks28, i.e. using the species-specific information on habitat

171

utilization documented through the IUCN Red List Habitat Classification Scheme.37 For each

172

of the 22,386 species of these three taxa, the IUCN Habitat Classification Scheme provides

173

information on the human land use types which a particular species utilises. This information

174

is collected on each species through field observation and literature synthesis by the Red List

175

assessors of the species in question (see Section 2.4.1, page 24, of IUCN’s documentation

176

standards and consistency checks for species accounts).38

177

For example, the House Sparrow (Passer domesticus) utilises all broad human land use types.

178

The frog Adenomera martinezi utilises agriculture and pasture but not in degraded forests,

179

urban areas and plantations. The primate Arctocebus calabarensis utilises degraded forests

180

only and cannot survive in any other human land use.37 The species utilisation of a particular

181

land use type depends upon its biophysical and life history traits such as diets, body size,

182

temperature range etc.

183

We obtained the list of all species of mammals, birds and amphibians present in each of the

184

804 ecoregions from WWF Wildfinder database.17 Next, for each terrestrial ecoregion j, we

185

counted the number of species of each taxon g for a particular land use type i (e.g. cropland)

186

was coded as suitable in the Habitat Classification Schemes of the IUCN Red List.37 The

187

number of species for which a given habitat was suitable were divided by total number of

188

species of that taxon occurring in the ecoregion, to obtain the fractional species richness 8 ACS Paragon Plus Environment

Page 8 of 32

Page 9 of 32

Environmental Science & Technology

189

(+,, ). The affinity of taxonomic group g to the broad land use type i is then calculated as the

190

proportion of all species that can survive in it (fractional richness), raised to the power 1⁄

191

(see Pereira et al):16 !" ℎ,,

 ,,, =   ,,





= .+,, /

192

Equation 3

193

Taxon affinity to land use intensity types

194

The IUCN Habitat Classification Scheme only provides data on species utilisation of five

195

broad human land uses - plantations, agriculture, pasture, urban areas, and degraded forests.

196

In order to derive taxon affinity to land use intensity classes (ℎ0,,1 ), we therefore combined

197

the calculated affinity to broad land use types (Eq. 3) with information from existing meta-

198

analysis studies26,35 documenting species richness in different land use intensity levels relative

199

to each broad land use class (fractional relative richness, 233 ). The taxon affinity to a land use

200

intensity type (ℎ0,,1 ) is calculated from the affinity to broad land use type as: !" ℎ,, = ℎ,, ∙ .233,, /



201

Equation 4

202

Here the fractional relative richness .233,, / is the local species richness in a land use

203

intensity type (e.g. intense use pasture) divided by the average local species richness in a

204

broad land use type (e.g. pasture). The mean and 95% confidence interval for the 233,, value

205

for minimal, light and intense use pasture, plantations, cropland and urban land uses was

206

obtained from Newbold et al.26 The data for fractional relative richness of the three managed

207

forest intensity classes (reduced-impact logging, selective logging and clear-cut) was

208

obtained from the meta-analysis of Chaudhary et al.35

209

For example, the forestry meta-analysis of Chaudhary et al.35 showed that the average

210

fractional species richness of amphibians in the broad land use class ‘managed forests’ is 0.89 9 ACS Paragon Plus Environment

Environmental Science & Technology

Page 10 of 32

211

while their fractional species richness in forests managed under clear-cut (intense use) regime

212

is 0.544. Therefore the fractional relative richness (233 ) value is equal to: 0.544⁄0.89 = 0.6. In

213

other words, the local species richness in ‘intense use’ (clear-cut) forests is just 60% of the

214

species richness value in broad land use ‘managed forests’.

215

For reptiles and plants, the IUCN habitat classification scheme37 does not yet provide data for

216

all species, and therefore we derive the affinity values solely using the fractional species

217

richness from Newbold et al.26 The affinity of taxa to ‘regenerating secondary vegetation’

218

(under no or minimal use) was also taken directly from Newbold et al.26

219

As an illustration, again consider the ecoregion Biak-Numfoor Rain Forests [WWF ecoregion

220

code: AA0103] that hosts three amphibian species: Papurana papua, Cornufer papuensis and

221

Nyctimystes infrafrenatus. According to IUCN Habitat Classification Scheme37, out of these

222

three species, only one (Nyctimystes infrafrenatus) utilises pastureland. Being an island, the

223

z-value for this ecoregion is 0.44. Therefore, the affinity of amphibians to broad land use

224

!" ‘pasture’ in this ecoregion according to Eq. 3 is: ℎ!9(:

!,(!$; ,90%). In the biodiversity context, selective logging>reduced impact logging. The

441

differences between the CFs for agriculture and pasture intensity classes were less

442

pronounced. These trends reflect the relative order of taxon affinities to different land use

443

intensity types (Table 2).

444

We found that for a given land use type, the new CFs are in general higher than the CFs

445

calculated by Chaudhary et al.4 mainly because their taxon affinity estimates were higher than

446

those calculated in this study. Importantly, the 95% confidence intervals of our new CF

447

values are significantly narrower than the CFs of Chaudhary et al.4 which ranged from

448

negative to positive. Finally, unlike the CFs of Chaudhary et al.4, the new CFs enable

449

comparison of products coming from three different management intensity levels (minimal,

450

light and intense).

451

In order the newly calculated CFs can be operationalized into the LCA softwares, we also

452

provide a look up table that matches the land use intensity categories used in this study with

453

the land use intensity classes recommended by the UNEP-SETAC for use in LCA (Koellner

454

et al.52) and also used in LCA inventory database such as ecoinvent v3 (supplementary Table

455

S7).

456

Application Example

457

Figure 2 shows that the new characterization factors (CFs) are better able to discern the

458

impact of different forest management types as compared to the old CFs provided by

459

Chaudhary et al.4 It can be seen that using the Chaudhary et al.4 CFs for each taxa, the 21 ACS Paragon Plus Environment

Environmental Science & Technology

calculated 95% confidence interval of the potential species loss values under all four regimes

461

overlap with each other as well as cross the CF=0 horizontal line. This means that one cannot

462

conclude if the calculated damage is significantly different from zero (statistically) or if the

463

wood sourced from any one regime is less damaging than the other.

Potential species loss/m 3 (*10 -10)

460

464 465

Figure 2. Biodiversity damage caused by 1m3 of wood sourced from four different forest

466

management regimes in India calculated by using a). New characterization factors derived in this

467

study, b). Characterization factors derived by Chaudhary et al.4 Except timber plantations, the other

468

three management regimes fall into ‘extensive forestry’ category of Chaudhary et al.4 See

469

supplementary Table S8 for detailed calculations (including plants) and data sources.

470

In contrast, the potential species loss for any given taxa, calculated using the new CFs clearly

471

show that wood sourced from plantation forest causes the most damage followed by that

472

sourced from clear-cut, selective logging and reduced impact logging (RIL) regimes because

473

the 95% confidence interval under a particular regime neither overlaps with the other three

474

nor intersects the CF = 0 line (except RIL).

22 ACS Paragon Plus Environment

Page 22 of 32

Page 23 of 32

Environmental Science & Technology

475

For example, using the new CFs, the calculated impact due to land occupation associated

476

with producing a 1 m3 of wood in a selectively logged forest and a forest managed using RIL

477

is equal to 2.04 × 10-10 and 7.09 × 10-12 potential species loss respectively (Figure 2;

478

Supplementary Table S8). In other words, the land occupation impact is 30 times lower if the

479

wood is sourced from RIL forest compared to selectively logged forest. However, the

480

calculated impact using the old CFs of Chaudhary et al. is equal to 5.19 × 10-11 for both

481

selectively logged and RIL forest. This is because both these forestry regime fall under the

482

‘extensive forestry’ definition of Chaudhary et al.4

483

Therefore, the simple example above shows that new CFs are able to differentiate the impact

484

caused by different land use intensity classes, something that was not possible with old CFs.4

485

Note that here we only compared the land occupation impacts of the four forestry regimes but

486

the ranking of total impact might change if other factors such as land use change, fuel use etc.

487

associated with each regime are taken into account.

488

DISCUSSION

489

The study presents first attempts towards addressing the gaps identified by UNEP-SETAC

490

life cycle initiative regarding calculation of land use biodiversity impacts within life cycle

491

assessment (LCA).14 First and most importantly, by using the global land use and land use

492

intensity maps20,32,33, we derived the new characterization factors (CFs) that enables one to

493

discern the damage caused by three levels of intensity classes within a particular broad land

494

use type, something not possible through the CFs of Chaudhary et al.4

495

Note that the UNEP-SETAC recommended to both expand the number of land use classes as

496

well as include different management regimes in the biodiversity model (see page 139 of

497

UNEP report).14 While we do not improve upon the number of land use classes, we do

498

include three management regimes (minimal, light and intense use) for each broad land use

499

class that were absent from Chaudhary et al.4 Second, we derive CFs for plants. 23 ACS Paragon Plus Environment

Environmental Science & Technology

500

Third, for calculating taxon affinity to different land use intensity types, we applied a novel

501

hybrid approach of leveraging Habitat Classification Scheme from the IUCN Red List37 and

502

combining it with empirical data on species richness comparisons from literature.26,35 This

503

enabled us to bring down the uncertainty range of the CFs within a factor of two (Fig. 2;

504

Table S3-S6). Previous CFs from Chaudhary et al.4 had huge uncertainty range varying over

505

several orders of magnitude (from positive to negative) and thereby limiting their practical

506

utility.

507

Fourth, model validation and goodness of fit metrics (Table 3) showed that extinction

508

projections by our parametrized model compare very well with the number of species

509

documented as threatened with extinction by the IUCN Red List44, and much better than the

510

model parametrized by Chaudhary et al.4

511

Finally, we demonstrate the application of new CFs through a case study on comparing the

512

land occupation biodiversity impacts of 1m3 of wood sourced from four different forest

513

management regimes (Fig. 2).

514

The innovations presented here do not resolve all issues identified by the UNEP14, but rather,

515

are steps forward in advancing the field. The main novelty of this study is the hybrid

516

approach of combining the Hoskins et al.20 map with a land use intensity datasets32,33 to

517

calculate the areas of five broad land use types under three intensity levels (Eq. 2, Table 2)

518

and combining the IUCN Habitat Classification Scheme37 with data from meta-analyses26,35

519

on fractional richness to calculate affinity of taxon to the land use types under different

520

intensity levels (Eq. 3, 4). Previous studies (Chaudhary & Brooks28) used the Hoskins et al.20

521

map to calculate the areas of four broad human land use types per ecoregion and the IUCN

522

Habitat Classification Scheme37 to calculate the affinity of mammals, birds and amphibians to

523

these four land use types only.

24 ACS Paragon Plus Environment

Page 24 of 32

Page 25 of 32

Environmental Science & Technology

524

Our intermediate results on the calculation of taxa affinity to different land use intensity types

525

(Table 2) and projected extinctions per ecoregion (Fig. 1a, b) provide important insights for

526

biodiversity conservation and policy making. For example, we found that taxon affinity

527

values differ substantially across the three intensity levels for managed forest land use (Table

528

2). This implies that shifting from minimal use forestry regimes (reduced impact logging) to

529

light/intense use forestry (selective logging/clear-cut) would cause substantial additional

530

losses. This has implications for regions planning to intensify their production systems, or,

531

conversely, to shift them towards biodiversity-friendly management practices.

532

Despite significant improvements, substantial uncertainties in accurately predicting the

533

biodiversity damage caused by a particular land use in different parts of the world remain due

534

to several data gaps, model and parameter uncertainties and the value choices involved in the

535

characterization approaches. Where available, such as in case of z-values19 and fractional

536

relative richness26,35 (fRR, Eq. 4), we propagated the uncertainty in model parameters (Table

537

1) into the characterization factors. However, in many cases it was not possible because the

538

underlying datasets do not report the uncertainty intervals around their estimates such as in

539

the case of the land use and land use intensity maps we used to derive areas of different land

540

use types under minimal, light and intense use20,32,33.

541

We acknowledge that our approach of combining two different land use datasets might add

542

uncertainty to the results that could not be quantified by us. We could only consider five

543

broad land use classes and their three intensity levels but future studies should attempt to

544

include further land use classes (e.g. annual, permanent crops, organic farms etc.) Although

545

rough procedure do exist, such as adopted by Chaudhary et al.4, to translate country level

546

proportion39 of annual and permanent crops to ecoregion level proportion, these are bound to

547

add further unknown amount of uncertainty into the results and therefore not applied here.

25 ACS Paragon Plus Environment

Environmental Science & Technology

548

New global high-resolution, harmonized and validated land use maps with more land use and

549

intensity classes are needed to reduce the uncertainty in characterization factors.

550

Next, the Habitat Classification Scheme of the IUCN37 is not yet available for all species of

551

plants and reptiles and is only available for six broad human land uses for mammals, birds

552

and amphibians. We therefore still had to rely on fractional relative richness factors

553

(233,, , RS. 4) from published meta-analyses26,35 of limited field studies (Eq. 4). A big

554

limitation of the fractional relative richness factors that we used is that they are taxa and

555

region generic. It is possible that one taxa is more sensitive to an increase in intensity of land

556

use than others or that an increase in intensity affects the species differentially in different

557

ecoregions. This might be the reason that the calculated CFs do not differ much across the

558

three intensity classes for cropland and pasture land use (Table S3). For example, according

559

to Newbold et al.26, the average fractional species richness (+,, , Eq. 3) in ‘minimal use’ and

560

‘intense use’ cropland is 0.73 and 0.64 respectively, i.e. a difference of just 9%.

561

Consequently, the calculated global CFs for cropland vary