오소리의 적합 서식지 모형과 개체군 유전학에 관한 연구 = A Study on Habitat Suitability Modeling and Population Genetics of Meles amurensis （Mustelidae）

다국어 초록 (Multilingual Abstract)

1. The biology and natural histroy of Meles amurensis inhabiting the Korean Peninsula was comprehensively reviewed in terms of taxonomy, morphology, ecology, genetics, conservation, and management in terms of taxonomy, morphology, ecology, genetics, and conservation and management. Meles amurensis is a medium-sized mustelid distributed broadly across the Northeastern Asia, including the Korean Peninsula. Although Eurasian badgers were long treated as a single species （M. meles）, successive molecular phylogenetic and morphometric studies led to the recognition of M. leucurus and M. anakuma as distinct species. Puzachenko et al. （2025） proposed elevating the Northeastern Asian population of badger to full species rank as M. amurensis, based on cranial morphology demonstrating clear differentiation from M. leucurus. The species is characterized by a stocky body, short limbs, and narrow black facial stripes through the eyes, weighing 4–10 kg. It is distributed throughout the Korean Peninsula including Jeju Island, with core habitats concentrated along the Baekdudaegan mountain range. Reproduction involves embryonic diapause （delayed implantation）, with parturition occurring in late winter to early spring. Behavioral ecology is predominantly nocturnal, with peak activity concentrated after sunset. Diet is broadly omnivorous and opportunistic, varying with local prey availability. Ecologically, the species functions as a mesopredator, seed disperser, and ecosystem engineer. Genetic studies indicated a low level of genetic diversity within the Korean population but confirmed the Korean badger as a distinct phylogeographic lineage. Given the combination of habitat loss, road mortality, historical overharvesting, and the recent reclassification of the taxon, an independent conservation status assessment and species-specific management plan are urgently warranted.

2. Medium and large mammals on the Korean Peninsula are experiencing habitat loss and fragmentation under the compounding pressures of climate change and land-use transformation. This study quantified the current potential habitat distribution of Meles amurensis and identified the key environmental determinants of its habitat suitability using an ensemble SDM framework. Occurrence records collected from 2021 to 2025 through the National Institute of Biological Resources, National Institute of Ecology, and GBIF were integrated with 17 environmental predictor variables （VIF 〈 7） encompassing WorldClim 2.1 bioclimatic variables, topographic indices derived from a digital elevation model （DEM）, the Human Footprint Index （HFI）, forest structural attributes, and soil depth. Using the biomod2 package in R, a weighted mean ensemble model （EMwmean） was constructed by combining eight algorithms—GAM, GLM, MARS, ANN, GBM, RF, CTA, and FDA—and evaluated through 10-fold cross-validation. The ensemble model achieved an AUC of 0.925 and a TSS of 0.678, demonstrating high predictive performance. Among individual algorithms, Random Forest attained the highest AUC （0.935） but exhibited evidence of overfitting, as indicated by a calibration AUC of 1.000. Applying the TSS-optimized threshold of 0.527, the total potential suitable habitat was estimated at approximately 48,700 km², representing 48.53% of the study area.Variable importance analysis consistently identified the Human Footprint Index （HFI; 20.3%） and distance to forest （dist_forest; 18.4%） as the dominant predictors across all algorithms, followed by isothermality （BIO3; 7.4%）, soil depth （6.2%）, and maximum temperature of the warmest month （BIO5; 6.5%）. Response curve analysis revealed a nonlinear threshold response to HFI, with presence probability declining sharply beyond a critical disturbance level, consistent with a nonlinear threshold response. Distance to forest exhibited the highest suitability at 0 m, decreasing steeply with increasing distance from forest edges. These findings indicate that badger distribution is primarily constrained by anthropogenic disturbance avoidance and forest continuity, with climatic stability and burrowing-suitable topographic and edaphic conditions playing complementary roles. The habitat suitability maps generated in this study provide a quantitative basis for prioritizing ecological corridor design and spatially explicit habitat management, while offering a reproducible methodological framework applicable to broader conservation planning for medium-sized mammals facing complex environmental pressures on the Korean Peninsula.

3. This study predicted the spatiotemporal dynamics of potential habitat suitability for Meles amurensis in the Korean Peninsula from the 2020s to the 2100s, using the MaxEnt algorithm in conjunction with five CMIP6 global climate models （GCMs: ACCESS-CM2, BCC-CSM2-MR, GISS-E2-1-G, INM-CM5-0, and MRI-ESM2-0） and four Shared Socioeconomic Pathways （SSP1-2.6, SSP2-4.5, SSP3-7.0, and SSP5-8.5）. Nineteen bioclimatic variables from WorldClim 2.1 and SRTM-derived elevation data were used as environmental predictors. Model performance, evaluated through 15-fold cross-validation, was satisfactory, with a mean training AUC of 0.758 and test AUC of 0.745. Variable contribution analysis identified elevation （Elev） as the dominant determinant of habitat suitability （68.84%）, with climatic variables including the maximum temperature of the warmest month （Bio5）, isothermality （Bio3）, and temperature seasonality （Bio4） playing supplementary roles. Under current climate conditions, potentially suitable habitat covered 20.9% of the total analysis area （30,190 grid cells）, concentrated along the Baekdudaegan mountain range encompassing the Seoraksan, Odaesan, Sobaeksan, and Jirisan massifs. Across all GCMs and scenarios, a consistent near-term expansion of potential habitat was projected for 2021–2040. However, trajectories diverged markedly from mid-century onward: low-emission pathways （SSP1-2.6） maintained habitat area near current levels, while high-emission scenarios （SSP3-7.0 and SSP5-8.5） showed progressive habitat contraction in the latter half of the century. Substantial inter-model uncertainty was evident; INM-CM5-0 projected habitat increases across all scenarios and periods, whereas BCC-CSM2-MR projected a 74.92% reduction under SSP5-8.5 by the end of the century. Gain–Stable–Loss analyses revealed that even scenarios with limited net change involved considerable spatial redistribution of suitable habitat, underscoring the importance of maintaining habitat connectivity. These findings highlight the role of high-elevation Baekdudaegan zones as potential climate refugia and the need to secure elevational movement corridors linking lowland and highland habitats, providing a scientific basis for climate-adaptive conservation strategies for M. amurensis in the face of ongoing climate change.

4. Despite being an important component of temperate forest ecosystems in Northeast Asia, genome-scale population genetic studies of Meles amurensis have been lacking. In this study, double-digest restriction site-associated DNA sequencing （ddRADseq） was applied to characterize genome-wide single nucleotide polymorphism （SNP） variation in Korean populations of M. amurensis. A total of 67 wild individuals were sequenced. Samples were obtained from wildlife rescue centers, endangered species restoration facilities, and national park authorities throughout South Korea. Following library preparation and sequencing, SNP calling was conducted using the Stacks2 pipeline and subsequently filtered using TASSEL 5.0. A final dataset comprising 5,883 high-quality SNP markers from 53 individuals was retained for downstream analyses. Genetic diversity analyses revealed moderate levels of variation, with mean expected heterozygosity （He = 0.243） and observed heterozygosity （HO = 0.254） showing comparable values. The overall mean inbreeding coefficient （FIS = -0.027） indicated a slight excess of heterozygotes and no evidence of substantial inbreeding. Pairwise FST values among regional populations were generally low to moderate, suggesting limited genetic differentiation across most of the Korean Peninsula. Bayesian clustering analysis implemented in STRUCTURE identified K = 3 as the most likely number of genetic clusters. However, genetic differentiation among regions was weak, and no clearly defined geographic population structure was detected. Most individuals exhibited varying degrees of admixture, with particularly pronounced admixture observed in northern Gyeongsangbuk-do. These patterns suggest ongoing gene flow among populations across the Korean Peninsula, while also indicating that geographic features such as the Baekdudaegan mountain range may have partially influenced historical dispersal and connectivity. The present study provides the first genome-wide assessment of genetic diversity and population structure in Korean M. amurensis and establishes an important baseline for future conservation planning, long-term genetic monitoring, and landscape genetic research.

5. Wildlife laundering—the practice of disguising illegally captured wild animals as legally captive-bred individuals through established breeding and trade systems—is increasingly recognized as a major challenge to wildlife conservation and regulatory enforcement worldwide. This study examined institutional vulnerabilities associated with wildlife laundering in South Korea using three representative case studies: the Asiatic black bear （Ursus thibetanus）, wild boar （Sus scrofa）, and Amur badger （Meles amurensis）. Across all three cases, common structural characteristics were identified, including government- supported wildlife farming policies, document-based administrative management systems, and the absence of scientific methods for verifying the origin of individual animals. To evaluate the genetic characteristics of farmed badgers, ddRADseq-based SNP analyses were conducted on a total of 77 individuals, including 10 captive individuals from commercial farms and 67 wild individuals collected throughout South Korea. Following quality filtering, 5,883 high-quality SNP markers from 61 individuals were retained for subsequent analyses. STRUCTURE and principal coordinates analyses （PCoA） showed that captive individuals did not form a genetically distinct cluster separate from wild populations. Likewise, no substantial differences in genetic diversity or overall population structure were detected between captive and wild groups, indicating a high degree of genetic similarity. Under long-term closed breeding conditions, captive populations are generally expected to exhibit measurable genetic differentiation resulting from founder effects, genetic drift, and inbreeding. The absence of such patterns in the present study suggests several non-exclusive explanations, including the relatively recent establishment of the captive population, the use of numerous wild founders during the early development of the farming industry, or continued gene flow between captive and wild populations. Although these findings do not constitute direct evidence of wildlife laundering, they demonstrate that clear genetic boundaries between captive and wild badger populations are currently lacking. The results further suggest that document-based administrative systems alone may be insufficient for objectively verifying the biological origin of captive individuals. The development of SNP-based reference databases, combined with genetic forensic approaches for origin assignment and traceability, would provide a valuable scientific framework for improving wildlife management, monitoring captive breeding systems, and supporting future investigations of illegal wildlife trade and wildlife laundering in South Korea.

6. By integrating the overall findings, this study demonstrated that populations of Meles amurensis on the Korean Peninsula are exposed to multiple threats, including habitat loss and fragmentation driven by climate change, as well as genetic risks associated with illegal capture, wildlife trade, and captive breeding systems. Based on these findings, an integrated conservation and management strategy is proposed. The increasing occurrence of badger in urban areas, particularly in the Seoul metropolitan region, together with associated human injury incidents, has placed the species in the dual role as both a conservation concern and a target of conflict management. We reviewed limitations of the current legal framework under the Wildlife Protection and Management Act and discussed the need for population monitoring, non-lethal mitigation measures, and evidence-based management strategies. SNP marker-based genetic forensic approaches demonstrated the potential for monitoring illegal wildlife trade and tracing population structure, supporting the establishment of a national genetic reference database. Climate change projections identified high-elevation areas along the Baekdudaegan mountain range as potential future climate refugia, emphasizing the importance of maintaining elevational movement corridors and conserving mature forest ecosystems. Given the ecological role of M. amurensis as an ecosystem engineer, conservation outcomes may be further strengthened through the integration of long-term scientific monitoring and citizen science participation. Despite these limitations, this study provides the first comprehensive framework integrating taxonomy, ecology, species distribution modeling, population genetics, wildlife forensics, and conservation management for Meles amurensis in South Korea.

국문 초록 (Abstract)

1. 본 연구는 한반도에 서식하는 오소리（Meles amurensisSchrenck, 1859）의 생물학적 특징을 분류학, 생태학, 유전학 및 보전·관리 측면에서 종합적으로 고찰하였다. Meles amurensis는 식육목 족제비과에 속하는 중형 포유류로 한반도를 포함한 동북아시아에 광범위하게 분포한다. 유라시아 오소리류는 오랜 기간 단일 종（M. meles）으로 취급되어 왔으나, 분자계통학 및 형태 계측학 연구의 축적을 통해 M. meles, M. leucurus, M. anakuma로 순차적으로 분리되었으며, 최근 Puzachenko et al. （2025）은 두개골 형태 변이 분석을 근거로 동북아오소리를 M. leucurus의 아종이 아닌 독립 종 M. amurensis로 재분류하였다. 본 종은 통통한 체형, 짧은 사지, 눈 주변의 좁은 흑색 줄무늬 등 뚜렷한 형태적 특징을 가지며, 체중은 4-10 kg 수준이다. 한반도 전역 및 제주도에 분포하며 백두대간을 중심으로 한 산림 지역이 핵심 서식지이다. 지연 착상 전략을 통해 이듬해 초봄에 출산하고, 야행성 활동 패턴을 유지하면서 일몰 이후 활동이 집중된다. 식성은 광범위한 잡식성으로 지역 먹이 가용성에 따라 유동적이며, 종자 산포, 토양 교란 등 생태계 엔지니어 기능을 수행한다. 유전학적으로는 한반도 개체군이 독자적인 계통지리학적 특성을 보유하나 유전적 다양성은 비교적 낮은 것으로 알려져 있다. 서식지 훼손, 로드킬, 과거의 포획 압력 등으로 인한 지역적 개체군 감소가 우려되며, 최근의 분류학적 재검토 결과를 반영한 보전 상태 재평가와 독립적인 관리 계획 수립이 요구된다.

2. 한반도 중·대형 포유류는 기후변화와 토지 이용 변화가 중첩되는 이중 압력하에 서식지 감소와 파편화를 경험하고 있다. 본 연구는 생태계 엔지니어이자 동북아 온대 산림의 지표종인 Meles amurensis를 대상으로, 앙상블 종 분포 모델링（Ensemble SDM）을 통해 현재 시점의 잠재적 서식지 분포와 이를 결정짓는 핵심 환경 요인을 정량화하였다. 국립생물자원관, 국립생태원 및 GBIF에서 수집한 2021–2025년 출현 기록과 WorldClim 2.1 생물기후 변수, 수치표고모형（DEM） 기반 지형 변수, 인위적 간섭지수（HFI）, 산림 구조 및 토양 관련 변수를 포함한 총 30개 변수 중 다중공선성을 고려하여 선별한 17개 환경변수（VIF 〈 7）를 사용하였다. R의 biomod2 패키지를 활용하여 GAM, GLM, MARS, ANN, GBM, RF, CTA, FDA 8종의 알고리즘을 결합한 TSS 가중 평균 앙상블 모델（EMwmean）을 구축하고, 10회 반복 교차 검증으로 성능을 평가하였다. 앙상블 모델은 AUC 0.925, TSS 0.678을 기록하여 높은 예측력을 확인하였다. 개별 알고리즘 가운데 RF가 AUC 0.935로 가장 우수한 성능을 보였으나, 보정 단계 AUC 1.000에서 확인된 과적합 경향이 지적되었다. TSS 최적 임계치（0.527） 적용 결과, 오소리의 잠재 서식 적합 면적은 약 48,700 km²（전체 연구면적의 48.53%）로 산출되었다. 변수 기여도 분석에서 인위적 간섭지수（20.3%）와 산림까지의 거리（18.4%）가 전체 알고리즘에 걸쳐 일관되게 최상위 예측 변수로 식별되었으며, 등온성（7.4%）, 유효토심（6.2%） 최난월 최고기온（6.5%）이 보조적 기여를 나타냈다. 반응 곡선 분석에서 HFI는 일정 임계치를 초과하는 시점에서 출현 확률이 비선형적으로 급감하는 생태적 문턱값 반응을 보였으며, 산림까지의 거리는 0 m에서 최고 적합도를 보이다가 이격 거리 증가에 따라 급격히 감소하였다. 이상의 결과는 오소리의 분포가 인위적 교란 회피와 산림 연속성에 의해 일차적으로 제한되며, 기후 안정성과 굴 조성에 적합한 지형, 토양 조건이 이를 보완하는 다층적 구조로 형성됨을 시사한다. 본 연구에서 도출된 서식 적합성 지도는 생태 통로 설계 우선순위 설정과 서식지 기반 공간 관리 계획의 과학적 근거로 활용될 수 있으며, 복합적 환경 압력에 노출된 한반도 중·대형 포유류 보전 전략 수립을 위한 재현 가능한 방법론적 틀을 제공한다.

3. MaxEnt 알고리즘과 CMIP6 기반 5개 전 지구 기후 모형（GCM: ACCESS- CM2, BCC-CSM2-MR, GISS-E2-1-G, INM-CM5-0, MRI-ESM2-0） 및 4개 공통 사회경제 경로（SSP1-2.6, SSP2-4.5, SSP3-7.0, SSP5-8.5）를 적용하여 Meles amurensis의 미래 잠재 서식지 변화를 2020년대부터 2100년대까지 시공간적으로 예측하였다. 환경 변수로는 WorldClim 2.1의 19개 생물 기후 변수와 SRTM 기반 고도 데이터를 활용하였으며, 15-fold 교차 검증 결과 훈련 AUC 0.758, 검증 AUC 0.745로 양호한 모형 성능이 확인되었다. 변수 기여도 분석에서 고도가 68.84%로 서식지 적합성의 지배적 결정 인자로 확인되었으며, 최난월 최고기온, 등온성, 기온 계절성 등 기후 변수가 보조적 역할을 수행하였다. 현재 잠재 서식지는 전체 분석 면적의 20.9%（30,190 grid cells）로 나타났으며, 백두대간을 축으로 한 설악산·오대산·소백산·지리산 권역의 주요 산악지대에 집중적으로 분포하였다. 미래 예측 결과, 모든 GCM과 시나리오에서 2021년에서 2040년 사이에는 잠재 서식지가 현재 대비 확대되는 공통된 경향이 나타났다. 그러나 세기 중·후반에는 시나리오별 분기가 뚜렷하여, 저탄소 경로（SSP1-2.6）에서는 현재와 유사한 수준이 유지되는 반면, 고배출 경로（SSP3-7.0, SSP5-8.5）에서는 서식지 면적이 감소하는 경향이 나타났다. 다만 GCM 간 예측 편차가 매우 커, INM-CM5-0은 전 기간 · 전 시나리오에서 현재 대비 서식지 증가를 예측한 반면, BCC-CSM2-MR은 SSP5-8.5 세기말 기준 74.92% 감소를 예측하였다. Gain·Stable·Loss 분석에서는 총량 변화가 제한적인 시나리오에서도 기존 서식지의 공간적 재배치가 확인되어, 서식지 연결성 유지의 중요성이 부각되었다. 분석 결과는 백두대간 고산 지대의 기후 피난처 기능과 저지대-고지대 간 수직 이동 통로 확보의 필요성을 시사하며, 기후 위기 대응형 오소리 서식지 보전 전략의 과학적 근거를 제공한다.

4. Meles amurensis는 동북아시아 온대림 생태계의 핵심 구성 종임에도 불구하고 유전체 수준의 대량의 분자 마커를 활용한 집단유전학적 연구는 수행된 바 없었다. 본 연구에서는 ddRADseq 기법을 적용하여 전국 야생동물구조센터, 멸종위기종복원센터, 국립공원공단 등을 통해 수집한 야생 오소리 67개체의 유전체 라이브러리를 시퀀싱 하였다. Stacks 2 파이프라인과 TASSEL 5.0을 이용한 SNP 스크리닝 결과 최종 53개체로부터 5,883개의 SNP 마커를 확보하였다. 유전적 다양성 분석 결과 평균 기대 이형접합도（He = 0.243）와 관찰 이형접합도（Ho = 0.254）는 유사한 수준을 보였으며, 전체 평균 근친교배 계수（FIS = -0.027）는 0에 근접한 음의 값을 나타내어 개체군이 Hardy-Weinberg 평형에 가까우면서도 약한 이형접합 과잉 경향을 보이는 것으로 해석되었다. Pairwise FST 값은 대부분의 개체군 조합에서 0.03-0.18 범위로 낮은 수준의 분화를 보였다. STRUCTURE 분석에서 최적 군집 수는 K = 3임에도 불구하고, 지역별 뚜렷한 유전적 구조는 확인되지 않았다. 대부분의 개체는 admixture 패턴을 나타냈으며, 특히 경북 북부（GB_B） 개체군에서 이러한 경향이 두드러졌다. 이는 오소리 개체군 간 유전자 흐름이 지속적으로 이루어지고 있음을 시사한다. 다만 일부 지역에서 관찰된 유전적 조성의 차이는 백두대간과 같은 지형적 요인이 개체군 간 유전적 교환에 부분적인 영향을 미쳤을 가능성을 제기한다.
5. 야생동물 세탁（wildlife laundering）은 불법으로 포획된 야생 개체를 합법적인 사육·번식 체계를 통해 합법 개체로 위장하여 유통하는 행위로, 전 세계적으로 야생동물 보전과 관리 체계의 신뢰성을 저해하는 주요 문제로 인식되고 있다. 오소리（Meles amurensis） 사례를 중심으로 국내 야생동물 세탁의 구조적 특성과 제도적 취약성을 분석하였다. 분석 결과, 국가 주도의 사육 장려 정책, 서류 중심의 행정 관리 체계, 그리고 과학적 기원 검증 수단의 부재가 결합하여 야생 개체가 사육 개체로 편입될 수 있는 구조적 취약성을 공유하는 것으로 나타났다. 특히 국내 오소리 사육 산업의 유전적 특성을 평가하기 위하여 농장 오소리 10개체와 야생 오소리 67개체를 포함한 총 77개체를 대상으로 ddRADseq 기반 SNP 분석을 수행하였다. 분석 결과 총 5,883개의 SNP가 확보되었으며 STRUCTURE 및 PCoA 분석에서 사육 개체군은 야생 개체군과 명확하게 구분되는 독립적인 유전 집단을 형성하지 않았다. 또한 사육 집단과 야생 집단 간 유전적 다양성 및 유전 구조에서도 뚜렷한 차이가 확인되지 않아 두 집단 간 높은 유전적 유사성이 나타났다. 이러한 결과는 국내 사육 오소리 집단이 장기간 폐쇄적으로 유지된 집단이라기보다 비교적 최근에 형성되었거나, 과거 다수의 야생 개체가 창시 집단으로 활용되었을 가능성을 시사한다. 또한 관찰된 낮은 유전적 분화는 초기 야생 개체의 대규모 도입, 짧은 격리 기간, 또는 지속적인 유전자 흐름 등 다양한 요인으로 설명될 수 있으며 사육 집단과 야생 집단 간 경계가 명확하지 않음을 보여준다. 본 연구는 국내 야생동물 관리 체계가 문서 기반의 행정 관리만으로는 개체의 실제 기원을 검증하는 데 한계가 있음을 제시하며, 향후 사육 개체 이력 관리의 신뢰성 확보를 위해 SNP 기반 유전자 프로파일링 데이터베이스 구축과 유전적 법과학을 활용한 과학적 검증 체계 도입이 필요함을 시사한다.

6. 전체 연구 결과를 종합하여, 본 연구는 한반도 내 Meles amurensis 개체군이 기후변화에 따른 서식지 소멸과 불법 포획·유통에 의한 유전적 오염이라는 복합적 위협에 처해 있음을 밝히고, 이에 기반한 통합적 보전 및 관리 전략을 제언하였다. 수도권을 중심으로 도심 출몰과 부상 피해 사례가 잇따르면서 오소리는 보전 대상인 동시에 갈등 관리 대상 종이라는 이중적 지위에 놓이게 되었다. 현행 「야생생물 보호 및 관리에 관한 법률」 체계의 한계를 검토하고, 유해야생동물 지정 논의, 비치사적 관리 기법 도입, 개체군 현황 조사의 필요성을 제시하였다. SNP 마커 기반 유전적 법과학 접근을 통해 불법 포획·유통 모니터링 및 개체군 구조 추적이 가능함을 확인하였으며 전국 단위 유전정보 데이터베이스 구축을 제언하였다. 기후 변화 시나리오 분석 결과, 백두대간 고산 지대가 미래 기후 피난처로 기능할 것으로 예측되었으며, 수직적 이동 통로 확보 및 노령림 보전을 통한 생태 네트워크 강화가 요구된다. 아울러 오소리의 생태계 엔지니어로서의 기능적 가치를 고려할 때 과학적 모니터링과 시민 과학 참여를 결합한 사회적 인식 제고 전략이 장기 보전 성과에 기여할 것으로 판단된다. 본 연구는 출현 자료의 공간적 편향 유전 분석 표본의 지리적 제한, 경관 유전학적 통합 분석의 부재 등 방법론적 한계를 내포하며, 이를 보완하는 후속 연구가 요구된다.

목차 (Table of Contents)

• Ⅰ. 오소리의 생물학적 특징 1
• 1. 분류학적 고찰 1
• 2. 형태적 특징 3
• 3. 분포 7
• 4. 화석 기록 10
• 5. 번식 11
• 6. 생태학적 특징 12
• 1） 서식지 12
• 2） 행동 13
• 3） 먹이 15
• 4） 질병과 기생충 17
• 7. 유전학적 연구 18
• 8. 오소리의 분류학적 재검토 19
• 9. 보전 및 관리 21
• Ⅱ. 오소리의 현재 분포와 서식지 모형 23
• 1. 서론 23
• 2. 연구 방법 27
• 1） 분석 개요 27
• 2） 출현 자료 27
• 3） 환경 변수 구성 28
• （1） 환경 변수 구축 및 추출 28
• 4） 앙상블 종분포모형 구축 30
• （1） 개별 알고리즘 선정 및 적합 30
• 3. 결과 32
• 1） 변수 스크리닝 및 후보 모델 선택 32
• 2） 개별 모델 평가 결과（Model Evaluation） 34
• 3） 모델링 분석 결과 및 변수 기여도 36
• （1） ANN （Artificial Neural Network） 36
• （2） CTA （Classification Tree Analysis） 42
• （3） FDA （Flexible Discriminant Analysis） 48
• （4） GAM （Generalized Additive Model） 54
• （5） GBM （Generalized Boosted Models） 60
• （6） GLM （Generalized Linear Model） 66
• （7） MARS （Multivariate Adaptive Regression Splines） 72
• （8） RF （Random Forest） 78
• 4） 앙상블（Ensemble） 모델 구축 및 분석 84
• （1） 모델 선정 및 앙상블 구축 과정 84
• （2） 모델의 타당성 및 신뢰도 검증 및 결합 방식 비교 86
• （3） 앙상블 참여 알고리즘별 변수 기여도 분석 89
• （4） 앙상블 모델 서식 적합성 분석 결과 91
• （5） 앙상블 모델 환경 변수 기여도 분석 94
• （6） 앙상블 모델과 개별 알고리즘의 비교 분석 96
• 4. 고찰 102
• 1） 인위적 간섭과 공간적 파편화의 결정적 제약 102
• 2） 지형적 복잡성과 은신처 확보의 생태적 기작 103
• 3） 기후적 안정성과 생리적 적응성 104
• 4） 앙상블 모델의 방법론적 타당성 및 보전적 의의 106
• 5） 결론: 보전 전략의 패러다임 전환 107
• Ⅲ. 오소리의 미래분포 예측 110
• 1. 서론 110
• 2. 연구 방법 112
• 1） 기후 시나리오 및 전 지구 기후 모형 112
• 2） 환경 변수 선정 및 데이터 구축 114
• 3） MaxEnt 모형 설정 및 최적화 116
• 4） 변화 분석 117
• 5） 변수 중요도 및 생태적 해석 117
• 3. 결과 118
• 1） MaxEnt 모델 성능 평가 및 변수 기여도 검증 118
• 2） 기후변화 시나리오별 미래 서식지 공간 투영 결과 120
• 3） 미래 시나리오 분석 결과 122
• （1） ACCESS-CM2 모델 기반 미래 서식지 변화 예측 122
• （2） BCC-CSM2-MR 모델 기반 미래 서식지 변화 예측 127
• （3） GISS-E2-1-G 모델 기반 미래 서식지 변화 예측 132
• （4） INM-CM5-0 모델 기반 미래 서식지 변화 예측 137
• （5） MRI-ESM2-0 모델 기반 미래 서식지 변화 예측 142
• 4. 고찰 147
• 1） MaxEnt 모형 성능의 의의와 구조적 제약 147
• 2） 핵심 환경 변수의 생태적 해석 148
• （1） 고도의 지배적 중요도 148
• （2） 기후 변수의 보조적 역할 148
• 3） 시나리오별 서식지 변동 양상의 해석 149
• （1） 세기 초반 서식지 확장의 의미 149
• （2） 탄소 배출 경로와 서식지 변동의 방향성 149
• （3） 기후모델 간 차이 149
• 4） 서식지 변동의 경관 생태적 함의 151
• （1） 서식지 총량 변화와 내부 전환 구조의 분리 해석 151
• （2） 고배출 시나리오 하에서의 서식지 압착 151
• （3） 기후 피난처로서의 고산 산림 151
• 5） 관리 함의 152
• （1） 오소리의 생태적 기능과 관리적 고려 152
• （2） 기후 변화 대응형 서식지 관리 방향 152
• （3） 서식지 확장기의 갈등 관리 152
• 6） 연구 한계 및 향후 연구 153
• Ⅳ. ddRADseq을 활용한 오소리의 개체군 유전학 154
• 1. 서론 154
• 2. 연구방법 158
• 1） 시료 수집 및 DNA 추출 158
• 2） ddRADseq 라이브러리 제작 및 시퀀싱 160
• 3） SNP 분리 및 스크리닝 161
• 4） 유전적 다양성 및 개체군 구조 분석 162
• 3. 결과 163
• 1） 유전적 다양성 163
• 2） 개체군 간 유전적 분화도 165
• 3） 개체군 유전 구조 167
• 4. 고찰 172
• 1） 유전적 다양성 및 근친교배 양상 172
• 2） 개체군 유전 구조 및 연결성 174
• 3） 관리 전략 수립을 위한 제언 175
• 4） 연구 한계 및 향후 연구 176
• Ⅴ. 우리나라 야생동물과 유전체 연구를 활용한 오소리 사례 177
• 1. 야생동물 세탁의 역사적 배경 177
• 1） 야생동물 세탁의 정의 177
• 2） 국제 사회의 대응 178
• 3） 다른 나라의 야생동물 세탁의 사례 179
• 2. 우리나라 야생동물 세탁 181
• 1） 반달가슴곰（Ursus thibetanus） 181
• 2） 멧돼지（Sus scrofa） 184
• 3. 오소리（Meles amurensis） 사례 연구 186
• 1） 연구방법 188
• 2） 결과 190
• （1） 개체군 유전 구조 190
• 3） 오소리 사례에 대한 고찰 195
• 4. 고찰 197
• Ⅵ. 오소리의 보전 및 관리 200
• 1. 도심 출몰과 인명 피해: 유해야생동물 지정 논의 200
• 2. 법적 지위 체계 및 관리 제도 정비 202
• 3. 유전적 법과학 기반 모니터링 구축 203
• 4. 기후 적응형 서식지 관리 및 생태 네트워크 강화 205
• 1） 미래 기후 피난처의 예비 보호구역 지정 205
• 2） 생태 통로의 전략적 배치 205
• 3） 노령림 보전을 통한 미세 기후 안정성 확보 206
• 4） 인간-야생동물 갈등 완화 206
• 5. 생태계 엔지니어 가치 확산 및 사회적 인식 제고 207
• 6. 연구의 한계와 향후 과제 208
• Ⅶ. 종합 제언 209
• 참고문헌 210
• 영문초록 284