Article Evaluation of Forestry Ecological Efficiency: A Spatiotemporal Empirical Study Based on China’s Provinces

Shuai Chen and Shunbo Yao *

Center of Economy and Environmental Management Research, College of Economics and Management, Northwest A & F University, Xianyang 712100, China; [email protected] * Correspondence: [email protected]; Tel.: +86-029-8708-1291

Abstract: Forests play a very important role in carbon dioxide emissions and climate change, and the development of China’s forestry is of great significance to our citizens. However, it is an arduous task for us to improve forestry output at a high and good level while taking environmental factors into account. In this paper, the non-expected super-efficiency SBM (slacks-based measure) model was used to measure the forestry ecological efficiency (FEE) of 31 provinces in China from 2004 to 2018, and the spatial and temporal evolution of FEE in different regions of China was analysed by using spatial econometric method. Tobit regression and random forest algorithm were selected to analyze the influence on FEE. The results showed that, firstly, the average annual increase of the national total factor change of China’s forestry was 1.2%, and that the average annual increase of the national total factor productivity change in the eastern region was lower than that in the central and western regions. Secondly, the distribution of China’s FEE of the northeast and the south was higher, and FEE of China’s middle regions was relatively lower in 2004, but then the FEE in Northeast China has decreased, and the FEE has increased gradually from north to south in 2018. The agglomeration of high-tech industries in most regions of China had obvious positive spatial correlation characteristics in 2018. Thirdly, there was a negative correlation between forestry fixed  assets investment and FEE, environmental regulation was an important factor affecting the ecological  efficiency of forestry in China, and the level of economic development and industrial structure also Citation: Chen, S.; Yao, S. Evaluation had a certain impact on FEE. of Forestry Ecological Efficiency: A Spatiotemporal Empirical Study Keywords: forestry ecological efficiency; spatial distribution; super-efficiency SBM; spatiotemporal Based on China’s Provinces. Forests difference; influencing factors 2021, 12, 142. https://doi.org/ 10.3390/f12020142

Received: 12 December 2020 1. Introduction Accepted: 21 January 2021 Published: 26 January 2021 Forestry is not only an important foundation of China’s national economy, but also which is an important part of public welfare undertakings. The ownership of the whole

Publisher’s Note: MDPI stays neutral people and collective is the main form in terms of ownership of forests in china, ownership with regard to jurisdictional claims in by the whole people refers to state-owned forestry accounting for more than 70% of the published maps and institutional affil- total forest area in China, including state-owned forestry enterprises and state-owned iations. forest farms. China’s forestry industrial structure and output value have been constantly optimized and improved in recent years. Specifically, China’s forest coverage rate has increased from 18.2% in 2004 to 23% in 2018, and the total forestry output value has increased from 93.65 billion yuan in 2000 to 575.57 billion yuan in 2019. The growth rate of forestry was more than five times, and it has a continuous upward trend. Although Copyright: © 2021 by the authors. Licensee MDPI, Basel, Switzerland. the output value of forestry has increased rapidly, the quality and cost of growth were This article is an open access article easily ignored by this simple quantity of growth, and people have had to face the ecological distributed under the terms and imbalance, crisis, environmental pollution, and other negative problems. Since conditions of the Creative Commons the 19th National Congress of the Communist Party of China, the State Council and the Attribution (CC BY) license (https:// people’s government have focused on forestry development sustainably. creativecommons.org/licenses/by/ FEE is to measure the quantity and quality of forestry products under the premise 4.0/). of minimizing environmental pollution and resource consumption [1]. At present, more

Forests 2021, 12, 142. https://doi.org/10.3390/f12020142 https://www.mdpi.com/journal/forests Forests 2021, 12, 142 2 of 17

and more scholars have turned to forestry efficiency. Some scholars used data envelop analysis (DEA) model to calculate the efficiency value generally. Zhang et al. [2] used two-stage DEA to analyze the data of 30 provinces in China and found that the ecological efficiency gradually decreased from the east to the west, and the eastern coastal areas had high economic efficiency but low environmental efficiency. Wang [3] used the model to calculate the efficiency of China’s green economy and found that the green efficiency value of developed regions, such as Shanghai, Beijing, and Tianjin, was at a high level, and the relationship between environmental regulation and green economic efficiency was inverted “U”. Lin and Ge [4] adopted the three-stage DEA model to adjust regional forest carbon by introducing economic and environmental factors. Chen et al. [5] studied the effective number of DEA efficiency provinces increased from 2003 to 2017, however, which still was lower than half of all provinces. The efficiency value of most provinces increased evidently in the southern forest area while the efficiency value of the northern forest area decreased. In addition, some scholars used spatial econometric method to analyze the efficiency problem. Chen et al. [6] analyzed the impact of technological innovation on urban ecological efficiency by using the spatial Durbin model. They found that 273 prefecture level cities in China showed strong spatial heterogeneity, and higher innovation ability can greatly improve urban ecological efficiency. Ren et al. [7] adopted the s-ebm mixed distance model of bad output to find that the positive spatial autocorrelation of ecological efficiency in China has gradually increased, and gradually decreased from relatively developed areas to developing areas. Zheng et al. [8] concluded the externalities brought by industrial agglomeration would promote the ecological efficiency of forestry industry in econometric model regression. There were also some scholars evaluated the forestry industry. Xiong [9] used the stochastic frontier analysis method to find that the forestry production efficiency of the six northwest provinces was in a state of decline from 2005 to 2015, and the spatial differences of the forestry production efficiency of the provinces were gradually reduced, Moreover, the collective forest right reform had a negative impact to forest right reform. Yang et al. [10] calculated the total factor productivity of state-owned key forestry enterprises in northeast, southwest and northwest regions by using Malmquist index. They found that the average growth rate of total factor productivity in Northwest China was higher than that in northeast and southwest of China. China’s western development has made more and more contributions to the economic development of the western region. Foreign scholars have also made good achievements in forestry research. One of the most frequently used models to measure forestry efficiency was DEA through the study of foreign literature. Viitala and Hanninen [11] studied the efficiency of regional forestry committees funded by 19 countries. They found that there were great differences in the efficiency of the Forestry Commission, and the investment saving potential was about 20%. Sporocic et al. [12] evaluated the efficiency of the basic organizational unit of forest office in Croatia by using DEA, they studied that DEA was a powerful multi criteria decision-making tool, and it played a significant part in forest management assessment. Similarly, Sowlati [13] proved that DEA had a good applicability in forestry efficiency measurement in his study. Toma et al. [14] adopted a nonparametric method to evaluate the efficiency performance of Italian forestry companies, which was stratified according to the environmental risk measurement. They found that the input and output oriented efficiency performance is higher for forestry companies belonging to medium and high environmental risk categories in Italy. However, forestry companies facing low environmental risks have shown greater progress in improving efficiency. In this context, this paper attempts to measure the total factor productivity of forestry from 2004 to 2018 to study the dynamic changes of forestry production. The super-efficiency SBM model is used to measure the ecological efficiency of forestry in three periods in recent years to reveal the distribution characteristics and evolution law of FEE in different periods in China. In addition, the internal relationship between various factors and FEE was Forests 2021, 12, 142 3 of 17

analyzed in the study. These results provide a reference framework for macro policy making to coordinate the development of forestry.

2. Materials and Methods 2.1. Indicator Description and Data Sources The input elements of FEE should at least include land, labor and capital [15], this pa- per calculates the ecological efficiency of forestry from the perspective of green production. The aim is to achieve the maximum forestry output at the cost of less resource reduction and environmental pollution. Considering the characteristics of forestry industry and the availability of data, the area of forestry land is used as the index of land production factors, the number of employees in the forestry system at the end of each year is taken as the labor factor, and capital element is the annual fixed asset investment. These elements are the three basic input elements of forestry industry investment. In addition, energy investment is another input element and it’s index is obtained by calculating the proportion of forestry output value to industrial output value and industrial energy investment. These indicators can be seen from Table1.

Table 1. Input and output indicators for measuring forestry ecological efficiency.

The Index Type Level Indicators The Secondary Indicators Unit Input Land The forest area Ten thousand km2 The number of employees in the forestry Labor People system at the end of each year Capital Fixed asset investment Ten thousand Yuan Energy Forest energy investment Hundred million Yuan Desirable output Economic benefit Total output value of forestry Hundred million Yuan Direct benefit Timber felling Ten thousand cubic meters Undesirable output Wastewater pollution Forestry wastewater discharge Ten thousand cubic meters Air pollution Waste gas emission Ten thousand ton

The expected output includes total output value of forestry and timber felling, and the non-expected output mainly include forestry wastewater discharge and waste gas emission [1]. The data are calculated from the ratio of forestry output value to total output value and the discharge amount of waste gas and wastewater. Based on the existing research results [16–20], this paper enriched and perfected the influencing factors of FEE. The specific indicators were as follows: EDL was the level of economic development by per capita GDP. IEC was the investment in forestry ecological construction, which was specifically expressed by the investment in forestry fixed assets [21]. ERL was used to represent environmental regulation level, the investment in wastewater treatment and the investment in waste gas treatment were selected, and then entropy method was used to calculate the two weights, as a result, the environmental regulation level of each region can be calculated. The industrial structure was represented by ISL, and it was expressed by calculating the proportion of regional output value of forestry to total forest output value. FRE was used to express resource endowment, and forest resource coverage was used to reflect it. The dependent variable was forestry ecological efficiency, which was expressed by FEE. The all indicators are from China Statistical Yearbook, China Forestry Statistical Yearbook, China Energy Statistical Yearbook and China Environment Statistical Yearbook, the missing data are estimated by simple linear extrapolation method.

2.2. Malmquist Index Malmquist productivity index model was first proposed by Malmquist in 1953 [22], and then Fare et al. [23] redefined it as a multi-stage dynamic production efficiency evalua- tion method guided by production efficiency. Malmquist is a total factor productivity index based on panel data, it evaluates the efficiency of input and output on the basis of distance Forests 2021, 12, 142 4 of 17

function. This model gives the distance function from two different angles of input and output, input vector can reduce the level of production frontier and evaluate the efficiency of production technology [24]. The change rate of technical efficiency can be divided into scale efficiency change rate and pure technical efficiency change rate, the change index of scale efficiency represents the change of forestry production level, and the pure technical efficiency change index only represents the change of allocation and utilization level of forestry production factor resources [25]. If the change rate of technical efficiency is more than 1, it means that the forestry production efficiency is improved in the period from t to t + 1. If the change rate of technical efficiency is less than 1, it means that the forestry production efficiency has decreased in this period; if the change rate of technical efficiency is equal to 1, it means that the forestry production efficiency has not changed in this period. When the returns to scale remain unchanged, the Malmquist productivity index formula in the period from t to t + 1 can be written as follows as in Equation (1).

" # 1 Dtxt+1, yt+1 Dt+1xt+1, yt+1 2 t+1 = o o × o o Mo t t t t+1 t t (1) D (xo, yo) D (xo, yo)

t t t  In Equation (1), D xo, yo represents the technical efficiency level of the current period; t t+1 t+1 D xo , yo represents the technical efficiency level of phase t + 1 with the technology of t+1 t t  phase t; D xo, yo represents the technical efficiency level of phase t with the technology t+1 t+1 t+1 of phase t + 1; D xo , yo represents the technical efficiency level of phase t + 1 with the technology of phase t + 1.

2.3. SBM Model Slacks-based measure (SBM) model is one of the DEA models which was first proposed by Tone [26]. The DEA methods are mainly applied to traditional models, such as constant return to scale (CCR) model and variable return to scale (BBC) models [27], and the output of these models is mostly based on the expected output without considering the problem of unexpected output. The SBM model can take the unexpected output into account, and it has two measurement methods including radial and non-radial [28]. SBM model not only has the characteristics of the traditional CCR model, but also has higher accuracy in the case of multiple inputs. Because the model takes both slack input and slack output into account, so it is more effective to evaluate the actual efficiency level [29]. Decision making unit (DMU) is an operational entity that can transform certain input into corresponding output [30]. Some DMU may be calculated relatively effective in SBM model, that means the efficiency value of these DMU is 1. Therefore, it is impossible to compare the production efficiency of these effective DMU. Therefore, this paper combines super efficiency with SBM model, that is super efficiency SBM model [31]. The efficiency value of the effective DMU can be calculated more than 1 in this model, and the gap between the efficiency values measured by this model is larger, so it is better to campare the efficiency values of each DMU [32]. Therefore, this paper uses the super efficiency SBM model to calculate and analyze the FEE of 31 provinces in China.

2.4. Moran’s I Index This paper analyzes the FEE of 31 provinces and cities in China from 2004 to 2018 by using spatial econometrics, the spatial autocorrelation of FEE is studied to find out the specific distribution law of FEE in China. Generally speaking, Moran’s I index is used to judge the existence of spatial correlation. Moran’s I can be divided into global Moran’s I and local Moran’s I, global Moran’s I is used to test global spatial correlation, while local Moran’s I is used to measure local spatial correlation [33]. The expression for Moran’s I is as follows: Forests 2021, 12, 142 5 of 17

N N N ∑ ∑ Wij(yi − y)(yj − y) i=1 j=1 Moran0s I = , (i 6= j) (2) N N N 2 ( ∑ ∑ Wij) ∑ (yi − y) i=1 j=1 i=1

In Equation (2), yi and yj represent the observed values of i and j areas respectively, N is the number of spatial units, which is the mean value of FEE in each region, and Wij is the spatial weight matrix. If the province i and j is adjacent, then Wij = 1; if not, then Wij = 0. The range of Moran’s I index is between [−1, +1]. If Moran’s I value is significantly greater than 0, it means that there is a spatial positive correlation; if Moran’s I value is significantly less than 0, then there is a negative correlation. If Moran’s I value is equal to 0, it is spatial uncorrelated [34]. Local spatial autocorrelation is used to determine whether a certain attribute value of each space has spatial correlation locally. The local Moran’s I index is presented below in Equation (3): N(y − y) N 0 = i ( − )2 Moran s Ii 2 ∑ Wij yi y (3) S j=1

N 2 1 2 S = ∑ (yi − y) (4) N i=1 In Equation (3), if Moran’s I is positive, it means that there is similar spatial agglom- eration around the regional unit; if Moran’s I is negative, it means that there is spatial agglomeration with non-similar values around the regional unit. The higher the Moran’s I index is, the higher the degree of closeness is [35]. The biggest drawback of global Moran’s I is that it can’t analyze the spatial agglomeration characteristics between local areas, how- ever, local Moran’s I just makes up for this vacancy. It can reflect greatly whether there is spatial dependence between a certain area and its surrounding areas [36]. Therefore, this paper uses the local Moran’s I index to analyze the local spatial correlation of China’s FEE. The results of local Moran’s I can be divided the local spatial characteristics into four types: the first quadrant is the high-high (H-H) aggregation area, which means that the level of FEE in this region and its adjacent areas is high, and the spatial correlation is high-level area. The second quadrant is the low-high (L-H) aggregation area, which means that the ecological efficiency of forestry in this region is low, but the adjacent areas are high, so the spatial correlation is in the development stage. The third quadrant is the low- low (L-L) aggregation area, which means that the ecological efficiency of forestry in this region and its adjacent areas is low. Therefore, the spatial correlation is low efficiency area. The fourth quadrant is the high-low (H-L) cluster area, which means that the ecological efficiency of forestry in this region is higher than that of adjacent areas, and the spatial performance is spillover effect [37]. Generally, H-H type and L-L type mean that there is positive spatial correlation between regions, while L-H type and H-L type indicate that there is spatial negative correlation between regions.

2.5. Tobit Regression Tobit model is a model in which dependent variables satisfy certain constraints. Tobit model includes two kinds of models: discrete dependent variable model and restricted dependent variable model [38]. The restricted dependent variable is the case that the dependent variable is restricted in practice [39]. In this case, the sample data obtained from a subset of the population may not fully reflect the population. The basic definition model of Tobit model is Equation (5).

 FEE, i f FEE > 0 (5) 0, i f FEE ≤ 0 Forests 2021, 12, 142 6 of 17

The value of FEE in this paper is such a dependent variable, and the review regression model of the restricted dependent variable model is used for regression analysis. The econometric model is shown in Equation (6).

FEE = β0 + β1 EDL + β2 IEC + β3 ERL + β4 ISL + β5 FRE + ε (6)

In Equation (6), β0 is the intercept term, β1, β2, β3, β4, β5 are the independent variable coefficients of the parameter to be estimated, and ε is the random interference term. EDL was the level of economic development by per capita GDP, IEC was expressed by the investment in forestry fixed assets, ERL represents environmental regulation level, ISL means the industrial structure, FRE was used to express resource endowment level. The dependent variable FEE represents forestry ecological efficiency.

2.6. Random Forest Generally, ensemble learning is a method that combines multiple learners to improve the performance of solving a problem. In the 1980s, Breiman et al. [40] came up with a classification tree algorithm after a long time of research. The data can be classified precisely through this classification tree. Later, Breiman combined the classification trees previously into random forest [41], that is taking the decision tree as the basis, the learner uses bagging to do simple voting on the classified objects based on decision tree, so it is a classifier combined with decision tree based on learners as an algorithm in ensemble learning [42]. In order to make the effect of processing data better, the individual learner’s discrimination effect should not be too bad. Random forest is random in the selection of attributes, and then select an optimal attribute to divide these samples. The principle of random forest is as follows: k samples are extracted from the sample set D (bootstrap is put back for sampling) for training, and then the k decision trees are combined into random forests, so that these decision trees can judge the classification results of each sample separately [43]. If there are many voting results of a certain classification, it will be divided into this category first. Random forest algorithm has many excellent characteristics in application. Generally speaking, random forest has high prediction accuracy in dealing with classification and regression problems [44]. When the data are missing or unbalanced, the random forest algorithm has good adaptability and the effect is not bad. When we increase the number of learners, the prediction error in random forest will gradually decrease and finally converge to a certain level [45]. Moreover, it can judge the importance of thousands of explanatory variables, it performs well on the data set, and when multiple randomness is introduced, it will not easily lead to over-fitting [46]. In addition, random forest is usually used to deal with classification problems, but it can also be used to deal with regression problems [47]. In general, linear regression model analysis in this study, there are explanatory variables and explained variables. In random forest regression model, each factor can be regarded as explanatory variable. In this paper, EDL, IEC, ERL, ISL, and FRE are used as the explanatory variables, and FEE is the explanatory variable. In random forest model, the importance of variables is evaluated by IncMSE and IncNodePurity. The IncMSE index is actually the abbreviation of increase in MSE, it represents the increasement of mean square error (MSE). When variable in the sample is assigned arbitrarily, if the variable is very important, the prediction error would increase. Therefore, the greater the value of IncMSE is, the more important it is [47]. IncNodePurity index is the abbreviation of increase in node purity. It refers to the error value of leaf nodes. In dealing with regression problems, the increase of node purity is equivalent to the decrease of Gini index, and the larger the IncNodePurity value is, the more important the factor is [48], which these variables should be focused on. Forests 2021, 12, 142 7 of 17

3. Results 3.1. Analysis of Total Factor Productivity of China’s Forestry In this paper, Malmquist productivity index model is used to measure the total factor productivity change (TFPC) of China’s provincial forestry production in 14 periods from 2004 to 2018. In order to better analyze the regional characteristics of TFPC in China, China is divided into three regions according to the National Bureau of Statistics. The eastern region includes 12 provinces, municipalities and autonomous regions, including Beijing, Tianjin, Hebei, Liaoning, Shanghai, Jiangsu, Zhejiang, Fujian, Shandong, Guang- dong, Guangxi, and Hainan while the central region includes Shanxi, Inner Mongolia, Jilin, Heilongjiang, Anhui, Jiangxi, Henan, Hubei, and Hunan, and the western region includes Chongqing, Sichuan, Guizhou, Yunnan, Tibet, Shaanxi, Gansu, Ningxia, Qinghai, and Xinjiang. According to Table2, the average growth rate of Malmquist index of TFPC from 2004 to 2018 was 1.012, and the average growth values of ecological efficiency in eastern, central and western regions were 0.994, 1.011 and 1.038, which showed that the average TFPC in the western region was greater than that of the eastern and central regions. However, the overall trend was not stable enough, and the fluctuation was obvious. Noticeably, the TFPC of China, the east and the central regions from 2008 to 2009 reached the lowest point, which were 0.845, 0.783, and 0.834, that showed the TFPC in this year has been greatly deccreased compared with the previous period, and the forestry industry has suffered a huge strike. However, the TFPC of China’s forestry showed a fluctuating growth trend from 2012 to 2018, and the values were all greater than 1, which indicated that forestry was developing in a benign direction. In the western region during 2012–2013, the TFPC of west forestry reached the maximum value during 2015–2016, but then began to decline. Generally speaking, China’s forestry TFPC was in an upward trend and the fluctuation range of TFPC in the eastern and western regions was greater than that in the central region.

Table 2. Total factor productivity change (TFPC) of forestry in China from 2004 to 2018.

Time All East Central West 2004–2005 0.942 0.926 0.927 0.986 2005–2006 1.058 1.066 1.045 1.020 2006–2007 1.042 1.099 1.055 0.925 2007–2008 1.011 0.962 1.052 1.007 2008–2009 0.845 0.783 0.834 0.952 2009–2010 1.033 1.010 1.044 0.957 2010–2011 0.941 0.857 0.910 1.117 2011–2012 0.983 0.980 1.042 0.988 2012–2013 1.078 1.075 1.090 1.111 2013–2014 1.039 1.035 1.067 1.064 2014–2015 1.026 1.021 0.996 1.054 2015–2016 1.050 0.997 0.988 1.265 2016–2017 1.056 1.040 1.046 1.148 2017–2018 1.090 1.124 1.088 0.997 Average 1.012 0.994 1.011 1.038

3.2. Analysis of China’s FEE 3.2.1. Descriptive Analysis of China’s FEE There is the static evaluation of FEE of 31 provinces and cities in China in 2004, 2011 and 2018. Which was measured in super efficiency SBM model with DEA-SOLVER Pro 5.0 software. The FEE value and ranking of each region were shown in the Table3 below: Forests 2021, 12, 142 8 of 17

Table 3. Forestry ecological efficiency values of 31 provinces and cities in China in 2004, 2011 and 2018.

Region 2004 FEE Rank 2011 FEE Rank 2018 FEE Rank Beijing 1.000 17 1.000 14 1.123 9 Tianjin 5.105 1 1.649 3 1.669 3 Hebei 0.112 26 0.072 25 0.103 23 Shanxi 0.030 31 0.013 29 0.041 27 InnerMongolia 1.427 6 0.135 23 0.048 26 Liaoning 0.302 21 0.195 21 0.293 18 Jilin 1.089 11 1.059 11 0.129 22 Heilongjiang 1.056 12 0.235 19 0.102 24 Shanghai 3.489 2 1.647 4 1.587 4 Jiangsu 1.045 13 1.001 13 1.014 14 Zhejiang 1.243 7 0.306 18 0.380 17 Anhui 1.132 8 1.194 6 1.158 8 Fujian 1.443 5 1.097 10 1.162 7 Jiangxi 1.127 9 1.001 12 1.016 13 Shandong 0.248 22 0.402 15 0.455 16 Henan 1.018 16 0.325 16 0.265 19 Hubei 0.203 24 0.311 17 0.249 20 Hunan 1.039 14 1.182 8 1.027 11 Guangdong 0.520 20 1.190 7 1.300 5 Guangxi 1.031 15 1.540 5 1.721 2 Hainan 3.458 4 4.219 2 3.651 1 Chongqing 0.068 29 0.094 24 1.021 12 Sichuan 0.126 25 0.155 22 0.224 21 Guizhou 0.225 23 0.204 20 1.042 10 Yunnan 1.102 10 1.164 9 1.209 6 Tibet 3.484 3 6.555 1 0.996 15 Shaanxi 0.088 28 0.050 27 0.011 28 Gansu 0.031 30 0.013 28 0.007 29 Qinghai 0.999 18 0.010 30 0.003 31 Ningxia 0.999 19 0.001 31 0.006 30 Xinjiang 0.100 27 0.055 26 0.064 25

In terms of year, the highest value of FEE in 2004 was Tianjin, and the value was 5.105. The ecological efficiency of forestry in Shanghai and Tibet came next, the values were 3.489 and 3.484. The three regions of the lowest FEE were Chongqing, Gansu and Shanxi. The values were 0.068, 0.031, and 0.030, it can be seen that the gap of FEE in different regions was very large in 2004. Tianjin was ranked first in forestry ecology in 2004 has fallen behind to the third place in 2011, and the FEE ranking of Tibet and Hainan Province has increased, their values were 6.555 and 4.219 respectively. The three regions of the lowest FEE were Shanxi, Qinghai and Ningxia. In 2018, the FEE of Hainan province has risen to the first place, and the Guangxi with the rapid increase of FEE was ranked second. The regions of the lowest FEE were Shaanxi, Gansu, Ningxia, and Qinghai. In terms of region, Tianjin, Shanghai and Hainan were the areas with high level of FEE. The FEE of Tibet ranked third in 2004 and first in 2011, but the FEE has been far lower than other regions in 2018, and its ranking was behind 15. Hebei, Shanxi, Shaanxi, Gansu, Ningxia were at the bottom of the list in the three periods. The regions in the middle level of forest ecological efficiency in the three periods included Liaoning, Jilin, Heilongjiang, Jiangsu, Shandong, Henan, Hubei, Hunan, Sichuan, Guizhou, etc. The areas where the relative efficiency of forestry ecology has been greatly improved included Beijing, Guangdong, Chongqing, and so on.

3.2.2. Spatial Distribution Analysis of China’s FEE In order to study the spatial characteristics of FEE in China, this paper made a spatial analysis of the above ecological efficiency values of forestry industry in 2004, 2011, and 2018 by using ArcGIS software. The spatial distribution is presented in Figures1–3. Forests 2021, 12, x FOR PEER REVIEW 9 of 18

In terms of year, the highest value of FEE in 2004 was Tianjin, and the value was 5.105. The ecological efficiency of forestry in Shanghai and Tibet came next, the values were 3.489 and 3.484. The three regions of the lowest FEE were Chongqing, Gansu and Shanxi. The values were 0.068, 0.031, and 0.030, it can be seen that the gap of FEE in different regions was very large in 2004. Tianjin was ranked first in forestry ecology in 2004 has fallen behind to the third place in 2011, and the FEE ranking of Tibet and Hainan Province has increased, their values were 6.555 and 4.219 respectively. The three regions of the low- est FEE were Shanxi, Qinghai and Ningxia. In 2018, the FEE of Hainan province has risen Forests 2021, 12, 142 to the first place, and the Guangxi with the rapid increase of FEE was ranked second.9 The of 17 regions of the lowest FEE were Shaanxi, Gansu, Ningxia, and Qinghai. In terms of region, Tianjin, Shanghai and Hainan were the areas with high level of FEE. The FEE of Tibet ranked third in 2004 and first in 2011, but the FEE has been far lower than otherIn Figure regions1, China’s in 2018 FEE, and of its the ranking northeast was and behind the south 15. Hebei, was higher, Shanxi, and Shaanxi, FEE of China’sGansu, Ningxiamiddle regionswere at the was bottom relatively of the lower. list in The the FEE three of periods. Shanghai The and regions Tianjin in ranked the middle top inlevel the ofeast forest region, ecological the reason efficiency was thatin the forestry three periods area and included forestry Liaoning, practitioners Jilin, had Heilo nongjiang, serious Jiangsuredundancy., Shandong, Moreover, Henan, the Hubei, total forestry Hunan, output Sichuan, value Guizhou, of Shanghai etc. The reached areas 13.14 where trillion the relativeyuan, the efficiency timber outputof forestry of Tianjin ecology was has about been 40greatly thousand improved cubic included meters. TheBeijing, FEE Guang- in Tibet dong,was the Chongqing, highest, becauseand so on. the emissions of waste gas and wastewater in Tibet were the lowest in the whole country, and there was also no redundancy of forestry practitioners. 3.2.2.Therefore, Spatial the D FEEistribution in this regionAnalysis was of relativelyChina’s FEE high. The high forestry efficiency of Hainan province was related to the low energy investment and the high output value of forestry. In order to study the spatial characteristics of FEE in China, this paper made a spatial The central regions of low ecological efficiency were Hubei, Hunan, Shaanxi, and Hebei, analysis of the above ecological efficiency values of forestry industry in 2004, 2011, and because there was redundancy in various input elements including forestry planting area, 2018 by using ArcGIS software. The spatial distribution is presented in Figures 1–3. energy input, emissions of waste gas and wastewater.

Forests 2021, 12, x FOR PEER REVIEW 10 of 18

Figure 1. Distribution of forestry ecological efficiency (FEE) in China in 2004. Figure 1. Distribution of forestry ecological efficiency (FEE) in China in 2004.

Figure 2. Distribution of forestry ecological efficiency (FEE) in China in 2011. Figure 2. Distribution of forestry ecological efficiency (FEE) in China in 2011.

Figure 3. Distribution of forestry ecological efficiency (FEE) in China in 2018.

In Figure 1, China’s FEE of the northeast and the south was higher, and FEE of China’s middle regions was relatively lower. The FEE of Shanghai and Tianjin ranked top in the east region, the reason was that forestry area and forestry practitioners had no seri- ous redundancy. Moreover, the total forestry output value of Shanghai reached 13.14 tril- lion yuan, the timber output of Tianjin was about 40 thousand cubic meters. The FEE in Tibet was the highest, because the emissions of waste gas and wastewater in Tibet were the lowest in the whole country, and there was also no redundancy of forestry practition- ers. Therefore, the FEE in this region was relatively high. The high forestry efficiency of Hainan province was related to the low energy investment and the high output value of forestry. The central regions of low ecological efficiency were Hubei, Hunan, Shaanxi, and

Forests 2021, 12, x FOR PEER REVIEW 10 of 18

Forests 2021, 12, 142 10 of 17 Figure 2. Distribution of forestry ecological efficiency (FEE) in China in 2011.

Forests 2021, 12, x FOR PEER REVIEW 11 of 18

Hebei, because there was redundancy in various input elements including forestry plant- ing area, energy input, emissions of waste gas and wastewater. Figure 3. Distribution of forestry ecological efficiency (FEE) in China in 2018. FigureIn F3.igure Distribution 2, the spatialof forestry differentiation ecological efficiency feature ( FEE)of FEE in China in China in 20 was18. high in southeast and south but low in northwest except Tibet. Compared with 2004, the ecological effi- In Figure2, the spatial differentiation feature of FEE in China was high in southeast ciency of forestry in Northwest China remained unchanged except Qinghai. The low value andIn south Figure but low1, China’s in northwest FEE of except the northeast Tibet. Compared and the with south 2004, wa thes higher ecological, and efficiency FEE of of FEE in Qinghai was mainly due to the increase of input from various factors and the China’sof forestry middle in Northwest regions wa Chinas relatively remained lower unchanged. The FEE except of Shanghai Qinghai. and The Tianjin low value ranked of FEEtop output has not increased significantly, which leaded to the decrease of FEE. The ranking inin the Qinghai east region, was mainly the reason due wa to thes that increase forestry of area input and from forestry various practitioners factors and ha thed no output seri- of FEE in Tibet, Shanghai, Tianjin, and Hainan remained unchanged, all of which ranked oushas redundancy. not increased Moreover significantly,, the whichtotal forestry leaded tooutput the decrease value of of Shanghai FEE. The reache rankingd 13.14 of FEE tril- in at the top. lionTibet, yuan Shanghai,, the timber Tianjin, output and of Hainan Tianjin remained was about unchanged, 40 thousand all of cubic which meters. ranked The at theFEE top. in As can be seen from Figure 3, China’s FEE was basically higher in the south in 2018. Tibet Aswas can the be highest, seen from because Figure the3, China’semissions FEE of waswaste basically gas and higher wastewate in ther southin Tibet in were 2018. It was worth noting that Tibet, which has been greatly reduced due to the reduction of theIt was lowest worth in the noting whole that country, Tibet, and which there has wa beens also greatly no redundancy reduced due of f toorestry the reduction practition- of timberers.timber Therefore, production. production. the BecauseFEE Because in this of ofthe region the dec decrease reasewas relativelyof of timber timber high.output output The and andhigh the the forestryincrease increase efficiency of of forestry forestry of fixedHainanfixed assets assets province investment investment was relatedin in the the three to three the northeast northeastlow energ provinces, provinces,y investment the the forestry and forestry the production high production output efficiency efficiencyvalue of wasforestrywas relatively relatively. The creduced.entral reduced. region Shanghai, Shanghai,s of low Tianjin, Tianjin,ecological and and efficiencyHainan Hainan still stillwere maintained maintained Hubei, Hunan, a aleading leading Shaanxi, position position and inin terms terms of of FEE. FEE.

3.2.3.2.3.3. Characteristi Characteristicscs of of Spatial Spatial Scatter Scatterofof FEE FEE in in China China Moran’sMoran’s I scatter I scatter plot plot is isan an important important tool tool to to analyze analyze local local spatial spatial correlation. correlation. The The FigFigureure 44 show showeded the the local local Moran’s Moran’s I scatter I scatter diagram diagram of of China’s China’s FEE FEE agglomeration agglomeration in in 20 2018.18.

Figure 4. Moran’s I scatter chart of provincial forestry ecological efficiency (FEE) cluster in 2018. Figure 4. Moran’s I scatter chart of provincial forestry ecological efficiency (FEE) cluster in 2018.

The result showed Moran’s I index of ecological efficiency of China’s 31 provinces was 0.222, which indicated that the ecological efficiency of 31 provinces in China had a significant positive autocorrelation in spatial distribution, and most regions were distrib- uted in the H-H cluster area and the L-L cluster area. As depicted in Figure 5, the H-H cluster area of China’s FEE in 2018 included Hainan, Guangxi, Guangdong, Yunnan, Fujian, Guizhou, Hunan, Anhui, Jiangsu, Jiangxi, and Shanghai. These provinces were mostly located in the southern coastal areas with a more developed development level, and the areas had reasonable industrial layout and struc- ture level, superior location advantages, and more advanced environmental protection concepts. So, the ecological efficiency has been continuously improved, which had a great radiation and diffusion effect on the surrounding areas.

Forests 2021, 12, 142 11 of 17

The result showed Moran’s I index of ecological efficiency of China’s 31 provinces was 0.222, which indicated that the ecological efficiency of 31 provinces in China had a sig- nificant positive autocorrelation in spatial distribution, and most regions were distributed in the H-H cluster area and the L-L cluster area. As depicted in Figure5, the H-H cluster area of China’s FEE in 2018 included Hainan, Guangxi, Guangdong, Yunnan, Fujian, Guizhou, Hunan, Anhui, Jiangsu, Jiangxi, and Shanghai. These provinces were mostly located in the southern coastal areas with a more developed development level, and the areas had reasonable industrial layout and structure Forests 2021, 12, x FOR PEER REVIEW 12 of 18 level, superior location advantages, and more advanced environmental protection concepts.

So, the ecological efficiency has been continuously improved, which had a great radiation and diffusion effect on the surrounding areas.

Figure 5. Clustering map of China’s provincial regional forestry ecological efficiency (FEE) in 2018. Figure 5. Clustering map of China’s provincial regional forestry ecological efficiency (FEE) in 2018. The L-H cluster area included Zhejiang, Henan, Hubei, and Xinjiang. These four regionsThe wereL-H cluster relatively area scattered included with Zhejiang, the obvious Henan, location Hubei, character and Xinjiang. of being These diffused, four re- the gionsecological were efficiencyrelatively wasscattered greatly with improved, the obvious and thelocation speed character of improvement of being wasdiffused, fast. the ecologicalThe L-Lefficiency cluster was area greatly included im theproved, economically and the s underdevelopedpeed of improvement northwest was fast. and north- eastThe provinces L-L cluster in 2018, area including included Gansu, the economic Heilongjiang,ally underdeveloped Qinghai, Ningxia, northwest Jilin, Liaoning, and northeastInner Mongolia, provinces Shaanxi, in 2018, Shanxi,including Shandong, Gansu, Heilongjiang, and Sichuan. Qinghai, We found Ningxia, that Jilin, L-L clusterLiao- ning,area Inner was also Mongolia, the agglomeration Shaanxi, Shanxi, area Shand with lowong, FEEand bySichuan. looking We at found Figures that3 andL-L 5cluster. The areanorthwest was also was the rich agglomeration in resources, area however, with low the traditionalFEE by looking forestry at Figure cultivations 3 and mode 5. The and northwestthe environmental was rich pollutionin resources, brought however, by forestry the traditional investment forestry made thecultivation ecological mode efficiency and theof forestryenvironmental in these pollution areas became brough thet by state forestry of L-L inv agglomeration.estment made These the ecological factors made effi- it ciencydifficult of forestry to change in thethese state areas of L-Lbecame concentration the state of in L- aL shortagglomeration. time. These factors made it difficultThe to H-L change cluster the area state mainly of L-L includedconcentration Tianjin in a and short Beijing, time. Tibet, Chongqing and Hebei.The BeijingH-L cluster and Tianjinarea mainly have included made remarkable Tianjin and achievements Beijing, Tibet, in Chongqing social and economicand He- bei.development Beijing and as Tianjin the regions have withmade high remarkable openness achiev in theements north, andin social they haveand economic a high level de- of velopmentFEE in 2018. as the However, regions due with to high the lackopenness of effective in the regionalnorth, and cooperation they have mechanism a high level and of FEEincomplete in 2018. industrialHowever, chain due to layout the lack with of surrounding effective regional areas, there coop waseration no effectivemechanism radiation and incompleteand driving industrial effect on chain the surrounding layout with areas,surroundin whichg resultedareas, there in the was phenomenon no effective thatradia- the tionecological and driving efficiency effect of on the the regions surrounding was high areas, and which the adjacent resulted areas’ in the ecological phenomenon efficiency that thewas ecological low. efficiency of the regions was high and the adjacent areas’ ecological effi- ciencyIt was can low. be seen that more than 70% of the regions were distributed in the first and third quadrants,It can be and seen only that lessmore than than 30% 70% of of the the provinces regions were were distributed distributed in the in the first second and third and quadrants, and only less than 30% of the provinces were distributed in the second and fourth quadrants. This showed that most of China’s high-tech industrial agglomeration in 2018 had obvious positive spatial correlation characteristics, and a few provinces and cit- ies had a negative spatial autocorrelation.

3.3. Analysis of Influencing Factors Traditional forestry only considered the output of forestry in the process of produc- tion, but it was unreasonable that the influence of resources, economy, environment and other factors was ignored. In order to study the improvement of ecological efficiency in forestry, it was necessary to study the factors affecting the ecological efficiency of forestry.

Forests 2021, 12, 142 12 of 17

fourth quadrants. This showed that most of China’s high-tech industrial agglomeration in 2018 had obvious positive spatial correlation characteristics, and a few provinces and cities had a negative spatial autocorrelation.

3.3. Analysis of Influencing Factors Traditional forestry only considered the output of forestry in the process of production, Forests 2021, 12, x FOR PEER REVIEWbut it was unreasonable that the influence of resources, economy, environment and13 other of 18 factors was ignored. In order to study the improvement of ecological efficiency in forestry, it was necessary to study the factors affecting the ecological efficiency of forestry.

3.3.1.3.3.1. AnalysisAnalysis ofof TobitTobit RegressionRegression InIn thisthis paper,paper, thethe influencinginfluencing factorsfactors ofof FEEFEE ofof 3131 regionsregions werewere discussed.discussed. ThereThere werewere 465465 samplessamples andand thethe FEEFEE valuevalue waswas calculatedcalculated byby thethe supersuper efficiencyefficiency SBMSBM modelmodel fromfrom 20042004 toto 2018.2018. TheThe resultresult isis shownshown inin FigureFigure6 6..

Figure 6. Comparison of regional forestry ecological efficiency (FEE) in China from 2004 to 2018. Figure 6. Comparison of regional forestry ecological efficiency (FEE) in China from 2004 to 2018. It can be seen that the fitting results of this model was good in Table4, and the It can be seen that the fitting results of this model was good in Table 4, and the cor- correlation was significant at the level of p = 0.05, which indicated that EDL, IEC, ERL, ISLrelation and FREwas affectedsignificant the at FEE the to level a certain of P = extent. 0.05, which The detailed indicated analysis that EDL, of each IEC, factor ERL, was ISL asand follows: FRE affected the FEE to a certain extent. The detailed analysis of each factor was as follows: Table 4. Regression results of influencing factors. Table 4. Regression results of influencing factors. Variabe Intercept EDL IEC ERL ISL FRE Variabe Intercept EDL IEC ERL ISL FRE −3 −2 −4 −5 0 −3 Coeffificient −6.695 × 10 Coeffificient2.197 × 10 −6.695*** ×10−−3 3.2532.197××1010−2 ***** −3.2532.146××10−104 ** ***2.146 ×6.39910−5 ***× 106.399*** ×100 ***3.967 3.967× 10×10−3*** *** Note:Note: **, and **, *** and denote *** statisticaldenote stati signifificancestical signifificance at the 5% and at 1%the levels, 5% and respectively. 1% levels, respectively.

TheThe resultresult showed showed that that the the economic economic development development level level was was positively positively correlated correlated with thewith FEE. the TheFEE. higher The higher the level the of level economic of economic development development was, the was, higher the the higher value the of value forestry of ecologicalforestry ecological performance performance was. With was. the expansionWith the expansion of economic of economic scale and thescale improvement and the im- ofprovement people’s incomeof peopl level,e’s income people level, had peopl highere had demand higher for demand forestry for production forestry production efficiency. Therefore,efficiency. speedingTherefore, up speeding economic up developmenteconomic development was an important was an important factor to improvefactor to theim- ecologicalprove the efficiencyecological ofefficiency forestry. of forestry. ThereThere waswas aa negativenegative correlationcorrelation betweenbetween thethe investmentinvestment inin forestryforestry ecologicalecological con- structionstruction and the the FEE FEE value. value. This This showed showed that that the the higher higher the theinvestment investment in forestry in forestry eco- ecologicallogical construction construction was, was, thethe lower lower thethe value value of FEE of FEE was. was. It was It was true true that that forestry forestry ecolog- eco- logicalical cons constructiontruction investment investment has has provided provided impetus impetus for for th thee growth growth of of forestry forestry in industry,dustry, butbut redundancyredundancy of of forestry forestry ecological ecological construction construction investment investment in forestry in forestry production production would would reduce the FEE, ecological construction would also increase the environmental bur- den. When the pollution exceeded the of the environment, environmen- tal pollution would inevitably affect the improvement of forestry ecological productivity and economic development. Therefore, the coefficient of investment in forestry ecological construction was negative. The relationship between environmental regulation factors and forestry ecology was very significant. The FEE calculated in this paper took the environmental pollution factors into account. The purpose of environmental regulation was to reduce the unexpected out- put such as wastewater and waste gas, because reducing the pollution caused by forestry

Forests 2021, 12, 142 13 of 17

Forests 2021, 12, x FOR PEER REVIEW 14 of 18 reduce the FEE, ecological construction would also increase the environmental burden. When the pollution exceeded the carrying capacity of the environment, environmental pollution would inevitably affect the improvement of forestry ecological productivity and economic development. Therefore, the coefficient of investment in forestry ecological production can improve the ecological efficiency of forestry. Therefore, the local govern- construction was negative. ments should increase the investment in environmental governance to improve the eco- The relationship between environmental regulation factors and forestry ecology was logical efficiency of forestry. very significant. The FEE calculated in this paper took the environmental pollution factors There was a strong positive correlation between the industrial structure and the eco- into account. The purpose of environmental regulation was to reduce the unexpected logical efficiency of forestry, and the correlation coefficient was as high as 6.399, indicating output such as wastewater and waste gas, because reducing the pollution caused by forestrythat the productionindustrial structure can improve of forestry the ecological had a pos efficiencyitive impact of forestry. on improvi Therefore,ng the ecological the local governmentsefficiency of forestry. should increase Therefore, the on investment the one hand, in environmental the government governance accelerated to improve the upgrad- the ecologicaling of industrial efficiency structure of forestry. and vigorously developed the forest industry. −3 ThereThe correlation was a strong coefficient positive between correlation FRE between and FEE the was industrial 3.967×10 structure. Generally and speaking, the eco- logicalthe better efficiency the resource of forestry, endowment and the correlationwas, the higher coefficient the FEE was was. as high However, as 6.399, the indicating resource thatendowment the industrial cannot structure be changed of forestry in a short had time a positive. impact on improving the ecolog- ical efficiency of forestry. Therefore, on the one hand, the government accelerated the upgrading3.3.2. Importance of industrial Analysis structure of Factor ands vigorously developed the forest industry. TheThe correlationre was an analysis coefficient of the between importance FRE and of FEEfactors was affec 3.967ting× 10forestry−3. Generally ecologica speak-l effi- ing,ciency the by better using the random resource forest. endowment The importance was, the higher of factors the FEE can was. be ranked However, by theIncMSE resource and endowmentIncNodePurity cannot, and be the changed results inare a shown short time. in Figure 7 and Table 5.

3.3.2.Table Importance5. Analysis results Analysis of IncMSE of Factors and IncNodePurity.

There was an analysisEDL of the importanceIEC of factorsERL affecting forestryISL ecologicalFRE effi- ciency by%IncMSE using random forest.47.837 The importance53.186 of factors46.328 can be ranked51.520 by IncMSE33.890 and IncNodePurity,IncNodePurity and the results9.442 are shown12.972 in Figure 7 and12.413 Table 5. 11.270 5.969

FigureFigure 7.7.Scatter Scatter diagramdiagram ofof IncMSEIncMSEand and IncNodePurity.IncNodePurity.

TableI 5.t canAnalysis be concluded results of IncMSE that IEC and had IncNodePurity. the highest value in both IncMSE and IncNodePu- rity, which indicated that the variable had an important impact on FEE. Then, ERL and EDL IEC ERL ISL FRE ISL came next. It can be seen that these two factors were a little important in affecting the ecological%IncMSE efficiency 47.837 of forestry. In 53.186 contrast, FRE 46.328 was the lowest 51.520 important factor 33.890 among theseIncNodePurity indicators, and9.442 the values 12.972of IncMSE and 12.413 IncNodePurit 11.270y were 33.89 5.969and 5.969, which showed that forestry resource endowment was not particularly significant to FEE. It can be concluded that IEC had the highest value in both IncMSE and IncNodePurity, which4. Discussion indicated that the variable had an important impact on FEE. Then, ERL and ISL cameThe next. result It can of the be seenMalmquist that these index two of China’s factors wereforestry a little efficiency important only was in affecting 1.012, whic theh ecologicalreflected the efficiency growth ofrate forestry. of China’s In contrast, forestry FREtotal wasfactor the is lowestslow, t importanthe result is factorin accordance among thesewith indicators,the study by and Liu the and values Li [ of49 IncMSE]. However, and IncNodePuritythe TFPC of China’s were 33.89 forestry and 5.969,from 2008 which to showed2009 reached that forestry the lowest resource value endowment in recent years. was Due not particularlyto severe low significant temperature, to FEE. rain, snow and freezing disaster occurred in southern China at the beginning of 2008, and Wenchuan earthquake occurred in May at the same year. What’s worse was that the international

Forests 2021, 12, 142 14 of 17

4. Discussion The result of the Malmquist index of China’s forestry efficiency only was 1.012, which reflected the growth rate of China’s forestry total factor is slow, the result is in accordance with the study by Liu and Li [49]. However, the TFPC of China’s forestry from 2008 to 2009 reached the lowest value in recent years. Due to severe low temperature, rain, snow and freezing disaster occurred in southern China at the beginning of 2008, and Wenchuan earthquake occurred in May at the same year. What’s worse was that the international financial crisis continued to spread in the second half of the year, the forest product market and export fluctuated significantly, and the TFPC of the whole region of China reduced greatly. The results of this paper also reflected China’s FEE was basically higher in the south and lower in the northwest in recent years. The result was roughly the same as Luo et al.’s study [50]. Similarly, the H-H cluster area of China’s FEE in 2018 was in the south, which the L-L cluster area included northwest regions. It was not conducive to the sustainable development of forestry in northwest regions. Our results also indicated that there was a negative correlation between the investment in forestry ecological construction and the FEE value. This is because the redundancy of investment in forestry fixed assets would also increase the environmental burden. On the other hand, the increase of forestry fixed assets would lead to the increase of tax, and the increase of tax can lead to the decrease of forestry efficiency. This conclusion was based on study of Ghebremichael and Potter-Witterverify [51], because they confirmed tax incentives did stimulate capital formation and total factor productivity growth. So, it was conducive to adjust the fixed assets of forestry in Northwest China. What’s more, IEC was the most important factor in random forest regression, so more attention should be paid on investment in forestry ecological construction. This paper results suggested the governments should promote the regional coordi- nated development of forestry industry. The results indicated that the ecological efficiency of forestry development was different in all regions, and the eastern and southern FEE was high, therefore, the advanced technology, experience and capital of the eastern region to the central and western regions were the focus of investment consideration. At the same time, the rich forestry resources in the northwest can be put into the forestry construction in the south, only in this way can it greatly promoted the optimization of the forestry industry in the whole region of China. Of course, the local government needed to innovate forestry science and technology means and improved the forestry management system. From the analysis of the factors of forestry ecological efficiency, EDL was a moderating important factor. However, in the study of Ester and Sebastian [52], GDP was a vital factor to the efficiency of the forest sector. The reasons for the different conclusions the 28 EU/EFTA countries taken as a subject, and GDP played an important role in interna- tional efficiency evaluation. It was also shown in the above results that ERL was other important factors affecting the ecological efficiency of forestry. The result was in accord with study of Chand et al. [53] that forest product benefits and environmental benefits were complementary to each other. Forestry enterprises need to strengthen the research on clean production technology, such as clean energy development, waste recycling, wastewater, and waste gas treatment technology in the production process, as governments need to im- prove the environmental management level of forestry industry and meet the requirements of forestry ecological security construction. In terms of the study of input and output efficiency between different regions, some scholars used conventional DEA, which likely lead to the situation that the calculated value of many units is 1, and the efficiency value of some units has little difference, so it is not better to compare the efficiency values of different regions. This result is the same as study of some scholars [32,54–61]. In contrast, the super efficiency SBM used in this paper can solve this problem well, which is also an advantage of efficiency comparison in this paper. However, there are also some shortcomings in this paper, the regional characteristics of FEE are very obvious, which is affected by regional climate, topography, water source, Forests 2021, 12, 142 15 of 17

and other factors, and these factors are not taken into account in this paper. In the later analysis of influencing factors of FEE, although reading a lot of relevant literature, some factors affecting FEE are not taken into account because the data are restricted, such as the implementation of forestry industry policies by local governments, the willingness of farmers for forest production, and the application of new forestry technology. These contents will be considered in the next step of forestry production efficiency research.

5. Conclusions Based on selected forests industry statistics from 31 Chinese provinces in China, the methods of super-efficiency SBM model and Malmquist productivity index to measure forestry efficiency has been applied in this paper from the static and dynamic angles. Then, to better study the distribution and evolution of FEE in China, this paper made a spatial analysis of the above ecological efficiency values of forestry industry by using ArcGIS software. Last, tobit regression and random forest were applied to explore the driving factors for FEE, and the conclusions were as follows. From dynamic angles, China’s forestry productivity change was in an upward trend from 2004 to 2018. The utmost average growth rate of TFPC from 2004 to 2018 is western regions, the TFPC of all Chinese regions from 2008 to 2009 decreased significantly. From static angles, Tianjin, Shanghai, and Hainan were the areas with high level of FEE and Hebei, Xinjiang, Gansu, and Shaanxi were the areas with high levels of FEE in these three periods. The result of spatial analysis showed that China’s FEE of the northeast and the south was basically higher, which FEE of China’s middle regions was relatively lower in 2004. However, as time went on, the distribution of FEE gradually appeared higher in the south and lower in the north. The H-H cluster appeared in south regions of China, these areas were economically developed and possessed high-tech industries. The regions of the L-L cluster area appeared the northeast and most of the northwest, which was consistent with the lower distribution FEE. Through the factor analysis of forestry ecological efficiency, there was a negative correlation between forestry fixed assets investment and FEE, but it was the most important factor to affect FEE. ERL and ISL were slightly important factors affecting FEE and they had a positive correlation with FEE.

Author Contributions: Conceptualization, S.C. and S.Y.; methodology, S.C.; software, S.C.; valida- tion, S.C. and S.Y.; formal analysis, S.C. and S.Y.; investigation, S.C.; sources, S.C. and S.Y.; data curation, S.C.; writing—original draft preparation, S.C.; writing—review and editing, S.C. and S.Y.; visualization, S.C.; supervision, S.Y. project administration, S.Y. funding acquisition, S.Y. All authors have read and agreed to the published version of the manuscript. Funding: This research was funded by the National Natural Science Foundation of China, grant number 71473195, 71773091; Special Fund for Scientific Research of Forestry Commonwealth Industry, grant number 201504424. Institutional Review Board Statement: Not applicable. Informed Consent Statement: Not applicable. Data Availability Statement: The data presented in this study are available in article. Acknowledgments: I sincerely thank my tutor for the guidance of my thesis. I am also grateful to the students in Center of Resource economy and Environmental Management Research for their help to me. Conflicts of Interest: The authors declare no conflict of interest. Forests 2021, 12, 142 16 of 17

References 1. Chen, K.; Li, X.; Piao, H. Analysis on Forestry Eco-Efficiency and Ecological Productivity Measure in China. For. Econ. 2016, 36, 115–120. [CrossRef] 2. Zhang, J.; Yang, J.; Cui, H. Research on the regional ecological efficiency evaluation based on the network DEA model in China. Soft Sci. 2016, 30, 15–19. [CrossRef] 3. Wang, R. The Influence of Environmental Regulation on the Efficiency of China’s Regional Green Economy Based on the GMM Model. Pol. J. Environ. Stud. 2020, 29, 2395–2402. [CrossRef] 4. Lin, B.; Ge, J. Carbon sinks and output of China’s forestry sector: An ecological economic development perspective. Sci. Total. Environ. 2019, 655, 1169–1180. [CrossRef][PubMed] 5. Chen, N.; Qin, F.; Zhai, Y.; Cao, H.; Zhang, R.; Cao, F. Evaluation of coordinated development of forestry management efficiency and forest ecological security: A spatiotemporal empirical study based on China’s provinces. J. Clean. Prod. 2020, 260, 121042. [CrossRef] 6. Chen, W.; Si, W.; Chen, Z.-M. How technological innovations affect urban eco-efficiency in China: A prefecture-level panel data analysis. J. Clean. Prod. 2020, 270, 122479. [CrossRef] 7. Ren, Y.; Fang, C.; Li, G. Spatiotemporal characteristics and influential factors of eco-efficiency in Chinese prefecture-level cities: A spatial panel econometric analysis. J. Clean. Prod. 2020, 260, 120787. [CrossRef] 8. Zheng, Y.; Gao, C.; Lei, G. Empirical Analyses of Forestry Industry Agglomeration and Eco-Efficiency:Based on Panel Data Test on 15 Provinces in China. Econ. Geogr. 2017, 10, 136–142. [CrossRef] 9. Xiong, L.; Wang, F.; Cheng, B.; Yu, C. Identifying factors influencing the forestry production efficiency in Northwest China. Resour. Conserv. Recycl. 2018, 130, 12–19. [CrossRef] 10. Yang, H.; Yuan, T.; Zhang, X.; Li, S. A Decade Trend of Total Factor Productivity of Key State-Owned Forestry Enterprises in China. Forests 2016, 7, 97. [CrossRef] 11. Viitala, E.; Hanninen, H. Measuring the efficiency of public forestry organizations. For. Sci. 1998, 44, 298–307. [CrossRef] 12. Sporocic, M.; Martinic, I.; Landekic, M.; Lovric, M. Measuring Efficiency of Organizational Units in Forestry by Nonparametric Model. Croat. J. For. Eng. 2009, 30, 1–13. [CrossRef] 13. Sowlati, T. Efficiency Studies in Forestry using Data Envelopment Analysis. For. Prod. J. 2005, 55, 49–57. [CrossRef] 14. Toma, P.; Miglietta, P.P.; Morrone, D.; Porrini, D. Environmental risks and efficiency performances: The vulnerability of Italian forestry firms. Corp. Soc. Responsib. Environ. Manag. 2020, 27, 2793–2803. [CrossRef] 15. Li, L.; Hao, T.; Chi, T. Evaluation on China’s forestry resources efficiency based on big data. J. Clean. Prod. 2017, 142, 513–523. [CrossRef] 16. Katrin, S.E.; Francisco, J.A.; Zuray, M. A three-stage DEA-SFA efficiency analysis of labour-owned and mercantile firms. J. Ind. Manag. Optim. 2017, 7, 573–592. [CrossRef] 17. Mcdonald, J. Using least squares and tobit in second stage DEA efficiency analyses. Eur. J. Oper. Res. 2009, 197, 792–798. [CrossRef] 18. Upadhyay, T.P.; Shahi, C.; Leitch, M.; Pulkki, R. An application of data envelopment analysis to investigate the efficiency of lumber industry in northwestern Ontario, Canada. J. For. Res. 2012, 23, 675–684. [CrossRef] 19. Miroslav, K. Efficiency of the Slovak forestry in comparison to other European countries: An application of Data Envelopment Analysis. Cent. Eur. For. J. 2018, 64, 46–54. [CrossRef] 20. Zadmirzaei, M.; Limaei, S.M.; Olsson, L.; Amirteimoori, A. Assessing the impact of the external non-discretionary factor on the performance of forest management units using DEA approach. J. For. Res. 2017, 22, 144–152. [CrossRef] 21. Zhang, Y.; Xiong, X. An empirical analysis of ecological efficiency evaluation and influencing factors of forestry in China under the background of green development: Based on DEA analysis perspective. J. Cent. South Univ. For. Technol. 2020, 40, 149–158. [CrossRef] 22. Malmquist, S. Index numbers and indifference surfaces. Trab. Estada 1953, 4, 209–242. [CrossRef] 23. Färe, R.; Grosskopf, S.; Lovell, C.A.K.; Grifelltatje, E. Biased Technical Change and the Malmquist Productivity Index. Scand. J. Econ. 1997, 99, 119–127. [CrossRef] 24. Färe, R.; Grosskopf, S.; Yaisawarng, S.; Li, S.; Wang, Z. Productivity growth in Illinois electric utilities. Resour. Energy 1990, 12, 383–398. [CrossRef] 25. Favero, C.A.; Papi, L. Technical efficiency and scale efficiency in the Italian banking sector: A non-parametric approach. Appl. Econ. 1995, 27, 385–395. [CrossRef] 26. Tone, K. A slacks based-measure of efficiency in data development analysis. Eur. J. Oper. Res. 2001, 130, 498–509. 27. Golany, B.; Roll, Y. An application procedure for DEA. Omega 1989, 17, 237–250. [CrossRef] 28. Froese, D.S.D.; Dobson, C.; White, A.; Wu, X.; Padovani, D.; Banerjee, R.; Haller, T.; Gerlt, J.; Surette, M.; Gravel, R. Sleeping beauty mutase (sbm) is expressed and interacts with ygfd in Escherichia coli. Microbiol. Res. 2009, 164, 1–8. [CrossRef] 29. Huang, H.; Wang, T. The Total-Factor Energy Efficiency of Regions in China: Based on Three-Stage SBM Model. Sustainability 2017, 9, 1664. [CrossRef] 30. Andersen, P.; Petersen, N.C. A Procedure for Ranking Efficient Units in Data Envelopment Analysis. Manag. Sci. 1993, 39, 1261–1264. [CrossRef] Forests 2021, 12, 142 17 of 17

31. Li, H.; Shi, J.-F. Energy efficiency analysis on Chinese industrial sectors: An improved Super-SBM model with undesirable outputs. J. Clean. Prod. 2014, 65, 97–107. [CrossRef] 32. Tone, K. A slacks-based measure of super-efficiency in data envelopment analysis. Eur. J. Oper. Res. 2002, 143, 32–41. [CrossRef] 33. Jong, P.; A Sprenger, C.J.; Veen, F. On Extreme Values of Moran’s I and Geary’s c. Geogr. Anal. 2010, 16, 17–24. [CrossRef] 34. Tiefelsdorf, M.; Boots, B. The Exact Distribution of Moran’s I. Environ. Plan. Econ. Space 1995, 27, 985–999. [CrossRef] 35. Lin, G.; Zhang, T. Loglinear Residual Tests of Moran’s I Autocorrelation and their Applications to Kentucky Breast Cancer Data. Geogr. Anal. 2007, 39, 293–310. [CrossRef] 36. Tiefelsdorf, M.; Boots, B. A Note on the Extremities of Local Moran’s Iis and Their Impact on Global Moran’s I. Geogr. Anal. 2010, 29, 248–257. [CrossRef] 37. Bivand, R.; Müller, W.G.; Reder, M. Power calculations for global and local Moran’s I. Comput. Stat. Data Anal. 2008, 53, 2859–2872. [CrossRef] 38. Liang, K.-Y.; Zeger, S.L. Longitudinal data analysis using generalized linear models. Biometrics 1986, 73, 13–22. [CrossRef] 39. Wooldridge, J.M. Simple solutions to the initial conditions problem in dynamic, nonlinear panel data models with unobserved heterogeneity. J. Appl. Econ. 2005, 20, 39–54. [CrossRef] 40. Gordon, A.D.; Breiman, L.; Friedman, J.H.; Olshen, R.A.; Stone, C.J. Classification and Regression Trees. Biometrics 1984, 40, 874. [CrossRef] 41. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [CrossRef] 42. Liaw, A.; Wiener, M.; Liaw, A. Classification and Regression with Random Forest. R News 2002, 23, 18–22. 43. Zarei, A.R.; Mahmoudi, M.R.; Shabani, A. Investigating of the climatic parameters effectiveness rate on barley water requirement using the random forest algorithm, Bayesian multiple linear regression and cross-correlation function. Paddy Water Environ. 2021, 19, 137–148. [CrossRef] 44. Islam, S.; Amin, S.H. Prediction of probable backorder scenarios in the supply chain using Distributed Random Forest and Gradient Boosting Machine learning techniques. J. Big Data 2020, 7, 1–22. [CrossRef] 45. Svetnik, V.; Liaw, A.; Tong, C.; Culberson, J.C.; Sheridan, R.P.; Feuston, B.P. Random Forest: A Classification and Regression Tool for Compound Classification and QSAR Modeling. J. Chem. Inf. Comput. Sci. 2003, 43, 1947–1958. [CrossRef] 46. Manilich, E.; Zsoyo, L.; Meral, Z.; Trubachev, V.; Radivoyevitch, T. Classification of large microarray datasets using fast random forest construction. J. Bioinform. Comput. Biol. 2011, 9, 251–267. [CrossRef] 47. Seo, D.K.; Kim, Y.H.; Eo, Y.D.; Park, W.Y. Learning-Based Colorization of Grayscale Aerial Images Using Random Forest Regression. Appl. Sci. 2018, 8, 1269. [CrossRef] 48. Zhang, L. Research on Importance of Influencing Factors of Resident Outbound Tourism in China Based on Random Forest Model. Resour. Dev. Mark. 2017, 33, 711–716. [CrossRef] 49. Liu, T.; Li, J. Temporal-Spatial Differentiation and Its Influencing Factors of Green Total Factor Productivity of Forestry in China. World For. Res. 2020, 33, 56–61. [CrossRef] 50. Luo, X.; Li, Z.; Li, R.; Xue, L. Temporal and regional variation of forestry production effienciency in China. J. Arid Land Resour. Environ. 2017, 31, 95–100. [CrossRef] 51. Ghebremichael, A.; Potter-Witter, K. Effects of tax incentives on long-run capital formation and total factor productivity growth in the Canadian sawmilling industry. For. Policy Econ. 2009, 11, 85–94. [CrossRef] 52. Ester, G.; Sebastián, L. Cross-country comparison of the efficiency of the European forest sector and second stage DEA approach. Ann. Oper. Res. 2020.[CrossRef] 53. Chand, N.; Kerr, G.N.; Bigsby, H. Production efficiency of forest management in Nepal. For. Policy Econ. 2015, 50, 172–179. [CrossRef] 54. Jha, R.; Gundimeda, H.; Andugula, P. Assessing the Social Vulnerability to Floods in India: An Application of Superefficiency Data Envelopment Analysis and Spatial Autocorrelation to Analyze Bihar Floods. Econ. Eff. Nat. Disasters 2021, 559–581. [CrossRef] 55. Aydin, U.; Karadayi, M.; Ulengin, F. How efficient airways act as role models and in what dimensions? A superefficiency DEA model enhanced by social network analysis. J. Air Transp. Manag. 2020, 82, 101725. [CrossRef] 56. Park, C. A Study on the Efficiency of Fishing-Ports Based on Super-SBM. J. Fish. Bus. Admin. 2010, 41, 129–151. 57. Gokgoz, F.; Erkul, E. Investigating the energy efficiencies of European countries with super efficiency model and super SBM approaches. Energy Effic. 2019, 12, 601–618. [CrossRef] 58. Ma, F.; Wang, W.; Sun, Q.; Liu, F.; Li, X. Integrated Transport Efficiency and Its Spatial Convergence in China’s Provinces: A Super-SBM DEA Model Considering Undesirable Outputs. Appl. Sci. 2018, 8, 1698. [CrossRef] 59. Azadi, M.; Saen, R. Developing an Output-Oriented Super Slacks-Based Measure Model with an Application to Third-Party Reverse Logistics Providers. J. Multi-Criteria Decis. Anal. 2012, 18, 267–277. [CrossRef] 60. Zimková, E. Technical Efficiency and Super-efficiency of the Banking Sector in Slovakia. Procedia Econ. Financ. 2014, 12, 780–787. [CrossRef] 61. Zanboori, E.; Rostamy-Malkhalifeh, M.; Jahanshahloo, G.; Shoja, N. Calculating Super Efficiency of DMUs for Ranking Units in Data Envelopment Analysis Based on SBM Model. Sci. World J. 2014, 2014, 1–7. [CrossRef][PubMed]