key: cord-0802988-rqanyhb0 authors: Ma, Ning; Yin, GuangWei; Li, Huajiao; Sun, WenLi; Wang, Ze; Liu, Gang; Xie, Dan title: The optimal industrial carbon tax for China under carbon intensity constraints: a dynamic input–output optimization model date: 2022-03-12 journal: Environ Sci Pollut Res Int DOI: 10.1007/s11356-022-19162-6 sha: 4696319dc1b2e8afcf3a24b4fbdf7099ab01440f doc_id: 802988 cord_uid: rqanyhb0 To reduce carbon emissions, the Chinese government is considering introducing a differentiated industrial carbon tax on enterprises outside the carbon trading market in the future. An efficient carbon tax must consider not only how carbon taxes impact the current economy but also how the size of the tax should be adjusted across time due to external changes. To calculate the optimal industrial carbon tax for China which is subject to certain constraints, this paper investigates the economic and environmental effects of four possible industrial carbon tax rate models under carbon intensity constraints from 2021 to 2030 by a dynamic input–output optimization model. The results show that the dynamic tax rate model leads to larger fluctuations in GDP growth than the other tax models, with a low initial tax rate in the beginning and a high tax rate exceeding ¥180/t in 2030. Second, a large quantity of capital stock is distributed across the energy-intensive industries, which leads the existing capital investment structure to be path-dependent. This offsets the performance of carbon taxes. Third, indirect energy-intensive industries such as construction and transport are insensitive to the industrial carbon tax. Finally, comparing the impacts of the four tax rate models, the optimal industrial carbon tax for China is found to be a fixed differentiated tax rate, in which energy-intensive sectors are taxed ¥75/t and low-carbon sectors are taxed ¥50/t. SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s11356-022-19162-6. As a responsible country that actively reduces CO 2 emissions, the Chinese government promised to reach the peak of CO 2 emissions by 2030 and achieve carbon neutrality by 2060 . The goal of carbon neutrality is not an issue of energy substitution and technological breakthroughs but puts forward new and hard constraints on China's sustainable development (Chen and Lin 2021; Zhang 2021) . To promote the enterprises to manage CO 2 emissions and transform to a green production model, China has launched a national carbon emission trading system (namely, ETS) covering all large enterprises in key energy-intensive industries in 2017. However, the high fluctuation of the carbon price in the ETS market causes the expected return on investments in low-carbon technology upgrading to be uncertain, which is unfavorable to carbon neutrality. Given the extremely large amount of carbon emissions accompanied by vast diversity across different regions and industries, the government is considering introducing a differentiated industrial carbon tax on enterprises outside the carbon trading market in the future (National Development and Reform Commission Responsible Editor: Roula Inglesi-Lotz * Huajiao Li hli@cugb.edu.cn 1 2018). Because more than 70% of energy consumption and carbon emissions come from industries in the production field, choosing an appropriate industrial carbon tax rate for China is important for policymakers to reduce carbon emissions from industries (Weng et al. 2021) . Furthermore, as China is the largest CO 2 emitter in the world, the analysis of the optimal industrial carbon tax for China can provide valuable information for emission abatement across other regions in the world. Designing an efficient tax rate is a vital task. Regarding the carbon tax rate model, the existing literature mainly proposes two dimensions: uniform (Li et al. 2019; Jia and Lin 2020) vs. differentiated carbon taxes (Tan et al. 2021; Weng et al. 2021 ) and dynamic vs. fixed carbon taxes (Cao et al. 2019; Chien et al. 2021 ). The former reflects whether a uniform tax rate is applied to different taxpayers, and the latter reflects whether the tax rate changes over time. Thus, there are four basic carbon tax rate models: the uniform tax rate, the differentiated tax rate, the fixed tax rate, and the dynamic tax rate. Most of the existing literature explores the environmental policy performance by a fixed uniform carbon tax (Alegoz et al. 2021; Carroll and Stevens 2021) . However, the fixed uniform carbon tax neglects the obvious difference between taxpayers and fails to reflect the pollutant-pays principle. That will cause inequalities and offset the policy efficiency of carbon taxes . Meanwhile, an efficient carbon tax must consider not only how carbon taxes impact the current economy but also how the size of the tax should be adjusted across time due to external changes from technological improvements, environmental constraints, etc. (Engström and Gars 2015) . Therefore, it is necessary to explore the effects of different carbon tax rate models deeply. To comprehensively reflect the mix of the carbon tax rate model, the sub-models include the fixed uniform tax rate, the fixed differentiated tax rate, the dynamic uniform tax rate, and the dynamic differentiated tax rate. Our analysis is based on the sub-model. The CGE model is the most widely applied method in the existing literature (Tang et al. 2015) . However, it can only simulate the impacts of a given tax on the economy and therefore fails to calculate a reasonable tax rate under certain constraints (Jiang and Zhang 2004) . Against this background, by adopting a dynamic optimal tax module, we construct a dynamic input-output (namely, I-O) optimization model based on the general equilibrium theory. The model is mainly constructed based on balanced relationships in the I-O table. Meanwhile, to fully describe the interactive relationships among all agents, some modules in the dynamic CGE model are concluded as follows: (1) the module of income and expenditure of agents and (2) the module of savings and investment and (3) the dynamic block. The dynamic I-O optimization model is widely applied to water management , carbon tax (Li et al. 2019) , pollution management , etc. Therefore, it is suitable for environmental policy analysis. As to the carbon intensity constraints, China's carbon peak policy strengthens the controls on total carbon emission and carbon intensity, while China's carbon neutrality strategy addresses the overall transformation of both production mode and lifestyle (Zhang 2021) . Considering the existing industrial structure and technical level, emission reduction is an arduous task, which cannot be accomplished in one stroke (Commentator of Economic Daily 2022) . Therefore, in the short term, given the coal-based energy structure, the country is guided to make major technological improvements and energy transition through market tools, such as ETS or carbon tax (Chen and Lin 2021; Song et al. 2021 ). In the long run, the carbon neutrality strategy should vigorously promote a comprehensive and profound transformation of production mode and lifestyle under ecological civilization (Zhang 2021) . That means that China's sustainable development should change the existing energy-intensive development model to a green one and balance the relationship between development and emission reduction. This paper focuses on the carbon intensity constraints in the short term. This paper aims to investigate the optimal industrial carbon tax rate for China from 2021 to 2030 under given carbon intensity constraints via a dynamic I-O optimization model. By comparing the effects of the two basic options for setting tax rates, this paper comprehensively explores strategies for designing an appropriate industrial carbon tax rate. Our contributions are as follows. First, by adopting a dynamic optimal tax module, we construct a dynamic I-O optimization model based on the theory of equilibrium. This model ensures the carbon tax policy analysis considers not only how carbon taxes impact the current economy but also how the size of the tax should be adjusted across time due to external changes. Second, the two dimensions of the tax rate model are comprehensively investigated in one analysis, including whether a uniform tax rate should be applied for different taxpayers and whether the tax rate should change over time. The expansion is helpful to fully understand the effects of different tax rate models. Third, this study extends the policy analysis of the carbon tax to the industry level and explores the impacts of a carbon tax shock on different industries. It provides the necessary information for setting the tax rate within a differentiated industrial carbon tax in the future. The remainder of this paper is structured as follows. Section 2 presents the literature review. Section 3 describes the methodologies, data sources, and scenarios. Section 4 shows the results and discussion. Section 5 concludes and provides policy implications. According to the basic approach to calculating the optimal carbon tax rate used, the existing literature on the optimal carbon tax can be divided into two categories: the Pigouvian tax approach and the standards and pricing approach. The former approach seeks to calculate a damage function in which the marginal net private profit equals the marginal external cost. In a first-best world, where climate externalities are the only source of inefficiency, equilibrium is achieved by imposing a tax equal to the difference between that profit and cost. This tax rate is the optimal Pigouvian tax (Pigou 1932; Engström and Gars 2015) . However, the recycled revenue effect and the tax interaction effect cause the optimal carbon tax to deviate from the Pigouvian tax (Bovenberg and Mooij 1994; Bovenberg and Goulder 1996; Bovenberg and Ploeg 1998) . Because of the challenges in the calculation of the cost of environmental damage, the Pigouvian tax is difficult to identify in practice. Against this background, Baumol and Oates (1971) proposed an alternative approach called the standards and pricing approach. This approach involves choosing environmental standards first and using environmental taxes to bring the level of environmental damage down to that accepted in the standards. Under the framework of standards and pricing approach, extensions of this literature have established integrated assessment models by introducing climate-related temperature equations into neoclassical growth models (Nordhaus 2007; Hassler and Krusell 2012; Engström and Gars 2015; Wesseh and Lin 2018) . To date, all carbon taxes have been implemented against certain environmental or economic standards rather than based on the Pigouvian tax (Ekins and Barker 2001) . This paper adopts this approach. Regarding the research on the carbon tax rate model for China, scholars have conducted much research, and the main conclusions are as follows. First, a dynamic uniform carbon tax rate is advocated by some researchers (Fan et al. 2016; . These scholars believe that the optimal carbon tax should increase with the GDP growth rate over time. A mild initial tax rate can offset the negative impacts of a carbon tax on the economy (Liao and Zheng 2016) . Second, some scholars believe that a dynamic differentiated carbon tax rate is suitable for addressing the diversity across regions and sectors. Zhang et al. (2019) proposed that the optimal carbon tax model for China is a dynamic differentiated one, following the principle of a low initial rate that gradually increases. Li et al. (2012) noted that a differentiated carbon tax is better than a uniform carbon tax rate for steel because of the lower abatement costs. Third, some scholars have focused on the fixed carbon tax model in different scenarios (Hu et al. 2021; Ma et al. 2021a; Tan et al. 2021; Weng et al. 2021) . Because the assumptions and methodology vary significantly, the optimal carbon tax rate is found differentiated accordingly. Generally, most studies agree that a lower carbon tax rate cannot effectively lead to cleaner energy consumption, while a higher carbon tax rate on energy-intensive industries can achieve greater environmental improvements (Lin and Jia 2018) . In general, the existing literature on the optimal carbon tax for China has mainly focused on the uniform carbon tax given certain environmental standards, and few studies have paid attention to the regional differentiated carbon tax. However, the research on industrial carbon taxes is insufficient. Second, few studies have comprehensively investigated the two dimensions of the carbon tax rate in one analysis: uniform carbon tax vs. differentiated carbon taxes and dynamic carbon taxes vs. fixed carbon taxes. Third, compared with the research in the fixed carbon tax, the studies of the dynamic carbon tax are scarce. As a result, the key information for the design of a carbon tax is relatively inadequate. This study adopts the method proposed by Li et al. (2019) and Song et al. (2019) . The model is constructed based on balanced relationships in the I-O table, including the horizontal balance and the vertical balance. The interactive relationships among all agents are shown in Fig. 1 . All industries are engaged in production, consumers provide the factor inputs and obtain income, the government collects taxes from producers and consumers, and consumers and the government spend their income on consumption and investment. A balance exists between the supply and demand of materials and values. CO 2 emissions are caused by production and consumption activities. This model includes two parts: objective functions and constraints. There are three objectives: the maximization of GDP, minimization of carbon emissions, and minimization of the carbon tax burden. Given the balanced relationships in the I-O table, the following constraints exist: balance in material flows, balance in value flows, balance between the income and expenditure of institutions, balance between investment and savings, the target reduction in carbon intensity, and requirements for GDP growth, among others. This paper has extended the reference model as follows: (1) an optimal dynamic tax block is added to explore the optimal carbon tax rate problem subject to a given environmental standard and (2) in the dynamic block, technological improvements and capital flow across industries are consolidated. The assumptions of this model are as follows: (1) a fixed production structure for industries, which means the direct requirement coefficient is fixed and (2) constant returns to scale, and (3) the carbon cost can be fully transferred to the downstream industries. The terms can be found in Table. 8 in Appendix 1. The economic objective of the government is to maximize economic growth ( GDP t ) subject to certain economic and environmental constraints. The formula can be written as follows: where v j denotes the value-added coefficient of industry j , P j,t denotes the price rate of industry j , X j,t denotes the output of industry j . t represents the year. 2) Minimization of carbon emissions. To achieve sustainable economic development, the environmental objective of the government is to minimize the total amount of carbon emissions ( TCO t ) pro- duced for a given level of output. The formula can be written as follows: where SCO i,t denotes the carbon emissions of industry i ,ECO i,t denotes the carbon emissions of end-users. 3) Minimization of the carbon tax burden. According to the principle of profit maximization, the optimal carbon tax rate is defined as the tax rate that makes each industry achieve its emission reduction target with the minimum extra carbon tax cost. The carbon tax burden ( ETB t ) is defined as the ratio of the total carbon tax ( TCTAX t ) to GDP, according to the study of Wang and Liu (2005) . The formula can be written as follows: (1) Balance in material flows. In line with the horizontal balance in the I-O table, the demand for commodities and services must equal the supply. For industry i , total output equals or is greater than the sum of intermediate inputs and final demand. where A denotes the direct requirement coefficient matrix, t denotes the factor-productivity coefficient, HD i,t denotes the household consumption of the i-th good, GD i,t denotes the government consumption of the i-th good, INV i,t denotes the capital formation for the i-th good, EM i,t denotes the exports from the i-th good, and IM i,t denotes the imports to the i-th good. (2) Balance in value flows. In line with the vertical balance in the I-O table, the revenue of each industry must equal its costs. For industry j , total output and subsidies from the government ( subsidy j,t ) are equal to the sum of intermediate inputs, labor compensation, fixed asset depreciation, direct taxes, and the carbon tax. where tacu j,t denotes the carbon tax rate per unit of output in industry j , ratehk j denotes the income rate for industry j , depr j,t denotes the rate of depreciation in industry j , K j,t denotes the capital stock of industry j , tind j,t denotes the indirect tax rate in industry j , and SCTAX j,t denotes the carbon tax of industry j. (3) Balance between the income and expenditure of households and the government. a. Income and expenditure of households. Household income is the sum of labor remuneration and the operating surplus from each industry (Eq. 6). Total household income ( YHT t ) is spent on direct taxes ( GHTAX t ) (Eq. 7), household savings ( SAVH t ) (Eq. 8), and household consumption (Eq. 9). The share ( conh i ) of all goods in household consumption is fixed according to the I-O structure. Total household income equals total household expenditure (Eq. 10). (4) where t h denotes the direct tax rate and s h denotes the household savings rate. b. Income and expenditure of the government. Government income ( YGT t ) consists of indirect taxes ( GINDTAX j,t ), direct taxes, and the total carbon tax (Eq. 11). It is spent on government savings ( SAVG t ) (Eq. 12), emission abatement subsidies (Eq. 13), and government savings (Eq. 14). Total government income must equal total government expenditure (Eq. 15). where s g denotes the government savings rate, tcsr denotes the abatement subsidy rate, and cong i denotes the share of the i-th good in government expenditure. (4) Balance between investment and savings. Total savings is sourced from household savings, government savings, and fixed asset depreciation. The fixed asset depreciation of one industry is linearly related to its capital stock through the rate of depreciation. Total savings is spent on net exports, inventory, and new investments. (5) Carbon emissions and carbon tax. a. Carbon emissions. Carbon emissions, including industrial and household carbon emissions, are caused by fossil fuel consumption (Eq. 17). The industrial carbon emissions coefficient is the rate of industrial carbon emissions per unit of output ( ECI i ). For details on the calculation process, please refer to the study of Ma et al. (2021b) . The industrial carbon emission coefficient is affected by improvements in energy efficiency ( ). National carbon intensity ( N_intensity t ) is the rate of total carbon emissions per unit of GDP (Eq. 18). b.Carbon tax. According to a study conducted by the Ministry of Finance, the initial taxation stage for the carbon tax in China is proposed to be set at the production level. Then, for industry j , the carbon tax can be written as in Eq. 19. (6) Carbon intensity. According to the relevant national policies, the reduction in carbon intensity must meet the set target. where u denotes the upper bound on the national carbon intensity and l denotes the lower bound on the national carbon intensity reduction. (7) Other constraints. a. Steady growth in imports and exports. Because the proportion of imports and exports to total output is relatively low, this paper assumes that it is linearly related to total output. The formulas can be written as follows: where emi i denotes the rate of export for the i-th good and imi i denotes the rate of import for the i-th good. b. Steady growth in GDP. The yearly GDP growth rate ( rategdp t ) should be maintained above a certain level. c. Steady growth in output. The yearly output growth rate should remain within a certain range. where ratex u denotes the upper bound on the output growth rate and ratex l denotes the lower bound on the output growth rate. d. The range of the acceptable optimal carbon tax. The yearly optimal dynamic carbon tax rate should remain within a certain acceptable range. where u is the upper bound on the optimal carbon tax rate and l denotes the lower bound on the optimal carbon tax rate. (8) Dynamic block. a. Technological improvement. The technological improvement coefficient must increase at a certain rate ( ). b. Capital stock and capital flow. For industry j , the capital stock is linearly related to its total output through the capital production coefficient ( h j ) (Eq. 27). Net investments and the capital stock in year t collectively accumulate as the capital stock in year t + 1 (Eq. 28). Capital is assumed to flow from industries with a lower return on capital (ROI) to industries with a higher ROI until the ROI is constant across industries (Eq. 29 and Eq. 30). The amount of new investment in industry i depends on total new investments and the dynamic share of the investment distribution in that industry ( i,t ) (Eq. 31). where ARK t denotes the average return on capital, R t denotes the return on capital factor, kdist i,t denotes the coefficient of capital distortion, and i denotes the capital flow coefficient for industry i. This paper adopts the NSGA-II to solve this multiobjective optimization model, which is one of the most widely used and most stable evolutionary algorithms. Multi-objective optimization usually results in a Pareto solution set with numerous nondominated solutions instead of only one best solution as in a single-objective optimization problem. Therefore, the multicriteria tournament decision (MTD) method is used to select the most appropriate solution from the numerous Paretooptimal solutions in the set according to the preferences of the decision-maker. MTD is widely used to ascertain the final solution from a Pareto-optimal set for multiobjective decision problems (Parreiras and Vasconcelos 2007; Yu et al. 2015) . Further details on the MTD can be found in Appendix 2. Each objective receives equal weight. Sensitivity testing for the proposed model can be found in Appendix 3. The industries are separated into two groups according to the Chinese energy conservation regulations. Eight industries are identified as energy-intensive industries, such as coal-mining and fuel processing products, while the remaining 23 industries are identified as low-carbon industries. The detailed industry classification is shown in Table. 9 in Appendix 4. The data are mainly derived from the Chinese I-O Lou (2015) . The industrial carbon emission coefficient ( ECI i ) is taken from the research of Ma et al. (2021b) . The range of industrial output changes [ ratex l , ratex u ] is set to [0, 20%] on the basis of the research of Yu et al. (2015) and Jiang (2020) . The abatement subsidy rate ( tcsr ) is set to 20%, following the research of Tang et al. (2015) and Xiao et al. (2020) . Since most scholars suggest a moderate carbon tax at the beginning and few studies propose a carbon tax rate that exceeds ¥200/t (Yi Li et al. 2021a; Ma et al. 2021a; Weng et al. 2021) , the range of acceptable optimal carbon tax rates [ l , u ] is set to [¥20/t, ¥200/t], taken from the study of Zhao and Yin (2016) . The coefficients u and l are set in line with the 13th Five-Year Plan of China. By 2018, national carbon intensity in China had fallen 45.8% since 2005 (Ministry of Ecological and Environment PRC 2019). To meet the target reduction in carbon intensity, national carbon emissions intensity should further decrease by 19.2-24.2% from its 2005 level. Accordingly, the annual average rate of decrease is set to be within (1.6%, 2.02%). According to the two dimensions for setting tax rates, the sub-models include the fixed uniform tax rate, the fixed differentiated tax rate, the dynamic uniform tax rate, and the dynamic differentiated tax rate. The scenarios are designed based on the sub-models. For detailed scenario descriptions, please refer to Table 1 . The baseline case without an industrial carbon tax is the BAU scenario. In line with the relevant previous literature, the rate of technological improvement ( ) is set at 2% (Lu and Chen 2020) , and the yearly energy efficiency improvement rate ( ) is set at 1.5% (Wei and Ma 2015; Zhang and Ma 2020) . The forecasted GDP growth rate from 2021 to 2030 ( rategdp t ) is taken from the research of CASS (2020). Based on related government reports and the previous literature, the fixed uniform carbon tax rates are set to ¥40/t (Li and Yao 2020), ¥50/t (Yi Li et al. 2021a ) and ¥60/t (Dong 2020) . The differentiated carbon tax in the European countries that have implemented a carbon tax is usually set to 50% across industries (Su et al. 2009 ). Therefore, under the fixed differentiated tax rate, we increase the carbon tax rate by 50% for energy-intensive industries. Thus, the tax rates for low-carbon industries are set to ¥40/t, ¥50/t, and ¥60/t, while the tax rates for energy-intensive industries are set to ¥60/t, ¥75/t, and ¥90/t. The optimal dynamic uniform and differentiated tax The optimal dynamic uniform and differentiated carbon tax rates under given carbon intensity constraints from scenarios C1 and D1 for 2021-2030 are shown in Fig. 2 . From 2021 to 2022, there is no tax because this stage still allows the national carbon intensity reduction target to be met. Starting in 2023, the optimal tax rates in both scenarios increased over time. The optimal dynamic uniform tax rate increases from ¥35.84/t to ¥166.54/t, while the optimal dynamic differentiated tax rate for energy-intensive and low-carbon industries increases from ¥42.21/t and ¥28.14/t, to ¥182.57/t and ¥121.71/t, respectively. For the dynamic tax rate, the tax implementation date in both scenarios is relatively late, and the initial tax rate level is low, beginning in the range [¥28/t, ¥60/t] in 2023-2024. This leads to weak carbon intensity constraints on industries, which could not effectively guide industries to begin consuming cleaner energies. As a result, carbon emission levels of all industries are high in the initial stage. Beginning in 2026, the dynamic tax rate rises rapidly, even becoming several times higher than the fixed tax rate by 2030. This demonstrates that if the industrial carbon tax is promulgated later together with a lower initial tax rate, although GDP growth can be maintained at a higher level in the beginning, such a tax rate level cannot effectively guide industries to choose cleaner energy consumption patterns. Therefore, to meet the target reduction in carbon intensity in a short time, more of the pressure on industries is transferred to a later stage. All this causes a relatively higher tax rate to become necessary. The optimal dynamic uniform tax rates are found to be higher than the optimal dynamic differentiated tax rates in most cases from 2021 to 2030. The average carbon tax burdens in the two scenarios are 1.0978% and 0.9102%, respectively (see Table 2 for details). This is due to the differentiated carbon tax model emphasizing that different polluters should pay carbon taxes according to their damages, and high-emission industries should bear additional tax costs. In this way, negative externality spillovers can be prevented, and tax equity can be reflected to some extent. The total tax burden of the differentiated carbon tax model is smaller. As Fig. 3 shows, total carbon emissions in China continue to rise in all scenarios from 2021 to 2030. The carbon emissions in scenario BAU increase from 9.69 billion tons to 19.96 billion tons, with an average annual growth rate of 7.55%. Given that the industrial carbon tax was introduced, the increase in carbon emissions is effectively slowed down in all scenarios. It is shown that there are lag effects in the industrial carbon tax and the energy transition cannot complete in a short time. In the long run, the policy impacts of the industrial carbon tax have been gradually enlarged, which leads to carbon emission reductions being more significant. In addition, it is found that the carbon emissions level of the dynamic tax is higher than that of the fixed tax before 2026, which can be explained by the lower tax rate in the beginning. The carbon emission reductions during 2021-2030 in all scenarios are illustrated in Table 3 . From 2021 to 2030, the carbon emissions' reduction in all scenarios gradually rises over time. The carbon emission reductions in the differentiated tax are larger than that of the uniform tax. That means the more tax is levied on energyintensive industries, the more carbon emissions reduction there is. This result is in line with the findings of Lin and Jia (2018) . Although the tax rates of the dynamic carbon tax are the highest since 2026, the overall carbon emissions are the smallest. This implies that the environmental effect of the dynamic carbon tax is inferior to that of the fixed carbon tax. Figure 4 illustrates the variation in the industrial carbon emissions reductions in all scenarios compared with the BAU scenario in 2030. In all scenarios, the top three contributors to emissions abatement are electricity production, metal smelting, and chemical industry. Carbon emission reductions are mainly driven by output changes and the carbon emission coefficient. The top three industries are all basic industrial commodity providers with high energy consumption and have the largest carbon emissions in the baseline year. Therefore, their carbon emission reductions are larger than those of other industries. The shares of coal consumption in the top three sectors reach 60%, 16.96%, and 6.84%, respectively. For such energy-intensive industries, the industrial carbon tax would lead to high carbon tax costs and provoke a decline in supply. Thus, the carbon emission reductions in these three industries are relatively higher than those of the others. Carbon intensity in all scenarios in 2030 is shown in Table 4 . Although the amount of carbon emissions continues to rise over time, carbon intensity declines. In scenario BAU, carbon intensity falls by 12.27% from 2021 to 2030, decreasing from 1.55 to 1.37 t/¥10,000. When the industrial carbon tax is introduced, the decreases in carbon intensity accelerate. The decreases in scenarios B2, B3, C1, and D1 are the greatest, accounting for 19.25%, 19.35%, 19.20%, and 19.22%, respectively. Each of the above four scenarios meets the minimum target for carbon intensity reductions, while the rest of the scenarios fail to do so. Meanwhile, the decreases in the carbon intensity of the dynamic carbon tax are still not significant. Similar to the effects of carbon emissions, the carbon intensity reductions in the differentiated tax are better than that of the uniform tax. The variation in the carbon intensity reductions of industries in all scenarios relative to the BAU scenario in 2030 is shown in Fig. 5 . In all scenarios, the top three industries with the largest carbon intensity reductions are electricity production, fuel processing, and metal smelting. In scenarios B2, B3, C1, and D1, which meet the carbon intensity reduction target, the top three industries contribute approximately 18%, 17%, and 16%, respectively. However, the decreases in carbon intensity among some of the energy-intensive industries are even smaller than those among some of the low-carbon industries, such as construction and transportation. This is mainly because these two industries can be categorized as indirectly energy-intensive industries, which consume a large amount of energy-intensive intermediate inputs, while their direct energy consumption is smaller. As a result, the direct carbon tax costs are lower, and the direct impacts of the industrial carbon tax are smaller. The average annual GDP growth rate during 2021-2030 is shown in Table 5 . The average annual growth rate of GDP is 5.97% in scenario BAU. All other scenarios cause GDP losses to some extent. In general, the GDP losses from the dynamic tax rates are greater than those from the fixed tax rates, while the GDP losses from the differentiated tax rates are greater than those from the uniform tax rates. Nevertheless, even with the greatest GDP losses in scenario D1, the magnitude of these GDP losses is only − 0.67%. Thus, the negative impacts of the industrial carbon tax are relatively mild. Greater carbon tax rates are found to be associated with more fluctuations in the GDP annual growth rate (see Fig. 6 for details). Among all scenarios, the fluctuations in scenarios C1 and D1 are the largest, at 3.30% and 3.22%, respectively. Because of the late taxation date and the low initial tax rates, the dynamic tax rate exchanges environmental losses for GDP benefits in the early stages, when high levels of emissions are accumulated. To reduce the carbon intensity to the target level in a short period, a high tax rate needs to be imposed on industries. Thus, the high tax rate has large negative impacts on the economy, resulting in the largest fluctuations in GDP growth. In this way, to ensure stable economic growth, the industrial carbon tax should be promulgated with moderate rates and as early as possible, which should encourage all industries to reduce their carbon emissions to an appropriate level. As shown in Fig. 7 significantly. The industries with the largest declines in output are coal mining and electricity production. In scenario C1, the decline in output for coal mining reaches 38.73% in 2030. Obviously, the industrial carbon tax can effectively reduce the coal consumption of industries and guide them to choose cleaner energies. However, except for the above two energy-intensive industries, the declines in output for the remaining energy-intensive industries are not significant. This is due to the large quantity of capital stock distributed across the energy-intensive industries, which is caused by China's extensive development in the past. The overall share of the capital stock in eight energy-intensive industries in the base year is 39.81%. More specifically, the share of the capital stock in transportation reaches 26.26%, and the coefficient of the capital distortion in the construction industry is 38.29. As a result, investments are path-dependence and still flow into energy-intensive industries, which offset the impacts of the industrial carbon tax to some extent. In contrast, the decreases in the output of some low-carbon industries exceed decreases in some energy-intensive industries, such as other manufacturing, repair services, and furniture manufacturing. This is mainly because such industries have the lowest share of investment, accounting for 0.05%, 0.02%, and 0.50% in the base year, respectively. Moreover, the capital production coefficients are as high as 3.86, 3.63, and 2.40. Therefore, a small change in investment could lead to a large change in output. Given this fact, more attention needs to be given to the problem of path-dependence in investments when implementing an industrial carbon tax. As shown in Fig. 8 , the variation in industry prices in all scenarios relative to the BAU scenario in 2030 is significant. In all scenarios, the energy-intensive industries have the greatest price changes. The industry with the largest price change is electricity production, reaching 4.91% in scenario D1. As the third-largest energy provider, electricity production has a wide influence on other industries, so its price change is the largest. Energy-intensive industries, such as fuel processing, metal smelting, and coal mining, provide basic industrial commodities and energy, each of which plays important role in the industrial chain. Thus, their corresponding prices increase significantly. The price changes in the energy-intensive industries under the dynamic carbon tax are the largest; those under the fixed differentiated carbon tax are next, and those under the fixed uniform carbon tax are the smallest. For example, the price increase for electricity production in scenario D1 is 4.78 times higher than that in scenario A1, 3.82 times higher than that in A2, and 3.17 times higher than that in A3. The reason is that the tax rate in scenario D1 exceeded ¥150/t in 2030, a very large amount. We can conclude that the higher is the carbon tax rate, the higher is the commodity price. This implies that the tax burden of the dynamic carbon tax is the largest. The prices of the indirect energy-intensive industries are insensitive to the carbon tax. The price changes in construction and transportation average only 0.02% and 0.26%, respectively. In contrast, some low-carbon industries which have close production relationships with energy-intensive industries have large price increases. For example, the price increase of P&G extraction and metal mining reach 0.49% and 0. 69%, respectively. This paper simulates the environmental and economic effects from 2021 to 2030 of four different tax rate models. In Table 6 , a comprehensive policy evaluation table is established that includes the four dimensions of cumulative GDP losses, cumulative carbon tax burden ratio, cumulative carbon emissions reductions, whether the carbon intensity reduction target is met, and whether the carbon neutrality target is met. The fixed differentiated tax rate model results in greater cumulative losses than the other three models. This model places higher tax rates on energy-intensive industries and achieves relatively greater environmental improvements. This reflects the pollutant-pays principle to an extent and effectively encourages energy transitions. Regarding the dynamic tax rate model, scenario C1 results in the highest cumulative carbon tax burdens, reaching up to 1.0978%. Although the late tax implementation date and low initial carbon tax rate prevent more GDP losses, the high level of carbon emissions in the early years results in a higher tax rate in the later years. Furthermore, the higher carbon tax rate needed to reduce carbon intensity to a target level in the short term causes large fluctuations in GDP growth. The policy benefits of scenarios C1 and D1 are only 1.2108% and 1.3956%, respectively. According to the study of Liu and Hu (2022) , the carbon sink of China in 2030 will increase to 38.23 billion t, which can absorb and store approximately 47.41-67.67 billion tons of CO 2 . Comparing the carbon emissions with carbon sink, all scenarios fail to meet the carbon neutrality target. Regarding the overall policy benefits, we conclude that the differentiated tax rate model is superior to the uniform tax rate model, and the fixed tax rate model is superior to the dynamic tax rate model. Since carbon emissions are mainly produced by energy-intensive industries, a higher carbon tax rate on those industries can not only achieve greater emission reductions but also provide a large degree of tax fairness. On the other hand, compared with the dynamic tax rate model, the fixed tax rate model has milder negative impacts on the economy and imposes a smaller carbon tax burden. Therefore, the fixed differentiated tax rate model has the advantages of both differentiated tax rates and fixed tax rates, so it is the best tax rate model. Compared with all the other scenarios, the overall policy benefits of scenario B2, which impose a tax of ¥50/t on lowcarbon industries and ¥75/t on energy-intensive industries, are the largest. s1 s2 s3 s4 s5 s6 s7 s8 s9 s10 s11 s12 s13 s14 s15 s16 s17 s18 s19 s20 s21 s22 s23 s24 s25 s26 s27 s28 s29 s30 s31 Carbon intensity reduction (tons/10 thousand yuan) Industry A1 A2 A3 B1 B2 B3 C1 D1 comparisons with other studies are listed in Table 7 . We compared the critical indicators of models, including carbon emissions, carbon intensity reduction, and GDP. All three results of this research are smaller than those of Li et al. (2019) because we place limitations on the GDP growth rate according to the study of Lou (2015) . The volume of carbon emissions and carbon intensity in 2030 in our research are greater than the studies of Lin and Jia (2018) and Niu et al. (2020) . These differences are because that there is no limitation on the labor and capital factors supply, which results in a greater output. This limitation will be improved in future research. Regarding the dynamic optimal carbon tax rate, the carbon tax burdens in scenarios C1 and C2 in 2030 are 1.73% and 1.82%, respectively, which are similar to 2% in the study of Fan et al. (2016) . To investigate the optimal industrial carbon tax rate for China, this paper simulates the environmental and economic effects of different scenarios from 2021 to 2030 via an integrated dynamic I-O model. The main conclusions are as follows. (1) Given the same carbon intensity constraints, the optimal tax rate model is the fixed differentiated tax rate. The fixed differentiated carbon tax can not only lead to greater environmental improvements but also provide a large degree of tax fairness. The largest policy benefits in all scenarios are from scenario B2, 2.1256%, where the tax is set to ¥50/t for low-carbon industries and ¥75/t for energy-intensive industries. (2) The existing investment structure is path-dependence and offsets the policy effects of a carbon tax to some extent. A large amount of capital stock is distributed across the energy-intensive industries due to China's energy-intensive development. As a result, investments are path-dependent and still flow to energy-intensive industries, which offset the impacts of the industrial carbon tax. (3) Indirect energy-intensive industries such as construction and transport are insensitive to the industrial carbon tax. Such industries consume a large amount of energy-intensive intermediate inputs, while their level of direct energy consumption is low. As a result, the output and price changes in indirectly energy-intensive industries are not significant. Based on the above results and discussions, we present the following policy suggestions for setting the industrial carbon tax rate: (1) As a supporting environmental tool of Chinese ETS, the fixed differentiated carbon tax is proposed to be introduced at a moderate rate step by step. The high fluctuation of the carbon price in the ETS market causes the expected return on investments in lowcarbon technology upgrading to be uncertain, which is unfavorable to carbon neutrality. In contrast, the fixed differentiated carbon tax is beneficial to a low overall carbon tax burden and the stability of economic growth. Furthermore, it can reflect the pollution-pays principle to the greatest extent. It is suggested that the fixed differentiated carbon tax should be promulgated from loose to tight and adjusted according to the policy performance by stages. Some energy-intensive industries, such as electricity production and metal smelting, can be selected as test pilots for the fixed carbon tax. Based on experimental experiences, the fixed differentiated carbon tax can be promoted to the whole country. (2) More attention should be paid to improving the existing investment structure to a cleaner one. More capital should be guided to the equipment investments on the key infrastructure of emerging information industries, such as cloud computing, big data, mobile internet, and artificial intelligence. This will improve the production efficiency and product quality of the existing production system. Meanwhile, to attract more long-term capi-tal to enter into the high-tech industries, the tolerance of investment failures in high technology in financial institutions should be improved. (3) For indirect energy-intensive industries, technology standards other than carbon taxes should be adopted for regulation. For construction, to reduce energy inputs from the source, high durable and green materials are first to be used to prolong the service life of buildings. In addition, whole life-cycle management is suggested to cover the stages of planning, design, production, and operation of buildings, which can reduce the waste in reworking. For transport, it is suggested to speed up changes in transportation models, from the highway to railway and from the highway to water in mediumand long-distance transportation; on the other hand, the application of new energy should be accelerated. The limitations of this article are as follows: (1) Because there are no constraints on the labor and capital factors, the short-term shock of Covid-19 on the labor supply and household consumption cannot be captured. (2) The heterogeneity of key industries with priority development in future national strategy is not concluded. All the above limitations will be improved in future research. Appendix 1 The multicriteria tournament decision method The MTD is a widely used decision method that provides a ranking of all alternatives in a Pareto-optimal set by a tournament function, according to the preferences of the decision-makers (Parreiras and Vasconcelos 2007; Yu et al. 2015) . The detailed process is as follows. First, for a given objective in the model, the number of times that alternative a wins the tournament against each other solution b from the Pareto-optimal set U can be described as follows: where objective s is a tournament. fit s (a) represents the value of the objective function when alternative a is chosen. fit s (b) represents the value of the objective function when alternative b is chosen. Second, the tournament function Ts (a,U) provides the ratio of the number of times that alternative a wins the tournament against other solutions in the Pareto-optimal set U. The formula can be written as follows: where U represents the entire Pareto-optimal set. Finally, to aggregate the scores into a global ranking function, all criteria and their weights are considered. The aggregation function can be evaluated using the following formula: where R(a) denotes the total score for alternative a in all tournaments. w s denotes the weight of objective s , with w s > 0, and ∑ o s=1 w s = 1 . o denotes the total number of objectives. Sensitivity testing is necessary, as some key parameters in the model may significantly affect the robustness of the simulations. A sensitivity analysis is performed to (33) T s (a, U) = ∑ ∀b∈U,a≠b t s (a,b) (U−1) determine how different values of certain key independent variables affect the simulation results. This model represents the multi-objective optimization problem of setting the industrial carbon tax rate subject to certain carbon intensity reduction constraints. Thus, the carbon intensity reduction constraints and the weights on the objectives are the key parameters of the model. Since most of the other parameters adopted in the model have already been tested in the existing literature, only the carbon intensity reduction constraints and objective weights are tested in this paper. Given that the simulation results mainly focus on the environmental and economic effects of the industrial carbon tax, the key indicators are CO 2 emissions reduction, carbon intensity reduction, and the GDP growth rate. Based on the above, the sensitivity analysis investigates the effect of the carbon intensity reduction constraints and the objective weights on the above three key indicators. Following the study of Hosoe et al. (2010) , the following testing standards for the sensitivity analysis are proposed: (1) A variable attribute is considered to be stable if the sensitivity analysis results are stable, that is, if there is no change from a positive to negative coefficient.⋅⋅ (2) There should be no significant change in the order of importance within the same group of variables. The changes in the key indicators when the carbon intensity reduction constraints change by + 0.1% or − 0.1% are shown in Fig. 9 . The changes in the key indicators when the objective weights change by + 10% or − 10% are shown in Fig. 10 . None of the key indicators in either test changed from positive to negative and all showed strong stability. We can conclude that the key parameters pass the sensitivity testing and that the settings for these key parameters are reasonable.⋅⋅ The changes in the key indicators from the simulation due to different carbon intensity constraints -0.30% -0.20% -0.10% 0.00% 0.10% 0.20% Carbon intensity (yuan/ton) Rate of change -10% 10% Fig. 10 The changes in the key indicators from the simulation due to different weights on the objectives Table 9 Industry classification, industry codes Industry classification Industry code Industry Industry classification Industry code Industry Energy-intensive industries S2 Coal mining Low-carbon industries S10 Education-related products S11 Fuel processing S15 Metal products S12 Chemical industry S16 General machinery S13 Nonmetallic products S17 Special machinery S14 Metal smelting S18 Transport equipment S25 Electricity production S19 A comparison of pure manufacturing and hybrid manufacturing-remanufacturing systems under carbon tax policy Dynamic simulation for comprehensive water resources policies to improve water-use efficiency in coastal city The use of standards and prices for protection of the environment. The economics of environment Optimal environmental taxation in the presence of other taxes: general-equilibrium analyses Environmental levies and distortionary taxation Consequences of environmental tax reform for unemployment and welfare China's emissions trading system and an ETS-carbon tax hybrid The short-term impact on emissions and federal tax revenue of a carbon tax in the U.S. electricity sector Research on the potential growth of the chineses economy in the next 15 years and the main goals and indicators of economic and social development during the 14th five-year plan period Towards carbon neutrality by implementing carbon emissions trading scheme: policy evaluation in China The effects of green growth, environmental-related tax and ecoinnovation towards carbon neutrality target in the US economy Commentator of Economic Daily (2022) Overall plan to steadily promote carbon peak and carbon neutralization Comparison of carbon emission reduction between carbon tax and energy efficiency enhancement-based on the simulation analysis of CGE model Carbon taxes and carbon emissions trading Optimal taxation in the macroeconomics of climate change Externalities of dynamic environmental taxation, paths of accumulative pollution and long-term economic growth Economics and climate change integrated assessment in a multi-region world Textbook of computable general equilibrium modeling: programming and simulations A comparative study on the environmental and economic effects of a resource tax and carbon tax in China: analysis based on the computable general equilibrium model Rethinking the choice of carbon tax and carbon trading in China Tax included computable nonlinear dynamic input-output model and determination of its parameters The industrial structure optimization for minimizing carbon emissions based on the global industrial chain network perspective. Dissertation The study of carbon dioxide emission intensity abatement mechanism of iron and steel industry based on two-stage game model Assessment of the carbon emissions reduction potential of China's iron and steel industry based on a simulation analysis Research on "Double Dividend" of the carbon tax policy based on CGE model. The Theory and Pract of Fin and Econ A quantitative roadmap for China towards carbon neutrality in 2060 using methanol and ammonia as energy carriers. iScience Can energy supply-side and demand-side policies for energy saving and emission reduction be synergistic?-a simulated study on China's coal capacity cut and carbon tax Environment, life expectancy and economic development: the optimal environmental tax The energy, environmental and economic impacts of carbon tax rate and taxation industry: a CGE based study in China China builds the world's largest forestry carbon sink market (2020-2050) The theory and application of Chinese economy-energyenvironment-tax dynamic computable general equilibrium model The scenarios of China reaching the peaking of carbon emission before 2030 and its macroeconomic implications The nexuses between energy investments, technological innovations, emission taxes, and carbon emissions in China The shortterm price effects and transmission mechanism of CO 2 cost passthrough in China: a partial transmission model Opinions on innovation and improvement of price mechanism for promoting green development Can China achieve its 2030 carbon emissions commitment? Scenario analysis based on an improved general regression neural network To tax or not to tax: alternative approaches to slowing global warming A multiplicative version of Promethee II applied to multiobjective optimization problems Exploring potential pathways towards fossil energy-related GHG emission peak prior to 2030 for China an integrated input-output simulation model Unraveling economic and environmental implications of cutting overcapacity of industries: a city-level empirical simulation with input-output approach Carbon emission control policy design based on the targets of carbon peak and carbon neutrality Performance evaluation of remanufacturing strategies under uniform carbon tax and discriminated carbon tax Carbon emissions trading scheme exploration in China: a multi-agent-based model Judgment on the level of macroeconomic tax rate in China Production and emissions reduction decisions considering the differentiated carbon tax regulation across new and remanufactured products and consumer preference The impact of SO 2 emissions trading mechanism and sulfur tax on air pollution governance Impact of differentiated industrial carbon taxes on China's economy and environment Optimal carbon taxes for China and implications for power generation, welfare, and the environment The impacts of carbon tax on China's macro-economy and the development of renewable energy power generation technology: based on CGE model with disaggregation in the electric power sector Can China realise its energysavings goal by adjusting its industrial structure? Why should carbon neutrality be incorporated into the construction of ecological civilization: a theoretical explanation and its policy implications Research on enterprise carbon emission decision considering the substitution effect of stepped carbon tax and carbon trading Carbon emission reduction investment for coal-fired power plants and policy analysis considering carbon price floor Emission reduction effect and emission reduction cost of haze by environmental tax-a dynamic multi-region CGE analysis Bi-level programming model of carbon emission reduction decision-making under the coexistence of carbon trading and carbon tax Author contribution Ning Ma: conceptualization, methodology, software, writing-original draft preparation. Guangwei Yin: investigation, writing-reviewing and editing. Huajiao Li: conceptualization, writing-original draft preparation. Wenli Sun: resources, supervision. Ze Wang: resources, investigation. Gang Liu: resources, supervision. Dan Xie: writing-reviewing and editing.Funding This research is supported by grants from the National Natural Science Foundation of China (Grant No. 41871202, No.71991481 and No.71991480), the Fundamental Research Funds for the Central Universities (Grant No. 3-7-6-2021.Data availability All data generated or analyzed during this study are included in this published article. Ethics approval and consent to participate Not applicable. The authors declare no competing interests. Supplementary Information The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s11356-022-19162-6.Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.