서 론

최근 전 세계적으로 가속화되는 기후변화로 인해 다양 한 사회·환경 문제들이 대두되고 있다. 기후변화에 관한 정부 간 협의체 (IPCC, Intergovernmental Panel on Climate Change)에서 발간한 지구온난화 1.5°C 특별보고서에서는 최근 평균온도가 산업화 이전 대비 약 1.0°C가 증가하였으 며, 2052년에는 1.5°C가 상승할 것이라 제시하였다 (IPCC 2018). 최근에 발간된 IPCC 제6차 보고서에서는 2040년 무렵 1.5°C가 증가하여 앞서 발간한 특별보고서에 비해 10 년 이상이 감축되어 지구온난화로 인해 기후변화가 급격 하게 진행되고 있다고 언급하였다 (IPCC 2021).

우리나라의 경우 2022년의 사례만 보아도 많은 기후재 해가 발생하였다. 충청 및 영남지방의 가뭄과 수도권 및 강 원지역의 집중호우, 울진의 산불 등이 인간과 자연환경에 큰 영향을 주었다. 특히, 가뭄과 홍수 등 수문학적 재해는 산불과 산사태로 이전되어 산림생태계에 많은 교란을 야 기하고 있다 (Lim et al. 2019;Lim and Kim 2022). 폭염과 국 지적 한파, 지속적 온난화는 수종별 생장이나 서식처 질에 도 많은 영향을 준다 (Yoon et al. 2022). 즉, 넓은 토지공간 을 차지하여 관리와 적응이 제한적인 산림생태계는 기후 변화에 완전하게 노출되어 있다고 볼 수 있다. 장기적으로 도 극한기상현상의 증가는 산림생태계에 교란 요인으로 작용하여 미래의 산림생태계에 큰 위협이 될 것으로 예상 된다 (Lee 2022). 현재는 아고산 지역의 침엽수림이 산림생 태계의 기후변화 영향을 가장 적나라하게 보여주고 있는 데 (Lim and Kim 2020;Yoo et al. 2020), 기후변화가 장기화 될수록 그 영향은 한반도 전역으로 확대될 것으로 전망되 고 있다 (Lim et al. 2018).

우리나라는 산림분야에서도 기후변화 대응을 위한 다양 한 정책을 수립하였다. 그중 6차 산림기본계획에서는 산림 자원 및 산지관리체계 고도화를 통한 산림 역할 강화를 제 시하고 있다. 특히, 국가 온실가스 감축 목표 달성을 위해 경제림 육성단지를 거점으로 탄소 흡수력이 높은 수종으 로 갱신하고, 국산 목재 이용의 확대를 계획하고 있다 (KFS 2018). 수목을 가구재, 건축재 등 목재로 이용하면 탄소를 고정하는 역할을 하게 되며 이는 온실가스 감축에 기여할 수 있다 (Zhang et al. 2020;Lim and Choi 2021). 이에 6차 산림기본계획에서는 목재를 활용한 기후변화 대응을 위해 2030년까지 2015년도 목재생산량보다 70% 이상 증가하는 것을 목표로 삼고 있다.

다만, 기후변화의 산림생태계 영향이 목재를 생산하는 용재수종에도 유효하다는 점을 고려해야 한다. 우리나라 평균기온이 1°C 상승할 경우 용재수종으로 사용되는 졸 참나무, 서어나무, 개서어나무의 점유율은 증가하지만 잣 나무와 신갈나무는 감소할 것이라 전망하였다 (Ahn et al. 2009). 특히, 현재 아고산 침엽수림의 서식공간이 급격히 감소하는 것처럼, 장기전망에서는 현재 중부지방에 식재되 는 침엽수림의 서식공간이 크게 감소할 것으로 전망되고 있다 (Cui et al. 2016). 목재생산량을 증대시키기 위해서는 변화하는 기후에 적합한 용재수종을 선정하고 식재하는 과정이 필요한 것이다.

이에 본 연구에서는 국내에서 목재로 생산하고 있는 용 재수종들을 대상으로 신기후변화 시나리오에 따른 서식지 적합도를 평가하고자 한다. 기계학습기반의 종분포모델 앙 상블 기법에 기후·지형인자 및 수종별 위치정보를 활용하 여 기후변화로 인한 서식지 적합도의 변화를 예측한다. 특 히 현재까지 중점적으로 조림되던 용재수종과 난대·아열 대 용재수종을 함께 모의하여 앞으로의 수종대체 가능성 을 확인하고자 한다. 이러한 접근은 기존에 조림되고 있던 용재수종의 난대·아열대 수종으로 대체 가능성을 평가하 고 앞으로의 용재수종 선정에 대한 과학적 정보를 제공할 수 있을 것이다.

재료 및 방법

1. 연구대상지 및 범위

본 연구는 우리나라 전체 국토를 대상으로 수행하였다. 우리나라는 북위 33°~39°, 동경 124°~131° 사이에 위치하 며 전체 면적은 약 100,412 km²다. 전체 면적 중 약 63.3% 가 산림이며 대부분 산은 동부와 북부에 위치한다. 최근 우 리나라의 연 평균기온은 약 13.3°C, 연 강수량은 약 1,300 mm, 평균 상대습도는 68%이다 (KFS 2021). 다만, 기후변 화에 따라 최근에도 기온의 상승과 강수계절성 변화가 관 측되어 산림서식처에 영향을 주고 있으며 (Lim and Kim 2020), 미래 기후전망에서는 연 평균기온이 최대 7°C까지 상승할 수 있는 것으로 예상되고 있다 (KMA 2021).

분석에 사용한 데이터는 기후자료의 공간해상도를 고려 해 전체 공간해상도는 1 km²로 설정하였으며, 대상 시기는 2000~2019년을 현재 시기로, 2041~2060년을 근미래로, 2081~2100년을 먼미래의 기후변화에 따른 영향의 시기로 선정하였다.

2. 수종 위치정보

본 연구에서는 기후변화로 인한 우리나라 산림의 수 종별 서식지역 변화에 대응하기 위하여 세 가지 용재수 종을 활용하였다. 현재 기후에 적합한 수종의 경우 전국 적으로 조림 및 분포하며 조림되고 있는 잣나무 (Pinus koraiensis)와 미래 용재수종으로는 제주도와 남해안 도서 지역에서 주로 조림 및 분포하고 있는 삼나무 (Cryptomeria japonica)와 편백 (Chamaecyparis obtusa)을 선정하였다 (Fig. 1). 세 수종의 위치정보 추출을 위해 전국 약 4,500개 의 표본점을 조사하는 국가산림자원조사 (NFI; National Forest Inventory) 자료 중 가장 최근 자료인 제7차 NFI (2016~2020) 자료에서 기후적합도를 고려하여 천연림에 서만 잣나무와 삼나무, 편백의 위치정보를 추출하였다. 잣 나무는 침엽수종 중 소나무와 편백, 낙엽송과 더불어 많이 조림되는 침엽수종이며 (KMS 2021) 기후변화의 영향으로 피해를 입을 가능성이 높은 기후변화 취약 산림식물종 100 종에 포함되어 있다 (Lim et al. 2016). 편백은 피톤치드로 인한 독특한 향과 영향으로 인해 목재와 건축재로 수요가 많아 산지 소유주들이 조림을 희망하는 인기가 많은 용재 수종 중 하나이며 ( Joung et al. 2020), 삼나무 또한 피톤치드 의 함유량이 높으며 방풍림으로도 사용되고 있는 수종이 다 (Kim et al. 2013;Kim et al. 2017). 세 수종은 모두 상록침 엽교목이며 산림청에서 선정한 권장 용재수종이라는 공통 점을 가지고 있다.

3. 환경 변수

본 연구에서 기후변수로는 19가지 생물기후변수와 WI (Warmth Index)를 이용하였다. 생물기후변수를 제작하기 위해선 최고기온, 최저기온, 누적 강수량에 대한 자료가 필 요하다. 현재시기의 기후정보인 기상청에서 제작한 MKPRISM v2.1 자료와 근미래와 먼미래의 기후변화 시나리 오인 SSP (Shared Socioeconomic Pathways)를 기상청이 운 영하는 기후정보포털을 통해 제공받았다. MK-PRISM v2.1 자료는 종관기상관측장비 (Automated Synoptic Observing System)와 방재기상관측장비 (Automatic Weather Station) 의 기상관측자료를 PRIDE (PRISM based Dynamic downscaling Error correction) 모델에 적용하여 2000~2019년의 일별 및 월별, 연간의 격자 단위 기후자료를 제공한다. MKPRISM (Modified Korea - Parameter-elevation Regressions on Independent Slopes Model) 기법은 고도, 거리, 지향면, 해양도의 영향을 고려하여 고해상도의 격자 자료를 산출 하는 PRISM (Parameter-elevation Regression on Independent Slopes Model)을 대한민국에 맞게 수정한 PRISM이 며 PRIDE 모델은 관측 해상도보다 더 고해상도로 기후변 화 시나리오를 통계적으로 상세화하는 모델이다 (Lim and Kim 2022). IPCC 5차 평가보고서에 사용된 RCP (Representative Concentration Pathways) 시나리오는 온실가스 농 도의 변화만 고려하여 제작되었지만 6차 평가보고서에 기 반한 SSP 시나리오는 인구통계, 경제발달, 복지, 생태계 요 소, 자원, 제도, 기술 발달, 사회적 인자, 정책 등 미래사회경 제와 온실가스 농도 전망을 함께 반영하여 보다 미래 전망 에 대해 현실성을 개선한 시나리오이다. SSP 시나리오 데 이터는 영국 기상청 해들리센터의 전지구기후모델 (Had GEM2-AO)을 기반으로, 지역기후모델인 HadGEM4-RA, WRF, CCLM, GRIMs, RegCM4의 앙상블 결과를 사용하였 다. 본 연구에서는 재생에너지 기술의 발달로 온실가스 감 축이 성공적으로 수행되어 화석연료 사용이 최소화되고 친환경적이고 지속가능한 경제 성장을 가정하는 SSP1-2.6 시나리오와 온실가스 감축이 미이행되어 화석연료 사용 이 증가하고 도시 위주의 개발이 확대될 것을 가정하는 SSP5-8.5의 시나리오를 사용하였다. SSP1-2.6과 SSP5-8.5 시나리오는 온실가스 감축과 사회변화에 대해 가장 큰 차 이를 볼 수 있는 시나리오이다. 특히 SSP1-2.6 시나리오 는 IPCC가 2100년까지 지구 평균온도 상승폭을 1.5°C로 제한하기 위해서는 2050년까지 전 지구적으로 탄소중립 을 실현해야 한다고 제시하였으므로 (IPCC 2021), 2050 탄소중립이 실현되었을 경우와 가장 가까운 시나리오이다.

위의 기후정보를 활용하여 19가지 생물기후변수를 제작 하였다. 생물기후변수는 식물의 분포 및 생장에 관계성이 높은 기후요소로 구성된 변수들로, 이를 활용하여 수종의 분포와 서식처를 예측하는 많은 선행 연구가 존재하였다 (Chun et al. 2015;Lee et al. 2021). R프로그래밍 언어 (Version 3.6.1)에서 R Studio의 ‘dismo’ 패키지를 이용하여 19가 지 생물기후변수를 제작하였으며 변수 간 다중공선성 (Multicollinearity)을 제거하기 위해 Pearson 상관관계분석 을 시행하였다. 변수 간 연관성 (상관계수 0.9 이상)이 높은 변수들은 모델 적용에서 제외하였고, 상관계수가 높은 두 변수 중 선행 연구에서 식물의 분포에 영향이 높은 것으로 제시된 변수를 선정하였다. 선정된 생물기후변수는 Table 1 과 같다.

추가적으로 WI를 분석에 포함하였다. WI는 월평균기온 5°C 이상인 달에 대하여 월평균기온과 5°C와의 차를 1년 동안 합한 값으로, 식물 생장에는 일정량 이상의 열이 필 요하다는 전제로 제작되었다. 기온을 이용하여 식물의 분 포 및 이동을 설명할 수 있는 대표적 지수이다 (Park et al. 2012).

지형 자료의 경우, Worldclim (www.worldclim.org)에서 제공하는 수치표고모델을 토대로 지형위치지수 (Topographic Position Index)와 지형습윤지수 (Topographic Wetness Index)를 제작하여 사용하였다. 토양 자료는 ISRIC (data. isric.org)에서 제공하는 토양 pH를 사용하였다.

4. 종 분포 예측 모델 및 앙상블

본 연구에서는 기계학습 알고리즘을 토대로 용재수종의 출현지점과 생물기후지수를 학습시켜 적합 서식지를 예측 하였다. 분석에 사용한 기계학습 알고리즘은 BIOMOD2 모형에서 제공하는 통계 기반의 Generalized Linear Models (GLM), Multivariate Adaptive Regression Splines (MARS)와 기계학습 기반의 Random Forest (RF), Maximum Entropy (Maxent), Generalized Boosted Models (GBM), Artificial Neural Networks (ANN), Classification Tree Analysis (CTA), Flexible Discriminant Analysis (FDA)로 총 8가지의 모델을 분석에 이용하였다. MaxEnt 모델을 제외한 나머지 7개의 모델은 모두 종의 출현 - 비출현 형태의 자료가 필요한 모 형이다. 실제 출현 - 비출현좌표를 이용하는 것이 정확도 가 높다는 선행 연구에 따라 NFI의 출현좌표 중 대상 수종 을 제외한 출현좌표를 비출현좌표로 선정하였다 (Seo et al. 2008;Lee et al. 2021). 이후 선정된 비출현좌표를 출현좌표 와 비출현좌표의 비율이 3대 7이 되도록 무작위 추출하였 다. 8개의 모델을 10회 반복 시행하여 AUC 값이 0.7 이상 인 결과와 CV (Committee Averaging) 앙상블 방법론을 기 반으로 앙상블을 진행하였다. 앙상블은 종 분포 단일 모형 의 오차 및 불확실성을 보완할 수 있는 방법론으로 많은 선 행 연구에서 이용된 검증된 방법론이다 (Park et al. 2016;Noce et al. 2017;Khan and Verma 2022). 종 분포 모델과 앙상블 모델에 대한 정확도를 ROC (Receiver Operating Characteristic) 곡선의 AUC (Area Under the Curve) 값을 통해 측정하였다. AUC 값은 0에서 1 사이의 값으로, 1에 가 까울수록 정확도가 높으며 AUC가 0.7 이상이면 모델의 설 명력이 좋은 것을 의미한다 (Lee et al. 2018). 모형 검증은 훈련자료 (training data)와 시험자료 (test data)를 80대 20으 로 나누어 3회 반복 시행하였다. 단일 모델 및 앙상블한 결 과는 서식 적합성이 0과 1 사이의 값으로 나타나며 적합성 을 동일 간격 (equal interval)을 통해 5개의 등급으로 구분 하였다. 0.2 이하는 서식지로 부적합한 지역, 0.2 초과 0.4 이 하는 서식지 적합도가 낮은 지역, 0.4 초과 0.6 이하는 서식 지로 적합한 지역, 0.8 초과 1 이하는 서식지로 매우 적합한 지역으로 분류하였다 (Li et al. 2020).

결과 및 고찰

1. 모델의 정확도 및 인자별 중요도 평가

8가지 모델의 10회 반복 시행 결과는 아래 Table 2와 같 다. 모델들 중 비교적 CTA와 ANN이 낮게 평가되고 RF, MaxEnt, GBM, GLM이 높게 평가되었다. 8가지 모델을 앙 상블한 모델의 경우 기존 모델들보다 AUC가 높게 나타나 단일 모형의 오차 및 불확실성을 보완했다. 특히, 잣나무의 경우 8가지 모델의 AUC 평균은 0.737이었지만 앙상블한 결과 AUC가 0.898까지 증가하게 되어 오차 및 불확실성을 보완하였다고 볼 수 있다. 잣나무의 경우 출현자료가 중부 지방에서 넓게 분포하여 상대적으로 AUC가 낮게 나타났 으며, 삼나무와 편백의 경우 출현자료가 남부지역에 국한 되어 높은 AUC를 보였다.

수종별 서식지 분포에 영향을 미치는 인자의 중요도를 파 악한 표는 Table 3과 같다. 인자의 중요도는 8종의 모델을 10 번 반복 시행한 중요도를 평균으로 나타낸 값으로 상대적인 값을 나타낸다. 잣나무는 연간 온도 범위 (BIO7), WI, 가장 추운 달의 최소 온도 (BIO6) 순으로 높은 중요도를 나타냈 으며 온도에 민감하게 반응하였다. 삼나무는 가장 추운 분기 의 강수량 (BIO19), 가장 따뜻한 분기의 강수량 (BIO18), 강 수 계절성 (BIO15) 순으로 중요도가 높게 나타났으며 편백 은 가장 추운 분기의 강수량 (BIO19), 강수 계절성 (BIO15), 가장 따뜻한 분기의 강수량 (BIO18) 순으로 높은 중요도를 나타냈다. 난대·아열대 수종인 삼나무와 편백은 온대 수종 인 잣나무에 비해 강수량에 민감하게 반응하였다.

2. 수종별 서식지 분포 변화

1) 잣나무

잣나무의 잠재 서식지의 예측 결과는 Fig. 2와 같다. 잣 나무의 현재 잠재 서식지는 내륙지방 중 북부에 서식에 적 합한 지역이 분포하고 있으며 북부를 제외한 지역은 적합 도가 낮게 나타났다. SSP1-2.6과 현재를 비교하면 근미래 에는 북부지방을 중심으로 매우 적합 서식지가 확산하지 만 먼미래에는 확산되었던 ‘매우 적합’ 서식지가 축소되 고 남부지역의 서식지가 감소하면서 경기도, 강원도 지역 에만 적합 서식지가 분포한다. SSP5-8.5와 현재를 비교하 면 근미래에는 매우 적합한 서식지의 면적이 줄어들게 되 며 서식지로 부적합한 지역이 중부에서 북상하여 경기도 와 강원도 부근에만 적합 서식지가 분포하게 된다. 먼미래 에는 서식지로 매우 적합한 지역은 대부분이 사라지게 된 다. 미래로 갈수록 침엽수의 적합 서식지가 대부분 사라지 고 강원도 고산지대 일부에만 분포가 가능하다는 선행 연 구 (Yoo et al. 2020)와 유사한 결과가 나타났다. 잣나무의 시 나리오에 따른 시기별 면적 변화는 Table 4와 같다. SSP1- 2.6에서 적합 서식지의 면적은 현재의 적정 면적과 비교해 근미래에 약 1.447배가 증가하지만 먼미래에 약 0.803배로 감소할 것으로 예상되었다. 반면, SSP5-8.5의 적합 서식지 면적은 근미래에 약 0.665배이며 먼미래는 0.008배로 급격 하게 감소할 것으로 전망되었다.

2) 삼나무

삼나무의 잠재 서식지의 예측 결과는 Fig. 3과 같다. 삼나 무의 현재 잠재 서식지는 한라산 고산지대를 제외한 제주 도와 남부 도서지역, 내륙의 남부지역에 매우 적합한 지역 이 분포하고 있으며, 내륙 중부지방부터는 서식지로 부적 합한 지역이 넓게 분포한다. SSP1-2.6과 현재를 비교하면 근미래에는 남서쪽을 중심으로 적합지역이 지역이 북상하 고 있으며 먼미래에는 서식에 적합한 지역이 서해안과 동 해안을 중심으로 북상하였다. SSP5-8.5 시나리오에서 근미 래에는 서해안과 동해안을 중심으로 서식 적합지역이 확 장하여 강원도 해안지역까지 북상하였으며, 내륙지방으로 도 확산되었다. 먼미래에는 서식지 적합도가 높은 지역은 근미래와 비슷한 면적이지만 적합도가 낮은 지역이 크게 감소하여 경기·강원 산간지역 일부에서만 분포하였다. 삼 나무의 시나리오에 따른 시기별 면적 변화는 Table 4와 같 다. SSP1-2.6에서 적합 서식지의 면적은 현재의 적합 면적 과 비교하여 근미래에 약 0.734배로 감소하지만 먼미래에 는 약 2.765배로 급격하게 증가한다. SSP5-8.5의 적합 서식 지 면적은 근미래에 약 2.856배이며 먼미래는 2.570배로 근미래에 비해 면적이 감소할 것으로 전망되었다. SSP5-8.5 먼미래에 근미래보다 적합 서식지의 면적이 적어지는 이 유는 수분 조건이 양호할 때 적합 서식지로 판정될 가능성 이 높다고 평가한 선행 연구 (Chun et al. 2014)와 강수량의 중요도가 높게 나타난 결과를 비추어 볼 때 제주도를 비롯 한 남부 해안지역의 기온이 지속적으로 상승함에 따라 열 대화가 진행되면서 삼나무의 수분 조건에 부정적 영향을 미쳐 적합 서식지의 면적이 감소한 것으로 보인다.

3) 편백

편백의 잠재 서식지 예측 결과는 Fig. 4와 같다. 편백의 현재 잠재 서식지는 한라산 고산지대를 제외한 제주도와 남부 도서지역, 내륙의 남부지역에 매우 적합한 지역이 분 포하고 있으며 내륙의 중부지방부터 서식지로 부적합하다. SSP1-2.6과 현재를 비교하면 근미래에는 매우 적합한 지역 이 북상하였으며, 부적합지역이 넓었던 중부지방 또한 낮 은 수준 적합도로 상당 지역 변화하였다. 먼미래에는 중부 지방의 적합도도 매우 높거나 높은 수준으로 증가할 것으 로 나타났다. SSP5-8.5 시나리오에서 근미래에는 강원도 와 경기 북부를 제외한 대부분의 지역이 매우 적합한 지역 으로 확산되었다. 먼미래에는 일부 고산지역을 제외한 우 리나라 대부분 지역에서 서식 적합도가 높게 전망되었다. 편백의 시나리오에 따른 시기별 면적 변화는 Table 4와 같 다. SSP1-2.6의 적합 서식지의 면적은 현재의 적정 면적 과 비교해 근미래에 약 1.636배가 증가하며 먼미래에는 약 2.602배로 급격하게 증가한다. SSP5-8.5의 적합 서식지 면 적은 근미래에 약 3.460배이며 먼미래는 4.705배로 SSP1- 2.6에 비해 넓은 면적이 크게 증가할 것으로 전망되었다.

3. 난대 용재수종의 잣나무 대체 가능성

기존 전국적으로 활용되던 용재수종을 난대·아열대 수 종이 대체할 수 있는지 확인하기 위해 현재 잣나무의 잠재 서식지와 난대·아열대 수종의 서식 적합지역을 중첩하여 대체 가능성을 분석하였다 (Fig. 5). 또한, 현재 잣나무 적합 서식지 면적에서 미래 난대·아열대 수종의 서식 적합지를 통해 대체할 수 있는 면적을 확인하였다 (Table 5). 현재와 SSP1-2.6의 잣나무와 삼나무의 적합 서식지 중 편백과 중 복되는 지역이 거의 존재하지 않았으며 SSP1-2.6의 먼미래 에는 삼나무와 편백이 수도권 지역 중 일부에 서식이 적합 해지면서 현재 잣나무의 적합 서식지 중 약 5.22%를 대체 할 수 있는 것으로 확인되었다. SSP5-8.5 시나리오의 근미 래에는 SSP1-2.6의 먼미래에 비해 상대적으로 수도권에서 삼나무와 편백의 적합 서식지가 증가하였으며 현재 잣나 무의 적합 서식지 중 약 24.30%를 대체할 수 있을 것으로 나타났다. SSP5-8.5의 먼미래에서는 경우 백두대간 부근을 제외하고 난대·아열대 수종이 대체 가능하였으며 현재 잣 나무의 적합 서식지 중 약 85.39%를 대체할 수 있는 것으 로 전망되었다.

4. 함의 및 한계

용재수종의 조림에는 일정량의 기간이 필요하며 최근 기온이 급격하게 변하는 상황에서 조림수종의 선정에 유 의할 필요가 있다. SSP 시나리오 중 온실가스 배출량이 적 은 경우인 SSP1-2.6과 온실가스 배출량이 많은 경우인 SSP5-8.5 두 가지 시나리오로 분석을 진행함으로써 온실가 스 배출량에 따라 용재수종 및 적합지를 분석하였다. 특히, SSP1-2.6에서는 잣나무의 대체 필요성이 높지 않았으나, SSP5-8.5에서는 서식 적합지역이 크게 감소하며, 대체 요 구도가 높았다. 이는 전 세계의 탄소중립과 온실가스 감축 에 따라 용재수종의 대체 여부도 고려해야 함을 시사한다.

수종별 예측 결과에서 잣나무의 경우 SSP1-2.6의 근미래 를 제외하고는 적합 서식지가 현재에 비해 감소하고 있으 며 다양한 선행 연구들 역시 기후변화로 인해 침엽수종의 서식적합이 지속적으로 감소하고 있다는 경향을 보이고 있다 (Park et al. 2014;Lim and Kim 2020;Yoo et al. 2020). 이는 21세기 내에 기존의 주요 용재수종 중 일부는 더 이 상 국내에 조림이 어렵다는 것을 의미한다. 반면, 미래에 대 체 가능한 용재수종으로 제시한 편백과 삼나무의 경우, 서 식 적합지 북상으로 기존 용재수종의 적합지 중 80% 이상 을 대체가 가능할 것으로 확인된 것이다. 편백과 삼나무의 서식 적합지역 확대는 유사한 방향성을 보이나 적합도에 서 차이를 보였다. 편백은 SSP5-8.5의 먼미래에 전국의 대 부분에서 서식이 적합하게 확인되었으나, 삼나무 경우 우 리나라 중부지역까지만 적합지가 높게 나타났다. 두 가지 난대·아열대 용재수종에서도 차별성을 발견한 것이며, 온 실가스 감축과 미래 기후변화에 따라 유사 수종군에서도 적정 용재수종의 선정이 필요할 것을 시사한다.

다만, 이번 분석이 가지는 몇 가지 한계점이 존재한다. 현 재 분석에 이용한 출현좌표는 NFI의 자료를 이용하였다. NFI의 자료는 표본점 조사를 통한 자료이므로 실제 분포 하고 있는 수종의 모든 좌표를 반영할 수 없으며 실제 생 물의 분포 경향과는 차이점이 존재한다. 또한, 분석에 사용 한 자료는 비생물적 요인인 기후요소와 지형요소만 적용 되었다. 경쟁, 포식, 공생 또는 종간 관계와 같은 생물적 요 인은 식물의 분포에 영향을 주므로 분석 결과는 실제 분포 와 다소 다른 결과로 나타났을 가능성이 있다. 이번 연구에 서는 기존의 주요 용재수종 중 한 수종에 대한 분석만 진행 하여 우리나라의 모든 기후대의 용재수종에 대한 반영이 되지 못했으므로 추후 다양한 기후대의 용재수종의 적합 지 분포를 분석할 필요가 있다. 또한, SSP 시나리오 중 기존 RCP에서 국제적으로 BAU (Business As Usual)로 판단하는 SSP2-4.5와 SSP3-7.0을 이용한 예측을 진행하지 못하였으 므로 추후 모든 SSP 시나리오를 이용한 용재수종의 적합지 분포를 분석할 필요가 있다.

적 요

본 연구는 국내 주요 용재수종인 잣나무와 삼나무, 편백 에 대한 종 분포 예측 모델의 결과를 앙상블하여 기후 시나 리오에 따라 현재, 근미래, 먼미래의 서식 적합지를 예측하 였고 잣나무와 삼나무, 편백의 기후변화 시나리오별 분포 적합지를 분석하였다. 특히, 잣나무를 삼나무와 편백이 대 체할 수 있는지 평가하였다. 기준연도 (현재) 잣나무의 매 우 적합한 서식지는 전국의 약 13.87%를 차지하지만 SSP5- 8.5 하의 먼미래에서는 약 0.11%까지 낮아진다. 삼나무의 경우 기준연도의 서식 적합지는 약 7.08%이며 SSP5-8.5 하 의 먼미래에서는 약 18.21%까지 증가한다. 편백의 경우 기 준연도의 서식 적합지는 약 19.32%이며 SSP5-8.5 하의 먼 미래에서는 약 90.93%까지 차지하는 것으로 예측되었다. 전국적으로 조림하던 잣나무는 기후변화의 영향으로 서식 처가 점차 북상하여 우리나라에서 적합한 서식처가 크게 감소하였으므로 21세기 중반 이후에는 국내에서 용재수종 으로 조림하기에는 부적합하며 높은 수준의 서식 적합도 를 갖는 편백이나 삼나무로 대체될 필요가 있다고 전망된 다. 특히, 편백은 대부분의 영역에서 잣나무를 대체 가능하 며, 삼나무의 경우 남해안과 중부지방 일부분을 대체할 수 있다고 평가되었다. 결론적으로 미래에는 조림하는 용재수 종의 변화가 생길 것이며 다양한 수종을 대상으로 한 연구 를 통해 기후변화에 대응하는 방안이 마련되기를 기대한다