본 연구에서는 콕스 비례위험모형을 이용하여 산불피해 소나무의 고사에 영향을 미치는 인자를 밝히고자 하였다. 지표화 피해 소나무를 대상으로 고사 영향 인자를 조사하고 산불 발생 7년 차까지 고사 발생 모니터링을 수행하였다. 수집된 자료를 기반으로 생존분석을 실시하였다. 분석 결과, 고사 위험성을 증가시키는 변수는 dNDVI(delta Normalized Difference Vegetation Index), dNBR(delta Normalized Burn Ratio), 경사, 나무에 남겨진 그을음의 상대적인 면적과 평균적인 높이를 나타내는 수피 그을음 지수(Bark Scorch Index, BSI)와 수피 그을음 높이(Bark Scorch Height, BSH)로 나타난 반면, 음의 관계를 가지는 변수는 고도, 흉고직경, 수관층 수분스트레스 변화를 나타내는 수분스트레스지수(dleta Moisture Stress Index, dMSI)로 나타났다(p<0.001). 콕스 비례위험모형의 유의성을 확인하기 위한 변수별 비례위험가정 검증에서는 사면방향을 제외한 모든 인자가 모형에 적합하며 고사 발생에 유의미한 영향을 미치는 것으로 나타났다. 생존 곡선 분석에서 가장 큰 생존율 차이를 보인 변수는 BSI였으며(p<0.0001), 원격탐사를 통해 얻어진 환경변화 인자들(dNDVI, dNBR, dMSI) 역시 큰 생존율 차이를 나타내었다(p<0.0001). 이러한 결과는 산불 이후 소나무의 잠재적인 고사 위험성을 고려한 복원계획 수립을 위한 기초자료로 활용될 것으로 기대된다.

In this study, we aimed to identify the factors influencing post-fire mortality in Korean red pine (Pinus densiflora) using Cox’s proportional hazards model and analyze the impact of these factors. We monitored the mortality rate of fire-damaged pine trees for seven years after a forest fire. Our survival analysis revealed that the risk of mortality increased with higher values of the delta normalized difference vegetation index (dNDVI), delat normalized burn ratio (dNBR), bark scorch index (BSI), bark scorch height (BSH) and slope. Conversely, the risk of mortality decreased with higher elevation, greater diameter at breast height (DBH), and higher value of delta moisture stress index (dMSI) (p < 0.01). Verification of the proportional hazards assumption for each variable showed that all factors, except slope aspect, were suitable for the model and significantly influenced fire occurrence. Among the variables, BSI caused the greatest change in the survival curves (p < 0.0001). The environmental change factors determined through remote sensing also significantly influenced the survival rates (p < 0.0001). These results will be useful in establishing restoration plans considering the potential mortality risk of Korean red pine after a forest fire.