기후변화에 따른 대한민국의 표범장지뱀 서식지 변화 예측 = Projecting Habitat Shifts of Eremias argus under Future Climate Change Scenarios in the Republic of Korea

국문 초록 (Abstract)

본 연구는 기후변화에 따른 대한민국 표범장지뱀（Eremias argus）의 현재 및 미래 잠재 서식지 변화를 정량적으로 평가하고, 복수의 미래 조건에서 서식지 감소가 반복적으로 나타나는 지역을 도출하는 것을 목적으로 하였다. 이를 위해 2007–2018년의 출현자료 23개 지점과 기후, 지형, 토지피복·식생 및 거리 기반 환경변수를 구축하였으며, 다중공선성 검토를 거쳐 14개의 예측변수를 선정하였다. BIOMOD2에서 8개 알고리즘을 적용한 뒤 AUC 0.70 이상의 모형을 이용하여 앙상블 가중평균 모형을 구축하였다. 미래 잠재 서식지는 SSP2–4.5와 SSP5–8.5 시나리오를 중기 미래（2041–2070년）와 후기 미래（2071–2100년）로 구분하여 예측하였다.
현재 서식적합도는 강원 영동과 경북 북부의 일부 산악지역에서 상대적으로 높게 나타났으며, 서부 저지대·평야지역과 남부 해안지역에서는 낮게 나타났다. 앙상블 모형의 ROC 기반 AUC는 0.978로 나타났으며, 이는 본 연구의 자료 구성과 무작위 분할 조건에서 나타난 내부 판별 성능으로 해석하였다. 변수 중요도 분석에서는 자연지역까지의 거리（d_NAT）, 지형습윤지수（TWI）, 해양까지의 거리（d_OCE） 및 강수 계절성（BIO15）이 상대적으로 높은 중요도를 나타냈다.
현재 8,095 km²였던 적합 서식지 면적은 SSP2–4.5 중기 4,271 km², 후기 4,596 km², SSP5–8.5 중기 3,181 km², 후기 4,428 km²로 감소하였다. 두 시나리오 모두 후기 적합면적이 중기보다 일부 증가하였으나, 모든 미래 조건에서 현재보다 작은 면적을 보여 전반적인 감소 경향은 유지되었다. 공간적으로는 서부 및 남부의 저지대·평야지역, 해안 인접 지역과 현재 적합지 주변부에서 상실이 주로 나타난 반면, 강원 영동과 경북 북부의 일부 산악지역에서는 상대적으로 높은 적합성이 유지되었다.
합의 기반 핫스팟 분석 결과, 4개 미래 조건 중 3개 이상에서 공통적으로 확인된 서식지 상실 핫스팟은 3,867 km², 감소강도 핫스팟은 140 km²로 나타났다. 본 연구는 표범장지뱀의 미래 적합성 유지지역과 서식지 감소 위험지역을 국가 규모에서 구분하여 제시하였다는 점에서 의의가 있다. 도출된 결과는 실제 점유분포나 보전 우선지역을 확정하는 자료라기보다, 향후 현장조사, 장기 모니터링 및 보전·관리 대상지역을 검토하기 위한 공간적 기초자료를 제공한다는 점에서 의의를 갖는다.

다국어 초록 (Multilingual Abstract)

This study projected current and future changes in the potential habitat of the Mongolian racerunner （Eremias argus） in the Republic of Korea and identified areas where habitat decline was consistently predicted under multiple future climate conditions. Twenty-three occurrence records collected between 2007 and 2018 and 14 climate, terrain, land-cover, vegetation, and proximity variables were used to construct an ensemble species distribution model in BIOMOD2. Eight algorithms were applied, and models with an area under the receiver operating characteristic curve （AUC） of 0.70 or higher were combined using an ensemble weighted-mean approach.
Future habitat suitability was projected under the SSP2–4.5 and SSP5–8.5 scenarios for 2041–2070 and 2071–2100. The ensemble model produced an AUC of 0.978. Current habitat suitability was relatively high in parts of Gangwon Yeongdong and northern Gyeongsangbuk-do. Distance to natural areas, the topographic wetness index, distance to the ocean, and precipitation seasonality showed relatively high variable importance.
The current suitable habitat area was estimated at 8,095 km². Under SSP2–4.5, it decreased to 4,271 km² in 2041–2070 and 4,596 km² in 2071–2100. Under SSP5–8.5, it decreased to 3,181 km² and 4,428 km², respectively. Although suitable habitat increased slightly in the late-century period compared with the mid-century period, all future estimates remained below the current level.
Consensus-based analysis identified 3,867 km² of habitat-loss hotspots and 140 km² of strong-decline hotspots occurring under at least three of the four future conditions. Habitat loss was mainly projected around the margins of current suitable habitats and in western, central, lowland, and coastal areas, whereas parts of Gangwon Yeongdong and northern Gyeongsangbuk-do retained relatively high suitability.
These results provide national-scale spatial information for identifying areas requiring field surveys, long-term monitoring, and conservation review. The projected areas represent potential environmental suitability rather than confirmed occupancy or population persistence.

목차 (Table of Contents)

• 제1장 서론 1
• 제1절 연구 배경 및 목적 1
• 1. 연구 배경 1
• 2. 연구 필요성 3
• 3. 연구 목적 4
• 4. 소결 6
• 제2장 선행연구 8
• 제1절 대상종 특성 및 관련 연구 동향 8
• 1. 표범장지뱀의 생태·분포 특성 및 보전 필요성 8
• 2. 기후변화에 따른 파충류 서식지 변화와 취약성 10
• 3. 종분포모형의 개념, 적용 및 한계 12
• 4. 기후변화 시나리오 기반 잠재서식지 예측 연구 동향 15
• 5. 선행연구 종합 18
• 제2절 연구의 차별성 20
• 제3장 연구 방법 23
• 제1절 연구 흐름도 23
• 제2절 연구 범위 26
• 1. 공간적 범위 26
• 2. 시간적 범위 27
• 제3절 출현자료 및 환경변수 구축 28
• 1. 출현자료 구축 및 전처리 28
• 2. 환경변수 구축 및 변수선정 30
• 제4절 SDM 모형 구축 및 검증 35
• 1. 개별 종분포모형 구축 35
• 2. 앙상블 모형 구축 38
• 3. 모형 검증 39
• 제5절 미래 기후 시나리오 적용 41
• 제6절 서식지 변화 분석 방법 43
• 1. 이진화 및 등급화 43
• 2. 적합 서식지 면적 변화 분석 44
• 3. 서식지 상실 및 감소강도 핫스팟 분석 46
• 4. 서식지 변화 양상 분석 47
• 제4장 연구 결과 50
• 제1절 서식적합도 및 고적합지역 분석 결과 50
• 1. 현재 서식적합도 및 고적합지역 분석 결과 50
• 2. 모형 성능 및 변수 중요도 분석 결과 51
• 제2절 미래 시나리오별 서식지 변화 분석 54
• 1. 미래 시나리오별 서식지 변화 54
• 2. 서식지 상실 및 감소강도 핫스팟 분석 결과 58
• 3. 서식지 변화 양상 분석 결과 60
• 제5장 고찰 62
• 1. 현재 고적합지역 분포의 생태적 의미 62
• 2. 주요 환경변수의 생태적 의미와 해석상의 한계 63
• 3. 미래 시나리오별 서식지 변화 양상 해석 67
• 4. 서식지 감소 위험지역의 공간적 의미와 취약지역 해석 69
• 5. 미래 적합성 유지지역과 보전·관리를 위한 제언 71
• 6. 연구의 한계와 향후 과제 74
• 제6장 결론 78
• 참 고 문 헌 81
• ABSTRACT 88