Current Issue

지진은 단층에 누적된 응력이 일시적으로 해소되는 단층파열(rupture) 현상의 결과로 발생한다. 내륙 지진과 관련된 단층파열은 주로 취성변형 조건의 상부지각에서 개시 및 전파되는데, 소-중규모 지진은 파열면이 지표까지 확장되지 못하고 지하에서 제한된다. 반면 일반적으로 규모 6 이상의 중-대규모 지진은 지하에서 전파되는 단층파열이 지표까지 확장되어 지표면의 변형을 야기하는데, 이를 지진지표파열(earthquake surface rupture)이라 한다. 고지진(paleoearthquake)에 의한 지표파열 기록은 고지진 발생 이전에 형성된 지형(면)에서 비교적 쉽게 인지할 수 있다. 그러나 지형 변위나 변형의 보존 여부 및 그 정도는 지표파열을 동반하는 지진의 재발주기(지구조적 요소)와 해당 지역의 침식/퇴적률(기후적 요소)에 의해 좌우되며, 도심 개발과 같은 인위적 요소가 영향을 미치기도 한다. 고지진에 수반된 지표파열 기록을 탐지하여 그 분포를 파악하고 변위정보를 분석하는 것은 미래에 재발할 수 있는 지진의 특성을 평가하는데 활용된다(Wallace, 1968; Sieh, 1978; Zielke et al., 2015).

고지진 연구에서는 지표파열을 탐지 및 분석하기 위한 방법으로 인공위성이나 항공기를 이용한 원격탐사 기법을 유용하게 활용해왔다(e.g., Zielke et al., 2015; Johnson et al., 2014; McCalpin and Gath, 2021; Ha et al., 2022). 원격탐사는 관측장치(센서)와 이를 상공에 전달하는 플랫폼으로 구성되며 활용 목적에 따라 센서와 플랫폼의 조합은 다양한데, 고지진 지형 분석에는 주로 입체사진측량 기법(Mikhail et al., 2001)과 라이다(Light Detection And Ranging, LiDAR) 기법(Kim et al., 2014)이 활용된다. 입체사진측량은 자연광원을 활용하는 수동적 센싱 기법으로, 두 장 이상의 2차원 광학사진이 서로 겹쳐지는 중복도를 갖도록 촬영한다. 우리가 사물을 두 눈으로 볼 때 입체감을 느끼듯 동일한 피사체를 서로 다른 각도에서 촬영했을 때 맺히는 상의 차이(시차차; parallax difference)로부터 표고를 계산하는 원리이다. 식생이나 건물과 같은 지표 장애물이 적은 지형조건에서 주로 활용되며, 최근 도심지 개발 이전에 촬영된 과거 항공사진을 이용해서 지형분석을 수행하기도 한다(e.g., Talebian et al., 2016; Kim, T. et al., 2023; Kim and Lee, 2023). 라이다는 인공적으로 레이저 펄스를 발생시키는 능동적 센싱 기법으로, 레이저 펄스의 속도와 반사되어 돌아오는데 걸린 시간 측정을 통해 지표지물까지의 거리를 계산하는 방식이다. 이때 식생 및 건물 등 피복물에서 반사되어 돌아온 데이터는 제거하고 지면에서 반사되어 돌아온 데이터를 보간하면 피복물이 제거된 순수지형을 모델링할 수 있다.

드론 기반의 지형분석은 위성과 항공기에 비해 데이터를 취득할 수 있는 면적은 제한적이지만, 낮은 운영고도에서 느린 속도로 높은 점밀도의 센싱을 통해 1 m 이상 급의 고해상도 지형 데이터를 취득할 수 있다(Oh, 2019). 특히 현장 운용이 상대적으로 용이하며, 결정적으로 위성이나 항공기에 비해 적은 비용으로 데이터를 취득할 수 있기 때문에 최근 고지진 연구에도 드론 활용이 시도되고 있다(e.g., Bemis et al., 2014; Rao et al., 2020; Li et al., 2021; Zou et al., 2022). 본 연구는 지진에 의해 변위된 지형을 대상으로 드론에 의해 취득된 수치표고모델을 활용한 고지진 지표파열 분석을 수행하였다. 연구대상은 현재의 지진지구조환경에서 고지진 지표파열 기록이 보고된 바 있는 한반도 동남부의 양산단층과 몽골 중부의 모고드단층으로 선정하였다. 양산단층 연구는 선행연구(Lee et al., 2023)에 의해 고지진 지표파열이 확인된 경주 삼릉지역의 단층애(fault scarp)에 관한 것으로 식생 제거 효과를 위해 드론 라이다를 활용하였고, 모고드단층 연구는 변위하천에 관한 것으로 건조한 기후로 인해 식생의 간섭이 작아 드론 기반의 입체사진측량 분석을 수행하였다. 연구결과를 토대로 각 지역의 고지진 운동특성을 추론하고 향후 지질재해 원격탐사로의 드론 활용성을 토의하였다.

지표파열은 지진을 일으킨 단층파열이 지표까지 진행되어 단층면이 지표에 노출된 결과, 즉 단층면과 지표면이 교차하는 선형의 지질구조로, 강력한 지진동, 산사태, 침하 등 다른 원인에 의한 지표면 붕괴와 구별된다(그림 1a). 지진지표파열의 가장 일반적인 형태는 선형의 지형학적 낙차를 보이는 단층애로, 이는 단층 양측의 서로 다른 암상에서 비롯되는 차별침식 등 고지진이 아닌 다른 요인에 의해 기반암에 내재된 단층이 인지되는 단층선애(fault-line scarp)와 구별된다(Giano and Schiattarella, 2014). 또한 단층파열면의 미끌림 운동(slip)에 의해 단층 양측 지괴가 서로 반대방향으로 어긋나는 변위가 발생하고, 이러한 변위에 의해 지형이 변형된다. 일반적으로 변위 정보(운동감각, 변위량 등)를 추론할 수 있는 지형을 지형학적 변위지시자(geomorphic offset marker)라 하며, 이는 고지진의 운동감각 또는 단층의 종류에 따라 전형적인 지형 특성을 보이기도 한다. 이 연구에서 다루는 주향이동지진(strike-slip earthquake)은 주로 수평변위를 동반하므로 단층을 가로지르는 선형의 지형이 변위지시자 역할을 한다. 가장 대표적인 예로 변위하천(offset channel)과 하천을 절단하여 배수를 차단하는 능선인 셔터리지(shutter ridge)를 들 수 있으며, 국부적인 수직변위 성분에 의한 새그폰드(sagpond) 또는 프레셔리지(pressure ridge) 등이 있다(그림 1b).

Fig. 1. 
(a) The down-dip rupture width of full- and partial-rupture earthquake along a fault in relation to the extent of seismogenic depth (from Zielke et al., 2015). (b) Geomorphology of strike-slip fault zone (from Burbank and Anderson, 2012). (c) A photo showing co-seismic offset on the surface rupture caused by the 2023 Mw 7.7 Kahramanmaraş, Turkey, earthquake on the Eastern Anatolian Fault (Credit: Prof. Dr. Hasan Sözbilir; Source: https://www.veryansintv.com/prof-dr-sozbilir-en-az-3-fay-segmenti-kirilmis-oldu/). (d) An aerial photo of Wallace Creek and San Andreas Fault (Photograph copyright by David Lynch).

일반적으로 지형변위는 한 번의 지진에 의해 형성되기 보다는 반복적으로 재발하는 지진에 의해 변위가 누적되어 형성된다. 지진 발생과 동시에 변위가 발생하는 것을 지진동시성 지표변위(co-seismic offset)라 하며, 한 차례의 지진에 의해 발생된 변위라는 의미를 내포한다(그림 1c). 지진동시성 변위가 여러 번 발생되어 변위가 누적된 것을 누적변위(cumulative offset)라 한다(그림 1d). 여기서 변위가 누적된다는 것은 지진지표파열이 동일한 단층 구간이나 지점에서 반복되어 발생한다는 의미를 내포하는데, 이는 대부분의 중-대규모 지진이 이미 존재하고 있는 활성단층에서 지진성 거동이 다시 발생한 결과이기 때문이다. 따라서 과거에 여러차례 발생했던 고지진의 흔적을 연구하여 해당 단층(구간)에서의 지진규모를 추정하고 지진재발주기를 분석하는 것이 가까운 미래의 지진 발생 특성을 추론할 수 있는 방법이다(Brune, 1968; Wallace, 1970; Youngs and Coppersmith, 1985; Mason et al., 2006). 한편, 고지진학 접근은 당시의 지진지표파열에 의한 지형변위 흔적이 지금까지 훼손되거나 유실되지 않고 잘 보존되어 있어야 그 연구가 가능한데, 훼손되지 않은 지질학적 기록을 찾는 것은 쉽지 않다. 특히 지형학적 기록은 지형변위가 형성된 후 침식 및 퇴적 작용의 영향을 받아 변형되기 때문에 온전히 보존되기 어렵다.

양산단층과 모고드단층은 모두 판 경계부와 멀리 떨어진 지역에 위치하는 주향이동단층으로, 반복적인 중대형 지진에 의해 지형변위가 누적되어 있다. 그러나, 두 단층을 대상으로 지금까지 수행된 고지진 연구는 주로 지표지질의 층서학적 변형기록을 토대로 고지진 이력을 해석해왔으며, 지형학적 변형기록에 대한 정밀분석은 부족한 상황이다. 한반도 동남부에 위치한 양산단층은 활발한 풍화와 침식작용으로 인하여 지형변위지시자 및 변위 측정 기준점이 보존되기 어렵고, 고속도로를 비롯한 사회 기반 시설물과 울창한 식생 등으로 지형이 피복되어 있어서 지형학적 분석이 상당히 제한된다. 이번 연구의 대상인 양산단층 경주 삼릉구간은 소나무 숲이 발달하기 때문에 지형이 식생으로 완전히 피복되어 있어 광학영상을 활용한 지형분석이 불가능하고, 변위를 측정하기 위한 지형변위지시자를 특정하기 어렵다. 이에 따라 식생 피복물을 제거하고 정밀한 지형분석을 하기 위해 경주 삼릉구간에서는 라이다 기법을 사용하였다. 반면 몽골 모고드단층은 건조한 기후로 인해 풍화와 침식작용이 더디게 진행되며, 더불어 비교적 최근인 1967년 Mw 7.1의 지진으로 지표파열이 발생한 바 있어 지진으로 인해 변위된 지형이 훼손되지 않고 잘 보존되어 있다. 또한 지형 정보를 가리는 식생 및 인공지물 등의 피복물이 없어 라이다 센서의 필요성이 부각되지 않는다. 따라서 모고드단층 사례연구는 데이터 획득과 처리의 용이성 및 경제성을 고려하여 입체사진측량 기법을 사용하였다.

지진의 운동감각은 활성지구조 환경하에서 작동하는 응력 상태와 지각운동 특성을 이해하는데 중요한 정보로, 주로 측지자료, 지진자료 등을 통해 추론 가능하다. 관측 장비가 개발되기 이전에 발생한 고지진의 경우, 단층운동감각과 응력 상태는 고지진에 수반된 지표변형 정보를 통해 추론할 수 있다. 나아가 지형변위지시자가 보이는 누적 변위와 해당 지형변위지시자의 형성 연대를 분석하여 해당 기간 동안의 변위률을 파악할 수 있다. 이는 고지진 굴착조사 정보(지진 발생 횟수, 지진 재발주기, 최후기 운동시기 등)와 더불어 미래지진의 시기, 규모 등을 추론하는 소위 지진재발모델(earthquake cycle model)를 평가하는데 중요하다. 이번 장에서는 두 사례연구를 통해 도출한 고지진 운동특성에 관해 토의하였다.

양산단층 경주 삼릉 구간에 대해 지형단면도에서 인지되는 수직 낙차를 측정한 결과 단층을 따라 변화의 폭이 비교적 큰 수직 낙차를 확인하였다. 고위면의 경우, 비록 모든 측선에서 단층의 동측 지괴가 상대적으로 융기한 것으로 확인되나, SR-C6을 포함한 북쪽은 3 m 이상의 큰 수직 낙차를 보이는 반면 SR-C7 측선부터 남쪽으로 갈수록 변위량이 감소하는 양상을 보인다. 한편 저위면의 경우 비록 수직 낙차는 모두 1 m 미만으로 작지만, 단층의 동측부는 남쪽으로 갈수록 고도가 낮아지는 반면 서측부는 남쪽으로 갈수록 고도가 높아지는 양상을 보인다. 이는 단층을 중심으로 서로 다른 경사방향을 보이는 지형면이 배치되어 있음을 의미하고, 단층의 주향을 따라 상대적으로 융기한 지괴가 달라지는 결과를 초래한다. 단층을 따라 관찰되는 수직 낙차의 다양성은 지진지표파열이 수직변위보다는 주로 수평 변위를 수반한 결과로 해석된다. 고위면과 저위면에서 관찰되는 수직 낙차의 크기 차이는 해당 지형면들이 서로 다른 시기에 형성되었기 때문에 서로 다른 누적 변위를 보이는 것으로 판단된다. 또한 고위면에서 동측지괴가 일관되게 상승한 이유는 수평변위성분보다 그 크기가 작지만 여러 번의 지진에 수반된 수직변위성분이 누적된 결과(Lee et al., 2023)로 판단된다. 한편 저위면의 지형변화가 고지진이 아닌 인위적 요인(왕릉 조성)에 의해 발생했을 가능성을 완전히 배제할 수 없어, 굴착조사와 같은 추가 고지진 연구가 필요하다.

그림 6은 동쪽에서 서쪽으로 경사진 능선과 하천이 반복되는 지형에 대하여 동측 지괴를 고정하고 서측 지괴를 다양한 운동감각의 지진지표파열로 변위 시켰을 때의 수직 낙차 양상을 모식도로 표현한 것이다. 서측 지괴의 지형이 낮아지는 수직변위만 발생한 경우는 단층대를 따라 비교적 일정한 수직변위가 나타나며, 측선 a-a’와 측선 b-b’에서 모두 급경사의 단층애(fault scarp) 지형을 보인다(그림 6a). 서측 지괴의 지형이 높아지는 수직변위만 발생한 경우는 단층대를 따라 일정한 정도의 수직변위가 나타나는 반면, 측선 c-c’와 측선 d-d’에서 모두 경사급변점 전면에 완경사 내지 역경사를 보이는 계단 상의 벤치(bench) 지형을 볼 수 있다(그림 6b). 이와 달리 우수향 변위만 발생한 경우, 측선 e-e’에서는 단층애 지형이 나타나는 반면 측선 f-f’에서는 벤치 지형이 나타나는 등 단층대를 따라 일관되지 않은 정도의 수직변위가 나타난다(그림 6c). 이는 양산단층 경주 삼릉구간에서 측정된 수직변위의 양상을 잘 반영한다. 삼릉 지역의 양산단층은 동측이 융기한 수직변위성분을 보이는데 이는 시추조사에서 관찰된 실트층의 변위와 지형면의 낙차 등을 통해 잘 확인된다(Lee et al., 2023). 반면, 이번 연구에서 제안한 수평변위성분의 경우 변위량과 운동감각을 파악하기 쉽지 않다. 다만 지형학적 변위 분석을 통해 수평변위성분이 수직변위성분에 비해 우세했던 것은 확실해 보인다.

Fig. 6. 
Simplified geomorphic offset model. (a) Model simulating a case where vertical separation occurs as the high-altitude eastern-block uplifts due to vertical displacement movement. (b) Model simulating a case where vertical separation occurs as the low-altitude western-block uplifts due to vertical displacement movement. (c) Model simulating a case where along-fault variations of vertical separation occur due to right-lateral displacement.

국내에서 선행된 고지진 연구에서도 운동감각이 추론된 바 있으나, 주로 지형변위가 아닌 제4기 퇴적층을 변위시킨 단층면상에서 관찰되는 단층조선의 선주각을 이용해 도출되었다. 이러한 이유는 2개 정도로 요약되는데, 활발한 침식과 삭박 작용으로 누적 변위가 지형학적으로 잘 남아있지 않은 조건과 주로 단면노두나 굴착단면 관찰을 통한 고지진 조사가 수행되었다는 점이다. 연구지역인 삼릉구간 또한 시추를 통해 수직변위가 분석된 바 있으나 주향이동 운동감각에 대한 접근은 부족한 실정이다. 최근 지진지표파열에 관한 여러 연구결과는 지표에서의 변형 분할이나 분산 등에 의해 지표파열 복잡성이 비교적 높다는 사실을 잘 보여준다(Klinger et al., 2005; Choi et al., 2018). 가장 대표적인 예로 사교이동지진(oblique-slip earthquake)에 의한 변위의 수평변위성분과 수직변위성분이 지표 및 천부에서 분기된 서로 다른 단층에 구분되어 발생하는 점(Toda et al., 2016; Kim, N. et al., 2023)과 지진지표파열면 상에서 서로 다른 각도나 휘어지는 단층조선이 관찰되는 점(Kearse and Kaneko, 2020; Macklin et al., 2021)을 들 수 있다. 이는 하나의 노두나 굴착조사에서 얻어진 구조 요소와 더불어 특정 단층 구간에 대한 지형학적 접근을 통한 고지진 운동감각 연구가 필요함을 의미한다.

몽골의 경우, 비교적 최근인 1967년 지진으로 지표파열이 발생함은 물론 지형학적 누적 변위가 잘 보존되는 조건으로 인해 수평이동 운동감각이 우수향임을 추론하고 수평변위량(33 m)을 측정하는 것은 어렵지 않다. 여기서는 2차원 단면 분석에서는 동측 지괴가 상승한 반면 3차원 지형복원에서는 서측 지괴가 상승한 것으로 분석된 수직변위에 대해 토의하고자 한다.

몽골 모고드단층에 의해 변위된 하천을 기준으로 누적 수직변위는 약 2 m 로 누적 수평변위의 약 6%에 해당된다. 이는 모고드단층이 주로 거의 수평에 가까운 우수향지진을 반복적으로 겪었음을 지시하며, 현재 변형된 상태의 지형이 주로 주향이동운동에 의한 산물임을 암시한다. 변위하천 및 주변 지형을 3차원 복원한 결과, 서측 지괴(상류)가 단층 운동에 의해 융기하였음을 확인하였으나, 이는 동측 지괴가 더 높은 고도로 나타나는 야외관찰 및 단면도 분석 결과와 다르다. 우리는 하천이 발달하는 지형면이 지진지표파열을 따라 북쪽이 상대적으로 높고 남쪽이 낮은 경사를 보이는 점에 주목하였다. 남쪽으로 경사진 지형면의 우수향 수평 변위를 고려하면, 단층의 양측 지괴가 서로 다른 고도를 보이는 수직 낙차를 발생시킬 수 있다.

그림 7은 고정된 서측 지괴를 기준으로 동측 지괴를 동일한 정도의 변위로 이동시켰을 때 사면의 경사 방향에 따라 수직 낙차가 달라짐을 표현한 것으로, 3개의 경우 모두 동측 지괴가 낮아지는 방향의 수직변위와 우수향 수평 변위를 부여했다. 경사가 거의 없는 평지의 경우, 수직변위와 유사한 정도의 수직 낙차가 발생했으며, 수직 낙차는 동측 지괴가 낮아진 것을 지시한다(그림 7a). 남동에서 북서로 경사한 지형의 경우에는 동측 지괴의 지형이 낮아 보이는 수직 낙차가 형성되었으나, 수직변위보다 큰 정도의 수직 낙차를 보인다(그림 7b). 반면 북서에서 남동으로 경사한 지형의 경우에는 동측 지괴의 지형이 낮아지는 수직변위가 발생하였으나 지표에 형성된 수직 낙차는 동측 지괴의 지형이 높은 방향으로 낙차가 형성되었다(그림 7c). 이는 지형 단면도 수직 낙차 측정 및 LaDiCaoz 분석 결과를 잘 반영하며, 해당 모델은 경사진 지형면이 수평 운동을 겪을 때 지표 지형에서 관찰되는 수직 낙차가 실제 변위 방향과 다른 겉보기 변위(위변위; apparent offset)일 수 있음을 보여준다.

Fig. 7. 
Simplified geomorphic offset model. All models triggered the same degree of vertical and right-lateral displacement based on the fixed western block. (a) Case of a flat surface with no slope. (b) A case of a slope from southeast to northwest, a larger vertical separation occurs than the actual vertical offset. (c) A case of a slope from southeast to northwest, a smaller vertical separation occurred than the actual vertical displacement, and the direction was opposite.

지형변위는 규모 ~ 6 이상 중-대규모 지진의 고유 산물 중 하나로, 이번 연구는 드론에 의해 취득된 수치표고모델을 활용해서 고지진 지형 분석 사례연구를 수행하였다. 드론 기반 데이터는 항공기 및 물리탐사 등 기존에 활용했던 데이터의 분석 결과와 부합하였으며, 특히 양산단층 경주 삼릉 구간에 대해 항공라이다 기반 지형데이터에서 인지되지 않았던 선상지 저위면의 지형 변위가 드론라이다 기반의 지형데이터에서 추론되었다는 점을 주목할 필요가 있다. 경애왕릉이 위치하는 저위면의 수직 낙차에 대해 인위적 정지(整地), 암상 경계부의 차별침식, 형산강 범람 등 다양한 가능성이 존재하기에 수직 낙차의 형성 원인에 대한 추가 조사가 필요하다. 그럼에도 불구하고 양산단층의 자취를 따라 수직 낙차가 분포하기 때문에 고지진에 의해 변위된 결과로 추론할 수 있으며, 수 cm 정확도의 드론 기반 데이터를 통해 인지되지 않았던 변형을 탐지하는 것이 가능했다는 측면에서 그 의미가 있다. 더불어 지형변위 분석 수행 시 수평 운동에 의한 겉보기 수직변위를 주의해야 하며, 운동특성 추론을 위해서는 지형을 3차원적으로 해석하는 것이 중요하다는 점을 시사한다. 특히 비교적 주향이동지진이 우세한 우리나라의 지진지구조환경을 고려할 때, 단층대를 따라 일관되지 않은 수직 낙차가 인지되었을 때 수평이동 감각이 우세한 주향이동단층일 가능성을 고려할 수 있다.

고지진 연구는 단일 데이터만으로는 완전한 해석을 하기 어려우며, 지형분석, 물리탐사, 시추 및 굴착조사, 연대측정 등 다학제적 분석을 통하여 종합적으로 해석하는 것이 중요하다. 양산단층 경주 삼릉구간의 고지진학적 해석은 이번 연구의 정밀 지형 분석과 앞선 연구(Lee et al., 2023)의 시추 결과 등을 종합하여 도출되었다. 이번 연구는 고지진 기록을 추적함에 있어, 고해상도의 자료를 짧은 주기로 취득할 수 있는 드론이 지진 지표변형 분석에서 그 역할을 유용히 수행할 수 있음을 보여준다. 모고드 단층에 대한 해석 또한 이번 연구의 지형분석 결과와 더불어 물리탐사 및 굴착조사 등의 결과가 추가된다면 그 해석의 신뢰도를 더욱 높일 수 있을 것으로 생각된다.