Next Article in Journal
Drivers of Wildfire Occurrence Patterns in the Inland Riverine Environment of New South Wales, Australia
Next Article in Special Issue
The Use of Satellite Information (MODIS/Aqua) for Phenological and Classification Analysis of Plant Communities
Previous Article in Journal
The Importance of Moisture for Brown Rot Degradation of Modified Wood: A Critical Discussion
Previous Article in Special Issue
A Holistic View of Soils in Delivering Ecosystem Services in Forests: A Case Study in South Korea
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Forests
Volume 10
Issue 6
10.3390/f10060523
forests-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Assessing Forest Ecosystems across the Vertical Edge of the Mid-Latitude Ecotone Using the BioGeoChemistry Management Model (BGC-MAN)

by
Cholho Song
1,2,
Stephan A. Pietsch
2,
Moonil Kim
2,3,
Sungeun Cha
1,2,
Eunbeen Park
1,
Anatoly Shvidenko
2,
Dmitry Schepaschenko
2,
Florian Kraxner
2 and
Woo-Kyun Lee
1,*
1
Department of Environmental Science and Ecological Engineering, Korea University, Seoul 02841, Korea
2
Ecosystems Services and Management, International Institute for Applied Systems Analysis (IIASA), Schlossplatz 1, A-2361 Laxenburg, Austria
3
Environmental GIS/RS Center, Korea University, Seoul 02841, Korea
*
Author to whom correspondence should be addressed.
Forests 2019, 10(6), 523; https://doi.org/10.3390/f10060523
Submission received: 31 May 2019 / Revised: 18 June 2019 / Accepted: 21 June 2019 / Published: 23 June 2019
(This article belongs to the Special Issue Cool Forests at risk? The Critical Role of Boreal and Mountain Ecosystems for People, Bioeconomy, and Climate)
Browse Figures
Versions Notes

Abstract

:
The mid-latitude ecotone (MLE)—a transition zone between boreal and temperate forests, which includes the regions of Northeast Asia around 30°–60° N latitudes—delivers different ecosystem functions depending on different management activities. In this study, we assessed forest volume and net primary productivity changes in the MLE of Northeast Asia under different ecological characteristics, as well as various current management activities, using the BioGeoChemistry Management Model (BGC-MAN). We selected five pilot sites for pine (Scots pine and Korean red pine; Pinus sylvestris and P. densiflora), oak (Quercus spp.), and larch forests (Dahurian larch and Siberian larch; Larix gmelinii and L. sibirica), respectively, which covered the transition zone across the MLE from Lake Baikal, Russia to Kyushu, Japan, including Mongolia, Northeast China, and the Korean Peninsula. With site-specific information, soil characteristics, and management descriptions by forest species, we established their management characteristics as natural preserved forests, degraded forests, sandy and cold forest stands, and forests exposed to fires. We simulated forest volume (m3) and net primary productivity (Mg C ha−1) during 1960–2005 and compared the results with published literature. They were in the range of those specified in previous studies, with some site-levels under or over estimation, but unbiased estimates in their mean values for pine, oak, and larch forests. Annual rates of change in volume and net primary productivity differed by latitude, site conditions, and climatic characteristics. For larch forests, we identified a high mountain ecotype which warrants a separate model parameterization. We detected changes in forest ecosystems, explaining ecological transition in the Northeast Asian MLE. Under the transition, we need to resolve expected problems through appropriate forest management and social efforts.

1. Introduction

Forest ecosystems have experienced changes owing to climate change and anthropogenic pressure [1,2], and these changes are particularly rapid in areas with a sensitive ecosystem structure and high population density. In addition, forest degradation, deforestation, and desertification are the main issues that threaten sustainable terrestrial ecosystems, therefore maintaining ecological productivity is important [3]. The northern mid-latitude region (MLR), which is located between 30° and 60° N latitudes with various land covers and a lot of human activities, including the mid-latitude ecotone (MLE), is a transition zone between forests and drylands [4,5,6]. In Eurasia, dense vegetation is distributed along the MLE in regions such as Northeast Asia, Siberia, and Europe, however, drylands also exist in Mongolia, Central Asia, and the Mediterranean region [7,8].
Forests in the MLE consist of boreal and temperate forests, and form a belt connecting the Far East and Europe. This forest belt has important roles in maintaining ecological sustainability and preventing the expansion of drylands. The Northeast Asian MLE, including Korea, China, Japan, and Mongolia, has faced environmental changes during the past decades and will be affected by future climate change [9]. Through the forests of Mongolia and China, the boreal forests of Russia near Lake Baikal and Siberia are connected to the temperate forests in China and the Korean Peninsula [10]. Along this connection, each nation has different forest management issues and environmental conservation regimes.
Temperate forests in South Korea and Japan were restored during the 1960s to 1970s, and now both countries have successful restoration outcomes, with the aim of sustainable forest management [11]. As a case in point, the forest growing stock volume (GSV) in South Korean forests was 10.5 m3/ha in 1952, and it reached 100.4 m3/ha in 2007 [12]. However, forests in North Korea have suffered from deforestation and forest degradation due to logging for food and fuel production [13,14]. There are some differences in reports on the extent of deforestation in North Korea; the deforested areas there are estimated to have ranged between 1.2 million ha and 2.5 million ha during the late 1980s to 2000s [15]. These differences are also represented in the forests’ net primary productivity (NPP) as a negative carbon budget in North Korea [16]. Temperate forests were changed to boreal forests in China, Northern Mongolia, and Russia near Lake Baikal. North China and Inner-Mongolia are widely known as regions facing ongoing deforestation and land degradation [17], but the Chinese and Mongolian governments strive to recover their forest ecosystems to secure ecological benefits [18,19,20], whereas, the forests near Lake Baikal and around the MLE have suffered natural or anthropogenic forest fire [21]. There is great uncertainty in measuring the effect of fire in boreal forests; one estimate revealed that 2.9% of Russian boreal forest loss was caused by forest fires [22].
Therefore, forest management regimes in the Northeast Asian MLE can be summarized as: (1) well-restored and managed areas in South Korea and Japan; (2) land degradation in North Korea; (3) semi-arid and new plantation in China and Mongolia; and (4) cold temperature with forest fire in Russia. The complexity of the management regimes of different countries leads to different forest statuses. NPP and GSV can represent the effects of these differences, and many modeling approaches were applied in the MLE of Northeast Asia using Flux data, Moderate Resolution Imaging Spectroradiometer-Leaf Area Index (MODIS-LAI), Boreal Ecosystem Productivity Simulator (BEPS), Vegetation Integrative Simulator for Trace gases (VISIT), MAPSS-CENTURY1 (MC1), and Biome-BioGeoChemistry (Biome-BGC) [15,16,23,24,25,26,27]. Forest fluxes were calculated from data obtained from field surveys, but such data were procured at different times and locations [10,24]. Therefore, the use of forest fluxes to understand different management regimes has limitations. Applying models can be the alternative to overcome the spatio–temporal limits of the fluxes. Both satellite-based and process-based ecosystem models were applied to simulate productivities in the Northeast Asian MLE [15,16,28,29,30]. However, these modeling approaches had limitations in reflecting management details because of different parametrizations and representations of ecological processes and spatio–temporal data limitations.
Our synthetic study elucidates the process of applying the BioGeoChemistry Management Model (BGC-MAN) in the MLE of Northeast Asia. The BGC-MAN, a modified mechanistic biogeochemical model, extends simulations into forest management options in certain forest stands. The model was tested in European forests but has not yet been tested in the MLE of Northeast Asia [31,32,33,34]. We assessed past changes in NPP and GSV in the MLE under different environmental characteristics, as well as under various current management activities, using the BGC-MAN. We selected representative forest stands in the MLE of Northeast Asia from literature and assessed changes to the ecosystem and model uncertainties from the simulated results.

2. Materials and Methods

2.1. Applying BGC-MAN

BGC-MAN is a process-based ecosystem model, and the model also has the modules to apply forest management activities, so different forest ecosystem and management activities were analysed by BGC-MAN [31,32,33,34]. This process-based model requires more than 50 inputs, with daily weather data, bio-physical information, species characteristics of the forest stand, and site information. The model calculates different combinations by the climatic variables (e.g., daily weather data, solar radiation, vapor pressure, precipitation), bio-physical status (e.g., aspects, elevation, and soil and nitrogen characteristics), and species characteristics (e.g., leaf area index (LAI), stomatal conductance, carbon nitrogen ratio). From photosynthesis, the self-initialized model generates gross primary productivity (GPP), NPP, and the carbon allocation of the vegetation’s carbon pools, including respirations.
The working steps for applying BGC-MAN in the MLE of Northeast Asia start from selecting study sites. Under the limited carbon flux sites and temporal availability, the study sites were selected based on the literature review. From the literature, we can figure out dominant tree species as well as the bio-physical status and management history of forest stands. To run the model, missing data and parameters were collected and adjusted by geo-spatial analysis and previous BGC-MAN studies in Europe. In addition, we adopted daily weather data from 1951 to 2005 to find current changes of the forest ecosystem. After that, we simulated forest volume (m3) and NPP (Mg C ha−1). We checked the results with different literature and general trends among the study regions.

2.2. Study Sites

To cover the MLE in Northeast Asia, we selected five pilot sites from each of the southern parts of the selected sites, including oak (Quercus spp.), pine (Pinus sylvestris and Pinus Densiflora), and larch (Larix Mongolica and Larix Sibirica) stands [35,36,37,38] (Figure 1). The study sites were distributed along latitudinal (vertical) transects. They were also located near the borders of forests with other land covers and different forest densities [7,39]. In addition, various human effects due to different management of forest activities were found from literature [6,40]. Therefore, this region has ecological and socio-economic transitions that vary according to each tree species and nation.
Oak forests are distributed in the Korean Peninsula, Japan, and Far East Russia. The Kasuya stand, located in Kyushu, Japan, is one of 34 long-term monitoring forest sites of JaLTER (Japan Long Term Ecological Research Network), so forest disturbance has been very limited here in comparison with the other forests [35,41]. Based on the literature review, its forest age is around 140 years. The Gwangneung stand in South Korea has been well-preserved as it is a national arboretum, therefore it has a very old forest stand age, around 80–200 years [16,36,42]. In addition, many field studies by the National Institute of Forest Science (NIFoS), as well as individual previous studies, have monitored forest changes. The Seolmacheon stand in South Korea was also recently added to the AsiaFlux network [43,44]. This forest stand is young, aged around 20–40 years, and some forest hydrology experiments have been conducted here. Detailed forest characteristics of the Sohung stand in North Korea, one of the degraded and partially restored areas, were assumed based on a previous case study for Reducing Emission from Deforestation and land Degradation (REDD+) [45]. Bikinsky is located in Far East Russia, near the border of China, and forest fires are common in the surrounding forests [46].
Pine forests are widespread around the Korean Peninsula, Northeast China, and near Lake Baikal. Gwangneung also has old protected pine forests for royal timber production, and presently has well-preserved natural forests [38,47,48]. Its forest age is more than 80 years. Wonsan, which is located on the east coast of North Korea, is one of the focus areas for restoring deforested North Korean forests [11,49]. Honghuaerji forest is located in a sandy area with a relatively young forest stand age. This site was established by Chinese government projects to restore the sandy lands and reduce soil desertification [50,51]. One of the Russian forest stands, named Oblast, is located on the western side of Lake Baikal near Irkutsk, and the other, named Pribaykalsky, is located on its eastern side. Both forest stands experience forest fires and they are some of the coldest among the pine tree sites [46].
The larch trees in Laoshan, China were planted in 1969 as an experimental forest, and various studies have been conducted in this site [26,37]. In Mongolia, the Tsenkhermandal and Khentei sites were selected for our simulation. Tsenkhermandal is located in the high mountains with a 110-year-old forest, and the Khentei site represents typical Mongolian larch forests near Ulaanbaatar. There have been different estimates of its forest age from previous studies, but it is estimated to be as old as 64–146, according to tree ring surveys [52,53,54]. Yeravninsky is a 100-year-old forest and is located on the eastern side of Lake Baikal across the mountains there. Buryatia is the northernmost site and comprises a 100-year-old high mountain forest. Both forest stands experience forest fires [46]. From these previous studies, a site information table, especially including information on species, forest age, and location with management history, was constructed (Table 1).
To add more details, site-specific management information and mortality rates were adopted from the previous studies. For the pine forests, a site-specific mortality rate of 1.4% was applied for the Gwangneung forest [38], but the average mortality rate was 3% in the other Korean young forests [55]. The average mortality of oak forests in Korea was 1.7–2.0%, however, it was as high as 4.2% in some oak stands [55]. We assumed that the mortality of old forest stands was lower than average but allowed a higher mortality rate for young forests. In the case of Russian forests, we applied an average mortality rate of 0.6% through calculations from data on growth and productivity [46]. Similarly, the mortality rate of larch forests from the Korean Peninsula was applied to the young forests in China. This value was 5.0–5.6% in L. kaempferi. For the other forests, a mortality rate of 0.2–0.3% was applied.
We set management characteristics as: (1) well-restored and managed areas in South Korea and Japan; (2) land degradation in North Korea; (3) semi-arid and new plantation in China and Mongolia; and (4) cold temperature with forest fires in Russia. Following these differences, the model adopted options of planting, clear-cutting, natural preservation, and logging. From the forest age, we initialized the planting and clear-cutting options for modeling to fit the respective management history. For example, many sites in South Korea and Japan were planted and naturally preserved if these were their historical management strategies, but we assumed clear-cutting in North Korea, following their history of deforestation.

2.3. Data Processing

To generate sufficient bio-physical data, which is required for running the BGC-MAN (Table 2), we supplemented the data with global scale GIS data, such as the Shuttle Radar Topography Mission (STRM) and Harmonized World Soil Database (HWSD). For the nitrogen deposition values and pre-industrialization and current values, we adopted global modeling results [56,57]. Parameters of nitrogen fixation were applied from several previous studies, which suggested different nitrogen fixation values of mature oak forests (12.1–14.0 N kg ha−1 year−1), pine (0.3–3.2 N kg ha−1 year−1), and mixed larch hardwood (0.6–4.9 N kg ha−1 year−1) [58,59,60]. For the larch forest, we followed previous studies, with consideration of the nitrogen fixation value from a mixture of western white pine and western larch [60,61]. We adopted parameters for species characteristics from Europe (the west side of the MLE), but we also adjusted some of the carbon and nitrogen allocation, as well as night time freezing temperature values [10,33,62].
For the daily weather data, we used global climate data at a grid resolution of 0.5° from the Environmental Policy Integrated Model (EPIC) and modified them through reflecting temperature lapse and horizontal angle of solar radiations using the Mountain Microclimate Simulation Model (MTCLIM) [63,64,65,66]. Because the study sites were vertically distributed, climatic characteristics of this region can be described by comparing the northernmost site (Buryatia, Russia) and the southernmost site (Kasuya, Japan) (Figure 2). In Russia, the temperature dropped to −34.5 °C in the winter, and climbed to 21.04 °C in the summer, whereas the temperature varied from −11.31 °C to 27.38 °C in Japan. The climatic variability of the MLE also could be explained by the warmth index (WI) which was used to present the distribution and suitability of vegetation [67,68]. The distribution of the WI in Kasuya, Gwangneung, Seolmacheon, Sohung, Bikinsky, Wonsan, Oblast, Pribaykalsky, Laoshan, Tsenkhermandal, Khentei, Yeravninsky, and Buryatia were 97.04, 123.36, 102.09, 142.98, 70.90, 84.12, 68.06, 62.89, 57.90, 85.57, 46.84, 63.56, 55.43, and 23.26, respectively. In addition, the accumulated precipitation in the sites in Russia and Japan was 30.41 cm and 157.29 cm, respectively. These differences contributed to different ecological transitions along this the region [5,9].

3. Results

3.1. GSV and NPP

From the BGC-MAN simulations, the estimated values of the GSVs of oak forests in Kasuya (Japan), Gwangneung (South Korea), Seolmacheon (South Korea), Sohung (North Korea), and Bikinsky (Russia) were 410, 511, 154, 65, and 153 m3 ha−1, respectively, in 2005. The GSVs of pine forests in Gwangneung (South Korea), Wonsan (North Korea), Honghuaerji (China), Oblast (Russia), and Pribaykalsky (Russia) were 392, 145, 147, 282, and 124 m3 ha−1, respectively, during the same period. The GSVs of larch forests in Laoshan (China), Tsenkhermandal (Mongolia), Khentei (Mongolia), Yeravninsky (Russia), and Buryatia (Russia) were 127, 156, 150, 227, and 116 m3 ha−1, respectively (Figure 3).
Along the MLE in Northeast Asia, the GSV ranges were wide in the order of oak, pine, and larch forests. Protected forests showed much higher GSV, but GSVs in this region converged around 154 m3 ha−1. Following the same order of sites among all forest types as for the GSV results, NPP values in oak forests were also estimated to be 5.05, 5.39, 5.22, 4.18, and 3.70 Mg C ha−1, respectively, in 2005. The NPP values of pine forest sites were 6.39, 6.75, 2.68, 3.49, and 2.50 Mg C ha−1, and the NPPs of larch forests were 5.24, 1.12, 1.35, 1.27, and 0.49 Mg C ha−1, respectively (Figure 4).

3.2. Comparing with the Literature

The estimated GSV and NPP values from the self-initialized model were compared with the previous literature. Many field studies, as well as modeling results for GSV and NPP, were available for some of our study sites, so we could compare our results with them directly. We obtained GSV and NPP values from previous reports if the data were available, however, we also assumed the values from diameter at breast height (DBH), gross primary productivity (GPP), and stand height in order to overcome the data limitation.
For the oak forests, the forest volume of Kasuya (Japan) was estimated at around 109–516 m3 ha−1, when considering some old forest stands. To estimate the value, we obtained the DBH and stand density from JaLTER and Inoue et al. (2008) [35,69], and applied the biomass expansion equation from the National Institute of Forest Science (NIFoS), Republic of Korea [70]. The tree volume of oak forests in Gwangneung (South Korea) was around 550 m3 ha−1 [36], and estimated annual NPP was 4.30–6.70 Mg C ha−1 [42,71]. In Seolmacheon (South Korea), we could estimate the tree volume from the DBH on National Forest Inventory (NFI) data, which was around 78–159 m3 ha−1 [43]. In addition, using the hydrology model named RHESSys (Regional Hydro-Ecological Simulation System), we calculated GPP as 11.93–13.96 Mg C ha−1 [72]. Sohung (North Korea) is one of the reforestation areas in North Korea, and a previous study assumed some forest characteristics, however, it estimated a forest volume of around 56–66 m3 ha−1 during 1989–2013 [45]. Bikinsky (Russia) is located in Far East Russia near the Chinese border. The forest volume here was estimated at around 154 m3 ha−1, and annual NPP at 5.05 Mg C ha−1.
The tree volume of the Gwangneung (South Korea) pine forest was estimated from a previous DBH survey at around 378–645 m3 ha−1, and a National forest inventory (NFI) plot located near the site indicated 270 m3 ha−1 [38]. The predicted annual NPP was 5.37–7.20 Mg C ha−1 [48]. Wonsan (North Korea), which was deforested and subsequently reforested, has an annual NPP range of 4.30–6.03 Mg C ha−1 based on the VISIT model simulation [16]. The forest volume of Honghuaerji (China) was estimated as 117–173 m3 ha−1 [50]. In the case of Oblast (Russia), the forest volume was estimated at around 282 m3 ha−1, and annual NPP was estimated at around 3.06 Mg C ha−1. The GSV of forests in Pribaykalsky (Russia) was reported to be around 102 m3 ha−1, and this forest’s NPP was predicted to be 2.72 Mg C ha−1.
For the larch trees of Laoshan (China), which are part of an experimental forest, GSV was estimated at around 117.86 m3 ha−1 in 2001, therefore a GSV value exceeding 133.36 m3 ha−1 would be expected in 2005 [73]. Other model simulations estimated the 300-year cumulated NPP to be 215.4–236.6 kg m−2, which can be converted to 7.18–7.89 Mg C ha−1 [26], but these estimations have big variations, as these were young forest stands. The GPP of the forest in Khentei (Mongolia) was reported to be around 5.62–6.10 Mg C ha−1 [24], and the GSV was estimated to be 64–146 m3 ha−1 via tree ring survey [52,53]. The forest volume of Yeravninsky was reported to be 227 m3 ha−1, and annual GPP as 3.38 Mg C ha−1. In Buryatia (Russia), the forest volume was estimated to be 117 m3 ha−1, and annual GPP as 2.81 Mg C ha−1 (Table 3).
Furthermore, the simulated NPPs were compared with the Moderate Resolution Imaging Spectroradiometer (MODIS)’s NPP (MOD17A3H) product from the National Aeronautics and Space Administration (NASA), the United States [74]. Since the MOIDIS NPP product data are available from 2000 and our study has the BGC-MAN results until 2005, the six years of overlapping NPP results were compared to find the relationship (Figure 5).
The NPP from the BGC-MAN simulation with climatic variations and the MODIS NPP from different characteristics of spectral reflectance may vary from time to time, but overall, they showed a positive correlation. The BGC-MAN results were underestimated compared to the NPPs from the MODIS NPP product, and the variations were shown in some larch forests because there was under estimation of the NPP from BGC-MAN.

4. Discussion

We identified changes to the forest ecosystems in the Northeast Asian MLE, and these results can explain some ecological transitions there. Under these transitions, we need to resolve expected problems through appropriate forest management and social efforts. The ecological transitions in the MLE of Northeast Asia are apparent from the annual rates of change of the forest volume and NPP. These values were estimated for a total of 15 sites under different environmental conditions and different forest management activities, and the simulation results were quite different according to latitude, site conditions, and climatic variations. At the species level, the overall BGC-MAN results from pine and oak forest sites were quite acceptable, even though we used parameterizations from European pine and oak forests to simulate P. sylvestris, Q. robur, and Q. petraea agg. characteristics here [31,33,34,62,75], however, there were some instances of under or over estimation. These BGC-MAN parameters were effective for P. densiflora and the other oak tree species in the MLE of the Northeast Asia. The GSVs of larch forests from Laoshan, Tsenkhermandal, Khentei, and Buryatia were quite similar to previous reports, but the estimated NPP was very low. This implies that, although our simulations were already based on modified biomass allocation parameters, the European parameters of L. decidua need to be adjusted further for the simulation of L. Mongolica and L. Sibirica [10,76]. This variation might be caused by some ecotypes and ecoregions on the high mountains [9,52,77]. With the heavy wind and cold temperature, larch forests in Mongolia and Russia were differentiated in stomatal conductance and the tree root system in comparison with European larch forests and larch forests in China and the Northern Korean Peninsula.
Another reason for these variations in larch forests could be the over estimation of forest fires in the self-initialized model. We assumed values for depicting forest fires, GSV, and NPP based on an inventory [46], and there might be differences in spatial resolution between the published reports and modeling stands. This means that the statistics aggregate the effect of forest disturbance and management, but the model was directly affected by the initialized parameters. When we use a process-based ecosystem model, the model usually simplifies environmental changes [10,78]. Furthermore, the model was self-initialized, so the generalization of established parameters would be limited in specific ecosystems [62]. Therefore, the model requires more adjustments in order to be suitable for application in specific ecosystems. However, the historical information of the sites was limited, and this gave rise to difficulties in the first initialization of the model. Therefore, more detailed management and human intervention data in a site-specific manner will be required in future research to apply the model to specific sites with greater accuracy.
For the study sites subjected to simulations in this study, the reported NPP values along the Northeast Asian MLE were 7.21 Mg C ha−1 for Mt. Jumbong, South Korea [79]; 5.16–8.04 Mg C ha−1 around the Korean Peninsula [16]; 2.04–6.20 Mg C ha−1 along the border of the Korean Peninsula and China [80]; 3.23 Mg C ha−1 for the boreal forests (L. gmelinii) of China [81]; 2.45 Mg C ha−1 for the larch forests of Eastern Siberia [82]; and 6.34 Mg C ha−1 among the East Asian regions, including our study area [23]. These estimated values also showed similar trends to our simulation results, which indicated a lower NPP at higher latitudes, although there were some limitations for larch forests. In particular, because the forest changed from temperate to boreal forest, these climatic influences are important for estimating ecosystem productivity [83,84]. In addition, there were many reports that MODIS NPPs were overestimated in Mid-Latitude region, so variations in NPP can be understood as general states [16]. In some forest stands without deforestation, such as in North Korea, the actual forest ecosystem’s productivity will be more positive than our simulations, for which we assumed dramatic deforestation and forest degradation due to logging for energy and food production [13,15]. Therefore, we found that the BGC-MAN was effective in explaining not only ecological and social transitions of forests in the MLR of Northeast Asia but could also explain future estimations under various scenarios of afforestation, if the parameters and other site-specific information were available.

5. Conclusions

Through our synthetic study, we aimed to elucidate the process of applying the BGC-MAN in the MLE of Northeast Asia, and we applied the model to oak, pine, and larch forests. There were some limitations, but the parameterizations of oak and pine species from European species were acceptable in similar ecosystems in the Northeast Asian MLE when we compared the simulated results of forest NPP and GSV with literature values. However, the parameterizations for the larch species, which were distributed over the entire study region, need to be modified to be as specific as possible for the different forest sites. Overall, our simulations were representative of the differences in NPP and GSV by latitude, parameter, and climatic variations, although there were limitations due to the self-initialized model and insufficient data in this region. Under this complexity, the results explained the changes in ecosystem along the latitudinal (vertical) ecological transects as well as the effects of differences in forest management under the model uncertainties. Therefore, we can conclude that the BGC-MAN is suitable for application in this region, especially for oak and pine forests, and it also can represent the different management regimes among the different regions studied. Therefore, our results can contribute to future research attempting to apply the BGC-MAN and process-based ecological models in the MLE of Northeast Asia.

Author Contributions

C.S., S.A.P., A.S. developed overall concept and methodology. C.S., S.A.P conducted analysis, M.K., S.C., E.P processed data. D.S., F.K., W.K.-L. supervised overall research. F.K., W.-K.L. supported fund for the research. All authors gave comments to improve and read the final manuscript.

Funding

This study was funded by the Ministry of Environment (MOE) of the Republic of Korea (Project Number: 2014001310008) and was carried out with the support of ‘R&D Program for Forest Science Technology (Project No. 2017044B00-1919-BB01 and Project No. 2018133B10-1820-BB01)’ provided by Korea Forest Service (Korea Forestry Promotion Institute).

Acknowledgments

This study was developed in the Young Scientists Summer Program (YSSP) at the International Institute for Systems Analysis (IIASA), Laxenburg, Austria with support from the National Research Foundation (NRF), Republic of Korea.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Soja, A.J.; Tchebakova, N.M.; French, N.H.F.; Flannigan, M.D.; Shugart, H.H.; Stocks, B.J.; Sukhinin, A.I.; Parfenova, E.I.; Chapin, F.S.; Stackhouse, P.W. Climate-induced boreal forest change: Predictions versus current observations. Glob. Planet. Chang. 2006, 56, 274–296. [Google Scholar] [CrossRef]
  2. Cowie, A.L.; Orr, B.J.; Castillo Sanchez, V.M.; Chasek, P.; Crossman, N.D.; Erlewein, A.; Louwagie, G.; Maron, M.; Metternicht, G.I.; Minelli, S.; et al. Land in balance: The scientific conceptual framework for Land Degradation Neutrality. Environ. Sci. Policy 2018, 79, 25–35. [Google Scholar] [CrossRef]
  3. Bai, Z.G.; Dent, D.L.; Olsson, L.; Schaepman, M.E. Proxy global assessment of land degradation. Soil Use Manag. 2008, 24, 223–234. [Google Scholar] [CrossRef]
  4. Neilson, R.P. Transient ecotone response to climatic change: Some conceptual and modelling approaches. Ecol. Appl. 1993, 3, 385–395. [Google Scholar] [CrossRef] [PubMed]
  5. Moon, J.; Lee, W.K.; Song, C.; Lee, S.G.; Heo, S.B.; Shvidenko, A.; Kraxner, F.; Lamchin, M.; Lee, E.J.; Zhu, Y.; et al. An Introduction To Mid-Latitude Ecotone: Sustainability and Environmental Challenges. Sib. J. For. Sci. 2017, 53, 41–53. [Google Scholar]
  6. Song, C.; Kim, S.J.; Moon, J.; Lee, S.J.; Lee, W.; Kim, N.; Wang, S.W.; Lee, W.K. Classification of global land development phases by forest and GDP changes for appropriate land management in the mid-latitude. Sustainability 2017, 9, 1342. [Google Scholar] [CrossRef]
  7. Lamchin, M.; Lee, W.K.; Jeon, S.W.; Wang, S.W.; Lim, C.H.; Song, C.; Sung, M. Long-term trend and correlation between vegetation greenness and climate variables in Asia based on satellite data. Sci. Total Environ. 2018, 618, 1089–1095. [Google Scholar] [CrossRef] [PubMed]
  8. Shvidenko, A.; Buksha, I.; Krakovska, S.; Lakyda, P. Vulnerability of Ukrainian forests to climate change. Sustainability 2017, 9, 1152. [Google Scholar] [CrossRef]
  9. Choi, Y.; Lim, C.H.; Ryu, J.; Jeon, S.W. Bioclimatic classification of Northeast Asia reflecting social factors: Development and characterization. Sustainability 2017, 9, 1137. [Google Scholar] [CrossRef]
  10. Ueyama, M.; Ichii, K.; Hirata, R.; Takagi, K.; Asanuma, J.; MacHimura, T.; Nakai, Y.; Ohta, T.; Saigusa, N.; Takahashi, Y.; et al. Simulating carbon and water cycles of larch forests in East Asia by the BIOME-BGC model with AsiaFlux data. Biogeosciences 2010, 7, 959–977. [Google Scholar] [CrossRef] [Green Version]
  11. Park, M.S. Inter-Korean forest cooperation 1998-2012: A policy arrangement approach. Sustainability 2015, 7, 5241–5259. [Google Scholar] [CrossRef]
  12. Bae, J.S.; Joo, R.W.; Kim, Y.S. Forest transition in South Korea: Reality, path and drivers. Land Use Policy 2012, 29, 198–207. [Google Scholar] [CrossRef]
  13. Lee, S.G.; Choi, H.A.; Yoo, H.; Song, C.; Cha, S.; Bae, S.W.; Son, Y.; Lee, W.K. Restoration plan for degraded forest in the democratic people’s Republic of Korea considering suitable tree species and spatial distribution. Sustainability 2018, 10, 856. [Google Scholar] [CrossRef]
  14. Kang, S.; Choi, W. Forest cover changes in North Korea since the 1980s. Reg. Environ. Chang. 2014, 14, 347–354. [Google Scholar] [CrossRef]
  15. Kim, D.; Lim, C.H.; Song, C.; Lee, W.K.; Piao, D.; Heo, S.; Jeon, S. Estimation of future carbon budget with climate change and reforestation scenario in North Korea. Adv. Sp. Res. 2016, 58, 1002–1016. [Google Scholar] [CrossRef]
  16. Cui, G.; Lee, W.K.; Kim, D.; Lee, E.J.; Kwak, H.; Choi, H.A.; Kwak, D.A.; Jeon, S. Estimation of forest carbon budget from land cover change in South and North Korea between 1981 and 2010. J. Plant Biol. 2014, 57, 225–238. [Google Scholar] [CrossRef]
  17. Feng, L.; Jia, Z.; Li, Q. The dynamic monitoring of aeolian desertification land distribution and its response to climate change in northern China. Sci. Rep. 2016, 6, 1–10. [Google Scholar] [CrossRef]
  18. Tsogtbaatar, J. Deforestation and reforestation needs in Mongolia. For. Ecol. Manag. 2004, 201, 57–63. [Google Scholar] [CrossRef]
  19. Wenhua, L. Degradation and restoration of forest ecosystems in China. For. Ecol. Manag. 2004, 201, 33–41. [Google Scholar] [CrossRef]
  20. Huang, L.; Xiao, T.; Zhao, Z.; Sun, C.; Liu, J.; Shao, Q.; Fan, J.; Wang, J. Effects of grassland restoration programs on ecosystems in arid and semiarid China. J. Environ. Manag. 2013, 117, 268–275. [Google Scholar] [CrossRef]
  21. Schaphoff, S.; Reyer, C.P.O.; Schepaschenko, D.; Gerten, D.; Shvidenko, A. Tamm Review: Observed and projected climate change impacts on Russia’s forests and its carbon balance. For. Ecol. Manag. 2016, 361, 432–444. [Google Scholar] [CrossRef]
  22. Potapov, P.; Hansen, M.C.; Stehman, S.V.; Loveland, T.R.; Pittman, K. Combining MODIS and Landsat imagery to estimate and map boreal forest cover loss. Remote Sens. Environ. 2008, 112, 3708–3719. [Google Scholar] [CrossRef]
  23. Matsushita, B.; Tamura, M. Integrating remotely sensed data with an ecosystem model to estimate net primary productivity in East Asia. Remote Sens. Environ. 2002, 81, 58–66. [Google Scholar] [CrossRef]
  24. Takagi, K.; Hirata, R.; Ide, R.; Ueyama, M.; Ichii, K.; Saigusa, N.; Hirano, T.; Asanuma, J.; Li, S.G.; Machimura, T.; et al. Spatial and seasonal variations of CO2flux and photosynthetic and respiratory parameters of larch forests in East Asia. Soil Sci. Plant Nutr. 2015, 61, 61–75. [Google Scholar] [CrossRef] [Green Version]
  25. Ueyama, M.; Kai, A.; Ichii, K.; Hamotani, K.; Kosugi, Y.; Monji, N. The sensitivity of carbon sequestration to harvesting and climate conditions in a temperate cypress forest: Observations and modeling. Ecol. Model. 2011, 222, 3216–3225. [Google Scholar] [CrossRef]
  26. Kondo, M.; Ichii, K.; Ueyama, M.; Mizoguchi, Y.; Hirata, R.; Saigusa, N. The role of carbon flux and biometric observations in constraining a terrestrial ecosystem model: A case study in disturbed forests in East Asia. Ecol. Res. 2013, 28, 893–905. [Google Scholar] [CrossRef]
  27. Churkina, G.; Running, S.W.; Schloss, A.L. Comparing global models of terrestrial net primary productivity (NPP): The importance of water availability. Glob. Chang. Biol. 1999, 5, 46–55. [Google Scholar] [CrossRef]
  28. Li, B.; Wang, W. Trade-offs and synergies in ecosystem services for the Yinchuan Basin in China. Ecol. Indic. 2018, 84, 837–846. [Google Scholar] [CrossRef]
  29. Merganičová, K.; Pietsch, S.A.; Hasenauer, H. Testing mechanistic modeling to assess impacts of biomass removal. For. Ecol. Manag. 2005, 207, 37–57. [Google Scholar] [CrossRef]
  30. Choi, S.; Lee, W.K.; Kwak, H.; Kim, S.R.; Yoo, S.; Choi, H.A.; Park, S.; Lim, J.H. Vulnerability assessment of forest ecosystem to climate change in Korea using MC1 model. J. For. Plan. 2011, 16, 149–161. [Google Scholar]
  31. Pietsch, S.A. Modelling Ecosystem Pools and Fluxes. Implementation and Application of Biogeochemical Ecosystem Models; Institute of Silviculture Department of Forest and Soil Sciences University of Natural Resources and Life Sciences (BOKU): Vienna, Austria, 2014. [Google Scholar]
  32. Thornton, P.E.; Law, B.E.; Gholz, H.L.; Clark, K.L.; Falge, E.; Ellsworth, D.S.; Goldstein, A.H.; Monson, R.K.; Hollinger, D.; Falk, M.; et al. Modeling and measuring the effects of disturbance history and climate on carbon and water budgets in evergreen needleleaf forests. Agric. For. Meteorol. 2002, 113, 185–222. [Google Scholar] [CrossRef]
  33. Pietsch, S.A.; Hasenauer, H.; Thornton, P.E. BGC-model parameters for tree species growing in central European forests. For. Ecol. Manag. 2005, 211, 264–295. [Google Scholar] [CrossRef]
  34. Akujärvi, A.; Shvidenko, A.; Pietsch, S.A. Modelling the impacts of intensifying forest management on carbon budget across a long latitudinal gradient in Europe. Environ. Res. Lett. 2018, 14, 034012. [Google Scholar] [CrossRef]
  35. Ishihara, M.I.; Suzuki, S.N.; Nakamura, M.; Enoki, T.; Fujiwara, A.; Hiura, T.; Hoshino, D.; Homma, K.; Hoshizaki, K.; Ida, H.; et al. Forest stand structure, composition, and dynamics in 34 sites over Japan. Ecol. Res. 2011, 26, 1007–1008. [Google Scholar] [CrossRef]
  36. Chae, N. Annual variation of soil respiration and precipitation in a temperate gorest (Quercus serrata and Carpinus laxiflora) under East Asian monsoon climate. J. Plant Biol. 2011, 54, 101–111. [Google Scholar] [CrossRef]
  37. Wang, H.M.; Saigusa, N.; Zu, Y.G.; Wang, W.J.; Yamamoto, S.; Kondo, H. Carbon fluxes and their response to environmental variables in a Dahurian larch forest ecosystem in northeast China. J. For. Res. 2008, 19, 1–10. [Google Scholar] [CrossRef]
  38. Noh, N.J.; Kim, C.; Bae, S.W.; Lee, W.K.; Yoon, T.K.; Muraoka, H.; Son, Y. Carbon and nitrogen dynamics in a Pinus densiflora forest with low and high stand densities. J. Plant Ecol. 2013, 6, 368–379. [Google Scholar] [CrossRef] [Green Version]
  39. Crowther, T.W.; Glick, H.B.; Covey, K.R.; Bettigole, C.; Maynard, D.S.; Thomas, S.M.; Smith, J.R.; Hintler, G.; Duguid, M.C.; Amatulli, G.; et al. Mapping tree density at a global scale. Nature 2015, 525, 201–205. [Google Scholar] [CrossRef]
  40. Kummu, M.; Varis, O. The world by latitudes: A global analysis of human population, development level and environment across the north-south axis over the past half century. Appl. Geogr. 2011, 31, 495–507. [Google Scholar] [CrossRef]
  41. Tran, H.; Iida, S.; Inoue, S. Initial sprouting of some tree species in natural forests following 9-month cutting in Kasuya Research Forest of Kyushu University Forest, Japan. J. Fac. Agric. Kyushu Univ. 2005, 50, 471–480. [Google Scholar]
  42. Lim, H.J.; Lee, Y.H.; Kwon, H. Evaluation of Community Land Model version 3.5-Dynamic global vegetation model over deciduous forest in Gwangneung, Korea. Korean J. Agric. For. Meteorol. 2010, 12, 95–106. [Google Scholar] [CrossRef]
  43. Kwon, H.; Lee, J.; Lee, Y.; Lee, J.W.; Jung, S.W.; Kim, J. Seasonal variations of evapotranspiration observed in a mixed forest in the Seolmacheon catchment. Korean J. Agric. For. Meteorol. 2009, 11, 39–47. [Google Scholar] [CrossRef]
  44. Shin, H.J.; Park, G.A.; Park, M.; Kim, S.J. Projection of forest vegetation change by applying future climate change scenario MIROC3.2 A1B. J. Korean Assoc. Geogr. Inf. Stud. 2012, 15, 64–75. [Google Scholar] [CrossRef]
  45. Piao, D.; Zhu, Y.; Piao, D.; Lee, W.; Zhu, Y.; Kim, M.; Song, C. Assessment of forest degradation and carbon storage for REDD+ project in North Korea. Korean J. Environ. Biol. 2016, 9999, 1–7. [Google Scholar] [CrossRef]
  46. Shvidenko, A.; Schepaschenko, D.; Nilsson, S.; Bouloui, Y. Tables and Models of Growth and Productivity of Forests of Major Forest Forming Species of Northern Eurasian; Federal Agency of Forest Management; International Institute for Applied Systems Analysis: Moscow, Russia, 2008. [Google Scholar]
  47. Noh, N.J.; Son, Y.; Lee, S.K.; Yoon, T.K.; Seo, K.W.; Kim, C.; Lee, W.K.; Bae, S.W.; Hwang, J. Influence of stand density on soil CO2efflux for a Pinus densiflora forest in Korea. J. Plant Res. 2010, 123, 411–419. [Google Scholar] [CrossRef] [PubMed]
  48. Eum, S.; Kang, S.; Lee, D. A simulation study to investigate climatic controls on net primary production (NPP) of a rugged forested landscape in The Mid-western Korean Peninsula. Korean J. Agric. For. Meteorol. 2005, 7, 67–78. [Google Scholar]
  49. Hanns Seidel Foundation German-North Korean Cooperation in the Field of Forestry. Available online: https://korea.hss.de/en/news/detail/german-north-korean-cooperation-in-the-field-of-forestry-news3931/ (accessed on 12 January 2019).
  50. Zhu, J.; Zhi-ping, F.; Zeng, D.; Jiang, F.; Matsuzaki, T. Comparison of stand structure and growth between artificial and natural forests of Pinus sylvestiris var.mongolica on sandy land. J. For. Res. 2003, 14, 103–111. [Google Scholar]
  51. Chen, Z.; Zhang, X.; Cui, M.; He, X.; Ding, W.; Peng, J. Tree-ring based precipitation reconstruction for the forest-steppe ecotone in northern Inner Mongolia, China and its linkages to the Pacific Ocean variability. Glob. Planet. Chang. 2012, 86–87, 45–56. [Google Scholar] [CrossRef]
  52. Khishigjargal, M.; Dulamsuren, C.; Leuschner, H.H.; Leuschner, C.; Hauck, M. Climate effects on inter- and intra-annual larch stemwood anomalies in the Mongolian forest-steppe. Acta Oecol. 2014, 55, 113–121. [Google Scholar] [CrossRef]
  53. Dulamsuren, C.; Khishigjargal, M.; Leuschner, C.; Hauck, M. Response of tree-ring width to climate warming and selective logging in larch forests of the Mongolian Altai. J. Plant Ecol. 2014, 7, 24–38. [Google Scholar] [CrossRef]
  54. Dulamsuren, C.; Hauck, M.; Leuschner, H.H.; Leuschner, C. Climate response of tree-ring width in Larix sibirica growing in the drought-stressed forest-steppe ecotone of northern Mongolia. Ann. For. Sci. 2011, 68, 275–282. [Google Scholar] [CrossRef] [Green Version]
  55. Kim, M.; Lee, W.K.; Choi, G.M.; Song, C.; Lim, C.H.; Moon, J.; Piao, D.; Kraxner, F.; Shvidenko, A.; Forsell, N. Modeling stand-level mortality based on maximum stem number and seasonal temperature. For. Ecol. Manag. 2017, 386, 37–50. [Google Scholar] [CrossRef] [Green Version]
  56. Holland, E.A.; Dentener, F.J.; Bobby, H.; Sulzman, J.M. Contemporary and pre-industrial global reactive nitrogen budgets. Biogeochemistry 1999, 46, 7–43. [Google Scholar] [CrossRef]
  57. Vet, R.; Artz, R.S.; Carou, S.; Shaw, M.; Ro, C.U.; Aas, W.; Baker, A.; Bowersox, V.C.; Dentener, F.; Galy-Lacaux, C.; et al. A global assessment of precipitation chemistry and deposition of sulfur, nitrogen, sea salt, base cations, organic acids, acidity and pH, and phosphorus. Atmos. Environ. 2014, 93, 3–100. [Google Scholar] [CrossRef]
  58. Granhall, U.; Lindberg, T. Nitrogen fixation in some coniferous forest ecosystems. Ecol. Bull. 1978, 26, 178–192. [Google Scholar]
  59. Todd, R.L.; Meyer, R.D.; Waide, J.B. Nitrogen Fixation in a Deciduous Forest in the South-Eastern United States. Source Ecol. Bull. 1978, 26, 114–121. [Google Scholar]
  60. Son, Y. Non-symbiotic nitrogen fixation in forest ecosystems. Ecol. Res. 2001, 16, 183–196. [Google Scholar] [CrossRef]
  61. Shugalei, L.S.; Vedrova, E.F. Nitrogen pool in northern-Taiga larch forests of central Siberia. Biol. Bull. 2004, 31, 200–208. [Google Scholar] [CrossRef]
  62. White, M.A.; Thornton, P.E.; Running, S.W.; Nemani, R.R. Parameterization and sensitivity analysis of the BIOME–BGC Terrestrial Ecosystem Model: Net primary production controls. Earth Interact. 2000, 4, 1–85. [Google Scholar] [CrossRef]
  63. Thornton, P.E.; Running, S.W. An improved algorithm for estimating incident daily solar radiation from measurements of temperature, humidity, and precipitation. Agric. For. Meteorol. 1999, 93, 211–228. [Google Scholar] [CrossRef]
  64. Warszawski, L.; Frieler, K.; Huber, V.; Piontek, F.; Serdeczny, O.; Schewe, J. The Inter-Sectoral Impact Model Intercomparison Project (ISI–MIP): Project framework. Proc. Natl. Acad. Sci. USA 2014, 111, 3228–3232. [Google Scholar] [CrossRef] [PubMed]
  65. Hempel, S.; Frieler, K.; Warszawski, L.; Schewe, J.; Piontek, F. A trend-preserving bias correction –The ISI-MIP approach. Earth Syst. Dyn. 2013, 4, 219–236. [Google Scholar] [CrossRef]
  66. Rosenzweig, C.; Elliott, J.; Deryng, D.; Ruane, A.C.; Müller, C.; Arneth, A.; Boote, K.J.; Folberth, C.; Glotter, M.; Khabarov, N.; et al. Assessing agricultural risks of climate change in the 21st century in a global gridded crop model intercomparison. Proc. Natl. Acad. Sci. USA 2014, 111, 3268–3273. [Google Scholar] [CrossRef] [PubMed]
  67. Kira, T. A New Classification of Climate in Eastern Asia as the Basis for Agricultural Geography; Horticultural Institute, Kyoto University: Kyoto, Japan, 1945. [Google Scholar]
  68. Choi, S.; Lee, W.; Kwak, D.; Lee, S.; Son, Y.; Lim, J.; Saborowski, J. Predicting forest cover changes in future climate using hydrological and thermal indices in South Korea. Clim. Res. 2011, 49, 229–245. [Google Scholar] [CrossRef] [Green Version]
  69. Inoue, T.; Enoki, T.; Tashiro, N.; Sakuta, K.; Inoue, S. Effects of topography and planted trees on the distribution of naturally regenerated broad-leaved trees in a 140-year-old Cryptomeria japonica plantation in northern Kyushu, Japan. J. For. Res. 2008, 13, 365–371. [Google Scholar] [CrossRef]
  70. Son, Y.M.; Kim, R.H.; Lee, K.H.; Pyo, J.K.; Kim, S.W.; Hwang, J.S.; Lee, S.J. Carbon Emission Factors and Biomass Allometric Equations by Species in Korea; National Institute of Forest Science: Seoul, Korea, 2014. [Google Scholar]
  71. Lim, J.; Shin, J.H.; Jin, G.Z.; Chun, J.H.; Oh, J.S. Forest stand structure, site characteristics and carbon budget of the Kwangneung Natural Forest in Korea. Korean J. Agric. For. Meteorol. 2003, 5, 101–109. [Google Scholar]
  72. Shin, H.J.; Park, M.J.; Kim, S.J. Evaluation of forest watershed hydro-ecology using measured data and RHESSys Model 2014for the Seolmacheon catchment. J. Korea Water Resour. Assoc. 2012, 45, 1293–1307. [Google Scholar] [CrossRef]
  73. Shi, F.; Chen, X.; Wang, W.; Takagi, K.; Akibayashi, Y.; Sasa, K.; Uemura, S. Vegetation characteristics of a Larch-dominant Site for CO2 flux monitoring study at the laoshan experimental station in northesast China. Eurasian J. For. Res. 2001, 3, 55–67. [Google Scholar]
  74. MOD17A3 MODIS/Terra Gross Primary Productivity Yearly L4 Global 1 km SIN Grid; NASA EOSDIS Land Processes DAAC. 2011. Available online: https://lpdaac.usgs.gov/products/mod17a3v055/ (accessed on 10 February 2019).
  75. Pietsch, S.A.; Hasenauer, H. Evaluating the self-initialization procedure for large-scale ecosystem models. Glob. Chang. Biol. 2006, 12, 1658–1669. [Google Scholar] [CrossRef]
  76. Kajimoto, T.; Matsuura, Y.; Sofronov, M.A.; Volokitina, A.V.; Mori, S.; Osawa, A.; Abaimov, A.P. Above- and belowground biomass and net primary productivity of a Larix gmelinii stand near Tura, central Siberia. Tree Physiol. 1999, 19, 815–822. [Google Scholar] [CrossRef]
  77. Kharuk, V.I.; Ranson, K.J.; Im, S.T.; Dvinskaya, M.L. Response of Pinus sibirica and Larix sibirica to climate change in southern Siberian alpine forest-tundra ecotone. Scand. J. For. Res. 2009, 24, 130–139. [Google Scholar] [CrossRef]
  78. Kang, S.; Lim, J.H.; Kim, E.S.; Cho, N. Modelling analysis of climate and soil depth effects on pine tree dieback in Korea using BIOME-BGC. Korean J. Agric. For. Meteorol. 2016, 18, 242–252. [Google Scholar] [CrossRef]
  79. Kang, S.; Lee, D.; Lee, J.; Running, S.W. Topographic and climatic controls on soil environments and net primary production in a rugged temperate hardwood forest in Korea. Ecol. Res. 2006, 21, 64–74. [Google Scholar] [CrossRef]
  80. Wang, P.; Sun, R.; Hu, J.; Zhu, Q.; Zhou, Y.; Li, L.; Chen, J.M. Measurements and simulation of forest leaf area index and net primary productivity in Northern China. J. Environ. Manag. 2007, 85, 607–615. [Google Scholar] [CrossRef] [PubMed]
  81. Gower, S.T.; Krankina, O.; Olson, R.J.; Apps, M.; Linder, S.; Wang, C. Net primary production and carbon allocation patterns of boreal forest ecosystems. Ecol. Appl. 2001, 11, 1395–1411. [Google Scholar] [CrossRef]
  82. Sato, H.; Kobayashi, H.; Iwahana, G.; Ohta, T. Endurance of larch forest ecosystems in eastern Siberia under warming trends. Ecol. Evol. 2016, 6, 5690–5704. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Shuman, J.K.; Shugart, H.H.; O’Halloran, T.L. Sensitivity of Siberian larch forests to climate change. Glob. Chang. Biol. 2011, 17, 2370–2384. [Google Scholar] [CrossRef]
  84. Jiang, H.; Apps, M.J.; Zhang, Y.; Peng, C.; Woodard, P.M. Modelling the spatial pattern of net primary productivity in Chinese forests. Ecol. Model. 1999, 122, 275–288. [Google Scholar] [CrossRef]
Forests 10 00523 g001 550
Figure 1. Study sites (oak, pine, and larch stands) along the mid-latitude ecotone (MLE) of Northeast Asia.
Figure 1. Study sites (oak, pine, and larch stands) along the mid-latitude ecotone (MLE) of Northeast Asia.
Forests 10 00523 g001
Forests 10 00523 g002 550
Figure 2. Climatic characteristics of the northernmost site (Buryatia, Russia) and the southernmost site (Kasuya, Japan) in 2005.
Figure 2. Climatic characteristics of the northernmost site (Buryatia, Russia) and the southernmost site (Kasuya, Japan) in 2005.
Forests 10 00523 g002
Forests 10 00523 g003 550
Figure 3. Growing stock volume (GSV) estimation through BGC-MAN simulations from 1960 to 2005.
Figure 3. Growing stock volume (GSV) estimation through BGC-MAN simulations from 1960 to 2005.
Forests 10 00523 g003
Forests 10 00523 g004 550
Figure 4. Net primary productivity (NPP) estimation through BGC-MAN from 1960 to 2005.
Figure 4. Net primary productivity (NPP) estimation through BGC-MAN from 1960 to 2005.
Forests 10 00523 g004
Forests 10 00523 g005 550
Figure 5. Relationship between the net primary productivities (NPPs) from MOD17A3H and BGC-MAN.
Figure 5. Relationship between the net primary productivities (NPPs) from MOD17A3H and BGC-MAN.
Forests 10 00523 g005
Table
Table 1. Site-specific information, including mortality rates, fire fractions, and applied management options.
Table 1. Site-specific information, including mortality rates, fire fractions, and applied management options.
ClassesOak Forest StandsPine Forest StandsLarch Forest Stands
SpeciesQuercus acutaQuercus serrataQuercus variabilisQuercus mongolicaQuercus mongolicaPinus densifloraPinus densifloraPinus sylvestrisPinus sylvestrisPinus sylvestrisLarix gmeliniiLarix sibiricaLarix sibiricaLarix gmeliniiLarix gmelinii
TypesEvergreenDeciduousDeciduousDeciduousDeciduousEvergreenEvergreenEvergreenEvergreenEvergreenDeciduousDeciduousDeciduousDeciduousDeciduous
LocationKasuyaGwang-neungSeolma-cheonSohungBikinskyGwang-neungWonsanHong-huaerjiOblastPribay-kalskyLaoshanTsenkher-mandalKhenteiYeravnin-skyBuryatia
Nation 1JPNROKROKDPRKRUSROKDPRKCHNRUSRUSCHNMNGMNGRUSRUS
Longitude (°)130.550127.149126.955126.032134.373127.162127.320119.993105.334107.985127.340109.068107.316111.411111.728
Latitude (°)33.65037.75037.93938.39146.77037.74839.14948.20052.18552.50445.20047.83549.10052.39754.960
Elevation (m)355340293169160425481836730620370179710219221612
Aspect (%)ENWNWNSSWNSWNESWSWNNWNWN
Slope (%)151415211519211262562483234
Sand (%)425052421550428965553731453438
Silt (%)36132536484536625374536313331
Clay (%)223723223752251081833242321
Soil depth (m)1.001.001.001.001.601.001.001.000.800.951.000.901.001.000.60
East horizon (°)5072632003163531616185
West horizon (°)601805410560101024802
Forest age14080–20020–40-8080+-40–60956040–5011060–150100100
Mortality (%) 21.7–2.01.7–2.02.0+2.0+1.7–2.01.4+1.4–3.01.4–3.00.70.8–2.05.0–5.60.20.2–0.30.2–0.30.2–0.3
Fire fraction (%)----0.005---0.250.15-0.15-0.080.08
Options 3P, NP, NC, P, NC, P, LP, NC, P, NC, P, LP, NP, NP, NC, P, NP, NP, NP, NP, N
1 Country code ROK, DPRK, JPN, CHN, MNG, and RUS represent South Korea, North Korea, Japan, China, Mongolia, and Russia, respectively. 2 Mainly based on Shvidenko et al. (2008) and Kim et al. (2017), but also followed site-specific information. 3 Options represent historic intervention in certain forest stands (P: Planting, C: Clear-cutting, N: Naturally preserved, L: Logging).
Table
Table 2. Species characteristics for running BioGeoChemistry Management (BGC-MAN).
Table 2. Species characteristics for running BioGeoChemistry Management (BGC-MAN).
ParametersUnitsOakPineLarch
Annual leaf and fine root turnover fractionyear−11.00.31.0
Annual live wood turnover fractionyear−10.70.30.7
Annual whole-plant mortality fractionyear−10.0080.1800.050
Annual fire mortality fractionyear−10.0170.7000.056
New fine root C:new leaf Cratio1.00.5230.8~1.2
New stem C:new total wood Cratio1.292.51.2~2.2
New live wood C:new total wood Cratio0.1200.0590.1
New coarse root C:new stem Cratio0.2500.2900.23
Current growth proportionratio0.50.50.5
C:N of leaveskgC/kgN26.933.125.8~27
C:N of leaf letterskgC/kgN63.3132.0111.9
C:N of fine rootskgC/kgN73.538.042.0
C:N of live woodkgC/kgN63.550.042.0
C:N of dead woodkgC/kgN450.01400.0442.0
Leaf litter labile proportionDIM0.2000.2570.390
Leaf litter cellulose proportionDIM0.5600.4930.440
Leaf litter lignin proportionDIM0.2400.2500.170
Fine root labile proportionDIM0.3400.2520.300
Fine root cellulose proportionDIM0.4400.4930.450
Fine root lignin proportionDIM0.2200.2530.250
Dead wood cellulose proportionDIM0.7040.7100.760
Dead wood lignin proportionDIM0.2960.2900.240
Canopy water interception coefficient1/LAI/d0.0380.0510.041
Canopy light extinction coefficientDIM0.5400.5100.500
All-sided to projected leaf area ratioDIM2.02.62.6
Canopy average specific leaf aream2/kgC35.013.020.3
Ratio of shaded SLADIM2.02.02.0
Fraction of leaf N in RubiscoDIM0.08800.04570.0750
Maximum stomatal conductancem/s0.00180.00100.0022
Cuticular conductancem/s0.000040.0000140.00001
Boundary layer conductancem/s0.0050.090.008
Leaf water potentialMPa(−0.1)–(−3.5)(−0.5)–(−2.2)(−0.7)–(−2.6)
Vapor pressure deficitMPa200–255050–2500800–3200
Night time freezing temperatureC(−1)–(−8)(−2)–(−8)(−8)–(−20)
Table
Table 3. Comparing BGC-MAN simulations and previous literature.
Table 3. Comparing BGC-MAN simulations and previous literature.
TypesSimulationLiteratureSourcesSpeciesLocations
GSV
(m3 ha−1)		410109–516Estimation [35,69,70]OakKasuya
511550Chae (2011) [36]OakGwangneung
15478–159Kwon et al. (2009) [43]OakSeolmacheon
6556–66Piao et al. (2016) [45]OakSohung
153154Shvidenko et al. (2008) [46]OakBikinsky
391378–645Noh et al. (2013) [38]PineGwangneung
146117–173Zhu et al. (2003) [50]PineHonghuaerji
282282Shvidenko et al. (2008) [46]PineOblast
124102Shvidenko et al. (2008) [46]PinePribaykalsky
127118–133Shi et al. (2001) [73]LarchLaoshan
227227Shvidenko et al. (2008) [46]LarchYeravninsky
116117Shvidenko et al. (2008) [46]LarchBuryatia
NPP
(Mg C ha−1)		4.284.30–6.05Lim et al. (2003) [71]OakGwangneung
5.394.00–6.55Lim et al. (2010) [42]OakGwangneung
3.705.05Shvidenko et al. [46]OakBikinsky
6.395.37–7.20Eum et al. (2005) [48]PineGwangneung
6.754.30–6.03Cui et al. (2014) [16]PineWonsan
3.493.06Shvidenko et al. (2008) [46]PineOblast
2.502.72Shvidenko et al. (2008) [46]PinePribaykalsky
5.237.18–7.89Kondo et al. (2013) [26]LarchLaoshan
1.283.38Shvidenko et al. (2008) [46]LarchYeravninsky
0.492.81Shvidenko et al. (2008) [46]LarchBuryatia
GPP8.7211.93–12.06Shin et al. (2012) [72]OakSeolmacheon
(Mg C ha−1)2.855.62–6.10Takagi et al. (2015) [24]LarchKhentei

Share and Cite

MDPI and ACS Style

Song, C.; Pietsch, S.A.; Kim, M.; Cha, S.; Park, E.; Shvidenko, A.; Schepaschenko, D.; Kraxner, F.; Lee, W.-K. Assessing Forest Ecosystems across the Vertical Edge of the Mid-Latitude Ecotone Using the BioGeoChemistry Management Model (BGC-MAN). Forests 2019, 10, 523. https://doi.org/10.3390/f10060523

AMA Style

Song C, Pietsch SA, Kim M, Cha S, Park E, Shvidenko A, Schepaschenko D, Kraxner F, Lee W-K. Assessing Forest Ecosystems across the Vertical Edge of the Mid-Latitude Ecotone Using the BioGeoChemistry Management Model (BGC-MAN). Forests. 2019; 10(6):523. https://doi.org/10.3390/f10060523

Chicago/Turabian Style

Song, Cholho, Stephan A. Pietsch, Moonil Kim, Sungeun Cha, Eunbeen Park, Anatoly Shvidenko, Dmitry Schepaschenko, Florian Kraxner, and Woo-Kyun Lee. 2019. "Assessing Forest Ecosystems across the Vertical Edge of the Mid-Latitude Ecotone Using the BioGeoChemistry Management Model (BGC-MAN)" Forests 10, no. 6: 523. https://doi.org/10.3390/f10060523

APA Style

Song, C., Pietsch, S. A., Kim, M., Cha, S., Park, E., Shvidenko, A., Schepaschenko, D., Kraxner, F., & Lee, W. -K. (2019). Assessing Forest Ecosystems across the Vertical Edge of the Mid-Latitude Ecotone Using the BioGeoChemistry Management Model (BGC-MAN). Forests, 10(6), 523. https://doi.org/10.3390/f10060523

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop