Next Article in Journal
Are Aquavoltaics Investable? A Framework for Economic and Environmental Cost-Benefit Analysis
Next Article in Special Issue
Interaction of Filth Flies and Epigeal Arthropods with Soil Nitrogen and Gas Emissions in Grazing Systems under a Legacy of Low Fertilization
Previous Article in Journal
Influence of Climatic Factors on Yields of Pistachio, Mango, and Bananas in Iran
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fine-Scale Species Distribution Modeling of Abies koreana across a Subalpine Zone in South Korea for In Situ Species Conservation

National Institute of Ecology, 1210 Geumgang, Seocheon 33657, Republic of Korea
*
Author to whom correspondence should be addressed.
Sustainability 2023, 15(11), 8964; https://doi.org/10.3390/su15118964
Submission received: 8 May 2023 / Revised: 26 May 2023 / Accepted: 30 May 2023 / Published: 1 June 2023
(This article belongs to the Special Issue Biodiversity, Biologic Conservation and Ecological Sustainability)

Abstract

:
Severe declines in the population of Abies koreana, a conifer native to the subalpine regions of South Korea, are likely a consequence of climate change. However, local-scale modeling of the species’ spatial distribution has seen limited application to in situ conservation policies. Therefore, we intended for this study to examine the applicability of fine-scale species distribution modeling of A. koreana in the Mt. Jiri National Park area in S. Korea in order to explore candidate areas for its in situ conservation. We simulated the potential habitat of the species in the area with four separate models using different dominance patterns, then created an index based on habitability probabilities and residual durations to determine priority conservation areas. Under the highest sensitivity of potential habitats to temperature (spatially downscaled based on geomorphological characteristics), rapid habitat reduction occurred under climate warming in all experiments. At the regional scale, hydrological characteristics such as precipitation and slope characterized different secondary habitat distributional patterns among the experiments. Final conservation priority sites were identified based on specified criteria for the designed index. Our results suggest that a fine-scale modeling system with adequate spatial resolution of environmental inputs is advantageous in representing local habitat characteristics of A. koreana and can be applied to in situ conservation strategies.

1. Introduction

Globally, isolated ecosystems such as mountains and islands are expected to be affected most by climate change [1,2]. The adaptability of species within most ecosystems can be diagnosed based on the movement of their habitats in response to the increase in average temperatures. However, regions such as mountains and islands, regardless of the mobility of species, could lose such opportunities for adaptation due to their isolated habitats [3].
The subalpine ecosystem in South Korea is representative of such isolated ecosystems that are vulnerable to the effects of climate change [4,5,6]. The subalpine zone in South Korea is mainly distributed in mountainous areas such as Mts. Seorak, Jiri, and Halla [6,7]. Subsequently, this subalpine zone is located in regions that are at the highest elevation of the temperate climate zone. This zone mainly comprises specific conifer species such as the Korean needle fir (Abies holophylla) and Korean fir (Abies koreana) [8]. These species have a short growing season and are sensitive to climate factors, such as continuous snowfall, dry winds, extreme temperatures, soil moisture, and evapotranspiration stress [9]. In particular, the decline of A. koreana, a unique conifer species in Korea, has been observed prominently in Mts. Jiri and Halla because of climate change [10,11]. Recent investigations [7,10] have revealed a rapid increase in the decline of A. koreana in the Mt. Jiri area over the past five years.
The causes of coniferous forest decline in high-altitude and subalpine zones worldwide are ascribed to climate change and, in particular, the increasing occurrences of winter heatwaves and spring droughts [12]. However, opinions differ as to the cause of the decline of A. koreana in the subalpine zone of South Korea. For example, several studies have suggested that the relatively shallow root system of A. koreana makes it more sensitive to changes in soil moisture conditions, and the spring droughts caused by strong winds that dry out the soil make the species more susceptible to decline than the other species [10,13,14]. Other studies have proposed the excess water supply derived from increased precipitation as the cause of the decline [15,16]. These differing opinions highlight the possibility that the cause of decline in the subalpine zone could vary according to the regional characteristics of the climate environment, even within a relatively small territory. For instance, the representative regions of coniferous forest decline in Korea, namely Mts. Jiri and Halla, are located approximately 2° apart in latitude (Mt. Jiri 35.33° N, Mt. Halla 33.36° N). However, the climate characteristics of these regions differ markedly, with Mt. Jiri having a continental climate, whereas Mt. Halla has a maritime climate.
To mitigate the decline of subalpine coniferous species, authorities follow two main approaches, namely ex situ and in situ conservation. In ex situ conservation, the surviving members of the species are moved to new habitats with similar environmental conditions. In situ conservation implies identifying the last surviving areas within the current habitat that are affected the least by climate change and implementing protective measures to limit the influence of external factors. In both approaches, species distribution models are typically employed for predicting potential habitats within a given area.
Species distribution modeling techniques have been developed with various approaches such as statistical and machine learning types that define the relationship between observed species distribution and environmental features. Due to their own strength and weakness in their simulating performances and accuracy, ensemble simulation systems using multi-model combination are also popular approaches that intended to reduce technical uncertainties. However, in applying advanced modelling techniques, the suitability of input data into experimental designs such as spatial scale is also critical in affecting modelling accuracy. In Korea, studies on potential habitats mainly employ low-resolution data (spatial resolution ~1 km) for the entire country [17,18,19]. However, the results of these studies could vary depending on the spatial resolution of the climate data and the information on species distribution, particularly for A. koreana, which has an extremely limited habitat range. For instance, Koo et al. (2016) [18] simulated potential habitats for A. koreana across the country at a resolution of 1 km, using only 14 observations. The results of applying future climate change scenarios in their study showed an expansion of potential habitat for A. koreana. This finding could be ascribed to the limited amount of observation data and low-resolution climate data (1 km), which failed to properly reflect that the actual distribution showed a high correlation with elevation.
Some previous experiments utilized higher resolutions to predict potential habitats. For example, using a climate change factor at a resolution of 1 km and other topographic information at a resolution of 100 m, Park et al. (2015) [20] predicted an overall decrease in the potential habitat of A. koreana. Compared with the results of Koo et al. (2016) [18], this finding suggested that a difference in spatial resolution could affect the accuracy of the prediction results. Notably, Manzoor et al. (2018) [21] demonstrated through comparative experiments conducted with environmental climate information at various spatial resolutions that increasing the resolution could not always produce prediction accuracy. However, only the input data corresponding to topographic information were used at a high resolution in their study, and no consideration was given to whether there were adequate observation data for a given species at the particular resolution. Consequently, sufficient observation data for a given species are considered a requisite for reliable prediction results at higher resolutions.
Alternatively, studies employing low-resolution data for potential habitat predictions may be suitable for exploring ex situ alternative habitat candidate areas in response to the effects of climate change in a relatively large area. However, prediction studies for in situ conservation planning in a relatively narrow area require precise predictions. More precise terrain and climate information is required to predict microhabitats that respond to microclimate environmental factors. Accordingly, grid-type high-resolution climate and environmental data are required for such research areas. However, obtaining climate data at a resolution of more than 1 km is challenging. For example, current regional-scale weather or climate prediction models have a maximum resolution of approximately 12.5 km. Although statistical techniques are being developed and applied to downscale data to a higher resolution [22], the current downscaling techniques have a limitation—that uncertainty increases as the resolution increases, owing to a dependence on the availability of ground observations. As the currently available downscaling techniques mainly consider altitude, aspect, slope, and other factors, it is worth applying such techniques to species, such as A. koreana, with a clear distribution pattern according to altitude.
In this study, we aimed to examine the applicability of fine-scale distribution modeling to A. koreana across Mt. Jiri National Park in South Korea, using the newly available high-resolution climate data from National Institute of Ecology (NIE). Our goal was to identify candidate areas for the in situ species conservation of A. koreana to combat the impact of climate change. The NIE has recently released fine-scale grid-type (100 m spatial resolution; the highest available in Korea) climate data for this area [23]. Our modeling experiments included various simulations of potential habitats of A. koreana under the current environmental and climate status and provided predictions for future changes in its distribution. Consequently, based on the simulated outputs, we attempted to develop an index for extracting priority information on in situ conservation areas for A. koreana.

2. Materials and Methods

2.1. Study Area

The study area constituted Mt. Jiri National Park, which covers an area of ~483 km2. Its highest peak, Cheonwangbong, is located on the eastern side of the park at an elevation of 1915 m. The distribution of A. koreana clusters in the park accounts for approximately 10% of the total park area. We gathered data from a survey of coniferous forest distributions in the subalpine zone conducted by the Korea National Park Research Institute [7] to determine the distribution of A. koreana clusters within the park. These data were acquired as a polygon shapefile with four levels of A. koreana dominance (10–25%, 25–50%, 50–75%, and >75%), expressed as the proportion of A. koreana stems per total stem count [6]. In this region, A. koreana occurs at an elevation of 800 m or higher but is concentrated at elevations above 1400 m. The species has a strong presence along the ridgeline of the highest peaks in the western Banyabong area (~1600 m) and the eastern Yeongsinbong and Cheonwangbong areas (~1700 m) (Figure 1). More than 75% of the total distribution area of A. koreana was observed at 1500–1800 m, with a slight decrease in distribution at an elevation above 1800 m (Figure 2a). The increase in A. koreana distribution along with elevation was observed in the high-density dominance category of over 75%, suggesting that the elevation, or temperature profile varying with the elevation, would be critical for the selection of environmental inputs for the modeling. In addition, this further suggests that contributions of climate and environmental input variables to the modeling performances may vary depending on the distribution density of species, and thus separate experiments according to the distribution density may be helpful to understand more about the nature of the species’ habitat. When extremely simple information is used in model applications on the presence/absence of a species, the results might not accurately reflect the differences in dominance by elevation for species such as A. koreana. Therefore, it is crucial for high-resolution prediction experiments in extremely narrow areas, such as our study area, to pay adequate attention to the existence of microhabitats for in situ species conservation.
Further, the distribution of A. koreana varied depending on the slope (Figure 2b). The maximum slope in the area was approximately 40° (based on a resolution of 100 m2), and the average slope of Mt. Jiri National Park was 17.39°. A comparison of the distribution of the species with the slope showed that the distribution area of A. koreana increased up to a certain level. For example, at an average slope of 35° or higher, the distribution area increased up to 20% within the same slope range. The distribution areas of all A. koreana dominance levels increased along with an increase in slope; however, the relationship between dominance and slope was unclear. Moreover, no clear relational characteristics were observed between aspects and dominance in the distribution of A. koreana.

2.2. Establishment of the High-Resolution Species Distribution Prediction System

2.2.1. Ensemble Prediction System Using Multiple Species Distribution Models

Species distribution models are used to predict the spatial distribution potential of a species, or its abundance, based on the relationship between the observed presence (or occurrence) and absence (non-occurrence) of a species and environmental variables (such as climate, topography, and hydrology). Such models evaluate the effects of various changes in the environment, including climate change, on the distribution of species [24].
Generally, the accuracy of these models is evaluated using a portion of the occurrence input data used to drive the model. Typically, a ratio of 7:3 is used, with 70% of the data used for model training (fitting) and prediction, and the remaining 30% for evaluating the accuracy of the predicted results. However, this internal validation method has limitations for organisms with relatively few occurrence records, as it inevitably reduces the number of information samples [25,26,27]. Moreover, such little data might not reflect the representativeness of environmental characteristics (climate, terrain, etc.) in the corresponding spatial range [28]. To overcome this limitation, various external validation approaches have been proposed, including independent temporal or polygon data that could reflect the distribution range of a species [29]. Various models are available for predicting the spatial distribution of organisms, broadly classified into statistical techniques (linear regression analysis, multiple regression analysis, etc.) and machine learning techniques (artificial neural network, decision tree technique, etc.). In addition to the uncertainty owing to the physical limitations of species survey information, further uncertainty arising from the technical limitations of the techniques used in each model significantly affect the reliability of the prediction results [30]. Therefore, to minimize the uncertainty of each model, ensemble techniques are widely used to evaluate the performance indicators for multiple model results [31].
Therefore, in this study, we used an ensemble platform for species distribution modeling, BIOMOD species (biomod2 v3.5.1), which allows selective ensemble prediction of ten species distribution models, including three statistical models (generalized linear models, generalized additive models, and multivariate adaptive regression splines) and seven machine learning models (classification tree analysis, flexible discriminant analysis, artificial neural networks, generalized boosted models, surface range envelope, random forest, and maximum entropy algorithm). The BIOMOD ensemble algorithm is based on the evaluation of the predictive performance of each model. After each model is run individually, weights are assigned based on the accuracy evaluation of the prediction results of each model, and the results are combined according to these weights [32].
Training or fitting of models within BIOMOD is carried out individually by each model operating within the system. However, performance diagnosis metrics that determine the optimization of the models are equally applied to all models, such as the area under the curve (AUC) of the receiver operating characteristic (ROC) curve and true skills statistics (TSS). The AUC value is calculated by comparing and validating the predicted distribution probabilities and observed values, representing a curve using a true positive rate (TPR) and false positive rate (FPR) for all threshold values that distinguish predicted results into the presence or absence of a species. The cumulative area distribution value under the ROC curve is used as the final model prediction performance metric, where closer to 100% or 1 represents the ideal prediction performance [33].
TSS is a metric developed to overcome the disadvantage of ROC-AUC, which is influenced by the size of the species distribution area (tendency to decrease as the distribution area widens) [34]. TSS utilizes the matching and non-matching information between the observations used for validation, and it is used in the ensemble of each model in BIOMOD. The weights are assigned and averaged based on the TSS values in the ensemble, after which the ensemble model data are generated. In addition, the TSS value that maximizes the value is also used as the threshold (criterion value). The TSS value ranges from −1 to 1, and the higher the value, the more accurate and reliable is the model.
The future potential habitat distribution is predicted by applying the altered future climate conditions based on the relationship between the distribution of the species derived from the analysis of the potential habitat of current climatic conditions and climate environmental variables. However, due to the characteristics of the species distribution model, only the potential habitat distribution results are derived at a specific point in time. Therefore, to predict the distribution changes over time, the future climate environmental variables are divided into multiple time periods, and the potential habitat is predicted for each period.

2.2.2. Building Input Data

Bioclimatic variables (BIOCLIM) were used as the major environmental model inputs in this study. These variables consider the ecological characteristics of the biological species by combining monthly maximum and minimum temperatures and precipitation values [35]. Other environmental variables could be considered based on habitat characteristics, such as altitude, aspect, slope, and connectivity of habitat patches, as well as hydrological features, depending on the habitat requirements of other species. However, when inputting data for all environmental variables, including climate and environmental factors, a process of selecting input variables based on the correlation analysis results between variables is followed. This process eliminates over- or underestimation of species distribution predictions, which are ascribed to the high correlation between variables, such as between spatially averaged temperature distribution and altitude. This implies that one of the variables is excluded in cases of high correlation between variables. However, this process inevitably involves the researcher’s subjective opinion. Furthermore, even if two variables showed high correlation, one might not be excluded depending on its ecological importance, which could be another source of uncertainty in the model prediction results.
In this study, to select input variables for the model, we analyzed the cross-correlation between 19 bioclimatic variables at a resolution of 100 m in the study region [23], as well as geomorphological variables, such as altitude, slope, and aspect in the same area. The final selected variables consisted of four bioclimatic variables related to temperature, and two bioclimatic variables related to precipitation, slope, and aspect (Table 1). Owing to a high correlation with temperature distribution (BIO1), the elevation data were excluded. Future BIOCLIM data were predicted according to the representative concentration pathway (RCP) 4.5 and 8.5 scenarios of the Intergovernmental Panel on Climate. These pathways were employed to predict the potential habitat of A. koreana using 20-year mean monthly climatology data for four periods. These periods were 2040s (2021–2040), 2060s (2041–2060), 2080s (2061–2080), and 2100s (2081–2100) [23].
The current average annual temperature in the Mt. Jiri subalpine region is ~9.4 °C, which is approximately 2.5 °C lower than the national average temperature in South Korea. According to the RCP 4.5 scenario, the temperature is expected to increase at a rate of 0.031 °C/year, reaching an average of 11.9 °C by 2100. According to the RCP 8.5 scenario, the temperature is expected to increase at a rate of 0.056 °C/year, reaching an average of 13.8 °C by 2100 [23].
Precipitation over this region is predicted to increase by approximately 100–200 mm annually from the current (2020) average, which is approximately 10–20% of the mean annual precipitation of 2000 mm. However, precipitation varies over an extensive range, unlike the gradual range of temperature increase. In particular, under the RCP 8.5 scenario, rainfall is expected to decrease sharply to 1850 mm until the 2030s and, subsequently, increase by more than 10% to approximately 2050 mm in the 2050s.

2.3. Experimental Design for Potential Habitat Prediction

We investigated various experimental designs for predicting the high-resolution distribution of the BIOMOD model. The predictive performance of each model result was compared and analyzed by dividing the distribution of A. koreana into four dominance categories. This was carried out to determine how the differences in the altitude distribution of A. koreana by dominance (discussed earlier) would affect the potential habitat prediction within the model. For this purpose, the polygon-formatted A. koreana dominance data were divided into areas of 75% or more, 50% or more, 25% or more, and 10% or more. The corresponding grids in the set dominance range were converted into grid-formatted data at 100 m intervals. For example, if the species dominance of 75% covers majorly in a given grid cell, the cell is assigned as species “presence” of 75% dominance. Then, the central coordinates of those cells are input as the A. koreana presence location information. Eventually, 382 locations with a dominance of 75% or more, 1201 locations with a dominance of 50% or more, 2941 locations with a dominance of 25% or more, and 4101 locations with a dominance of 10% or more were input into the model. In addition, for areas outside the A. koreana occurrence area, random points were extracted and applied to the model as absence data. Finally, the prediction was performed by using all ten species distribution models within BIOMOD and repeated five times to produce the final ensemble prediction result.

3. Results and Discussion

3.1. Prediction Results of Potential Habitats for A. koreana Based on Current Climate Conditions

The prediction results of the potential habitat distribution of A. koreana based on current climate conditions by dominance category indicated wider potential habitat distribution areas with lower dominance. The model performance evaluation (AUC = 0.990–0.993, standard deviation = 0.0015; TSS = 0.896–0.908, standard deviation = 0.0054) was consistently high in all experiments (Figure 3). Although TSS was the highest in Experiment 4 (EXP4), the differences in TSS between each experiment were insignificant.
The percent contribution (PC) and permutation importance (PI) of each input variable for each predictive experiment are shown in Table 2. Contribution refers to the degree to which each variable contributed to the optimal prediction results during the model training process using the observed data. Contribution also represents the relative contribution of variables that could vary depending on the model algorithm [36]. Variable importance measures to what extent a variable influences the dependent variable relative to other variables [36]. In the species distribution prediction experiment for A. koreana, PC and PI showed different patterns, but the climate/environmental variables (BIOCLIM) generally had the highest scores. Specifically, variables related to temperature (BIO1 and BIO4) showed the most influence on the prediction performance in both PC and PI. This was particularly observed in experiments with high dominance, where the weight of precipitation-related variables (BIO14 and BIO15) increased in PC. However, in PI, the model was only sensitive to temperature variables, but in EXP4 the weight of BIO6 and BIO9 increased notably. These variables represent the seasonal characteristics of temperature and precipitation, where BIO6, BIO9, and BIO14 represent the winter climate characteristics of the Mt. Jiri area. The implication was that the dominance of A. koreana in areas with high dominance was closely related to the temperature and precipitation characteristics of winter, which are in turn closely related to the water stress of the biological species. These findings supported the results of previous studies [37] that identified the cause of the population decline of A. koreana in the area. In contrast, the contribution of the topographical variable (slope) decreased notably in areas with high dominance. As the mean temperature variable (BIO1) has a high correlation with altitude, the overall distribution of A. koreana was observed to vary according to mean temperature or altitude. Dominance decreased along with a decrease in altitude, but the species distribution tended to be concentrated in areas with gentle slopes. However, areas with high altitudes and relatively low mean temperatures and where the dominance of A. koreana was high were more sensitive to precipitation (particularly during the dry spring seasons).
The average predicted distribution probability of potential habitats by altitude was observed to increase as the altitude increased (Figure 4). However, the pattern of probability distribution values varied according to altitude between each experiment and the points where the probability converged differed. EXP1 and EXP2 showed peaks around 1500 m, EXP3 around 1600 m, and the probability value of EXP4 continued to increase after 1800 m. However, in all experiments, probability converged at points where distribution probability was 90% or higher. EXP1 and EXP2 converged around 1400 m, EXP3 around 1300 m, and EXP4 around 1500 m. Generally, the height of the distribution probability is considered as the reference point for a potential habitat, which was 75% at around 1100 m on average in EXP1, 2, and 3, and around 1250 m in EXP4.
Furthermore, when the distribution probability was divided into quartiles (25%, 50%, and 75% or higher) for each experiment and the total area of the corresponding potential habitat was examined (Table 3), the potential habitat area decreased along with an increase in the dominance information used. In the experiment where only the information from areas with a dominance of 75% or higher was used, compared with the current distribution of A. koreana, the potential habitat area with a distribution probability of 75% or higher was the closest to the total distribution area of A. koreana (4400 ha, areas with dominance of 10% or higher). This result showed the possibility of overestimating the potential habitat when information from areas with low A. koreana dominance is also used. Therefore, further analysis is required to determine which experiment among the four dominance information utilization experiments conducted in this study could best represent the current distribution of A. koreana. Furthermore, such analysis should explore how the threshold point could be set in the predicted distribution probability for interpreting the potential habitat.
For further examination of the predictive performance of each experiment, we analyzed which experiment best represented the actual distribution of A. koreana by dominance (Table 4). For example, we considered the accuracy of each experiment in predicting the potential habitat in areas with a dominance of 75%. The experiment with the highest accuracy was EXP3. However, interestingly, EXP4 showed the lowest accuracy. This low accuracy could be ascribed to the inability of EXP4 to adequately represent the decrease in A. koreana distribution above an altitude of 1800 m.

3.2. Changes in Potential Habitat Distribution of A. koreana under Future Climate Conditions

The predicted potential habitat distribution results under future climate change scenarios showed a rapid decline in the potential habitat over time in the RCP 4.5 and RCP 8.5 scenarios (Appendix A). These results indicated that A. koreana was extremely sensitive to temperature increases, with the Cheonwangbong area in the eastern high-altitude region of the study area generally predicted to be the last surviving habitat. Considering only the predicted potential habitat results with a distribution probability of 75% or higher, complete extinction was predicted in all the experiments at the end of the 21st century, with the timing of extinction varying depending on the experiment (Figure 5). The rate of potential habitat decline was generally faster in the RCP 8.5 scenario than that in the RCP 4.5 scenario, and this rate was also faster when using information from areas with high dominance, compared to those with low dominance. The timing of extinction obviously differed depending on the distribution probability used as the judgment criterion for potential habitat; however, the phenomenon of rapid decline in potential habitat over time was observed clearly in all areas.
This study, however, had a limitation, namely that only species presence information was employed, regardless of the existence of dead individuals. In the Mt. Jiri subalpine zone, a rapid decline in A. koreana was recently observed [10]. As our study area included the area where this decline was observed, it was necessary to consider the inclusion of information related to this decline in the species distribution modelling. Additionally, the geographical distribution of biological species is directly related to the climatic and environmental factors of their habitat, but the indirect relationship with climatic and environmental factors associated with competition among different species is also an important consideration. Owing to the significant expansion of Sasa borealis in the Mt. Jiri and Mt. Halla subalpine zones, interference with the growth of coniferous trees and the decrease in subsequent tree growth opportunities have become crucial considerations [38]. Accordingly, it was necessary to consider the impact of competition from other species on the potential habitat changes of A. koreana.
The significance of this study is that it was a novel attempt to predict the species distribution of A. koreana using a species distribution model for limited local areas. This was made possible by the increase in spatial resolution of the modeling system to 100 m, which is approximately 100 times higher than that of previous studies. This increase in resolution facilitated prediction experiments for smaller areas. Thereby, more realistic predictions were calculated by reflecting the distribution of A. koreana by altitude, which might be unclear in some previous studies [18]. Our study showed that the distribution pattern of A. koreana is much more sensitive to temperature when detailed altitude characteristics are included in the model, leading to a significant effect of the potential habitat decrease under future temperature increases.

3.3. An Index for Extracting In Situ Conservation Areas for A. koreana

Our modeling results could probably be applied for the identification of conservation areas for A. koreana in the subalpine zone of Mt. Jiri National Park. However, as a rapid decrease in the potential habitat of this species was predicted in the area based on future climate scenarios, exploring alternative habitats in other areas outside Mt. Jiri should be considered. Nevertheless, delaying the extinction of A. koreana in its current habitat for as long as possible would be a valuable conservation measure for preserving ecosystems. Therefore, this study proposes a simple conservation grading method for identifying areas for the in situ conservation of A. koreana in the subalpine region. First, the most rational starting point would be selecting the experiment with the highest prediction performance, considering the minimal uncertainty for the predicted results. Therefore, the results of EXP3 could be applied to the conservation grading of A. koreana, as these best represent the dominance of A. koreana in areas with more than 75% representation. Second, the grading method for conservation could be employed using future prediction results. Collectively, the aim would be to delay the complete extinction of a species occurring in a particular area for as long as possible. Our proposed grading was conducted based on the probability of future distribution presented in the model and information about the areas that could survive the future climate scenarios. Subsequently, starting with areas with the highest grades, employing artificial measures to increase the adaptability of A. koreana would be a possibility. Our method proposes a way to determine the areas with the highest probability for the survival of A. koreana in response to the impacts of climate change. We attempted to create a conservation grading system based on conservative criteria (Table 5) to increase the success rate of A. koreana conservation. The RCP 4.5 scenario was selected to represent the climate change trajectory, and the distribution probability minimum value for potential habitat was set at 20%, with areas with a distribution probability of 50% or more being considered the best fit. Employing the matrix method, the final selected conservation area was classified into four levels: strong recommendation (SR), high recommendation (HR), medium recommendation (MR), and low recommendation (LR). The SR area should be given the highest priority for implementing conservation measures, as it is the area where the survival of A. koreana is most likely to last the longest.
The spatial distribution of each conservation grade is illustrated in Figure 6. Notably, the Cheonwangbong area (42 ha), with elevation >1600 m, was designated with the SR grade. Compared with the current distribution of A. koreana, the SR areas included 31 ha (86% of the entire SR area) of distribution areas with a dominance of 50% or more, as well as some areas where A. koreana was absent (constituting approximately 14% of the entire SR area).
These priority conservation zones derived from the indexing method proposed in this study can be considered a reasonable approach to aid policy makers in the formulation of a species conservation plan for A. koreana. For example, they are useful to set strategic management sites to retard the habitat reduction due to temperature rising. The criteria for distinguishing the recommendation classes, however, can vary according to the importance of policies to conserve A. koreana and the effectiveness of policy implementation. If in situ conservation of the species is of policy importance with very high ecological value, the criteria may need to be set more conservatively. On the contrary, if maximizing the effectiveness of policy implementation is the top priority, then the criteria may need be adjusted to lead to more limited conservation areas.

4. Conclusions

We conducted a model-based study to develop in situ conservation strategies for the native conifer Abies koreana in the subalpine zone of Mt. Jiri that is experiencing severe decline owing to climate change. We constructed a high-resolution species distribution prediction system at a scale of 100 m for precise species distribution prediction in a narrow area, and we performed potential habitat prediction experiments using A. koreana distribution information according to their dominance in Mt. Jiri. The predicted results for potential habitats according to current climatic conditions well represented the distribution characteristics according to temperature and consequently, altitude. The experimental results according to dominance showed that the degree of influence of climatic and environmental factors varied depending on the A. koreana distribution situation, with differences in the contribution or importance of input variables. For example, in areas where A. koreana dominance was high, the altitude of the distribution area was generally high. Furthermore, temperature and precipitation in winter and spring appeared to have a significant influence. Differences in prediction accuracy were observed among the dominance-based experiments. These findings indicate that when a species distribution model is applied in a relatively small area, such as our study area, considering species dominance information could reduce the uncertainty of the predicted results compared with simply using occurrence/non-occurrence information. The predicted A. koreana potential habitat under current climate conditions was slightly wider than the actual distribution area; however, the decrease in potential habitat owing to future climate change was extremely rapid. These findings were ascribed to the potential habitat of A. koreana being susceptible to rising temperature, and the potential habitat with the highest capacity to survive the effects of climate change was predicted to be located mainly in high-altitude areas. Finally, through an indexing method based on the simulation results, this study suggested spatial information about the recommended areas for A. koreana conservation at least to reduce species disappearing over the study region due to warming climate conditions.

Author Contributions

Conceptualization, K.L. and S.H.; Methodology, S.H.; Data curation, J.C.; Formal analyses, K.L., S.H. and D.K., Writing—Original Draft Preparation, K.L. and S.H.; Writing—Review and Editing, K.L. and S.H.; Visualization, S.H. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the projects, “Assessment of climate change risks and adaptation strategy research for ecosystem in Korea”, grant numbers NIE-B-2019-11 and NIE-B-2020-11 and “Study of adaptation capacity to climate change risk for ecosystem”, grant numbers NIE-B-2021-35, NIE-B-2022-35, and NIE-B-2023-35.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Acknowledgments

We would like to thank the Korea National Park Service and Korean National Park Institute for providing the field observation data on Abies koreana from Mt. Jiri.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The following images show the future predictions of A. koreana habitat distribution. EXP1–4, Experiments 1–4; RCP, representative concentration pathway.
Sustainability 15 08964 i001

References

  1. Horikawa, M.; Tsuyama, I.; Matsui, T.; Kominami, Y.; Tanaka, N. Assessing the potential impacts of climate change on the alpine habitat suitability of Japanese stone pine (Pinus pumila). Landsc. Ecol. 2009, 24, 115–128. [Google Scholar] [CrossRef]
  2. IPCC. Summary for policymakers. In Climate Change 2014 Impacts, Adaptation, and Vulnerability. Part A: Global and Sectoral Aspect; Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2014; pp. 1–32. [Google Scholar]
  3. Kong, W.S.; Lim, J.H. Disjunctive distribution of Vaccinium vitis-idaea and thermal condition. J. Korean Geog Soc. 2008, 43, 495–510, (In Korean with English abstract). [Google Scholar]
  4. Kong, W.S.; Kim, K.; Lee, S.; Park, H.; Cho, S.H. Distribution of high mountain plants and species vulnerability against climate change. J. Environ. Impact Assess. 2014, 23, 119–136, (In Korean with English abstract). [Google Scholar] [CrossRef]
  5. Kim, E.S.; Lee, J.S.; Park, G.E.; Lim, J.H. Change of subalpine coniferous forest area over the last 20 years. J. Korean Soc. For. Sci. 2019, 108, 10–20, (In Korean with English abstract). [Google Scholar]
  6. Park, H.C.; Lee, H.Y.; Lee, N.Y.; Lee, H.; Song, J.Y. Survey on the distribution of evergreen conifers in major national parks—A case study on Seoraksan, Odaesan, Taebaeksan, Sobaeksan, Deogyusan, Jirisan National Park. J. Natl. Park. Res. 2019, 10, 224–231. [Google Scholar]
  7. Korea National Park Research Institute. Needleleaf Species Distributions over Subalpine Areas of Korean National Parks; Korea National Park Research Institute: Wonju, Republic of Korea, 2016; p. 79. (In Korean) [Google Scholar]
  8. National Institute of Ecology. Climate Change and Vegetation Colony (Subalpine Coniferous Colony) in Korea I; National Institute of Ecology: Seocheon, Republic of Korea, 2014; p. 134. [Google Scholar]
  9. Millar, C.I.; Rundel, P.W. Subalpine forest. In Ecosystems of California, 1st ed.; Zavaleta, E., Mooney, H., Eds.; University of California Press: Berkeley, CA, USA, 2016; pp. 579–662. [Google Scholar]
  10. Park, H.C.; Moon, G.S.; Lee, H.; Lee, N.Y. A study on the spatial information and location environment of dead coniferous trees in subalpine zone in Jirisan National Park. Korean J. Environ. Ecol. 2020, 34, 42–54, (In Korean with English abstract). [Google Scholar] [CrossRef]
  11. Kim, J.K.; Koh, J.G.; Yim, H.T.; Kim, D.S. Changes of spatial distribution of Korean fir forest in Mt. Hallasan for the past 10 years (2006, 2015). Korean J. Environ. Ecol. 2017, 31, 549–556, (In Korean with English abstract). [Google Scholar] [CrossRef]
  12. Campbell, D.R. Early snowmelt projected to cause population decline in a subalpine plant. Proc. Natl. Acad. Sci. USA 2019, 116, 12901–12906. [Google Scholar] [CrossRef] [PubMed]
  13. Koo, K.A.; Park, W.K.; Kong, W.S. Dendrochronological analysis of Abies koreana W. at Mt. Halla, Korea: Effect of climate change on the growths. Korean J. Ecol. 2011, 24, 281–288, (In Korean with English abstract). [Google Scholar]
  14. Seo, J.W.; Kim, Y.J.; Choi, E.B.; Park, J.H.; Kim, J.H. Investigation of death years and inter-annual growth reduction of Korean firs (Abies koreana) at Yeongsil in Mt. Halla. J. Korean Soc. Environ. Restor. Technol. 2019, 22, 1–14, (In Korean with English abstract). [Google Scholar]
  15. Ahn, U.S.; Kim, D.S.; Yun, Y.S.; Ko, S.H.; Kim, K.S.; Cho, I.S. The inference about the cause of death of Korean fir in Mt. Halla through the analysis of spatial dying pattern—Proposing the possibility of excess soil moisture by climate changes. Korean J. Agric. For. Meteorol. 2019, 21, 1–28, (In Korean with English abstract). [Google Scholar]
  16. Ahn, U.S.; Yun, Y.S. Causes of decline in the Korean fir based on spatial distribution in the Mt. Halla region in Korea: A meta-analysis. Forests 2020, 11, 391. [Google Scholar] [CrossRef]
  17. Kim, N.S.; Han, D.; Cha, J.Y.; Park, Y.S.; Cho, H.J.; Kwon, H.J.; Cho, Y.C.; Oh, S.H.; Lee, C.S. A detection of novel habitats of Abies koreana by using species distribution models (SDMs) and its application for plant conservation. J. Korean Soc. Environ. Restor. Technol. 2015, 18, 135–149, (In Korean with English abstract). [Google Scholar] [CrossRef]
  18. Koo, K.A.; Kim, J.; Kong, W.S.; Jung, H.; Kim, G. Projecting the potential distribution of Abies koreana in Korea under the climate change based on RCP scenarios. J. Korean Soc. Environ. Restor. Technol. 2016, 19, 19–30, (In Korean with English abstract). [Google Scholar] [CrossRef]
  19. Yun, J.H.; Nakao, K.; Tsuyama, I.; Matsui, T.; Park, C.H.; Lee, B.Y.; Tanaka, N. Vulnerability of subalpine fir species to climate change: Using species distribution modeling to assess the future efficiency of current protected areas in the Korean Peninsula. Ecol. Res. 2018, 33, 341–350. [Google Scholar] [CrossRef]
  20. Park, H.C.; Lee, J.H.; Lee, G.G.; Um, G.J. Environmental features of the distribution areas and climate sensitivity assessment of Korean fir and Khinghan fir. J. Environ. Impact Assess. 2015, 24, 260–277, (In Korean with English abstract). [Google Scholar] [CrossRef]
  21. Manzoor, S.A.; Griffiths, G.; Lukac, M. Species distribution model transferability and model grain size—Finer may not always be better. Sci. Rep. 2018, 8, 7168. [Google Scholar] [CrossRef]
  22. Kim, M.K.; Han, M.S.; Jang, D.H.; Baek, S.G.; Lee, W.S.; Kim, Y.H.; Kim, S.J. Production technique of observation grid data of 1 km resolution. J. Clim. Res. 2012, 7, 55–68, (In Korean with English abstract). [Google Scholar]
  23. Hong, S.B.; Lee, K.E.; Kim, M.K.; Sang, J. High-Resolution Bioclimatic Prediction Data Production Technology Guide for the Subalpine Zone—Based on the RCP Climate Change Scenario; National Institute of Ecology (NIE): Maseo-myeon, Republic of Korea, 2022. (In Korean) [Google Scholar]
  24. Austin, M.P. Spatial prediction of species distribution: An interface between ecological theory and statistical modelling. Ecol. Model. 2002, 157, 101–118. [Google Scholar] [CrossRef]
  25. Marmion, M.; Luoto, M.; Heikkinen, R.K.; Thuiller, W. The performance of state-of-the-art modelling techniques depends on geographical distribution of species. Ecol. Model. 2009, 220, 3512–3520. [Google Scholar] [CrossRef]
  26. Buisson, L.; Thuiller, W.; Casajus, N.; Lek, S.; Grenouillet, G. Uncertainty in ensemble forecasting of species distribution. Glob. Chang. Biol. 2010, 16, 1145–1157. [Google Scholar] [CrossRef]
  27. Dobrowski, S.Z.; Thorne, J.H.; Greenberg, J.A.; Safford, H.D.; Mynsberge, A.R.; Crimmins, S.M.; Swanson, A.K. Modeling plant ranges over 75 years of climate change in California.; USA: Temporal transferability and species traits. Ecol. Monogr. 2011, 81, 241–257. [Google Scholar] [CrossRef]
  28. Corcoran, J.; Knight, J.; Pelletier, K.; Rampi, L.; Wang, Y. The effects of point or polygon based training data on RandomForest classification accuracy of wetlands. Remote Sens. 2015, 7, 4002–4025. [Google Scholar] [CrossRef]
  29. Hirzel, A.H.; Lay, G.L.; Helfer, V.; Randin, C.; Guisan, A. Evaluating the ability of habitat suitability models to predict species presences. Ecol. Model. 2006, 199, 142–152. [Google Scholar] [CrossRef]
  30. Araújo, M.B.; New, M. Ensemble forecasting of species distributions. Trends Ecol. Evol. 2007, 22, 42–47. [Google Scholar] [CrossRef]
  31. Elith, J.; Graham, C.H.; Anderson, R.P.; Dudik, M.; Ferrier, S.; Guisan, A.; Hijmans, R.J.; Huettmann, F.; Leathwick, J.R.; Lehmann, A.; et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography 2006, 29, 129–151. [Google Scholar] [CrossRef]
  32. Hao, T.; Elith, J.; Guillera-Arroita, G.; Lahoz-Monfort, J.J. A review of evidence about use and performance of species distribution modeling ensembles like BIOMOD. Divers. Distrib. 2018, 25, 839–852. [Google Scholar] [CrossRef]
  33. Landis, J.R.; Koch, G.G. The measurement of observer agreement for categorical data. Biometrics 1977, 33, 159–174. [Google Scholar] [CrossRef]
  34. Allouche, O.; Tsoar, A.; Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 2006, 43, 1223–1232. [Google Scholar] [CrossRef]
  35. Hijmans, R.J.; Cameron, S.E.; Parra, J.L.; Jones, P.G.; Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. Int. J. Climatol. 2005, 25, 1965–1978. [Google Scholar] [CrossRef]
  36. Baldwin, R.A. Use of maximum entropy modeling in wildlife research. Entropy 2009, 11, 854–866. [Google Scholar] [CrossRef]
  37. Koo, K.A.; Kim, D.B. Review forty-year studies of Korean fir (Abies koreana Wilson). Korean J. Environ. Ecol. 2020, 34, 358–371, (In Korean with English abstract). [Google Scholar] [CrossRef]
  38. Cho, S.; Kim, Y.; Choung, Y. Distribution and synchronized massive flowering of Sasa borealis in the forests of Korean National Parks. J. Ecol. Environ. 2018, 42, 37. [Google Scholar] [CrossRef]
Figure 1. (a) Study area and distribution map of A. koreana according to dominance level (source: Korea National Park Service), photographs of showing (b) the landscape of study area and (c) standing dead trees of A. koreana.
Figure 1. (a) Study area and distribution map of A. koreana according to dominance level (source: Korea National Park Service), photographs of showing (b) the landscape of study area and (c) standing dead trees of A. koreana.
Sustainability 15 08964 g001
Figure 2. Distribution area of A. koreana in the Mt. Jiri National Park area by dominance according to (a) altitude (the black line is the total area of A. koreana) and (b) slope.
Figure 2. Distribution area of A. koreana in the Mt. Jiri National Park area by dominance according to (a) altitude (the black line is the total area of A. koreana) and (b) slope.
Sustainability 15 08964 g002
Figure 3. Prediction results of potential habitats of (A) the experiment 1(EXP1), (B) EXP2, (C) EXP3, and (D) EXP4 for A. koreana based on current climate conditions.; AUC, area under the curve; TSS, true skills statistic.
Figure 3. Prediction results of potential habitats of (A) the experiment 1(EXP1), (B) EXP2, (C) EXP3, and (D) EXP4 for A. koreana based on current climate conditions.; AUC, area under the curve; TSS, true skills statistic.
Sustainability 15 08964 g003
Figure 4. Average distribution probability by experiment according to altitude. EXP1–4, Experiments 1–4.
Figure 4. Average distribution probability by experiment according to altitude. EXP1–4, Experiments 1–4.
Sustainability 15 08964 g004
Figure 5. Distribution probability average of potential habitats by time period. EXP1–4, Experiments 1–4; RCP, representative concentration pathway.
Figure 5. Distribution probability average of potential habitats by time period. EXP1–4, Experiments 1–4; RCP, representative concentration pathway.
Sustainability 15 08964 g005
Figure 6. Map of conservation grades of A. koreana that are color-coded to represent strong, highly, medium, and low recommendations: SR, HR, MR, and LR, respectively.
Figure 6. Map of conservation grades of A. koreana that are color-coded to represent strong, highly, medium, and low recommendations: SR, HR, MR, and LR, respectively.
Sustainability 15 08964 g006
Table 1. Selected BIOCLIM variables that were used as the model inputs.
Table 1. Selected BIOCLIM variables that were used as the model inputs.
BIOCLIMDescriptionUnit
BIO1Annual mean temperature°C
BIO4Temperature seasonality
BIO6Mean minimum temperature of the coldest month°C
BIO9Mean temperature of the driest quarter°C
BIO14Precipitation in the driest monthmm/month
BIO15Precipitation seasonality
BIOCLIM, Bioclimatic variables.
Table 2. Variable contribution of the current potential habitat prediction results of A. koreana according to dominance.
Table 2. Variable contribution of the current potential habitat prediction results of A. koreana according to dominance.
BIOCLIM
Variables
EXP1 (≥10%)EXP2 (≥25%)EXP3 (≥50%) EXP4 (≥75%)
PCPIPCPIPCPIPCPI
BIO16877.121.4722777.516.322
BIO40.78.827.513.646.31.72228.6
BIO60.1100009.122.9
BIO9000.20.70.96.69.218.9
BIO140.48.60.88.21.29.926.13.5
BIO150.10.29.43.924.13.215.22.9
Slope30.63.920.71.70.401.20.9
Aspect00.100000.80.4
EXP1–4, experiments 1–4; BIOCLIM, Bioclimatic variables; PC, percentage contribution; PI, permutation importance.
Table 3. Total area of A. koreana distribution estimated in each experiment according to inhabitable probabilities.
Table 3. Total area of A. koreana distribution estimated in each experiment according to inhabitable probabilities.
Inhabitable ProbabilityEXP1EXP2EXP3EXP4
≥25%16,061 ha14,415 ha13,032 ha11,013 ha
≥50%12,806 ha11,175 ha10,693 ha8094 ha
≥75%9316 ha8045 ha8234 ha5449 ha
EXP1–4, Experiments 1–4.
Table 4. Average inhabitable probabilities in the areas of the actual distributions of A. koreana.
Table 4. Average inhabitable probabilities in the areas of the actual distributions of A. koreana.
Density of
A. koreana
Average Inhabitable Probability
EXP1EXP2EXP3EXP4
75–100%89.9688.9192.1488.29
50–75%89.6088.8291.3384.82
25–50%86.6585.4485.3377.33
10–25%85.1081.7582.3971.97
EXP1–4, Experiments 1–4.
Table 5. Matrix approach for the grading of A. koreana conservation to mitigate the effects of climate change.
Table 5. Matrix approach for the grading of A. koreana conservation to mitigate the effects of climate change.
Inhabitable ProbabilityResidual PeriodConservation Grade
2100208020602040GradeDescription
>50%SRHRHRMRSRStrong recommendation
40–50%HRHRMRMRHRHigh recommendation
30–40%HRMRMRLRMRMedium recommendation
20–30%MRMRLRLRLRLow recommendation
SR, strong recommendation; HR, high recommendation; MR, medium recommendation; LR, low recommendation.
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

Lee, K.; Kim, D.; Cha, J.; Hong, S. Fine-Scale Species Distribution Modeling of Abies koreana across a Subalpine Zone in South Korea for In Situ Species Conservation. Sustainability 2023, 15, 8964. https://doi.org/10.3390/su15118964

AMA Style

Lee K, Kim D, Cha J, Hong S. Fine-Scale Species Distribution Modeling of Abies koreana across a Subalpine Zone in South Korea for In Situ Species Conservation. Sustainability. 2023; 15(11):8964. https://doi.org/10.3390/su15118964

Chicago/Turabian Style

Lee, Kyungeun, Daeguen Kim, Jaegyu Cha, and Seungbum Hong. 2023. "Fine-Scale Species Distribution Modeling of Abies koreana across a Subalpine Zone in South Korea for In Situ Species Conservation" Sustainability 15, no. 11: 8964. https://doi.org/10.3390/su15118964

APA Style

Lee, K., Kim, D., Cha, J., & Hong, S. (2023). Fine-Scale Species Distribution Modeling of Abies koreana across a Subalpine Zone in South Korea for In Situ Species Conservation. Sustainability, 15(11), 8964. https://doi.org/10.3390/su15118964

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