Changes in ecosystem services and an analysis of driving factors for China's Natural Forest Conservation Program

Abstract China's Natural Forest Conservation Program (NFCP) is aimed at improving the fragile and unstable ecological environment and has become one of the largest ecological restoration programs in the world because of its enormous investment and effects. It is important to work out and strengthen new measures to overcome difficulties to promote more ecosystem services and human well‐being in the NFCP. This study focused on how to evaluate the ecosystem services change brought about by implementing the NFCP. Taking the key state‐owned forest areas in the Northeast and Inner Mongolia as the study area, we provide a basic overview of development and construct an evaluation index system and a distributed calculation method for the NFCP to analyze the implications of the NFCP on ecosystem services combined with multi‐source data coupling. An evaluation index system for NFCP ecosystem services was constructed. The system includes five ecological service functions and 12 evaluation indices. The trade‐off and/or synergistic analysis of ecosystem services were carried out. The regional characteristics and changes in the NFCP ecosystem services were emphasized. Although it has not been implemented for a long time, the NFCP has had a great impact on ecosystem services because it reduces soil and water losses, increases soil fertility, strengthens the forest carbon sink and helped the forest accumulate nutrients and purify the atmosphere. Socioeconomic factors affect the NFCP ecosystem services, such as the implementation area of NFCP, investment amount of NFCP, area ratio of nature reserves, and yield of tree stock volume. Policy drivers of the NFCP, changes in the economic structure and reductions in forest yield are the main factors affecting the change in NFCP ecosystem services. Although the NFCP has positively affected the society, the economy, and the ecological environment, it has also generated some problems, such as the improper management of forest resources, shortage of capital investment, staff transfer, etc. The social and economic problems will be transient with implementation of the NFCP, and the structural changes in forestry and agriculture may eventually benefit the forestry workers and other stakeholders.


| INTRODUC TI ON
Ecosystem services are constituted by the close association between environmental conditions and human well-being (Tallis, Goldman, Uhl, & Brosi, 2009). However, the increased social and economic activities of humans have led to degradation of the ecological environment, resulting in many challenges for global sustainable development (Mee, Dublin, & Eberhard, 2008), and nearly 60% of the degradation is unsustainable in ecosystem services (MEA (2005)). Human activities, particularly agricultural expansion and deforestation have covered half of the Earth's land area (Kareiva, Watts, McDonald, & Boucher, 2007; Vitousek, Mooney, Lubchenco, & Melillo, 1997). As forests have functions of soil and water conservation, carbon sequestration, and biodiversity conservation (Foley et al., 2005;Hansen et al., 2013), their degradation will become the most harmful of human activities (Crowther et al., 2015). Under such a challenging background, ecological restoration and protection have been widely recognized and implemented and are based on various forms of ecological restoration projects (Allan et al., 2013;Bullock, Aronson, Newton, Pywell, & Rey-Benayas, 2011). The China's Natural Forest Conservation Program (NFCP) was launched after China suffered a great flood in 1998 and was a huge project for protecting and restoring the forests. Natural forest resources are effectively protected to prohibit commercial felling of natural forests and forest coverage and ecosystem services have improved (NDPRC, 2009). In addition, the implementation of NFCP has a certain global significance, although it was initially cope with the declining ecological environment in China. However, if implemented reasonably and continuously, it will have an important impact on China and the rest of the world. This study focuses on how to quantitatively evaluate the implications of the ecosystem services.
In recent years, the research of natural forest resources protection projects has emerged in the world. The main research direction can be summed up in four aspects: the practical research on the conservation project of natural forest resources in China, the past experience and future focus of natural forest conservation projects in the world, study on the relationship between conservation of natural forest resources and biodiversity and ecological protection in developing countries (Chen, Cao, & Su, 2017). In the past decades, there has been no interruption in the research of natural forest resources protection program. With the continuous deepening of people's understanding of it, the focus of its research is also changing. In the years after China implemented the natural forest protection project in 1998, due to its large scale, long span and huge investment, it has become the focus of academic research. The research in this period is mainly focused on the policy background, afforestation measures, management methods, and effects of natural forest resources protection projects in China (Lewis, 2005;Lu & Zhang, 2003). As the research goes deep, relevant scholars have studied the protection of natural forest resources for biodiversity conservation and effective land use, and gradually turned from macro to microcosmic (Hannon & McCallum, 2006;Khanina, Bobrovsky, Komarov, & Mikhajlov, 2007). Then, people realize that the important role of this program in carbon fixation has become a hot spot of research (Ryan et al., 2010;Soares-Filho et al., 2010). Later, the ecosystem service value of this project and the restoration of forest vegetation have become the focus of attention (Santopuoli, 2012;Zhou, Svensson, Yan, Chen, & Li, 2014). Consequently, more people are paying close attention to how to quantitatively assess the impact of the NFCP on ecosystem services, how implementing the NFCP has changed ecosystem services, and the key drivers of the impact (Yu & Han, 2016). But for many reasons, such as lack of unified evaluation index system, complexity of social economy, and lack of reference materials, making such a quantitative assessment very challenging.
National key ecological restoration projects have had an important impact on China's degraded ecosystems, and are the focus of current research. Studies have shown that terrestrial ecosystems in China are important carbon sinks, and that any increase in carbon storage is largely attributed to climate change, ecological restoration projects, and agricultural land management (Fang, Yu, Liu, & Chapin, 2018). For example, ecological restoration projects have contributed to an increase in carbon storage of 74 Tg C/year (Lu et al., 2018). From 1978 to 2015, China has invested more than USD 370 billion across 6.2 million km 2 (65% of the country's total area) of land through 16 major ecological restoration programs. More than 500 million people have participated in them. Their implementation has resulted in very positive impacts: for example, forests now cover over 22% of China's land area, grasslands have expanded and regenerated, and soil erosion and sedimentation have been significantly reduced. The sediment load of the Yellow River has dropped by 90%, agricultural productivity has grown by 5% annually on average, and millions of farmers have been lifted out of poverty . The experiences and achievements gained by many major ecological restoration programs in China help improve future sustainable development projects and policies, as well as to provide examples and lessons for sustainable development in other countries .
In order to reveal the impact of ecological restoration projects on regional ecosystem services, further strengthening of regional ecological management will be generated based on the Northeast and Inner Mongolia key state-owned forest region as the study area and the status quo of the NFCP, according to the People's Republic of China Forestry Standard: Forest Ecosystem Service Valuation Norms (LY/ T1721-2008), which was issued by the State Forestry Administration, P.R. China in 2008(SFAC, 2008. The following research objectives were addressed in this study: (a) analyze the changes in ecosystem services of the NFCP from multivariate data; (b) demonstrate the key drivers for change in NFCP ecosystem services; and (c) discuss the impact of drivers on ecosystem service functions, to clarify the existing problems and challenges, and to provide suggestions for overcoming the difficulties and enhancing their potential. The results will provide a reference for the development, management, and evaluation of similar ecological restoration projects in China and other countries worldwide.

| Study area
The key state-owned forest areas in the Northeast and Inner Mongolia are mainly distributed in Heilongjiang, Jilin and east of the Inner Mongolia Autonomous Region Province (Figure 1), located at 117°06′-135°05′E, 41°25′-53°23′N. The main geomorphic types included mountains, hills, plains, and plateaus, from 30 to 1,600 m.
The region is located in the eastern part of the Eurasian continent, most of which belongs to the central temperate zone. Only a few parts of the north belong to the cold temperate zone. The soil types are dominated by dark brown soil, black soil, and brown forest soil.
The NFCP covers the key state-owned forest areas in Northeast China and Inner Mongolia provinces with a total area of 3,589.97 × 10 4 ha, accounting for 3.74% of the national land area. Since implementation of the NFCP, forest resources have received effective protection and restoration to the NFCP requirements, and commercial logging has stopped; therefore, the forests area and volume have increased from 2,673.42 × ha10 4 ha to 2,788.72 × ha10 4 ha and 18.96 × 10 8 m 3 to 23.30 × 10 8 m 3 , respectively (Table 1).

| Ecological function monitoring and layout of the assessment area
The main indicators of ecological function monitoring and the assessment area according to the differences in natural conditions of the natural forest protection projects in key state-owned forest areas F I G U R E 1 Layout of the forest ecological stations belongs to CFERN of the ecosystem services and the ecological functions monitoring and the assessment area of the NFCP in the key state-owned forest areas in Northeast and the eastern of Inner Mongolia. CFERN: Chinese Forest Ecosystem Research Station Network. NFCP: Natural Forest Conservation Program. The following the same. I-A-a-1: Micro wind erosion area of deciduous coniferous forest in cold temperature zone; I-A-a-2: Non erosive of deciduous coniferous forest in cold temperature zone; I-A-a-6: Intensity water erosion area of deciduous coniferous forest in cold temperature zone; II-A-b-2: No erosion area in humid temperate coniferous and broad-leaved mixed forest; II-A-b-3: Micro water erosion area of coniferous and broad-leaved mixed forest in middle temperate; II-A-b-4: Mild water erosion area in humid temperate coniferous and broad-leaved mixed forest; II-A-b-5: Moderate water erosion area in coniferous and broad-leaved mixed forest of moderate temperate zone; II-A-B-6: Intensity of water erosion area in humid temperate coniferous and broad-leaved mixed forest; II-A-c-6: Moderate humid grassland in moderate intensity water erosion area; II-B-c-1: Micro wind erosion region in mid temperate semi humid grassland; II-B-c-2: No erosion area in mid temperate semi humid grassland; II-B-c-3: Micro water erosion region in mid temperate semi humid grassland. The following the same in Northeast and Inner Mongolia are temperature and water (Zheng, 2008), an index of vegetation regionalization (Zhang, 2007), and soil erosion indices (Li, Zhao, & Yang, 2014). Geometric rectification and vectorization were carried out by defining the projection uniformly in ArcGIS10.0, and then the NFCP implementation scope layer, different temperature and moisture distribution layer regions of China, Chinese vegetation distribution layer, and China soil erosion distribution layer were obtained. Finally, the corresponding layers of Northeast China and Inner Mongolia area were defined, and the ecological function monitoring and evaluation regionalization of the study area was generated ( Figure 1).

| Forest ecology observation station
Some of the current forest ecosystem stations have regarded NFCP ecosystem services as one of the main monitoring objectives (Wang, Cui, & Yang, 2004). The network layout for the forest ecology station was guided by typical sampling, based on the water and heat distribution and the conditions at the forest site and typical, representative areas were selected to complete the forest ecological network layout (Guo, Wang, & Niu, 2015). Figure 1 shows the forest ecology station responsible for monitoring the ecosystem services of the NFCP in the key state-owned forest areas in Northeast China and eastern Inner Mongolia. At least 10% of the counties or forestry bureaus in the NFCP area are under long-term observations to ensure observational accuracy and collect sufficient observational data. There are eight forest industry groups, 41 counties/cities, 100 forestry bureaus, and 18 forest ecology stations (Niu, Wang, & Wei, 2013) in the study area.

| Data sources
The data sources for the changes in physical measurements of the  Tables 2 and 3). Finally, the quality and quantity of the NFCP ecosystem services in the study area were evaluated.  the forest ecosystem (Wang et al., 2003). The index system included five ecological service functions and 12 evaluation indices ( Figure 2). The five ecological service functions are water conservation (water quantity regulation and water purification), soil conservation (decrease soil loss and protection of soil fertility), carbon fixation and oxygen release, nutrients accumulation in the forest, atmospheric environmental purification (negative ion supply, pollutant absorption, total suspended particulate(TSP) absorption, PM 10 absorption, and PM 2.5 absorption; Wang et al., 2011). Table 2 shows the model, formulas, and definition of the 12 evaluation indicators for NFCP ecosystem services (Wang et al., 2003).

| Trade-off and/or synergies analysis of ecosystem services for the NFCP
Zoning statistics of ecosystem services were carried out based on the unit of the forestry bureau in the NFCP area. The values were standardized, and the changes over the whole research period (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015) were calculated. Finally, the changes were imported into SPSS for bivariate correlation analysis, and correlation coefficients between water conservation, soil conservation, carbon sequestration and oxygen release, tree nutrient accumulation, and air purification were obtained (Table 6).
TA B L E 2 Models of the ecosystem services evaluation of NFCP

Evaluation index Formula Parameter descriptions
Water quantity regulation and purification G WAC = 10A·(P-E-C)·F G WAC is the quantity of water regulated (or water quality purification) for NFCP per year (m 3 /a); P is precipitation in open land (mm/a); E is evapotranspiration of trees and other plants in the stand (mm/a); C is surface runoff (mm/a); and A is the area of NFCP (ha); F is forest ecological function correction coefficients (FEFCC).The following is the same Decrease of soil loss G DSL = A·(X 2 −X 1 )·F G DSL is the quantity of fixed soil per year (t/a); X 1 and X 2 are the modulus of soil erosion of forestlands and non-forestlands for NFCP (t [hm −2 a −1 ]).The following the same Protection of soil fertility G PSF = A·(N + P + K + M)·(X 2 −X 1 )·F G PSF is the quantity of fertility maintained for NFCP per year (t/a); N is the concentration of nitrogen in forest soil of NFCP (%); P is the concentration of phosphorus in forest soil of NFCP (%); K is the concentration of potassium in forest soil of NFCP (%); M is the concentration of organic matter in forest soil of NFCP (%) Carbon fixation G CAF = A·(1.63·R c ·B a + F sc )·F G CAF is the quantity of carbon fixation for NFCP per year (t/a); R c is the carbon content in carbon dioxide, it is 27.2%; B a is the net primary productivity of the forest for NFCP (t·(hm −2 ·a −1 )); F sc is the carbon sequestrated in the soil per unit area for NFCP per year(t·(hm −2 ·a −1 )). The following the same Oxygen released G OXR = 1.19·A·B a ·F G OXR is the quantity of oxygen released for NFCP per year (t/a).

Nutrients accumulation
G NUA is the quantity of N, P, and K accumulated in a stand type of NFCP (t/a); N f , P f and K f are the concentrations of N, P, and K in dry matter (%) Negative ion supply G NIS = 5.256.10 15 ·Q a ·A·H·F/L G NIS is the number of negative ion produced for NFCP per year (N/a); Q a is the concentration of anions in the atmosphere in a stand type (ge/cm 3 ); H is the height of trees for NFCP (m); L is the lifetime of a negative ion in the air for NFCP (min) G POA is the quantity of pollutants absorbed for NFCP per year (t/a); Q SO2 is the quantity of sulfur dioxide absorbed for NFCP per unit area per year (kg·(hm −2 a −1 )); Q F is the quantity of fluoride absorbed for NFCP per unit area per year (kg·(hm −2 a −1 )); Q NOx is the quantity of nitrogen oxides absorbed for NFCP per unit area per year (kg·[hm −2 a −1 ]) TSP absorption G TSP = Q d ·A·F/1,000 G TSP is the quantity of TSP absorbed for NFCP per year (t/a); Q d is the quantity of dust absorbed for NFCP per unit area per year (kg·[hm −2 a −1 ]). The following the same PM 10 absorption G PM10 = 10·Q PM10 ·A·n·F·LAI G PM10 is the quantity of PM 10 absorbed for NFCP per year (t/a); Q PM10 is the measured absorption of PM 10 in the unit area of the stand (kg/a); LAI is leaf area index. The following the same PM 2.5 absorption G PM2.5 = 10·Q PM2.5 ·A·n·F·LAI G PM2.5 is the quantity of PM 2.5 absorbed for NFCP per year (t/a); Q PM2.5 is the measured absorption of PM 2.5 in the unit area of the stand (kg/a) Notes. CAF: Carbon fixation, it is the sum of plant carbon fixation and soil carbon fixation; DSL: Decrease of soil loss; NUA: Nutrients accumulation, it is the sum of accumulation of nitrogen, phosphorus, potassium; NIS: Negative ion supply; OXR: Oxygen release; POA: Pollutants absorption, it is the sum of absorption of sulfur dioxide, fluoride and nitrogen oxides; PM 10 : PM 10 absorption; PM 2.5 : PM 2.5 absorption; PSF: Protection of soil fertility, it is the sum of soil nitrogen, phosphorus, potassium and organic matter content; TSP: Total suspended particulate absorption; WAC: Water conservation. The following the same.

| Analysis of factors influencing NFCP ecosystem services based on the DPSIR-mDSS model
Driving force-Pressure-State-Impact-Response (DPSIR) is a conceptual framework used to reveal the causal relationship between the state of the environment and human activities, and has been widely used in the study of resource management, population, and sustainable development (Benini, Bandini, Marazza, & Contin, 2010;Burkhard & Müller, 2008;Fassio, Giupponi, Hiederer, & Simota, 2005;Malekmohammadi & Jahanishakib, 2017). Because of the different evaluation objects and emphasis, the DPSIR conceptual framework often has great differences in the evaluation indexes and methods of the ecological environment response assessment, but its core idea is that natural forces and human activities (driving forces) constitute a huge "pressure" on the ecosystem, shifting the ecological environment into a deteriorated "state". This state produces an unsustainable "impact," which human beings must take certain measures (responses) to reduce and thus achieve sustainable development (EEA, 2005;OECD, 1991). In a project of sustainable use of water resources at the river basin scale, the European Union tried to combine the multi-project integrated decision support system (mDSS) with the DPSIR conceptual framework, and achieved good results. issues have been innovatively addressed. Therefore, the combi- (1) The objective weight of the entropy method is as follows: where W j is the weight of the j index, k is the index number, g j is the j index difference coefficient, e j is the j index entropy, R ij is the pro- where Q j and Q i represent the composite values of the j and i schemes, q jk and q ik represent the composite values of the index k in the j and i schemes, and w k is the weight of index k.

| Evaluation of ecosystem services during implementation of the NFCP
The evaluation of ecosystem services during implementation of the NFCP in the key state-owned forest areas in Northeast China and Inner Mongolia is shown in Table 5. The differences in the ecosystem services between 2000 and 2015 were mainly manifested in the following aspects.
Soil and water losses were reduced. As shown in Table 5, after implementing the NFCP, the quantity of water conservation (WAC) in the study area increased by 157.57 × 10 8 m 3 /a, which was equivalent to 40.09% of the maximum capacity of the Three Gorges reservoir (393 × 10 8 m 3 ) and a 30.21% increase compared with prior to the NFCP.
The decrease in soil loss (DSL) increased by 2.38 t/a, equivalent to 44.65% of the total amount of soil erosion (5.34 × 10 8 t, Bulletin of Soil and Water Conservation in China, 2014; MWRPRC, 2014) of China's 11 major rivers in 2014 and 21.14% higher than before the NFCP. Table 5, protection of soil fertility (PSF) in the study area increased by 1,747.78 × 10 4 t/a, equivalent to 29.02% of China's fertilizer application rate (6,022.6 × 10 4 t,
Enhanced forest carbon sequestration and accumulated nutrients. As shown in Table 5, after implementing the NFCP, the quantity of carbon sequestration (CAS) in the study area increased by 1,573.03 × 10 4 t/a, equivalent to 0.52% of the total carbon released by China's standard coal consumption (43 × 10 8 t, China Statistical Yearbook, 2016;NBSPRC, 2016) in 2015, and 26.98% higher than before the NFCP. Nutrients accumulation (NUA) increased by 127.83 × 10 4 t/a, which was 39.99% higher than before the NFCP.
Purify the atmospheric environment. As shown in Table 5

| Trade-off and/or synergies of ecosystem services for the NFCP
Correlation analysis showed that there were strong trade-offs between water conservation and soil conservation, water conservation and carbon fixation and oxygen release, water conservation and nutrients accumulation, soil conservation and carbon fixation and oxygen release, soil conservation and nutrients accumulation (the highest correlation coefficient was 0.912). However, there is a weak balance between atmosphere environmental purification and the four functions of water conservation, soil conservation, carbon fixation and oxygen release, and nutrients accumulation (the maximum absolute value of correlation coefficient is only 0.226) (Table 6). Specifically, water conservation and soil conservation, water conservation and carbon fixation and oxygen release, water conservation and nutrients accumulation are positively correlated, showing a mutually beneficial synergistic relationship (correlation coefficients were 0.873, 0.682, and 0.578, respectively). Water is one of the key elements in structuring ecosystems and controlling the flow of ecological processes and functions (Fu, Wang, Su, & Forsius, 2013), stable water supply is the basis for maintaining many ecosystem services (Maes et al., 2012). There were positive correlations between soil conservation and carbon fixation and oxygen release, and between soil conservation and nutrients accumulation (correlation coefficients were 0.842 and 0.912, respectively), which showed that stable soil conservation provided a guarantee for carbon fixation of ecosystem and nutrients accumulation. Atmosphere environmental purification is negatively correlated with water conservation, soil conservation, carbon fixation and oxygen release, and nutrients accumulation, which is a trade-off relationship between growth and decline (correlation coefficients were 0.105, 0.145, 0.094, and 0.226, respectively).

| Key factors affecting the change in ecosystem services in the NFCP
Some studies have shown that the socioeconomic (policy, population, economic development) and ecological environments (climate, soil, vegetation) were the two major factors affecting ecosystem services during ecological restoration. This is somewhat similar to the results of this study (Gao, Niu, Wang, & Zheng, 2015;Liu, Li, Ouyang, Tam, & Chen, 2008). A sensitivity analysis (Figure 4) showed that the lateral axis represents the difference between the composite values, and the current situation represents the composite value difference in the contrast scheme. The more comprehensive value difference was less than zero, indicating that the stronger the ability of this factor to change the comprehensive value, the higher the sensitivity.
The weights of one factor were changed sequentially, and the weights of other factors remained unchanged. The sensitivity of TA B L E 5 Annual ecosystem services of different ecological function monitoring and assessment areas of the implementation period for NFCP in key state-owned forest areas in Northeast and Inner Mongolia

| Comments on the DPSIR-mDSS model
The DPSIR-mDSS model uses the DPSIR conceptual model to construct the ecosystem service evaluation index system for ecological restoration to explore the causes of ecosystem service status (state), direct factors (pressure), indirect factors (driving force), and ecosystem service status (state) caused by the changes in ecosystem service function (impact), and the relevant measures for the adverse effects (response; Giupponi et al., 2004). In this process, the complex relationship between human behavior and the ecological system becomes a simple cause and effect relationship, making the evaluation process clearer and easier to understand. Because it is a causal relationship, a quantitative index for the interaction between them is possible (Durairaj, Sudharsun, & Swamynathan, 2013).
Notably, ecological restoration needs to go through the process within a certain time interval (Celentano & Rousseau, 2016

| Problems and recommendations for the NFCP
The sensitivity analysis based on the model showed that policy and capital input for the NFCP played a pivotal role in ecosystem service changes. The NFCP has had positive impacts on the social, economic, and ecological environments but also faces some problems (Ren et al., 2015;Xu, Yin, Li, & Liu, 2006). The accuracy of ecosystem services valuation is closely related to ecosystem services monitoring, ecosystem services management can provide support for relevant policies, and wages and population density vary with staff turnover. As a factor that influences the change in NFCP ecosystem services, forest yield has a direct relationship with forest resource management and capital investment.
Therefore, this study combined factors that impacted ecosystem services change for the NFCP to reveal the problems with the NFCP from the following aspects and puts forward some suggestions.

| Special monitoring of ecosystem services
Currently, there are obvious deficiencies in the basic data for the evaluation of the NFCP ecosystem services (Ke, Chen, Robson, Wen, & Tian, 2015;Yin et al., 2010). First, the forest resources data of the basic implementation units need to be determined in the NFCP area. In terms of current technical level and database, it is feasible and necessary to make use of the spatial analysis technology to analyze the shared high resolution remote sensing image data, and determine the NFCP area and the implementation of unit space boundary (Lü et al., 2015;Viña, McConnell, Yang, Xu, & Liu, 2016). Second, no NFCP ecological monitoring station has been established. In this study, only the CFERN and auxiliary observation sites in nearby areas or similar habitats were replaced, which reduces evaluation accuracy. Therefore, the next step is to establish a national ecological monitoring station for the natural forest protection project area, so that the data can be used. In the future, more and more measured data can be obtained with increases in the number of monitoring stations and determining the layout of natural forest protection project.

| Ecosystem services management
A comprehensive and systematic simulation and estimation of ecosystem services is a prerequisite for ecosystem services management and related policy formulation. The spatial heterogeneity of geographical environment will lead to uncertainty in the trade-off and coordination of ecosystem services. This paper concludes that the trade-off and synergy between various ecosystem services in the study area is basically equivalent, and the synergies are more significant, which also provides a good way for the region to formulate ecological restoration and protection policies. This indicates F I G U R E 4 Results of sensitive analysis for the NFCP that the implementation of the NFCP has played a certain effect. In addition, there is a staggered form of trade-offs and/or synergies between various ecosystem services, and its internal mechanism and mechanism of action need further exploration and in-depth analysis in the future.

| Forest resource management
Forestry limited to commercial forest areas and limited cutting areas for forest annual tending is not conducive to the protection of forest resources. However, a large number of forests in the forbidden cutting area need to be tended and felled sanitarily.
Without tending and sanitary felling, the growth of trees is poor, and some of the dead trees are attacked by pests and diseases. In

| Staff shunt transfer
The biggest problem facing enterprises stopped after cutting is the reorienting workers (Ke et al., 2015;Yin et al., 2010). Therefore, forest industry enterprises should actively encourage employees through forest tending, management and development of the forest economy, multi-channel development of entrepreneurial opportunities and jobs (Huang, Liu, & Tian, 2014;Zang, Zhang, & Hao, 2014).
The first is to promote the sustainable management of forestry, strengthen the cultivation of reserve resources, and enhance the ability to create employment and absorb extra workers. The second is a multi-channel placement of surplus staff and workers (Dawley, Marshall, Pike, Pollard, & Tomaney, 2014;MacKinnon, 2017), promote the reorientation of workers shunt. Third is to encourage independent entrepreneurship, entrepreneurship to promote employment, the establishment of entrepreneurship and employment "two wheeled" drive mechanism (Baumgartner, Pütz, & Seidl, 2013;Parker Harris, Renko, & Caldwell, 2014).

ACK N OWLED G M ENTS
The study was subsidized by the Integration and application dem-

CO N FLI C T O F I NTE R E S T
None declared.

AUTH O R CO NTR I B UTI O N S
Bing Wang designed the study. Xiang Niu and Peng Gao performed the analysis. Qingfeng Song helped collect the data. All authors contributed to the interpretation of the results. Longsheng Huang wrote the initial draft and refined the manuscript. Special thanks are given to the referees and the editors for their instructive comments, suggestions, and editing for the manuscript.

DATA ACCE SS I B I LIT Y
Upon article acceptance, original data including evaluation results of ecosystem services and social public will be archived in the public repository Dryad.