Carbon Sequestration Total Factor Productivity Growth and Decomposition: A Case of the Yangtze River Economic Belt of China

To find out whether carbon sequestration is both effective at mitigating climate change and promoting economic growth, in this paper, by adopting a stochastic frontier panel model with translog production function, carbon sequestration is incorporated into endogenous variables to establish estimation model of carbon sequestration total factor productivity (CSTFP) and examine CSTFP growth and its drivers decomposition of the Yangtze River Economic Belt (YREB) of China in three estimations. The result shows that, (1) compared to traditional TFP growth, CSTFP growth in YREB is improved by 26.74 percentages (from −26.55% to 0.20%), contributed by three positive drivers of technical efficiency change (28.59%), technological progress change (18.55%), and scale efficiency change (3.99%); (2) different CSTFP growth exists in three watershed segments of YREB, which firstly is the upper reaches (0.62%), then the lower reaches (0.11%) and the middle reaches (−0.14%). Improved CSTFP growth owes to carbon sequestration’s harmonious symbiosis where natural ecosystems and human activities are naturally blended while insufficient synergies are bottleneck for promotion of CSTFP growth in YREB. Related policy suggestions are provided in the end. The proposed analysis framework is efficient to disclose CSTFP growth in YREB, and can also be applied to similar analysis on CSTFP in regions and extended to multi-country/region analysis.


Introduction
Facing challenges of global climate change, carbon sequestration can endogenously drive both economic development and mitigate climate change. CO 2 emissions reduction and carbon sequestration are both effective measures for dealing with dilemmas between promoting economic development and mitigating climate change [1,2]. Compared to CO 2 emissions reduction in reality, carbon sequestration can fundamentally promote win-win strategy of economic development and mitigate climate change by two ways. One is carbon sinks and reservoirs of CO 2 through photosynthesis, such activities as sustainable forest management practices, afforestation and reforestation, and sustainable forms of agriculture [2]. Because, photosynthesis can make the plant convert CO 2 and water into organic material and release oxygen through biochemical process under visible light irradiation, with light and dark reaction and the use of photosynthetic pigment. As the basis of the biosphere's survival, photosynthesis is the main driver of the global carbon cycle, closely linked to the climate system through numerous feedback processes [3,4], which are associated with extreme weather events, such as heat waves, droughts, and floods [5,6]. More importantly, photosynthesis can increase crop yields and plant products production for the needs of population and economic growth [7,8]. The other way is artificial carbon capture and use through activities of technical innovation and clean development, such as enhancement of energy efficiency in relevant sectors of the national economy, new and renewable forms of energy, research and development (R&D), and increased use of CO 2 sequestration technologies and of advanced and innovative environmentally sound technologies [2].
Carbon sequestration is efficient for improvement of productivity growth. On the aspect of natural reproduction, carbon sequestration creates gross primary production, the absorption of atmospheric CO 2 by plants through photosynthesis, raises ecosystem productivity and contributes to mitigating global warming [9,10]. On the aspect of human activity, carbon sequestration stimulates human productivity and increases output value by offering abundant natural resources for human production activity and economic development. Forest carbon sinks, for example, promote efficiency of ecological economy and green total factor productivity of forests by supply of forest coverage, fresh air, green environment and wood, medicine, mushrooms, and other human production and living necessities [1,11,12]. Especially under the challenge of the green paradox where carbon tax or subsidies to carbon capture, use and storage (CCUS) technology cause higher global emissions through stimulating fossil fuel production [13][14][15], to disclose the endogenous characteristics of carbon sequestration to total factor productivity (TFP) growth and its driver decomposition is extremely inevitable and valuable for sustainable economic development and environmental conservation.
Carbon sequestration as an endogenous variable to TFP growth is in need of urgent discussion. For research of source of productivity growth, scholars used classical production function of the Solow residuals to explore TFP growth by introducing capital, labor, and gross domestic product (GDP) into consideration [16][17][18][19][20]. In the face of increasingly environmental problems caused by CO 2 emissions, regarding CO 2 emissions as negative residual output, together with positive output such as GDP, some researchers took energy consumption and environmental impacts such as CO 2 emissions as endogenous variables into the production function, constructed estimation model to examine green TFP growth [21][22][23][24][25][26][27][28][29][30][31][32][33][34][35]. Basically, two lines of researches on green TFP growth exist in literature. One line is mainly focused on estimating the industrial green TFP growth of fossil energy-intensive industries and sectors, such as estimation of green TFP of Swedish and Canadian pulp and paper industries [21], investigation of green TFP of China's industrial sectors [22][23][24][25][26][27][28], and measurement of ecological economic efficiency and TFP of forests in China [1]. The other line is focused on green TFP growth at the national or regional level. Examples are the exploration on environmentally sensitive productivity growth of 41 developing and developed countries [29], green and environmental TFP measure of China's provinces and regional economies [30][31][32][33], the drivers of TFP change with an illustrative example of 14 EU countries [34], and aggregate green productivity growth estimation at the aggregate level in countries from the Organization for Economic Co-operation and Development (OECD) [35]. Carbon sequestration's endogenous feature to TFP growth is ignored, although a few scholars discussed the functions of forest carbon sinks on TFP of forest and economic growth [1,11,36], however, at the aggregate level, carbon sequestration is not introduced as an impact factor to measure the TFP growth in the literature; and again, most studies on TFP are based on perspective of administrative area, without consideration of watershed geographic correlations of TFP growth in natural river basin regions.
Since carbon sequestration is effective for dealing with dilemmas between economic growth and environmental protection both by carbon sinks and reservoirs of CO 2 through photosynthesis and by artificial carbon capture and use through activities of technical innovation and clean development, it is of great necessity to discuss the role of carbon sequestration as an endogenous variable of TFP growth especially in natural river basin regions by estimation of carbon sequestration total factor growth (CSTFP) growth.
For assessment of the CSTFP growth, the stochastic frontier analysis (SFA) is adoptable. Up to now, abundant literature has sought various approaches to open the "black box" of the TFP growth; the stochastic frontier analysis (SFA) is one of the main approaches in the application of estimating TFP growth. SFA was originally proposed by Aigner et al. and Meeusen and van den Broeck [37,38], and has been broadly applied by many researchers into study [39][40][41][42][43][44][45][46][47][48][49][50][51][52][53][54] due to the following fundamental advantages. Firstly, it can explain more accurately the theoretical definition of a production function which estimates the maximum output attainable from a given set of inputs, and secondly, the technical efficiency of the sample may be estimated by comparing the observed output with the predicted (or attainable) output [44,46,[55][56][57]. Thus, SFA is appropriate for the investigation of CSTFP growth.
Drivers decomposition of CSTFP growth needs consideration on integration of technology progress, technical efficiency change, factor allocation efficiency change, and scale efficiency change. On the drivers of TFP growth, technical effect and technology progress at the national level are discussed but the structural effect and scale effect at the regional level are ignored in most literature [58][59][60][61][62][63][64][65][66]. Actually, even if all regions are efficient, it is still required to optimize resource utilization for the TFP growth at the regional level through coordination for improvement of input-output mixes [67]. Recognizing the shortage of previous research, some scholars have tried to improve analysis on the structural effect or/and scale effect in decomposition of TFP growth [68][69][70][71], such as a decomposition of TFP indicator into components of technical change, technical inefficiency change, and scale inefficiency change on the agricultural sector at the state-level in the U.S. [68] and measurement of green productivity evolution for 30 OECD countries with decomposition for green productivity growth changes into technological progress, technical efficiency change, and structural efficiency change [35]. But, integration of technology progress, technical efficiency change, factor allocation efficiency change, and scale efficiency change has yet to be discussed in literature. Considering carbon sequestration's harmonious symbiosis where natural ecosystems and human activities are naturally blended in the same regions, drivers decomposition of CSTFP growth requires examination on the base of integration of technology progress, technical efficiency change, factor allocation efficiency change, and scale efficiency change.
Based on the above research findings and shortcomings, this paper adopts the stochastic frontier approach with translog production function to estimate the CSTFP growth and discusses its drivers decomposition with the case study of the Yangtze River Economic Belt (YREB) of China. The main innovative contributions to the literature is as follows.
First, although some scholars investigated TFP growth by introducing energy consumption, carbon emissions, and other variables in the production functions, or discussed the interaction between forest carbon sinks and economic growth, there still lacks an analysis on the impact of carbon sequestration on TFP growth. In this paper, carbon sequestration as an endogenous variable is introduced into the production function to investigate CSTFP growth from the perspective of input-output mixes and discloses its impact on TFP growth by comparison of results in three estimations.
Second, most literature focus on analysis of the industrial TFP or regional TFP growth from the level of administrative areas but neglect the watershed geographic correlations of TFP growth in the same river basin. In this paper, the focus on carbon sequestration's characteristics of collaborative symbiosis between economic development and environment optimization, the difference of CSTFP growth with different watershed characteristics of the same river basin, are disclosed.
Third, decomposition of TFP growth in the current research is analyzed by only the use of three components of technical progress, technical efficiency change and scale efficiency change, or those of technical progress, technical efficiency change, and factor allocation efficiency change. By decomposition of CSTFP growth in this paper, four components of technical progress, technical efficiency change and scale efficiency change and factor allocation efficiency change are integrated together to disclose the drivers of CSTFP growth and their harmonious coexistence in YREB.
The structure of this paper is as follows. In Section 2, the conceptual framework and model are presented. Section 3 presents the study area and data, Section 4 is empirical estimation, Section 5 delivers results, and Section 6 is the discussion. Section 8 is the conclusion with policy implications.

Conceptual Framework
As Figure 1 shows, the conceptual framework of CSTFP growth estimation is set up based on features caused by carbon sequestration's harmonious symbiosis correlations between economic development and environment optimization as follows.
decomposition of CSTFP growth in this paper, four components of technical progress, technical efficiency change and scale efficiency change and factor allocation efficiency change are integrated together to disclose the drivers of CSTFP growth and their harmonious coexistence in YREB.
The structure of this paper is as follows. In Section 2, the conceptual framework and model are presented. Section 3 presents the study area and data, Section 4 is empirical estimation, Section 5 delivers results, and Section 6 is the discussion. Section 7 is the conclusion with policy implications.

Conceptual Framework
As Figure 1 shows, the conceptual framework of CSTFP growth estimation is set up based on features caused by carbon sequestration's harmonious symbiosis correlations between economic development and environment optimization as follows. Firstly, compared to CO2 emissions as undesirable output, carbon sequestration is desirable output which produces watershed harmonious symbiosis where natural ecosystems and human activity circle are naturally blended in economic regions. On the environmental sphere, carbon sequestration creates gross primary production, the absorption and decomposition of atmospheric CO2 by plants through photosynthesis, raises ecosystem productivity, and contributes to mitigating global warming in natural ecosystem of economic regions [9,10]. On the economic sphere, carbon sequestration promotes human productivity and makes output value added by offering abundant natural resources for human activities of production and life in the human activity circle of economic regions [1,11,36].
Secondly, as an endogenous variable with integration of other relative variables such as labor, capital, energy consumption, CO2 emissions, industrial structure, open-up, R&D, and urbanization, carbon sequestration's synergy effect on CSTFP growth is showed from the perspective of inputoutput mixes. On the end of input, carbon sequestration stimulates factor investment such as labor, Human Activity Natural Ecosystem Variables Figure 1. The conceptual framework of carbon sequestration total factor productivity (CSTFP) growth estimation.
Firstly, compared to CO 2 emissions as undesirable output, carbon sequestration is desirable output which produces watershed harmonious symbiosis where natural ecosystems and human activity circle are naturally blended in economic regions. On the environmental sphere, carbon sequestration creates gross primary production, the absorption and decomposition of atmospheric CO 2 by plants through photosynthesis, raises ecosystem productivity, and contributes to mitigating global warming in natural ecosystem of economic regions [9,10]. On the economic sphere, carbon sequestration promotes human productivity and makes output value added by offering abundant natural resources for human activities of production and life in the human activity circle of economic regions [1,11,36].
Secondly, as an endogenous variable with integration of other relative variables such as labor, capital, energy consumption, CO 2 emissions, industrial structure, open-up, R&D, and urbanization, carbon sequestration's synergy effect on CSTFP growth is showed from the perspective of input-output mixes. On the end of input, carbon sequestration stimulates factor investment such as labor, capital, and technology into primary productive activities of afforestation, agriculture, forestry, and ecological restoration of rivers, grasslands, and wetlands, as well as human productive activities of manufacturing and services, resulting in the superposition of relevant elements of input. On the end of output, carbon sequestration makes CO 2 absorbed through photosynthesis of nature to release oxygen or clean production to capture and use CO 2 , providing diversified products such as timber, food, and other raw materials for the human society, promoting integration of natural reproduction and social reproduction synchronously.
Thirdly, focus on the target of green and low-carbon sustainable growth of economic regions, the growth of CSTFP is estimated by the stochastic frontier panel model with translog production function with introduction of relative core variables and exogenous variables, and the four drivers of technical progress, technical efficiency change, scale efficiency change, and factor allocation efficiency change are integrated to drive the CSTFP growth for the sustainable development of economic growth and environmental conservation. Thus, is disclosed that the growth of CSTFP and its four-driver linkage naturally promote the harmonious co-prosperity of the human activity circle and natural ecosystem in economic regions.
Based on the conceptual framework of CSTFP growth estimation, the relative methodology is discussed in detail as follows.

Stochastic Frontier Analysis
For estimating CSTFP growth in YREB, stochastic frontier analysis (SFA) is introduced in this paper. Actually, SFA is originally used to study production optimization [37][38][39]. The models of SFA adopt the compound residual measure, and can simply be expressed as follows: where Equation (1) is stochastic frontier production function; i means the provincial region, i = 1, 2, . . . , n. t means time, t = 1, 2, . . . , n; y it is Logarithm of actual local gross domestic product for province i in period t; f [X it (t), β] is Frontier production function; X it (t) is Production factor combination; β is Coefficient to be estimated; (υ it − µ it ) is Compound error structure, supposed µ it and υ it is independent, and it is assumed that these two kinds of errors are not related to X it (t); if (υ it − µ it ) > 0, the sum of noise effect and technical efficiency is positive, indicating that the output is above the determined value of production frontier. Conversely, below the determined value of production frontier. υ it is random error terms assumed to be independent and identical distribution. Equation (2) shows that the distribution of random error term µ it follows a semi-normal distribution and is a non-negative random variable. Equation (3) indicates the distribution of the technical invalidity rate term. Equation (4) is the mean value equation of technical inefficiency. z it is a set of exogenous explanatory variables affecting technological inefficiency. δ is coefficient of explanatory variable of the technical inefficiency equation. If coefficient is positive, the variable has a positive effect on the inefficiency term. Equation (5) is technical invalidity rate term variance equation, θ is the value to be estimated of the coefficient of variance of technical inefficiency term. Equation (6) is stochastic error term variance equation, λ is the value to be estimated of the coefficient of variance of random error term.
Where σ 2 = σ 2 υ + σ 2 µ represents the variance of the compound residual term which is (υ it − µ it ). Battese [42,43] set variance parameter as γ = σ 2 µ /σ 2 (γ ∈ [0.1]) in order to verify the reasonableness of adopting SFA. If γ is close to 1, it is reasonable and necessary to adopt SFA model. If γ is close to 0, the least square estimation can be used.

Translog Production Function Setting
When constructing the SFA model to measure the CSTFP growth, the appropriate approach will be used to test it. Meanwhile, the production function should be satisfied with conditions that any one of the four input variables as labor input (L), capital input (K), energy consumption (E), and carbon sequestration (CS) is affected by the other three input variables, as well as the formation of scale effect and factor allocation effect. Compared with the Cobb-Douglas production model, Leontief's production model and the constant elastic production model, translog production function not only eliminate substitution effect and interaction between input elements, but also reflect influence of time variation, so it is more general [37][38][39][40][41][42][43][44][45][46]. The errors caused by the function can be effectively avoided [37][38][39][40][41][42]. Therefore, translog production function with time variable as specific form of the frontier production function is adopted in this paper, that is:

Technical Inefficiency Function Setting
For the purpose of investigating influence of factors of second industry ratio of (CI) to (SI), urbanization rate (UE), R&D (RD), opening outside (OP) on CSTFP growth, the mean value equation of technical inefficiency is introduced on the basis of (7) as the following [42,43]: where µ it is a random error term.

Decomposition of CSTFP Growth
Based on the relative analysis on TFP decomposition [46,56,57], in this paper, the growth rate of CSTFP can be decomposed into four components: technical progress (∆TP), technical efficiency change (∆TE), scale efficiency change (∆SE), and factor allocation efficiency change (∆AE).

Rate of Technical Progress
The rate of technological progress refers to the change rate of output varying over time with that the relevant input is constant. It shows moving up or down of the production frontier curve along with change of frontier technological progress, as follows [46,56,57]:

Rate of Technical Efficiency Change
It is the change rate of the ratio of actual output to frontier output with time under premise that the cutting-edge technology level is constant, and that of actual output to frontier output with given factors input. It reflects the change in proportion of output over time that the actual efficiency of cutting-edge technology has failed to reach. It measures capacity of producers to use technology, as follows [46,56,57]:

Rate of Scale Efficiency Change
It refers to the change of output corresponding to the increase of factor input with expansion of production scale, which can be in three situations: increasing rate, decreasing rate, and keeping it constant, as follows [46,56,57]: ε l +ε k +ε e +ε cs × ∆K K + ε e ε l +ε k +ε e +ε cs × ∆E E + ε cs ε l +ε k +ε e +ε cs × ∆CS where ε j ( j = l, k, e, cs) represents separately output elasticity of labor, capital, energy consumption, or carbon sequestration; ε is the overall factor scale elasticity, which is the sum of output elasticity of labor, capital, energy consumption, and carbon sequestration. If ε = 1, this part obviously does not have to be added; if ε 1, when the scale efficiency ratio is more than 1, the increase rate of productivity will be raised by increasing the input quantity of factors. But if the scale efficiency is less than 1, the scale of output can only be reduced. .
x j is the change rate of factor input:

Rate of Factor Allocation Efficiency Change
The rate of factor allocation efficiency change mainly examines the contribution of structural change of factor allocation to CSTFP growth in the process of production with given frontier technology level, as follows [46,56,57]: where s j represents the proportion of a certain factor per unit cost, E j E − s j means there is difference between the ratio of output elasticity of a certain factor to unit cost and that of unit cost. When the value is positive, it shows the factor allocation efficiency has been improved, and conversely, the factor allocation efficiency has been worsened.
On the basis of the methodology discussed above, next, let us come to the study area of YREB and related data analysis.

The Yangtze River Economic Belt (YREB) of China
As Figure 2 shows, the YREB is an economic belt with watershed features of harmonious symbiosis where natural ecosystems and human activities are naturally blended. It spans the diversified topography of the western high mountains, the central hilly valley, and the eastern plain delta, it belongs to the subtropical and warm temperate climatic belt and the intersection of the southeast ocean warm current and the northwest plateau cold current with four distinct seasons and rainfall. It is abundant with lush forests, diversified agriculture, modern manufacturing, and a prosperous economy. It is a bridge linked with China's eastern coastal leading economic zone, the western great developing zone, and the central rising economic zone.
The YREB covers 9 provinces of Sichuan, Yunnan, Guizhou, Hubei, Hunan, Jiangxi, Anhui, Jiangsu, and Zhejiang and 2 municipalities of Chongqing and Shanghai. For the sake of comparative analysis, it is divided into three watershed segments: the upper reaches, the middle reaches, and the lower reaches (Yangtze River Delta). The upper reaches include 3 provinces and 1 municipality of Sichuan, Guizhou, Yunnan, and Chongqing; the middle reaches include 4 provinces of Anhui, Jiangxi, Hunan, and Hubei; and the lower reaches include 2 provinces and 1 municipality of Zhejiang, Jiangsu, and Shanghai. The YREB covers 9 provinces of Sichuan, Yunnan, Guizhou, Hubei, Hunan, Jiangxi, Anhui, Jiangsu, and Zhejiang and 2 municipalities of Chongqing and Shanghai. For the sake of comparative analysis, it is divided into three watershed segments: the upper reaches, the middle reaches, and the lower reaches (Yangtze River Delta). The upper reaches include 3 provinces and 1 municipality of Sichuan, Guizhou, Yunnan, and Chongqing; the middle reaches include 4 provinces of Anhui, Jiangxi, Hunan, and Hubei; and the lower reaches include 2 provinces and 1 municipality of Zhejiang, Jiangsu, and Shanghai.
Recently, the YREB is the biggest natural river basin economic belt, accounting for 40% of China's population and GDP as well as CO2 emissions. Since 2000, as shown in Figure 3, carbon sequestration and carbon productivity have increased while energy intensity decreased, along with growth of GDP per capita, economic development correlates ever closer with energy consumption and carbon sequestration. So, it is very valuable to investigate the CSTFP growth in the natural river basin economic belt of YREB.

Variables and Data
In this paper, the growth of CSTFP in YREB is measured based on data of 2000-2016 from statistical yearbooks of 11 provinces and municipalities of the YREB [72,73]. For meeting the needs of Recently, the YREB is the biggest natural river basin economic belt, accounting for 40% of China's population and GDP as well as CO 2 emissions. Since 2000, as shown in Figure 3, carbon sequestration and carbon productivity have increased while energy intensity decreased, along with growth of GDP per capita, economic development correlates ever closer with energy consumption and carbon sequestration. So, it is very valuable to investigate the CSTFP growth in the natural river basin economic belt of YREB. The YREB covers 9 provinces of Sichuan, Yunnan, Guizhou, Hubei, Hunan, Jiangxi, Anhui, Jiangsu, and Zhejiang and 2 municipalities of Chongqing and Shanghai. For the sake of comparative analysis, it is divided into three watershed segments: the upper reaches, the middle reaches, and the lower reaches (Yangtze River Delta). The upper reaches include 3 provinces and 1 municipality of Sichuan, Guizhou, Yunnan, and Chongqing; the middle reaches include 4 provinces of Anhui, Jiangxi, Hunan, and Hubei; and the lower reaches include 2 provinces and 1 municipality of Zhejiang, Jiangsu, and Shanghai.
Recently, the YREB is the biggest natural river basin economic belt, accounting for 40% of China's population and GDP as well as CO2 emissions. Since 2000, as shown in Figure 3, carbon sequestration and carbon productivity have increased while energy intensity decreased, along with growth of GDP per capita, economic development correlates ever closer with energy consumption and carbon sequestration. So, it is very valuable to investigate the CSTFP growth in the natural river basin economic belt of YREB.

Variables and Data
In this paper, the growth of CSTFP in YREB is measured based on data of 2000-2016 from statistical yearbooks of 11 provinces and municipalities of the YREB [72,73]. For meeting the needs of

Variables and Data
In this paper, the growth of CSTFP in YREB is measured based on data of 2000-2016 from statistical yearbooks of 11 provinces and municipalities of the YREB [72,73]. For meeting the needs of analysis on CSTFP growth and decomposition as in Figure 1, the total variables and data are firstly discussed as the following. The regional actual GDP is used as the index of output variable, the labor (L), capital (K), energy consumption (E), and carbon sequestration (CS) are taken as input variables; and the carbon emission intensity of (CI), the second industry ratio (SI), urbanization rate (UE), R&D level (RD), and open-door level (OP) are the explanatory variables of technical inefficiency. The statistical description of the variables selected in this paper is shown in Table 1.  [72,73], here we take RMB¥ as the unit of value in data processing.

Output and Input Variables
(1) Total output variable (unit: 100 million RMB¥). GDP is selected as total output of each province in YREB. The nominal GDP is converted into real GDP on the base of consumer price index of the year 2000. As total output variable, GDP is recorded as Y.
(2) Labor input variables (unit: 10,000 people). The year-end number of employees is selected as labor input of each province in YREB, and is recorded as L. The nominal average salary data of employees is converted into real average salary data based on consumer price index of the year 2000 as cost of labor input.
(3) Capital stock input variables (unit: 100 million RMB¥). The method of perpetual storage is used here that the formula is K it = I it /P it + K it−1 (1 − ρ), I it represents total fixed capital formation, P it is fixed asset formation price index, ρ is fixed assets depreciation rate (here is 10.96%) [74]. The formula for calculating the initial capital stock is K it = I it /(g it + ρ), among which, g it is mean value of the growth rate of the total fixed capital formation in the provinces and municipalities from 2000 to 2004. Through the above methods to obtain the real stock of capital in each province, is recorded as K.
(4) Energy input variables (unit: 10,000 tons of equivalent coal). In research of this paper, data on energy consumption and fossil energy equivalent coal coefficient from 2000 to 2016 was obtained from China Statistical Yearbooks and China Energy Statistical Yearbooks. Using annual energy consumption of each province as energy input, the consumption of coal, coke, crude oil, gasoline, kerosene, diesel, fuel oil, natural gas, and other energy is converted into a unified unit due to their different categories of energy consumption in each province. "Tonnes of equivalent coal" (see Table 2), the total energy consumption is recorded as E. The nominal unit price of equivalent coal, coming from the nationally announced coal-fired power generation grid benchmark price, is converted into real unit price of equivalent coal on the base of consumer price index of the year 2000. In 2016, the coal-fired power generation benchmark price is derived from the website of China National Development and Reform Commission. Table 2. Equivalent coal converted coefficient of fossil energy and carbon emissions factor [75]. (5) Carbon sequestration input variables (unit: 10,000 tons of carbon dioxide). Carbon sequestration mainly refers to the amount of carbon dioxide absorbed and stored by nature such as forest vegetation, agricultural crops, wetlands, and offshore water. The factors and indicators of carbon sequestration are shown in Table 3. Abide by the characteristics of natural resources in YREB, carbon sequestration of forests, agricultural crops, wetlands, and offshore water is measured as an indicator of carbon sequestration, which is recorded as CS. In this paper, its cost is recorded as zero. Table 3. Carbon sequestration conversion coefficient.

Impact Variables of Technical Inefficiency Terms
(1) Total carbon emissions in region (unit: 100 million tons). It is mainly calculated from consumption of coal, coke, crude oil, gasoline, kerosene, diesel, fuel oil, natural gas, and other energy sources in YREB. The carbon emission coefficient here (see Table 2) comes from the Institute of Environmental Development of the Chinese Academy of Agricultural Sciences, the Ecological Institute of the Chinese Academy of Forestry, and the Climate Impact Center of the China Institute of Environmental Sciences, which is originally based on the national standards of the People's Republic of China "General rules for calculation of comprehensive energy consumption"(GB/T2589-2008) and modified by the Climate Division of China National Development and Reform Commission "Guidelines for the preparation of provincial greenhouse gas inventories".
(2) Carbon emission intensity. Carbon emission intensity index refers to carbon dioxide emissions per unit of GDP, calculated on the basis of carbon emissions and the real GDP data of the provinces in YREB in the past years, recorded as CI.
(3) Ratio of secondary industry to GDP. It is directly derived from the Statistical Yearbooks of China and 11 provinces and municipalities in YREB, recorded as SI.
(4) Urbanization rate. Urbanization rate is a measure of the level of urbanization, generally using demographic indicators, that is, the proportion of urban population to the total. This paper uses the proportion of non-agricultural population to the total population at the year-end as the urbanization rate, here recorded as UE.
(5) R&D. R&D is represented by the ratio of the total R&D expenditure of each province to the regional GDP, which is recorded as RD.
(6) Opening-up to the outside world. It is expressed as the ratio of total imports and exports by province to regional GDP, recorded as OP.
With support of related variables and data analyzed above, empirical estimation of CSTFP growth in YREB can be accordingly carried out next.

Model Test
For the most suitable model estimation of frontier production function on CSTFP growth, the following original hypothesis (H 0 ) is established to test the model.
(1) Random frontier model production function form test: if H 0 is true, there is no interaction between variables, so it is possible to use the ordinary Cobb-Douglas production function without using the transcendental logarithmic form.
(2) Technical inefficiency feature information test: if H 0 is true, then µ obeys semi-normal distribution rather than truncated normal distribution, and the technical inefficiency term has no time-varying effect [71].
(3) Inspection of neutral technological progress in frontier production function: if H 0 is true, the form of technological progress is Hicks neutral technological progress that does not change the ratio of marginal product of capital to labor.
(4) Whether there is a factor of technological progress in frontier production function: The degree of freedom of constraint conditions is the number of constraint conditions, and the test results of each hypothesis are shown in Table 4.  According to the test results in Table 4, the four null hypotheses are rejected, which means that the ordinary Cobb-Douglas production function is not applicable here, there is a technical inefficiency item, technological progress exists, and it is not Hicks neutral technology progress. This shows it is reasonable to take translog production function and carry out random frontier analysis. The technical inefficiency term exists objectively [55][56][57].

Parameter Estimation
In this paper, we use Stata 14.0 to perform a random frontier regression for translog production function on CSTFP growth in YREB. For comparing the different impact of factors such as carbon sequestration, energy consumption, and carbon intensity on the estimation of CSTFP growth, regressions are carried out through different estimations: Estimation I, the two-factor estimation, representing the traditional TFP estimation with two input factors of labor and capital; Estimation II, the three-factor estimation, representing the green TFP estimation with energy consumption and carbon as well as labor and capital; Estimation III, the four-factor estimation, representing the CSTFP estimation with carbon sequestration added on the basis of three-factor estimation. The overall estimation results are shown in Table 5.  Note: (1) In the brackets in the table, the morning is the standard deviation, *, **, and *** represent the significance level of 10%, 5%, and 1% respectively; (2) γ is the technical efficiency term variance in the whole random The proportion of the error variance, the importance of the reaction inefficiency term relative to the whole random error; (3) When the model is established, the energy consumption and carbon sequestration are embedded in the frontier production function, and the carbon emission intensity is embedded in the technical inefficiency term.
As shown in Table 5, the value of Wald chi2 (13) is large and prob > chi2 = 0, indicating that the overall effect of the model is better. Most explanatory variables are significant at the level of 1% significant. γ value obtained by calculation of three estimations [γ = σ 2 µ /(σ 2 υ + σ 2 µ )] are close to 1, indicating the variance of the technical inefficiency term te for the overall random error term (υ − µ). The variance is the main one, the inefficiency term exists objectively, and the analysis method using the random frontier is reasonable.

The CSTFP Growth Is Obviously Improved Owing to Carbon Sequestration's Synergy Effect
As shown in Figure 4, the rate of TFP growth in YREB under the three estimations fluctuate with tendency of up-down-up, increasing yearly before 2005, then falling to 2009, and rebounding up after 2014. In total, the changes in the rate of TFP growth in YREB under the three estimations are very small, but the four-factor estimation with carbon sequestration (0.20%) are significantly better than three-factor estimation (−3.76%) and traditional two-factor estimation (−26.55%). Compared with the estimation results of three-factor estimation and two-factor estimation, the estimation results of four-factor estimation with carbon sequestration is further improved, suggesting that carbon sequestration produce two-way synergy effect on input and output for the growth of TFP in YREB. Therefore, the CSTFP growth in YREB is obviously improved. However, due to the alternation of the two driving forces of resource-driven economic growth dominated by energy consumption and innovation-driven economic growth improved by technological progress, the rate of CSTFP growth in YREB is deadlocked in the range of 0.05%-0.3%. As shown in Figure 4, the rate of TFP growth in YREB under the three estimations fluctuate with tendency of up-down-up, increasing yearly before 2005, then falling to 2009, and rebounding up after 2014. In total, the changes in the rate of TFP growth in YREB under the three estimations are very small, but the four-factor estimation with carbon sequestration (0.20%) are significantly better than three-factor estimation (-3.76%) and traditional two-factor estimation (-26.55%). Compared with the estimation results of three-factor estimation and two-factor estimation, the estimation results of fourfactor estimation with carbon sequestration is further improved, suggesting that carbon sequestration produce two-way synergy effect on input and output for the growth of TFP in YREB. Therefore, the CSTFP growth in YREB is obviously improved. However, due to the alternation of the two driving forces of resource-driven economic growth dominated by energy consumption and innovationdriven economic growth improved by technological progress, the rate of CSTFP growth in YREB is deadlocked in the range of 0.05%-0.3%.

Three of Four Components Drives Positively the CSTFP Growth
As shown in Table 6, compared to that by two-factor estimation and three-factor estimation, the driving factors estimated by four-factor estimation reveal better results in decomposition of CSTFP growth in YREB (average annual growth of 0.2%). Among them, the rate of technological progress reaches 1.72% with a contribution rate of 18.55% to CSTFP growth, the rate of technical efficiency change reaches 2.65% with a contribution rate of 28.59%, and the rate of change in scale efficiency reaches 0.37% with a contribution rate of 3.99%. Only the rate of factor allocation efficiency change is negative (-4.53%) with contribution rate of -48.87% to CSTFP growth in YREB. Note: a) CR represents contribution rate. b) the data derives from the annual average data from 2000 to 2016. Figure 5a shows the rate of change in technical efficiency in three estimations. In comparison, the rate of technical efficiency change measured by four-factor estimation fluctuates above zero (2.65%), and it is significantly two percentage points higher than the other two estimations in each

Three of Four Components Drives Positively the CSTFP Growth
As shown in Table 6, compared to that by two-factor estimation and three-factor estimation, the driving factors estimated by four-factor estimation reveal better results in decomposition of CSTFP growth in YREB (average annual growth of 0.2%). Among them, the rate of technological progress reaches 1.72% with a contribution rate of 18.55% to CSTFP growth, the rate of technical efficiency change reaches 2.65% with a contribution rate of 28.59%, and the rate of change in scale efficiency reaches 0.37% with a contribution rate of 3.99%. Only the rate of factor allocation efficiency change is negative (−4.53%) with contribution rate of −48.87% to CSTFP growth in YREB.  Figure 5a shows the rate of change in technical efficiency in three estimations. In comparison, the rate of technical efficiency change measured by four-factor estimation fluctuates above zero (2.65%), and it is significantly two percentage points higher than the other two estimations in each period. Technical efficiency growth is the primary contributor with the rate of 28.59% to CSTFP growth in YREB, which is supported by the four following aspects. Firstly, it is release of policy effects. The promotion of strategic measures such as "Ecological priority for green development" and "Elimination of smog for the blue sky" actively guide green development with clean production technology to increase real productivity and improve technical efficiency. Secondly, it is effective for introducing technology from the outside world (the inefficiency term OP is −0.138). By introduction, digesting and absorption of foreign technology for meeting domestic demand of development, the rate of technical efficiency change is improved. Finally, it is the result of urbanization for urban newcomers in use of technology (inefficient item UE is −0.975). The urban newcomers continue to improve their ability of employment by learning applicable technology and management experience, and this helps to accumulation of knowledge and application of technology innovation, and improvement of technological efficiency. But overall, this rate of technological efficiency change in the growth of CSTFP in YREB will be unsustainable if there is no steady increase in technological progress.

Technical Efficiency Growth Is the Primary Contributor to the CSTFP Growth
period. Technical efficiency growth is the primary contributor with the rate of 28.59% to CSTFP growth in YREB, which is supported by the four following aspects. Firstly, it is release of policy effects. The promotion of strategic measures such as "Ecological priority for green development" and "Elimination of smog for the blue sky" actively guide green development with clean production technology to increase real productivity and improve technical efficiency. Secondly, it is effective for introducing technology from the outside world (the inefficiency term OP is -0.138). By introduction, digesting and absorption of foreign technology for meeting domestic demand of development, the rate of technical efficiency change is improved. Finally, it is the result of urbanization for urban newcomers in use of technology (inefficient item UE is -0.975). The urban newcomers continue to improve their ability of employment by learning applicable technology and management experience, and this helps to accumulation of knowledge and application of technology innovation, and improvement of technological efficiency. But overall, this rate of technological efficiency change in the growth of CSTFP in YREB will be unsustainable if there is no steady increase in technological progress.
(a) Technical efficiency change.
(b) Scale efficiency change.

Scale Efficiency is Stabilized While Factor Allocation Efficiency is Lowered
Carbon sequestration is conducive to promoting the stable growth of scale efficiency of CSTFP growth in YREB. Figure 5b shows the different trends of the scale efficiency changes in the three estimations. The rate of scale efficiency change calculated by four-factor estimation fluctuates within a narrow range of 0%-0.4%, which is relatively stable. The change rate of scale efficiency calculated by three-factor estimation has upward tendency after 2006, indicating that the change of scale efficiency brought by energy conservation, CO2 emission reduction and low-carbon development, is significant. However, the change rate of return to scale calculated by two-factor estimation has been decreasing all the time and becomes negative after 2012, indicating that ignoring energy consumption and CO2 emissions is not conducive to the improvement of scale efficiency. In comparison, the fourfactor estimation with introducing carbon sequestration and energy consumption integrates the estimation results of the other two estimations, showing that it stabilizes the steady growth of scale efficiency of CSTFP in the YREB. However, as Figure 5c shows, among the tendency of the change rate in factor allocation efficiency in the three estimations, the results of three-factor estimation with energy consumption and CO2 emissions (-4.75%) and four-factor estimation with carbon sequestration added (-4.53%) are significantly lower than that of traditional two-factor estimation (-1.36%), indicating that the introduction of energy consumption, CO2 emissions, and carbon sequestration lowers the efficiency of factor allocation of CSTFP in YREB.

Scale Efficiency is Stabilized While Factor Allocation Efficiency is Lowered
Carbon sequestration is conducive to promoting the stable growth of scale efficiency of CSTFP growth in YREB. Figure 5b shows the different trends of the scale efficiency changes in the three estimations. The rate of scale efficiency change calculated by four-factor estimation fluctuates within a narrow range of 0%-0.4%, which is relatively stable. The change rate of scale efficiency calculated by three-factor estimation has upward tendency after 2006, indicating that the change of scale efficiency brought by energy conservation, CO 2 emission reduction and low-carbon development, is significant. However, the change rate of return to scale calculated by two-factor estimation has been decreasing all the time and becomes negative after 2012, indicating that ignoring energy consumption and CO 2 emissions is not conducive to the improvement of scale efficiency. In comparison, the four-factor estimation with introducing carbon sequestration and energy consumption integrates the estimation results of the other two estimations, showing that it stabilizes the steady growth of scale efficiency of CSTFP in the YREB. However, as Figure 5c shows, among the tendency of the change rate in factor allocation efficiency in the three estimations, the results of three-factor estimation with energy consumption and CO 2 emissions (−4.75%) and four-factor estimation with carbon sequestration added (−4.53%) are significantly lower than that of traditional two-factor estimation (−1.36%), indicating that the introduction of energy consumption, CO 2 emissions, and carbon sequestration lowers the efficiency of factor allocation of CSTFP in YREB.

Distance on CSTFP Growth in Three Watershed Segments is Significant
Carbon sequestration stimulate more the growth of CSTFP in the upper reaches of YREB. As Figure 6 shows, the average growth rate of CSTFP in the upper reaches of YREB is reaching 0.62%, which is the highest, three times higher than that in YREB (0.2%), while it is 0.11% in the lower reaches and -0.14% in the middle reaches. It indicates that the great capacity of carbon sequestration provided by the abundant ecological resource pool originating from the Qinghai-Tibet plateau and Yunnan-Guizhou plateau has a significant impact on the growth of CSTFP in the upper reaches of YREB, while lack of natural resource conditions in the middle reaches and the lower reaches of YREB, there is no such influence of carbon sequestration on the growth of CSTFP in those regions. Carbon sequestration stimulate more the growth of CSTFP in the upper reaches of YREB. As Figure 6 shows, the average growth rate of CSTFP in the upper reaches of YREB is reaching 0.62%, which is the highest, three times higher than that in YREB (0.2%), while it is 0.11% in the lower reaches and -0.14% in the middle reaches. It indicates that the great capacity of carbon sequestration provided by the abundant ecological resource pool originating from the Qinghai-Tibet plateau and Yunnan-Guizhou plateau has a significant impact on the growth of CSTFP in the upper reaches of YREB, while lack of natural resource conditions in the middle reaches and the lower reaches of YREB, there is no such influence of carbon sequestration on the growth of CSTFP in those regions. Contributors to the growth rate of CSTFP in each watershed segment of YREB are significantly different as Table 7 shows.
(a) In the upper reaches of YREB, the growth rate of CSTFP reaches 0.62%, three times higher than that of YREB (0.20%). The major contributors are technological progress with growth rate of 3.86% and contribution rate of 23.84%, and technical efficiency with growth rate of 4.06% and the contribution rate of 25.07%. While the weakness is factor allocation efficiency with growth rate of -7.80% and contribution rate of -48.09%.
(b) In the middle reaches of YREB, the growth rate of CSTFP is -0.14%, 2.43 times lower than that of YREB (0.20%). Its first positive contributor is technical progress with growth rate of 2.38% and contributing rate of 20.96%. Next, is technical efficiency with growth rate of 3.0% and contributing rate of 26.35%, while the foot-dragging driver is factor allocation efficiency with the growth rate of 5.75% and contribution rate of -50.62%.
(c) In the lower reaches of YREB (the Yangtze River Delta), the growth rate of CSTFP reaches 0.11%, which is 0.55 times lower than that of YREB (0.20%). The main contributor is factor allocation efficiency with growth rate of 1.44% and the contribution rate of 34.9%. While the most weakness is technological progress with growth rate of -2.01% and contribution rate of -48.66%.  Contributors to the growth rate of CSTFP in each watershed segment of YREB are significantly different as Table 7 shows.  (a) In the upper reaches of YREB, the growth rate of CSTFP reaches 0.62%, three times higher than that of YREB (0.20%). The major contributors are technological progress with growth rate of 3.86% and contribution rate of 23.84%, and technical efficiency with growth rate of 4.06% and the contribution rate of 25.07%. While the weakness is factor allocation efficiency with growth rate of -7.80% and contribution rate of −48.09%.
(b) In the middle reaches of YREB, the growth rate of CSTFP is −0.14%, 2.43 times lower than that of YREB (0.20%). Its first positive contributor is technical progress with growth rate of 2.38% and contributing rate of 20.96%. Next, is technical efficiency with growth rate of 3.0% and contributing rate of 26.35%, while the foot-dragging driver is factor allocation efficiency with the growth rate of 5.75% and contribution rate of −50.62%.
(c) In the lower reaches of YREB (the Yangtze River Delta), the growth rate of CSTFP reaches 0.11%, which is 0.55 times lower than that of YREB (0.20%). The main contributor is factor allocation efficiency with growth rate of 1.44% and the contribution rate of 34.9%. While the most weakness is technological progress with growth rate of −2.01% and contribution rate of −48.66%.
The growth rate of technological progress driver is relatively the highest in the upper reaches compared to those in the middle reaches and the lower reaches, as Figure 7a shows. It indicates that the frontier production function in the lower reaches has been moving downward on the bottom due to its slow technological progress keeping input factors and production structure unchanged since 2003, while the leading production function stays on the top and the rate of technological progress changes better in the middle reaches and the upper reaches year by year.
The growth of technical efficiency driver in the lower reaches leads the change of technical efficiency drivers in the middle reaches and the upper reaches. Figure 7b shows that the highest change rate of technical efficiency is in the lower reaches, while the lowest rate is in the upper reaches. The change rate in the middle reaches and upper reaches approaches convergent to that of the lower reaches, indicating that the growth of technical efficiency of the middle reaches and the upper reaches is led by that of the lower reaches. Note: a) CE represents contribution rate. b) The data is from 11 provinces and municipalities from 2000 to 2016, taking the average time of each region.
The growth rate of technological progress driver is relatively the highest in the upper reaches compared to those in the middle reaches and the lower reaches, as Figure 7a shows. It indicates that the frontier production function in the lower reaches has been moving downward on the bottom due to its slow technological progress keeping input factors and production structure unchanged since 2003, while the leading production function stays on the top and the rate of technological progress changes better in the middle reaches and the upper reaches year by year.   The growth of technical efficiency driver in the lower reaches leads the change of technical efficiency drivers in the middle reaches and the upper reaches. Figure 7b shows that the highest change rate of technical efficiency is in the lower reaches, while the lowest rate is in the upper reaches. The change rate in the middle reaches and upper reaches approaches convergent to that of the lower reaches, indicating that the growth of technical efficiency of the middle reaches and the upper reaches is led by that of the lower reaches.
Difference in the rate of scale efficiency drivers in three watersheds of YREB was reversed by year of 2009. As Figure 7c shows, the rate of change in scale efficiency in the upper and middle reaches fluctuated around zero from 2001 to 2016. Compared to that in the upper and middle reaches, the trend of change in scale efficiency in the lower reaches showed a higher tendency before 2007 but dropped down and lower than that in the middle reaches and the lower reaches after 2009, indicating that the simple expansion of scale has not promoted further economic growth for the lower reaches of YREB. Figure 7d shows contribution of factor allocation driver to CSTFP growth in three segments of YREB. Overall, the change rate of factor allocation efficiency in each watershed segment shows a volatility. Among them, after 2003, the change rate of factor allocation efficiency in the lower reaches is above zero, indicating that the factor allocation is relatively reasonable. While the change rate of factor allocation efficiency in the middle and upper reaches is below zero all the time, reflecting that Difference in the rate of scale efficiency drivers in three watersheds of YREB was reversed by year of 2009. As Figure 7c shows, the rate of change in scale efficiency in the upper and middle reaches fluctuated around zero from 2001 to 2016. Compared to that in the upper and middle reaches, the trend of change in scale efficiency in the lower reaches showed a higher tendency before 2007 but dropped down and lower than that in the middle reaches and the lower reaches after 2009, indicating that the simple expansion of scale has not promoted further economic growth for the lower reaches of YREB. Figure 7d shows contribution of factor allocation driver to CSTFP growth in three segments of YREB. Overall, the change rate of factor allocation efficiency in each watershed segment shows a volatility. Among them, after 2003, the change rate of factor allocation efficiency in the lower reaches is above zero, indicating that the factor allocation is relatively reasonable. While the change rate of factor allocation efficiency in the middle and upper reaches is below zero all the time, reflecting that the factor allocation is not reasonable enough.

Technological Progress Is Improved by Synergistic Effect of Energy Inputs and Carbon Sequestration
As Figure 5d shows, the rate of technological progress measured by three-factor estimation introducing energy consumption and four-factor estimation introducing carbon sequestration is converging and almost overlapping together in each period, and it is 25% higher than that by the results of the two-factor estimation, indicating that energy consumption has a positive correlation with carbon sequestration (co-efficient is 0.335) and this stimulates the improvement of technological progress. Because, energy consumption and industrial restructuring under the constraints of energy conservation and emission reduction (inefficiency item SI decreased from −0.923 to −1.252) stimulate related technology advance and reduce CO 2 emissions (inefficiency term CI has been reduced from 0.56 to 0.524). Meanwhile, carbon sequestration supplies clean and renewable energy to replace the fossil energy and makes energy consumption structure transited from CO 2 emissions-intensive energy to more clean and renewable energy. However, on the whole, the rate of technological progress in YREB has been declining year by year or even falling into a negative growth trend because of insufficient investment in research and development (inefficiency term R&D is stabilized from 0.109 to 0.124), and this will weaken sustainable ability of economic development in YREB.

Technical Efficiency Growth Is Benefited from Harmonious Symbiosis Influence of Carbon Sequestration
Carbon sequestration promotes harmonious symbiosis correlations with economic development and environment optimization, and this is beneficial for technical efficiency growth of CSTFP. On one end, carbon sequestration stimulates forest vegetation and ecological restoration of mountains, rivers, grasslands, wetlands, and other ecological activities, CO 2 is absorbed through photosynthesis of nature to release oxygen, and this makes ecological productivity increased and environmental improvement. On the other end, carbon sequestration creates natural resources and conditions such as carbon-rich soil, clean air, and water for agriculture, afforestation, animal husbandry, and fishery production. Under the input of labor, capital, and technology use, diversified products such as timber, food, raw materials are produced as output and create economic value for society. The function of carbon sequestration for resulting in the superposition of multitudinous input elements promotes technical efficiency and integration of natural reproduction and economic reproduction.

Inefficient Resource Allocation Is Mainly Due to Conflicts of Interest and Insufficient Coordination among Regions in YREB
In the CSTFP growth of YREB, although the concept of "ecological priority for green development" has been deeply rooted in the hearts of the people, and factor allocation efficiency has presented upward trend in total, but local governments encounter conflicts of interest in aspects of industrial development transformation, adjustment of energy structure, responsibility sharing of CO 2 emission reduction, ecological investment, and environmental governance, and each region has its own administrative boundary/territory, acts independently, lack of compromise of interests, and mutual cooperation. This results in negative growth of resource allocation efficiency (average output elasticity of capital is −1.66). Meanwhile, the unbalanced benefits being allocated between the contributors and beneficiaries of afforestation, improved forestry, or agricultural practices and ecological management (average output elasticity of CS is −1.491) deviate or restrict allocation efficiency of relevant factors improved better than before.

CSTFP Growth Is Inadequate and Uncoordinated among 11 Provinces and Municipalities of YREB
As Figure 8 shows, the upper reaches and its provinces and municipalities of Guizhou, Chongqing, and Sichuan are the leaders of the CSTFP growth of the whole basin, but inadequacies and incongruities of CSTFP growth are embodied in the middle reaches and lower reaches and relative provinces and municipalities of YREB. For example, the rate of CSTFP growth in Guizhou, Chongqing, and Sichuan reaches 3.22%, 0.87%, and 0.68%, respectively, which is higher than that in Shanghai (0.52%) and Jiangsu (0.17%). Meanwhile, the rate of CSTFP growth is negative in Yunnan (−2.3%), Hubei (−0.12%), Jiangxi (−1.09%), and Zhejiang (−0.36%), respectively. The gap of CSTFP growth rate between the highest and lowest is 5.52 percentages among different provinces and municipalities of YREB, which is driven relatively by differential growth levels of technical progress, technical efficiency, scale efficiency, and factor allocation efficiency. Shanghai (0.52%) and Jiangsu (0.17%). Meanwhile, the rate of CSTFP growth is negative in Yunnan (-2.3%), Hubei (-0.12%), Jiangxi (-1.09%), and Zhejiang (-0.36%), respectively. The gap of CSTFP growth rate between the highest and lowest is 5.52 percentages among different provinces and municipalities of YREB, which is driven relatively by differential growth levels of technical progress, technical efficiency, scale efficiency, and factor allocation efficiency.  As Figure 9a shows, the higher rate of change in technical progress in the upper and middle reaches and their provinces and municipalities support the technical progress of the whole basin in YREB. Guizhou (5.85%) and Yunnan (5.17%) in the upper reaches have a higher rate of technological progress than Shanghai (−2.82%) and Jiangsu (−2.73%) in the lower reaches, where there is a gap of 8.64 percent to close because of the coefficient of intersection with time β et and β cst , which is 0.0187 and 0.0269, respectively, indicating that, with increase of carbon sequestration and reduction of energy consumption, provincial level of relevant technical progress has disintegrated, multi-speed green and low-carbon technical progress growth has produced based on differential natural resources and environmental advantage in provinces and municipalities of YREB. Shanghai (0.52%) and Jiangsu (0.17%). Meanwhile, the rate of CSTFP growth is negative in Yunnan (-2.3%), Hubei (-0.12%), Jiangxi (-1.09%), and Zhejiang (-0.36%), respectively. The gap of CSTFP growth rate between the highest and lowest is 5.52 percentages among different provinces and municipalities of YREB, which is driven relatively by differential growth levels of technical progress, technical efficiency, scale efficiency, and factor allocation efficiency. As Figure 9a shows, the higher rate of change in technical progress in the upper and middle reaches and their provinces and municipalities support the technical progress of the whole basin in YREB. Guizhou (5.85%) and Yunnan (5.17%) in the upper reaches have a higher rate of technological progress than Shanghai (-2.82%) and Jiangsu (-2.73%) in the lower reaches, where there is a gap of 8.64 percent to close because of the coefficient of intersection with time et β and cst β , which is 0.0187 and 0.0269, respectively, indicating that, with increase of carbon sequestration and reduction of energy consumption, provincial level of relevant technical progress has disintegrated, multi-speed green and low-carbon technical progress growth has produced based on differential natural resources and environmental advantage in provinces and municipalities of YREB.
(a) Technical progress growth rate. (b) Technical efficiency growth rate.  . CSTFP growth drivers among 11 provinces in YREB (label is the same as Figure 8).
As shown in Figure 9b, the higher rate of change in technical efficiency in the upper and middle reaches and their provinces and municipalities leads the technical efficiency growth of the whole basin in YREB. The growth rate of technical efficiency in Shanghai (0.52%), Jiangsu (0.47%), and Zhejiang (-0.11%) in the lower reaches is significantly lower than that of the provinces and municipalities in the upper reaches and middle reaches. As Table 7 shows, based on the four-factor estimation, technical efficiency growth contributes 25.07% and 26.35%, respectively, and becomes the main source of relative CSTFP growth in the upper reaches and middle reaches of YREB, owing to the speeding up of urbanization (inefficient UE coefficient is -0.795), the second industry development (inefficient SI coefficient is -1.225), and the improvement of open outside to the world (inefficient OP coefficient is -0.138) jointly promote the technical efficiency of the upper reaches and middle reaches of YREB.
Zhejiang and Jiangsu are leaders at the growth of scale efficiency among the provinces and municipalities of YREB. As Figure 9c shows, except for Guizhou (-0.12%), Anhui (-0.47%), and Shanghai (-1.29%), the rate of change in scale efficiency in most of provinces and municipalities of YREB is positive. For example, the rate of change in scale efficiency is relatively higher in Jiangsu (1.01%) and Zhejiang (1.45%) in the lower reaches with positive contribution of 20.54% and 44.44%, respectively, to the growth rate of CSTFP, but the rate of change in scale efficiency in Shanghai is the lowest (-1.29%) with negative contribution (-14.78%) to the growth rate of CSTFP. The gap between the highest and the lowest is 2.74 percentages. Due to constraints of resources and environment, economic growth brought about by blind expansion of economic scale has been unsustainable, especially in Shanghai, Guizhou, and Anhui. Therefore, it is urgent to eliminate surplus production capacity and nurture new kinetic energy for transforming regional CSTFP growth into a high-quality mode.
Negative growth is a great challenge in the change rate of factor allocation efficiency in the provinces and municipalities of YREB. As shown in Figure 9d, only Shanghai (4.12%) and Jiangsu (1.07%) have positive change rate of factor allocation efficiency, contributing 47.15% and 21.66% to the growth rate of CSTFP, respectively. The change rate of factor allocation efficiency in other provinces and municipalities is negative, which poses a great challenge to the growth of CSTFP and high-quality development in YREB on the whole.

Total Synergies Are Not Being Unleashed Enough for Promotion of the CSTFP Growth in YREB
As Table 5 shows, labor input has a positive impact (3.278) on economic growth while capital investment has a negative impact (-0.571) in two-factor estimation with ignorance of impact of energy consumption and CO2 emissions on TFP growth. Introducing energy consumption into endogenous variables in three-factor estimation, the impact of energy and labor on TFP growth is positive with respectively average output elasticity of 2.366 and 1.863, while that of capital on TFP growth is Figure 9. CSTFP growth drivers among 11 provinces in YREB (label is the same as Figure 8). Figure 9b, the higher rate of change in technical efficiency in the upper and middle reaches and their provinces and municipalities leads the technical efficiency growth of the whole basin in YREB. The growth rate of technical efficiency in Shanghai (0.52%), Jiangsu (0.47%), and Zhejiang (−0.11%) in the lower reaches is significantly lower than that of the provinces and municipalities in the upper reaches and middle reaches. As Table 7 shows, based on the four-factor estimation, technical efficiency growth contributes 25.07% and 26.35%, respectively, and becomes the main source of relative CSTFP growth in the upper reaches and middle reaches of YREB, owing to the speeding up of urbanization (inefficient UE coefficient is −0.795), the second industry development (inefficient SI coefficient is −1.225), and the improvement of open outside to the world (inefficient OP coefficient is -0.138) jointly promote the technical efficiency of the upper reaches and middle reaches of YREB.

As shown in
Zhejiang and Jiangsu are leaders at the growth of scale efficiency among the provinces and municipalities of YREB. As Figure 9c shows, except for Guizhou (−0.12%), Anhui (−0.47%), and Shanghai (−1.29%), the rate of change in scale efficiency in most of provinces and municipalities of YREB is positive. For example, the rate of change in scale efficiency is relatively higher in Jiangsu (1.01%) and Zhejiang (1.45%) in the lower reaches with positive contribution of 20.54% and 44.44%, respectively, to the growth rate of CSTFP, but the rate of change in scale efficiency in Shanghai is the lowest (−1.29%) with negative contribution (−14.78%) to the growth rate of CSTFP. The gap between the highest and the lowest is 2.74 percentages. Due to constraints of resources and environment, economic growth brought about by blind expansion of economic scale has been unsustainable, especially in Shanghai, Guizhou, and Anhui. Therefore, it is urgent to eliminate surplus production capacity and nurture new kinetic energy for transforming regional CSTFP growth into a high-quality mode.
Negative growth is a great challenge in the change rate of factor allocation efficiency in the provinces and municipalities of YREB. As shown in Figure 9d, only Shanghai (4.12%) and Jiangsu (1.07%) have positive change rate of factor allocation efficiency, contributing 47.15% and 21.66% to the growth rate of CSTFP, respectively. The change rate of factor allocation efficiency in other provinces and municipalities is negative, which poses a great challenge to the growth of CSTFP and high-quality development in YREB on the whole.

Total Synergies Are Not Being Unleashed Enough for Promotion of the CSTFP Growth in YREB
As Table 5 shows, labor input has a positive impact (3.278) on economic growth while capital investment has a negative impact (−0.571) in two-factor estimation with ignorance of impact of energy consumption and CO 2 emissions on TFP growth. Introducing energy consumption into endogenous variables in three-factor estimation, the impact of energy and labor on TFP growth is positive with respectively average output elasticity of 2.366 and 1.863, while that of capital on TFP growth is negative (−1.368). And further as shown in four-factor estimation with addition of carbon sequestration to endogenous variables, the average output elasticity of energy consumption (0.581) and labor (3.42) is positive, but that of carbon sequestration (−1.491) and capital(−1.66) is negative, indicating that synergies of carbon sequestration on CSTFP growth are not being unleashed enough due to the deadlock intertwined and knotted by the rival game of resources-driven economic growth dominated by fossil energy consumption and environmental-friendly economic growth driven by technological progress.

Conclusions
For making it clear whether carbon sequestration is both effective at mitigating climate change and promoting economic growth, in this paper, based on carbon sequestration's feature of harmonious symbiosis correlations between economic development and environment optimization, we firstly design an conceptual framework of CSTFP growth by taking carbon sequestration into endogenous variables from the perspective of input-output mixes, and then set up an estimation model of CSTFP growth by adopting the stochastic frontier panel model with translog production function. Finally, we estimate and decompose CSTFP growth in YREB with 17-year data of 9 provinces and 2 municipalities under three estimations and in three watershed segments. The results show: Compared to the result of two-factor estimation on traditional TFP growth in YREB, the CSTFP growth in YREB by four-factor estimation is obviously improved by 26.74 percentages (from −26.55% to 0.20%), which owes to carbon sequestration's synergy effect on input and output mixes as an endogenous variable in the CSTFP growth.
Three of four components positively drive the CSTFP growth in YREB, of which technical efficiency change is the primary driver with contribution share of 28.59%, then technological progress change with contribution share of 18.55%, and scale efficiency change with contribution share of 3.99%, while factor allocation efficiency is the negative driver with reverse contribution share of −48.87% in the CSTFP growth.
The difference on the CSTFP growth in three watershed segments of YREB is significant, of which the highest is that of the upper reaches (0.62%) efficiently driven by technical efficiency change (25.07%) and technological progress change (23.84%), then that of the lower reaches (0.11%) efficiently driven by factor allocation efficiency (34.9%), and scale efficiency change (9.43%), and finally that of the middle reaches (−0.14%) efficiently driven by technical efficiency change (26.35%) and technological progress change (20.96%).
From the above analysis, technological progress is improved by synergistic effect of energy inputs and carbon sequestration, technical efficiency growth is benefited from harmonious symbiosis influence of carbon sequestration, meanwhile, inefficient resource allocation is mainly due to conflicts of interest and insufficient coordination among regions in YREB and overall synergies not being unleashed enough for driving CSTFP growth in YREB.

Overall Synergies of Carbon Sequestration on Input-Output Mixes Should Be Given Full Play to Driving CSTFP Growth in YREB
The human activity circle is parasitic in the natural ecosystem, and both are coupled into a living community of carbon sequestration. Therefore, the following strategies should be taken to stimulate overall synergies of carbon sequestration on input-output mixes. Firstly, afforestation is the direct strategy for the increase of carbon sequestration's synergy effect on active use of photosynthesis of nature, including barren hills and slopes greening, forests and grasslands protection, forestry and animal husbandry production, and development of ecological agriculture and farmland protection. Secondly, river pollution control and ecological restoration is the major and urgent strategy for self-restraint of carbon sequestration's synergy effect on water resources protection, such as pollution prevention and control and ecological restoration activities of rivers, lakes, reservoirs, wetlands, and marine protection. Thirdly, fossil fuels energy low-carbon transformation is the primary strategy for the improvement of carbon sequestration's synergy effect on CO 2 emissions reduction and energy economic development by replacing coal power with hydro-power, wind power, and nuclear power for energy restructuring. Finally, low-carbon technologies R&D and industrialization is the fundamental strategy to improve carbon sequestration's synergy effect on technical progress and efficiency by active promoting CCUS technology R&D and industrialization, clean production and green manufacturing, harmless disposal of garbage, resource recycling, and circular economic development.

Regional Close Cooperation and Harmonious Symbiosis Development Should Be Strengthened among Provinces and Municipalities of YREB
The upper reaches of YREB are rich in coal and hydro-power resources, with high forest coverage and outstanding natural carbon sink. With Chongqing and Chengdu as the strategic hubs, regional and inter-provincial cooperation should be improved for afforestation, hydro-power and wind power development, and coal cleaning. The cooperation mechanism needs to be continuously deepened jointly on development of the western region, development of the YREB, and construction of the "Belt and Road" economy and environmental conservation, for actively promoting the optimal allocation of natural ecological and economic diversified resources, and striving to improve resources allocation efficiency for CSTFP growth in YREB.
The middle reaches of YREB distribute natural rivers networks, clusters of manufacturing industry, prosperous agriculture, developed transportation, and energy industrialization. Taking the core cities of Wuhan and Nanchang as bases, focus on central rising strategy, the cross-provincial consultation, and cooperation should be improved on aspects of river pollution control and ecological restoration, fossil fuels energy transition, and low-carbon agriculture and low-carbon technologies R&D and industrialization, for continuously guiding local departments and universities, institutions, enterprises, and industry alliances actively to participate in cross-regional inter-provincial exchanges and cooperation for information sharing, joint innovation, industrial collaboration, and standard docking, and enhance resource allocation efficiency and economies of scale for CSTFP growth.
The lower reaches of YREB locate at the mouth of the sea, blessed with developed information and a strong capacity to gather resources. Setting up Shanghai as the strategic leader, the lower reaches should actively use the national strategy of the Yangtze River Delta rise to improve regional coordination and cooperation mechanism of YREB, strengthen innovation-leading, and increase local original technology while maintaining the leading edge of foreign technology absorption and digestion. R&D and industrialization should be focused on low-carbon technologies and clean energy and environmental technology for increasing the growth rate of technological progress and technological efficiency and radiating throughout the whole of YREB, and comprehensively enhancing the growth of CSTFP and modern development capacity of YREB.