Carbon stock in Japanese forests has been greatly underestimated

An accurate estimate of total forest carbon (C) stock and C uptake is crucial for predicting global warming scenarios and planning CO2 emission reductions. Forest inventory, based on field measurements of individual tree sizes, is considered the most accurate estimation method for forest C stock. Japan’s national forest inventory (NFI) provides stand-scale stem volume for the entire forested area based on (1) direct field measurements (m-NFI) and (2) prediction using yield tables (p-NFI). Here, we show that Japanese national and local forestry agencies and some research studies have used p-NFI and greatly underestimated the Japanese forest C stock (58–64%) and net annual C uptake (41–48%). This was because approximately 10% of the forest area was not counted in p-NFI and because the yield tables in p-NFI, which were constructed around 1970, were outdated. For accurate estimation of the forest C stock, yield tables used in p-NFI should be reconstructed or ideally field measurement campaigns for m-NFI should be continued. In the future, appropriate forest management plans are necessary to effectively use the high CO2 absorption capacity of Japanese forests and these should be compared with other industries’ CO2 reduction plans from a cost-benefit perspective.

Forests accumulate large amounts of Carbon (C) from the atmosphere by incorporating photosynthesis into long living woody biomass 1,2 , thereby greatly affecting the rate of increase of atmospheric CO 2 concentration and global warming 3 . The Marrakesh Accords determined that net C uptake, which reflects the temporal variation in forest C stocks, can be counted in national CO 2 reduction plans if the forest is properly managed 4 ; this ruling remains active 5 . Therefore, an accurate estimate of the total forest C stock and net C uptake at national to global scales is crucial for predicting future global warming scenarios and planning future CO 2 emission reductions.
The national forest inventory (NFI) is a national forest survey that was historically aimed to estimate timber productivity. It generally includes information on tree species composition and stand-scale or individual tree-scale stem/timber volume in targeted forests. The Nordic countries began using the NFI from the 1910s 6 . Although each country has its own survey standards 7,8 , the NFI generally provides an accurate and valid estimation of national forest C stocks and net C uptake because it is based on multipoint observations and the uncertainties of the estimated values can be calculated [9][10][11] . Using NFI data, many studies have clarified the forest C stocks and net C uptake from the past to the present 6,9,[11][12][13][14][15][16] and predicted their future changes with forest growth models [17][18][19] .
In Japan, two types of NFI have existed (Fig. 1). One (hereinafter referred to as p-NFI) is based on predictions from yield tables, which are empirical relationships between stand age and the stand stem volume. These yield tables were constructed using actual observation data obtained around 1970. The other is based on direct field measurements (hereinafter referred to as m-NFI). Using 10,000 survey points m-NFI surveys were carried out in 1961, 1966 and, following a long gap, 1999, based on revised sampling methods 20 (Fig. 1). Despite its broader use, p-NFI tends to underestimate the actual stem volume 21 . The m-NFI conducted in 1999-2003 was also observed to underestimate the stem volume 22 . However, for the m-NFI in 2009-2013, experts checked the accuracy and surveyors underwent technical training. Thus, the m-NFI in 2009-2013 was the most reliable NFI in Japan.
In fact, few studies with scientific accuracy and of academic value have evaluated Japan's nationwide C stock. In this study, we assumed that even the few studies recognized as academically valid have committed a fallacy: some studies estimated the C stock using p-NFI, which resulted in an apparent underestimation 23,24 , and some examined the climate and environmental change effects on C uptake by mistaking p-NFI for m-NFI, which caused them to misunderstand the effect [25][26][27][28] . These studies might have misled policymakers in their Japanese forest management policy and legislation. Hence, our main objective in this study was to offer an accurate estimation

Results and Discussion
C stocks estimated from m-NFI in 1961 and 1966 were 851.7 ± 23.0 (the plus-minus sign represents the 95% confidence interval) and 834.0 ± 22.5 TgC, respectively, which lie between values from p-NFI in 1956 (774.2 TgC) and 1975 (905.8 TgC) (Fig. 2). Until the 1970s, C stocks estimated from p-NFI were expected to be equivalent to those estimated from m-NFI, because the yield table for each p-NFI was constructed using the field data at that time 23 . Thus, we can confirm that the total forest C stocks from the 1940s to the 1970s were not more than 1000 TgC and were somewhat constant (that is, the annual net C uptakes were small).
The m-NFI-derived C stock in 2011 was estimated using two types of BEF (BEF 1 and BEF 2 ; see the Methods) to give 3016.2 ± 26.9 and 2696.4 ± 25.4 TgC, respectively (Fig. 2). The p-NFI-derived C stock in 2012 was 1750.0 TgC, which was considerably smaller than the m-NFI values. Additionally, considering the validation surveys since 2010, the m-NFI probably underestimates the stem volume by 4-8% 29,30 . This suggested that the difference between the p-NFI-and m-NFI-derived C stocks would be even greater than the present study showed. Therefore, we can conclude that the p-NFI-derived C stock greatly underestimated the actual forest C stock in Japan. Year Figure 1. Survey years of National Forest Inventory (NFI) in Japan. p-NFI and m-NFI denote NFI based on the predictions from yield tables and the direct field measurements, respectively. Open circles mean that the surveys were done with different methods to the current ones. Each m-NFI survey has a 5-year field measurement period; therefore, m-NFIs that began in 1999 are plotted at their central years. Temporal changes in total forest carbon stock in Japan. p-NFI and m-NFI denote National Forest Inventory (NFI) based on the predictions from yield tables and the direct field measurements, respectively. The carbon stocks until 1995 were estimated by Fang et al. (2005) 25 . BEF 1 and BEF 2 represent two different biomass expansion factor (BEF) types. BEF 1 and BEF 2 were used for all calculations and a calculation from m-NFI in 2009-2013, respectively (see Methods section). We did not show the 95% confidence interval because it was invisibly small.
Such a discrepancy between p-NFI-and m-NFI-derived C stocks resulted in a large difference in the average net C uptakes between p-NFI and m-NFI: 19.9 TgC year −1 for 1966-2012 using p-NFI and 48.5 ± 0.8 and 41.4 ± 0.8 TgC year −1 for 1966-2011 using m-NFI with BEF 1 and BEF 2 , respectively. The m-NFI-derived net C uptake was more than twice that of the p-NFI value. The area-based m-NFI C uptakes with BEF 1 and BEF 2 were 1.83 ± 0.03 and 1.56 ± 0.03 MgC ha −1 year −1 , respectively, and that of p-NFI was 0.84 MgC ha −1 year −1 (Here, we assumed the forest area in 1966 was equal to that in 2012 p-NFI, 23.7 million ha). The net ecosystem productivity (NEP) observed using a micrometeorological method (the eddy-flux measurement) at seven forest sites across Japan ranged from 1.23 to 3.88 MgC ha −1 year −1 31 . The NEP estimated using a carbon cycle model and eddy-flux data was more than 1 MgC ha −1 year −1 in most areas of Japan 32 . Additionally, the NEP of Japanese deciduous broadleaf forest estimated from MODIS sensor and eddy-flux data was 3.47 ± 2.88 MgC ha −1 year −1 33 . These NEP values were comparable with the net C uptakes estimated from m-NFI, suggesting the validity of m-NFI-derived net C uptake and C stock.
We also estimated C stock and the area-based C stock (C density) for major forest types in Japan using m-NFI with BEF 1 and BEF 2 (Fig. 3). Cryptomeria japonica, Chamaecyparis obtusa, Fagus and Quercus have almost the same level of C density (ca. 120 MgC ha −1 ), which implies their relatively higher productivity (Fig. 3a). Cr. japonica and Quercus (famil) in Japan had the dominant C stocks in those forest types (Fig. 3b). In terms of an estimation error caused by the adoption of BEF 1 or BEF 2 , we found that the BEF 1 estimations were greatly larger than the BEF 2 ones for Cr. japonica (5.6%), Lalix (5.8%), Fagus (5.5%), Quercus (33.9%), other conifers (15.6%) and deciduous broadleaf trees (8.3%), and greatly smaller for Pinus (16.1%; Fig. 3). Due to the difference of the C content in the calculations with BEF 1 and BEF 2 (Eqs. 2 and 3), the C stocks with BEF 2 were 2% higher and 4% lower than those with BEF 1 for conifers and broadleaf forests, respectively. Even considering the effects of the C contents, there were still large differences between the C stocks estimated with BEF 1 and BEF 2 in some forest types. In particular, the C stocks of Quercus using BEF 1 and BEF 2 were estimated as 676.1 ± 14.1 and 447.1 ± 9.7 TgC, respectively (Fig. 3b). The difference between the two accounted for most of the difference between the total forest C stocks estimated with BEF 1 and BEF 2 (ca. 300 TgC; Fig. 2). Because the C stock of Quercus seems comparable to that of Cr. japonica, successfully adopting the BEF method for Quercus might be critical for the estimating Japan's total forest C stock.
Only two tree species, Cr. japonica and Ch. obtusa, the major plantation tree species for timber production in Japan, comprise 34.5% of Japan's total forest C stock (Fig. 3). We estimated the Cr. japonica and Ch. obtusa C stocks in 1990, 1995 and 2005 using m-NFI with BEF 1 (see Methods section) and compared to those in two previous studies 23,28 . One previous study used a simple growth function with parameters derived from raw data for national stand density diagrams in the 1970s 23 and the other simply used p-NFI 28 . The estimated C stocks were 452.1, 507.0 and 616.7 TgC for Cr. japonica and 236.8, 265.5 and 323.0 TgC for Ch. obtusa, in 1990, 1995 and 2005, respectively, which were much higher than those in the previous studies (Fig. 4). Again, these results demonstrated that the past C stock calculations with p-NFI have greatly underestimated the actual forest C stock. www.nature.com/scientificreports www.nature.com/scientificreports/ There are several possible factors for the large discrepancy between p-NFI and m-NFI. One potential factor is the difference in forested areas measured in p-NFI and m-NFI. Although the area categorized as forest in m-NFI in 2009-2013 was 25.9 million ha, the total area of planted and natural forests in p-NFI was 23.7 million ha in 2012 34 . This difference results in a C stock underestimation of approximately 10%. In p-NFI, 1.6 million ha of forested area was reported as the sum of treeless land and bamboo forest, which would be a major cause of the difference (ca. 2.2 million ha) between p-NFI and m-NFI. The majority of this 1.6 million ha was treeless land, because bamboo forest is estimated to occupy only 1% of the total forest area (~0.3 million ha; See Methods section). Fundamentally, when harvesting or planting forests in Japan, submitting a report to the local government is required. If no notification of planting is made after harvesting, the area will be treated as treeless land. It is quite possible there exist forests in some treeless land because of the natural growth of trees. Additionally, the forested area may have increased slightly because of conversion from agricultural land to forest. Only if the owner reported the change, the area could be reflected in p-NFI.
Another potential factor is that the yield tables of p-NFI, which had been constructed around 1970, were outdated with the following reasons.
(1) The underestimation of older forests' growth. There were few older forests in Japan around 1970 because many forests were clear-cut after World War II to meet increased timber demand 35 . Furthermore, the remaining older forests at that time might have been so poor growth that they were not worth harvesting. Thus, older forests' growth characteristics were not reflected in the yield tables, which caused C stock underestimation for the older forests and the total forests. Although the latest studies found a higher growth rate than expected in older forests 36,37 , few studies on older forests' growth have been conducted in Japan 38,39 . Therefore, intensive observations of older forests should be performed in the future and reflected in the yield tables. (2) The increased forest growth due to changes in environmental factors. Many experimental studies have shown that rising air temperatures can enhance vegetation growth 40,41 . Besides, the net primary production of Japanese forests was strongly constrained by air temperature and solar radiation 31,42,43 . Therefore, it is highly probable that the rate of Japanese forest growth increased with the recent rise in air temperature. Furthermore, Nitrogen (N) deposition has increased in Japan since the 1980s because of emissions from China 44 . Of the non-N limitation conditions, increasing atmospheric CO 2 concentration possibly enhanced forest productivity [45][46][47][48] . (3) The lack of thinning. A main purpose of the yield table is to predict stem/timber production volume for a target commercial tree species on condition that common forestry practices are conducted in a given environment. However, a decline in timber prices in Japan over many years has reduced the motivation of forestry practices, and in particular, the lack of thinning has become a serious problem 35,49 . Some previous studies reported that the stand-level forest stock of non-thinned plots was comparable with 50 or significantly higher than that of thinned plots [50][51][52] , even decades after thinning. Therefore, the lack of thinning in Japanese forests might have led to an increase in forest C stock. Year C stock (TgC) www.nature.com/scientificreports www.nature.com/scientificreports/ concluding Remarks Our study shows that p-NFI, which Japanese public forestry agencies have believed to be authentic, has underestimated the m-NFI-derived actual total C stock in Japanese forests. The main causes are the underestimation of the forest area of 2.2 million ha and the underestimation of the forest biomass due to the outdated yield tables. Our results have implied that reconstructing the yield tables in p-NFI and ideally changing the timber production estimation system completely from p-NFI to m-NFI is necessary for accurate estimation of forest C stock. In principle, m-NFI is more suitable for accurate estimation of forest C stock and decision-making than p-NFI because m-NFI estimates forest area more accurately than p-NFI and calculates uncertainties of the results. We should encourage Japanese forest management policymakers and legislators to adopt the revised timber production/ forest C stock estimation system.
The present m-NFI system began in 1999; the fourth survey of 2014-2018 is now complete and is currently being analyzed. The accuracy check in m-NFI is ongoing, but the latest result shows an annual mean difference between the actual and estimated stem volume for 2014-2016 of 5-7% 30 , which suggests a successful effort to improve the accuracy of m-NFI estimation. At present, only stand-area-based stem volume data are available in the m-NFI data supply system. Therefore, we applied the stand-area-based BEF for a single dominant tree species in a given stand to the estimation of forest stand C stock, which is a potential source of error. The wider availability of individual tree-scale configuration data such as stem volume, tree height and diameter at breast height (DBH) will allow us to more accurately estimate forest stand C stock. Furthermore, these data can describe the status of Japanese forests in detail, revealing the factors that yield tables underestimating current forest stocks. Elucidation of these factors is one possible method for accurately correcting past Japanese forest C stocks and predicting future changes. It could help policy makers in their decision-making, and therefore, is an important future task.
As previous studies noted 53 , BEF use can cause errors when estimating the forest C stock (Figs. 2 and 3). In particular, we found a significant difference in C stock estimation between BEF 1 and BEF 2 for the Quercus family, which is a major forest type and has significant biomass in Japan (Fig. 3). Again, the non-availability of raw data is a major obstacle. Integrating the data used for BEF 1 and BEF 2 would enable us to obtain BEF more accurately.
Biomass and C dynamics in forest ecosystems in warm-cool temperate zones, which have the largest population and the most active human life in the world, have received increasing attention. Therefore, various studies relating to the C uptake process, such as eddy-covariance measurements, remote sensing and terrestrial ecosystem process modelling, have been conducted. We emphasize that forest inventory based on direct field measurements conducted at numerous plots is the strongest tool to validate results from such studies. The present study can help the researchers of these related fields to correctly understand and properly use the current NFI in Japan and will contribute to the progress of C dynamics studies in warm-cool temperate zones.

Methods
National forest Inventory based on direct field measurements: m-NFI. In one campaign of the current m-NFI survey, direct field measurements were conducted for 5 years at more than 10,000 plots throughout Japan 20 . These plots formed a 4-km grid interval across Japan. Measurements were conducted at given plots including forests. The third m-NFI survey was conducted in 2009-2013 (Fig. 1); the data can be downloaded from the Japan Forestry Agency website (http://www.rinya.maff.go.jp/j/keikaku/tayouseichousa/). The number of planned plots was initially 14,838; actual forest plots were then selected, which resulted in 14,719 plots (14,522, if eliminating bamboo plots). Practically, however, the number of the target plots was 13,357 because of their accessibility.
Each plot comprises triple concentric circles with areas of 0.01, 0.04 and 0.1 ha. In the largest (0.1 ha), medium (0.04 ha) and smallest (0.01 ha) circular areas, the DBHs of trees ≥18, 5 and 1 cm DBH, respectively, are measured. Tree height is measured for each of the selected 20 trees in the plot and tree heights of the remaining trees are estimated using the DBH-tree height relationships derived from the 20 tree measurements. Stem volume of a single tree is calculated from the empirical relationship between DBH and tree height for each tree species and its growing region 54,55 . Note that the publicly available data do not include individual tree scale information as mentioned above, but rather the plot scale totalized information.
This study used only the published data, that is, stem volume per hectare (Vs; m 3 ha −1 ), name of dominant tree species and forest age for a given plot. Here, the dominant tree species denotes a tree species with the largest total cross-section area at breast height in the plot. The total number of plots used for the survey was 12,457 because we excluded 703 plots in which not all of the three types of the data were identified and 197 bamboo plots for which the C stocks were difficult to estimate. Measurement validation commenced in 2010. The annual mean differences in Vs between the practical measurements and the validation surveys were 4%-8% in 2010-2013 29 .

National forest inventory based yield table predictions: p-NFI. The yield tables used in p-NFI
show empirical relationships among stand-scale stem volume, age class (ranked from 1 at 5-year intervals) and land-productivity index for each major tree species. Although the land-productivity index is critical for estimating the stem volume in a given stand, no clear standard existed for many target forest stands on how the land-productivity index was determined. As far as we investigated, no one in national and local forestry agencies knows how to revise the land-productivity index, which can be affected by various environmental changes, which are a large source of error for the stem volume estimation. Note that the currently-used yield tables were constructed using the actual measurements from around 1970.
Yield tables for Japan differ depending on the type of forest-ownership. For national forests (ca. 30% of the total forested area), yield tables created by the national government have been used. However, for local-government-owned and private forests (ca. 70%), each of 47 prefectural governments to which the given forests belong created its own yield tables. Surprisingly, most yield tables have not been published and even the (2020) 10:7895 | https://doi.org/10.1038/s41598-020-64851-2 www.nature.com/scientificreports www.nature.com/scientificreports/ national forestry agency, that is, the Japan Forestry Agency, cannot know all the information on their yield tables. The Japan Forestry Agency has aggregated the data submitted by administrators for national and prefectural forests and published estimations of the stem volume with some relevant information (http://www.rinya.maff. go.jp/j/keikaku/genkyou/index1.html).
The forests were categorized into two types in p-NFI: plantations and natural forests. Data for plantations comprises stem volume and forest area, which are further categorized into 10 forest types, 47 prefectures and 19 age classes. Data for natural forests include only stem volume categorized into 11 forest types and 47 prefectures. Notably, no information on forest area for each forest type is included for natural forests. carbon density estimation using biomass expansion factor (Bef). BEF is the coefficient derived from allometry measurements, which enables us to convert Vs to the stand scale C stock. BEF has several definitions, and one type of BEF calculates total biomass (sum of stems, branches, leaves and roots) and another type calculates aboveground biomass (sum of stems, branches and leaves). Previous studies in Japan have used both types of BEF; here, we refer to these as BEF 1 and BEF 2 . BEF 1 was suggested in a study examining long-term changes in Japanese forest C stocks 25 We used the past forest C stock values that were calculated with BEF 1 in a previous study 25 ; therefore, BEF 1 is indispensable for the present study. BEF 2 was used in the official reports from the Japanese government 56 and is probably a more commonly used BEF in Japan than BEF 1 . Therefore, we estimated all forest C stocks with BEF 1 to clarify the temporal changes in forest C stock and used both BEFs for m-NFI in 2009-2013 to quantify the uncertainty derived from the difference in BEF. The conversion methods using these BEFs are shown below.
BEF 1 (Mg m −3 ) is for the conversion from Vs to total biomass per hectare and was presented for 10 plantation-forest types in p-NFI. According to the evidence that BEF 1 tends to increase with decreasing Vs, the following equation is reached: 1 where a (Mg m −3 ) and b (Mg ha −1 ) are empirical parameters obtained for a given forest type. In the case that Vs is extremely small in Eq. 1, BEF 1 becomes unrealistically large. Therefore, we set an upper limit of 1.5 for BEF 1 according to data obtained in a previous study 25 . In addition, in previous studies using BEF 1 , C content was defined as 0.5 25 . As a result, forest stand scale C density (Cd; MgC ha −1 ) is given as follows: 1 BEF 2 (unitless) is the conversion factor to expand stem biomass to the aboveground biomass (that is, stems, branches and leaves) and is presented for two age classes (below and above 20 years old) of 32 tree species, some alien species, other conifers and other broadleaf trees (see Table S1). Cd using BEF 2 can be written as follows: 2 where WD is woody density (Mg m −3 ) and R (unitless) is the ratio of root biomass to aboveground biomass for a given tree species. Stem biomass ( × V WD s ) is converted into aboveground biomass by BEF 2 , then multiplied by ( + R 1 ) to calculate total biomass. α is C content: 0.51 and 0.48 for conifers and broadleaf trees, respectively 56 .
carbon stock estimation using m-nfi. The m-NFI-derived C stocks in 1961 and 1966 have been already reported in a previous study 25 ; therefore, we estimated C stock for m-NFI in 2009-2013. Again, the numbers of actual forest and bamboo plots were 14,522 and 197, respectively. We calculated Cd at 12,457 of the 14,522 forest plots using both BEF 1 and BEF 2 , which were determined for the dominant tree species at the plot. We then obtained plot-averaged Cd for each of the 11 forest types categorized in p-NFI. The total forest C stock was estimated as a total of each forest type C stock, which was calculated by multiplying the forest type averaged Cd and forest type area for the 14,522 forest plots. Surveyors could not reach or tree species could not be identified in 2,065 of the 14,522 forest plots. We assumed the 2,065 plots to be natural forests but could not apply a forest type (using the p-NFI categories). This means that while we could estimate forest type areas for 12,457 observation plots, those for 2,065 plots are unknown, which resulted an incorrect total forest type area in the 14,522 total forest plots (see Table S2).
Therefore, it is necessary to estimate the forest type areas for these 2,065 plots and the total for the 14,522 plots (Table S2). Assuming that the areas of each forest type in the p-NFI natural forest category are proportional to each forest type's estimated stem volume in the natural forest category in the 2012 p-NFI, we assigned 11forest types to the 2,065 plots. The sum of the assigned plots for each forest type in the 2,065 non-observation (natural forest) plots and the 12,457 observation plots can be assumed to represent the area ratio of each forest type in the 14,522 total forest plots (Table S2). The forest area of each forest type can be calculated by multiplying Japan's total forest area and this area ratio. Using aerial photographs, satellite images and field surveys, Japan's total forest area was estimated to be 25.86 million ha (http://www.rinya.maff.go.jp/j/keikaku/tayouseichousa/tikuseki.html). The ratio of the bamboo forest plots to the actual forest plots (~1%) was subtracted from the total forest area, leaving 25.6 million ha as the total forest area. carbon stock estimation using p-nfi. The p-NFI-estimated C stocks until 1995 have been already reported in previous studies 25,28 ; therefore, we estimated C stock for the subsequent three p-NFIs in 2002, 2007 and 2012 using the same method as the previous studies. We calculated Vs at plantations by dividing the p-NFI-estimated stem volume by the forest area for each factor, that is, prefecture, forest type and age class. Then, the Cd was obtained by using the BEF 1 calculation with the Vs (see Eqs. 1 and 2). The C stocks at plantations were estimated by multiplying the Cd and the forest area for each factor and summing the products. (2020) 10:7895 | https://doi.org/10.1038/s41598-020-64851-2 www.nature.com/scientificreports www.nature.com/scientificreports/ No information is given in the p-NFI natural forest category about forest area; therefore, we estimated the C stock in natural forests without using the Vs. We first obtained plantations' BEF 1 values for each prefecture and forest type and used these in place of those for natural forests. Then, the C stock in natural forests was estimated by multiplying the BEF 1 and the total stem volume for each prefecture and forest type in natural forests and summing the products. Japan's total forest C stock is determined by the sum of the C stocks in plantations and natural forests.
Calculation of 95% confidence interval (CI). We calculated the 95% CI for m-NFI-derived C stocks. No information on the C stock uncertainty was available in the previous study, which used the 1961 and 1966 m-NFIs to calculate the forest C stocks 25 . Thus, we assumed that the widths of the 95% CIs of the total forest C stocks calculated from the m-NFI of 1961 and 1966 were equal to that of total stem volume in the 1961 m-NFI, i.e., 5.4% of the mean value 57 . For m-NFI in 2009-2013, the 95% CI of the area-based C stock was calculated for each of the 11 forest types (Table S2) as follows: i i i where ε i , σ i , and n i represent the half width of 95% CI (MgC ha −1 ), standard deviation (MgC ha −1 ), and sample size of area-based C stock of forest type i, respectively. Then the half width of 95% CI of total forest C stock (ε T ; MgC) can be calculated with the forest area of forest type i (A i ; ha) as follows: Reproducing past carbon stocks of Cr. japonica and Ch. obtusa. Assuming that for a given forest type i (such as Cr. japonica) the net forest C uptake flux (F i ; MgC ha −1 year −1 ) is proportional to the area-based C stock (V i /A i , in which V i and A i are C stock (MgC) and forest area (ha), respectively), the net C uptake of the forest type i (F i × A i ; MgC year −1 ) may be approximated as follows: where F T and V T are net total forest C uptake (MgC year −1 ) and total forest C stock (MgC), respectively. Thus, the past C stock of i at the year Y 1 (Cst Y1, i ; MgC) can be calculated as: where Y 0 is the reference year and Cst Y0, i is the C stock (MgC) at Y 0 .