Incorporating carbon sink of harvested wood products into ecological footprint accounting: model and case study

ABSTRACT Ecological footprint (EF) accounting system is widely used in assessing environmental sustainability at various spatial scales. To improve the EF performance, a new EF calculation framework was established in this study, in which HWP is considered as carbon sinks and the model was then applied to the Yangtze River basin of China. Results showed that (1) carbon sink in the basin is about 2.99 Tg per year and it varies substantially across subbasins, which has a nonlinear relationship with net primary productivity; (2) available biocapacity calculated in the whole basin increases from 0.5068 to 0.5759 national ha (nha) per capita when incorporating the carbon sinks. In a few subbasins, the added biocapacity is even larger than the biocapacity estimated from the conventional approach; (3) the EF of the basin is 5.598 nha per capita and it is much larger than its biological carrying capacity (BC), indicating that ecological deficit is ubiquitous across the basin. If carbon sinks are considered, the ecological deficit would decrease by 3.2–14.3% for different regions in the basin. It is necessary to incorporate HWP and other potential terrestrial carbon sinks in the EF accounting, and more efforts are needed in improving carbon sink accuracy.


Introduction
With the rapid development of the global population and economy, human's impact on the planet is increasing considerably and has exceeded the limit of its carrying capacity in many aspects (Rockstrom et al. 2009). To mitigate deteriorating environmental sustainability, measuring anthropogenic influence on the natural environment is the fundamental step to manage both humans and nature (Wackernagel, Onisto, and Bello et al. 1999). Various environmental models have been suggested to quantify human's impact on the natural environment to assess environmental sustainability, such as net primary productivity (NPP) (Vitousek, Ehrlich, and Ehrlich et al. 1986), material flow analysis (Fischer-Kowalski and Haberl 1998), pressure-state-response model , emergy analysis (Bonilla et al., 2019), life cycle assessment (Castellani and Sala 2012), and ecological footprint (EF) (Wackernagel, Onisto, and Bello et al. 1999). Among them, EF is one of the most widely used approaches (Mostafa 2010) and has been proved to be an effective tool to assess the impact of anthropogenic activity on the natural ecosystem . The EF is an important model for measuring environmental sustainability, which was first introduced by William Rees (1992) and developed by Wackernagel, Onisto, and Bello et al. (1999), and was applied in many countries. The EF accounting system measures environmental sustainability by comparing the bio-productive land areas that are consumed and produced in a region (Venetoulis and Talberth 2008). As a composite model of quantifying anthropogenic intervention and natural resources, the EF approach is consisted of a lot of calculating items, parameters, and assumptions that may lead to uncertainty in EF calculation results (Venetoulis and Talberth 2008;Beynon and Munday 2008;Lenzen and Murray 2001). To decrease model uncertainty and improve its performance, various modified models of EF have been established by including new calculation items (Lenzen and Murray 2001) or presenting new methods of calculating parameters such as the equivalence factor and the yield factor (YF) (Barrahmoune, Lahboub, and El Ghmari 2019;Li, Zhang, and Chang et al. 2018). In the EF case study of Australia, a revised model which considered greenhouse gases other than CO 2 and regional productivity other than the world average was used in the calculation (Lenzen and Murray 2001). To improve parameter estimation in the model, NPP was suggested to estimate the YF and equivalent factor (EQF) (Liu, Fu, and Jiang et al. 2019). Some assumption in the model such as the calculation of YF and EQF for the built-up land has also been assessed and revised. In a recent study (Li, Zhang, and Chang et al. 2018), the assumption that built-up land is converted from cropland is modified from historic land-use change, concluding that in some urbanized regions, the assumption has a significant influence on the EF. New EF components such as water EF are included in the calculation in a recent study . The revised form of EF accounting system has improved the model performance and enhanced its ability to deal with different environmental and socioeconomic situations.
In the conventional EF accounting, biologically productive land used for assimilating carbon emission from fossil combustion (or carbon footprint) often accounts for the largest source of EF (Li, Zhang, and Chang et al. 2018;Yilanci, Gorus, and Aydin 2019;Wackernagel, Schulz, and Deumling et al. 2002;Lin et al. 2018). Therefore, accurate estimation of biological capacity for carbon assimilation is critical to estimate EF. The current EF accounting system assumes that biologically productive land is mutually exclusive (Wackernagel, Onisto, and Bello et al. 1999). Under this assumption, forest is only calculated as biological land use for producing forestry products mainly including roundwood and fuelwood, while its function for absorbing carbon dioxide is not calculated as biocapacity, and sequestration area is set as zero (Wackernagel, Onisto, and Bello et al. 1999;Lin et al. 2018). Although forest has long been identified as an important terrestrial carbon sink (CS), in the current EF calculation, only ocean is taken as a CS for assimilating carbon dioxide from fossil combustion, while no other natural CSs are considered (Wackernagel, Onisto, and Bello et al. 1999). With the development of CS study, terrestrial ecosystem such as forest and pasture is identified as an important CSs (Sarmiento and Gruber 2002;Pugh, Lindeskog, and Smith et al. 2019;Jmo and Hall 2010), which are excluded from the existing EF accounting, leading to underestimation of ecological carrying capacity and overestimation of ecological overshoot or ecological deficit. As pronounced land-use types of terrestrial ecosystem, forest and pasture make up a significant part of CS. Globally, tropical grasslands and savannas may act as essential CSs with a magnitude of 0.5 Pg C per year (Jmo and Hall 2010), and forest contributed about 2.15 Pg C per year across 2001-2010 (Sarmiento and Gruber 2002). Furthermore, land carbon sequestration may be accelerated with elevated atmospheric CO 2 concentration (Liu, Piao, and Gasser et al. 2019). In China, forest and grasslands are also significant CSs, taking up 75.2 Tg C per year and 20.94-30.94 Tg C per year, respectively (Fang, Guo, and Piao et al. 2007). Therefore, incorporating terrestrial ecosystem CSs into the EF accounting is useful to reflect real carbon adsorption in the ecosystem, which may improve the model performance and reflect the real state of sustainability. Some publications have noticed the shortcoming of overlooking the carbon absorption capacity of the terrestrial ecosystem in the EF study ), but it was only used in calculating the EF of energy consumption, and no approach for estimating its role in biocapacity has been suggested. However, it is worth noting that there are double-counting flaws if all terrestrial vegetable products are directly used to estimate CSs. Some terrestrial vegetable products such as crop, fruits, and firewood cannot act as CSs for a long time, and most carbon will return to atmosphere through rapid biological and chemical process within a year and thus not considered as important CSs (Fang, Guo, and Piao et al. 2007). Carbon sequestered in harvested wood products (HWP) usually has a much longer half-life time and acts as an important CS . It is true that carbon cannot be permanently stored in HWP and only acts as a temporary CS, but the stored carbon works as an important sink for a certain period because of lag release effect, which is a significant countermeasure for mitigating greenhouse effect ). Due to the rapid increase of HWP and its relatively slow release of carbon, the role of HWP in carbon sequestration attracts substantial attention, and it is reported that its net sequestration of carbon is increasing rapidly in many regions . With our knowledge, no attempts have been made to incorporate HWP CS in the EF accounting system to assess its impact on the biocapacity. Based on the study of HWP CS in the territorial ecosystem mentioned above, HWP is an increasingly important CS, and it should be included in the calculation of EF.
The objective of this study is to establish a new EF accounting system incorporating CS from HWP and to assess its applicability at regional scale. The established model was then applied to the Yangtze River basin to assess its impact on regional EF estimation. First, HWP sinks in the watershed were calculated based on historic statistic data. Second, the estimated CS was converted into carbon assimilation lands to estimate biological capacity in the region. Third, biological carrying capacity with and without HWP CS in the watershed was calculated and compared. Finally, the EF of the region was computed and compared to the biological carrying capacity. The revised model may renew our knowledge about the role of HWP CS in EF calculations and supply a novel insight for calculating the framework of EF.

Study case
The Yangtze River is the longest river in China with a drainage area of 1.78 million km 2 (Figure 1), ranking the third longest and most water-abundant river in the world. In 2015, the total population in the basin was 450 million, making up 32% of China. The Yangtze River basin is also an important economic region in China. Its gross demotic product was 28.8 trillion RMB, contributing 35% of China. Land use in the basin ( Figure S1) is dominated by forest (41.7%), cropland (26.9%), and pasture (23.3%). Human population, land use, and economic intensity diversify across different subbasins. Nine main subbasins were discerned in the region based on China watershed classification according to China Water Resource Second Class Zone (http:// www.resdc.cn/Default.aspx). In recent decades, increased riverine nutrient loads and significant altered land use in the basin have been found (Liu, Beusen, and Van Beek et al. 2018;Cheng, Cheng, and Wu et al. 2017), and human pressure in the basin is intensifying.

Ecological footprint analysis
EF model is consisted of EF and biocapacity. The consumption demand of the population of the region can be compared to the biological regeneration of the region by comparing the relative magnitude of the two components.
where RE is the ratio of EF to accessible BC in national hectares which is total biologically productive land minus the land reserved for biodiversity protection. When the value of RE is greater than 1, it indicates ecological deficit, while the value less than 1 indicates ecological surplus, respectively; BC is the total existing bioproductive areas in national hectares (Lenzen and Murray 2001), indicating regional biological carrying capacity. The constant 0.12 is the percent of land suggested for protecting biodiversity on the Earth from the total existing biocapacity (Wackernagel, Onisto, and Bello et al. 1999).

Ecological footprint
In the EF accounting framework, human's requirement of biologically productive space derives from six human activities (Wackernagel, Schulz, and Deumling et al. 2002): growing crops on cropland, grazing animals on grassland, harvesting timber on forest, fishing in freshwater and marine water, accommodating infrastructure on built-up land, and burning fossil fuel that requires forest to absorb greenhouse gases, which can then be categorized as three components: biological resources, built-up land, and fossil fuel emissions.
where EF is the total ecological footprint from all anthropogenic activities; EF b , EF u , and EF f refer to ecological footprint from consumption of biological resources, built-up land, and fossil fuel emissions, respectively.
The agricultural product is harvested from cropland, forest, pasture, and fishing ground. Their EF can be estimated by the following formula: where i is the kind of agricultural products humans consumed; BR is the quantity of biological products consumed (t/year or m 3 round wood/year for forest); YF is the yield factor, which can be estimated from local yield and national average yield (nha/ha); EQF is the equivalence factor; LY is the regional yield of the ECOSYSTEM HEALTH AND SUSTAINABILITY product (t/ha/year or m 3 round wood/ha/year for forest).
EF from occupied built-up land can be derived from actual land-use image and then multiplies the YF and equivalence factor of cropland (Wackernagel, Onisto, and Bello et al. 1999). Since land uses other than cropland may also contribute to the converted built-up area, the YF and equivalence factor of built-up land in this study were estimated from actual land-use conversion. Thus, the EF of built-up land can be presented as follows (Li, Zhang, and Chang et al. 2018): where UA is the built-up land area (ha); PA i is the percentage of the converted area from land use i to the total land area converted to built-up land; and YF and EQF are the yield factor (nha/ha) and equivalence factor of land use.
Carbon dioxide emitted from anthropogenic activity poses an important threat to the natural ecosystem. Its EF is quantified as the additional land area that can absorb CO 2 . Since forest is usually the largest CO 2 consumer in terrestrial ecosystems, fossil energy lands are treated as the forest area needed for taking up CO 2 . The EF f can be described as follows: where FE is the total CO 2 emitted from fossil combustion (t C/year), which is usually estimated from burning coal, oil, and natural gas; S o is the proportion of CO 2 absorbed by oceans, which is about 35.2% (Sarmiento and Gruber 2002); FS is the global average forest CO 2 sequestration (t C/ha/year), which is 0.45 t C/ha/year estimating from the reference Pugh, Lindeskog, and Smith et al. (2019); and EQF f is the equivalence factor of forest.

Biocapacity and its revised form from terrestrial carbon sinks
Natural resources are the basis for calculating carrying capacity in a specific region. Various land use can supply a lot of ecosystem services such as food and absorb carbon dioxide for human society. Since the productivity of land use is different across the world, land-use area will firstly be transformed into national or global land-use area by multiplying its YF. Different land use will then be aggregated to national or global hectares by multiplying its equivalence factor.
where i is the set of five land-use types; LA is the area of land use (ha); EQF is the equivalence factor; and YF is the yield factor (nha/ha). Different from traditional form (6), CS in the region is expressed as an important item of carrying capacity (7).
where i is the set of five land-use types; LA is the area of land use (ha); EQF is the equivalence factor; YF is the yield factor (nha/ha); CS is carbon sink of HWP in the region (t/year); FS is the CO 2 sequestration of forest (t C/ha/year); and EQF f is the equivalence factor of forest.
The ocean and the terrestrial biosphere act as the two significant sinks of carbon in nature (Sarmiento and Gruber 2002). The magnitude of CS is usually calculated as the difference in carbon stocks at different periods, which can be estimated by atmospheric inverse model (Cheng, Cheng, and Wu et al. 2017) or NPP (Fang, Guo, and Piao et al. 2007). Due to data scarcity and model uncertainty, estimating CSs at local scales is rare in the current scientific research. Most studies of CSs have been carried out at global or national scales (Pugh, Lindeskog, and Smith et al. 2019;Fang, Guo, and Piao et al. 2007). In China, CS of forest and pasture is found to be highly related to NPP (Fang, Guo, and Piao et al. 2007), based on which regional CS could be calculated (Fang, Guo, and Piao et al. 2007;Yang, Ma, and Zhu et al. 2016). It is worth noting that CS can only be estimated when NPP value is within the extent between 114 and 536 g C·m −2 ·year −1 ; otherwise, negative value will be expected, which is a great limitation for this model. What is more, the estimated CS is the total amount of products from forest and pasture ecosystem, which may cause double counting for forestland. To avoid the unrealistic estimation of CSs in the region, CS of HWP is estimated from the following model revised from references Yang, Ji, and Yang et al., 2013).
where CS is carbon sequestration rate of HWP in the region (t/year); CSS indicates the carbon stock of HWP in a specific year, CSS (1990) = 0 (t/year); k is the decay coefficient representing the annual rate at which the product goes out of use, citing from reference ; TP is the timber production (m 3 ); D refers to the density of timber (t/m 3 ); CBC is the carbon content coefficient for timber; and i refers to the year studied. Different from the current HWP carbon storage assessment, the CS estimated from formulas 1 and 2 is CS contributed from local forest production; thus, no export or import is considered here.

Calculation of YF and equivalence factor based on NPP
YF and EQF are two key factors affecting EF calculation, and they can be estimated from detailed agricultural production of local and national or global agriculture. However, it is usually difficult to get a full inventory of agricultural goods at different spatial scales, and the data vary a lot across years. In this study, NPP is used to estimate the YF and EQF for the basin following the methods suggested by references Venetoulis and Talberth, 2008;Liu, Fu, and Jiang et al., 2019).
where LNPP i is the local NPP of type i biological production area (g C·m −2 ·year −1 ); NNPP i is the national NPP of type i biological production area (g C·m −2 ·year −1 ); and NA i is the national area of type i biological production area (ha).

Data sources
Land use at 30 m resolution for the basin in 2015 was obtained from Tsinghua University Land Use Internet Database (http://data.ess.tsinghua.edu.cn/). The landuse data were then aggregated on the ArcGIS (10.2) Platform to six land-use types (Table S1). The annual NPP data at 1 km resolution from 2000 to 2010 of China are calculated by the GLM_PEM model and were downloaded from the Resource and Environmental Data Cloud Platform (http://www.resdc.cn/Default. aspx). Human population, crop productions, and other socioeconomic data at city unit are derived from associated provincial and agricultural statistical yearly book (2016), and the data were then aggregated to the basin based on area-weighted methods. The human population at 250 m resolution accessed through Resource Watch (http://www.resourcewatch.org). Forest productions and wood-based panels from 1990 to 2017 were obtained from China Forestry Statistical Yearbook and China statistical Yearbook. Land-use and socioeconomic data are two key input data for estimating CSs, EF, and biocapacity in this study. Main socioeconomic data are obtained in the year 2016 in this study. However, since land-use data cannot be obtained annually, we selected the land-use data in 2015 which is the nearest year available to calculate related index in this study. Although land use may change over time, especially at small spatial scales, the change for 1 year at about 1.8 million square kilometers (the area of the Yangtze River basin) is supposed to be negligible.

Estimating HWP carbon sinks in the Yangtze River basin
Based on timber production from the 19 regions in the Yangtze River basin, carbon sequestration of HWP was estimated at subbasin scales. In 2016, HWP in the Yangtze River basin is 2.99 Tg C/year, and there is great variety across different subbasins (Figure 2). HWP in the Lake Dongting watershed contributed 19% of the total CS in the basin, ranking the largest contributed region. However, in terms of CS ability per area, Jialingjiang ranks the lowest, while Lake Tai ranks the highest. The two regions are located at upstream and downstream the basin. It is interesting to find that  ECOSYSTEM HEALTH AND SUSTAINABILITY from upstream to downstream, CS per area is decreasing, indicating that forest protection measures are more strict in upstream regions.

Calculating YF and EQF of the Yangtze River basin
Based on multiple-year average of NPP (Table S2), YFs and EQFs for the Yangtze River basin were estimated by comparing regional NPP to national NPP (Table S3). On average, the YFs of cropland, pasture, and fishing ground in the basin are higher than 1.0, indicating that productivities of these land uses in the basin are higher than averages of China (Table 1). However, the YFs for different productive land used vary from each other. Pasture has the highest YF value, while forest has the lowest. Besides, the YFs vary greatly across different subbasins, and the YFs in the upstream subbasins generally are smaller than the downstream subbasins because a large area of snow and ice land use locating in the Tibet and Qinghai plateau ( Figure  S1). Different from the traditional method of estimating the YF for built-up land which assumed that built-up land is entirely derived from cropland, the YFs of built-up lands are estimated from actual land-use conversion. Therefore, the YFs of built-up lands in the subbasins are not the same as cropland.
The EQFs are estimated from actual NPP of China that averages over the 2000s and the land use in 2015. The EQFs of cropland, forest, pasture, fishing ground, and built-up land in the basin are 1.15, 1.83, 0.31, 0.77, and 1.20, respectively. The results of EQFs estimated in this study are different from some studies (Liu and Li 2009) in which the EQF of cropland is 1.74, and the difference is from the estimated NPP for cropland. In this study, NPP for the cropland is calculated as 499 g C m −2 year −1 , while it was calculated as 741 g C m −2 year −1 in the reference. Since only 1 year (2001) of NPP was used in the reference, the interannual variety is overlooked which may bring uncertainty in the calculation. Our estimated NPP for the cropland is similar to the result based on averages of 2001-2010 which is 412 g C m-−2 year −1 (Wang, Shi, and Wang et al. 2017), indicating that multiple-year average is important for estimating EQF by the NPP method.

Biological carrying capacity of the Yangtze River basin
Based on the actual land use in 2015, biocapacity for the Yangtze River basin is estimated as 0.5759 nha per capita (Figure 3). The BCs across subbasins diversify substantially due to the difference in land use and human density. Jinshajiang enjoys the largest biocapacity (1.225), while Lake Tai ranks the least (0.091). Apart from differences in land use across the subbasins, population density can well explain the distribution of biocapacity among subbasins ( Figure S3). In the Lake Tai watershed, human population density is 1706 capita per km 2 , ranking the highest across the subbasins. By contrast, in Jinshajiang, there is only 50 capita per km 2 .  Forest contributed the most biocapacity in the basin, making up 65.1% of the total for the whole Yangtze River basin. However, the percent of forest varies a lot across the watershed, ranging from 24.3% to 74.3%, and only in Lake Tai, the cropland outnumbered the forest. Since a certain proportion of biological carrying capacity has to be reserved for protecting biodiversity, the accessible biocapacity is only 88% of the total biocapacity ( Figure 4). Incorporating CS in the biocapacity, the revised biocapacity for the whole basin increased from 0.5068 nha per capita to 0.5338 nha per capita, increasing by 5.3%. For different subbasins, the proportion of CSs to accessible biocapacity varies from 3.3% to 16.7%. In Lake Tai and Wujiang, incorporating HWP in the CSs can increase their biocapacity by more than 10%, indicating the importance of CSs in carrying capacity assessment.

Comparison of EF between conventional and revised models
Based on food and energy consumption in the basin, EF was estimated for the whole basin and its subbasins ( Figure 5). The EF of the Yangtze River basin is 5.598 nha per capita, and it ranges from 3.904 to 9.337 nha capita −1 . Compared to the biocapacity, EF in the basin is much larger than the BC. If CSs are not considered, the RE is 11.0, indicating that EF is 10 times overshoot of the BC. The largest overshoot is found in Lake Tai basin, the  ECOSYSTEM HEALTH AND SUSTAINABILITY population of which is mostly dense in the basin. Although the RE of Lake Poyang basin is the smallest, it is still in heavily overshoot status. When HWP CSs are considered, the revised RE is smaller than the original values due to the significant increase in BC ( Figure 6). The ratio of EF to BC decreases from 11.0 to 10.5 for the whole basin, and the ratios range from 3.2% to 14.3% across different watersheds.

Estimated HWP and total carbon sinks in the Yangtze River basin
To invest the role of HWP in the whole terrestrial ecosystem CS, the carbon sequestration in the Yangtze River basin is calculated by the regression model from 1 -km resolution NPP (Fang, Guo, and Piao et al. 2007) ( Figure S2). The results showed that CS in the whole basin is 0.125 Mg C/ha/year ( Figure 5). Considering the area in the basin that can take up carbon, the total amount of CS in the basin is 29.30 Tg C/year. The terrestrial CS in China is about 96.1-106.1 Tg/year (Fang, Guo, and Piao et al. 2007). CS in the Yangtze River basin accounts for about 27.6-30.5% of China. According to land use in 2015, forest and pasture in the Yangtze River basin make up 30.7% and 13.0% of China, respectively. The CS in this study is well consistent with the spatial distribution of land uses. However, there is great uncertainty in CS calculation in China, and some land use may shift between sinks and sources of carbon (Piao, Fang, and Huang 2010). Carbon sequestration varies a lot across different subbasins. Based on CS ability per area, Jialingjiang ranks the lowest, while Wujiang ranks the highest (Figure 7). Although CS is well correlated to NPP, their relationship is nonlinear. Therefore, it is difficult to explain the distribution of CS merely from NPP. The lowest NPP in the subbasins is located in Jinshajiang watershed (296 g C m −2 year −1 ) and the highest is in Lake Poyang (958 g C m −2 yr −1 ), but they both are not the extremes of CSs. It is interesting to find that no significant linear correlation is observed between carbon sequestration from that of HWP and that of forest and pasture (p > 0.1). The CS estimated from HWP makes up 18% of forest and pasture and varies from 5% to 45% across subbasins, leaving more than 50% uncalculated CSs. Although great uncertainty exists in the CS calculation of forest and pasture (Fang, Guo, and Piao et al. 2007), other potential CSs need to be investigated and assessed in the future.

Limitation of this study
In this study, only HWP are considered as new CSs, but other CSs such as pasture and soil sinks may also act as significant CSs (Sarmiento and Gruber 2002;Fang, Guo, and Piao et al. 2007;Mccauley, Robles, and Woolley et al. 2019), and it is important to incorporate other CSs in the EF calculation when data are available. Traditional EF calculation assumes that biologically productive areas are mutually exclusive (Wackernagel, Onisto, and Bello et al. 1999), based on which carbon dioxide absorption is excluded from the functions of forest and pasture. Since it is difficult to extract firework from the total areas due to lack of data, the BC results in this study may be overestimated. The YFs and EQFs are key factors influencing EF calculation, and they are calculated from NPP. The NPP may vary significantly across different years (Zhu, Zhao, and Zhu et al. 2017), but only a short period (11 years) in this study is used, which may bring uncertainty to the estimation. The YFs and EQFs derived from NPP are different from other methods (Venetoulis and Talberth 2008;Barrahmoune, Lahboub, and El Ghmari 2019) and are not even consistent with the results from NPP due to diverse calculation methods of NPP and land-use classification (Liu and Li 2009), implying that the NPP-based YF and EQF estimation is still challenging. CS estimation based on NPP can only be used within a small range of NPP values, and it is of great need to develop a reliable approach for estimating CS at various scales.

Conclusions
Biologically productive space for sequestrating CO 2 is one of the most essential components of EF for a region, usually ranking the biggest contributor. The classic approach of the EF accounting system assumes that biologically productive land uses are mutually exclusive and forest is excluded from CS calculation, leading to underestimation of carrying capacity. In this study, a new calculating framework incorporating HWP CSs to biocapacity was established and was applied to the Yangtze River basin. Results showed that incorporating carbon sinks can significantly improve the biocapacity of the basin, and in some subbasins, CSs are enhanced by 10% compared to their original calculated biocapacity. Therefore, the role of CSs from HWP on EF accounting needs to be fully considered. The model established in this study supplies a primary method to incorporate CSs from HWP in EF accounting, which proved to be very important in intense forestry production areas. To accurately estimate the impact of CSs on biocapacity, more efforts should be done on developing a reliable approach to estimate various CSs at regional scales and to tackle the problem of exclusive assumption of biologically productive areas in EF.