Growing environmental footprint of plastics driven by coal combustion

Research on the environmental impacts from the global value chain of plastics has typically focused on the disposal phase, considered most harmful to the environment and human health. However, the production of plastics is also responsible for substantial environmental, health and socioeconomic impacts. We show that the carbon and particulate-matter-related health footprint of plastics has doubled since 1995, due mainly to growth in plastics production in coal-based economies. Coal-based emissions have quadrupled since 1995, causing almost half of the plastics-related carbon and particulate-matter-related health footprint in 2015. Plastics-related carbon footprints of China’s transportation, Indonesia’s electronics industry and India’s construction sector have increased more than 50-fold since 1995. In 2015, plastics caused 4.5% of global greenhouse gas emissions. Moreover, 6% of global coal electricity is used for plastics production. The European Union and the United States have increasingly consumed plastics produced in coal-based economies. In 2015, 85% of the workforce required for plastics consumed by the European Union and the United States was employed abroad, but 80% of the related value added was generated domestically. As high-income regions have outsourced the energy-intensive steps of plastics production to coal-based economies, renewable energy investments throughout the plastics value chain are critical for sustainable production and consumption of plastics. Despite concerns about plastics in the environment, not enough attention is paid to the impacts of the various stages of the plastics value chain globally. This study finds that most environmental and socioeconomic impacts from plastics are due to their growing production in coal-based economies.

T he global demand for plastics has quadrupled over the past four decades 1 and is projected to further increase in the future, intensifying the impacts on the environment and human health [2][3][4][5] . Strategies for the sustainable production and consumption of plastics require information on the value chain of plastics [6][7][8][9] , such as fossil resource extraction and processing, resin production, manufacturing into plastics products, plastics use and end-of-life treatments. Many studies have addressed the environmental impacts of (micro)plastics pollution 2,4-6,8,10-14 and plastics incineration 7,15,16 . Less attention has been paid to plastics production, which also has substantial environmental impacts, such as those caused by the release of greenhouse gas (GHG) 17-20 emissions. Furthermore, plastics production induces health impacts, such as through the release of particulate-matter (PM) emissions, and socioeconomic impacts, such as by employing a workforce and creating value added.
Since the value chain of plastics spans the entire globe, plastics are often produced in a different country from the one in which they are ultimately consumed 21 . Therefore, the environmental, health and socioeconomic impacts resulting from one country's plastics consumption can occur elsewhere around the globe. Multiregional input-output (MRIO) analysis allows these impacts to be assessed along the global value chain [22][23][24][25][26][27][28] . However, the accuracy of results from standard MRIO analysis has been limited when analysing the cumulative impacts of materials such as plastics due to double counting [29][30][31][32] . For example, when assessing cumulative GHG emissions (including upstream emissions) of primary plastics production and plastics recycling, double counting occurs because some primary plastics are ultimately recycled. In standard MRIO analysis 32 , the emissions of these primary plastics are counted again as upstream emissions in plastics recycling (a detailed explanation is provided in Cabernard et al. 31 ).
In this article, we apply an enhanced method based on MRIO analysis that prevents double counting [29][30][31]33 to assess the environmental impacts of global plastics production from 1995 to 2030 and extend this method to evaluate the role of coal combustion (Methods). We assess GHG emissions that occur in the global plastics value chain, called the carbon footprint of plastics. In this context, the study highlights the importance of plastics production, including resin production, manufacturing into plastics products and related upstream activities. Moreover, we analyse fossil resources used as a fuel and feedstock for plastics production, called here the fossil resource footprint of plastics. To evaluate the role of trade, the link between plastics-producing and consuming regions is mapped. Further, we analyse the future evolution of the global carbon footprint of plastics assuming that the world follows the International Energy Agency's (IEA's) projection for a 2 °C or 6 °C scenario 34,35 . Finally, this study analyses the PM health impacts, the workforce employed and value added created in the global plastics production chain to provide an overview of the health and socioeconomic impacts (Methods).
(6%, year 2015, Extended Data Fig. 1). Consequently, coal-based emissions caused almost half of the carbon footprint of global plastics production in 2015, due mainly to electricity and heat supply from coal for resin production and manufacturing into plastics products (Fig. 1a,b). In total, 6% of global coal electricity was used for plastics production in 2015.
Due to the increased reliance on coal, the fossil resource footprint of plastics, including fossil resources used as fuel and feedstock for plastics production, has tripled since 1995 (Fig. 2a). Fossil fuels combusted for global plastics production released a total of 1.7 GtCO 2 e in 2015 (Fig. 1b, 88%). The carbon contained in fossil resources used as feedstock for plastics production accounted for another 890 MtCO 2 e (meaning this amount would be released if all plastics produced in 2015 were combusted without credits from energy recovery). Thus, twice as much fossil carbon is combusted as fuel for plastics production (1.7 GtCO 2 e) than contained as feedstock in plastics (890 MtCO 2 e). Our results further indicate that if all plastics produced in 2015 were incinerated, this would increase the annual carbon footprint of plastics by 19% (350 MtCO 2 e, subtracting credits from energy recovery 36 ). While the GHG emissions of plastics incineration are commonly known 7,8,15 , our results show that even in a worst-case scenario in which all plastics were combusted, the major share of GHG emissions would still occur in the production phase.   a-e, The sum of each horizontal bar of the flow chart refers to the carbon footprint of global plastics production in 2015 (1.9 GtCO 2 e, 100%) and allocates it to the different perspectives in the global value chain: the sectors where GHG emissions are released (a); the processes that release GHG emissions, that is, combustion (88%) and non-combustion (7%) processes of fossil fuels (left side: 95% in total) and other processes (right side: 5 % in total): biogenic emissions (2.7%), hydrofluorocarbon emissions (1.7%) and cement production (0.8%) (b); the regions where GHG emissions are released (production perspective) (c); the regions where plastics are finally used (consumption perspective) (d); and the end products or sectors where plastics are finally used: plastics packaging material is allocated to the end product or sector where it is used (for example, food packaging in the food sector) (e). The flows show the linkages between the perspectives. To enhance clarity, linkages contributing to less than 0.1% of the plastics-related carbon footprint are not shown. The small graphs on the right show the temporal evolution of the carbon footprint of global plastics production for each perspective (a-e) over the past two decades  and in the future (2020-2030) if the world follows the IEa's projection for a 2 °C or 6 °C scenario 34,35 (but not accounting for the decrease in global GHG emissions in 2020 due to the COVID-19 pandemic, as in Liu et al. 66 ). The colours of the graphs on the right correspond to the bars of the flow chart. The red line in each graph represents the global share. The end-of-life stages of plastics (recycling, incineration and landfills) contributed another 120 MtCO 2 e in 2015, not illustrated in this figure. a detailed value-chain analysis of the carbon footprint of global plastics recycling is shown in Extended Data Fig. 2.
Growth in plastics production in coal-based emerging economies, such as China, Indonesia and South Africa, was the major driver of the increasing carbon footprint of plastics ( Fig. 1). Since 1995, China's plastics-related carbon footprint has more than tripled from both a production and a consumption perspective ( Fig. 1b-d). In 2015, 40% of the global plastics-related carbon footprint and more than 60% of the related coal-based emissions were caused in China. In Indonesia, coal mining for plastics production has increased by a factor of 300 since 1995 (Fig. 2b). In 2015, 15% of Indonesia's totally mined coal was used for plastics production, either domestically or abroad. More than 10% of Indonesia's total domestic GHG emissions were attributed to plastics production in 2015. South Africa's plastics-related carbon footprint has increased tenfold since 1995, and 95% of it was caused by domestic coal consumption in 2015. The reason for this is that South Africa uses coal not only to supply electricity and heat but also as a feedstock for plastics production 37 . Almost half of the plastics produced in China and South Africa were exported, such as to the European Union and to the United States (Fig. 1c,d).
High-income regions, such as the European Union and the United States, contributed substantially to the increasing global carbon footprint of plastics by their rising demand for plastics produced in lower-income regions, particularly in coal-based economies such as China (Figs. 1 and 3 and Extended Data Fig. 3). Due to the outsourcing of plastics production to lower-income regions, the European Union's plastics-related carbon footprint increased although their domestic plastics-related GHG emissions decreased (Extended Data Fig. 4). In 2015, two-thirds of the European Union's plastics-related carbon footprint was emitted abroad, mainly in lower-income regions with less-stringent environmental policies. The fraction of the plastics-related carbon footprint caused abroad was even higher for Australia, Canada and the United States (>80% in 2015, Fig. 3b). Since 1995, the United States' plastics-related carbon footprint generated abroad quadrupled, while their plastics-related domestic GHG emissions decreased (Extended Data Fig. 5). Consequently, the fraction of the United States' plastics-related carbon footprint induced abroad increased from 39% in 1995 to 78% in 2015. One-third of the United States' plastics-related carbon footprint in 2015 occurred in China (Fig. 1c, Supplementary Fig. 2 for results on the carbon intensity of plastics resin production per region).
In addition to increased plastics exports, emerging economies contributed to the rising global carbon footprint of plastics by their growing plastics demand, due mainly to their growing infrastructures, transportation systems and digitalization. Since 1995, the plastics-related carbon footprint of China's transportation system, Indonesia's electronics industry and India's construction sector has increased more than 50-fold. In 2015, 15% of the global carbon footprint of plastics was attributed to plastics used for construction, and almost half of these emissions were attributed to China's construction (Fig. 1d,e). In addition, plastics are responsible for 15% of the  carbon footprint of the global automotive industry, and more than one-third of these GHG emissions were attributed to China's automotive industry. While China also manufactured most of the plastics embodied in electronics (65%), a smaller portion (42%) was used by China itself, and the majority (58%) was exported to other regions.
Climate change scenarios. If the world follows the IEA's projection for a 6 °C scenario 34,35 , the global carbon footprint of plastics would grow by 31% from 2015 to 2030, and thus almost at the same speed as projected plastics production (+40%, Fig. 1, see Extended Data Fig. 6 for per capita projections). The strongest increase is expected for coal-based economies, such as China, India and South Africa, while the largest international customers of plastics produced in these coal-based economies continue to be the European Union and the United States. Following the IEA's projection for a 2 °C scenario 34,35 would reduce the global carbon footprint of plastics by 10% from 2015 to 2030, while plastics production would increase by 40%. This would be attributed to investment in renewable energy production, mainly clean electricity, and improved energy efficiency in resin production. However, coal combustion would still contribute more than a third of the carbon footprint of plastics in 2030 if the world followed the IEA's projection for a 2 °C scenario 34,35 . This highlights the potential of a rapid phase-out of coal to further reduce the carbon footprint of plastics in the future, as keeping global warming below 1.5 °C is critical for preventing major climate-related hazards 38 . Despite the growing demand for plastics in emerging economies, the carbon footprints of plastics remain distinctly higher in high-income regions on a per capita level (Fig. 3a), and this imbalance is projected to persist into the future (Extended Data Fig. 6). Taking into account the projected reduction in plastics-related carbon footprints per region if the IEA's measures for a 2 °C scenario 34,35 were implemented, high-income regions in particular would have a high saving potential at the per capita level, which would also reduce income-related differences in the carbon footprints of plastics.
Health and socioeconomic impacts of plastics production. Since 1995, the global PM health footprint of plastics has increased by 70%, causing the loss of 2.2 million disability-adjusted life years (DALY) and representing 2.8% of global PM health impacts in 2015 (based on EXIOBASE3 22 and Cabernard et al. 31 , Extended Data Figs. 7 and 8). Similar to GHG emissions, plastics production accounted for the majority (96%) of the plastics-related PM health footprint, half of which was attributed to coal combustion (year 2015, Fig. 4a,b). Moreover, the majority (75%) of PM health impacts were caused in China, India, Indonesia and other Asia, while high-income regions increasingly consumed plastics produced in these regions ( Fig. 4c,d). In 2015, the majority of the European Union's (80%) and the United States' (91%) plastics-related PM health footprint was caused abroad, mainly in Asia. Thus, only 8% of the global plastics-related PM health footprint effectively occurred in high-income regions (Extended Data Fig. 9c), although one-third of global plastics was consumed in high-income regions.
The number of employed workers for global plastics production has increased by 54% since 1995. At the same time, the value added created by global plastics production, including workforce compensation (48%), operating surplus (40%) and taxes (12%), has more than doubled. In 2015, 106 million workers (in full-time equivalents (FTE); 2.3% of the global workforce) were employed in plastics production, while the associated value added reached €1.56 trillion (2.6% of global gross domestic product (GDP)). The number of workers and the created value added are unevenly distributed around the globe: although more than 90% of the plastics-related workforce was employed in low-and middle-income regions, more than half of the plastics-related value added was generated in high-income regions (Extended Data Fig. 9d,e). China and other Asia made up the majority of the plastics-related workforce (Fig. 5a), while the European Union generated most of the plastics-related value added (Fig. 5b). In 2015, 70% of the workforce required for the European Union's plastics consumption was employed abroad, but 80% of the related value added was generated domestically. The reliance on the low-paid foreign workforce was even higher for the United States' plastics consumption, where 90% of the workforce was occupied abroad. This means that high-income regions, such as the European Union and the United States, outsource the low-paid steps in the plastics production value chain to lower-income regions and focus on the valuable steps of manufacturing plastics into finished products.

Discussion
Comparison with literature. Other than previous estimates [39][40][41] , which indicated equal amounts of oil used as fossil fuel and feedstock for plastics production, this study showed that twice as much fossil carbon is combusted as fuel for plastics production than contained as feedstock in plastics. The reason for this difference is that our approach accounts not only for the amount of oil, but for GHG emissions released by combusting all types of fossil fuels along the plastics production chain as well, taking the increased reliance on coal into consideration. In this context, coal also contributes to the carbon footprint of plastics due to its increasing use as a feedstock for plastics production, such as in China 42 and South Africa 37 .
The carbon footprint of plastics obtained in this study is higher (+16%) compared with Zheng and Suh 17 , who used bottom-up life-cycle analysis to assess the global carbon footprint of plastics in 2015 but did not account for double counting (see Extended Data Fig. 10 for a comparison of the results). In our study, double counting would have overestimated the carbon footprint of plastics by 29% but was avoided by means of the method applied 31 . The reason the carbon footprint of plastics obtained here is still higher compared with Zheng and Suh 17 is that we considered the regionalized fuel-specific energy mix, such as the increased reliance on coal, while Zheng and Suh 17 calculated with the average global energy mix.
The importance of considering the energy mix is also reflected in the results of the scenarios: in contrast to previous estimates 8, 15 , the plastics' share in total global GHG emissions is projected to  processes of fossil fuels (left side: 88% in total) and other non-combustion processes (right side: 12% in total): metals (7.6%), chemicals (1.8%), cement production (1.5%) and biogenic emissions (1.2%). c, regions where PM health impacts are caused (production perspective). d, regions where plastics are consumed (consumption perspective). e, End products or sectors where plastics are used. Plastics packaging material is allocated to the end product or sector where it is used (for example, food packaging in the food sector).
decrease in the scenarios used here (Fig. 1). One reason for this is the assumption of increased future investment in renewable energy technologies and improved energy efficiency in the plastics production value chain, which can together exercise substantial leverage for reducing the global carbon footprint of plastics. Another reason could be that the growth in plastics production may be underestimated in the scenarios of this study. This highlights the importance of improved scenarios that take both changes in the energy mix and plastics production into consideration.
Policy implications. In addition to commonly known issues about (micro) plastics pollution 2-6,8,10-14 and plastics incineration 7,15,16 , this study highlights the need for improved policy measures to reduce the increasing carbon footprint of plastics production, which bears the major share of the plastics-related GHG emissions (even in a worst-case scenario where all plastics would be incinerated). On the one hand, our results underscore the importance of ongoing initiatives to reduce primary plastics production by avoiding, reusing and recycling plastics as discussed in the context of circular economy 8,9,39,43 . However, a general ban on plastics is counterproductive as alternative materials often have higher environmental impacts 44 . On the other hand, this study highlights the particularly strong leverage in the plastics production chain itself to reduce the carbon footprint of plastics. Efficient measures include phasing out coal, transitioning to renewables and improving the energy efficiency in the plastics production process. As shown here for the past and future, decreasing the emissions in high-income regions as specified in the Paris Agreement is not sufficient. Such an approach even fosters a shift of plastics production to emerging regions with less-stringent environmental policies and limited economic power to implement state-of-the-art low-carbon technology. Thus, it is important that high-income regions invest in clean energy production throughout the supply chain.
Since renewable energy investments are currently hampered by the lack of economic incentives 45-47 , a key measure required is for the government to implement a carbon price, such as carbon taxes, cap-and-trade emissions schemes and renewable energy subsidies. Once a reasonable carbon price is in place on a producer and consumer level, the economic opportunities to decarbonize (plastics-related) supply-chain emissions will be particularly attractive for consumer-facing companies in high-income regions 45 . One reason for this is that consumer-facing companies can often reduce many more GHG emissions by decarbonizing their (plastics-related) supply chain than by focusing exclusively on their direct emissions. Another reason is that decarbonization of the (plastics-related) supply chain is less costly for consumer-facing companies in high-income regions than for plastics-producing companies in emerging markets 45 . This is because plastics represent only a small portion of end-use prices, and consumer-facing companies usually generate higher value per emissions produced (as shown in this study for high-income regions; Figs. 1c,d and 5).
In addition to carbon pricing, creating transparency with regard to (plastic-related) supply-chain emissions is a key measure required to incentivize consumer-facing companies to reduce their supply-chain emissions 45 since it would allow companies to benefit from increased consumer demand for green products. In this study, transparency with regard to plastics-related supply-chain emissions has been improved on a global and country level using the methodology of Cabernard et al. 31 . Thus, future research is crucial to provide company-level data and to thereby close research information gaps.
One incentive for emerging economies to implement a carbon price to foster renewable energy investments is that such a measure can also reduce PM health impacts and thus local health costs 48 (because unlike the global climate impact of GHG emissions, health impacts depend on the location where they occur). This might be of particular interest to China, which bears the burden of most  plastics production (including resin production, manufacturing into plastics products and related upstream activities) from production and consumption perspectives in 2015. The production perspective indicates the region where the plastics-related workforce is employed and the value added is created. The consumption perspective allocates the plastics-related workforce and value added to the region where plastics were finally consumed. Coloured flows refer to plastics consumed in another region than produced due to international trade. The value added of global plastics production includes workforce compensation (48%), operating surplus (40%) and taxes (12%).
plastic-related PM health impacts but has limited emission reduction potential from further flue gas treatment improvements 49 . By contrast, the installation of advanced flue gas treatment could substantially reduce plastic-related PM health impacts in India, Indonesia and other Asian countries 49 . Another incentive for emerging economies to implement a carbon price is the technical and economic opportunities that decarbonization presents in the long term, such as the future competitive advantage to be gained in view of the energy transition and the growth of new industries and services 46,47,50 . Consequently, decarbonization efforts might also mitigate the income-related regional imbalance between plastics-related PM health impacts, workforce and value added shown in this study.

Methods
Background. Our methodology is based on MRIO analysis. In MRIO analysis, the world economy is aggregated into a given number of regions and industrial sectors, whose transactional, environmental and socioeconomic accounts are captured for a given period 22,51,52 . This allows the environmental and socioeconomic impacts of these regions and sectors to be assessed along the global value chain [23][24][25][26]53 . Standard MRIO analysis is accurate for assessing the impacts of plastics end products such as tableware directly purchased by households [25][26][27][28]54,55 . However, most plastics are not directly used by households but intermediately by the industry, such as plastics for food packaging, for construction materials or in electronic products. The impacts of these intermediately used plastics could not previously be properly addressed by standard MRIO methodology 25-28,54,55 due to double counting [29][30][31][32]56 . This methodical issue of double counting existed not only for plastics but for any materials, such as minerals, fossil resources and biomass, when assessing their cumulative impacts, including cumulative upstream and direct impacts. For example, when assessing cumulative GHG emissions of steel and coke production, double counting occurs because some coke is used for steel production. This means that the emissions related to coke production are counted again as upstream emissions in steel production. Similarly, when assessing cumulative GHG emissions of primary plastics production and plastics recycling, double counting occurs because some primary plastics are ultimately recycled. In standard MRIO analysis 32 , the emissions of these primary plastics are counted again as upstream emissions in plastics recycling (a detailed explanation is provided here in ref. 31 ).
Recently, a method was developed that prevents double counting 29,30 . It was extended to MRIO analysis 31 to assess the cumulative impacts of any sector or regions and to track them upstream and downstream the global value chain (without double counting). It was applied to assess GHG emissions 31,33 , PM-related health impacts 31 and other environmental and socioeconomic impacts of global material production 31 . An explanation of how double counting is prevented and how this affects the environmental footprint of global material production is provided by Cabernard et al. 31 .
Overview and scope. We used this method 31 to assess fossil resource extraction, GHG emissions (carbon footprint), PM health impacts and socioeconomic impacts (value added and workforce) along the global plastics value chain. We extended this method to assess the processes of GHG and PM emissions release and the type of fossil fuels used at the different life-cycle stages in the plastics value chain. Further, we compared the amounts of fossil resources used as fuels and feedstock for plastics production. To project how the global carbon footprint of plastics will evolve until 2030, we applied our method on a forward-looking MRIO database 35 , covering the 2 °C and 6 °C scenarios of the IEA 34 . This study includes cumulative upstream and direct impacts of resin production, further manufacturing, recycling, incineration and landfills of all types of plastics and rubber, called footprint of plastics in what follows. Impacts related to plastics waste in the environment, such as microplastics pollution, were not addressed.
Database and terminology. We used the global MRIO database EXIOBASE3 (version 3.7) 22 , which aggregates the global economy into 163 industrial sectors for 44 countries and 5 rest-of-the-world regions (49 regions). This results in 7,987 sector-region combinations, whose economic flows are recorded in the coefficient matrix A and the final demand matrix Y for each year from 1995 to 2015 in constant prices 22 . The coefficient matrix A indicates the monetary input of each sector-region combination per monetary output of each sector-region combination (7,987 rows × 7,987 columns, in Euro Euro -1 ). The final demand matrix Y indicates the final demand of each region for each sector-region combination (7,987 rows × 49 columns, in Euro). For each sector-region combination and year, the EXIOBASE3 database provides a set of environmental and socioeconomic accounts (i), which was extended to GHG emissions and PM health impacts 31 , and which is called the impact coefficient matrix (D, i rows × 7,987 columns, impact/Euro). In this context, each type of GHG emission listed by the impact coefficient matrix of EXIOBASE3 (CO 2 , CH 4 , N 2 O, hydrofluorocarbons and perfluorinated compounds) was multiplied with the respective global warming potential to estimate the amount of emitted CO 2 e. To calculate PM health impacts, each type of PM emission listed by the impact coefficient matrix of EXIOBASE3 (PM 2.5 , NO x , SO x and NH 3 ) was multiplied with sector-specific characterization factors (CFs) to measure the human burden of disease in DALY, as done by Cabernard et al. 31 . These CFs were derived from the United Nations Environment Programme (UNEP) and Society of Environmental Toxicology and Chemistry (SETAC) 57 , who provided CFs for PM 2.5 , NO x , SO x and NH 3 emitted to different compartments (outdoor urban, outdoor rural, indoor urban and indoor rural) and CFs for PM 2.5 emitted at different heights (ground level, low stack, high stack, very high stack). Since we cannot distinguish emission locations by population density in EXIOBASE3, we assumed rural outdoor emissions (with ground-level emissions for PM 2.5 ) for all extracting sectors, we averaged rural and urban outdoor emissions for the manufacturing sectors and we assumed the average of urban outdoor emissions for the rest of the economy, as done by Cabernard et al. 31 . The indicators 'fossil resource extraction' (in tons of extracted fossil resources), 'value added' (in Euros) and 'workforce' (in number of people working full-time) were directly adopted from the impact coefficient matrix of EXIOBASE3 22 .
Environmental and socioeconomic footprints of plastics. The method from Cabernard et al. 31 can be applied to any industrial sector(s) and region(s) of any MRIO system. The first step is to define the sectors and regions of interest, called target sectors and target regions. Out of the 163 sectors covered by EXIOBASE3, we defined five sectors as target sectors: 'plastics resin production' , 'manufacturing of plastics and rubber products' , 'plastics recycling' , 'incineration of plastics' and 'landfills of plastics' . Since we analysed the global plastics value chain, we defined all 49 regions as target regions. This resulted in 245 target-sector regions (T, 5 target sectors × 49 target regions) referring to the global plastics economy, and 7,742 non-target-sector regions (O) referring to the global non-plastics economy. Both target-sector regions and non-target-sector regions represent the global economy (all, 7,987 sector regions). Following the method from Cabernard et al. 31 , we calculated the footprints (FP) of the entire plastics value chain for the indicators (i) fossil resource extraction, GHG emissions, PM health impacts, value added and workforce: Equation (1) describes the cumulative upstream and direct impacts of resin production, further manufacturing, recycling, incineration and landfills of all types of plastics and rubber (called footprint of plastics in this study) but excludes the impacts related to the use phase of plastics as well as the impacts of plastics disposed to the environment. In equation (1) ) and the direct final demand for the outputs of the non-plastics economy (Y o-all ). This includes, for example, plastics used for electronics consumed by households. Diagonalizing the impact coefficients (diag(D i-all )) allocates the global plastics footprints to the sector and region where the fossil resources are extracted, the emissions are released, the impacts are caused, the value added is generated and the workforce is employed. A detailed explanation of equation (1) and a list of all terms is provided in Cabernard et al. 31 by way of the example of global material production.

Division by plastics sectors.
To divide the environmental footprints of plastics between the plastics sectors plastics resin production, manufacturing of plastics and rubber products, plastics recycling, incineration of plastics and landfills of plastics' , we divided equation (1) into two terms referring to direct and cumulative indirect impacts, respectively: The first term refers to the direct impacts of the plastics sectors, where D i-T equals the impact coefficients of the plastics sectors (T) and L T-T is an excerpt of the Leontief inverse, which indicates the cumulated inputs of the plastics economy into the plastics economy. The second term indicates the cumulative indirect impacts of the plastics sectors caused in the upstream chain by the non-plastics economy:  (2) is the same as equation (1), but diagonalizing the cumulative upstream impact coefficients (diag (D i−O × L O−T )) allocates the cumulative indirect impacts (caused in the upstream chain by the non-plastics economy) to the plastics sector (while equation (1) allocates the cumulative indirect impacts to the sector in the upstream chain). To track the footprints along the global value chain, we diagonalized the different terms of equations (1) and (2), as explained by Cabernard et al. 31 .
Distinguishing different processes of GHG and PM emissions. The impact coefficient matrix of EXIOBASE3 indicates the GHG and PM emissions per Euro of output for each sector-region combination and for different processes. These processes include combustion of fossil fuels, non-combustion processes of fossil resources, minerals extraction and processing, and biogenic GHG emissions due to agriculture for bio-based plastics. To allocate the carbon and PM health footprint of global plastics production to each of these processes, we applied the process-specific impact coefficients (D iprocess ) to equations (1) and (2).
Types of fossil fuels causing combustion-related emissions. For the carbon and PM health footprint related to fossil fuel combustion, we further investigated which types of fossil fuels, such as coal, petroleum, natural gas and coke, are combusted in the global plastics value chain. Since this is not indicated by EXIOBASE3, we weighted the combustion-related carbon and PM health footprint of global plastics production with fuel-specific contribution matrices (C GHG fuel−all and C PM fuel−all , 4 rows × 7,987 columns): C GHG fuel−all and C PM fuel−all indicate the contribution (in %) to GHG emissions and PM health impacts, respectively, by the combustion of coal, petroleum, natural gas or coke (4 rows) in each sector-region combination of the EXIOBASE3 database (7,987 columns). Since each column sums up to 100%, the resulting footprints are the same as in equations (1) and (2). However, the application of the fuel-specific contribution matrices allows the types of fossil fuels combusted at the different life-cycle stages of the global plastics value chain to be distinguished, as illustrated in Figs. 1b and 4b. We compiled fuel-specific GHG emissions factors (e GHG/phy fuel , kgCO 2 e t -1 ) for coal, petroleum, natural gas and coke 58 . For coal, we weighted the emission factor of anthracite, bituminous coal, sub-bituminous coal, lignite and peat with the extracted volumes indicated by the United Nations International Resource Panel Material Flows Database 59 . Fuel and region-specific PM health impact factors per amount of fuel (e PM/phy fuel , DALY t -1 ) were calculated on the basis of a global coal emission inventory 49 (with a few outliers and the small units below 50 MW with lower data quality being removed) in combination with a highly resolved regionalized PM health impact assessment methodology 60 . Emission factors for petroleum and natural gas were obtained from the European Environment Agency 61 . Health impact factors were approximated on the basis of Oberschelp et al. 49,60 and an average electrical efficiency of 40%, assuming the same spatio-temporal emission distribution patterns as for coal power generation. The health impact factors indicate the DALY per amount of combusted coal, petroleum, natural gas and coke for each of the 49 regions covered by EXIOBASE3 for the year 2015. For coke, we assumed the same region-specific PM health impact factors as for coal.
Fossil resources used as fuels and feedstock. Fossil resources are used both as fuels and feedstock for plastics production. While fossil resources used as fuels release GHG emissions during plastics production and related upstream activities, fossil resources used as feedstock release GHG emissions only if plastics are incinerated or composted. We calculated the GHG emissions of fossil resources used as fuels for plastics production (FP resin+manu i,comb ) by applying equations (1) and (2) to the combustion-related GHG emissions indicated by EXIOBASE3 separately. To derive the GHG emission potential of fossil resources used as feedstock for plastics resin production (E resin GHG,feedstock , in kgCO 2 e), we subtracted the GHG emissions directly released through plastics resin production (E resin GHG, direct , in kgCO 2 e) from the GHG emission potential of the fossil resource input into plastics resin production (E resin GHG,input , in kgCO 2 e) for the year 2015: We derived the GHG emission potential of the fossil resource input into plastics resin production (E resin GHG,input ) from the element-wise product of the monetary transaction matrix of EXIOBASE3 (T mon fuel−resin , in Euros), the inverted price vector of fossil resources ( p phy/mon fuel , in t Euro -1 ) and the fuel-specific GHG emission factor vector (e GHG/phy fuel , in kgCO 2 e t -1 ) for the year 2015: Thereby, FP plastics GHG,prod refers to the carbon footprint of global plastics production (based on equation (1)), and s refers to GHG emissions savings through heat and electricity recovery during plastics waste incineration. According to the ecoinvent 3.4 36 life-cycle database, the incineration of 1 kg of plastics generates on average 3.92 MJ of electricity and 7.66 MJ of useful heat, which we multiplied with the life-cycle GHG emissions of the global electricity (0.544 kgCO 2 e MJ -1 ) 42,62 and heat mix (0.073 kgCO 2 e MJ -1 ) 63 as well as global plastics production in 2015 to assess carbon savings through energy recovery, as done by Zheng and Suh 17 . This estimation assumes a scenario in which all plastics produced in 2015 end up as waste. In reality, some plastics end up in stocks (for example, infrastructure), while some plastics produced in previous years add to the waste. Further, heat and electricity savings (s) are based on average values 36 .

Climate change scenarios.
To evaluate how the carbon footprint of plastics evolves until 2030 under different climate change scenarios, we applied the procedure described in equations (1)-(7) on a forward-looking version of EXIOBASE3 35 , which is based on climate change scenarios of the IEA 34 . The forward-looking MRIO database involves different measures required to pursue the 2 °C and 6 °C scenarios until 2100 as projected by the Energy Technology Perspectives of the IEA 34 . The most important measures include increased investment in renewable energy technologies and improvements in material and energy efficiency on the basis of the growing world population and GDP. An overview of the measures implemented is provided in Table 1 in Wiebe et al. 35 . Note that the United States' current boom in shale gas extraction may increase their domestic plastics production and related GHG emissions in the future, 64,65 which was not accounted for in the scenarios used here 34,35 .
Reporting Summary. Further information on research design is available in the Nature Research Reporting Summary linked to this article.

Data availability
The study is based on the multiregional input-output (MRIO) database EXIOBASE3 22 and the impact assessment methodologies recommended by UNEP-SETAC 57 . A forward-looking version of EXIOBASE3 from Wiebe et al. 35  This study's global carbon footprint of plastics production is 16% higher compared to Zheng and Suh 17 . The reason for this might be that the MrIO-based approach applied here accounts for the fuel-specific energy mix, such as the increased reliance on coal, while Zheng and Suh 17  Last updated by author(s): Oct 8, 2021 Reporting Summary Nature Research wishes to improve the reproducibility of the work that we publish. This form provides structure for consistency and transparency in reporting. For further information on Nature Research policies, see our Editorial Policies and the Editorial Policy Checklist.

Statistics
For all statistical analyses, confirm that the following items are present in the figure legend, table legend, main text, or Methods section.

n/a Confirmed
The exact sample size (n) for each experimental group/condition, given as a discrete number and unit of measurement A statement on whether measurements were taken from distinct samples or whether the same sample was measured repeatedly The statistical test(s) used AND whether they are one-or two-sided Only common tests should be described solely by name; describe more complex techniques in the Methods section.
A description of all covariates tested A description of any assumptions or corrections, such as tests of normality and adjustment for multiple comparisons A full description of the statistical parameters including central tendency (e.g. means) or other basic estimates (e.g. regression coefficient) AND variation (e.g. standard deviation) or associated estimates of uncertainty (e.g. confidence intervals) For null hypothesis testing, the test statistic (e.g. F, t, r) with confidence intervals, effect sizes, degrees of freedom and P value noted

Software and code
Policy information about availability of computer code Data collection Data were compiled from EXIOBASE3 (22) based on the impact assessment methodology recommended by UNEP-SETAC (57), and by applying and extending the method from Cabernard et al (31) to assess the type of fossil resources used as fuel and feedstock for plastics production (based on Oberschelp et al (49,60) for the PM health impact assessment).

Data analysis
Data were analyzed with matlab and illustrated with R-studio and tableau.
For manuscripts utilizing custom algorithms or software that are central to the research but not yet described in published literature, software must be made available to editors and reviewers. We strongly encourage code deposition in a community repository (e.g. GitHub). See the Nature Research guidelines for submitting code & software for further information.

Data
Policy information about availability of data All manuscripts must include a data availability statement. This statement should provide the following information, where applicable: -Accession codes, unique identifiers, or web links for publicly available datasets -A list of figures that have associated raw data -A description of any restrictions on data availability The study is based on the multiregional input-output (MRIO) database EXIOBASE322 and the impact assessment methodologies recommended by UNEP-SETAC57. A forward-looking version of EXIOBASE3 from Wiebe et al (35) was used for future scenarios, based on climate change scenarios from the International Energy Agency (34). Further results of this study are presented in the Extended Data Figures 1-10 and the SI (SupportingInformation_Plastics.pdf). All data are available upon reasonable request. Additional results on the plastics value chain and its environmental impacts, e.g., for specific countries, can be calculated using the tool from Cabernard et al (31).

nature research | reporting summary
April 2020 Field-specific reporting Please select the one below that is the best fit for your research. If you are not sure, read the appropriate sections before making your selection. All studies must disclose on these points even when the disclosure is negative.

Study description
The objective of this work was to assess the environmental and health impacts of global plastics production from 1995 to 2030, taking different climate change scenarios into consideration. In contrast to previous research, our approach accounts for the regionalized energy mix, such as the increased reliance on coal to produce plastics, which is shown to drive the growing carbon and health footprint of plastics. Due to the shift of plastics production to coal-based economies, mainly China, the carbon and health footprint of plastics has doubled since 1995. In 2015, plastics caused 4.5% of global greenhouse (GHG) emissions and 6% of global coal electricity was used for plastics production, which is higher than generally assumed owing to our methodology. Through enhanced multi-regional input-output analysis, we map the global plastics value chain, such as the link of plastics production and consumption. It shows that the rising impacts of plastics are driven by emerging economies' growing infrastructures, while highincome regions, such as the EU and USA, increasingly consume plastics produced in coal-based economies.
Research sample n/a Sampling strategy n/a Data collection Data were compiled from EXIOBASE3 (22) based on the impact assessment methodology recommended by UNEP-SETAC (57), and by applying and extending the method from Cabernard et al (31) to assess the type of fossil resources used as fuel and feedstock for plastics production (based on Oberschelp et al (49,60) for the PM health impact assessment).
Timing and spatial scale The study covers the period from 1995 to 2030, referring to the global production, consumption, and trade of plastics.

Data exclusions
No data were excluded.

Reproducibility
The matlab code to calculate the results of this study is available upon reasonable request and can be used to reproduce the results of this study.