Phenological Classiﬁcation Using Deep Learning and the Sentinel-2 Satellite to Identify Priority Afforestation Sites in North Korea

: The role of forests to sequester carbon is considered an important strategy for mitigating climate change and achieving net zero emissions. However, forests in North Korea have continued to be cleared since the 1990s due to the lack of food and energy resources. Deforestation in this country has not been accurately classiﬁed nor consistently reported because of the characteristics of small patches. This study precisely determined the area of deforested land in North Korea through the vegetation phenological classiﬁcation using high-resolution satellite imagery and deep learning algorithms. Effective afforestation target sites in North Korea were identiﬁed with priority grade. The U-Net deep learning algorithm and time-series Sentinel-2 satellite images were applied to phenological classiﬁcation; the results reﬂected the small patch-like characteristics of deforestation in North Korea. Based on the phenological classiﬁcation, the land cover of the country was classiﬁed with an accuracy of 84.6%; this included 2.6 million ha of unstocked forest and reclaimed forest. Sites for afforestation were prioritized into ﬁve grades based on deforested characteristics, altitude and slope. Forest area is expanded and the forest ecosystem is restored through successful afforestation, this may improve the overall ecosystem services in North Korea. In the long term, it will be possible to contribute to carbon neutrality and greenhouse gas reduction on the Korean Peninsula level through optimal afforestation by using these outcomes.


Introduction
With accelerating climate change, many countries have declared their goal of carbon neutrality. The success of mitigation initiatives is dependent on the long-term net balance between carbon gains and losses [1]. From a carbon-neutral perspective, it is very important to estimate carbon uptake and emissions based on land cover. It is possible to estimate carbon emissions and absorption by establishing a land cover map and tracking changes through continuous mapping. Forests play an important role in the global carbon cycle and are thus valued worldwide for the array of ecosystem services they provide to society [2]; this includes the sequestration or emission of carbon, one of the major greenhouse gases [3]. Forests take up atmospheric carbon dioxide during photosynthesis, using the gas to produce sugar and other organic compounds [4]; this process facilitates carbon sequestration function. The amount of carbon uptake in forests varies depending on the climate, and forest physiognomy and age. For example, forests that are 30 years of age in South Korea take up more than 10 t/ha of carbon annually [5]. The forested area of South Korea is approximately 6.3 million ha [6], sequestering 67.91 million t of carbon. North Korea, also imagery [24]. However, if a single-period satellite image is used, the classification result may differ depending on the time when the image was acquired, even at the same location [25]. For satellite imagery acquired during winter is problematic as it does not recognize the difference between agricultural land and barren. The Korean Peninsula, where North Korea is located, experiences four distinct seasons, and where most land cover consists of vegetation; as such, the phenological characteristics of vegetation are very well represented [26]. If these characteristics are used during the classification process, it enables the compensation of errors in the classification results that occur with the use of single-period satellite imagery [25]. To accurately capture the vegetation phenological characteristics, satellite images of at least four or more scenes must be used to reflect all characteristics across all seasons; this necessitates the use of a large amount of data in the analysis. For analysis, the satellite image classification may be effectively and accurately carried out using an artificial intelligence deep learning U-Net algorithm; this has recently attracted considerable research attention [27,28].
This study employed U-Net-based deep learning for land cover classification, including deforestation, for North Korea by reflecting the phenological characteristics of vegetation and utilizing Sentinel-2 satellite imagery with 10 m spatial resolution. The usage of Sentinel-2 satellite imagery allows to accurately analyze the small patch-like deforestation occurring in North Korea. Through the derivation of satellite imagery from across North Korea that reflects the phenological characteristics of vegetation, it is possible to effectively classify land cover. However, there is a large amount of available data from the phenological Sentinel-2 satellite imagery spanning North Korea, forming a type of Big Data. Deep learning has been hugely successful in Big Data analysis; as such, it may be a good candidate to capitalize on the potential of such complex and massive data [29]. The U-Net deep learning algorithm performs well in the classification of satellite imagery [30]. As the exact deforested area in North Korea has not yet been derived, and forest restoration has not been consistently effective, it is necessary to identify exact deforestation status and establish an evidence-based forest restoration plan for recovering ecosystem services and forest environment [31,32]. This means that optimal sites for afforestation in North Korea need to be derived by using accurate deforestation data. Effective forest restoration, based on graded afforestation target data, should be promoted to restore and manage forest resources in the country [33]. This data may form the basis forest cooperation between the international community and North Korea.

Study Area
North Korea is located in East Asia, between 37 • 41 00"N and 43 • 00 36"N (latitude), and 124 • 18 41"E and 130 • 41 32"E (longitude); it occupies the northern part of the Korean Peninsula ( Figure 1). North Korea spans an area of 12,054,000 ha, in which 80% of its territory is mountainous. The administrative districts of the country consist of nine provinces, one capital city, and three special cities. The country has an average altitude above sea level of~440 m, including plateaus above 2000 m in altitude. The terrain in North Korea may be characterized as being high in the east and low in the west. The country is located in a mid-latitude ecotone, characterized by a temperate climate and four distinct seasons; this has created clearly observable phenological characteristics in the vegetation of North Korea. In North Korean, rice planting and sowing is usually come in June to July. And then crops are harvested in September to October. Phenological characteristics of vegetation is appeared by distinct season and agricultural characteristics.
Previously, the mountainous regions of North Korea were forested, accounting for 60% of forested area in the entire Korean Peninsula; however, the current forest cover has been decreasing [34]. The principal reason for this deforestation is the lack of energy and food resources [35]. In North Korea, the lack of an established energy infrastructure forces the residents to use forest resources, including firewood, for energy. The lack of food also creates the demand for more agricultural land, in which forested land is reclaimed after felling timber for firewood. For this reason, North Korea has been undergoing gradual deforestation, beginning with the forests in low altitude areas. The Food and Agriculture Organization of the United Nations (FAO) has reported that since the 1990s, 12% of the forests in North Korea have been cleared [36]. North Korea has been attempting to resolve its energy and food issues to protect its forest resources [36]; as a part of this effort, they conducted 'Agroforestry Management,' to reclaim land with slopes < 15 • and afforest land on slopes > 15 • in the mountain region. Since 2012 when Kim Jong-un commenced his reign, North Korea has considered that forest-related policies have sufficient importance to consider these issues a "forest restoration battle;" as such, there has been emphasis for these issues to be included in various propaganda materials [35,36].
Although North Korea claims that afforestation is being implemented, its rate of deforestation is much faster than afforestation. In attempting to identify the results of these afforestation efforts using Google Earth Pro (GEP), this study identified a small area in Pyongyang special city ( Figure 2). The area of the afforested patch in Figure 2 is approximately 3 ha and its altitude ranges from 110 to 158 m; the surrounding area was barren, reclaimed forest.
x FOR PEER REVIEW 4 of 15 Previously, the mountainous regions of North Korea were forested, accounting for 60% of forested area in the entire Korean Peninsula; however, the current forest cover has been decreasing [34]. The principal reason for this deforestation is the lack of energy and food resources [35]. In North Korea, the lack of an established energy infrastructure forces the residents to use forest resources, including firewood, for energy. The lack of food also creates the demand for more agricultural land, in which forested land is reclaimed after felling timber for firewood. For this reason, North Korea has been undergoing gradual deforestation, beginning with the forests in low altitude areas. The Food and Agriculture Organization of the United Nations (FAO) has reported that since the 1990s, 12% of the forests in North Korea have been cleared [36]. North Korea has been attempting to resolve its energy and food issues to protect its forest resources [36]; as a part of this effort, they conducted 'Agroforestry Management,' to reclaim land with slopes < 15° and afforest land on slopes > 15° in the mountain region. Since 2012 when Kim Jong-un commenced his reign, North Korea has considered that forest-related policies have sufficient importance to consider these issues a "forest restoration battle;" as such, there has been emphasis for these issues to be included in various propaganda materials [35,36].
Although North Korea claims that afforestation is being implemented, its rate of de-  Previously, the mountainous regions of North Korea were forested, accounting for 60% of forested area in the entire Korean Peninsula; however, the current forest cover has been decreasing [34]. The principal reason for this deforestation is the lack of energy and food resources [35]. In North Korea, the lack of an established energy infrastructure forces the residents to use forest resources, including firewood, for energy. The lack of food also creates the demand for more agricultural land, in which forested land is reclaimed after felling timber for firewood. For this reason, North Korea has been undergoing gradual deforestation, beginning with the forests in low altitude areas. The Food and Agriculture Organization of the United Nations (FAO) has reported that since the 1990s, 12% of the forests in North Korea have been cleared [36]. North Korea has been attempting to resolve its energy and food issues to protect its forest resources [36]; as a part of this effort, they conducted 'Agroforestry Management,' to reclaim land with slopes < 15° and afforest land on slopes > 15° in the mountain region. Since 2012 when Kim Jong-un commenced his reign, North Korea has considered that forest-related policies have sufficient importance to consider these issues a "forest restoration battle;" as such, there has been emphasis for these issues to be included in various propaganda materials [35,36].
Although North Korea claims that afforestation is being implemented, its rate of deforestation is much faster than afforestation. In attempting to identify the results of these afforestation efforts using Google Earth Pro (GEP), this study identified a small area in Pyongyang special city ( Figure 2). The area of the afforested patch in Figure 2 is approximately 3 ha and its altitude ranges from 110 to 158 m; the surrounding area was barren, reclaimed forest.

Sentinel-2 Data Collection and Processing
Bi-monthly vegetation greenness data spanning North Korea is required for North Korea when using the phenological characteristics of vegetation; this data is made up of Sentinel-2 satellite imagery. This study uses the normalized difference vegetation index (NDVI) as the vegetation greenness indicator, in which the vegetation vitality was determined based on the difference between the red and near infrared (NIR) bands from the satellite imagery. As the spatial resolution of the red, green, blue, and NIR bands from the Sentinel-2 satellite is 10 m, this same spatial resolution was established for the NDVI derived from the Sentinel-2. The NDVI was calculated using Equation (1): To calculate the NDVI, basic data processing of Sentinel-2 satellite imagery was performed using the Google Earth Engine (GEE). The swath width of Sentinel-2 satellite imagery is 290 km, provided with swath shaped linear image, cut into squares. As one swath cannot cover full North Korea, several scenes were selected. Of the Sentinel-2 satellite imagery acquired by GEE after March 2018, the image with the least cloud cover was mosaiced. The NDVI image was calculated by applying the NDVI formula to the band 4 (red) and band 8 (NIR) of the selected satellite imagery. Then maximum value of NDVI is composited between selected satellite images for reflecting representative vegetation greenness. From this process, six bi-monthly NDVI images were generated from March, 2018 to February, 2019 using images with the fewest clouds ( Figure 3). Of the NDVI image results, some were excluded from analysis (i.e., Figure 3f), as the NDVI was not normally expressed in the snow-covered area. The NDVI outputs indicated differences in the seasonal NDVI values across North Korea. This is based on the difference in the seasonal NDVI characteristics of vegetation; Figure 3a,b are before rice planting and sowing, that some area has low value of NDVI. However, in Figure 3c, the general value of NDVI is high over the territory. This shows that the NDVI difference according to landcover is distinct within the scenes, during vegetation growth, sow, harvest, and leaf fall. Thus, Figure 3a-e were stacked into one multi-band NDVI image using the 'composite band' tool in ESRI ArcGIS software. Each band of result is comprised of Figure 3a-e. The multi-band NDVI image was able to reflect difference in the NDVI values for each scene with respect to the coordinates. Multi-band satellite image can be displayed via red-green-blue (RGB) color combination. Multi-band NDVI image also demonstrate color differences based on land cover according to NDVI value difference over time periods (Figure 4). In Figure 4, red area is supposed as agricultural land, because agricultural land has high value of NDVI only in Jul-Aug, 2018 band, after sowing. Land cover is identified using GEP by visually comparing Figure 4 and GEP satellite images according to color of Figure 4.

Classification Using U-Net Based Algorithm
The U-Net algorithm is a deep convolutional neural network (CNN), first introduced in biomedical image segmentation by Olaf et al. [37]; it has since been used in many segmentation tasks. U-Net, being one form of segmentation, can be used for land cover classification using satellite imagery. The U-Net algorithm is operated by supervised training, requiring ground truth labeled dataset for modeling. This study mapped the land cover of North Korea based on stepwise classification using the U-Net algorithm; this is a multiple, single-class segmentation.

Classification Using U-Net Based Algorithm
The U-Net algorithm is a deep convolutional neural network (CNN), first introduced in biomedical image segmentation by Olaf et al. [37]; it has since been used in many segmentation tasks. U-Net, being one form of segmentation, can be used for land cover classification using satellite imagery. The U-Net algorithm is operated by supervised training,

Classification Using U-Net Based Algorithm
The U-Net algorithm is a deep convolutional neural network (CNN), first introduce in biomedical image segmentation by Olaf et al. [37]; it has since been used in many seg mentation tasks. U-Net, being one form of segmentation, can be used for land cover clas sification using satellite imagery. The U-Net algorithm is operated by supervised training Five scenes of training data for U-Net were extracted from multi-band NDVI images, which were 20 × 20 km in size. Each scene is extracted from five different provinces to minimize data differences due to variations in latitude and longitude. Since there is no ground truth data, that could be referenced for generating labeled dataset, the 'Band thresholding' method was used to generate pseudo labeling data, which enables to use Remote Sens. 2021, 13, 2946 7 of 15 unlabeled data for training with a limited amount of labeled data. Among the five NDVI bands of multi-band NDVI images, one with binary NDVI distributions was used for thresholding (in the case of first classification, May-June 2018 layer was used, because that period is before sowing, so that NDVI difference between mountain region and nonmountain region is distinct). Threshold is set on the point of inflection or dividing the histogram into two according to Otsu's method [38]; then, the setting range was confirmed by visual interpretations on randomly distributed inspection points in the target area. Figure 5 presents the training data from Pyeongyang Province and the forest labeling process using L3Harris ENVI software; Figure 5a is the histogram of May-June 2018 layer of Figure 5b and band thresholding parameters, and Figure 5b,c are training data and labeled data. In Figure 5b, forests were distinguished by color as light yellow and fully covered by green color in Figure 5c. This process enabled U-Net model to learn seasonal NDVI characteristics of training data (in the case of Figure 5, NDVI characteristics of forest). ENVI software automatically separates input training data to 80% of training data and 20% of validation data and shuffle.
which were 20 × 20 km in size. Each scene is extracted from five different provinces to minimize data differences due to variations in latitude and longitude. Since there is no ground truth data, that could be referenced for generating labeled dataset, the 'Band thresholding' method was used to generate pseudo labeling data, which enables to use unlabeled data for training with a limited amount of labeled data. Among the five NDVI bands of multi-band NDVI images, one with binary NDVI distributions was used for thresholding (in the case of first classification, May-June 2018 layer was used, because that period is before sowing, so that NDVI difference between mountain region and nonmountain region is distinct). Threshold is set on the point of inflection or dividing the histogram into two according to Otsu's method [38]; then, the setting range was confirmed by visual interpretations on randomly distributed inspection points in the target area.  Figure 5c. This process enabled U-Net model to learn seasonal NDVI characteristics of training data (in the case of Figure 5, NDVI characteristics of forest). ENVI software automatically separates input training data to 80% of training data and 20% of validation data and shuffle. Although U-Net is a type of semantic segmentation that can perform multi-class classification, stepwise classification was conducted to improve the data processing efficiency. Figure 6 presents the stepwise classification process, in which the full scene is initially classified into mountain and non-mountain regions, and then further divided from nonmountain region to vegetation cover and non-vegetation cover. Vegetation cover may be categorized into agricultural region and reclaimed forest. The non-vegetation region was not clearly classified as the used area and water, as there were very small NDVI differences between these two categories. The forest was classified as either an evergreen or deciduous forest, in which the NDVI difference for leaf growth and winter season were elucidated using the Mar-Apr and Nov-Dec layers from Figure 3. Time series multi-band NDVI imagery were classified into seven scenes using a six-time classification and six different CNN models. Although U-Net is a type of semantic segmentation that can perform multi-class classification, stepwise classification was conducted to improve the data processing efficiency. Figure 6 presents the stepwise classification process, in which the full scene is initially classified into mountain and non-mountain regions, and then further divided from nonmountain region to vegetation cover and non-vegetation cover. Vegetation cover may be categorized into agricultural region and reclaimed forest. The non-vegetation region was not clearly classified as the used area and water, as there were very small NDVI differences between these two categories. The forest was classified as either an evergreen or deciduous forest, in which the NDVI difference for leaf growth and winter season were elucidated using the Mar-Apr and Nov-Dec layers from Figure 3. Time series multi-band NDVI imagery were classified into seven scenes using a six-time classification and six different CNN models.

Identifying Priority Afforestation Sites Using Topographic Data
One of key geographic features in North Korea is the distribution of highlands at an altitude of 2000 m. These highlands occur at altitudes higher than the tree limit layer, and alpine vegetation occur above the tree limit layer. A previous study reported that alpine vegetation is distributed at altitudes between 1900 and 2000 m above sea level (asl) [39]. To extract the alpine region, Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM) was used; this has a spatial resolution of 1-ARC (~30 m). The SRTM DEM imagery were provided by USGS, by 1 • tiles. All tiles on the Korean Peninsula were mosaiced and clipped using the Korean Peninsula polygon ( Figure 1); regions above 1900 m were classified as plateau vegetation as per Kong et al. [39] and identified using GEP.

Identifying Priority Afforestation Sites Using Topographic Data
One of key geographic features in North Korea is the distribution of h altitude of 2000 m. These highlands occur at altitudes higher than the tree l alpine vegetation occur above the tree limit layer. A previous study repor vegetation is distributed at altitudes between 1900 and 2000 m above sea To extract the alpine region, Shuttle Radar Topography Mission (SRTM) Di Model (DEM) was used; this has a spatial resolution of 1-ARC (~30 m). Th imagery were provided by USGS, by 1° tiles. All tiles on the Korean P mosaiced and clipped using the Korean Peninsula polygon ( Figure 1); regio m were classified as plateau vegetation as per Kong et al. [39] and identifie The priority afforestation areas were identified from the classification areas reflected the most effective afforestation sites in North Korea. Effectiv criteria differed based various factors including the soil type, weather and f National Institute of Forest Science of South Korea has been monitoring deforested areas in North Korea every 10 years, beginning in 1999. They afforestation sites as sites a slope > 8°, as per the Agricultural Land Act of So North Korea reported that areas with a slope >15° were regarded as affor sites based on the agroforestry policy data of North Korea, procured international workshop [7]. Yu et al. [40] demonstrated that the soils in N relatively dry; as such, applying the slope > 8° criteria for target sites is de afforestation of North Korea [40]. Afforestation areas with slopes > 8° may with the past materials. Thus, deforested areas with slopes > 8° were extracte those criteria from SRTM DEM data. Deforested areas with slopes > extracted to evaluate the agroforestry policy of North Korea. For e restoration, the priorities of afforestation sites were established; these prio at five grades, and additionally difficult areas for reforestation were restoration targets. The afforestation priority was initially graded based deforestation, followed by slope. Each slope standard was establi reforestation-related studies [40,41] in Korea. The priority afforestation areas were identified from the classification results; these areas reflected the most effective afforestation sites in North Korea. Effective afforestation criteria differed based various factors including the soil type, weather and forest type. The National Institute of Forest Science of South Korea has been monitoring forested and deforested areas in North Korea every 10 years, beginning in 1999. They define target afforestation sites as sites a slope > 8 • , as per the Agricultural Land Act of South Korea [7]. North Korea reported that areas with a slope >15 • were regarded as afforestation target sites based on the agroforestry policy data of North Korea, procured during a 2015 international workshop [7]. Yu et al. [40] demonstrated that the soils in North Korea are relatively dry; as such, applying the slope > 8 • criteria for target sites is desirable for the afforestation of North Korea [40]. Afforestation areas with slopes > 8 • may be compared with the past materials. Thus, deforested areas with slopes > 8 • were extracted by applying those criteria from SRTM DEM data. Deforested areas with slopes > 15 • were also extracted to evaluate the agroforestry policy of North Korea. For effective forest restoration, the priorities of afforestation sites were established; these priorities were set at five grades, and additionally difficult areas for reforestation were set as natural restoration targets. The afforestation priority was initially graded based on the type of deforestation, followed by slope. Each slope standard was established through reforestation-related studies [40,41] in Korea.

Land Classification Results and Accuracy Verification
North Korean land cover was classified into five macro-classes and eight sub-classes (Figure 7), based on the classification results using the U-Net algorithm and the phenological characteristics of vegetation. Reclaimed forest was included in the agricultural area, as it is essentially used as agricultural land. Furthermore, unstocked forest is included in area of forest, as it was forest, but deforested.
To accurately verify the results, 990 points were randomly scattered and visually interpreted and identified using GEP data and Figure 4. The confusion matrix was derived from the visual interpretation results in Table 1. The classification accuracy was 84.8% and the kappa value was 0.7740; these exceed the standards for securing a 75% classification accuracy [42].
Based on the spatial distribution of the classification results, the agricultural land was largely distributed in the west. The plateau vegetation, extracted using DEM, was distributed in the Ryanggang and Hamgyong provinces, where the Kaema Plateau is located. Notably, the deforested areas were distributed at the boundary between forest and agricultural lands at low altitude. This spatial distribution reflects the characteristics of Remote Sens. 2021, 13, 2946 9 of 15 deforestation caused by the human reclamation of forest. Table 2 shows that the total area of agriculture is the sum of agricultural land and reclaimed forest, which is approximately 35% of North Korean territory. By contrast, 18% of South Korea is agricultural land, which means that North Korea is using a large area for agriculture due to food shortages. To accurately verify the results, 990 points were randomly scattered and visually i terpreted and identified using GEP data and Figure 4. The confusion matrix was deriv from the visual interpretation results in Table 1. The classification accuracy was 84.8% an the kappa value was 0.7740; these exceed the standards for securing a 75% classificatio accuracy [42]. Based on the spatial distribution of the classification results, the agricultural land w largely distributed in the west. The plateau vegetation, extracted using DEM, was distri uted in the Ryanggang and Hamgyong provinces, where the Kaema Plateau is locate Notably, the deforested areas were distributed at the boundary between forest and ag cultural lands at low altitude. This spatial distribution reflects the characteristics of defo estation caused by the human reclamation of forest. Table 2 shows that the total area agriculture is the sum of agricultural land and reclaimed forest, which is approximate 35% of North Korean territory. By contrast, 18% of South Korea is agricultural land, whi means that North Korea is using a large area for agriculture due to food shortages.  The classification results were used to calculate the deforestation ratio according to administrative districts, as shown in Table 3. The deforestation ratio in the southwest region including the Pyongan and Hwanghae provinces, was higher (maximum of 61%), compared to the northeast region including the Hamgyong, Chagang, and Ryanggang provinces. As the southwest region is an area with a high proportion of agricultural land, forests may have been reclaimed to further expand agricultural land.
The agricultural land and forest areas differed when comparing the results of this study and the land cover map from the Korean Ministry of Environment (Table 2). However, the total agricultural land and forest areas were very close (i.e., 11.3 and 11.5 million ha). This error may have occurred as researchers used different criteria to classify the land cover; the FAO reported that the forested area of North Korea was~9 million ha in 1987 [34], while the Korean Rural Economic Institute (KREI) reported that the agriculture area of North Korea is 1.96 million ha in 2009 [16]. As such, in terms of the reports from the FAO and KREI on land cover change over time, this study includes deforested areas within forest area, and the classification results are convincing.

Comparison with Previous Studies for the Deforested Area
The area of priority afforestation sites was derived using the classification results and compared with the Statistics Korea and Jin et al. [7,14] results in Table 4. As a result of the derivation, the total area of deforestation was 2.6 million ha; this is the same as the Statistics Korea result [7]. However, the Statistics Korea result was derived by applying slopes > 8 • , while the results from this study were derived without slope. Therefore, slope is considered in this study and compared against the Statistic Korea results, the total area of deforestation from this study is 1.7 million ha versus the 2.6 million ha [7]. For slopes > 15 • as per the agroforestry policy data of North Korea, the deforestation area was 1.1 and 1.4 million ha for this study and from Statistics Korea, respectively [7]. There is a difference between the results of these two studies; this may be due to the differences from the spatial resolution of the satellite images and the data acquisition time for images used. Statistics Korea classified land cover using single period RapidEye satellite imagery with 5 m resolution; compared with Jin et al. [14], the deforestation area is significantly different. Although this study and Jin et al. [14] classified using the phenological characteristics of vegetation, there was a large difference in the spatial resolution (i.e., 10 versus 250 m) [14]. Therefore, the area of deforested land was overestimated as the deforested land in North Korea is distributed in small patches. Based on the spatial distribution of the classification results, the spatial distribution of deforestation is also very similar. When comparing the results of the three studies, the classification area varies depending on the resolution, as the results of this study and that from Statistics Korea were similar [7]. This emphasized the importance of using high-resolution satellite imagery to effectively analyze the deforestation area in North Korea.

Identification of Priority Afforestation Target Site
The priority afforestation target sites were graded in Table 5 and mapped in Figure 8 in order ensure effective restoration and identify effective sites in North Korea. Reclaimed forest and unstocked forest of Figure 7 are graded using slope and altitude data derived from STRM DEM. According to the results of a study by the KREI for afforestation in North Korea, it is difficult to proceed with the project due to poor accessibility in areas with an altitude exceeding 600 m [41]. As such, it is necessary to prevent the further destruction of forests by human activity and restore forests through the natural transition process in deforested lands at altitudes > 600 m. Among the deforested areas in North Korea, the area with the first grade for afforestation were the unstocked forests at altitudes < 600 m. Unstocked forests make up less than 30% of the crown area and are suitable for reforestation as it is essentially a forest physiognomy where vegetation grows in a natural state. Compared to reclaimed forests, unstocked forests have better soil physical and chemical conditions as crops are not continuously grown in these areas. Afforestation priorities two to five were prioritized based on the slope of reclaimed forests at altitudes < 600 m; the higher the slope, the greater the priority, and the more difficult it is for crop cultivation which increases the risk of soil erosion. In North Korea, it is desirable to continue using reclaimed forests used as agricultural land as an agroforestry policy until energy and food resources are secured; this is because it is the lack of energy and food security that is driving the deforestation in forests of North Korea. Areas with a slope >20 • (the second grade), have a high risk of soil erosion, as such, forest restoration through the Afforestation and Reforestation Clean Development Mechanism (A/R CDM) project is required [41]. The slope criteria for the third grade were set at 15 • , which is the lower limit of the 15 • to 30 • slope range; the standard as per the North Korean agroforestry policy. The rationale for setting the lower limit of 8 • for the fourth grade, was based on the research results of Yu et al. [40] and the Korean Agricultural Land Law. The last fifth grade area was an area with a slope < 8 • ; this area may be afforested when all energy and agricultural resource problems are resolved. As such, the restoration of forests in North Korea is a very complex problem that must consider the scientific aspects of forest ecosystems, and socio-economic issues such as energy and food resources [41]. The restoration of forests is extremely challenging for the North Korea to achieve independently; thus, forestry cooperation of the Republic of Korea and the international communities is necessary [9,29,30]. search results of Yu et al. [40] and the Korean Agricultural Land Law. The last fifth grade area was an area with a slope < 8°; this area may be afforested when all energy and agricultural resource problems are resolved. As such, the restoration of forests in North Korea is a very complex problem that must consider the scientific aspects of forest ecosystems, and socio-economic issues such as energy and food resources [41]. The restoration of forests is extremely challenging for the North Korea to achieve independently; thus, forestry cooperation of the Republic of Korea and the international communities is necessary [9,29,30].

Implication of Phenological Classification with Sentinel Satellite and Deep Learning Algorithm
To identify the deforestation area in North Korea and determine target sites for priority afforestation, the phenological characteristics of vegetation were used to reflect the deforestation characteristics in North Korea. The derived phenological characteristics of vegetation were reflected in five scenes that were learned, resulting in high-accuracy classification results. These results confirmed that phenological features according to vegetation cover may appear in time-series satellite imagery and can be learned by the classification algorithm. As such the Sentinel-2 satellite imagery, which can acquire sufficient number of images due to its short visit period, was suitable for this methodology. This imagery has a spatial resolution of 10 m, which indicates its ability to reflect the characteristics of the small, patch-like deforestation that occurs in North Korea. Those imagery enabled confirmation of the area and spatial distribution of deforested land, and the spatial derivation of effective afforestation target sites.
By using a deep learning algorithm and a stepwise classification methodology, we were able to overcome limitations posed by the lack of labeling data. If high-accuracy labeling data exists, the land cover of North Korea may be more effectively classified with deep learning; however, field investigations are difficult due to due to limited accessibility of these areas. As such, precise classification using deep learning and high-resolution satellite imagery was feasible. Due to the high resolution of Sentinel-2 satellite imagery, the images covering all of North Korea is very huge data; this allowed for very high classification accuracy as the deep learning algorithm is suitable for big data analysis. Thus, it was possible to map North Korean land cover with high accuracy. However, a major limitation is that unpaved roads, small farmland or green areas in cities are classified as agricultural land, forest, on used area.
This multi-series vegetation phenological classification methodology enables the mapping of land cover on an annual basis. As such, it is possible to effectively track changes in forests in the Korean Peninsula, including North Korea and regions with similar climatic zones; using this data, it is possible to estimate the carbon uptake of forests.

Conclusions
A thorough classification of land cover and deforested areas using phenological characteristics and the deep learning technique has significant implications for afforestation in North Korea. The land cover of North Korea was classified with an accuracy of 84.8%, producing results similar to those of previous studies. this study derived a total of 2.6 million ha of afforestation sites, which were graded to enable effective afforestation and forest restoration. The application of deep learning algorithms for big data analysis of shortperiod, high-resolution satellite imagery was very effective. Though U-Net algorithm demands labeled data for training, the histogram threshold supplementing the lack of labeling data, and the classification accuracy demonstrates the adequacy of this method. A major limitation is that objectivity is not secured; every data is remotely sensed and visually interpreted and identified. Using the classification methodology of this study, it was easy to classify land cover in temperate areas with clear seasonal changes and high vegetative cover. As the phenological characteristics of vegetation appear on a yearly basis, land cover maps may be updated annually, enabling the effective tracking of land cover change.
As deforestation has been ongoing for many years in North Korea, it is necessary to track changes in land cover and forested areas. Using the methodology present in this study, land cover maps may be established on an annual basis, enabling the tracking of changes. If afforestation is carried out in stages based on the priority results from this study, the effect of afforestation is maximized and the potentiation for re-deforestation of afforested sites is low. If the forest area is expanded and the forest ecosystem is restored through successful afforestation, this may improve the overall ecosystem services in North Korea.