Virtual CO2 Emission Flows in the Global Electricity Trade Network

May 8, 2018 - Virtual CO2 Emission Flows in the Global Electricity Trade Network. Shen Qu† ... *Phone: +1-734-763-8644; fax: +1-734-936-2195; e-mail...
0 downloads 0 Views 1MB Size
Subscriber access provided by UNIV OF NEWCASTLE

Energy and the Environment

Virtual CO2 emission flows in the global electricity trade network Shen Qu, Yun Li, Sai Liang, Jia-Hai Yuan, and Ming Xu Environ. Sci. Technol., Just Accepted Manuscript • DOI: 10.1021/acs.est.7b05191 • Publication Date (Web): 08 May 2018 Downloaded from http://pubs.acs.org on May 8, 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 22

Environmental Science & Technology

1

Virtual CO2 Emission Flows in the Global Electricity Trade Network

2

Shen Qu a, Yun Li b,a, Sai Liang c, Jiahai Yuan b, Ming Xu a,d,*

3 4

a

5 6

b

7 8

c

9 10

d

School for Environment and Sustainability, University of Michigan, Ann Arbor, Michigan 48109-1041, United States School of Economics and Management, North China Electric Power University, Beijing 102206, People’s Republic of China State Key Joint Laboratory of Environmental Simulation and Pollution Control, School of Environment, Beijing Normal University, Beijing, 100875, People’s Republic of China Department of Civil and Environmental Engineering, University of Michigan, Ann Arbor, Michigan 48109-2125, United States

11 12

*corresponding authors.

13

Phone: +1-734-763-8644; fax: +1-734-936-2195; e-mail: [email protected] (Ming Xu).

1 ACS Paragon Plus Environment

Environmental Science & Technology

Page 2 of 22

14

Abstract

15 16 17 18 19 20 21 22 23 24 25 26

Quantifying greenhouse gas emissions due to electricity consumption is crucial for climate mitigation in the electric power sector. Current practices primarily use production-based emission factors to quantify emissions for electricity consumption, assuming production and consumption of electricity take place within the same region. The increasingly intensified crossborder electricity trade complicates the accounting for emissions of electricity consumption. This study employs a network approach to account for the flows in the whole electricity trade network to estimate CO2 emissions of electricity consumption for 137 major countries/regions in 2014. Results show that in some countries, especially those in Europe and Southern Africa, the impacts of electricity trade on the estimation of emission factors and embodied emissions are significant. The changes made to emission factors by considering inter-grid electricity trade can have significant implications for emission accounting and climate mitigation when multiplied by total electricity consumption of the corresponding countries/regions.

27 28

TOC

29 30

2 ACS Paragon Plus Environment

Page 3 of 22

Environmental Science & Technology

31

1. Introduction

32 33 34 35 36 37 38 39 40 41

Electric power generation contributes significantly to global greenhouse gas (GHG) emissions. In 2014, over 40% of the global carbon dioxide (CO2) emissions were from the electric power sector 1. Mitigation initiatives, strategies, and policies related to the power sector have taken place at various scales, including the national, regional, organizational, and even individual scales. Underpinning such effects is the accurate and fair accounting for GHG emissions, for emissions from both electricity generation and consumption. Consumption-based accounting is particularly relevant to initiatives and policies at regional, organizational, and individual levels. Converting grid electricity consumption (or purchased electricity) into emissions from power generation requires the measurement and use of the emission factor, which is defined as the emission generated due to unitary electricity consumption.

42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61

Current practices mostly use production-based emission factors for estimating emissions driven by electricity consumption. For example, the Greenhouse Gas Protocol requires organizations to report Scope 2 emissions due to grid electricity consumption by using emission factors based on national electricity generation 2, 3, except for US where emission factors for eGRID subregions should be used 4. However, as electricity is purchased and consumed from interconnected grids, production-based emission factors lead to inaccurate measurements of emissions due to electricity consumption. In particular, power grids are connected and interdependent at the regional and even the global level. Indeed, global electricity trade has been steadily increasing in past decades. For example, electricity exports and imports of OECD countries have been growing by 4.5% and 4.3% annually from 1974, and reached 511 TWh and 510 TWh in 2015, respectively 5. Electricity trade brings about economic benefit, since it opens the opportunities to exploit region variations in natural resources, climate and load timing, reducing the surplus generation capacity needed 5, 6. However, similar to the fact that globalized supply chains distance production and consumption and render environmental responsibilities more “invisible” 7, 8, cross-border electricity trade furthers the separation between electricity generation and consumption. Such separation subsequently complicates the accounting for emissions due to electricity consumption, since consumed electricity may be generated by multiple regional grids where electricity generation may rely on specific energy sources, use various technologies, and produce different amounts of emissions for unitary electricity generated.

62 63 64 65 66 67 68 69

Given the complexity of interconnected grid and increasingly intensified electricity trade, a proper accounting framework for cross-border electricity trade is needed to correctly measure the emissions due to electricity consumption. Without such a framework, accounting exercises may lead to wrong estimations of GHG emissions by consuming entities, and even consequently undermine the promised or expected effects of mitigation efforts due to carbon leakage in the power sector 9-11. Data disclosure initiatives have begun to report emissions from power plants worldwide 12, 13. However, without controlling for emissions from electricity trade, countries/regions may well stabilize or even reduce direct emissions from their power sector 3 ACS Paragon Plus Environment

Environmental Science & Technology

Page 4 of 22

70 71

while importing embodied emission from other nations, and thus compromise the effectiveness of initiatives directly targeting electricity production.

72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87

In this study, we explore the impacts of cross-border electricity trade on country-/region- level average CO2 emission factors of electricity consumption. We include 137 countries/regions with complete data on electricity generation, cross-border electricity transfers, fuel mixes and emission factors for electricity generation. Over 90% of global electricity trade in 2014 occurred among these countries/regions. We employ and compare three different methods to estimate emission factors of electricity consumption (which will be detailed in the next section), including a network-based approach for estimating embodied emission flows in electricity trade networks 14-16 which, to our knowledge, has never been employed at the global scale before this study. Therefore, we are able to identify where and how electricity trade may change embodied CO2 emissions and fuel mix of electricity consumption. Even though for most world regions the effects of cross-border electricity trade on emission factors remain relatively marginal, the impacts can be significant when combined with the amounts of electricity consumed. In addition, we evaluate the virtual emission transfers across borders due to electricity trade and identify top net emission importers and exporters, which are hotspots for collective efforts among countries/regions to mitigate carbon leakage occurring in the power sector.

88 89

2. Method

90 91 92 93 94 95 96 97 98 99 100 101 102

Although the emission factor of electricity consumption is conceptually analogous to consumption-based emission accounting, there is a fundamental difference between exchanges in electricity and exchanges in other goods and services, which render the quantification of emission factors of electricity consumption unique and nontrivial 17. For example, at a particular moment, if grid A only imports electricity from grid B which in turn imports from grid C, in effect grid A is importing a blend of electricity (or electrical energy) from both B and C. In other words, there is a virtual emission flow generated in grid C due to electricity consumption in grid A, although there is no direct electricity transfer from grid C to A. The virtual flows of emissions should follow the pattern of inter-grid electricity transfers. Therefore, when n grids are generally interconnected and transfer electricity among themselves, to correctly estimate the emission factor of electricity consumption from a particular grid, ideally one should include all the n grids in the model, instead of only considering its direct trade partners. This becomes more crucial as inter-grid electricity trade intensifies.

103 104 105 106 107 108

Table 1 summarizes the methodologies, resolutions, and scales of previous research on embodied emissions in electricity consumption. We exclude studies that simply base their estimates on electricity generation within the considered geographical boundary and neglect electricity trade/transfers. There are mainly two methods to account for electricity trade. The first method considers only a grid’s direct electricity trade partners to quantify its emission factor (i.e., direct trade-adjustment approach). Embodied emissions in electricity consumption 4 ACS Paragon Plus Environment

Page 5 of 22

109 110 111 112 113 114 115 116

Environmental Science & Technology

in a grid is calculated as the emissions directly from electricity generation minus the emissions for producing exported electricity plus the emissions (generated outside the grid) for producing the imported electricity, with the (partly incorrect) assumption that the electricity a grid exports is entirely generated in this grid. The second method (i.e., network-based approach), improves upon the first one, enlarging the geographic boundary when evaluating every single grid’s embodied emissions to reflect all virtual emission flows in the electricity trade network. Recently we developed this method with extensive discussion on its relation to previous approaches 17, and applied it to China’s interprovincial electricity trade network 18.

117 118 119

Table 1. Studies that report embodied emissions in electricity consumption and take account of inter-grid electricity trade. Study

Method to account for electricity trade

Geographic coverage and resolution

Boundary for calculating embodied emissions for a grid

Marriott and Matthews (2005)19

Optimization

US and US states

States relatively close to the considered state

Kang et al. (2012)14

Network-based

China and subnational grids

Entire China

Song et al. (2013)20

Direct tradeadjustment

China and subnational grids

Sub-nation grids with direct electricity transfers

Lindner et al. (2013)21

Direct tradeadjustment

China and Provinces with direct Chinese provinces electricity transfers

Ji et al. (2015)16

Network-based, with only net electricity transfers

Eurasian continent and nations

Entire electricity trade network

Kodra et al. (2015)15

Network-based

US and US power control areas

Entire US

Zafirakis et al. (2015)22

Direct tradeadjustment

Europe and nations

Nations with direct electricity transfers

Colett et al. (2016)23

Nested approach

US and US power control areas

NERC regions surrounding the power control area

Qu et al. (2017)18

Network-based

China and Entire China Chinese provinces

120 121 122

Besides the above two approach, other approaches in Table 1 are not directly relevant to our study. Marriott and Matthews (2005) use an optimization method to estimate interstate 5 ACS Paragon Plus Environment

Environmental Science & Technology

Page 6 of 22

123 124 125 126 127

electricity transfers in the US without real-world data. Colett et al. (2016) adjusts for the electricity transfers between a power control area of the US and the surrounding North American Electric Reliability Corporation (NERC) regions to estimate GHG emission factors, considering two protocols based respectively on bilateral and net electricity trade which are essential based on the direct trade-adjustment method.

128 129 130 131 132 133 134 135 136 137 138 139 140 141 142

It is also worth mentioning that the electricity market has been systematically modelled in sophisticated life cycle inventory (LCI) databases such as ecoinvent 24. In particular, the supply mix of electricity to consumers is modelled as domestic production plus imports, according to Itten et al. (2012) 25. Conceptually, electricity imports should have the same composition as the supply mixes of the exporting countries, but in current practice of LCI database compilation , production mixes of the exporting countries were used for simplification, to avoid “virtual feedback loops between the electricity supply mixes of different countries” 25. However, modeling such feedback loops is not only theoretically valid, but may prove necessary as electricity trade intensifies (e.g., see ref 15). Recent methodological contribution by us and coauthors 17 has simplified and standardized this procedure, and clarified both its similarity to and difference from the Economic Input-Output (EIO) model. In the Discussion section, we will calculate the fuel compositions of electricity supply in different countries with the methods mentioned above. Although this paper focuses on CO2 emissions, the inferred changes in fuel composition (e.g. more hydropower and less fossil fuels) could have important life cycle implications.

143 144 145 146 147 148 149 150

The interconnected grid network. In this study, a national/regional grid is modelled as node. In the considered period (i.e., a year), we define total electricity inflow to the node and total electricity outflow from the node. The former equals to total electricity generation in the country/region plus total electricity import. The latter is total electricity consumption in the country/region plus total electricity export. These two terms are equal by conservation of energy, as in equation (1): n

n

j =1

j =1

xi = pi + ∑Tji =ci + ∑Tij

(1)

151 152

where xi denotes the total electricity inflow or outflow of grid i, pi is electricity generation, ci is electricity consumption, and Tij is the annual electricity transfer from grid i to grid j.

153 154 155 156 157

Based only on annual data on grid electricity generation, consumption, and inter-grid electricity transfers among the 137 countries, it is not possible to differentiate components of electricity inflow (i.e., imports from other grids and generation) by their destination (i.e., export to other grids or used for own consumption). We assume that, for each node, electricity generation and imports are first mixed, and then either consumed or exported to other grids.

6 ACS Paragon Plus Environment

Page 7 of 22

Environmental Science & Technology

158 159 160

Next we use a matrix representation of the electricity supply chains in an n-grid electricity trade network to compute the emissions embodied in each grid’s electricity consumption. We define the direct outflow matrix, B, and total outflow matrix, G, in equations (2) and (3).

161

B = xˆ −1T (2)

162

G = ( I − B)−1 = I + B + B2 +L (3)

163

where xˆ is the diagonal matrix with total electricity flow xi on the diagonal, and T is the inter-

164

grid electricity transfer matrix with T ij as the (i,j)th element. The right-hand-side of equation (3)

165 166

invokes the Taylor expansion of ( I − B)−1 , which is a commonly used technique in input-output economics 26.

167 168 169 170 171

In this accounting framework, the (i,j)th element of total outflow matrix G is the amount in the total electricity inflow to grid j that comes from unitary electricity generation in grid i. This transfer is enabled through infinite electricity supply chains, including self-supply (represented by the term I), direct transfers (B), transfers through one intermediate node (B2), transfers through two intermediate nodes (B3), and so on.

172 173 174 175

Note a known fraction of the total electricity inflow to a grid is consumed without further exporting. Based on the above results, we can immediately derive the generation-consumption matrix, H, where the (i,j)th element is the electricity consumed in grid j that corresponds to the unitary electricity generation in grid i:

176

ˆ ˆ −1 (4) H = Gcx

177 178 179

where cˆ is the diagonal matrix that contains electricity consumption of each grid. Since virtual emission flows are concurrent with electricity transfers, the (i,j)th element of H also represents the allocation of grid i’s production-based emissions to grid j’s electricity consumption.

180

184

Virtual emission flows and emission factors. Now we are ready to link the embodied emissions of electricity consumption in a particular grid to direct emissions of each grid in the network, and therefore to compute the emission factors of electricity consumption. Given emission factors of electricity generation efi G for each national/regional grid i, the total direct emission

185

from each grid is eiG = ef i G pi . We can now compute the matrix of inter-grid virtual emission

186 188

flows E C , the row vector of embodied emissions in grid electricity consumption eC , and the corresponding vector of network-based emission factors of electricity consumption ef C , Network , respectively, with equations (5)-(7).

189

EC = EG H (5)

190

eC = [1,L,1] E C (6)

181 182 183

187

7 ACS Paragon Plus Environment

Environmental Science & Technology

Page 8 of 22

191

ef C , Network = eC cˆ −1 (7)

192

where E G is diagonal matrix with eiG on the diagonal, E C is an n by n matrix with the (i,j)th

193 194

element being the direct and indirect emission flow from i to j, eC is the 1 by n vector with the ith element being the embodied emissions in country/region i’s electricity consumption, and

195

ef C , Network is the associated 1 by n vector containing emission factors for electricity consumption.

196 197 198 199 200 201 202 203

Adjustment for open systems. In practice, the n grids in the model may not constitute a closed system, because of incomplete data for all countries/regions in the global electricity trade network (see section 3.3). Some countries/regions, therefore, have to be left out, but they may still export electricity to or import electricity from the system in the real world. Therefore, the method presented in the previous section needs slight adjustments to reflect the open nature of the system in practice. In particular, we can keep all the procedures the same but simply rewrite equations (1) and (5) as follows.

204

xi = pi + ∑Tji + ∑T%ji =ci + ∑Tij + ∑T%ij

n

mi

n

li

j =1

j =1

j =1

j =1

205

E C = ( E G + E In ) H (5’)

206

where

(1’)

mi

∑ T%ji denote the mi electricity transfers to grid i from outside of the system, j =1

207 208 209

li

∑ T%

ij

j =1

denote the li electricity transfers from grid i to outside, and E In is the diagonalized matrix where each diagonal element is the emission generated from producing electricity that is imported into the corresponding country in the n-grid system from outside the system.

210 211 212

Direct trade-adjustment approach. The following equation lays out the direct trade-adjustment approach, which we used to compute emission factors of purchase electricity for comparison.

213

  ef i C , DirectAjust =  ef i G ⋅ pi + ∑ ef jG ⋅ T ji − ∑ ef i G ⋅ Tij  ci j ≠i j ≠i  

214

As described previously, the estimation of ef i C , DirectAjust is only based on its own emission factors

215

of generation, efi G , and values of ef jG of those direct electricity trade partners with grid i.

(8)

216

8 ACS Paragon Plus Environment

Page 9 of 22

Environmental Science & Technology

217

3. Data

218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237

Electricity trade. Global electricity trade data (in kWh) in 2014 are reported in the United Nations (UN) Comtrade Database 27. To address statistical discrepancies, we cleaned the data with the following procedures. First, if only export or import data exist between an exporterimporter pair, we simply use the existing data. Second, if both export and import data are reported but not equal, we use the quantity which is greater. This is because some might conceal trade volumes to evade regulations such as tariffs. Alternatively, one may argue that the import data should be used, since they are recorded with greater accuracy for regulation purposes. Therefore, we also calculate the results with import data if inconsistency occurs. The effect on emission factors are discussed in Section 4. Third, when the physical electricity trade volumes are unrealistically high, such as Ecuador importing 2.2×108MWh from Colombia, we adjusted these data and report the adjustment in Table S1. The principle is to rely on trade volumes reported in US dollars to correct data in kWh. Given the possibility of typos when physical data are entered, we multiply the unrealistic data by 10-n such that, after the adjustment, the electricity price falls in the reasonable range as reported by the International Energy Agency 5. Lastly, we delete electricity trade records that are not within the reasonable geographical ranges (Table S2). For example, it was reported that Japan sent 7,050 MWh electricity to the United Arab Emirates in 2014, which is possibly due to electricity purchased by cargo ships and evidently does not conform to the concept of electricity trade in our model. We simply delete such records. All the adjustments we made to the original UN Comtrade data are described in Table S1 and S2.

238 239 240

Furthermore, all electricity not transferred out of a country, including the extra amount of electricity to enable the transmission to another country, is assumed to be consumed in the country. In this way, we implicitly assign the transmission losses to the generating country 15.

241 242 243 244 245 246 247

Electricity generation and direct emissions. Average CO2 emission factors of electricity generation at the country level in 2014 are from IEA statistics 1. IEA has estimated these values according to the “Tier 1” approach in IPCC guidelines 28, using country-specific net calorific values for various types of fuels 1. These statistics provide the most comprehensive coverage, but could be substituted with more accurate data if available. National/regional electricity generation and fuel mixes divided into six categories (fossil fuels, nuclear, geothermal, hydro, solar/wind, and biofuels) for electricity generation are also from IEA 5.

248 249 250 251

Countries/regions included in this study. There are 137 countries/regions (listed in Table S4) for which data are complete for this study. The electricity trade among them comprises about 91% of total global electricity trade recorded in the UN Comtrade database. We include these 137 countries/regions in our final analysis for emission factors and virtual emission flows.

252 253 254 255

Some countries that are left out in this study may export electricity to and import electricity from the countries/regions included in this study, as reported by the UN Comtrade database. The quantities of these electricity transfers are not significant. For each of these outside countries, we made reasonable assumptions about its fuel mix and emission factor for 9 ACS Paragon Plus Environment

Environmental Science & Technology

256 257 258

Page 10 of 22

electricity generation (Table S3). For example, Laos, which is not among the 137 countries/regions, is from the “other Asia” category in IEA statistics 5. We then use the IEA data for “other Asia” on fuel mix and emission factor of electricity generation in equation (5’).

259 260

4. Results

261

Global electricity trade network in 2014

262 263 264 265 266 267 268

Figure 1 visualizes the global electricity trade network in 2014, where nodes represent electricity generation of the country/region and clockwise curves represent direct electricity transfers from origins to destinations. As previously revealed for the year 2011 29, there are four major communities in the global electricity trade network, namely, Eurasian, North and African, Central American, and South American. The insets of Figure 1A and 1B zoom in the African community and the European part of the Eurasian community, where electricity trade is more intensive and the network structure is more complex.

269 270 271 272 273

Obviously, limitations of transmission technology dictate the community structure of the crossborder electricity trade. National/regional policies, resource endowments, and generation technologies are also important factors underlying the structure of the global electricity trade network. This structure leads to the differences of emission factor estimates and patterns of virtual emission flows among countries/regions.

274

10 ACS Paragon Plus Environment

Page 11 of 22

Environmental Science & Technology

(A) 275 276 277 278 279 280

(B)

Figure 1. Global electricity trade network in 2014. Curves represent direct electricity transfers, which flow clockwise from origins to destinations. Curve widths indicate the amount electricity transfers, some of which are labeled in TWh. Node size represents electricity generation of countries/regions. Countries/regions are labeled with their ISO codes. The insets of (A) and (B) zoom in the African community and the European part of the Eurasian community respectively. Map images are from the GADM database of Global Administrative Areas 30.

281

11 ACS Paragon Plus Environment

Environmental Science & Technology

Page 12 of 22

282

CO2 emission factors of electricity consumption in countries/regions

283 284 285 286 287 288 289 290 291 292 293

Variation of CO2 emission factors of electricity consumption. At the global level, there is a wide range of CO2 emission factors of electricity consumption across countries/regions. Figure 2A illustrates the global distribution of emission factors of electricity consumption for each of the 137 countries/regions in our dataset with shades of color. For comparison, we also tabulate the specific estimates with different methods in Table S4. The places where consuming unitary electricity from the local grid is responsible for the most global CO2 emissions are found in East and Southeast Asia, Middle East, and Southern Africa. Emission factors of electricity consumption are highest in Iraq (1,169 g/kWh), Mongolia (1,148 g/kWh), South Africa (987 g/kWh), Turkmenistan (890 g/kWh), and India (813 g/kWh), due to the dominant uses of carbon-intensive fuels in their own electricity generation and/or their foreign electricity suppliers.

294 295 296 297 298 299 300 301 302 303 304 305 306 307 308

On the other hand, some countries have very low emission factors of electricity consumption mainly because of the energy composition for their electricity supply. Here the focus is on CO2 emissions from fossil fuel consumption, and the entire life cycle emissions of the generation technologies are beyond the scope of this study. Although it only offers an incomplete picture of the life cycle environmental impacts of electricity, our method can help improve the accuracy of a full LCA study with improved estimation of life cycle emissions from electricity consumption. In addition, our method is also consistent with some widely used standards such as the Scope 2 emission of the Greenhouse Gas Protocol 2. Iceland has zero CO2 emission from fossil fuels since it solely relies on hydro and geothermal energy (Table S4, Figure S1). Emission factors are also low in Ethiopia (1 g/kWh), Tajikistan (7 g/kWh), DR Congo (7 g/kWh), Norway (12 g/kWh), and Sweden (18 g/kWh). These countries mainly depend on hydropower, except for Sweden which consumes electricity from a mix of nuclear, hydro, and wind power and biofuels. In Section 5, we discuss how modelling cross-border electricity trade differently may change the fuel mixes of countries’ electricity consumption, which could have implications for relevant life cycle assessments.

309

12 ACS Paragon Plus Environment

Page 13 of 22

310 311

Environmental Science & Technology

Figure 2. (A) Emission factors of electricity consumption of countries/regions. Shades of color represent the values of emission factors ( efi C , Network for country/region i), while sizes of circles

315

represent total electricity consumption. White areas are where data are incomplete. (B) Differences between various accounting methods for country/region-level emission factors of electricity consumption. Blue area indicates that emission factors for electricity consumption ( efi C , Network ) is smaller than those for generation ( efi G ); red area indicates the opposite. For

316 317

hashed/meshed areas, enlarging the system boundary of accounting for electricity trade significantly changes estimates of emission factors, either downward (in hashed areas, where

312 313 314

13 ACS Paragon Plus Environment

Environmental Science & Technology

Page 14 of 22

318

efi C ,Network efi C ,DirectAdjust 105%).

319

Map images are from the GADM database of Global Administrative Areas 30.

320 321 322 323 324 325 326 327 328 329

Impacts of electricity trade on emission factors. Figure 2B illustrates the impacts of crossborder electricity trade on emission factors of electricity consumption. Cross-border electricity trade has the greatest impacts in emission factors of Mongolia, Argentina, and clusters of countries in Europe and Southern Africa, where cross-border electricity transfers are intensive and fuel mix differences among countries/regions are significant. For most countries/regions, the effect of indirect electricity transfers on emission factors is less significant than the effect of direct electricity transfers. But our results also show that, in some countries in Europe and Northern Africa, the former effect is already an important factor to consider when estimating emissions from electricity consumption.

330 331 332 333 334 335 336 337 338 339 340 341 342 343 344 345 346 347 348 349 350 351 352 353 354

Figure S2 plots specific values of emission factors of electricity consumption, for countries/regions where cross-border trade can change the results by at least 5%. We compare emission factors calculated based on three different methods, including the production-based approach, direct trade-adjustment approach, and network-based approach. In a couple of countries, cross-border electricity transfers already have important implications for emission factor estimates. Generally, if only adjusting for direct electricity trade raises (lowers) a country’s emission factor, further accounting for indirect electricity transfers with the networkbased approach tends to lower (raise) the estimate, although the change is usually smaller than the change from adjusting for direct trade. Take Estonia and Macedonia as examples. Their emission factors based on the network-based approach are respectively 81% and 86% of emission factors based on the production-based approach, and 120% and 118% of direct tradeadjustment emission factors. For Botswana and Zimbabwe, the emission factors from the network-based approach are respectively 50% and 92% of emission factors from the production-based approach, and 117% and 111% of direct trade-adjustment emission factors. Changes in the opposite direction occur for countries like Switzerland and Luxembourg, where the emission factors estimated with the network-based approach are 424% and 125% of the production-based emission factors, and 74% and 70% of direct trade-adjustment emission factors. For Mozambique and Namibia, emission factors of the network-based approach are respectively 920% and 4,441% of the production-based emission factors, and 61% and 3% of direct trade-adjustment emission factors. The drastic changes observed for Mozambique and Namibia are due to their pass-through roles in the cross-border electricity trade network: because the amounts of electricity imports and exports are much greater than local electricity generation and consumption, they serve as relaying nodes in the trade network. Therefore, it is virtually impractical to estimate their emission factors simply based on electricity generation or direct electricity trade.

14 ACS Paragon Plus Environment

Page 15 of 22

Environmental Science & Technology

355 356 357 358 359 360 361

Inconsistencies between electricity export and import data may also have significant impacts on emission factors estimations. Table S5 illustrate the changes of such estimations if we use electricity import data when export and import do not match. In 7 of the 137 countries, the results change by more than 5%. Also, for these countries, no matter what electricity trade data is used, emission factors for electricity consumption diverge significantly from those for production. Thus, greater international efforts to collect electricity trade data are needed for more accurate estimation of carbon footprint of electricity consumption.

362 363 364 365 366 367 368 369

Two general patterns in our results of different emission factor estimations are worth noting. First, more notable impacts of cross-border electricity trade occur in smaller countries. This is due to the fact that large countries are usually endowed with abundant and diverse energy sources and less reliant on electricity trade. However, if data on inter-grid electricity transfers in a large country are available with finer geographic resolutions, advanced methods that are able to properly account for direct and indirect electricity trade, such as the network-based approach used in this study, will be indispensable for understating the environmental impacts of electricity consumption (e.g., see Kodra et al.15 for inter-PCA emission flows in the US).

370 371 372 373 374 375 376 377 378 379

Second, as noted above, incorporating the effects of indirect electricity transfers tends to pull the emission factors back to the levels implied by local electricity generation. In retrospect, this is due to the general pattern that if a country/region imports electricity from another one with very carbon-intensive technologies, the exporting country/region tends to import less carbonintensive electricity from other grids and indirectly transfer the electricity to the original importing country/region. This result is consistent with previous findings 19 that accounting for inter-state electricity trade in the US tends to push the state emission factors to the country average level. The selection of geographic boundaries for emission factor estimation thus has important implications for the evaluation of emission responsibility on the consumer side, which we will further discuss in section 5.

380 381 382 383 384 385 386 387

Embodied CO2 emission flows in electricity trade. Cross-border electricity trade can lead to shifts in emission responsibilities or carbon leakages, whereby countries/regions may well decarbonize electricity generation inside their borders but at the same time import carbonintensive electricity. Among the countries included in this study, 177.9 Mt of CO2 emissions from electricity generation are embedded in cross-border electricity trade in 2014. This amount lies between the total CO2 emissions of the Netherlands (158 Mt) and Argentina (199 Mt) 31, and may increase in the future due to intensifying electricity trade.

388 389 390 391

Figure 3 and Table S6 show net virtual CO2 imports/exports from electricity trade of the studied countries/regions. The clusters of net virtual emission transfers mainly reside in Europe, Southern Africa, and North America (from Canada to the US). These virtual emission transfers coincide with those of the global electricity trade network (Figure 1). However, as the traded

15 ACS Paragon Plus Environment

Environmental Science & Technology

392 393

Page 16 of 22

electricity may be generated from different energy sources, the strengths of the specific virtual emission flows can be quite different from those simply implied by electricity trade flows.

394

395 396 397 398 399

Figure 3. Net virtual CO2 emission imports in electricity trade. Lines and Arrows indicate net flows of embodied emissions. The largest net virtual CO2 emission flows are listed in Table S6. Countries/regions are indexed with their ISO codes. Map images are from the GADM database of Global Administrative Areas 30.

400 401 402 403 404 405 406 407 408

When a country/region is a net virtual CO2 importer (or exporter) in the electricity trade network, the CO2 emission responsibility for its electricity consumption is greater (or less) than that for its electricity generation. In Europe, countries with the largest amounts of net virtual CO2 imports through electricity trade are Italy (11 Mt), Austria (7.7 Mt), Switzerland (5.5 Mt) and Hungary (5.1 Mt); and the most important net virtual CO2 exporters are Germany (32.4 Mt), Czech Republic (8.5 Mt) and Ukraine (5.4 Mt). In Africa, Mozambique and Botswana have the largest net virtual CO2 imports (5.5 Mt and 4.9 Mt, respectively), and South Africa is the most important net CO2 exporter (12.4 Mt).

409 410 411 412 413 414

For most of the above-mentioned countries, the net virtual CO2 imports through electricity trade constitute significant shares of embodied CO2 emissions of electricity consumption. For example, Austria’s net virtual CO2 imports are 45% of the total CO2 emissions embodied in its electricity consumption, while that percentage is 89% for Mozambique. These results again attest to the importance of virtual emission flows in cross-border electricity trade and consequently for environmental policies relevant to the electric power sector.

415 16 ACS Paragon Plus Environment

Page 17 of 22

Environmental Science & Technology

416

5. Discussion

417 418 419 420 421 422 423 424 425

Impacts of cross-border electricity trade on embodied CO2 emissions. In 2014, the impacts of cross-border electricity trade on emission factors of electricity consumption still remain marginal in most regions, and this is especially true for the impacts of indirect electricity transfers. However, for the key countries/regions identified in Figure 2B, virtual emission flows have become more important for CO2 emissions of their electricity consumption. Without properly accounting for such virtual emissions from electricity trade, the environmental implications of the electricity consumption can be significantly misunderstood. Moreover, electricity trade has been increasing 5, and may have even greater effects on embodied emissions embodied in the future.

426 427 428 429 430 431 432 433 434 435 436 437

For the shaded countries/regions in Figure 2B, we calculate how much CO2 emissions will be added to or deducted from their emission responsibilities due to electricity consumption, when we change the geographic boundary of estimation from the country/region itself or from the country/region with its direct electricity trade partners, to the entire electricity trade network, as shown in Figure S3. Red bars represent the effects of considering both direct and indirect trade, and green bars illustrate effects of accounting for indirect trade in addition to direct trade. Note that different methods use the same country-/region- level electricity consumption, but different emission factor estimates. Similar to the pattern we previously observed for emission factors, while directly adjusting for electricity trade can significantly change the country-/region- level embodied emissions, it may also overestimate the impact of electricity trade. Using the network-based approach generally leads to emission estimates closer to those based on country-/region- energy mixes for electricity generation.

438 439 440 441 442 443 444 445 446 447 448

At the country/region level, cross-border electricity trade has important implications for emission responsibilities of electricity consumption. Using emission factors fully accounting for cross-border electricity trade, in almost every country/region, the total embodied CO2 emissions in electricity consumption increase or decrease by at least one million tons (Figure S3, red bars). Such variation amounts to CO2 emissions of almost 150,000 homes’ annual electricity use and more than 200,000 passenger vehicles driven in one year, based on some rough estimations 32. Accounting for indirect cross-border electricity transfers in addition to direct transfers also results in increases in total CO2 emissions relate to electricity consumption over one million tons in Germany, Italy, South Africa, and Botswana, and decreases in total CO2 emissions over one million tons in Austria, Switzerland, Latvia, Slovak Republic and Mozambique (Figure S3, green bars).

449 450 451 452 453 454

Therefore, even with the current level of cross-border electricity trade, it is already important to properly account for emissions embodied in electricity trade when estimating emission factors of electricity consumption, especially for countries in Europe and Southern Africa. Moreover, global electricity trade has been intensifying over time, with the physical trade volume of OECD countries growing more than sixfold in the past four decades 5. In the future, the impacts of electricity trade on emissions accounting for the power sector can be much 17 ACS Paragon Plus Environment

Environmental Science & Technology

455 456

Page 18 of 22

more significant, and the network-based approach used in this paper will become more useful for understanding environmental implications of electricity consumption.

457 458 459 460 461 462 463 464 465 466 467

Fuel mix of electricity generation and consumption. With cross-border electricity trade, the energy sources virtually supporting a country/region’s electricity consumption may be quite different from those directly contributing to its electricity generation. As mentioned earlier, traditional LCI databases such as ecoinvent model the electricity supply to consumers in a country as coming from its generation plus imports, practically using the exporting countries’ generation mix to represent their electricity exports 24, 25. As electricity trade intensifies, even such adjustment may to a certain extent misestimate electricity fuel mix since it leaves out indirect electricity transfers. Figure S1 compares the fuel mixes of electricity generation, consumption and (production + import) (with the method in ref 25) for countries included in this study. In the following we discuss two notable patterns.

468 469 470 471 472 473 474 475 476 477

First, the fuel mix of the electricity consumption differs much from that of the electricity generation in a lot of countries/regions. For example, from the consumption side, Austria and Croatia respectively rely 18% and 24% less on hydropower, 10% and 5% more on fossil fuels, and 8% and 18% more on nuclear energy than from the production side. Due to cross-border electricity trade, consumed electricity usually has more diversified fuel mixes. For example, the fuel mix of electricity consumed by Mongolia and Botswana, inferred with the network-based approach, contains sizable shares of hydropower (4% and 36%), although their generated electricity almost entirely comes from fossil fuels (97% and 100%). Denmark generates electricity using fossil fuels, wind, and geothermal energy, but also consumes hydro and nuclear power virtually from other countries/regions.

478 479 480 481 482 483

Second, the (production + import) calculation 25 provides a good approximation for fuel mixes while for most countries/regions , in a few countries, the results are significantly different (Figure S1). The network approach reveals that Botswana consumes considerably more fossil fuels (from 51% to 62% of electricity consumption) and less hydropower (from 47% to 35%). Noticeable changes in fuel mixes also occur in Central and Eastern Europe, including Croatia, Hungary, Latvia, Montenegro and Slovenia.

484 485 486 487 488 489 490

The above differences not only can drive the CO2 emission factors for electricity consumption away from those for electricity generation, but may also have other important environmental implications since generation technologies are not innocuous to the environment 33-35. These results attest the importance of cross-border electricity trade in the electric power sector which may have rendered the environmental impacts of electricity consumption more complicated than ever. Thus appropriately and timely accounting for its effects is critical for understanding the environmental implications of electricity consumption.

491 492

It will be useful to integrate more detailed electricity data with the network method and its adaption. Note here the fuel mix for electricity consumption is derived using aggregate cross18 ACS Paragon Plus Environment

Page 19 of 22

493 494 495 496 497 498 499 500 501 502

Environmental Science & Technology

border electricity trade data which per se does not differentiate various types of energy sources to produce transferred electricity. The results are inferred from the fuel mixes of each country/region’s electricity generation and the structure of aggregate electricity trade. In reality, different types of electricity can be mixed in the transmission and distribution system. Therefore it is impractical to separate them in cross-border trade, except when there are records of specific contracts between suppliers and users 3. Recent progress in building multiregional input-output (MRIO) database have enabled separate international trade records for specific types of electricity such as hydro or solar for as late as 2007 36, 37. However, assembling MRIO databases entails a lot of assumptions and the resulted data are meant for economic systems analysis rather than to precisely represent individual trade flows.

503 504 505 506 507 508 509 510 511 512 513 514 515

Limitations and future research. The resolution of the results is at the country level. This is more satisfactory for small countries, where incorporating the impacts of electricity trade can be crucial for estimating emission factors of electricity consumption. However, for large countries such as the US and China 15, 18, regional heterogeneity within the country is relatively more important for understanding environmental impacts of electricity consumption. To better understand the environmental impacts of electricity consumption in large countries in the future research, data on electricity generation, consumption and transfers should be acquired at finer geographic resolutions within the country. These data, of course, need to be linked with cross-border electricity trade data. Furthermore, if specific contracts for purchased electricity are available, it is more appropriate to allocate emission responsibilities based on these contracts rather than on any type of grid-average emission factors. However, only a small portion of electricity can be traced in this way 2, 25.

516 517 518 519 520 521 522 523 524 525 526 527

This study quantifies average emission factors of electricity consumption from specific grids. This type of emission factor serves to estimate emissions due to the electricity consumption of existing activities, and underpins the accounting of corporate emission responsibilities. Alternatively, marginal emission factors are used to measure impacts of incremental changes of electricity demand 38, for example, the large-scale deployments of new technologies 39. There is no estimation for marginal electricity factors that incorporate the impacts of electricity trade 38, despite the obvious importance of doing so. This is primarily because estimating marginal emission factors are more data intensive. However, if data on electricity generation, consumption and transfers are available at finer temporal and geographic resolutions, statistical approaches 40 can be used together with the matrix formulation in this study to estimate marginal emission factors from electricity consumption, taking account of the role of the ever-growing electricity trade.

19 ACS Paragon Plus Environment

Environmental Science & Technology

528

Supporting Information

529

Tables and Figures supporting the main text.

Page 20 of 22

530 531

Competing financial interests

532

Authors declare no competing financial interests.

533 534

Reference

535 536 537 538 539 540 541 542 543 544 545 546 547 548 549 550 551 552 553 554 555 556 557 558 559 560 561 562 563 564 565

1. IEA, CO2 Emissions from Fuel Combustion 2016. IEA: 2016. 2. WRI, The greenhouse gas protocol: A corporate accounting and reporting standard (revised edition). 2004. 3. Sotos, M. GHG Protocol Scope 2 Guidance; World Resources Institute: 2016. 4. USEPA, The Emissions & Generation Resource Integrated Database (eGRID). In Washington DC, US: United States Environmental Protection Agency, 2017. 5. IEA, Electricity Information 2016. https://www.iea.org/statistics/relateddatabases/electricityinformation/, 2016. 6. Antweiler, W., Cross-border trade in electricity. Journal of International Economics 2016, 101, 42-51. 7. Hertwich, E. G.; Peters, G. P., Carbon Footprint of Nations: A Global, Trade-Linked Analysis. Environ Sci Technol 2009, 43, (16), 6414-6420. 8. Kanemoto, K.; Moran, D.; Hertwich, E. G., Mapping the Carbon Footprint of Nations. Environ Sci Technol 2016, 50, (19), 10512-10517. 9. Peters, G. P.; Minx, J. C.; Weber, C. L.; Edenhofer, O., Growth in emission transfers via international trade from 1990 to 2008. Proceedings of the National Academy of Sciences 2011, 108, (21), 8903-8908. 10. Malik, A.; Lan, J.; Lenzen, M., Trends in global greenhouse gas emissions from 1990 to 2010. Environ Sci Technol 2016, 50, (9), 4722-4730. 11. Peters, G. P.; Hertwich, E. G., CO2 embodied in international trade with implications for global climate policy. Environ Sci Technol 2008, 42, (5), 1401-1407. 12. CGD, Carbon Monitoring for Action. In http://carma.org/, 2017. 13. WRI, Global Power Watch. In http://globalpowerwatch.org/, 2017. 14. Kang, C.; Zhou, T.; Chen, Q.; Xu, Q.; Xia, Q.; Ji, Z., Carbon emission flow in networks. Sci Rep-Uk 2012, 2, 479. 15. Kodra, E.; Sheldon, S.; Dolen, R.; Zik, O., The North American Electric Grid as an Exchange Network: An Approach for Evaluating Energy Resource Composition and Greenhouse Gas Mitigation. Environ Sci Technol 2015, 49, 13692-13698. 16. Ji, L.; Liang, S.; Qu, S.; Zhang, Y.; Xu, M.; Jia, X.; Jia, Y.; Niu, D.; Yuan, J.; Hou, Y.; Wang, H.; Chiu, A. S. F.; Hu, X., Greenhouse gas emission factors of purchased electricity from interconnected grids. Appl Energ 2015, 184, 751-758.

20 ACS Paragon Plus Environment

Page 21 of 22

566 567 568 569 570 571 572 573 574 575 576 577 578 579 580 581 582 583 584 585 586 587 588 589 590 591 592 593 594 595 596 597 598 599 600 601 602 603 604 605 606 607

Environmental Science & Technology

17. Qu, S.; Wang, H.; Liang, S.; Shapiro, A. M.; Suh, S.; Sheldon, S.; Zik, O.; Fang, H.; Xu, M., A Quasi-Input-Output model to improve the estimation of emission factors for purchased electricity from interconnected grids. Appl Energ 2017, 200, 249-259. 18. Qu, S.; Liang, S.; Xu, M., CO2 Emissions Embodied in Interprovincial Electricity Transmissions in China. Environ Sci Technol 2017, 51, (18), 10893-10902. 19. Marriott, J.; Matthews, H. S., Environmental Effects of Interstate Power Trading on Electricity Consumption Mixes. Environ Sci Technol 2005, 39, (22), 8584-8590. 20. Song, R.; Zhu, J.; Hou, P.; Wang, H., Getting every ton of emissions right: An analysis of emission factors for purchased electricity in China. World Resources Institute, Working Paper. 2013. 21. Lindner, S.; Liu, Z.; Guan, D.; Geng, Y.; Li, X., CO2 emissions from China’s power sector at the provincial level: Consumption versus production perspectives. Renewable and Sustainable Energy Reviews 2013, 19, 164-172. 22. Zafirakis, D.; Chalvatzis, K. J.; Baiocchi, G., Embodied CO2 emissions and cross-border electricity trade in Europe: Rebalancing burden sharing with energy storage. Appl Energ 2015, 143, 283-300. 23. Colett, J. S.; Kelly, J. C.; Keoleian, G. A., Using Nested Average Electricity Allocation Protocols to Characterize Electrical Grids in Life Cycle Assessment. J Ind Ecol 2016, 20, (1), 29-41. 24. Treyer, K.; Bauer, C., Life cycle inventories of electricity generation and power supply in version 3 of the ecoinvent database—part II: electricity markets. The International Journal of Life Cycle Assessment 2016, 21, (9), 1255-1268. 25. Itten, R.; Frischknecht, R.; Stucki, M.; Scherrer, P.; Psi, I., Life cycle inventories of electricity mixes and grid. In 2012. 26. Miller, R. E., Input-output analysis : foundations and extensions / Ronald E. Miller and Peter D. Blair. Cambridge University Press: Cambridge, UK ; New York, 2009. 27. UNComtrade. UN Comtrade Database (http://comtrade.un.org/data); United Nations Statistics Division, 2016. 28. Eggleston, H.; Buendia, L.; Miwa, K.; Ngara, T.; Tanabe, K., IPCC guidelines for national greenhouse gas inventories. Institute for Global Environmental Strategies, Hayama, Japan 2006, 2, 48-56. 29. Ji, L.; Jia, X.; Chiu, A. S.; Xu, M., Global Electricity Trade Network: Structures and Implications. PloS one 2016, 11, (8), e0160869. 30. Hijmans, R.; Kapoor, J.; Wieczorek, J.; Garcia, N.; Maunahan, A.; Rala, A.; Alex, M., GADM database of Global Administrative Areas. In https://gadm.org/data.html, 2016. 31. European Commission, Joint Research Centre (JRC)/Netherlands Environmental Assessment Agency (PBL), Emission Database for Global Atmospheric Research (EDGAR), release version 4.2. In edgar. jrc. ec. europa. eu (last access: 2016), 2015. 32. USEPA Greenhouse Gas Equivalencies Calculator. https://www.epa.gov/energy/greenhouse-gas-equivalencies-calculator 33. Weber, C. L.; Jaramillo, P.; Marriott, J.; Samaras, C., Life Cycle Assessment and Grid Electricity: What Do We Know and What Can We Know? Environ Sci Technol 2010, 44, (6), 18951901.

21 ACS Paragon Plus Environment

Environmental Science & Technology

608 609 610 611 612 613 614 615 616 617 618 619 620 621 622 623 624 625

Page 22 of 22

34. Masanet, E.; Chang, Y.; Gopal, A. R.; Larsen, P.; Morrow III, W. R.; Sathre, R.; Shehabi, A.; Zhai, P., Life-cycle assessment of electric power systems. Annual Review of Environment and Resources 2013, 38, 107-136. 35. Cai, W.; Wang, C.; Jin, Z.; Chen, J., Quantifying baseline emission factors of air pollutants in China’s regional power grids. Environ Sci Technol 2013, 47, (8), 3590-3597. 36. Wood, R.; Stadler, K.; Bulavskaya, T.; Lutter, S.; Giljum, S.; de Koning, A.; Kuenen, J.; Schütz, H.; Acosta-Fernández, J.; Usubiaga, A., Global sustainability accounting—developing EXIOBASE for multi-regional footprint analysis. Sustainability 2014, 7, (1), 138-163. 37. Tukker, A.; de Koning, A.; Wood, R.; Hawkins, T.; Lutter, S.; Acosta, J.; Rueda Cantuche, J. M.; Bouwmeester, M.; Oosterhaven, J.; Drosdowski, T., EXIOPOL–development and illustrative analyses of a detailed global MR EE SUT/IOT. Econ Syst Res 2013, 25, (1), 50-70. 38. Ryan, N. A.; Johnson, J. X.; Keoleian, G. A., Comparative Assessment of Models and Methods To Calculate Grid Electricity Emissions. Environ Sci Technol 2016, 50, (17), 8937-8953. 39. Gingerich, D. B.; Sun, X.; Behrer, A. P.; Azevedo, I. L.; Mauter, M. S., Spatially resolved air-water emissions tradeoffs improve regulatory impact analyses for electricity generation. Proceedings of the National Academy of Sciences 2017, 114, (8), 1862-1867. 40. Siler-Evans, K.; Azevedo, I. L.; Morgan, M. G., Marginal Emissions Factors for the U.S. Electricity System. Environ Sci Technol 2012, 46, (9), 4742-4748.

626

22 ACS Paragon Plus Environment