Next Article in Journal
Study of the Impact of Rural Land Transfer on the Status of Women in Rural Households
Previous Article in Journal
Research on the Coupling Co-ordination between Quality of County-Level New Urbanization and Ecosystem Service Value in Shaanxi Province
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Establishing an Ecological Network to Enhance Forest Connectivity in South Korea’s Demilitarized Zone

Major of Bio & Environmental Technology, Seoul Women’s University, 621 Hwarang-ro, Nowon-gu, Seoul 01797, Republic of Korea
Land 2024, 13(1), 106; https://doi.org/10.3390/land13010106
Submission received: 18 December 2023 / Revised: 12 January 2024 / Accepted: 17 January 2024 / Published: 18 January 2024
(This article belongs to the Section Landscape Ecology)

Abstract

:
This investigation delineated an ecological network within South Korea’s Demilitarized Zone (DMZ) to enhance ecosystem functionality, with a focus on forest landscape connectivity. The study employed spatial autocorrelation analysis and the minimum cumulative resistance model to identify key source patches. These patches are vital for maintaining connectivity across various topographies and land uses in the DMZ. Morphological spatial pattern analysis in the DMZ’s forest area showed a variation in forest patch types. The core class, less affected by external influences, was more prevalent in the eastern part. In contrast, the edge class, influenced by different environmental boundaries, was dominant in the western part. A 452.36 km long ecological network was constructed, serving as an essential corridor linking habitats of endangered species. This network covers a total forest area of 730.17 km2, which includes an existing forest protection area of 165.84 km2 (22.7%). The study also identified 564.33 km2 (77.3%) of previously non-designated forest regions as potential conservation areas. This aims to significantly increase forest connectivity within the DMZ. The research highlights the effectiveness of spatial planning tools in promoting ecosystem integrity in politically sensitive and restricted-access areas.

1. Introduction

In recent decades, the development-centric spatial planning paradigm has inflicted severe damage upon forest ecosystems across the globe. This harm primarily stems from habitat fragmentation, where a single habitat gets divided into numerous isolated segments [1,2]. Habitat fragmentation alters the structure and functioning of ecosystems, causing changes in environmental conditions such as light and wind, as well as disruptions in biological and non-biological flow dynamics across habitat patches [3]. Extensive research indicates that these intricate alterations ultimately lead to a decrease in biodiversity and ecological productivity. Persistent fragmentation has led to a situation where, as of 2015, around 20% of global forests are located within 100 m of a forest edge, and approximately 70% are within 1 km, posing a significant threat to the biodiversity of the world’s forest ecosystems [1].
Ongoing research continues to uncover the adverse effects of habitat fragmentation on forest ecosystem structure and function [2,4,5,6]. Moreover, recognizing that most social and ecological processes impacting trees and forests occur within specific spatiotemporal landscape areas [7], forest management has shifted towards a landscape-focused approach [8]. This shift is evident in national and supranational efforts aimed at managing forests within the context of the larger landscape. For example, the United States and Europe have established national parks and the Natura 2000 network, respectively, as areas where forest ecosystems are legally protected and conserved. South Korea, for its part, designates 20 classes of protected areas, covering 17.2% of the country’s total terrestrial area [9], in line with the United Nations Environment Programme’s Aichi Target 11 [10]. This target recommends conserving at least 17% of terrestrial and inland water areas and 10% of coastal and marine areas, especially those of significant biodiversity importance. However, as South Korea is a signatory to the Aichi biodiversity convention, it must also define additional protected areas to meet Post-2020 Global Diversity Framework Target 3, which sets a goal of conserving 30% of global land and sea areas by 2030 [11]. To achieve this, well-connected systems of protected areas must be designed, necessitating consideration of landscape connectivity. Landscape connectivity measures interactions that affect species migration, survival, gene flow, and the structure of the landscape, all of which are vital ecological processes [12]. Protecting landscape connectivity is, thus, a fundamental goal of ecosystem conservation [13,14,15].
A landscape-based ecological network consists of ecologically significant landscape segments that can contribute to biodiversity conservation [16]. Specifically, the forest patches play a core role in an ecological network as the structural connectivity of the forest patch is associated with its functional connectivity, such as genetic diversity, brought about by the potential migration of specific species [17,18]. Hence, when the structural connectivity of landscape elements is low or non-existent, the consequent deterioration of cause–effect relationships can damage ecological processes [19]. Considering this, strategies that increase the connectivity between forest patches with high and low landscape connectivity to enhance the ecosystem quality throughout the target area are expected to be useful approaches for managing forest landscapes [20,21].
Various theoretical approaches, including graph theory [22,23], source–sink theory [24,25], circuit theory [25,26,27,28,29], and habitat suitability theory [30], have been employed to establish ecological networks based on landscape connectivity. Notably, a recent approach combining morphological pattern theory [15] and graph theory [31] identifies core landscape elements influencing ecosystem connectivity [32,33]. This approach determines the minimum-cost path required to connect these elements, creating well-defined ecological networks that offer greater interpretive clarity than other methods [34,35,36,37,38]. Its practical applicability is particularly high for decision-making and planning.
This study aimed to establish an ecological network that enhances ecosystem function by promoting forest landscape connectivity. In particular, unlike previous studies that conducted global connectivity analyses to derive ecological network pathways [34,35,36,37,38], this study aimed to derive ecological network pathways that are effective in securing ecosystem connectivity in areas that are narrow and long, with distinct topographic and land use gradients. The study area encompassed the border region between South Korea and North Korea, where the absence of civilian activity after the Korean War allowed the forest ecosystem to recover towards its natural state [39,40,41,42]. The demilitarized zone (DMZ) in North Korea was excluded due to a lack of spatial data for analyzing forest landscape connectivity. Additionally, the DMZ’s barbed-wire border restricts the movement of large terrestrial mammals. Our study identified critical forest source patches for preserving forest landscape connectivity, and we developed an effective ecological network by incorporating topographic and land use characteristics of the study area.

2. Materials and Methods

The analytical processes employed to establish the ecological network for forest landscape connectivity are summarized in Figure 1. Each analytical step is described in detail in the following sub-sections.

2.1. Study Area

The study area encompasses a portion of the border region between South Korea and North Korea, located within the geographical coordinates of 126°38′48.69″ E–128°24′42.6″ E and 37°49′40.31″ N–38°38′8.53″ N. This South Korean border region includes the Demilitarized Zone (DMZ), which is defined as the area situated 2 km south of the Military Demarcation Line (MDL), as per the Korean Armistice Agreement, effective since 27 July 1953. Additionally, the DMZ extends to within 10 km of the Southern Limit Line (SLL) of the DMZ, forming the Civilian Control Zone (CCZ). The study area comprises portions of two administrative districts, Gyeonggi-do and Gangwon-do, one city (Paju-si), and six counties (Yeoncheon-gun, Cheorwon-gun, Hwacheon-gun, Yanggu-gun, Inje-gun, and Goseong-gun). For geographical accuracy, the study area was divided into four distinct regions: Paju (PJ), Yeoncheon (YC), Cheorwon (CW), and Gangwon (GW), which includes four counties: Hwacheon-gun, Yanggu-gun, Inje-gun, and Goseong-gun (Figure 2).
The study area exhibits an elevation range from 0 to 1216 m, characterized by an elevation pattern that inclines eastward and descends westward. The annual temperature within the region averages around 10.4 °C, with an average annual precipitation of approximately 1344.8 mm. Land use data for the study area are summarized in Table 1 and classified as follows: forest area was predominant at 70.96%, followed by agricultural area (9.08%) and grassland (6.64%). This distribution correlates with the topographical features, with higher elevations in the east and lower elevations in the west, resulting in a higher percentage of forested land in the eastern regions and a greater prevalence of land allocated for human activities, such as agriculture and urban development, in the western regions (Figure 3).
Deforestation attributed to military activities and forest fires has led to the wide distribution of vegetation, including grass, during the initial stages of secondary succession. Consequently, the ecological landscape is characterized by a simple structure primarily dominated by oak trees. In contrast, parts of the high-elevation eastern regions, where human activity is limited, feature old-growth forests [42].
The presence of access restrictions and landmines in the border area has impeded in-depth investigations of its flora and fauna. Nonetheless, both the Korea Ministry of Environment and the National Institute of Ecology [43] have compiled published research on the area, reporting a total of 6373 species inhabiting the border region. Among these, 102 species are classified as endangered by the Korea Ministry of Environment, comprising 18 Class I species (e.g., Urus thibetanus ussuricus, Naemorhedus caudatus, Moschus moschiferus, Grus japonensis, and Hyla suweonensis) and 84 Class II species (e.g., Pteromys volans aluco, Aegypius monachus, and Grus vipio).

2.2. Identification of Ecological Sources

2.2.1. Morphological Spatial Pattern Analysis (MSPA)

To analyze the landscape morphology of the study area, a binary raster map at a 30 m spatial resolution classifying land use as foreground (forest) and background (non-forest) was generated using the 1:5000 m scale digital land vector cover map produced by the Korea Ministry of Environment in 2020 [44]. Map features were then reclassified into seven morphological categories: core (core habitat patch), islet (isolated small habitat patch), perforation (gap within a habitat patch), edge (boundary between different land covers), loop (ring-like habitat), bridge (connection between habitat patches), and branch (narrow habitat) [15] via the eight-neighborhood analysis method in the GuidosToolbox [45]. Lastly, core patches deemed critical for protecting landscape connectivity were extracted as ecological sources to analyze connectivity.

2.2.2. Evaluation of Landscape Connectivity

Regional patches critical for connectivity were extracted to prevent bias in selecting candidates for ecological sources and to adequately reflect the forest distribution with its higher density toward the east of the study area. The landscape connectivity for each core patch was calculated based on the delta of the probability of connectivity (dPC) index using the Conefor 2.6 software [46] for core patches extracted in the MSPA [31,36,37,47]. The following dPC equation was used [48]:
P C = i = 1 n j = 1 , i j n P i j * a i a j A 2 L ,   d P C = P C P C P C × 100 %
Here, ai and aj are the area of the i and j patches, respectively; AL is the total area of the study region consisting of forest and non-forest areas; and Pij is the link intensity indicating the direct probability of dispersal across patches i and j. The threshold value of the patch connectivity distance was established at 1000 m, and the connectivity probability was fixed at 0.5, as referenced in previous studies [36,49]. Regional patches critical for connectivity were extracted to prevent bias in selecting candidates for ecological sources and to adequately reflect the forest distribution with its higher density toward the east of the study area. The landscape connectivity for each core patch was calculated based on the delta of the probability of connectivity (dPC) index using the Conefor 2.6 software [46] for core patches extracted in the MSPA [31,36,37,47].

2.3. Constructing the Ecological Network

2.3.1. Area Clustering Using Spatial Autocorrelation Analysis

To construct an ecological network that can enhance forest landscape connectivity throughout the study area, a spatial autocorrelation analysis based on the landscape connectivity index was used to cluster forest areas and draw corridors to indicate paths connecting all clusters. Whether attribute values of a spatial variable at a specific location are significantly correlated with a nearby attribute could be verified with a spatial autocorrelation analysis [50]. The study area was divided into 500 m × 500 m meshes, and the mean connectivity index for the core forest patch in each mesh was estimated, with the aim of ensuring the structural stability of core forest areas unaffected by edge effects and securing the feasibility of spatial planning. In addition, the spatial aggregation characteristics of these meshes were analyzed by drawing a local indicator of spatial association (LISA) distribution map using the GeoDa software 1.20 [37,47,51,52]. LISA is an analytical technique employed to study localized spatial aggregation phenomena, whereby clusters of high attribute values (HH) and those of low attribute values (LL) can be identified for a set level of significance (p = 0.05). The HH outliers surrounded by LL clusters and LL outliers surrounded by HH clusters can also be identified [37,53]. Among the detected core patches in each HH and LL cluster identified from the analysis, those with the highest landscape connectivity index were selected as ecological source patches.

2.3.2. Constructing Ecological Cost and Ecological Resistance Surfaces

Raster-type ecological and resistance surfaces at a 30 m resolution were produced using methods suggested by Wanghe et al. [34]. The cost value of each patch was calculated according to the product of its weight and suitability score:
c o s t = 100 w e i g h t × s u i t a b i l i t y   s c o r e
The efficiency of ecological network construction was improved by assigning high-weight values and suitability scores to land use classes with high ecosystem naturalness and by estimating the ecological cost per class (Table 2) [34]. For this estimation, the resistance, weight, and suitability coefficients utilized in this study were scientifically estimated based on a comprehensive review of the literature [34,54,55] and expert judgment in National Institute of Ecology, reflecting the relative importance and ecological impact of different land use classes on forest connectivity. While previous studies have focused on creating networks of urban green infrastructure, this study assigns resistance values to a network of forest patches.

2.3.3. Establishing the Ecological Network and Selecting Forest Area Candidates for Protection

Based on the constructed resistance surface, the minimum cumulative resistance (MCR) model was used to draw efficient corridors between core source patches. The MCR model comprised the following formula [36,56,57]:
m c r = f m i n j = n i = m ( D i j × R i )
Here, mcr is the resistance of the potential corridor connecting two patches; f is an unknown positive function that represents the minimum resistance at a specific spatial location, the distance to all sources, and the positive correlation across the landscape-based surfaces; Dij is the distance from the source point j to the spatial unit i; and Ri is the resistance coefficient for the spatial unit i [36].
Among the candidate corridors extracted using the MCR model, the core corridor paths constituting the ecological network were selected using the gravity model tool of ArcGIS 10.7 [58]. The gravity model simulates interactions between nodes connected through the corridors: more interactions indicate a higher significance of corridor connectivity. The equation outlining the gravity model is as follows [34]:
G a b = N a N b D 2 a b = L 2 m a x × S a × S b L 2 a b × p a × P b
Here, Gab is the interaction between nodes a and b; N is the weight on the corresponding node; Dab is the standardized resistance value of the potential corridor for nodes a and b; and S and P indicate the area and resistance of the corresponding node, respectively. The resistance value of a patch is defined according to its land use class, which is associated with the level of difficulty in constructing the given ecological corridor. Lab is the cumulative resistance of the corridor for nodes a and b; Lmax is the maximum resistance across all corridors in the study area.
Lastly, from potential corridors connecting core source areas at each cluster derived from the MCR model, corridors with the highest interaction values were selected to construct an ecological network. In addition, forests in areas where a network passed but was yet to be designated as a forest conservation zone were selected as protection-area candidates.

3. Results

3.1. Ecological Sources

Figure 4 shows the distribution of each landscape class following reclassification in the MSPA. The ratio of core classes throughout the study area accounted for the widest area at approximately 62%. The core area ratio in each region increased toward the east, while that of the edge area increased toward the west (Table 3).
A total of 2058 core patches were identified in the study area, with a mean area per patch of 0.71 ± 8.93 km2. The number of core areas distributed in each region was the greatest in YC (n = 692), followed by PJ (n = 519), GW (n = 437), and CW (n = 400), with the mean area as follows: GW (1.92 ± 16.91 km2), CW (0.72 ± 7.30 km2), YC (0.14 ± 1.18 km2), and PJ (0.05 ± 0.19 km2), indicating an increase in core area toward the east (Table 3).
The distribution of dPC estimated for the core patches in each region is presented in Figure 5a. Compared to the connectivity estimated for the entire study area (Figure 5b), the patches critical for protecting connectivity throughout the study area were distributed evenly in each region. An analysis of the entire study area revealed a trend of a concentrated distribution of most high-connectivity core patches toward the east, where there is a high forest density. The mean dPC of the core areas of the entire study area was 0.45, while the highest connectivity for the core patches was 78.42 ± 2.16 standard deviation. The CW region had the highest mean connectivity index (0.54 ± 4.41), followed by GW (0.49 ± 4.30), PJ (0.47 ± 1.89), and YC (0.34 ± 3.02) (Table 4).

3.2. Clustering Area by Landscape Connectivity

The spatial autocorrelation analysis of forest landscape connectivity reveals a distinct distribution pattern of high–high (HH) and low–low (LL) connectivity clusters across the study area. In the PJ and YC regions, there is a notable gradient in landscape connectivity (Figure 6), with lower connectivity observed on the western side and higher connectivity on the eastern side. The CW region exhibits a contrasting pattern, with lower connectivity on the southern side and higher connectivity on the northern side. The GW region presents a more complex scenario where HH and LL clusters are interspersed, indicating a mixed-connectivity landscape.

3.3. Ecological Network and Proposed Protected Area

A comprehensive ecological network stretching 452.36 km was established, designed to interlink core ecological source points across various regions (Figure 7). This network meticulously weaves through the border area, creating vital ecological corridors for several endangered species native to these habitats. The network specifically supports species such as H. suweonensis, G. japonensis and G. vipio in the PJ region; it continues through the YC region sustaining populations of G. japonensis and G. vipio; it supports M. moschiferus in the CW region; and it provides critical habitat connectivity for U. thibetanus ussuricus and N. caudatus in the GW region (Table 5).
The forests along the developed network are presented in Figure 8. The total area of these forests was 730.17 km2, with the areas of protected forest covering 165.84 km2. Based on this, forests not yet designated as protected areas (564.33 km2) were selected as potential candidates.

3.4. Comparing the Results of Constructing an Ecological Network with Global Approach

In constructing an ecological network across the study area, the sectional (stepwise) approach (Figure 8) demonstrated superior efficacy over the global approach (Figure 9). The topography-driven forest distribution, with high connectivity forests in the east and low connectivity in the west, necessitated a tailored strategy. Our sectional method, considering regional land-use characteristics, yielded a more effective ecological network, enhancing connectivity where it is most needed. This approach’s alignment with local spatial management ensures its practical applicability and potential for successful implementation in South Korea’s forest conservation efforts.

4. Discussion

4.1. Constructing Effective Ecological Networks That Take into Account Terrain and Land Use Characteristics

The establishment of ecological networks plays a crucial role in facilitating the colonization of wide-ranging species, primarily through the enhancement of species mobility and genetic exchange. These networks also have broader ecological implications: they can reduce habitat fragmentation, augment biodiversity, support pollinator migration, and enable species to adapt to climate change [16,17,18]. When targeting a particular species, it becomes imperative to incorporate coefficients that accurately reflect its behavioral traits, including home range and seasonal movement patterns, to compute the least cost distance. Integrating these extensive ecological benefits significantly augments the methodology of the current study, rendering it highly effective for the development of species conservation strategies that improve habitat connectivity. However, in this study, due to the absence of specified target species and the focus on establishing ecological network pathways to enhance forest connectivity, coefficients from prior studies were adopted, especially pertinent in regions with significant land-use intensity gradients, as observed in the study area.
Owing to the topographic characteristics of the Korean peninsula, where the east presents a high altitude and steep slopes and the west a low altitude with milder slopes, most forests are found concentrated on the eastern side, while human land-use areas, such as developed and agricultural areas, tend to be concentrated in the west. Hence, in constructing a longitudinal forest landscape network, there is the potential for a broad approach to producing a network biased towards the forests in the east. Even within the eastern area, the forest patches with the highest connectivity were concentrated on the southern side, indicating that a global approach to constructing an ecological network would ultimately partially limit the scope to those concentrated areas. The approach is, thus, limited in enhancing the ecosystem quality throughout the forest landscapes by constructing an ecological network for the target area of spatial planning. The approach is, thus, limited in enhancing the ecosystem quality throughout the forest landscapes by constructing an ecological network for the target area of spatial planning. Informed by recent advancements in spatial planning [59,60], this study adopted a stepwise approach over the traditional global methodology [33,34,35,36,37,38], where paths of the ecological network were individually tailored for each region based on specific land use characteristics and topographical challenges of the Korean peninsula, culminating in a final network that connects these paths, proving more effective in achieving our conservation goals. Our ecological network approach connected the forest landscapes more effectively than connections produced through a global approach. Because spatial management planning for natural environments in South Korea is conducted by each local government, the developed network is anticipated to prove effective in practice owing to its regional focus [61]. However, for an enhanced applicability of the method in policymaking, it is imperative that research on optimization persists to confront the inherent challenge associated with the stepwise methodology. This approach occasionally overlooks forests possessing substantial absolute connectivity but comparatively lower regional connectivity significance. This issue arises in the process of identifying core connectivity areas within highly connected regions. A case in point is the southwestern forests of the GW region, as depicted in Figure 6, which, despite their high connectivity, may be overlooked due to their relative positioning in the connectivity hierarchy. Addressing this limitation is crucial for refining the selection of key areas that contribute to overall ecological network efficacy.

4.2. Value of Ecological Networking in the South Korean Border Area

The ecological axis from the south to the north of the Korean peninsula, centered on a series of high-altitude mountain ranges on the eastern side, allows for relatively distinct compartmentalization, thus maintaining ecological connectivity to a certain degree due to limited land use by humans. However, connectivity is difficult to maintain along the ecological axis from the east to the west due to the land use by humans focused in the low-altitude regions. In this light, ecosystems at the border area with North Korea serve as a core axis of biodiversity in the Korean peninsula due to their limited land use by humans and uninterrupted connectivity due to the lack of roads, railroads, and urbanization [42]. However, as urban sprawl increases due to the overpopulation of metropolitan cities in South Korea, the pressure of human-related development is increasing, even near the border area [62]. In addition, the installation of structures for military purposes and road blockades used in the cold war landscape [63] could isolate forest ecosystems. The ecological connectivity-based management of border-area forest ecosystems is, thus, essential for conserving biodiversity in the Korean peninsula. Furthermore, considering regional characteristics with a high frequency of forest fires for military purposes, followed by the recovery of the forest landscape, the study area is anticipated to serve as a valuable reference model area for spatial planning, aiming to restore and manage forest landscapes in fire-prone regions in future [64]. Because the border area is still subject to military activity between South Korea and North Korea, the risk of remnants of war, such as landmines, limits the ability to thoroughly plan for managing forest ecosystems based on field surveys. To resolve such limitations, the current distribution of forests was determined by analyzing KOMPSAT-3 satellite images, aerial orthophotos, and the use of a national land-cover map [44] with classifications based on spatial thematic maps constructed by government agencies, such as the national forest and digital topographic maps. The paths of a forest ecological network that ensure landscape connectivity were drawn using the active and recently developed GIS-based spatial analysis. The methods applied in this study may prove effective in spatial planning for landscape ecosystem conservation in areas where field surveys are limited due to the large scale of the area or restricted access. It should also be noted that the structural integrity of forest patches determined in field surveys should be considered when deriving core forest patches to further enhance the effects of the constructed ecological network [65,66]. In addition, various ongoing studies are estimating the forest conservation value using remote sensing techniques [67,68,69], and the findings should be considered for constructing effective ecological networks, even in areas with a limited scope for field surveys.

4.3. Establishing Protected Areas for Ecological Connectivity in the South Korean Border Area

With the goal of forest ecosystem conservation, the South Korean government has designated legally protected areas for land management, including natural parks, ecological and scenery conservation areas, forest conservation zones, Baekdu-daegan protection areas, and protection districts for wildlife [70]. The total protected area, excluding duplicated areas where there are multiple designated protected areas, is around 10,160 km2 [9], accounting for approximately 10% of the total terrestrial area of South Korea (100,341.8 km2 as of 2021) [71]. Except for the Baekdu-daegan protection areas along a series of mountain ranges on the eastern side of the Korean peninsula, most protected forest ecosystem areas are separated by independent mountains. Ultimately, such natural compartmentalization is effective for preventing habitat loss and fragmentation within each protected area, despite the possibility of isolation due to changes in land use outside the protected area [72], making it difficult for the benefits resulting from forest protection to spread to neighboring, possibly unprotected, forests. Notably, land use by humans, mainly in low-altitude regions, may limit natural ecological processes, such as the migration and dispersal of wildlife that inhabit adequately conserved forests [73], and designating protected areas without considering ecosystem connectivity will not lead to the enhanced quality of overall forest ecosystems in each region.
Furthermore, while the Baekdu-daegan protection areas may be considered useful for adequately managing the functional connectivity of forest ecosystems in the north and south, besides differences in climate due to latitude, these areas exhibit mostly similar topographic characteristics. Hence, the diversity of forest landscape connectivity cannot be guaranteed regarding the topographic characteristics on a macroscopic level on the Korean peninsula. On the contrary, to ensure the effective conservation of national forests, it is necessary to increase connectivity across various forest and water edge ecosystems that develop along the macroscopic topography gradient from the east to the west. However, this has not been considered when designating most protected areas in South Korea. Thus, an effective method for designating protected areas that can connect various forest landscape classes in the east-to-west direction is expected to be useful for promoting forest management. The Korean Ministry of Environment in South Korea has already constructed a GIS-based digital thematic map that presents the level of conservation of the national ecosystems, such as the ecological and natural maps and the urban ecological map. These maps could be used to select forests with high conservation values as protected areas. Conversely, it is difficult to ensure that protected forest ecosystem areas selected according to the conservation level will reflect ecosystem connectivity as the criterion for assessing the conservation level as indicated on the maps, especially when it comes to assessing the conservation of vegetation and the structural integrity of plant communities. The vegetation metrics consider plant distribution and rarity, plant nativity, presence or absence of potential natural vegetation, the vertical structural integrity of the plant community, the time required for the community to develop, and inhabitation by plant species critical to conservation ecology. Thus, alongside valuable data for defining the structural conservation value of forest ecosystems, ecological connectivity will be useful for designating protected areas to best reflect both ecosystem structure and function. In this respect, the findings of our study on developing a method to establish a forest ecological network for the South Korean border area are anticipated to contribute to spatial planning for designating and managing protected areas in South Korea in the future.

5. Conclusions

South Korea has a considerably high population density (25th on the global scale and 516/km2 as of 2020) [71], with mountains constituting approximately 70% of the terrestrial land [74], and a land mass is characterized by differential land use characteristics to the east and west. The border area is the only terrestrial ecological axis in the east-to-west direction that has not been significantly impacted by rapid urbanization, making it a critical region for biodiversity conservation. It is necessary to develop conservation measures against potential forest landscape damage due to the increasing developmental pressure on this region. This study was conducted to establish an ecological network aiming to enhance forest landscape connectivity in the border area in South Korea, and the following conclusions were drawn.
(1)
Considering the macroscopic topographic characteristics of the Korean peninsula, such as a high altitude in the east and a low altitude in the west, with consequently varied land use patterns, the study area was divided into four regions. We selected the core source forest patches for each region by analyzing the forest landscape connectivity in each region. The paths drawn for constructing an ecological network were connected to reflect all paths throughout the study area. This method allowed for more effective connections throughout the area than those derived from a global analysis; hence, our method is anticipated to enhance the ecological connectivity of the entire study area.
(2)
A forest network spanning 452.36 km2 to connect the core source forest patches in each region effectively connects the habitats of the main endangered wildlife species inhabiting the study area. The developed network is, thus, anticipated to serve as an ecological corridor. Furthermore, to enhance forest ecosystem connectivity in the border area, it is necessary to consider forests of 564.33 km2 in the vicinity of the network paths as candidate forest conservation zones, which would contribute to achieving Target 3 outlined in the Post-2020 Global Diversity Framework.
(3)
The method applied in this study will prove effective for constructing ecological networks of two areas with unique land use and topographic characteristics while remaining versatile owing to its simple, intuitive, and quantitative analytical process.
(4)
Because this study utilized limited spatial information, only ecological functions were considered in the selection of forest reserve sites. However, it is believed that a more effective selection of reserve sites can be achieved by considering spatial information that considers the structural value of forest patches. In addition, although this study focused on establishing a method for optimal network construction, it is necessary to calculate the optimal resistance value according to the regional characteristics and the purpose of network construction when utilized for designing policy.

Funding

This work was supported by a research grant from Seoul Women’s University (2023-0130). The funder provided support in the form of salaries for author L.C.H. but did not have any additional role in the study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available due to national security policeis for the DMZ region.

Acknowledgments

I thank Yong Chan Cho of the Korea National Arboretum for reviewing the manuscript.

Conflicts of Interest

The author declares no competing interests.

References

  1. Haddad, N.M.; Brudvig, L.A.; Clobert, J.; Davies, K.F.; Gonzalez, A.; Holt, R.D.; Lovejoy, T.E.; Sexton, J.O.; Austin, M.P.; Collins, C.D.; et al. Habitat fragmentation and its lasting impact on Earth’s ecosystems. Sci. Adv. 2015, 1, e1500052. [Google Scholar] [CrossRef] [PubMed]
  2. Liu, J.; Wilson, M.; Hu, G.; Liu, J.; Wu, J.; Yu, M. How does habitat fragmentation affect the biodiversity and ecosystem functioning relationship? Landsc. Ecol. 2018, 33, 341–352. [Google Scholar] [CrossRef]
  3. Bennett, A.F.; Saunders, D.A. Habitat fragmentation and landscape change. In Conservation Biology for All; Sodhi, N.S., Ehrlich, P.R., Eds.; Oxford University Press: Oxford, UK, 2010. [Google Scholar] [CrossRef]
  4. Fahrig, L. Effects of habitat fragmentation on biodiversity. Annu. Rev. Ecol. Evol. Syst. 2003, 34, 487–515. [Google Scholar] [CrossRef]
  5. Chetcuti, J.; Kunin, W.E.; Bullock, J.M. Habitat fragmentation increases overall richness, but not of habitat-dependent species. Front. Ecol. Evol. 2020, 8, 607619. [Google Scholar]
  6. Zambrano, J.; Cordeiro, N.J.; Garzon-Lopez, C.; Yeager, L.; Fortunel, C.; Ndangalasi, H.J.; Beckman, N.G. Investigating the direct and indirect effects of forest fragmentation on plant functional diversity. PLoS ONE 2020, 15, e0235210. [Google Scholar] [CrossRef] [PubMed]
  7. Fischer, J.; Lindenmayer, D.B. Landscape modification and habitat fragmentation: A synthesis. Glob. Ecol. Biogeogr. 2007, 16, 265–280. [Google Scholar] [CrossRef]
  8. Wu, J. A landscape approach for sustainability science. In Sustainability Science: The Emerging Paradigm and the Urban Environment; Weinstein, M.P., Turner, R.E., Eds.; Springer: New York, NY, USA, 2012; pp. 59–77. [Google Scholar] [CrossRef]
  9. Korea Database of Protected Areas. Available online: http://kdpa.kr/ (accessed on 19 September 2022).
  10. UNEP. Aichi Biodiversity Targets. 2020. Available online: https://www.cbd.int/sp/targets/ (accessed on 8 December 2022).
  11. UNEP. 1st Draft of the Post-2020 Global Biodiversity Framework. UNEP—UN Environment Programme. 2021. Available online: http://www.unep.org/resources/publication/1st-draft-post-2020-global-biodiversity-framework (accessed on 8 December 2022).
  12. Luque, S.; Saura, S.; Fortin, M.-J. Landscape connectivity analysis for conservation: Insights from combining new methods with ecological and genetic data. Landsc. Ecol. 2012, 27, 153–157. [Google Scholar] [CrossRef]
  13. Forman, R.T.T.; Godron, M. Landscape Ecology; Wiley: New York, NY, USA, 1986; 648p. [Google Scholar]
  14. Crooks, K.R.; Sanjayan, M. Connectivity conservation: Maintaining connections for nature. In Connectivity Conservation; Crooks, K.R., Sanjayan, M., Eds.; Cambridge University Press: Cambridge, UK, 2006; pp. 1–20. [Google Scholar] [CrossRef]
  15. Vogt, P.; Ferrari, J.R.; Lookingbill, T.R.; Gardner, R.H.; Riitters, K.H.; Ostapowicz, K. Mapping functional connectivity. Ecol. Indic. 2009, 9, 64–71. [Google Scholar] [CrossRef]
  16. Fortin, M.J.; Dale, M.R.T.; Brimacombe, C. Network ecology in dynamic landscapes. Proc. Biol. Sci. 2021, 288, 20201889. [Google Scholar]
  17. Hellmund, P.C.; Smith, D. Designing Greenways: Sustainable Landscapes for Nature and People, 2nd ed.; Island Press: Washington, DC, USA, 2013; 288p. [Google Scholar]
  18. Velázquez, J.; Gutiérrez, J.; Hernando, A.; García-Abril, A. Evaluating landscape connectivity in fragmented habitats: Cantabrian capercaillie (Tetrao urogallus cantabricus) in northern Spain. For. Ecol. Manag. 2017, 389, 59–67. [Google Scholar]
  19. Hilty, J.; Worboys, G.L.; Keeley, A.; Woodley, S.; Lausche, B.J.; Locke, H.; Carr, M.; Pulsford, I.; Pittock, J.; White, J.W.; et al. Guidelines for Conserving Connectivity through Ecological Networks and Corridors; IUCN: Gland, Switzerland, 2020. [Google Scholar]
  20. Kettunen, M.; Terry, A.; Tucker, G.; Jones, A. Guidance on the Maintenance of Landscape Features of Major Importance for Wild Flora and Fauna—Guidance on the Implementation of Article 3 of the Birds Directive (79/409/EEC) and Article 10 of the Habitats Directive (92/43/EEC); Institute for European Environmental Policy (Institute for European Environmental Policy): Brussels, Belgium, 2007; 114p. [Google Scholar]
  21. García-Feced, C.; Saura, S.; Elena-Rosselló, R. Improving landscape connectivity in forest districts: A two-stage process for prioritizing agricultural patches for reforestation. For. Ecol. Manag. 2011, 261, 154–161. [Google Scholar] [CrossRef]
  22. Cantwell, M.D.; Forman, R.T.T. Landscape graphs: Ecological modeling with graph theory to detect configurations common to diverse landscapes. Landsc. Ecol. 1993, 8, 239–255. [Google Scholar] [CrossRef]
  23. Minor, E.S.; Urban, D.L. A graph-theory framework for evaluating landscape connectivity and conservation planning. Conserv. Biol. 2008, 22, 297–307. [Google Scholar] [CrossRef] [PubMed]
  24. Chen, L.; Fu, B.; Zhao, W. Source-sink landscape theory and its ecological significance. Front. Biol. China 2008, 3, 131–136. [Google Scholar] [CrossRef]
  25. Hansen, A. Contribution of source–sink theory to protected area science. In Sources, Sinks and Sustainability; Morzillo, A.T., Liu, J., Wiens, J.A., Hull, V., Eds.; Cambridge Studies in Landscape Ecology; Cambridge University Press: Cambridge, UK, 2011; 544p. [Google Scholar]
  26. McRae, B.H.; Dickson, B.G.; Keitt, T.H.; Shah, V.B. Using circuit theory to model connectivity in ecology, evolution, and conservation. Ecology 2008, 89, 2712–2724. [Google Scholar] [CrossRef] [PubMed]
  27. Thorne, J.H.; Choe, H.; Boynton, R.M.; Lee, D.K. Open space networks can guide urban renewal in a megacity. Environ. Res. Lett. 2020, 15, 094080. [Google Scholar] [CrossRef]
  28. Choe, H.; Keeley, A.T.H.; Cameron, D.R.; Gogol-Prokurat, M.; Hannah, L.; Roehrdanz, P.R.; Schloss, C.A.; Thorne, J.H. The influence of model frameworks in spatial planning of regional climate-adaptive connectivity for conservation planning. Landsc. Urban Plan. 2021, 214, 104169. [Google Scholar] [CrossRef]
  29. Landau, V.A.; Shah, V.B.; Anantharaman, R.; Hall, K.R. Omniscape.jl: Software to compute omnidirectional landscape connectivity. J. Open Source Softw. 2021, 6, 2829. [Google Scholar] [CrossRef]
  30. Boudreau, M.R.; Gantchoff, M.G.; Ramirez-Reyes, C.; Conlee, L.; Belant, J.L.; Iglay, R.B. Using habitat suitability and landscape connectivity in the spatial prioritization of public outreach and management during carnivore recolonization. J. Appl. Ecol. 2022, 59, 757–767. [Google Scholar] [CrossRef]
  31. Saura, S.; Torné, J. Conefor Sensinode 2.2: A software package for quantifying the importance of habitat patches for landscape connectivity. Environ. Model. Softw. 2009, 24, 135–139. [Google Scholar] [CrossRef]
  32. Saura, S.; Vogt, P.; Velázquez, J.; Hernando, A.; Tejera, R. Key structural forest connectors can be identified by combining landscape spatial pattern and network analyses. For. Ecol. Manag. 2011, 262, 150–160. [Google Scholar] [CrossRef]
  33. Wei, J.; Qian, J.; Tao, Y.; Hu, F.; Ou, W. Evaluating spatial priority of urban green infrastructure for urban sustainability in areas of rapid urbanization: A Case study of Pukou in China. Sustainability 2018, 10, 327. [Google Scholar] [CrossRef]
  34. Wanghe, K.; Guo, X.; Wang, M.; Zhuang, H.; Ahmad, S.; Khan, T.U.; Xiao, Y.; Luan, X.; Li, K. Gravity model toolbox: An automated and open-source ArcGIS tool to build and prioritize ecological corridors in urban landscapes. Glob. Ecol. Conserv. 2020, 22, e01012. [Google Scholar] [CrossRef]
  35. An, Y.; Liu, S.; Sun, Y.; Shi, F.; Beazley, R. Construction and optimization of an ecological network based on morphological spatial pattern analysis and circuit theory. Landsc. Ecol. 2021, 36, 2059–2076. [Google Scholar] [CrossRef]
  36. Li, Y.-Y.; Zhang, Y.-Z.; Jiang, Z.-Y.; Guo, C.-X.; Zhao, M.-Y.; Yang, Z.-G.; Guo, M.-Y.; Wu, B.-Y.; Chen, Q.-L. Integrating morphological spatial pattern analysis and the minimal cumulative resistance model to optimize urban ecological networks: A case study in Shenzhen City, China. Ecol. Process. 2021, 10, 63. [Google Scholar] [CrossRef]
  37. Wang, S.; Wu, M.; Hu, M.; Fan, C.; Wang, T.; Xia, B. Promoting landscape connectivity of highly urbanized area: An ecological network approach. Ecol. Indic. 2021, 125, 107487. [Google Scholar] [CrossRef]
  38. Jiang, J.; Abulizi, A.; Abliz, A.; Zayiti, A.; Akbar, A.; Ou, B. Construction of landscape ecological security pattern in the Zhundong region, Xinjiang, NW China. Int. J. Environ. Res. Public Health 2022, 19, 6301. [Google Scholar] [CrossRef]
  39. John, K.H. The Korean DMZ: A fragile ecosystem. Science 1998, 280, 803. [Google Scholar] [CrossRef]
  40. Kim, K.-G.; Cho, D.-G. Status and ecological resource value of the Republic of Korea’s de-militarized zone. Landsc. Ecol. Eng. 2005, 1, 3–15. [Google Scholar] [CrossRef]
  41. Choi, S.A.; Park, E.-J.; Park, S.-H. Conservation values of major resources in the Korean DMZ and its vicinity. Policy Res. 2010, 1, 1–216. [Google Scholar]
  42. Cho, D. The Ecological Values of the Korean Demilitarized Zone (DMZ) and International Natural Protected Areas. MUN HWA JAE–Annu. Rev. Cult. Herit. Stud. 2019, 52, 272–287. [Google Scholar]
  43. Korea Ministry of Environment; National Institute of Ecology. 5929 Species of Wildlife, Including 101 Endangered Species, Inhabit the DMZ; Korea Ministry of the Environment Press Release: Sejong-si, Republic of Korea, 2010.
  44. Korea Ministry of Environment. Construction of Current Land Use Map in 2020 (No. KME-11-1480000-001737-01); Ministry of the Environment: Sejong-si, Republic of Korea, 2020; 138p.
  45. GuidosToolbox Software. Available online: https://forest.jrc.ec.europa.eu/en/activities/lpa/gtb/ (accessed on 19 March 2022).
  46. Conefor 2.6 Software. Available online: http://conefor.org (accessed on 23 May 2022).
  47. Hu, C.; Wang, Z.; Wang, Y.; Sun, D.; Zhang, J. Combining MSPA-MCR model to evaluate the Ecological Network in Wuhan, China. Land 2022, 11, 213. [Google Scholar]
  48. Saura, S.; Rubio, L. A common currency for the different ways in which patches and links can contribute to habitat availability and connectivity in the landscape. Ecography 2010, 33, 523–537. [Google Scholar]
  49. Liu, W.; Hughes, A.C.; Bai, Y.; Li, Z.; Mei, C.; Ma, Y. Using landscape connectivity tools to identify conservation priorities in forested areas and potential restoration priorities in rubber plantation in Xishuangbanna, Southwest China. Landsc. Ecol. 2020, 35, 389–402. [Google Scholar]
  50. Anselin, L. Local indicators of spatial association—LISA. Geogr. Anal. 1995, 27, 93–115. [Google Scholar]
  51. GeoDa Software. Available online: http://geodacenter.github.io (accessed on 23 May 2022).
  52. Martinho, V.J.P.D. Forest fires across Portuguese municipalities: Zones of similar incidence, interactions and benchmarks. Environ. Ecol. Stat. 2018, 25, 405–428. [Google Scholar]
  53. Alvarado-Serrano, D.F.; Hickerson, M.J. Spatially explicit summary statistics for historical population genetic inference. Methods Ecol. Evol. 2016, 7, 418–427. [Google Scholar]
  54. Pelletier, D.; Clark, M.; Anderson, M.G.; Rayfield, B.; Wulder, M.A.; Cardille, J.A. Applying circuit theory for corridor expansion and management at regional scales: Tiling, pinch points, and omnidirectional connectivity. PLoS ONE 2014, 9, e84135. [Google Scholar]
  55. Choe, H.; Thorne, J. Omnidirectional connectivity of urban open spaces provides context for local government redevelopment plans. Landsc. Ecol. Eng. 2019, 15, 245–251. [Google Scholar] [CrossRef]
  56. Knaapen, J.P.; Scheffer, M.; Harms, B. Estimating habitat isolation in landscape planning. Landsc. Urban Plan. 1992, 23, 1–16. [Google Scholar] [CrossRef]
  57. Su, Y.; Chen, X.; Liao, J.; Zhang, H.; Wang, C.; Ye, Y.; Wang, Y. Modeling the optimal ecological security pattern for guiding the urban constructed land expansions. Urban For. Urban Green. 2016, 19, 35–46. [Google Scholar] [CrossRef]
  58. Gravity Model Tool of the ArcGIS 10.7 Software. Available online: https://github.com/wanghekunyuan/Gravity-model-toolbox (accessed on 25 May 2022).
  59. Jiang, H.; Peng, J.; Zhao, Y.; Xu, D.; Dong, J. Zoning for Ecosystem Restoration Based on Ecological Network in Mountainous Region. Ecol. Indic. 2022, 142, 109138. [Google Scholar] [CrossRef]
  60. Xiao, H.; Guo, Y.; Wang, Y.; Xu, Y.; Liu, D. Evaluation and Construction of Regional Ecological Network Based on Multi-Objective Optimization: A Perspective of Mountains–Rivers–Forests–Farmlands–Lakes–Grasslands Life Community Concept in China. Appl. Sci. 2022, 12, 9600. [Google Scholar] [CrossRef]
  61. Bae, M.G. A study on environmental conservation plan based on spatialization method in local governments. Environ. Policy 2017, 25, 25–60. [Google Scholar]
  62. Kim, J.H.; Park, S.; Kim, S.H.; Lee, E.J. Long-term land cover changes in the western part of the Korean demilitarized zone. Land 2021, 10, 708. [Google Scholar] [CrossRef]
  63. Jung, G.; Han, M.; Kang, I.; Jeon, W. Humanistic Research on the Current Non-peace State of the Border Region at DMZ; Korea Institute for National Unification: Seoul, Republic of Korea, 2020; pp. 35–37. [Google Scholar]
  64. Urgenson, L.S.; Nelson, C.R.; Haugo, R.D.; Halpern, C.B.; Bakker, J.D.; Ryan, C.M.; Waltz, A.E.M.; Belote, R.T.; Alvarado, E. Social perspectives on the use of reference conditions in restoration of fire-adapted forest landscapes. Restor. Ecol. 2018, 26, 987–996. [Google Scholar] [CrossRef]
  65. Harvey, E.; Gounand, I.; Ward, C.L.; Altermatt, F. Bridging ecology and conservation: From ecological networks to ecosystem function. J. Appl. Ecol. 2017, 54, 371–379. [Google Scholar] [CrossRef]
  66. Guimarães, P.R. The structure of ecological networks across levels of organization. Annu. Rev. Ecol. Evol. Syst. 2020, 51, 433–460. [Google Scholar] [CrossRef]
  67. Jin, Y.; Jeong, S.; Jeong, S.; Lee, D. Assessment on the forest conservation value considering forest ecosystem services—The case of Gapyung-gun. J. Environ. Impact Assess. 2015, 24, 420–431. [Google Scholar] [CrossRef]
  68. Mikusiński, G.; Orlikowska, E.H.; Bubnicki, J.W.; Jonsson, B.G.; Svensson, J. Strengthening the network of high conservation value forests in Boreal landscapes. Front. Ecol. Evol. 2021, 8, 595730. [Google Scholar] [CrossRef]
  69. Munteanu, C.; Senf, C.; Nita, M.D.; Sabatini, F.M.; Oeser, J.; Seidl, R.; Kuemmerle, T. Using historical spy satellite photographs and recent remote sensing data to identify high-conservation-value forests. Conserv. Biol. 2022, 36, e13820. [Google Scholar]
  70. Kim, G.H.; Kong, S.J.; Kim, O.S.; Son, S.W.; Lee, E.J. A strategy on extracting terrestrial protected areas of the Republic of Korea under the convention on biological diversity. J. Assoc. Korean Geogr. 2017, 6, 407–423. [Google Scholar]
  71. National Index of South Korea. Available online: http://www.index.go.kr/potal/ (accessed on 19 September 2022).
  72. Hong, J.-P.; Shim, Y.-J.; Heo, H.-Y. A study on Aichi biodiversity target 11—Focused on quantitative expansion goals and qualitative improvement goals of protected areas. J. Korean Soc. Environ. Restor. Technol. 2017, 20, 43–58. [Google Scholar]
  73. Niebuhr, B.B.S.; Wosniack, M.E.; Santos, M.C.; Raposo, E.P.; Viswanathan, G.M.; da Luz, M.G.E.; Pie, M.R. Survival in patchy landscapes: The interplay between dispersal, habitat loss and fragmentation. Sci. Rep. 2015, 5, 11898. [Google Scholar]
  74. Choi, J.S.; Jin, J.H.; Shim, W.J.; An, Y.S.; Shin, H.S.; Lee, S.J.; Park, S.J. A study on the development of topographical variables and algorithm for mountain classification. J. Korean Geomorphol. Assoc. 2018, 25, 1–18. [Google Scholar]
Figure 1. Flowchart outlining the study.
Figure 1. Flowchart outlining the study.
Land 13 00106 g001
Figure 2. Map of the study area. The distribution per land use class is not shown on the map, according to the national security policy. The vertical cross-section of the study area was obtained from a shuttle-radar-topography-mission digital elevation model of a 30 m resolution (30 m SRTM DEM).
Figure 2. Map of the study area. The distribution per land use class is not shown on the map, according to the national security policy. The vertical cross-section of the study area was obtained from a shuttle-radar-topography-mission digital elevation model of a 30 m resolution (30 m SRTM DEM).
Land 13 00106 g002
Figure 3. Ratio of land use class area by region (a) and elevation (b). Area for each land use class was estimated based on the 1:5000 land-cover map produced by the Korean Ministry of Environment in 2021, and elevation data were obtained from the 30 m SRTM DEM.
Figure 3. Ratio of land use class area by region (a) and elevation (b). Area for each land use class was estimated based on the 1:5000 land-cover map produced by the Korean Ministry of Environment in 2021, and elevation data were obtained from the 30 m SRTM DEM.
Land 13 00106 g003
Figure 4. Distribution map of each morphological spatial pattern analysis (MSPA) class. Points (AC) represent randomly selected areas with different forest landscape distribution characteristics.
Figure 4. Distribution map of each morphological spatial pattern analysis (MSPA) class. Points (AC) represent randomly selected areas with different forest landscape distribution characteristics.
Land 13 00106 g004
Figure 5. Comparison of a map that combines the connectivity indices (dPC) analyzed for (a) each region and a map that analyses connectivity for (b) the entire study area at once. (a) A regionally even distribution of high-connectivity forests across the study area, denoted by green areas. In contrast, the figure in (b) reveals that these high-connectivity forests are predominantly confined to the eastern regions, as indicated by the concentration of green.
Figure 5. Comparison of a map that combines the connectivity indices (dPC) analyzed for (a) each region and a map that analyses connectivity for (b) the entire study area at once. (a) A regionally even distribution of high-connectivity forests across the study area, denoted by green areas. In contrast, the figure in (b) reveals that these high-connectivity forests are predominantly confined to the eastern regions, as indicated by the concentration of green.
Land 13 00106 g005
Figure 6. The local autocorrelation distribution of the connectivity index (dPC) for the core patches.
Figure 6. The local autocorrelation distribution of the connectivity index (dPC) for the core patches.
Land 13 00106 g006
Figure 7. Paths comprising the ecological network developed based on forest landscape connectivity and distribution of the main endangered species. Photos of endangered species were obtained from the National Institute of Ecology (https://www.nie.re.kr/nieEng/main/main.do) (accessed on 23 September 2022).
Figure 7. Paths comprising the ecological network developed based on forest landscape connectivity and distribution of the main endangered species. Photos of endangered species were obtained from the National Institute of Ecology (https://www.nie.re.kr/nieEng/main/main.do) (accessed on 23 September 2022).
Land 13 00106 g007
Figure 8. Forest conservation network in the study area. The map shows the established ecological network (red line) connecting forests (green) and currently designated forest conservation zone (orange), with key ecological source points (stars). Forests not currently designated as a forest conservation zone have been identified as candidates for additional protected areas.
Figure 8. Forest conservation network in the study area. The map shows the established ecological network (red line) connecting forests (green) and currently designated forest conservation zone (orange), with key ecological source points (stars). Forests not currently designated as a forest conservation zone have been identified as candidates for additional protected areas.
Land 13 00106 g008
Figure 9. Ecological network produced through a global analysis of landscape connectivity. The black line indicates central tendency, dispersion, and directional trends of ecological connectivity based on the forest patch connectivity index (dPC) derived using the Directional Distribution tool of the ArcGIS 10.7 software.
Figure 9. Ecological network produced through a global analysis of landscape connectivity. The black line indicates central tendency, dispersion, and directional trends of ecological connectivity based on the forest patch connectivity index (dPC) derived using the Directional Distribution tool of the ArcGIS 10.7 software.
Land 13 00106 g009
Table 1. Land use classes in the study area (percentage of total area). Non-investigated indicates areas where investigation is prohibited for military purposes or due to the North Korean control of the territory. The units of area and percentage are km2 and %, respectively.
Table 1. Land use classes in the study area (percentage of total area). Non-investigated indicates areas where investigation is prohibited for military purposes or due to the North Korean control of the territory. The units of area and percentage are km2 and %, respectively.
Land Use ClassArea (Percentage)
PJYCCWGWTotal
Developed area4.66 (2.27)3.98 (1.18)7.47 (1.82)8.59 (0.81)24.69 (1.23)
Agricultural area44.19 (21.55)39.50 (11.76)73.93 (18)24.89 (2.35)182.50 (9.08)
Forest73.34 (35.77)163.27 (48.59)250.77 (61.07)938.76 (88.72)1426.14 (70.96)
Grassland22.24 (10.85)38.91 (11.58)47.20 (11.50)25.18 (2.38)133.54 (6.64)
Wetland5.93 (2.89)2.75 (0.82)9.82 (2.39)7.67 (0.72)26.17 (1.30)
Bare ground6.50 (3.17)3.54 (1.05)3.14 (0.76)11.13 (1.05)24.30 (1.21)
Water12.73 (6.21)3.70 (1.10)6.62 (1.61)4.99 (0.47)28.05 (1.40)
Non-investigated35.45 (17.29)80.34 (23.91)11.67 (2.84)36.92 (3.49)164.38 (8.18)
Total205.03 (100)336.00 (100)410.61 (100)1058.13 (100)2009.77 (100)
Table 2. Coefficients applied to estimate ecological cost and resistance.
Table 2. Coefficients applied to estimate ecological cost and resistance.
Land Use ClassWeightSuitabilityCostResistance
Developed area0.0510095100
Agricultural area0.25608530
Forest0.66100340.1
Grassland0.380760.5
Wetland0.2560851
Bare ground0.2509010
Water0.110991
Table 3. Characteristics of each MSPA class. The background area includes non-forested and non-investigated areas. The units of area and percentage are km2 and %, respectively.
Table 3. Characteristics of each MSPA class. The background area includes non-forested and non-investigated areas. The units of area and percentage are km2 and %, respectively.
MSPA ClassArea (Percentage)
PJYCCWGWTotal
Background131.69 (64.23)172.77 (51.42)159.98 (38.96)119.78 (11.32)584.23 (29.07)
Core45.39 (22.14)120.21 (35.78)217.84 (53.05)864.09 (81.66)1247.54 (62.07)
Islet0.45 (0.22)0.80 (0.24)0.30 (0.07)0.85 (0.08)2.40 (0.12)
Bridge0.42 (0.20)1.90 (0.57)4.76 (1.16)24.20 (2.29)31.28 (1.56)
Edge19.37 (9.44)28.41 (8.46)20.02 (4.88)27.61 (2.61)95.41 (4.75)
Loop1.10 (0.54)2.35 (0.70)3.45 (0.84)13.63 (1.29)20.54 (1.02)
Perforation2.90 (1.41)3.74 (1.11)1.67 (0.41)4.10 (0.39)12.41 (0.62)
Branch3.72 (1.81)5.81 (1.73)2.59 (0.63)3.85 (0.36)15.97 (0.79)
Total205.03 (100)336.00 (100)410.61 (100)1058.13 (100)2009.77 (100)
Table 4. Description of patch area and connectivity for each region (Aver.: Average value, Max.: Maximum value, Std.: Standard deviation).
Table 4. Description of patch area and connectivity for each region (Aver.: Average value, Max.: Maximum value, Std.: Standard deviation).
RegionCountArea (km2)Connectivity (dPC)
Aver.Max.Std.Aver.Max.Std.
PJ5290.052.040.190.4718.911.89
YC6920.1420.721.180.3453.163.02
CW4000.72140.127.300.5478.424.41
GW4371.92220.8516.910.4959.104.30
Total20580.71220.858.930.4578.422.16
Table 5. Classification and habitat characteristics of major endangered species in the study area. Appendix I lists species that are the most endangered among CITES-listed animals and plants, and Appendix II lists species that are not necessarily now threatened with extinction but that may become so unless trade is closely controlled.
Table 5. Classification and habitat characteristics of major endangered species in the study area. Appendix I lists species that are the most endangered among CITES-listed animals and plants, and Appendix II lists species that are not necessarily now threatened with extinction but that may become so unless trade is closely controlled.
Species Name (Scientific)GDEWIUCNCITESGDRLMajor Habitat TypeMajor Behavior
Moschus moschiferusClass IVUAppendix IICRForest area, Mountainous regionSeasonal movement; Altitudinal migration
Hyla suweonensisClass IEN-ENWetland, Aquatic vegetation areaBreeding behavior;
Egg laying on rice paddy
Grus japonensisClass IENAppendix IENWetlandFeeding behavior;
Foraging in rice paddy and field
Grus vipioClass IIVUAppendix INTWetlandFeeding behavior;Foraging in rice paddy and field
Ursus thibetanus ussuricusClass IVUAppendix IENForest areaTerritoriality;
Wide homeranges with seasonal variations
Naemorhedus caudatusClass IVUAppendix IVUForest area, Mounatinous terrainMovement;
Agile climbers, adapted for steep and rugged terrain
GDEW: Government-designated endangered wildlife, IUCN: IUCN Red List, GDRL: Government-designated Red List, CR: Critically endangered, EN: Endangered, VU: Vulnerable, NT: Near threatened.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lim, C.H. Establishing an Ecological Network to Enhance Forest Connectivity in South Korea’s Demilitarized Zone. Land 2024, 13, 106. https://doi.org/10.3390/land13010106

AMA Style

Lim CH. Establishing an Ecological Network to Enhance Forest Connectivity in South Korea’s Demilitarized Zone. Land. 2024; 13(1):106. https://doi.org/10.3390/land13010106

Chicago/Turabian Style

Lim, Chi Hong. 2024. "Establishing an Ecological Network to Enhance Forest Connectivity in South Korea’s Demilitarized Zone" Land 13, no. 1: 106. https://doi.org/10.3390/land13010106

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