Quantifying Impacts of National-Scale Afforestation on Carbon Budgets in South Korea from 1961 to 2014

Abstract

Forests play an important role in regulating the carbon (C) cycle. The main objective of this study was to quantify the effects of South Korean national reforestation programs on carbon budgets. We estimated the changes in C stocks and annual C sequestration in the years 1961–2014 using Korea-specific models, a forest cover map (FCM), national forest inventory (NFI) data, and climate data. Furthermore, we examined the differences in C budgets between Cool forests (forests at elevations above 700 m) and forests in lower-altitude areas. Simulations including the effects of climate conditions on forest dynamics showed that the C stocks of the total forest area increased from 6.65 Tg C in 1961 to 476.21 Tg C in 2014. The model developed here showed a high degree of spatiotemporal reliability. The mean C stocks of the Cool forests and other forests increased from 4.03 and 0.43 Mg C ha⁻¹, respectively, to 102.43 and 73.76 Mg C ha⁻¹ at a rate of 1.82 and 1.36 Mg C ha⁻¹ yr⁻¹ during the same period. These results imply that, although the total Cool forest area of South Korea occupied only about 12.3% (772,788 ha) of the total forest area, the Cool forests play important roles in C balances and forest ecosystems in South Korea. Annual C sequestration totals are projected to decrease at a low rate in the near future because the overall growth rate of a mature forest decreases as the stand ages. Our results quantified forest C dynamics in South Korean forests before and after national reforestation programs. Furthermore, our results can help in development of regional and national forest management strategies to allow for sustainable development of society and to cope with climate change in South Korea.

1. Introduction

The success of climate change mitigation and adaptation and implementing sustainable development goals depends on a diverse array of ecosystem services provided by forests. Over the coming century, these and other forest-derived ecosystem services will need to be sustained to meet the increasing demands of a growing human population [1]. The difficulty of this task is exacerbated by the impacts of anthropogenic climate change [2]. In addition, deforestation is a major driver of both climate change and biodiversity loss [3]. The various worldwide problems and damage caused by deforestation and climate change have been well-documented. Therefore, many countries, including India [4], China [5], New Zealand [6], Finland [7], the UK [8], and South Korea [9], have formulated and implemented a number of policies and programs aimed at forest conservation, afforestation, and reforestation.

South Korea is a rare, exemplary model of full-scale reforestation [10]. After the Korean War, almost half of its forestland was destroyed. Since 1973, following this period of serious deforestation, the South Korean government has implemented national plantation programs to promote forest recovery. After about 40 years of these efforts, South Korean forests have successfully recovered, and the stocking stem volume increased from 8.2 m³ ha⁻¹ in 1954 to 150.0 m³ ha⁻¹ in 2015 [11]. Furthermore, this success has spread to several different sectors with positive outcomes in the realms of land restoration, flood prevention, biodiversity recovery, and increases in water supply and recreational forests [12,13,14,15]. One of the major management issues that the forest sector in South Korea is faced with is understanding the impact of climate change and optimizing the use of forest resources [16]. There are numerous impacts of climate change associated with social, environmental, and economic problems, and a number of countries are creating mitigation and adaptation plans for climate change [17]. Carbon storage and sequestration are essential attributes of forests and important components of mitigation and adaptation strategies. However, this in turn indicates that forests are also subject to impacts of climate change [18,19].

Studies on forest carbon (C) stocks and balances in South Korea over the past few decades have been conducted for many years. Previous research based on the national forest inventory (NFI) and theoretical models has shown that the biomass C stocks of South Korean forests increased gradually in recent decades [9,15,20,21,22,23,24]. Lee et al. [24] estimated the forest carbon storage in the period of 1990–1997 as a part of estimating the changes in annual CO₂ fluxes in Korea. Choi et al. [15] analyzed the temporal variations in forest biomass, forest carbon density, and forest carbon uptake rate in the period of 1954–2000 and reported that Korean forests had a higher carbon uptake rate (Mg C ha⁻¹ yr⁻¹) in the late 1990s than in earlier years. It was concluded that this high carbon uptake rate resulted from the high growth rate of young trees present as a consequence of the 30 year reforestation and forest management program. However, there were some limitations to that study, such as (1) a relatively short simulation period of not more than two decades, (2) the large spatial scale of the simulation units, (3) insufficient spatial validation, and (4) the lack of an annual climate impact on forest stand dynamics.

In addition, Cool forests, defined at the 2018 International Boreal Forest Research Association Conference as forests at elevations above 700 m, have recently needed more careful monitoring and assessment than other forests. This is because Cool forests have been assessed as highly sensitive and vulnerable to climate change in South Korea and worldwide [18,25,26,27,28]. Therefore, comparison of C dynamics between the Cool forests and other forests through systematic observation and spatial-temporal modeling is required.

The primary objective of this study was to estimate C stocks and their changes in South Korean forests, including climate impacts on forest growth during the post-war period (1961–2014) using a Korea-specific dynamic model. To estimate the spatial distribution of forest C stocks, we constructed a spatial data set for Korean forests using the fifth NFI and the fifth Forest Cover Map (FCM). The model was validated by comparing estimated stem volumes with the observational data in the fifth and the sixth NFI on multiple spatial scales. Finally, we compared C stocks and annual C sequestration of Cool forests with those of other forests in South Korea.

2. Materials and Methods

2.1. Spatial Data for South Korean Forests

2.1.1. Study Area

South Korea is a peninsula located in mid-latitude East Asia and has a temperate climate. Annual average precipitation is 1000–1800 mm, and 50–60% of the precipitation is concentrated in the summer (June to August) (Figure 1a). The average annual temperature is 10–15 °C. More than 60% of the country is mountainous, and the terrain is high in the east and low in the west (Figure 1b). Most of the forests in South Korea are located in mountainous areas. Forests cover about 63.7% (6,369,000 ha) of the total land area of South Korea, and evergreen forests (mainly Pinus densiflora), deciduous broad-leaved forests (mainly Quercus spp.), and mixed forests made up approximately 40.5%, 27%, and 29.3%, respectively, of the total forest area in 2017 [11].

Cool forests were defined at the 2018 International Boreal Forest Research Association (IBFRA) Conference. Boreal and mountain forests develop over many thousands of years in regions with cold climates. They make up more than one-third of global forests, forming the largest terrestrial vegetation ecosystem. They are found from the circumpolar belt in the northern hemisphere to high-elevation forests spread over the entire planet. Forests on permafrost show many similarities with those in boreal and high mountain ecozones—especially with respect to species and growth patterns and in response to climate exposure. The 2018 IBFRA defined a Cool forest in the temperate zone as a forest at an elevation above 700 m. When this definition is applied, the total Cool forest area in South Korea is an estimated 772,788 ha based on high spatial resolution (30 m × 30 m) digital elevation model data (Figure 1c).

The Korea Meteorological Administration provided climate data, including monthly mean temperature, mean daily minimum temperature, and total precipitation, from 75 weather stations across South Korea from 1961 to 2014. The data were interpolated with a 0.01° grid size (≈1 km) using kriging and inverse distance weighting, considering the absolute temperature and precipitation lapse rate by elevation [29]. These datasets were resampled into a 0.01° spatial resolution. Figure 1d shows the annual mean temperature and the annual precipitation of South Korean forests area in 1961–2014. The annual mean temperatures of Cool and other forest areas rose by 1.53 °C and 0.78 °C, respectively, in this time. In addition, the annual precipitation of Cool and other forest areas rose by 362.7 and 299.4 mm, respectively, during the same period. These data imply that there was tendency for Cool forest areas to warm more quickly than the other forest areas during the last five decades.

2.1.2. National Forest Inventory Data and Forest Cover Map

The fifth NFI was conducted for the entirety of South Korean forests from 2006 to 2010 [30]. The survey design consisted of systematic sampling at intervals of 4 km (longitude) × 4 km (latitude) across South Korea (Figure 2a). Four circular sample plots were located at the intersection of each 4 km × 4 km grid line. Each sample plot (16.0 m radius) covered 0.08 ha. Forest characteristics [tree species, age, height, diameter at breast height (DBH), site index, and number of trees] and topographical factors (coordinates, elevation, slope, and aspect) were measured at all sites [30]. The sixth NFI was conducted from 2011 to 2015. Each permanent plot was surveyed five years after its previous survey.

A forest cover map (FCM) was used to predict the present stem volume. This map (at a scale of 1:5000) was produced from visual interpretation of aerial photographs and NFI data, and it outlined forest stands classified by tree species, DBH, age class, and canopy closure [31]. Seventeen categories of tree species were specified. The DBH and the age of stands were expressed as “DBH class,” (6–16 cm, 16–30 cm, and over 30 cm), and 10 year “age class” ranges.

2.2. Investment in the National Forestation Program

2.2.1. Establishing a Time Series of Forest Stock Maps

Information about the changes in forest areas in the Republic of Korea during 1961–2014 was provided only as statistical data, and there are no official spatial data. In addition, the forest area in the land cover data from the 1960s to 2014 is similar to that in the forest statistics. In other words, the change in the forest stock is not reflected in the statistics, as only the total forest area is provided. For this reason, in this study, the concept of a forest cover map was applied [9]. The ages of forest areas were estimated using a different method.

Spatial data for the time series of forest areas stocked though the national forestation program in the five decades of concern was produced using the FCM and the fifth NFI. We determined stand age information in the 2000s using these two date sources. In the FCM, stand age was assumed to be 5, 15, 25, 35, 45, and over 55 for age classes I, II, III, IV, V, and VI, respectively (Figure 2c). First, the representative age class was determined as the age class occupying the largest area of each grid. Second, for grids that represented age classes I, II, III, or IV, the stand age of each grid cell was calculated as an area-weighted average value of each stand within the grid cell. Third, the stand ages of grids that were representative of age class V were determined within a range from 51 to 107 based on the statistics of the fifth NFI. The larger the area of age class V in each grid was, the older the stand age was assumed to be. The probability distribution of stand age for grids that represented age class V followed the distribution of stand age in the fifth NFI data (Figure 2e). The probability distribution of estimated stand age showed a trend similar to that of observed stand age in the NFI data on a national scale. As the results were similar, it could be inferred that the model successfully reflected the trends in stem volume of Korean forests, further suggesting successful reconstruction of age distribution at a national scale, although uncertainties remained for individual stands.

2.2.2. Preparation of Inventory Data on Past and Present Forest Status

Four independent Korea-specific models, including height growth, tree mortality, DBH growth, and growth modifier models, were used to simulate annual stem volume and carbon stock changes. The connection among these models did not necessitate complex methods. The overall process used to estimate stem volume through the four different models in this study is summarized as follows: (1) estimation of growth of stand mean height and dominant tree height based on age and site index of each stand using regression models, as suggested by the Korea Forest Service (KFS) [32]; (2) estimation of stand density (trees per ha) change using a stand-level mortality model [33]; (3) estimation of annual stand mean DBH growth reflecting topographic and climatic conditions using regression models and a growth modifier developed by Piao et al. [34] and Kim et al. [35]; and (4) using the information from the previous three steps to estimate forest volume and carbon stocks using regression models, basic wood density (BWD), and biomass extension factors (BEFs) developed by the National Institute for Forest Science (NIFoS) [36].

SI is defined by the KFS as the height of the dominant tree species at 30 years of age. Therefore, to estimate the dominant tree height (ho) of each stand, we used a functional formula (Equation (1)):

(1)$h{o}_i=SI·{\left(\frac{1-e^{a·A}}{1-e^{a·30}}\right)}^b$

where ho is dominant tree height (m), SI is site index, A is stand age (years), i is the identification number of the stand, and a and b are species–specific coefficients. Coefficients proposed by the KFS [31] are shown in Table S1 in Supplementary Material. Regression equations and coefficients developed by the KFS [32] for the tree species that dominated each stand were applied in Equation (2).

hm = f(ho)(2)

where hm is mean tree height (m). Dominant tree height (ho) is one of the most commonly used indicators of site productivity because there is a close correlation between stem volume and site index, and it is generally accepted that ho is minimally affected by competition. A Korea-specific self-thinning model was used to estimate change in stand density [33]. This model was derived from Sterba’s theory based on the competition density (C-D) effect and dominant tree height; it showed a high degree of reliability (R² = 0.55–0.81) and no obvious dependencies or patterns in residuals for major Korean tree species. Coefficients for each tree species are shown in Table S1 in Supplementary Material:

(3)$\Delta S{N}_i=a·e^{b·\frac{S{N}_i}{MS{N}_i}}·\left(MS{N}_i-MS{N}_{i+1}\right)$

where i is stand age (years), SN_i is stand density at time i, MSN_i is the maximum stem number at time i, ΔSN is the number of dying trees from i to i + 1, and a and b represent the self-thinning index.

The DBH of each stand was estimated by incorporating not only stand-level factors such as stand age, SI, and stand density but also topographic and climatic factors to reflect various conditions of forests nationwide. For this, two different growth models were modified and integrated. First, Korea-specific regression models were used to calculate the basic (standard) pattern of DBH development by SI and stand age. These models, as described by Piao et al. [34], were developed to estimate the growth of DBH. These regression models were developed based on stand-level factors such as age, site index, and stand density using the fifth NFI. They showed relatively good performance with high correlation (0.710–0.789) for surveyed stand-level DBH at the national scale. To reflect the diameter growth response of trees to topographic and annual climatic variation, annual diameter growth by SI and age was calculated based on the derivative of these regression models with respect to age. Second, to reflect the annual diameter growth response of trees to topographic and climatic variation in DBH development, the estimated standard growth (eSG) model, as described by Kim et al. [35], was used as a growth modifier. To calibrate the differences between the impact of annual in-situ climate conditions on diameter growth of each stand and average climate conditions for each tree species in forests nationwide, we applied the growth modifier used in Kim et al. [35]. The following equation describes the growth modifier used here (Equation (4)):

(4)$db{h}_{i+1}=db{h}_i+2·\Delta r_{ij}·\left(\frac{eS{G}_{ij}}{meS{G}_p}\right)$

where i is each year from 1961 to 2014, j is the identification number of the stand, p is the total period (1961–2014), Δr is the mean estimated radial change from a differential regression model, eSG is the estimated standard growth, and meSG is the mean estimated standard growth of each species during the total period. The eSG_ij can be normalized using meSG_p. Because the normalized eSG_ij is integrated into the radial growth model, radial growth increments can be estimated from topographic and annual climatic conditions.

The KFS regression model was used to calculate the stem volumes of individual stands (Equation (5)). The coefficient of each tree species for Equation (4) was derived by the NIFoS [36]. To estimate all forest biomass carbon stocks (CS) from stem volume, a biomass expansion factor (BEF, defined as the ratio of all stand biomass to growing stock volume), which converts timber volume to mass and accounts for noncommercial components (such as branches and leaves), and root/shoot ratio (RS) were applied (Equation (6)). As with the previous research [37,38], the annual carbon sequestration (CSE) in forests biomass was calculated based on the difference between the CS of this year and the CS of the previous year (Equation (7)).

(5)$V_{ij}=a·db{h}_{ij}^b·h{m}_{ij}^c·N_{ij}·BW{D}^{-1}/1000$

(6)$C{S}_{ij}=V_{ij}·BEF·RS·0.5$

(7)$CS{E}_{ij}=C{S}_{ij}-C{S}_{ij-1}$

where V is stem volume (m³ ha⁻¹); CS is forest biomass carbon stocks (Mg ha⁻¹); CSE is annual carbon sequestration in forests biomass (Mg ha⁻¹ yr⁻¹); dbh is the mean stand diameter (cm); N is the number of trees per ha; BWD, BEF, and RS are the basic wood density, biomass expansion factors, and root/shoot ratio, respectively, proposed by the NIFoS [36] (Table 1); i is the year; j is the identification number of each stand; and a, b, and c are species-specific coefficients.

For this study, Korean forests were divided into seven major forest types: red pine (Pinus densiflora), Japanese larch (Larix kaempferi), Korean pine (Pinus koraiensis), cork oak (Quercus variabilis), Mongolian oak (Quercus mongolica), mixed forest A (Mixed-A: red pine and cork oak), and mixed forest B (Mixed-B: red pine and Mongolian oak) based on FCM classification. For the purposes of this study, other coniferous tree species were classified as red pine, which comprise the highest proportion of tree species in South Korea. Other broad-leaved tree species in low (below 700 m) and high (over 700 m) altitude were defined as cork oak and Mongolian oak based on the fifth NFI [30].

2.3. Validation

In this study, the stem volume of each stand in South Korean forests from 1960–2014 was estimated by applying the integrated model then validated using the fifth and the sixth NFI, statistical data and other research. We compared model results for stem volume to those of the statistical yearbooks of forestry in South Korea between 1960 and 2014 and estimated R² and the Root Mean Squared Error (RMSE) to indirectly validate the estimated biomass C stocks.

In addition, we attempted to directly spatially validate the model results using NFI plot data for three spatial-scale approaches as follows: (1) the plot level—comparison between the estimated stem volume of a stand and the mean stem volume of one NFI plot, which is composed of one center plot and three sub-plots; (2) the watershed level—comparison between the mean estimated stem volumes of forest areas within each watershed area and the mean stem volumes of NFI plots in each watershed area; (3) the administrative level—comparison between the mean estimated stem volumes of forest area and the mean stem volumes of NFI plots within each administrative boundary (Table 2, Figure 3). Generally, the watershed area and the administrative boundary are considered the spatial units for planning national forest management in South Korea. Therefore, this multi-spatial scale approach to validation is expected to be useful for local governments in developing forest activity strategies.

3. Results and Discussion

3.1. Stem Volume, C stocks, and Annual C Sequestration

Changes in the C stocks, including stem, branch, and root stocks, of South Korean forests in the years 1961–2014 are depicted in Figure 4. The results of the C stock simulation in this study reflected climatic effects on annual forest dynamic growth on a plot level. As shown in Figure 4a–c, the modeled spatial distribution of the C stocks in forests showed that C stocks tended to increase gradually from the 1960s to the 2000s. The highest C stocks were in the northeast and the center of South Korea (dark pixels). Most of these areas are the Taebaek Mountain Range and the Sobaek Mountain Range. The Sobaek Mountain Range runs from the Taebaek Mountain Range in the northeast to the southwest.

There was an increase in the total C stock of biomass in South Korean forests, from 6.65 Tg C in 1961 to 476.21 Tg C in 2014 (Figure 4d). The effects of national reforestation programs on forest C budget are shown more clearly when C stock is mapped by decade (Figure 4a–c). The C stock of South Korean forests increased from 13.15 Tg C in the 1960s to 135.38 Tg C in the 1980s due to the success of the national-scale forestation program [10]. Averaged over the period from the 1960s to the 1980s, the annual increase in C stock was 6.11 Tg C yr⁻¹. Based on the national statistical data, around one million ha was restored in just six years (1973–1978) [39]. Our simulation estimated that by the 2000s, the C stock had increased drastically to 375.58 Tg at a rate of 12.01 Tg C yr⁻¹. This reflects the high growth rates of the young forest stands and the national sustainable forest management plan. The national focus has shifted since the 1980s from planting to tending trees. This contributes to increased forest values and decreased deforestation activities such as household fuelwood use, illegal logging, and slash-and-burn fields [39]. Although the total C stock of Cool forests has been lower than that of the other forests since the 1970s, the mean C stock (Mg ha⁻¹) and the annual C sequestration (Mg ha⁻¹ yr⁻¹) of Cool forests are still estimated to be higher than those of other forests. Changes in the mean C stock per ha (Mg ha⁻¹) of Cool and other forests in 1960–2014 are depicted in Figure 4e. There were increases in the C stocks of biomass in both Cool and other forests from 4.03 and 0.43 Mg C ha⁻¹, respectively, in 1960 to 102.43 and 73.76 Mg C ha⁻¹ in 2014. The annual C balances from biomass in Cool and other forests were 0.97 and 0.26 Mg C ha⁻¹ yr⁻¹ before the onset of reforestation programs (1960–1973). The annual C balances of biomass in both forest types increased at rates of 2.07 and 1.67 Mg C ha⁻¹ yr⁻¹, respectively, after the onset of those programs (1974–2014). These results imply that the Cool forests have played important roles in C balances, biodiversity, terrestrial ecosystems, and supply ecosystem services of South Korea [14].

The simulation showed that the annual differences in mean C sequestration (Mg ha⁻¹ yr⁻¹) between the Cool forests and other forests started to gradually decrease after the 2000s. In addition, the estimated C sequestration of both Cool and other forests started to gradually decrease after the 2000s. According to both statistical data and model results, the stand age of more than half of South Korean forests reached 40 to 50 after afforestation. This result reflected the decreased growth rates of the aging forests of South Korea. Our model results also show the changes in annual C sequestration of each forest type (Figure 4g). The variations of annual C sequestration in both Cool forests and other forests were estimated to grow steadily after the 1980s. This is linked to recent trends in temperature and precipitation in South Korea. According to previous research, the warming trend of the mean annual maximum and minimum temperatures at most observational stations in South Korea have increased more rapidly after the 1980s than in earlier decades [40,41]. In addition, the phenomenon of a steep decline in C sequestration was observed more frequently during the 2000s than in other times in the model simulation. This result could be partially explained by the annual precipitation used as input data. South Korea has experienced much more frequent and extreme droughts since the late 1990s and early 2000s [42,43]. Accordingly, the annual growth rate and the tree mortality during that period were probably underestimated and overestimated, respectively, in our model. These drought events had a critical impact on South Korea, which led to regional water shortages and influenced the use of water, including for agricultural and household activities [44]. In addition, natural ecosystems were damaged by the drought, and vegetation indices on the national scale were low [33,35,45]. However, in model results, the growth peaks were observed after 2000. This result was because sharply rising temperatures and precipitation were reflected in the model. Both positive (e.g., increased water use efficiency, increases in forest vigor and growth from CO₂ fertilization, and longer growing seasons) and negative effects (e.g., increases in stress and mortality due to the combined impacts of climate change and climate-driven changes) of climate change on forests have been observed globally [33,46,47].

The model generally estimated the annual variation in mean C sequestration (Mg C ha⁻¹ yr⁻¹) of the Cool forests as larger than that of other forests during the simulation period. There are two possible explanations for this result. First, based on observed data, the high-elevation areas have been warming faster than other regions in South Korea during the five decades of concern (Figure 1). The differences in warming between the Cool forests and other forests during the simulation period influenced the model results. The other possible reason was a difference in the tree species composition between the Cool forests and other forests. Based on the FCM data used as input, the red pine (Pinus densiflora), the Mongolian oak (Quercus mongolica), and the mixed forests occupied approximately 10.3%, 43.3%, and 26.9%, respectively, of the Cool forest areas. On the other hand, these tree species occupied approximately 36.4%, 9.9%, and 31.3%, respectively, of the other forest areas. According to previous studies for not only South Korea but also worldwide, the impacts of climate change on the patterns of tree growth, mortality, and distribution vary by tree species and region [19,33,35,46,48].

3.2. Model Validation

The estimated stem volumes simulated by the model were compared to the Statistical Yearbook of Forestry [11] and to previous research results from South Korean forests (Table 3). The time sequence of estimated stem volumes showed a trend similar to that of observed stem volumes on a national scale (R² = 0.96, RMSE = 9.13; Figure 5a). According to the Statistical Yearbook of Forestry [11], the stem volume in South Korean forests increased from 9.6 to 142.2 m³ ha⁻¹ between 1961 and 2014. The simulation showed that it increased from 3.3 to 145.4 m³ ha⁻¹ during that period. As the results were similar, it could be inferred that the model successfully reflected the volume trends of Korean forests, which further suggests successful reconstruction of age distribution at the national scale, while uncertainties remain for individual stands.

The results of this research are consistent with other studies showing a large increase in biomass C stocks after the onset of reforestation programs. However, the mean biomass C density here was slightly higher than in other studies (Table 3). There are two possible reasons for these differences. First, as shown in Figure 5a, the model estimate of stem volume in recent years was higher than the observed stem volume. This caused higher estimates of recent mean biomass C density and stocks and annual C balance. The other possible reason was a difference in the methods of biomass C stock estimation. We estimated the biomass C stocks with species-specific growth functions and the latest BEFs developed by NIFoS [36]. In contrast, Choi and Chang [15] estimated biomass C stocks by multiplying stemwood volume with forest type-specific (coniferous, deciduous, and mixed) and constant BEFs developed from global data. Variable BEFs can overestimate biomass C of a young forest compared to constant BEFs [49]. As South Korean forests are relatively young, the estimates of biomass C stocks of South Korean forests with variable BEFs could be higher than the estimates with constant BEFs. Kim et al. [9] considered clear cut and thinning activities under the KFS law, but our research did not include a certain amount of timber harvest. Therefore, it stands to reason that our result is relatively higher than the previous study [9]. Lee et al. [22] used only plot-level data of the fifth NFI to estimate C budget for the entire forest area. In addition, the impacts of climate change on forest growth were not considered.

To assess model reliability, we compared the model results to observed values at various spatial scales (plot, watershed, and first-tier administrative levels). The estimated R² for model results in each spatial scale varied greatly (Figure 5b,c). The R² for plot, watershed, and administrative levels of fifth NFI were estimated at 0.25, 0.64, and 0.76, respectively (Figure 6). The model results of sixth NFI showed a similar pattern but were slightly lower than the R² for each spatial scale of fifth NFI. Accuracy was poor at plot level but increased when aggregations of cells were evaluated. This is because South Korean forests have highly heterogeneous characteristics among stands, even within a small area. The SD (standard deviation) of four plots’ stem volumes (m³ ha⁻¹) in one permanent plot of the NFI was used as a measure of forest structure heterogeneity in this study. One permanent plot from the Korean NFI consists of four subplots, such as one center plot and three other subplots, with the center subplot in the center of the plot and the three other subplots located 50 m away at azimuths 0°, 120°, and 240° from the center subplot. SD of the stem volume ranged from 0.0 to 165.0 and averaged 36.4 in 2860 plots (Figure 5b,c). The numbers of plots with SD < 40, 40 ≤ SD < 80, 80 ≤ SD < 120, and SD ≥ 120 were 1858 (65.0%), 832 (29.1%), 149 (5.2%), and 21 (0.7%), respectively. There was also a difference in the spatial scale between the estimations and the observed data. The model results represented the mean stem volume of forested areas in 1 km × 1 km grid cells. However, the observed stem volume of each NFI plot was calculated from all trees with at least 6 cm DBH in 0.08 ha with a 16 m radius. This implies that our model results were not useful at the plot level but could be useful at watershed and administrative levels.

Figure 7 shows the ranges of stem volumes in the Cool forests and other forests of South Korea for two time periods. According to the fifth NFI, the sixth NFI, and the model results, the total stem volume of Cool forests was higher than that of other forests. The mean stem volume of the Cool forests and other forests in the fifth NFI were 146.0 (SD: 70.5) and 123.1 m³ ha⁻¹ (SD: 74.4), respectively. After five years, these increased to 170.7 (SD: 77.4, 4.95 Mg ha⁻¹ yr⁻¹) and 147.2 m³ ha⁻¹ (SD: 84.0, 4.82 Mg ha⁻¹ yr⁻¹), respectively, according to sixth NFI data. Our model similarly showed that the estimated mean stem volume of the Cool forests and other forests changed from 149.7 (SD: 57.4) and 122.4 m³ ha⁻¹ (SD: 54.2), respectively, in the period 2006–2010 to 172.5 (SD: 59.6) and 146.7 m³ ha⁻¹ (SD: 57.3) in the period 2011–2015. However, the ranges of forest stem volume differed. As seen in Figure 7, the ranges of the observed stem volume values were wider than the estimated values. One major reason for this result is that the stem volume is a result of complex interactions between numerous factors such as environmental conditions, natural disturbances (disease and insect pests, wind damage, landslides), and observation errors. Another cause is that this study used nationwide forest inventory data that included forests in atypical locations such as urban areas and islands. Therefore, this range discrepancy does not represent an important issue for the reliability of our model. Because the observations and the estimations were similar, it can be inferred that our model appropriately reflected the forest volume and C stock change according to climatic variables on regional and national scales.

3.3. Uncertainty and Implications for Forest Management

In order to measure the function of the forests as a carbon sink, BEFs were used to estimate the carbon content of a whole tree using stem volume as input data. The biomass expansion factors developed by NIFoS [36] are dependent on tree species, which provide a relatively firm basis for quantifying forest carbon stock. The NIFoS report estimated the average uncertainties of basic wood density, biomass expansion factor, root/shoot ratio at 3.6%, 7.0% and 14.9%, respectively. Each coefficient was designed specifically for different tree species and is not influenced by tree size or age. However, according to many previous studies, these coefficients may also depend on site [50], tree size, and age [51]. However, we assumed that all C emission factors had constant values. Our model results reflected these uncertainties. In this study, we assessed only the C stock and the C sequestration in the tree layer; the total estimated forest C stock would be considerably higher if soil C and C in litter and ground vegetation were included. In addition, to enable a more precise and comprehensive assessment of South Korean forest C cycles, important influences on C balance, such as atmospheric CO₂ fertilization [52], leaching [53], N deposition [54], and other disturbances [46,55,56], need to be considered. In particular, natural disturbances are integral drivers of forest dynamics and contribute to the diversity and the adaptive capacity of ecosystems. Globally, climatic changes have been identified as a key driver behind increased disturbances in forests, as knowledge about climate change impacts on forests is continuously expanding [57]. Increasing disturbances could strongly impact forest ecosystem services and are suggested to contribute to the recently observed ecosystem functions in forests. Therefore, further research and future forest policies will require a stronger focus on disturbance risk and resilience. These factors must be considered in further research involving predictions of forest C dynamics under climate change and socio-economic scenarios.

Moreover, if substitution effects (such as forest regeneration) from using forest biomass as raw material instead of fossil fuels and carbon stored in wood products had been taken into account, this would also have affected the estimated carbon balance [58]. Climate mitigation effects are dependent on the system boundaries, e.g., the temporal and the spatial scales considered as well as the extents of the life cycles of forest products included in the analysis. For instance, some recent studies indicate that substitution effects and storage of C in harvested wood products may result in higher overall C sequestration on a long-term basis than storing C in the forest would [59]. These issues should be considered in future studies to widen the scope of understanding the effects of forests and forestry on the C balance.

Changes in species distribution, tree range shifts, and amplified tree mortality phenomena observed throughout the world are associated with climate change. Such changes including mortality of indigenous trees, expansion of subtropical forests, and the interactions of drought and insects leading to an increase in tree mortality rates are also seen in Korean forests. Other challenges include forest cultivation and management practices. In the last five decades, forest management practices in South Korea have focused on restoration and conservation activities [39]. Despite the substantial progress made thus far, there is still room for improvement in terms of the disproportionate tree age structure, the decreased annual carbon sequestration, the low self-sufficiency rate for forest resources, and the regulation of the water table. Enhanced forest management strategies need to be set out in light of these challenges lying ahead.

4. Conclusions

The national-scale afforestation program in South Korea is a successful program that has greatly improved the functioning of forest ecosystems. The key findings of this study are summarized as follows. Changes in Cool forests and other forests of South Korea, as defined in the 2018 IBFRA, were estimated for 1961–2014 through a model developed using Korean NFIs, FCM, and climate data. Our results indicated that the developed model accurately described the growing stem volume and the forest biomass C stocks at the watershed and the administrative levels, which are basic units for building public forest management plan. Therefore, the developed model could be useful for building the adaptive and the risk-resilient sustainable forest management in South Korea to cope with climate change. Simulation results showed that the mean C stocks of the Cool forests and other forests changed from 4.03 and 0.43 Mg C ha⁻¹ to 102.43 and 73.76 Mg C ha⁻¹ during the period 1961–2014. These results imply that, although the Cool forests of South Korea were only about 12.3% (772,788 ha) of the total forest area, the Cool forests have played important roles in C balances and the forest ecosystem of South Korea. Another major finding of this study is that Cool forests in South Korea have been identified as more highly sensitive to climate change than other forests. In addition, annual carbon sequestrations are projected to decrease slowly after the 2000s because the overall growth rate of mature forests decreases as stands age. Therefore, the next level of forest management in South Korea should address the following issues: (1) adaptation to climate change through forest management strategies suitable to each region or forest type; (2) maintenance and promotion of annual C sequestration; and (3) focus on development of policies and systems for using forest biomass such as bioenergy production, using raw material instead of fossil fuels for their substitution effects, and storage of C in harvested wood products. The developed model and the results of this study can support the improvement of spatiotemporal forest management strategies to allow for sustainable development of society and to cope with climate change in South Korea.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4907/10/7/579/s1, Table S1: Coefficients for the dominant tree height equation and stand-level mortality model by the Korea Forest Service [31] and Kim et al. [33].

Author Contributions

Conceptualization, M.K., F.K., Y.S., S.W.J. and W.K.L.; Data curation, B.Y.H., C.H.L., C.S. and M.H.; Formal analysis, M.K.; Funding acquisition, W.K.L; Methodology, M.K., A.S. and W.K.L.; Supervision, F.K., A.S, D.S. and W.K.L; Validation, B.Y.H., C.H.L. and C.S.; Visualization, B.Y.H., C.H.L. and C.S.; Writing – original draft, M.K.; Writing – review & editing, F.K., Y.S., S.W.J, A.S, D.S. and W.K.L.

Funding

This research was supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education (NRF-2018R1A6A3A11051238), and was carried out with the support of ´R&D Program for Forest Science Technology (Project No. 2017044B00-1919-BB01) provided by Korea Forest Service (Korea Forestry Promotion Institute). This study was also supported by Korea University Grants, and the Ministry of Environment (MOE) of the Republic of Korea (Project Number: 2014001310008).

Acknowledgments

The authors would like to thank the collaborative research between the Korea University and IIASA-ESM-CLR (International Institute for Applied Systems Analysis-Ecosystem Services and Management-Center for Landscape Resilience and Management). The authors are grateful to the International Boreal Forest Research Association (IBFRA) for encouraging innovation and knowledge-sharing of Cool forests.

Conflicts of Interest

The authors declare no conflict of interest.

Figures and Tables

Figure 1. (a) Study area and its (b) elevation map and (c) forest area map. (d) Annual mean temperature (left) and annual precipitation (right) for South Korean forests. Cool forests were defined as forests at elevations above 700 m in the temperate zone based on the 2018 International Boreal Forest Research Association (IBFRA) Conference. T: annual mean temperature, P: annual precipitation.

Figure 2. (a) Locations of National Forest Inventory (NFI) plots from which data were obtained. (b,c) Estimated site index and stand age distribution in 2010 of South Korean forests using the fifth NFI data and 1:5000 forest cover map (FCM). The distribution of (d) age classes in the FCM and (e) estimated stand age using the FCM and fifth NFI.

Figure 3. (a) Administrative boundaries of South Korea and an example of administrative-level spatial data for model validation. (b) Watershed boundaries of South Korea and an example of watershed-level and plot-level spatial data for model validation.

Figure 4. (a–c) Estimated forest carbon stock distribution in South Korea during the 1960s, the 1980s, and the 2000s. (d,e) Time series of total (Tg C) and mean C stocks (Mg C ha−1) for three categories in South Korean forests during the simulation period. (f,g) Time series of total (Tg C) and mean annual C sequestration (Mg C ha−1) for three categories in South Korean forests during the simulation period.

Figure 5. (a) Time series of the estimated and the observed (statistical data) mean stem volumes in South Korean forests during the simulation period. Distribution of standard deviations in observed stem volume of sub-plots within permanent plots from the (b) fifth and the (c) sixth Korean National Forest Inventory (NFI) data.

Figure 6. Observed and estimated stem volumes (m3 ha−1) at plot, watershed, and administrative levels.

Figure 7. Distribution plots of estimated and observed forest stem volume (m3 ha−1) for Cool forests and other forests. The fifth and the sixth National Forest Inventory data corresponded to the estimated mean values in 2005–2010 and 2011–2015. For each box plot, the top circle is the 95th percentile, the lower circle is the fifth percentile, the top bar is the 90th percentile, the lower bar is the 10th percentile, the top of the box is the upper or the third quartile, the bottom of the box is the lower or the first quartile, and the middle bar is the median value.

Table 1

Coefficients for stem volume equations and carbon emission factors by tree species, developed by the National Institute for Forest Science [36]. The coefficients for each tree species were applied to both Cool and other forests.

Tree Species	Stem Biomass Equation			Carbon Equation
	a	b	c	Basic Wood Density	Biomass Expansion Factor	Root/Shoot Ratio
Pinus densiflora	0.034	1.734	1.025	0.472	1.413	0.254
Pinus koraiensis	0.046	1.732	0.896	0.408	1.812	0.283
Larix kaempferi	0.005	2.458	0.904	0.453	1.335	0.291
Quercus variabilis	0.053	1.810	0.881	0.721	1.338	0.324
Quercus mongolica	0.098	1.406	1.135	0.663	1.603	0.388

Table 2

Descriptive statistics of validation data by spatial level.

Factors	Spatial Level (Std. dev.)
	Plot	Watershed Boundary	Administrative Boundary
No. of sample	3892 (n/a)	102 (n/a)	17 (n/a)
Mean area (ha)	0.24 (n/a )	93,514 (61,244)	590,374 (620,349)
Mean stem volume (m3 ha−1) offifth NFI and sixth NFI	125.8 (74.8),149.9 (84.7)	118.3 (25.8),142.3 (28.3)	124.9 (15.7),149.1 (18.7)

Table 3

Comparison of living biomass carbon (C) densities of South Korean forests estimated in this study with those found in previous studies.

Category	Year or Period	Model	Estimate	Reference
Mean C density(Mg C ha−1)	2001	Statisticalmodel	34.4	Choi and Chang [15]
		KFSC	39.7	Lee et al. [22]
		CBM-CFS3	34.1	Kim et al. [9]
		This study	42.1	This study
	2010	CBM-CFS3	68.0	Kim et al. [9]
		This study	67.5	This study