Estimation of CO 2 Sequestration by the Forests in Japan by Discriminating Precise Tree Age Category using Remote Sensing Techniques

This study estimates CO2 sequestration by forests in Japan using Land Remote Sensing Satellite (Landsat) Operational Land Imager (OLI) and the Advanced Land Observing Satellite (ALOS) Phased Array L-band Synthetic Aperture Radar (PALSAR) remote sensing data for the in-depth retrieval of forest growth stages (tree age). Landsat imagery was used to develop a detailed forest cover map, while the PALSAR data were used to estimate the volume information. The volume was converted to tree age information for each of the three forest types in Japan. An estimation of CO2 sequestration values for each forest type and for each tree age from the forest inventory data was made. The forest cover map results in four classes, and the overall accuracy yields approximately 74%. For the volume estimation, Root Mean Square Error (RMSE) was computed with the ground reference information resulting in 105.58 m/ha. The final result showed that total CO2 sequestration in Japan based on tree age forest subclasses yields 85.0 Mt∙CO2 (coniferous), 4.76 Mt∙CO2 (evergreen broadleaf) and 21.61 Mt∙CO2 (deciduous broadleaf), which in total is 111.27 Mt∙CO2. Using remote sensing techniques to quantitatively estimate CO2 sequestration in Japanese forests has been shown both to have advantages and to offer further possibilities. OPEN ACCESS Remote Sens. 2015, 7 15083


Introduction
Forest resources are essential to life on Earth because forests play a large role in providing organic matter through photosynthesis, which is important for countless species of plants and animals, and in processing atmospheric gasses such as CO2, which is of crucial environmental importance and among the most important issues faced in the 21st century: Global Warming/Climate Change.Research is performed and conferences are held to investigate this issue and its associated problems, and many nations are working together on policies to prevent or mitigate the associated problems and impacts [1,2].Japan has been at the forefront of nations involved in the many activities related to this issue.The Kyoto Protocol and the Reducing emissions from deforestation and forest degradation mechanisms (REDD, REDD+) provided a new approach to finding solutions, giving the role of forests more importance, especially when focusing on how forests function as a carbon sink.Forest carbon sinks are highly important in the current research context, as 67% of the land in Japan is covered by forests [3].Japan considers forests to be the top candidate for its CO2 reduction process [4].Prefectural Governments and private companies in Japan highly value forest carbon sequestration research; however, some Prefectures have not yet been involved in such processes even when they have a great amount of forest cover [5], which may be because they lack the equipment or effective methodology to quantitatively estimate carbon sequestration.
The estimation implemented in Japan mostly focuses on identifying forests as either natural or planted.Hiroshima and Nakajima [6] focused on the potential carbon sinks in Japanese plantation forests during the first commitment period of the Kyoto Protocol.They estimated that planted forests were expected to sequester 8.16-8.87Mt•C/yr, depending on various scenarios of forest management based on silvicultural practices, and employed subsidies and forest workers' wages as predictor variables.Sasaki and Kim [7] estimated potential CO2 sequestration by forests in Japan and the eligible sequestration under the condition of the Marrakesh Accords.Using the land use model and the carbon stock growth model, it has been estimated that the forests in Japan are likely to sequester 20.1 Mt•C/yr (planted forests: 15.3 Mt•C/yr, natural: 4.8 Mt•C/yr), and under the conditions of the Marrakesh Accords, it is estimated at 10.2 Mt•C/yr (planted forests: 7.3 Mt•C/yr, natural: 2.9 Mt•C/yr).Recently, the National Institute for Environmental Studies (NIES) Center for Global Environmental Research (CGER) calculated a result for the IPCC reports on the issue of forest carbon sinks, estimating the CO2 sequestration by the forests of Japan [8].Using the forest registration data from each prefectural government, and differentiating the forests by tree species, they have calculated the total CO2 amount sequestered both aboveground and in litter and soil.Their results showed an estimate of 77.67 Mt•CO2 for the year 2012.
Utilizing satellite remote sensing data for quantitative C/CO2 estimations is challenging but well-known method for improving broad scale estimations compared to conventional ground-based observations [9].The techniques can be separated into two major methodologies.The first is the indirect measure of C/CO2, which calculates several parameters related to the function of the forest system, such as Leaf Area Index (LAI), Photosynthetically Active Radiation (PAR) and others, and then models multiple parameters to estimate either Gross Primary Product (GPP) or Net Primary Product (NPP) [10,11].The second method integrates detailed land cover information with ground observations of forest inventories [5].Alternatively, both methods may be integrated together [12].All these methods have advantages and disadvantages.For example, one popular product by the former method is the Moderate Resolution Imaging Spectroradiometer (MODIS) NPP [13], which can be collected globally with a fine temporal resolution.Depending on the need, such products can be unsatisfying because the spatial resolution of 1 km by 1 km is too coarse; when we want to focus more on a regional-scale analysis, a spatial resolution of 1 km cannot delineate the actual local environment accurately.The latter method is usually implemented in a regional scale analysis using fine resolution imagery, such as Landsat, to discriminate details in land cover/land use types to estimate accurate values for C/CO2 content.This can also be conducted by indirectly estimating modeling parameters, but it is often performed with the associated forest inventory data collected from other ground observation sources and then integrated with the land cover information for the estimation.In comparison with MODIS, Landsat provides a much finer spatial resolution (30 m), making it possible to allocate the heterogeneity of the natural environment; however, there are again issues that must always be faced, such as the fundamental problem of utilizing optical imagery due to cloud cover.Fine resolution imagery has lower temporal resolution, and it is sometimes difficult to collect cloud free imagery through only a few observations.Sanga-Ngoie et al. [5] has utilized synthetic aperture radar (SAR) data to extract the forest parameters related to CO2 sequestration.This method has shown the usefulness of implementing SAR to observe forest extent, which could be broadly usable in conditions where there is frequent cloud cover over the area of interest.This method made it clear that not only land cover type discrimination, but also tree age, or the growth stages of the forests, plays an important role.Taking these variables into account provides more accurate estimations for the overall result, because each forest growth stage differs in its sequestration rate [5,14].Utilizing both optical and microwave satellite imagery makes this possible.
Most forest carbon sequestration capacity estimates within Japan are performed at a national scale, and detailed analyses at a regional scale are still rare; therefore, such estimates are either rough estimations on a national scale or detailed analysis on a local scale.Unfortunately, local-scale studies cannot be expanded to a national scale; they are only for smaller regions and do not apply to the country as a whole.The studies mentioned usually use a methodology involving a modeling base, site observation and statistics.However, detecting and locating forest resources does not take place.Locating precise forest resources of both forest types and their growth stages requires not only consideration of realistic estimation of sequestration [5,10,[14][15][16] but also an understanding of the current status of the forest resources, which may be the most important information to gather for long-term sustainable plans to manage forests in a strategic way.Consequently, the objective of this research is to reassess the estimation of CO2 sequestration (i.e., NPP) by forest cover in Japan based on a precise evaluation of forest extent by both forest type and age, together with the use of remote sensing data, to quantitatively estimate the CO2 sequestration potential.

Study Area
Our area of study focuses on mainland Japan (except Okinawa and small islands) (Figure 1).Japan is made up of over 7000 contiguous islands from the northeast to the southwest, with a total length of approximately 3000 km.Total land area is approximately 377,900 km 2 .Over 60% of the land is recognized as mountainous (over 2000 m in altitude) and 12% as hilly.The total forested area makes up approximately 67% [3], the majority of which is located around the aforementioned mountain or hill areas.Climatic conditions in Japan vary widely due to its geographical location; its lands stretch from a southern latitude of 30°N (Okinawa omitted) to a northern latitude of 45°N.This range, together with the ocean currents, gives each region a different climate.Because of the variety in topography and climatic conditions, Japan has a great diversity of vegetation both horizontally and vertically.For example, there are over 6000 species of Pteridophytes, each of which has evolved and adapted to its specific location [17].The dominant coniferous tree types in Japan are the Japanese cedar (Cryptomeria japonica) and cypress (Chamaecyparis obtusa) trees, which occur widely from the north to the south.As a whole, coniferous and broadleaf forests are distributed in approximately a 50/50 ratio, while 80% of coniferous forests tend to be planted forests that were part of a large plantation plan that lasted from the 1950s to the 1970s.On the other hand, the broadleaf forests are primarily natural (over 90%), and few deciduous or evergreen broadleaf trees have been planted.Broadleaf hardwoods are dominant in the central region, while mixed deciduous evergreen forests or evergreen broadleaf forests can be observed in the temperate south zones.Remote sensing data for Chiba Prefecture was used for a detailed analysis of statistical data, so a brief explanation is given.Chiba is located on the east coast of Japan along the Pacific Ocean, just east of the Tokyo metropolitan area, where the peninsula juts out.Climatic conditions include a warm oceanic climate, which is a condition with high humidity and high precipitation in summer, and low humidity and low precipitation in the winter.However, comparing the southern and northeastern regions of Chiba, the climate is warm throughout the year in contrast to the inland area, which shows more diversity such as larger temperature drops in the winter [18].Annual precipitation is highest in the southern area, with more than 2000 mm, and is next highest in the northern area, which has approximately 1400 mm to 1600 mm.The difference in precipitation clearly shows the distribution of the forests, which are more dense in the south and sparse in the north [18].A characteristic of Chiba is that it is one of the flattest regions in Japan; the average land altitude is 46 m above sea level, and there are no mountains higher than 500 m in altitude.

Methods and Data
The overall flowchart of the methodology is shown in Figure 2. A three-step process is used: (1).Estimating the volume of the forests in Japan using SAR and converting these data to tree age information data; (2).Developing a detailed land cover (forest cover) map to allocate the precise forest types across the land area of Japan; and (3).Utilizing the obtained data from the preceding steps together with the CO2 sequestration rate per unit area to estimate the sequestration values by each forest type and age.

L-Band PALSAR for Volume Mapping
In the first step of the study, a simple relationship analysis was performed between the biophysical parameters collected from field observations by the Chiba Prefectural Government and the L-band backscattering intensity data collected by the Advanced Land Observing Satellite (ALOS) Phased Array L-band Synthetic Aperture Radar (PALSAR) provided by the Japan Aerospace Exploration Agency (JAXA).This step involved extracting the biophysical parameters (volumes) of the forested areas and mapping them for the whole of Japan.The use of SAR for detecting forest parameters works because the radar waves are longer than the optical waves, so the radar beam penetrates through the top layer of canopy and scatters from the stems and branches of the trees when using L-band radar [19][20][21][22][23].A number of studies have been conducted to develop models related to those studies' specific regions of interest.Those studies confirm that the backscattering signature correlates with forest parameters and that SAR images can be used to extract information related to forest biomass.While most of the papers present and acknowledge the possibility of extracting forest parameters using the L-band SAR information, some studies show difficulties in data interpretation [24], mostly due to distortions caused by local topography, which is also one of the challenges for this work.
To meet this topography challenge for the analysis, Chiba Prefecture was selected due to its relatively flat terrain, which limits this critical and fundamental topography problem when utilizing SAR imagery [25].
The 25 m PALSAR Global Mosaic (PGM) image used in this study was converted from the provided Digital Number (DN) values to the backscattering intensity information, also known as backscattering coefficients, using Equation (1).The PALSAR specification is in ascending Fine Beam Dual (FB) mode characterized by a 34.3 degree of off-nadir angle, calculated as: where σ 0 is the backscattering intensity represented in decibel units (dB) and CF is the calibration factor for the data obtained, depending on the observation period and polarization [26].The image is already processed for ortho-rectification and slope correction [27], so further analysis uses these data directly without any additional correction.
A statistical analysis for formulating relationships among the field observation data of the forest stand characteristics and the remotely sensed microwave satellite data was carried out.The process analyzes the relationships between SAR backscattering intensity (σ 0 ) and the stand characteristics (volume) for the 25 m pixel spacing PALSAR image.The PALSAR image used here is the mean value of the years 2009 and 2010, to minimize local variations from such things as precipitation.Using the vegetation continuous field data [14], the sample plots were selected from an area of less than 70% forest cover to avoid uncertainties associated with scattering, as Iizuka and Tateishi [28] explain that attenuation of the radar beam is caused by high dense forested areas, resulting in lower backscattering values compared to the peak areas (70%-75% forest cover), and causing difficulties when interpreting the actual forest trend.
Using the result, the model function was applied to the 50 m PALSAR mosaic product (which can be downloaded free from the JAXA web page [29]) to obtain volume information for the forests across Japan.The 50 m PALSAR data were selected because of the spatial coverage and availability (it is freely available and covers the whole of Japan).The equation obtained from the relationship analysis was assumed to be acceptable to apply to the downscaled data, because several papers confirm an increase in accuracy and higher correlations when the pixel spacing is reduced [30,31].
DN values of the data from 2007 and 2009 were averaged to avoid small changes from backscattering caused by local environmental factors at that particular observation time (e.g., precipitation) and then converted to dB units.The data were then slope corrected, because they had been previously ortho-rectified but not slope corrected, resulting in enhancements of the backscatter from the rugged hills facing towards the sensor direction [28].
A correction method proposed by Zhou et al. [32] was modified and applied to the 50 m PALSAR mosaic product.The correction of the SAR imagery is performed by applying a slope correction factor to the original SAR data.The equation for the correction is as follows: where σ 0 c is the corrected sigma naught backscattering value, σ 0 is the original sigma naught backscattering data, q is the diffuse scattering, which is polarization dependent, and SCF is the slope correction factor for the correction of the topography, expressed as: where cos i is the cosine of the incidence angle, α is the slope and β is the aspect.The HV polarization image of PALSAR was corrected using this method.The value for the diffuse scattering depends on the environment and was unknown in this case due to the large study area with differing forest structures, so here it has been assumed that q = 1.The volume map for the whole of Japan was obtained by applying the volume model equation to the corrected 50 m PALSAR mosaic product.

Land Cover Mapping of Japan
We selected and analyzed 35 scenes from the Land Remote Sensing Satellite (Landsat) Operational Land Imager (OLI) data published by the United States Geological Survey (USGS) to develop a land cover map of Japan that takes forest types into account.Each scene was processed with an Integrated Radiometric Correction [33] method to reduce effects from the atmosphere and topography, resulting in fewer errors throughout the classification process.Landsat OLI 5 bands (bands 2-6) along with the Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI) [34], and Green Red Vegetation Index (GRVI) [35] were used to perform the supervised classification using a maximum likelihood algorithm.This yielded a final land cover product with 12 classes.However, our interests here are in the forest classes only, so the number of classes was reorganized into three forest type classes (coniferous, deciduous broadleaf, and evergreen broadleaf) plus one additional class (others), resulting in a total of four classes.These were utilized for further procedures along with the accuracy assessment, which is fundamental for evaluating the resulting categorical map [36].The date of the acquired scenes included two from March, four from April, four from May, six from June, five from August, five from September and nine from October.Selection of the images was made by considering a general and critical issue for such optical images: that they were cloud-free.Not all scenes were completely cloud free.In such cases, selections were made from the clearest images that could be found.The specifications of the chosen data are listed in Table 1.
Table 1.Specification of the Landsat Imagery used.

Estimating CO2 Sequestration from the Forests of Japan
Estimation of the sequestration from Japanese forests was obtained by utilizing the two data sets described above (i.e., the volume map and land cover map), along with an additional data set of forest inventory data provided by two Prefectural Governments and the precipitation data published by the Ministry of Land, Infrastructure, Transport and Tourism (MLIT) of Japan.CO2 sequestration per unit area data published by the Tochigi Prefectural Government was used to estimate the sequestration values not only for each forest type but also for each tree age.This was achieved by transforming the volume map into ages among the different forest types, referring to the volume-tree age relationship data provided by the Oita Prefectural Government.The precipitation data were also used to help discriminate between growth rates for each region, which is an important step, because even with the same stem volume, the forests' per unit areas could have different tree ages due to local site productivity [37].
The stem volume-age relationship data provided by the Oita Prefectural Government is based on observations made on ground-level plots at multiple sites within Oita Prefecture.The curves were averaged for each growth region and computed as high, moderate and low growth rate areas.Each location within the study area that lies on a matching growth-rate region was determined to have a particular growth curve for that forest; the tree age conversion is based on that curve.To differentiate the growth rate of the study area, a relationship analysis was made using the precipitation data along with the total volume of growth made between the two periods of time.The sugi (cedar) and kunugi (sawtooth oak) data were used as representative of the coniferous and deciduous broadleaf forests, respectively.Unfortunately, there is a lack of data for the evergreen broadleaf forests, so the growth rate is considered only for the first two classes.Here, the precipitation data were used to determine whether precipitation would be one of the limiting factors for forest growth.The sequestration values were modified from the reference data published by the Tochigi Prefectural Government by taking a mean value of each tree type at each tree age range for deciduous and evergreen broadleaf forests, using the tree data for sawtooth oak, Japanese beech, Japanese zelkova and hornbeam for deciduous broadleaf forests, and using evergreen oak and chinquapin oak for the evergreen broadleaf forests.For the coniferous forests, a weighted mean value of the sugi (cedar) and hinoki (cypress) trees was taken as 7:3 because the distribution of those trees approximates that ratio for the whole of Japan.Together with the estimation, an in-depth analysis of the uncertainties of the estimated results was assessed by considering possible errors that could generate inaccuracies in the results.

Stem Volume Modeling and Mapping
Figure 3 shows the trends of the stem volume vs. the backscattering intensity (dB) for both HH and HV polarizations.An increasing trend can be observed with the increasing parameter and increasing backscattering intensity, which is consistent with many of the other related studies [19][20][21][22][23]30,31].According to these results, a second polynomial function from the HV polarization was selected and applied to obtain the stem volume information for the whole of Japan, because it has higher correlation compared to HH (HH: R 2 = 0.2123, HV: R 2 = 0.2198).The equation of the function is described as: where V is the stem volume and σ 0 HV is the backscattering intensity from the HV polarization.Validation of the backscattering model was assessed by using the remaining observation plot data that were not used in the modeling process.Using the plot data as ground truth information, the root mean square error (RMSE) between the estimated results and the reference data was calculated.The RMSE result is 105.58 m 3 /ha, and R 2 = 0.164 as shown in Figure 4.This model will be applied to obtain the stem volume information for the forested areas within Japan.Figure 5 shows the interpretation of the original and the corrected imagery of the PALSAR data.In addition to the interpretation, a quantitative evaluation of the performance of the correction was conductedby assessing the non-correlation between cos i and the corrected data.Comparing the correlation coefficients, the value shows r = 0.4 and r = 0.06 for non-corrected and corrected images, respectively.From these results, it is presumed that the correction of the PALSAR was performed successfully.The formula (Equation ( 4)) obtained through the model was applied to the corrected PALSAR imagery to obtain stem volume information for the entire forest area of Japan.By doing this, the spatial distribution of the estimated stand volume for Japanese forests can be shown.However, this model is specific to the coniferous stands in Japan, and the possibility that backscattering characteristics of different forest types would give different trends due to the structure of the stand and leaf size should be taken into account.Instead, for the broadleaf forests, a polynomial model developed by Wijaya [38] from the evergreen trees of Indonesia was applied to fit the situation (Equation ( 5)).Utilizing the two different models specific to those forest types, the stem volume distribution across Japan was obtained, as seen in Figure 6a,b.Using the model developed from the results of this study, the maximum stand volume of the coniferous forests reaches up to 1069 m 3 /ha; for the broadleaf forests, it reaches up to 414 m 3 /ha.Although the coniferous stand volume gives a high maximum value, most of the distributed areas lie near the 400-500 m 3 /ha range and rarely exceeds 700 m 3 /ha.The original stem volume information was distributed over all land cover types, so the exact location of each forest was extracted using the land cover map.The white areas in Figure 6 indicate non-coniferous (a) and non-broadleaf forests (b) areas. = 256.85− (65.458σ 0 ) − (6.8σ 0 2 ) (5)

Land Cover Map and Its Accuracy
The Landsat OLI images captured on multiple dates from years 2013 to 2014 were classified into the following 12 land cover types: water bodies, urban, paddy fields, crop lands, orchards, grasslands, silver grass, bare lands, coniferous forests, evergreen broadleaf forests, deciduous broadleaf forests and snow/clouds.These detailed classes were then reclassified into the three general forest classes and one other-than-forest class, resulting in a total of four land cover categories.Only the forest categories are of interest; the accuracy assessment was also performed with this map.The final product of the land cover map (before reclassification) of Japan is shown in Figure 7.The ground truth points obtained from Google Earth together with the reference vegetation map [39] were used in the resulting land cover map to perform an accuracy assessment by outputting an error matrix (Table 2).The final number of ground truth points yields 3550 points.However, the number of points was normalized for each class to 1000 points, to make calculations and interpretation easier in the accuracy assessment.The overall accuracy of the land cover (forest cover) map result was 73.93%, and the overall Kappa Index of Agreement (KIA) result was 0.65.KIA is a statistical measure adapted for accuracy assessment in the remote sensing fields by Congalton and Mead [40].KIA values range from 0 to 1, with zero representing completely random chance and 1 meaning a perfectly true agreement.In this study, KIA = 0.65, which could be said to show that the map is 65% likely to be the result of true agreement.The stem volume-age relationship data were modified from the original Oita data [5] to compute high, moderate and low growth curves for coniferous and deciduous broadleaf forests (Figure 8), which were used in turn to convert the so-obtained stem volume map (Figure 6) into tree age classes.Figure 9 shows the analysis between annual precipitation and tree growth for coniferous and deciduous broadleaf forests.
The results show a clear trend in the relationship between precipitation and tree growth.For the sugi trees (coniferous), increasing precipitation results in larger tree growth.On the other hand, for the sawtooth oak (deciduous broadleaf), the opposite trend can be observed, where tree growth increases with decreasing precipitation.The negative trends in the deciduous broadleaf forests might be caused by the precipitation data used.The precipitation data are an annual value, which means that it was not taken into account during the time when deciduous trees are likely to be photosynthesizing.If the precipitation falls mostly in the fall or winter, the trees may not be able to take advantage of it and thus would have lower than expected growth rates.The distribution of the precipitation classes and the thresholds is seen in Figure 10.For the coniferous forests, Class 1 is the area with a low growth rate, Class 2 is moderate, and Class 3 is a high growth rate.For deciduous broadleaf forests, Class 1 is high growth rate, Class 2 moderate, and Class 3 is low growth rate.Finally, Figure 11 shows a graph of the distribution of tree age overlaid with the forest cover data to extract the exact forest types.For the coniferous and deciduous broadleaf forests, the growth rates are considered for the age conversion, however for the evergreen broadleaf it was not, due to the lack of reference data, so that has been converted directly from the stem volume map using the stem volumeage relationship information provided by Sasaki et al. [41].The statistical distribution of forest tree ages can be seen in Figure 12.Coniferous forests can be widely distributed among the different tree ages.The majority lie in the 25-30 year age range for Site 1 (high growth rate regions), 30-35 years range for Site 2 (moderate growth rate regions) and 35-40 years range for Site 3 (low growth rate regions).We can also confirm some sizable areas over 75 years of age, which can be understood to be old mature forested areas left naturally and still growing without any management, while the areas with lower tree ages are places where plantations or forest management is in operation.However, future work is needed to validate the accuracy of this concept.Compared to coniferous forests, the majority of deciduous and evergreen broadleaf forests can be observed to have tree ages of over 50 years.These forests include few to no areas that are actively managed; they are just left to grow naturally, resulting in a majority of the areas being in highly matured age classes.

CO2 Sequestration by the Forests of Japan
Estimations of CO2 sequestration by forest type were obtained by multiplying the observed CO2 sequestration rates from the Prefectural data by the total area of each forest type's tree age category.Tables 3-5 show the total CO2 sequestration estimates for coniferous, deciduous broadleaf, and evergreen broadleaf forest age categories, respectively.
Figure 13 shows graphs of the sequestration values among each tree age for (a) coniferous and (b) broadleaf forests.Together with the calculated sequestration value, a mean sequestration value of each forest type from the reference by Tadaki and Hachiya [15] has been plotted.The coniferous forests show a peak value at approximately 15-20 years of age and then gradually decrease to a value similar to that of the 0-10 year range when they reach 75 years.Each site dependency for the growth rate clearly shows a difference in the total sequestration amount, where compared with trees in low growth rate areas, high growth rate area trees have approximately a 10 t•CO2/ha/yr difference at the peak value.For the broadleaf forests in this study, the results do not show a large difference between evergreen and deciduous broadleaf forests.The peak value can be observed at approximately 15-20 years of age, which gradually decreases, similar to the coniferous forests; however, the total sequestration value is much smaller than for coniferous forests.It is obvious that Japan has operated a large-scale plantation operation using sugi and hinoki trees because of their fast growth characteristics, which these values reflect.By using the detailed sequestration values for tree ages, sequestration overestimates that implicitly occur with the estimation method using averaged values [15] for coniferous trees of less than 10 years of age and for trees of more than 20 years of age, as well as the underestimation for the peak tree ages (15-20 years), can be avoided altogether.Similarly, for deciduous and evergreen broadleaf forests, compared to the averaged values for all ages, the values for the tree age categories show lower results.
Finally, Figure 14 shows the overall output of the final work, visually illustrating CO2 sequestration values for forests in Japan.

CO2 Estimates and Potential Errors
The total CO2 from the forests of Japan has been calculated by taking into consideration both forest type and tree age.It should be understood that determining the accuracy of this result is quite difficult, because there is no known "true value" for this information.It is also necessary to note that there are always sources of error when attempting to calculate such estimations.Possible errors that could generate inaccuracies in the results have been identified and listed here: 1. Accuracy of the land cover map (forest cover map) 2. Dependency on the backscattering model for stem volume estimation 3. Stem volume-age relationship curve 4. The sequestration value itself Number 1 refers simply to errors resulting from the misclassification of each forest type.As acknowledged in the work, the total area size of a forest type is an important factor in calculating total CO2 sequestration, because the sequestration value per unit area is multiplied by the total area of each forest type.Figure 15 illustrates the relationship between accuracy and errors for each forest type in reference to the error matrix (Table 2).For example, coniferous forests may have had an 85% accuracy rate in classification, while 4% were misclassified into evergreen broadleaf, and 6% were misclassified into deciduous broadleaf forests.The same can be said for the other two classes.Thus, an 85% accuracy rate for coniferous may be true, but there are 4% and 6% chances of errors due to misclassification into evergreen and deciduous broadleaf forests, respectively, in calculating the total CO2 amount.The same can be shown for the other classes, mutatis mutandis.The evergreen broadleaf class shows 20% and 25% chances of error from deciduous broadleaf and coniferous forests, respectively, while deciduous broadleaf shows 11% and 15% chances of error from evergreen broadleaf and coniferous forests, respectively.
For number 2, errors in the backscattering model can be explained by the RMSE calculated from the estimated stem volume versus the ground truth volume information (Figure 4).As has been explained in an earlier section, the RMSE of the backscattering model resulted in 105.58 m 3 /ha.This can lead to errors in transforming the stem volume into the tree age information, and thus become an error in the estimation of CO2 sequestration.For example, in the coniferous trees in high growth rate areas, a 100 m 3 /ha difference can miss the true tree age by ±6 years.Larger differences accrue to lower growth rate areas, for example at Site 3 (low growth rate region) 100 m 3 /ha can miss the tree age by ±13 years by looking at the stem volume-age relationship data from the Oita Prefectural Government.A 6-year difference might change the estimation values; for example, the forested coniferous areas at peak ages could be classified to earlier or older age ranges and vice versa: the older age classes could change to the peak tree age classes.Therefore, this error source must be acknowledged as one possible error in the estimation.
For number 3, the tree age transformation used the reference data provided by the Oita Prefectural Government and from Sasaki et al. [41].According to this growth curve, the stem volume information was reclassified for the whole of Japan.The issue here is whether this growth curve also has implications for other regions in Japan.This study used the reference information as a representative growth curve for each tree type, and using the precipitation data, it has been assumed that the other regions of Japan would similarly lie along one of those growth curves.To make sure, additional references were collected to compare the curve information to other regions of Japan [42,43].According to the tree growth curve information from the Kagoshima region (Kyushu district) and Aomori region (Tohoku district), there are differences in the growth curves depending on the different site productivity areas; however, comparing the two figures, a large difference between growth curves even in opposite environmental regions (the warmer Kagoshima region, and the cooler Aomori region) cannot be confirmed [42,43].It is known that there is a most suitable growth curve for each Prefecture in Japan, but this study considered that by making an attempt to normalize the differences of the growth curves by using the precipitation data and setting the threshold.Thus it can be said that the growth curve differences among each region are averaged out and reduced to the least possible error in this aspect.Still, this is one point that could be considered as a potential error in the estimation.
Finally, for number 4, the sequestration information was modified from the Tochigi Prefectural Government by taking the mean or weighted mean of each tree type and using that value as representative of each forest type tree age class.The Prefectural data for the CO2 sequestration is estimated using the guidelines commonly used for estimations in the forests of Japan, based on the following equation: Annual growth is measured by stem volume growth.The biomass expansion factor is used for expanding the stem volume into an estimate for the whole volume of the tree, including the volume of leaves, branches and roots together.This value usually differs among tree types and species.Bulk Density represents the weight per cubic meter of volume.Carbon content commonly uses the value of 0.5, which means that the wood is composed of approximately 50% carbon.Finally, to convert the carbon value to a CO2 sequestration value, the CO2 Equivalent Factor of 44/12 is multiplied, where 12 and 44 are the molecular masses of carbon and CO2, respectively.
An error that can occur from this is the sequestration value used in multiplying the total forest area, which would be different if it were actually broken down in further detail for each specific tree type and age.Some trees capture more or less CO2 compared to other tree types within the same forest type, but those values were averaged for this study.Another possible error source is the minimum and maximum sequestration value for each tree age class.Because it is difficult to discriminate between forests using a one-year age distribution, the tree ages were reclassified into 5-year interval ranges.This study used the mean value within that range, and there are minimum and maximum values for those categories, and those can cause errors in the sequestration estimation as well.Table 6 shows the minimum and maximum sequestration values for each forest type at each tree age class range.As shown, there are some classes that have high sequestration values at the peak age of coniferous trees (10-15 years) at site 1, which is both a little higher than the average value from the evergreen broadleaf forests of Tadaki and Hachiya [15] and also similar to the deciduous broadleaf forests.Looking at all the ages, the biggest differences between maximum and minimum can be observed in the peak aged forest classes.This point is one of the more probable errors.The CO2 sequestration amount by the forests in Japan was calculated after each error was taken into consideration as discussed in Section 4.4.1.Table 7 shows the summary of the final output of the total CO2 amounts when these errors are considered.
In Table 7, the summary of the total CO2 calculated from this study is listed when each of the error cases is considered.Four different cases of potential errors were examined: land cover misclassification, reliability of the backscattering model, dependency of the tree age curve, and the sequestration value.The tree age curve is listed as N/A because it is simply the errors caused by the growth trend of the trees, and giving an estimation by using other Prefectural Governments' tree age curve information will give results similar to those obtained using the Oita data.However, for the errors from land cover mapping misclassification, backscattering model and sequestration values, and after considering the backscattering model and sequestration value errors together, the estimation was calculated.
In the table, Avgval Minval and Maxval represent the average, minimum and maximum sequestration value used for the estimation, respectively.Mintot and Maxtot represent the minimum total and maximum total CO2 sequestration by each category, respectively.Thus, for backscattering model error cases in the coniferous category, the total CO2 amount was estimated using the average sequestration value (the value listed in Table 3).Values in brackets written after N/A are the total sequestration outputs listed in the tables of deciduous and evergreen broadleaf forests (Tables 4 and 5).Errors of misclassification were considered by using the percentage of the total area of each distributed forest class modified based on the actual misclassification percentage shown in Figure 15.For example, 4% and 6% of coniferous classes were misclassified to evergreen broadleaf and deciduous broadleaf forests, respectively.Thus, 4% and 6% of the total area from those forest classes, respectively, were subtracted and added to the coniferous class (this is determined as the maximum status).At that point, the deciduous and evergreen broadleaf status was determined as the minimum status.On the other hand, 25% and 15% were misclassified to the coniferous class when they were likely to be evergreen and deciduous broadleaf forests, respectively.Similarly, those percentages of the total area were subtracted from the coniferous forests and added to the evergreen and deciduous broadleaf forest classes (this is determined as the minimum status).Again, at this point, the maximum deciduous and evergreen broadleaf statuses were determined.The dependency of the classification results to errors in estimates of sequestration were assumed in this method.As a result, the minimum total sequestration was 75.80 Mt•CO2, while the maximum was 127.11Mt•CO2.Here, the sequestration value by tree age class for each forest type does not change, and only the total area of the forests fluctuates.We can see how the total size of the forested area is one important factor that contributes to total sequestration.
Errors of the backscattering model were considered from the developed RMSE model, which resulted in 105.58 m 3 /ha.According to the stem volume-tree age curve information, the difference of stem volume results in 6-year, 7-year and 10-year differences for the coniferous trees on site 1, site 2 and site 3, respectively.Using this information, the estimation of the CO2 in coniferous forests was re-assessed by shifting the age range by the above tree age year difference.The minimum and maximum then became 65.46 and 99.81 Mt•CO2, respectively.The total CO2 was then calculated using these values along with the total values of other forest types (deciduous broadleaf (DBL) and evergreen broadleaf (EBL)).The estimation was made only for the coniferous forests in this case because the actual error values of the model developed by Wijaya [38] were not found.
For the sequestration value errors listed in Table 6, the minimum and maximum values were applied, multiplying by the total forested area of each forest type.Overall total CO2 resulted in 92.73 Mt•CO2 and 137.69 Mt•CO2 for minimum and maximum totals, respectively.This case is easier to interpret compared to the other errors because all the forest classes can be estimated clearly.
Finally both the errors from the backscattering model and the sequestration value were added, and then the estimation was re-assessed.This was conducted only for the coniferous class, and resulted in 64.65 Mt•CO2 and 91.45 Mt•CO2 for minimum and maximum amounts, respectively, when the minimum sequestration value (Table 6) was used, and 75.03Mt•CO2 and 110.64 Mt•CO2 for minimum and maximum amounts, respectively, when the maximum sequestration value (Table 6) was used.

Stem Volume Empirical Modeling
In Section 4.1, an empirical backscattering model was developed to estimate stem volume information across the forests of Japan.The challenge is that backscattering can have multiple causes within a forested area, and can be affected by species, structures, soil moisture [44], density of forests [18], etc.A simple empirical modeling method was used by taking a correlation analysis between the backscattering intensity and the ground reference data.We have avoided most of the possible distracting factors affecting the difference in the backscattering to make a more accurate trend between the parameters and SAR image, which allowed a better interpretation when compared to many studies performed in similarly rugged areas [24].We could discuss about the limitation of L-band SAR for the estimation of biomass, due to its sensitivity which gives saturation at a certain biomass level [45], however the most important variable in our study is the topographic effect.Although various radiometric corrections have been developed to normalize the SAR imagery [19,26,27,46] the radiometric accuracies of those correction methods can always be challenged.However, without their application, the topographic effects caused by rugged terrain will remain a crucial issue, especially for the study area here in Japan, and can result in higher errors compared to not applying them.There might be solutions to improve the interpretation of the modeling, but in this first study on national-scale estimation using remote sensing, we believe that the result can serve as a standard.Note that during the relationship analysis, speckle filtering was not applied to the PGM data.This decision was made because the data were downscaled to 25 m resolution from the original, which averaged out the area, reducing the noise.Further filtering might lose the actual information from the surface objects.However, in the future, a different filtering method might be applied to evaluate whether the difference affects the results.

Land Cover Mapping of Japan
A land cover map was produced that not only defined the areas as simple forests and non-forested areas but also provided more detail (i.e., subclasses) such as forest types (coniferous, deciduous and evergreen broadleaf), so the information could be used in many different ways.The need for a detailed digital land cover (or digital vegetation model) map, especially for the vegetation categories, is important for obtaining more precise estimations of CO2 sequestration values.For this study, only single-scene data sets were used to develop the land cover map for each area.Sanga-Ngoie et al. [5] states the usefulness of applying temporal data sets to improve land cover accuracy due to the clear difference of radiance values that change during different seasons of the year.In the case of forest classes, similar spectral signatures would be observed throughout a season after blooming finishes, but clear and distinguishable differences can be observed for different forest types, such as trees whose leaves fall in the winter.Such methods were not applied in this case because few cloud-free products were captured, which made it difficult to find a set of scenes with two different seasons for each location.Operation of Landsat 8 is still in its third year, so its scene database is still relatively small compared to Landsat 5 or 7, so such methods are currently limited to very small areas or must wait for more data to be collected over time.
However, because of the lack of temporal data scenes, the classification process was performed not by increasing the number of images but by collecting as many detailed training sites as possible.This method led to the ability to discriminate between forest classes in more detail and with more precision, because constructing more detailed and precise land cover classes resulted in greater chances of unclear pixels being classified into the forest class to which they should correctly belong.For example, if the training sites of the forest classes are specified and classified only as "coniferous, deciduous broadleaf and evergreen broadleaf", deciduous and evergreen broadleaf forest classification drop inaccuracy by 33% and 36%, respectively, compared to when more detailed training classes of those forest types are constructed.Therefore, this methodology is effective when there are few scenes available for classifying a region.Constructing detailed training sites can improve the categorical map.In the future, we recommend reprocessing the classification of areas that had lower accuracy results or regions that were covered by snow/cloud to improve and correct the land cover types used here.

CO2 Sequestration Values and Trends
The conventional method [15], when used alone for estimating CO2 sequestration, results in large errors in the total amount compared to methods that give special attention to tree age.The stem volume method for discriminating tree ages within each forest type clearly shows the importance of this variable for avoiding over-or underestimations, depending on the growth stage of the forests examined.The use of L-band SAR data rather than optical data [5] made this discrimination possible, yielding sequestration rates by tree age for trees up to 75 years of age for coniferous and 50 years of age for deciduous and evergreen broadleaf forests.It is thus clearly established that sequestration rates of 13.5 t•CO2/ha/yr [15] are only valid for 10-15-year-old coniferous forests in high growth rate regions.In contrast, deciduous broadleaf forests reach just below the 8.7-t•CO2/ha/yr threshold only in the 15-20-year-old range.For evergreen broadleaf forests, this is a large issue when the conventional method is applied, because it shows approximately 3 times greater value than observed in the results of this study, meaning that there would be a large overestimation of the forest carbon sink capacity if such a value is applied.High sequestration values from evergreen broadleaf forests are difficult to confirm for many reasons; however, the amount of sequestration by evergreen trees in Japan is acknowledged to depend on total solar radiation and air temperature.Gross Primary Product (GPP) results from the total radiation, while the Net Primary Product (NPP) is set by the amount of respiration subtracted from the GPP [47].In this case, higher air temperature leads to higher respiration and decreases the NPP.On the other hand, when both total radiation and air temperature are low, respiration rates decrease but NPP becomes the same or even greater compared to higher radiation periods [37].Depending on the season and location of the observation, the average value of evergreen broadleaf forests will differ, and it might well be possible that the large value obtained by Tadaki and Hachiya [15] resulted from such causes.
It is very important to note that the averaged sequestration value for the coniferous forest type matches the sequestration value by tree ages obtained here for the 10-15-year-old coniferous trees at Site 1.This is also the category of the second fastest-growing period for coniferous trees at that site.The same is also true for deciduous broadleaf forests for 15-20-year-old trees; it is just a little lower than the average value.This strongly implies that an estimation based on the averaged values alone will round up the values for trees in different growth stages, causing an overall overestimation of the carbon sink capacity [5], both for natural and planted forests.The results emphasize the fact that sequestration differs among tree ages.This reality should be taken into consideration in all cases, otherwise, erroneous values of the forests sequestration capacity will be obtained, generating results that possibly have errors from using less accurate values.
The uncertainty assessment gave a clear view of how and why the estimations can fluctuate, where the estimated total CO2 lies between the minimum 75.80 Mt•CO2 and maximum 137.69 Mt•CO2.The minimum estimate results from land cover error and the maximum from sequestration value error.After examining the results presented in this paper, it can be assumed that these two parameters are the most important factors when estimating sequestration amounts.Focusing on the results where 111.27Mt•CO2 was shown, the fluctuation of the estimated amount is a difference of approximately 23 to 32%.Certainly, if the errors for the DBL and EBL classes for the backscattering model errors could be obtained, some change might be observed in the results; however, it may be determined that even if this occurred, it would not result in a very large difference from the current numbers.Through the uncertainty assessment, a better understanding of the fluctuation of the CO2 trend became clear, and thus this result leads to a strong confidence in the practical usage of this methodology.

Validity of the Method in Comparison
For further knowledge, other estimation methods were investigated and compared to the results of this study.Four main methods were examined and compared; one is from this study, the second is from the conventional method by Tadaki and Hachiya [15], the third is the work conducted by Sasaki and Kim [7], and the last is from the National Institute for Environmental Studies (NIES) [8] of the IPCC report for the forest carbon sink.The total CO2 estimation is summarized in Table 8 [7,8,15] for all four methods.For the conventional method, the sequestration value was used along with the total forest areas obtained from our work.[15] using the averaged value compared to the results in this study, a large overestimation of the sequestration can be observed, because the sequestration value dumps everything into one category, omitting tree age information and ignoring lower sequestration stages, leading to high overestimation of the result.Sasaki and Kim [7] and NIES [8] show similar estimation values, however for the former study, the result is obtained by using the land use model and carbon stock growth model and estimates the sequestration value classifying the forests as planted and natural forests.Planted forests showed an estimation of 56.1 Mt•CO2, while natural forests were estimated at 17.6 Mt•CO2, and then summed for a total result of 73.7 Mt•CO2.The most interesting result is from the NIES, whose approach used forest registration data from each prefectural government and differentiated the forests by tree type, collecting annual growth information for each, and then performing an estimation by multiplying the annual sequestration value by the total forested area of each tree type, which is similar to this work.They have not only made calculations for aboveground sequestration but also for the litter and soil.However, the majority of the sequestration results from the aboveground biomass (aboveground: 76.53 Mt•CO2; litter: 0.41 Mt•CO2; soil: 2.58 Mt•CO2; dead tree: −1.85 Mt•CO2).If we assume that the methods used by the NIES [8] study are likely to be the most accurate estimation of the natural environment, the remote sensing method applied in this study shows only a 43% overestimation in comparison.This difference can be explained from the errors discussed in the previous section.If methods are developed to reduce those potential errors, our results may become somewhat similar to the results by the NIES.However, even though the NIES results raise some questions, such as the uncertainty of the total forest area, which is important as mentioned, total CO2 is calculated by multiplying the sequestration value by the total forested areas.The data from prefectural government is sometimes unclear, especially in terms of the area that the forest covers, as researchers have noted (the uncertainty rate of NIES estimates shows 5% for forest areas without discriminating by forest type [8]), so there is a possibility that the NIES is underestimating forest capacity.What is significant about this work compared to the NIES is that this method can be performed throughout various seasons in the future in less time and with lower labor costs through the use of remote sensing data.The work performed by NIES is possible in countries where a large database of forest information is compiled for use in applications, but in countries where this is not the case, it will be very difficult to implement such a method.On the other hand, the method developed in this study can be applied in any country as long as the required data sets are available.Moreover, the uncertainty of the total forest area compared to NIES can be overcome using our methodology, because the correct amounts and locations of forests are easier to discern.Overall, this study explains why the result is much lower (better) than the conventional method, which seems to result in high overestimation, while when compared to NIES results, the estimation arrived at by this study does not show serious underestimation.The true result likely falls between the two as discussed in this paper.The results achieved using the stem volume method in this study show both an improvement and the potential for obtaining even better results than other methods by using remote sensing method.As a note, the term sequestration used in this study can be described as NPP, while the NIES result can be described as Net Ecosystem Product (NEP) because it also considers decomposers.In the future, improvements in the methodology will likely take this into consideration and obtain NEP values of the forests so we can estimate a more precise value not only for use in forest management and policy but also for use in modeling the biogeochemical cycle of terrestrial ecosystems.

Concluding Remarks
Rapid change in climatic conditions in recent years has raised many pressing issues and problems that need to be solved.These are of global interest among researchers, companies, politicians and others.Therefore, in response, various measures are being taken, such as the Kyoto Protocol, REDD+ and others, to mitigate the impacts of climate change.Although such policies and actions are being implemented around the world, the quantitative estimation of how well they have contributed to the natural environment, such as in terms of carbon sequestration, is still under discussion.Meanwhile, this work has shown a path forward in quantitative estimates on a national scale for CO2 sequestration by utilizing fine resolution satellite imagery through a local/regional scale analysis.Focusing on forest type and tree age information using the LANDSAT and PALSAR data, respectively, yields more accurate data.This methodology gives better interpretations and more precise estimations of CO2 sequestration, as it is clear that without both tree type and tree age discrimination, overestimations may well result.
The CO2 amount and other important information revealed in this study has provided important data.Do old mature trees sequestrate as much as younger trees?The answer is no when we see the trend of the sequestration as a function of tree age.Using this information, we can discuss policies such as the REDD+ procedure.Is planting trees for 100 years' time [48] for a credit really an effective and correct way to manage forest resources?Rather, it could be suggested that planting faster growing trees in a period of 20 to 30 years and then replacing them continuously might be more sustainable and effective, even in economic terms.The findings in this paper can thus suggest these types of ideas to policy makers for more sustainable development/conservation.
The methodology is not only of interest for forestry, stakeholders, and policy makers but also for researchers who address modeling of the earth's ecosystems.The knowledge of CO2 trends in the terrestrial biosphere developed from this work contributes to a better understanding of the biogeochemical cycle and thus sheds more light on the influence of the terrestrial ecosystem and its interaction with the atmosphere.This could lead to a step forward in understanding more about climatic systems.Utilizing fine resolution satellite imagery on a more local/regional scale can be expected to improve modeling accuracy and help better understand the trends of the earth's systems.
The in-depth quantification of forest sequestration can play an important role in many fields, helping them to improve how they set strategic actions into operation.Promoting such methods will certainly support further development and discussion to make a difference in solving or mitigating the issues faced in a more effective way.

Figure 1 .
Figure 1.Study Area, Japan and its eight major districts.

Figure 2 .
Figure 2. Overall flowchart of the methodology.

Figure 4 .Figure 5 .
Figure 4. Validation of the backscattering model.The X-axis is the estimated volume (m 3 /ha) using the model, and the Y-axis is the reference data from the observation plot.R 2 = 0.164 Root Mean Square Error (RMSE) = 105.58m 3 /ha.

Figure 6 .
Figure 6.Stem volume map of Japan for the (a) coniferous forests and (b) deciduous and evergreen broadleaf forests.Forest cover data were overlaid with the stem volume map to extract only the forest types needed for each stem volume map.

Figure 7 .
Figure 7. Land Cover map of Japan developed using 35 Landsat Operational Land Imager (OLI) scenes yielding 12 classes in total.

Figure 8 .Figure 9 .Figure 10 .
Figure 8. Mean growth curve of (a) coniferous and (b) deciduous broadleaf forests in each growth rate region, computed from the original Oita data.

Figure 11 .
Figure 11.Tree age map transformed from the stem volume map for the (a) coniferous, (b) deciduous broadleaf and (c) evergreen broadleaf forests.

Figure 12 .
Figure 12.Statistical distribution of forest area by tree age for (a) coniferous and (b) deciduous/evergreen broadleaf forests in Japan.

Figure 13 .
Figure 13.Trend of CO2 values by tree age classes for (a) coniferous and (b) deciduous/evergreen broadleaf forests.

Figure 15 .
Figure 15.Triangular relationship of the forest classes and their accuracy and errors.

Table 2 .
Error Matrix showing the reference data (ground truth) versus the image classification.Overall accuracy shows approximately 74% and Kappa Index of Agreement (KIA) as 0.65.

Table 3 .
CO2 sequestration of coniferous forests by tree age using the PALSAR data and the stem volume method.

Table 4 .
CO2 sequestration of deciduous broadleaf forests by tree age using the PALSAR data and the stem volume method.

Table 5 .
CO2 sequestration of evergreen broadleaf forests by tree age using the PALSAR data and the stem volume method.

Table 6 .
CO2 sequestration minimum and maximum values for each forest type and tree age class.

Table 7 .
Summary of total CO2 sequestration estimated after considering each possible error source.

Table 8 .
Summary of the total CO2 sequestration estimated by various methodologies in the different studies.The stem volume method used in this study, shows a result of 111.27 Mt•CO2 by considering both forest type and tree age information (coniferous: 85.0 Mt•CO2; deciduous broadleaf: 21.61 Mt•CO2; evergreen broadleaf: 4.76 Mt•CO2).With the conventional method