Deriving the Metal and Alloy Networks of Modern ... - ACS Publications

Feb 29, 2016 - Philip Nuss,. †,∥. Wei-Qiang Chen,*,†,§ and Thomas E. Graedel. †. †. Center for Industrial Ecology, Yale School of Forestry ...
3 downloads 0 Views 3MB Size
Subscriber access provided by Universiteit Utrecht

Article

Deriving the Metal and Alloy Networks of Modern Technology Hajime Ohno, Philip Nuss, Wei-Qiang Chen, and Thomas E. Graedel Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.5b05093 • Publication Date (Web): 29 Feb 2016 Downloaded from http://pubs.acs.org on March 15, 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 32

Environmental Science & Technology

1

2

Deriving the Metal and Alloy Networks of Modern Technology

3

The authors:

4

Hajime Ohno *

5 6

Center for Industrial Ecology, Yale School of Forestry & Environmental Studies, 195 Prospect St, New Haven, CT 06511, USA

7 8 9 10

Graduate School of Engineering, Tohoku University, 6-6-04, Aramaki Aza Aoba Aoba-ku, Sendai, Miyagi 980-8579, Japan [email protected] +81+22-795-8569

11 12

Philip Nuss

13 14 15 16

Center for Industrial Ecology, Yale School of Forestry & Environmental Studies, 195 Prospect St, New Haven, CT 06511, USA [email protected]

17

Wei-Qiang Chen **

18 19

Center for Industrial Ecology, Yale School of Forestry & Environmental Studies, 195 Prospect St, New Haven, CT 06511, USA Institute of Urban Environment, Chinese Academy of Sciences, 1799 Jimei Road, Fujian 361021, P.R.China [email protected] +86-592-6190-763

20 21 22 23 24 25

Thomas E. Graedel

26 27 28

Center for Industrial Ecology, Yale School of Forestry & Environmental Studies, 195 Prospect St, New Haven, CT 06511, USA [email protected]

29 30 31 32

* Corresponding Author: Tel.+81-22-795-5869, Email: [email protected] ** Co-Corresponding Author: Tel. +86-592-6190-763, Email: [email protected]

1 ACS Paragon Plus Environment

Environmental Science & Technology

33

Abstract

34

Metals have strongly contributed to the development of the human society. Today, large amounts

35

of and various metals are utilized in a wide variety of products. Metals are rarely used

36

individually but mostly together with other metals in the form of alloys and/or other

37

combinational uses. This study reveals the inter-sectoral flows of metals by means of

38

input-output (IO) based material flow analysis (MFA). Using the 2007 United States IO table, we

39

calculate the flows of eight metals (i.e., manganese, chromium, nickel, molybdenum, niobium,

40

vanadium, tungsten and cobalt) and simultaneously visualize them as a network. We quantify the

41

interrelationship of metals by means of flow path sharing. Furthermore, by looking at the flows

42

of alloys into metal networks, the networks of the major metals iron, aluminum, and copper

43

together with those of the eight alloying metals can be categorized into alloyed-,

44

non-alloyed-(i.e. individual), and both mixed. The result shows that most metals are used

45

primarily in alloy form and that functional recycling thereby requires identification, separation,

46

and alloy-specific reprocessing if the physical properties of the alloys are to be retained for

47

subsequent use. The quantified interrelation of metals helps us consider better metal uses and

48

develop a sustainable cycle of metals.

49

50

Key words: input-output analysis, material flow analysis, alloying metals, alloys

51

52

2 ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32

53

Environmental Science & Technology

TOC

54

55

3 ACS Paragon Plus Environment

Environmental Science & Technology

56 57

1. Introduction Metals have been indispensable in our lives since the first invention of metal goods during

58

the Stone age1. The use of bronze during the Bronze Age marks the starting point for the use of

59

metals in alloy form1. For example, raw steel (an iron-based alloy containing up to 2% carbon) is

60

oftentimes used with varying alloying elements (e.g., chromium, manganese, molybdenum, and

61

vanadium) to obtain various physical and chemical properties, such as hardness, workability,

62

corrosion resistance, and heat resistance.

63

Because modern society is highly dependent on the use of various minerals and metals,

64

various studies in recent years have investigated the topic of resource criticality2-6. The goal of

65

these criticality studies is to provide insights into the geological distribution of metals and related

66

geopolitical, social, and regulatory aspects in the sourcing countries, companionality (i.e., the

67

fraction of an element obtained as a byproduct with another metal)7, environmental implications8,

68

aspects of substitutability9, recycling potential10, and the overall importance of the resource to

69

the global economy. Increasing recycling of metals can support the sustainable use of metals.

70

However, recovering an individual metal from its alloyed form can be difficult under current

71

technical and economic conditions11, 12. Several studies on the thermodynamic behavior of metals

72

in pyro-metallurgical processes discuss the difficulties associated with metal recovery from

73

major metal-based alloys (e.g. steel or aluminum- and copper-based alloys)12-15. These are

74

related to a metal’s thermodynamic behavior in pyro-metallurgical processes, e.g., whether a

75

metal; remains in the metal matrix, is being oxidized to the slag phase, or vaporized to the gas

76

phase. Except metals vaporization (e.g. zinc in steel recycling)16, commercial base technologies

77

for recovering individual metals from alloys have not yet been well developed. Therefore, as

78

long as alloys are recycled for the purpose of mainly recovering base metals, alloyed elements

4 ACS Paragon Plus Environment

Page 4 of 32

Page 5 of 32

Environmental Science & Technology

79

are unintentionally dissipated as contaminants and/or by-products17. To prevent the dissipation of

80

alloying metals during recycling, recent research points out the importance of the

81

quality-oriented scrap allocation sorting scrap according to contents of metals for the sustainable

82

metals use, especially for aluminum alloys18, 19. However, in order to implement the appropriate

83

scrap allocation for various metals and alloys, understanding metals and alloys usages in a

84

society would be necessary.

85

A potential tool for the detailed study of inter-sectoral flows in an economy is the use of

86

input-output tables (IOT) which are constructed on a regular basis by national statistical

87

offices20. In such tables, the rows and columns indicate individual sectors of the economy and

88

the matrix elements reflect transactions between two sectors in a given year. The level of detail

89

differs from country to country – in the United States the IOT is given at 389 x 389 matrix

90

resolution in 200721. Although the use of some 400 sectors would seem more than adequate to

91

describe a national economy, the 2007 IOT of the United States contains independent sectors for

92

only three metals, i.e., iron, copper, and aluminum, and even for these metals the entries are in

93

monetary rather than physical units. Nonetheless, it is possible to convert the monetary

94

transactions in IOTs to physical flows of metals22, 23 and thereby generate an inter-sectoral

95

network for metals characterized by individual sectors (see our companion paper 24 where we

96

derive a metal network for aluminum in the 2007 US economy).

97

For metals without specific IOT sector entries (i.e., elements other than aluminum, iron,

98

and copper), a different approach is required. For those metals, the IOTs provide only aggregated

99

information, e.g., the “Nonferrous metals” sectors which comprise data about several metals

100

within a single sector.

5 ACS Paragon Plus Environment

Environmental Science & Technology

101

However, IO based material flow analysis (MFA) provides some methods to help with the

102

disaggregation of sectors in an IOT based on the waste input-output material flow analysis

103

(WIO-MFA) model22, 23. Specific examples include plastics25, steel alloying elements17, 26, and

104

automobile parts27. Of particular interest here is the work of Nakajima et al.26, which visualized

105

flows of chromium, nickel, and molybdenum of automobile production in Japan.

106

Based on the above work and the approach described in our companion paper24 for

107

aluminum, we derive the metal networks of eight “alloying” metals (i.e., manganese, chromium,

108

nickel, molybdenum, niobium, vanadium, tungsten and cobalt) in the 2007 US economy. The

109

results are visualized using a Sankey diagram covering the entire economy-wide supply chain of

110

these metals. The eight metals were chosen because they have large consumptive uses and

111

because they are typically related to the three major metals in the form of alloys28 (Some relevant

112

characteristics of these metals are given in Table S1 in the Supporting Information). As

113

Nakamura et al.29 pointed out, “Sankey diagrams may become too complex for effective

114

visualization” for several hundreds of sectors. An additional challenge is that we have multiple

115

metals to be simultaneously visualized in a Sankey diagram. However, by applying the approach

116

of Chen et al.24 (a companion paper) to aggregate and simplify the number of sectors as well as

117

distributing sectors according to their degree of fabrication, the readable simultaneous picture of

118

metal flow can be visualized. The diagram then provides a comprehensive view of major metal

119

uses in society. In addition, by regarding sectors and flows as nodes and edges (directed and

120

weighted) respectively, the flows of metals can be regarded as inter-sectoral networks. “Edge” is

121

common terminology in network analysis to express the flow of resources or information

122

between two nodes30. Here the nodes are sectors and the edges represent either monetary or mass

123

flows. As Nuss et al.31 (another companion paper) demonstrates, various measures of network

6 ACS Paragon Plus Environment

Page 6 of 32

Page 7 of 32

Environmental Science & Technology

124

analysis can then be applied to the IO-based material flow networks (IO-MFNs) of metals to gain

125

insights into sector “importance”, supply chain risk, and the overall network structure. This

126

information is therefore also relevant in the context of resource criticality analysis31. In this

127

paper, we derive and compare metal IO-MFNs in terms of the sharing of edges. Besides showing

128

the overall network diagrams, we quantitatively study to what degree the metals coexist with

129

each other in the economy by analyzing the edge sharing among metals. The coexistence of

130

metals measured by edge sharing then reflects the challenges of any subsequent separation of

131

metals and alloys in recycling operations.

132

133

2. Methodology

134

2.1 Conversion of the monetary IOT to a metal IOT

135

The multi-step processes of converting the monetary IOT to a metal IOT and finally

136

obtaining an IO-MFN are outlined in Figure 1. A very detailed description of applying these

137

processes to aluminum has been provided in Chen et al24. Here we provide a shorter and general

138

version for the convenience of readers. The first step involves the disaggregation of sectors and

139

unit conversion of IO table  from monetary to physical units to obtain a hybrid IO table  ∗ .

140

Details are provided in section 2.2. Secondly, an input coefficient matrix  is calculated20.

141

Thirdly, the input coefficient matrix is multiplied with two filtering matrices, i.e., (1) a physical

142

flow filter Φ which is a binary matrix removing non-physical flows (e.g., services) and physical

143

flows that do not incorporate the studied metal into the final product (e.g., process catalysts), and

144 145 146

(2) the yield loss filter matrix Γ that removes the mass of inputs that become process waste. The filtered input coefficient matrix  = Γ ⊗ (Φ ⊗ ) then consists of  (which represents input coefficients whose ( , )-element is the monetary or physical input of Product to 7 ACS Paragon Plus Environment

Environmental Science & Technology

147

Page 8 of 32

Product ) and  (which represents input coefficients whose ( , )-element is the physical

148

input coefficient of Material to Product ). The unit of  is tonne per monetary unit (but

149

see below for a discussion on IO-MFN units). In the fourth step, the waste input-output material

150

flow analysis (WIO-MFA) model23 is computed to calculate the matrix of material composition

151

152 153 154 155

of products ( ) by:

 =  ( −  )

(1)

As Nakamura et al.23 noted, the IOT representing inter-sectoral transactions can readily be converted into the corresponding flow of materials. Consequently, during the fifth step the

hybrid IO table can be simply transformed into a physical IO table   for metals and alloys by formulating the product of the diagonalized  of a material (i.e. metal)  and the IOT as: ∗ diag ( (,∙))  =   ∗ 

156 157 158 159 160

∗ Here,  ∗ represents the endogenous IO matrix obtained in the second step, consisting of  ∗ and  with the same classification as  (for the detailed definition, see Supporting

Information).  (,∙) represents the th row of the matrix  . Let

be the vector of the

net trade of products and !(,∙) the mass of waste material  generated in the supply chain. Then the IO-MFN of metal  in an economic system   can be expressed as: ∗ diag ( (,∙))  =  ∗ 

161

(2)

diag ( (,∙)) (, )

!(,∙)  0

We address the calculation sequence in more detail in the Supporting Information.

162 163

2.2 Sector disaggregation and arrangement for the model application to the U.S. IOT

164

2.2.1 Definition of sectors

8 ACS Paragon Plus Environment

(3)

Page 9 of 32

Environmental Science & Technology

165

Although a large variety of metals are used by modern society, their trade in the economy

166

is highly aggregated in the original IOT for the U.S. into only five sectors (the five sectors listed

167

in the first column of Table 1). The sector for “Iron and steel mills and ferroalloy manufacturing”

168

includes sub-sectors producing raw materials based on iron, many types of ferroalloys, crude

169

steels, and finished steels. Similarly, the sector “Primary smelting and refining of nonferrous

170

metals” aggregates industries producing many different metals. The sector, “Nonferrous metal

171

rolling, drawing, extruding, and alloying” aggregates the nonferrous processing industries. In

172

order to disaggregate flows from those sectors into flows related to a specific nonferrous metal,

173

these sectors were disaggregated into the sectors indicted in the second column of Table 1. As a

174

consequence of the disaggregation, the five metal-related sectors were disaggregated into 34

175

sectors related to individual metals, metal materials, scrap, and products. The disaggregation

176

made the iron flow clear as well because iron sources were separated from the original

177

aggregated sector for iron.

178

Sectors of metal raw materials such as ferroalloys and metals themselves were

179

categorized as "Material" sectors for the analysis. The ferroalloys in the “Material” category

180

were considered as sources of both iron and corresponding alloying elements according to their

181

chemical composition28. Consequently, flows of iron and an alloying element incorporated in a

182

ferroalloy were separately analyzed. In addition, sectors of three different grades of crude steels

183

and nonferrous alloys were defined as "Fabricated material", which is midway between Material

184

and Product in the model. The distinction of this intermediate category enables us to analyze the

185 186

flows of alloys defined as “Fabricated material” by setting their inputs as products in 

instead of inputs as “Material” in the model23. In other words, when we describe direct inputs of

9 ACS Paragon Plus Environment

Environmental Science & Technology

187 188

"Material" sectors in  , we can obtain individual metal flows; when we utilize "Fabricated material" sectors in  , alloy flows can be obtained.

189 190 191

2.2.2 Disaggregation of sectors The disaggregation of sectors is conducted for both rows and columns of the IOT. In

192

addition to the disaggregation, values in the rows are converted to physical units. An exception

193

to this approach relates to the sector “Nonferrous metal (except copper and aluminum) rolling,

194

drawing, extruding, and alloying”, which is disaggregated into two sectors but remains in

195

monetary units because the total production of nonferrous alloys and superalloys in physical

196

units could not be obtained from the available data (Supporting Information).

197

The approach to the disaggregation of sectors highly depends on the data availability in the

198

region for which the IOTs are provided. Nakajima et al.26 and Ohno et al.17, 27 used the Japanese

199

IOT 32. In this work we used the U.S. IOT21. These tables have the highest resolutions compared

200

to IOTs of other economies, and can be used together with data such as the year book of iron and

201

steel for Japan 33 and the minerals year book for the U.S.28, respectively. If only lower IOT

202

resolution is available, similar disaggregation of sectors could be challenging for other countries

203

or regions. The disaggregation of sectors provides many benefits for metal flow analysis, but

204

introduces uncertainties as well, because assumptions and estimations are necessary when

205

adjusting the available data into the form of an IOT having around 400 sectors. The detailed

206

operations of the disaggregation taken in this study are explained in the Supporting Information.

207 208

Because the assumption of unique sectoral prices34 for metals is generally not satisfied in the sectors of the IOT, the direct conversion of monetary unit to physical unit may introduce

10 ACS Paragon Plus Environment

Page 10 of 32

Page 11 of 32

Environmental Science & Technology

209

uncertainties into the analysis. For example, when the aggregated “Iron and steel mills” sector

210

sells 1 million USDs of products to two different downstream sectors such as pig iron and

211

stainless steel, the physical flow distribution will vary depending on price differences between

212

pig iron and stainless steel, respectively. This problem can be partly reduced by further

213

disaggregating the economic sectors. However, given the current resolution of the US IO table

214

used in this study (~400 sectors), direct conversion of monetary units to physical units using

215

homogeneous sectoral prices will likely cause some inconsistences between the analysis and

216

reality. In this regard, rather than expressing the transformed metal units (TMUs) in mass terms,

217

calculated units are formally designated in “hektes” as discussed in our companion paper24. This

218

new unit expresses the relative magnitudes of physical flows without implying that an actual

219

physical unit such as kilograms is appropriate.

220 221

2.3 Comparison of IO-MFNs

222

The resulting networks for the eleven metals and seven alloys can be studied by regarding

223

each sector in the IOT as a node and the input of each sector to another as an edge. Because the

224

networks were built in the same inter-sectoral structure, we can compare the IO-MFNs in terms

225

of edge sharing. One edge property involves the number of shared edges, in which the degree of

226

edge sharing among metals represents their potential coexistence in technological uses. For

227

example, when the networks of Fe and Cu (having 57,321 and 31,654 edges, respectively, see

228

Table 2), share 31,654 edges with each other, the degrees of edge sharing between Fe network

229

and Cu network are calculated as 55% and 100% for Fe and Cu respectively. A second network

230

property takes the total strength of shared edges into account, thus representing the significance

231

of the relationships among metals in the whole network. This measure takes into account the

11 ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 32

232

weighting of edges (i.e., the quantity of physical flow of metal between sector as expressed in

233

hektes). For example, Fe and Cu are found to have 307 kh (kilohektes) and 5.4 kh of total edge

234

strength, respectively. If shared edges have 114 kh of edge strength in the Fe network and 5.4 kh

235

of edge strength in the Cu network, the weighted degree of edge sharing will be 37% for Fe and

236

100% for Cu. In this case, 37% of the total mass of iron appears to be related to the network of

237

copper.

238

For clarity in examining the flows of metals, we applied a threshold cutoff of 1 hekte

239

(otherwise almost all edges are shared by all metals, because nearly all sectors employ tiny

240

amounts of all metals). The quantitative effect of applying this threshold to the networks is

241

summarized in Table 2. There is substantial coherence across the metals in terms of edge

242

strength, but not of total edges. This indicates that the preponderance of metal flows occurs

243

together from early to late sector flows, but that small amounts of metals such as cobalt and

244

tungsten see diverse use in sectors other than those involving the major metals.

245 246 247 248

 ∗ We utilize  = diag( (,∙)) in   which represents a metal IOT including

Product to Product of the metal flow (i.e. inter-sectoral supply chain) to compare the IO-MFN of metals with each other. From set theory, edge sharing between two different metal ( ∈ α, β) networks (N % , N & ) can be expressed as follows:

Edge sharing = N % ∩ N &

249 250 251 252

(4)

Hereafter, metal IO-MFNs themselves are expressed as N  to distinguish them from the supply  chains of metals,  .

To quantify the degree of edge sharing based on the IOT that was obtained, let V  be

 the edge detecting binary matrix for metal “m”, having the ( , )-element v01 defined as

12 ACS Paragon Plus Environment

Page 13 of 32

Environmental Science & Technology

 v01 =2

253 254 255 256

 1 ( , ) − element of  >1 0 otherwise

(5)

  > 1 represents the threshold for removing edges with strengths less than 1 hekte from the

network. (This table is also defined as the adjacency matrix or Boolean matrix20). By comparing matrices V  of each metal pair, the degree of edge sharing can be calculated. Here, ∑0 ∑1 v 

represents the numbers of edges in an IO-MFN for metal . Because V  contains only 1 and 0,

257

the result of subtraction of V  for metal α and metal β provides 1, 0, and -1 values defined as

258

follows: %& ∆v01

1 N> ∖ N@ BBBBBBBBBBB > ∪ N@ = = 0 N > ∩ N @ or N −1 N@ ∖ N>

259

By summing ∆V %& = C∆v01

260

Finally, the degree of edge sharing can be calculated:

%&

D∆v01

%&

S = G1 −

261 262

= 1E, the numbers of unshared edges can be determined.

∑0 ∑1 ∆V %& H × 100 ∑0 ∑1 ∆V %

 element-wise product of  with V  , representing metal flows over 1 hekte and ∆J %& the

element-wise product of ∆V %& with J  , including the weights of shared edges. The

264

weighted degree of edge sharing is then ∑0 ∑1 ∆J %& SL = G1 − H × 100 ∑0 ∑1 ∆J %

266

(7)

The degree of weighted edge sharing can be calculated as well. Let J  be the

263

265

(6)

(8)

We can distinguish the two types of coexistence by comparing the networks of metals with those of alloys. Because alloys contain metals, metals contained in an alloy share the same

13 ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 32

267

network as the related alloy. Thus, the supply chain of a metal in alloy forms, N  0M NOOPQ , can be

268

calculated as follows: N ∗  = diagR (S,∙)T

 0M NOOPQU 

269

W

N = V  (,∙) NX

(9)

(10)

N where S denotes a specific alloy and  represents the supply chain of alloy S. For metal

270

networks not contained in alloys, the network Y  MPMNOOPQ can be obtained based on the

271

supply chains of metals not contained in alloys: 

 MPMNOOPQ

272

 =  − 

 0M NOOPQ

(11)

These two networks and the original networks satisfy the following relationship: N  = N  0M NOOPQU ∪ N  MPMNOOPQ

(12)

273

According to this relationship, the edge sharing among metals can be categorized into three

274

groups; A: sharing edges only as alloy forms, either alloyed together or combined physically as

275

alloy forms (i.e., A = N  0M NOOPQU ∖ N  MPMNOOPQ ); B: sharing edges only in non-alloy form (i.e.,

276 277

B = N  MPMNOOPQ ∖ N  0M NOOPQU ); and C: sharing edges in both alloy and non-alloy forms (i.e., C = N 0M NOOPQU ∩ N  MPMNOOPQ ).

278

279

3. Results and Discussion

280

3.1 IO-MFN visualization

281

The derived IOTs are commonly presented as numerical spreadsheets. However, graphical

282

approaches are useful methods for illustrating some of the most interesting features of the tables.

283

An example is the IO-MFN for cobalt, shown in Figure 2. For all network diagrams in this paper,

14 ACS Paragon Plus Environment

Page 15 of 32

Environmental Science & Technology

284

sectors in the IOT were aggregated into 106 nodes (which is more than the 99 nodes in Chen et

285

al24 due to sector disaggregation). The layout of the nodes is according to their place in the

286

fabrication sequence, that is, sectors at the left side of the diagrams represent early stages of

287

fabrication while those at the right side reflect later stages. The U.S. did not mine cobalt in 2007

288

28

289

of the figure). The imported cobalt had two principal destinations. One was its use in Batteries

290

(the pale blue arrow from lower left to the right center of the network). The other was the

291

formation of non-ferrous alloys (largely nickel-based superalloys) (the light blue line from lower

292

left to left-center), followed by Forging & Stamping and then by the flow to the Aerospace sector

293

for use in aircraft engines. More careful examination reveals such features as imports of batteries,

294

flows of cobalt-containing products to the Machinery sector, the transfer and reuse of new scrap,

295

and a web of cobalt-containing flows to a very large number of sectors in the national economy.

296

Similar diagrams for manganese, chromium, nickel, molybdenum, niobium, vanadium, and

297

tungsten appear in the Supporting Information (Figure S2 to S8). The flows coming directly from

298

a metal node correspond to relatively well known flows of metals because their values are

, but imported a significant amount (the pale orange arrow from the top left to the bottom left

299

∗ included in the IOT based on use statistics (i.e.,  in equation (2)). However, flows to

300

metal-containing products could not be directly obtained from statistics or by process

301

based-MFA. Thus, Figure 2 demonstrates that more detailed flow paths and end-uses were

302

uncovered in this research than have traditionally been identified. This feature of the results,

303

which is present over all parts of the metal supply chain, will help to understand how metals flow

304

through industrial networks, and to anticipate where metals tend to accumulate as potential

305

“urban mines”.

15 ACS Paragon Plus Environment

Environmental Science & Technology

306

Similar diagrams can be constructed for the metals in alloys, providing that metal-specific

307

IO-MFNs have been derived for all of the alloying metals of interest. As pictured

308

diagrammatically in the TOC art, the task is to combine common matrix elements and flows for

309

each of the alloying metals, continuing until all matrix elements have been addressed. When this

310

is done, the results are as shown in Figure 3(a). This simultaneous visualization helps us

311

recognize the coexistences of metals in products and society, and reveals that most metals flow

312

together via alloy sectors. In the present set of eight alloying metals most of the nodes of the

313

network are involved, but some metals are clearly more important than others. Figure 3(b) shows

314

an enlarged version of the right side of Figure 3(a); it demonstrates that most of the metals are

315

involved in flow to the Motor Vehicle sector and thus to Road Transportation. The Construction

316

sector is also important in that regard, and both of those important sectors utilize metals

317

primarily in alloy forms. In contrast, Machinery consumes less alloyed tungsten provided by the

318

Metalworking Machinery sector: that sector largely employs tungsten as tungsten carbide for

319

metal-cutting tools 28. This result implies that the potential tungsten recovery from the

320

Metalworking Machinery sector will be in the form of tungsten carbide rather than as a

321

constituent of an alloy.

322

In addition to following each metal in an alloy group independently, the metal flows can be

323

summed to create a network for a small number of alloy groups (Supporting Information: Figure

324

S9). This approach provides information specific to each alloy group, such as stainless steel or

325

alloy steel. For stainless steel, for example, the two largest flows include the one to Architectural

326

Metals through to the Construction sector, and the other to Motor Vehicle Parts through to the

327

Motor Vehicle sector. However, smaller flows go to many other sectors, as can be better seen in

328

the enlarged version of the center of Figure S9(a) shown in Figure S9(b). The network diagram 16 ACS Paragon Plus Environment

Page 16 of 32

Page 17 of 32

Environmental Science & Technology

329

from the perspective of alloy groups implies that alloys coexist as such in products. However, the

330

important difference in the coexistences between co-existing metals and the alloys themselves is

331

that different metals flowing to the same sectors may not be chemically combined (i.e., alloyed)

332

with each other but rather are physically combined in products. Physically combined (e.g.,

333

bolted) metals would be expected to be more recyclable than those that are chemically

334

combined12.

335

3.2 Comparison of IO-MFNs

336

Because most of the metals addressed in this study flow together and share edges, metal

337

relations of N % ∩ N & is summarized in Figure 4 by quantifying the degree of edge sharing and

338

weighted edge sharing including IO-MFN of base metals obtained as a consequence of the sector

339

disaggregation. Edges are weighted in units of hektes (see above). Figure 4(a) clearly

340

demonstrates that the degree of edge sharing is related to the total consumption of metals. Major

341

metals’ networks tend to be shared by minor metals, but iron is a central component of most

342

metal networks because other metals share virtually all their edges with iron (i.e. N ]^ ⊃ N  ).

343

Looking at other combinations of metals, for example, nickel shares 3% of its edges with those

344

of cobalt, and cobalt shares 86% of its edges with those of nickel. This reflects that a large

345

portion of cobalt is jointly used with nickel as a superalloy. However, this joint use represents a

346

tiny part of total nickel use as the main use of nickel is for stainless steel, which generally does

347

not contain cobalt. In the weighted edge sharing results shown in Figure 4(b), we find that some

348

metals share large fractions of weighted edges with other metals, while sharing much smaller

349

fractions of numbers of edges. For example, iron shares 86% of edge weight (metal flow) with

350

vanadium, but just 2% of edges. Because vanadium contents in carbon steel and tool steel are so

351

small, the flows of vanadium tend to be cut off by the threshold whereas the edges are shared 17 ACS Paragon Plus Environment

Environmental Science & Technology

352

with iron. However, this result also demonstrates that some alloying metals are likely to be

353

dissipated by being involved in the flow of mass-produced metals.

354 355

The relationships in the network shown in Figure 4 include both chemically and physically combined coexistences of metals. Figure 5 shows the degree of edge sharing and weighted edge

356

sharing within the whole network N  for the three categories. Although most metals tend to

357

share edges in alloy-related networks (i.e., the categories A and C) rather than in the non-alloy

358

network, tungsten has 78% and 55% of edges and edge weight in its own non-alloy network.

359

Furthermore, tungsten shares fewer edges with other metals in the non-alloy network. This

360

reflects the fact that tungsten in modern technology exists in relatively pure forms without

361

significant contamination by other metals (e.g., tungsten carbide cutting tools). Because

362

relationships among metals in category B include only physical coexistence, this implies

363

enhanced potential for metal recovery in the recycling phase. Metals other than tungsten tend,

364

however, to be highly combined with each other in alloy-related forms, especially in category A.

365

In this case, any subsequent recycling must consider the metals as largely part of alloy

366

combinations. Recycling thereby requires identification, separation, and alloy-specific

367

reprocessing if the characteristics of the alloys are to be retained for subsequent use. If this is not

368

done, the result is likely to be the incorporation of the alloy into a composite recycling flow of

369

(for example) carbon steel, in which case the physical and chemical properties for which the

370

alloy was designed will be lost.

371 372

373 374

3.3 Future directions We have demonstrated in this work a methodology to use IOTs as a starting point from which to construct IOTs and IO-MFNs for those metals that lack sectors of their own in 18 ACS Paragon Plus Environment

Page 18 of 32

Page 19 of 32

Environmental Science & Technology

375

economic IOTs. We have presented the results for eight metals that are associated with iron,

376

copper, and aluminum, the only three metals with designated sectors in the U.S. IOTs. Together

377

with the work of Chen et al.24 , this results in eleven different metals for which U.S. 2007

378

networks are available. Similar results for other regions and years can be expected. We have also developed approaches to combine sets of metals into IO-MFNs for alloy

379 380

groups of interest. This information, which is complementary to the networks of individual

381

metals, permits the analysis of the principal and subsidiary flows of alloys in an economy. This

382

work also shows by the degree of edge sharing that all the metals analyzed herein are

383

predominantly used in alloy forms and thus must be recycled in those forms rather than as

384

individual metals. This result has important implications for the recovery of metals in discarded

385

products. In connection with efficient scrap treatment and sorting for end-of-life (EoL) vehicles19,

386

27

387

included in EoL products would seem to be required for the sustainable use of metals.

388

, the development of metal scrap recycling strategies that consider combinations of metals

More complex network analysis approaches can reveal additional features of the networks

389

generated in this work31. For example, we anticipate that triangulated physical IOTs35, 36 for

390

metals will reveal more complex interrelationships among both metals and industrial sectors. In

391

addition, as Kagawa et al.37 and Okamoto38 have demonstrated, the identification of important

392

industry clusters for metal uses will be a potential application as well. While we have carried out

393

this type of detailed clustering analysis for aluminum31 doing similar analysis for the metals

394

reported in the present remains a task for the future.

395

19 ACS Paragon Plus Environment

Environmental Science & Technology

396 397

Supporting Information The detailed methodologies for the conversion of monetary IO table to a metal IO table, the

398

disaggregation of sectors, the results for individual IO-MFN of metals, and the figure for the

399

network of alloys can be found in the Supporting Information. This material is available free of

400

charge via the internet at http://pubs.acs.org.

401

402

Author Information

403

Corresponding Authors

404

*Phone: +81-22-795-8569; email: [email protected]

405

*Phone: +86-592-6190-763; e-mail: [email protected]

406

Notes

407

The authors declare no competing financial interests.

408

409 410

Acknowledgements We would like to thank the Yale Criticality Project for funding this research.

411

This research was supported by Japanese Society for the Promotion of Science (Grant-in-Aid for

412

JSPS Fellows 258801).

413

20 ACS Paragon Plus Environment

Page 20 of 32

Page 21 of 32

Environmental Science & Technology

414 415

Table 1. Disaggregated sectors and their category definition for WIO-MFA Original sector

Iron and steel mills and ferroalloy manufacturing

Alumina refining aluminum production

and

primary

Primary smelting and refining of copper

Primary smelting and refining of nonferrous metal (except copper and aluminum)

Sector after disaggregation Pig iron and DRI (direct reduced iron) Iron and steel scrap Ferromanganese Ferrochromium Ferronickel Ferromolybdenum Ferroniobium Ferrovanadium Ferrotungsten Other ferroalloys Crude steel-Carbon steel Crude steel-Stainless steel Crude steel-Alloy steel Hot-rolled steel-carbon steel Hot-rolled steel-stainless steel Hot-rolled steel-alloy steel Primary aluminum Aluminum scrap Aluminum alloy Primary copper Copper scrap Copper alloy Manganese Chromium Nickel Molybdenum Niobium Vanadium Tungsten Cobalt Nickel scrap Stainless steel scrap Nonferrous metal alloy*1

Nonferrous metal (except copper and aluminum) rolling, drawing, extruding, Superalloy*1*2 and alloying *1 Expressed in monetary units *2 Analyzed with special treatment (see SI) 416 417

21 ACS Paragon Plus Environment

Category of sector Material Material Material Material Material Material Material Material Material Material Fabricated material Fabricated material Fabricated material Product Product Product Material Material Fabricated material Material Material Fabricated material Material Material Material Material Material Material Material Material Material Material Fabricated material Fabricated material

Environmental Science & Technology

Page 22 of 32

418 419 420 421

Table 2. Numbers of edges and edge strengths of the metal networks Fe

Cu

Mn

Cr

Ni

Mo

Nb

V

W

Co

Total edges

74,608

73,948

73,948

74,608

74,608

74,606

74,608

74,278

74,608

74,278

74,278

>1 h (threshold) Remaining edges ratio (%)

57,321

39,590

31,654

20,052

17,733

12,835

3,533

1,778

1,123

1,401

738

76.83

53.54

42.81

26.88

23.77

17.20

4.74

2.39

1.51

1.89

0.99

306,890

20,572

5,458

1,626

1,255

560

46.7

26.7

13.6

10.6

13.7

Total >1 h (kh) 306,886 Remaining edge weight ratio(%) 100.00

20,565

5,451

1,618

1,247

553

42.8

24.1

11.6

8.4

12.1

99.97

99.87

99.52

99.39

98.75

91.57

90.38

85.17

77.81

88.02

8,660

2,270

588

260

196

21

9.9

5.3

11

9.0

Total edge * weights (kh )

Consumption (kt)

422

Al

*

96,500

kh = kilohekte (see text)

423

22 ACS Paragon Plus Environment

Page 23 of 32

Environmental Science & Technology

424 425 426

Figure Captions

427 428

Figure 1. The process of filtering the network for metal M in an economic system defined by a sectoral input-output table.

429 430 431 432

Figure 2. The cobalt input-output network for the United States, 2007. The pale orange dots at the top of the diagram are origin locations for flows of cobalt or cobalt-containing products into the U.S. economy. Arrow widths indicate approximate percentages of flows (i.e., edge strengths). The units of flow are hekhtes, discussed in the text.

433 434 435

Figure 3. (a) Networks of eight metals commonly used in alloys in the U.S. economy, 2007. Colors represent the individual metals; arrow widths represent edge strengths; (b) Expanded view of (a) to show detail in the Motor Vehicle and Construction sectors.

436 437 438 439 440

Figure 4. Percentage of shared edges between networks of the eleven metals disaggregated in this work. (a) Shared edge percentages for networks of the individual metals [e.g., for non-alloy uses]; (b) Shared edge percentages for alloy networks. By reading numbers for the horizontal direction of the tables, the degree to which a metal shares edges with other metals can be observed.

441 442 443 444 445 446 447

Figure 5. (top row): Percentages of shared edges for metals either alloyed together or combined physically as alloy forms(left top panel), for metals sharing edges in both alloy and non-alloy forms (center top panel), or for metals sharing edges only in non-alloy form (right top panel); (bottom row): Percentages of shared edge strengths for metals either alloyed together or combined physically (left bottom panel), for metals sharing edge strengths in both alloy and non-alloy forms (center bottom panel), or for metals sharing edge strengths only in non-alloy form (right bottom panel).

448

23 ACS Paragon Plus Environment

Environmental Science & Technology

449 450

Figure 1

451

24 ACS Paragon Plus Environment

Page 24 of 32

Page 25 of 32

Environmental Science & Technology

Figure 2

25 ACS Paragon Plus Environment

Environmental Science & Technology

Figure 3 (a)

26 ACS Paragon Plus Environment

Page 26 of 32

Page 27 of 32

Environmental Science & Technology

Figure 3 (b) 27 ACS Paragon Plus Environment

Environmental Science & Technology

(a) Degree of edge sharing

(b) Weighted degree of edge sharing

% Fe Al Cu Mn Cr Ni Mo Nb V W Co

Fe Al Cu Mn Cr Ni Mo Nb V W Co 100 37 37 100 99 40 92 86 86 24 20 100 100 58 99 99 98 69 34 30 61 51 100 100 100 98 99 98 57 28 22 63 44 100 40 39 100 100 53 92 80 83 31 27 100 43 43 100 100 99 92 85 77 28 23 100 50 49 100 100 100 90 80 73 36 31 100 45 45 100 100 91 100 61 82 65 56 100 35 35 100 100 51 100 100 94 35 31 100 29 29 100 100 55 100 75 100 48 44 100 86 85 99 99 97 86 59 56 100 65 100 62 62 98 100 100 91 78 78 79 100

Figure 4 (a) and (b)

28 ACS Paragon Plus Environment

Page 28 of 32

Page 29 of 32

Environmental Science & Technology

Figure 5

29 ACS Paragon Plus Environment

Environmental Science & Technology

Page 30 of 32

References 1. Gleich, A. v.; Ayres, R. U.; Gössling-Reisemann, S., Sustainable metals management : securing our future - steps towards a closed loop economy. Springer: Dordrecht, 2006; p xvi, 607 p. 2. Graedel, T. E.; Barr, R.; Chandler, C.; Chase, T.; Choi, J.; Christoffersen, L.; Friedlander, E.; Henly, C.; Jun, C.; Nassar, N. T.; Schechner, D.; Warren, S.; Yang, M. Y.; Zhu, C., Methodology of metal criticality determination. Environ. Sci. Technol. 2012, 46, (2), 1063-70; DOI 10.1021/es203534z. 3. Nuss, P.; Harper, E. M.; Nassar, N. T.; Reck, B. K.; Graedel, T. E., Criticality of Iron and Its Principal Alloying Elements. Environ. Sci. Technol. 2014, 48, (7), 4171-4177; DOI 10.1021/Es405044w. 4. Graedel, T. E.; Harper, E. M.; Nassar, N. T.; Nuss, P.; Reck, B. K., Criticality of metals and metalloids. PNAS 2015, 112, (14), 4257-62; DOI 10.1073/pnas.1500415112. 5. National Research Council, Minerals, Critical Minerals, and the U.S. Economy. The National Academies Press: Washington, DC, 2008; p 264. 6. European Commission (EC) Report on Critical Raw Materials for the EU; EC: Brussels, Belgium,, 2014. 7. Nassar, N. T.; Graedel, T. E.; Harper, E. M., By-product metals are technologically essential but have problematic supply. Sci. Adv. 2015, 1, (3), e1400180; DOI 10.1126/sciadv.1400180. 8. Nuss, P.; Eckelman, M. J., Life cycle assessment of metals: a scientific synthesis. Plos One 2014, 9, (7), e101298; DOI 10.1371/journal.pone.0101298. 9. Graedel, T. E.; Harper, E. M.; Nassar, N. T.; Reck, B. K., On the materials basis of modern society. PNAS 2015, 112, (20), 6295-6300; DOI 10.1073/pnas.1312752110. 10. United Nations Environment Programme (UNEP). Recycling Rate of Metals-A Status Report, A report of the Working Group on the Global Metal Flows to the International Resource Panel; 2011. 11. Nakamura, S.; Kondo, Y.; Matsubae, K.; Nakajima, K.; Tasaki, T.; Nagasaka, T., Quality- and dilution losses in the recycling of ferrous materials from end-of-life passenger cars: input-output analysis under explicit consideration of scrap quality. Environ. Sci. Technol. 2012, 46, (17), 9266-73; DOI 10.1021/es3013529. 12. Nakajima, K.; Takeda, O.; Miki, T.; Matsubae, K.; Nagasaka, T., Thermodynamic Analysis for the Controllability of Elements in the Recycling Process of Metals. Environ. Sci. Technol. 2011, 45, (11), 4929-4936; DOI 10.1021/Es104231n. 13. Amini, S. H.; Remmerswaal, J. A. M.; Castro, M. B.; Reuter, M. A., Quantifying the quality loss and resource efficiency of recycling by means of exergy analysis. J. Clean. Prod. 2007, 15, (10), 907-913; DOI 10.1016/j.jclepro.2006.01.010. 14. Nakajima, K.; Takeda, O.; Miki, T.; Matsubae, K.; Nakamura, S.; Nagasaka, T., Thermodynamic Analysis of Contamination by Alloying Elements in Aluminum Recycling. Environ. Sci. Technol. 2010, 44, (14), 5594-5600; DOI 10.1021/Es9038769. 15. Hiraki, T.; Miki, T.; Nakajima, K.; Matsubae, K.; Nakamura, S.; Nagasaka, T., Thermodynamic Analysis for the Refining Ability of Salt Flux for Aluminum Recycling. Materials 2014, 7, (8), 5543-5553; DOI 10.3390/Ma7085543. 16. Nakajima, K.; Matsubae-Yokoyama, K.; Nakamura, S.; Itoh, S.; Nagasaka, T., Substance Flow Analysis of Zinc Associated with Iron and Steel Cycle in Japan, and Environmental Assessment of EAF Dust Recycling Process. ISIJ Int. 2008, 48, (10), 1478-1483; 30 ACS Paragon Plus Environment

Page 31 of 32

Environmental Science & Technology

DOI 10.2355/isijinternational.48.1478. 17. Ohno, H.; Matsubae, K.; Nakajima, K.; Nakamura, S.; Nagasaka, T., Unintentional Flow of Alloying Elements in Steel during Recycling of End-of-Life Vehicles. J. Ind. Ecol. 2014, 18, (2), 242-253; DOI 10.1111/jiec.12095. 18. Gaustad, G.; Olivetti, E.; Kirchain, R., Toward Sustainable Material Usage: Evaluating the Importance of Market Motivated Agency in Modeling Material Flows. Environ. Sci. Technol. 2011, 45, (9), 4110-4117; DOI 10.1021/Es103508u. 19. Lovik, A. N.; Modaresi, R.; Muller, D. B., Long-term strategies for increased recycling of automotive aluminum and its alloying elements. Environ. Sci. Technol. 2014, 48, (8), 4257-65; DOI 10.1021/es405604g. 20. Miller, R. E.; Blair, P. D., Input-output analysis : foundations and extensions. 2nd ed.; Cambridge University Press: Cambridge England ; New York, 2009; p xxxii, 750 p. 21. U.S. Bureau of Economic Analysis. Input-Output Accounts Data. http://www.bea.gov/industry/io_annual.htm (accessed June 3, 2014), 22. Nakamura, S.; Nakajima, K., Waste input-output material flow analysis of metals in the Japanese economy. Mater. Trans. 2005, 46, (12), 2550-2553; DOI 10.2320/matertrans.46.2550. 23. Nakamura, S.; Nakajima, K.; Kondo, Y.; Nagasaka, T., The waste input-output approach to materials flow analysis - Concepts and application to base metals. J. Ind. Ecol. 2007, 11, (4), 50-63; DOI 10.1162/jiec.2007.1290. 24. Chen, W. Q.; Graedel, T. E.; Nuss, P.; Ohno, H., Building the Material Flow Network of Aluminum in the 2007 U.S. Economy. Environ. Sci. Technol. 2016, Submitting. 25. Nakamura, S.; Nakajima, K.; Yoshizawa, Y.; Matsubae-Yokoyama, K.; Nagasaka, T., Analyzing Polyvinyl Chloride in Japan With the Waste Input-Output Material Flow Analysis Model. J. Ind. Ecol. 2009, 13, (5), 706-717; DOI 10.1111/j.1530-9290.2009.00153.x. 26. Nakajima, K.; Ohno, H.; Kondo, Y.; Matsubae, K.; Takeda, O.; Miki, T.; Nakamura, S.; Nagasaka, T., Simultaneous material flow analysis of nickel, chromium, and molybdenum used in alloy steel by means of input-output analysis. Environ. Sci. Technol. 2013, 47, (9), 4653-60; DOI 10.1021/es3043559. 27. Ohno, H.; Matsubae, K.; Nakajima, K.; Kondo, Y.; Nakamura, S.; Nagasaka, T., Toward the efficient recycling of alloying elements from end of life vehicle steel scrap. Resour. Conserv. Recycl. 2015, 100, (0), 11-20; DOI 10.1016/j.resconrec.2015.04.001. 28. U.S. Geological Survey (USGS), 2007 Minerals Yearbook. USGS: Virginia, 2008. 29. Nakamura, S.; Kondo, Y.; Matsubae, K.; Nakajima, K.; Nagasaka, T., UPIOM: A New Tool of MFA and Its Application to the Flow of Iron and Steel Associated with Car Production. Environ. Sci. Technol. 2011, 45, (3), 1114-1120; DOI 10.1021/Es1024299. 30. Newman, M. E. J., Networks : an introduction. Oxford University Press: Oxford ; New York, 2010; p xi, 772 p. 31. Nuss, P.; Chen, W. Q.; Ohno, H.; Graedel, T. E., Structural Investivation of Aluminum in the US Economy using Network Analysis. Environ. Sci. Technol. 2016, Submitting. 32. Ministry of Industrial Affairs and Communications (Japan) 2005 input-output table for Japan; Ministy of Industrial Affairs and Communications: Tokyo, 2009. 33. Ministry of Economy, Trade and Industry (Japan),; Yearbook of Iron and Steel, Non-ferrous Metals, and Fabricated Metals Statistics; Ministry of Economy, Trade and Industry: Japan, 2006. 34. Weisz, H.; Duchin, F., Physical and monetary input-output analysis: What makes the difference? Ecolog. Econ. 2006, 57, (3), 534-541; DOI 10.1016/j.ecolecon.2005.05.011. 31 ACS Paragon Plus Environment

Environmental Science & Technology

35. Simpson, D.; Tsukui, J., The Fundamental Structure of Input-Output Tables, An International Comparison. Rev. Econ. Statist. 1965, 47, (4), 434-446; DOI 10.2307/1927773. 36. Kondo, Y., Triangulation of Input–Output Tables Based on Mixed Integer Programs for Inter-temporal and Inter-regional Comparison of Production Structures. J. Econ. Struct. 2014, 3, (1), 2; DOI 10.1186/2193-2409-3-2. 37. Kagawa, S.; Okamoto, S.; Suh, S.; Kondo, Y.; Nansai, K., Finding environmentally important industry clusters: Multiway cut approach using nonnegative matrix factorization. Soc. Networks 2013, 35, (3), 423-438; DOI 10.1016/j.socnet.2013.04.009. 38. Okamoto, S., Analyzing instability of industrial clustering techniques. Environ. Econ. Pol. Stud. 2015, 17, (3), 389-406; DOI 10.1007/s10018-014-0086-x.

32 ACS Paragon Plus Environment

Page 32 of 32