Detection of Surface Changes by the 6th North Korea Nuclear Test Using High-resolution Satellite Imagery

고해상도 위성영상을 활용한 북한 6차 핵실험 이후 지표변화 관측

Abstract

On September 3rd 2017, strong artificial seismic signals from North Korea were detected in KMA (Korea Meteorological Administration) seismic network. The location of the epicenter was estimated to be Punggye-ri nuclear test site and it was the most powerful to date. The event was not studied well due to accessibility and geodetic measurements. Therefore, we used remote sensing data to analyze surface changes around Mt. Mantap area. First of all, we tried to detect surface deformation using InSAR method with Advanced Land Observation Satellite-2 (ALOS-2). Even though ALOS-2 data used L-band long wavelength, it was not working well for this particular case because of decorrelation on interferogram. The main reason would be large deformation near the Mt. Mantap area. To overcome this limitation of decorrelation, we applied offset tracking method to measure deformation. However, this method is affected by window kernel size. So we applied various window sizes from 32 to 224 in 16 steps. We could retrieve 2D surface deformation of about 3 m in maximum in the west side of Mt. Mantap. Second, we used Pleiadas-A/B high resolution satellite optical images which were acquired before and after the 6th nuclear test. We detected widespread surface damage around the top of Mt. Mantap such as landslide and suspected collapse area. This phenomenon may be caused by a very strong underground nuclear explosion test. High-resolution satellite images could be used to analyze non-accessible area.

2017년 9월 3일 북한에서 발생한 인공지진 신호가 기상청 지진관측망에 관측되었다. 진앙은 풍계리 핵실험 지역으로 지금까지의 실험 중 가장 강력한 실험이었다. 직접적 접근이 제한되는 상황에서 6차 핵실험에 의한 주변 지표변화 연구를 위해 고해상도 위성 자료를 활용하였다. 우선, 지표변위 관측을 위해 ALOS-2 위성 자료를 활용하여 레이더 간섭기법(InSAR: SAR Interferometry)을 적용하였다. 하지만 6차 핵실험 주변 지역의 대규모 지표변위와 강한 진동으로 인해 낮은 긴밀도(coherence) 값을 지니며 레이더 간섭도가 생성되지 않았다. 이는 강한 진동으로 인한 표면의 변화와 레이더 간섭도가 생성 가능한 최대 지표변위 관측 범위보다 큰 변위가 발생했기 때문으로 추정되며 이러한 한계를 극복하기 위해 오프셋 트래킹(Offset Tracking) 방법을 활용하였다. 오프셋 트래킹 방법은 6차 핵실험 전 후 위성 영상레이더의 강도 영상(Intensity)에 대한 교차상관기법(Cross-Correlation)을 이용하는 것으로 상관관계 추정을 위해 사용된 윈도우 크기에 따라 결과가 달라지는 단점이 존재한다. 본 연구에서는 32부터 224까지 16단계로 윈도우 크기를 변화시키고 그 결과의 통계적 처리 후 지표변위를 생성하였다. 그 결과, 6차 핵실험 장소를 기준으로 만탑산 서쪽 지역에서 최대 3 m의 지표변위를 관측하였다. 또한, 고해상도 광학 영상을 활용하여 6차 핵실험에 의한 산사태 및 함몰지역으로 추정되는 지역을 확인하였다. 이러한 현상은 매우 강력한 지하 핵실험에 의한 것으로 판단되며 기존 음파 및 지진파를 이용한 핵실험 분석뿐만 아니라 고해상도 위성영상을 활용하여 비접근 지역에 대한 보조 분석 자료로 활용이 가능 할 것이다.

1. 서론

핵실험은 강력한 인공지진을 발생시키기 때문에 지진파 또는 공중음파 관측 등의 방법을 이용하여 핵실험을 관측하여 왔다(Che et al., 2014). 핵실험은 강력한 진동으로 인해 지표의 변화를 발생시키기 때문에 이에 대한 연구가 필요하지만 대부분의 경우 접근이 제한되기 때문에 직접 접근을 통한 연구가 부족한 실정이다. 이와 같이 접근이 어려운 원격지에 대해 고해상도 위성자료를 이용한 연구가 이루어지고 있다. 특히 위성 영상레이더(SAR: Synthetic Aperture Radar) 자료로부터 레이더 간섭기법(InSAR: SAR Interferometry)을 이용하여 미국 네바다 주, 중국 등 핵실험에 의해 발생한 지표변위가 관측되기도 하였다(Vincent et al., 2003; Vincent et al.,2013).

북한은 과거 2006년부터 풍계리 지역에서 지속적인 핵실험을 수행하였으며 2016년 1월에 수행된 4차 핵실험의 경우 ALOS-2 (Advanced Land Observing Satellite 2) 자료를 이용하여 10 cm의 지표변화가 관측되었다(Wei,2017). 특히, 2017년 9월 3일에 수행한 6차 핵실험은 과거 핵실험에 비해 수 배 이상 강력한 에너지가 발생하였으며 X-band TerraSAR-X 위성 자료를 이용하여 최대 3.5 m의 지표변위가 관측되었다(Wang et al., 2018). Wanget al.(2018) 연구에서는 위상 비간섭(decorrelation)으로 인해 레이더 간섭기법으로 지표변위 관측이 어렵기 때문에 오프셋 트래킹 기법을 활용하여 지표변화를 관측하였다. 오프셋 트래킹 기법은 레이더 간섭기법과 달리 이벤트 전·후 위성 영상레이더의 강도영상 사이의 교차 상관기법을 이용하는 방법으로 상관관계를 계산할 때 사용되는 윈도우 크기에 따라 결과가 달라지기 때문에 이에 대한 고려가 필요하다. 따라서 본 연구에서는 파장 길이가 긴 L-band ALOS-2 위성 자료를 이용하여 6차 핵실험에 의한 지표변위 산출 연구를 진행하였다. Lband는 X-band에 비해 파장이 길며 비간섭 효과가 작은 장점이 있기 때문에 이를 통한 레이더 간섭기법을 사용하였으나 Wang et al.(2018) 연구와 같이 비간섭으로 인한 낮은 긴밀도(coherence) 때문에 간섭도가 생성되지 않았다. 따라서 간섭도 생성 조건을 바탕으로 긴밀도가 낮은 이유를 추정하고 오프셋 트래킹의 단점을 극복하기 위해 윈도우 크기의 변화를 통해 최종 지표변위를 산출하였다. 이와 더불어 위성 영상레이더뿐만 아니라 고해상도 광학영상을 활용하여 시각적 변화 분석을 진행하고 핵실험 이후 지표변화를 종합적으로 분석하고자 하였다.

2. 연구지역 및 사용된 데이터

풍계리 핵실험장 지역은 북한 함경북도 길주군에 위치하고 있으며 만탑산 주변에 갱도를 만들고 2006년부터 2017년 9월 3일까지 총 6차례의 핵실험을 진행하였다(Fig. 1). 만탑산은 고도 2000 m 성층화산으로 아래부분의 기반암은 단단한 화강암으로 되어 있으나 표면은 약 200 m 두께의 화산퇴적물이 쌓여있다고 알려져 있다(Pavian and Coblentz, 2017). 이와 같은 지질학적 조건으로 인해 강한 진동에도 기반이 무너지지 않기 때문에 핵실험에 의한 방사능 유출 위험이 적은 곳으로 알려져 있으며, 이 곳에서 총 여섯 차례의 핵실험을 진행되었다.

Fig. 1. Study area on hill shade map with satellite coverage. Black square in the upper left inset represented Punggye-ri area. Dotted line represents SAR image coverage and solid line is optical image coverage. Symbol with number indicates the location of nuclear test in right figure.

6차 핵실험은 북한 핵실험 중 가장 강력한 에너지가 방출되었으며 2016년 9월 9일에 수행된 5차에 비해 TNT (Trinitrotoluene: 강력 폭약) 기준으로 환산한 값(수율)이 약 9배 이상 큰 것으로 추정되기도 하였다(Stevensand O’Brien, 2018). 강력한 핵 실험에 의해 많은 지표변화가 생겼을 것으로 추측이 되며 직접적 접근이 제한되는 상황에서 지표변화 관측을 위해 일본의 ALOS-2 영상레이더와 프랑스 Pleiades 고해상도 위성 광학영상을 활용하였다(Table 1). 영상레이더와 광학영상 모두 태양 동기궤도로서 정지궤도와 같이 지속적인 촬영은 불가능하지만 고해상도 영상을 획득할 수 있기 때문에 핵실험과 같이 국지적 변화 탐지가 가능하다. 따라서 6차 핵실험 발생 날짜인 2017년 9월 3일을 포함하는 쌍으로 데이터를 수집하여 지표변화를 관측하였다. 두 영상 모두 2 m 내의 고해상도이며 영상레이더는 마이크로파 대역 주파수를 사용하며 측면(Side-Looking)으로 관측한다. 이와 비교하여 광학영상은 가시광선 대역을 사용하며 거의 수직(Nadir)에 가깝게 관측하기 때문에 기하왜곡이 적은 편이며 육안으로 직접적인 변화탐지 활용이 가능하다. 영상레이더의 경우 기하학적 왜곡이 심하며 육안 분석이 쉽지 않지만 이벤트 전·후의 위상값 차이를 이용하여 지진, 화산, 산사태 등으로 인한 지표변위 관측에 많이 활용되고 있는 자료이다(Lee et al., 2017; Lee et al., 2018). 본 연구에서는 서로 다른 특성을 지니는 두 종류의 영상을 이용하여 6차 핵실험 전·후의 지표 변화 탐지를 수행하였다.

Table 1. Characteristics of SAR and Optical imagery

* Ascending track, Right-looking

3. 연 방법 및 결과

1) 연구방법

(1) 위성레이더 긴밀도를 이용한 지표 변화 영역 추정

Wei (2017)는 4차 핵실험에 의한 지표변위 관측을 위해 레이더 간섭기법을 사용하였으며 성공적으로 지표 변위를 관측하였다. 레이더 간섭기법은 수 cm 단위의 정밀 지표변위 관측이 가능하기 때문에 비록 핵실험이 지하 수 백 미터에서 수행되어도 이를 활용하여 핵실험에 의한 지표 변화 관측이 가능하다. 따라서 본 연구에서는 6차 핵실험에 의한 지표변화 관측을 위해 레이더 간섭기법을 적용하였으나 레이더 간섭도 생성이 되지 않았다. 레이더 간섭도를 이용한 정확한 지표변위 산정을 위해서는 긴밀도(Coherence) 값이 중요하다. 사용된ALOS-2 영상은 L-band로 파장이 길기 때문에 X-, C-band와 같은 짧은 파장에 비해 긴밀도가 높게 유지되는 특징이 있음에도 불구하고 핵실험이 수행된 지역 주변에 매우 낮은 긴밀도가 관측되었다(Fig. 2(a)). 일반적으로 긴밀도 값이 0.3 이하인 경우 지표변위 산정의 신뢰도가 매우 떨어지기 때문에 정확한 지표변위 산출이 불가능하다(Baran et al., 2005).

긴밀도에 영향을 미치는 다양한 요소가 있지만 서로 다른 시기에 촬영된 위성의 거리차이에 의한 베이스라인 비상관(Baseline decorrelation), 시간이 지남에 따른 식생 또는 토양의 지표 변화로 인한 시간 비상관(temporaldecorrelation) 등이 크게 영향을 준다(Lee et al., 2015a). Fig. 2(a)는 6차 핵실험 기간을 포함하여 14일의 관측 시기 차이를 가진 긴밀도 지도이며, Fig. 2(b)는 6차 핵실험 이후의 14일 관측 시기 차이를 가진 긴밀도 지도로써 각각 -11 m, -9 m의 매우 짧은 수직베이스라인 거리를 지닌다. 두 그림에서 관측 기간 차이는 동일하며 수직베이스라인도 무시할 수 있을 정도의 차이이기 때문에 6차 핵실험에 의해 발생한 강한 진동으로 토양의 지표가 매우 광범위하게 변화하고 그 변화량도 매우 커져 긴밀도가 낮아진 것으로 예상된다. 긴밀도가 낮은 영역을 기반으로 핵실험에 의한 지표변위 영향 범위는 약 9.2 km²인 것으로 추정된다.

Fig. 2. Coherence map comparison near the Mt. Mantap area. (a) Coherence map of 20170829-20170912 pair (including 6th nuclear test date), (b) coherence map of 20170912-20170926pair (after 6th nuclear test date).

6차 핵실험 전·후 영상레이더에서 레이더 간섭도 생성을 위해서는 최대 위상변화량 이내의 값이 요구되며 최대 위상변화값은 다음과 같다(Hanssen, 2001).

\(f=\frac{B w}{R S F}\)       (1)

여기서, f는 최대 위상변화값(critical phase gradient:cycle/pixel)이며, Bw는 bandwidth, RSF는 Range Sampling Frequency를 의미한다.

최대 위상변화값을 변위값으로 변경하면 영상레이더에서 관측이 가능한 최대 지표변위량은 아래와 같다.

\(\mathrm{d}=\frac{f \times \lambda}{2 \times 1000}\)       (2)

여기서, d는 최대 지표변위 관측범위(mm/pixel) 이며, λ는 파장이다.

(1) 식을 이용하여 ALOS-2 위성 영상레이더 파라미터를 계산하면 간섭도가 생성될 수 있는 최대 위상변화 값은 0.7577(cycle/pixel)이며 (2)식에 의해 최대 지표변위량은 약 89 mm/pixel로 계산된다. 최대 지표변위량은 경사가 급하지 않는 지형을 가정한 상태이나 핵실험이 진행된 만탑산의 경우 산악지형이기 때문에 그 값은 더욱 감소하게 된다. 특히, 핵실험에 의한 영향 범위가 반경 2 km이내의 비교적 좁은 지역이지만 그 변위량은 최대 3.5 m 이상의 지표변위가 발생한 것으로 알려져 있다(Wang et al., 2018). 이와 같은 결과를 통해 6차 핵실험으로 인한 지표변위는 레이더 간섭기법을 활용하여 관측할 수 있는 최대 지표변위량보다 큰 변위가 생겼고 대규모의 지표변위와 강한 진동으로 인한 긴밀도 저하로 인해 기존의 레이더 간섭기법을 적용할 수가 없음을 알 수 있다.

(2) Multi-Kernel Offset Tracking을 이용한 영상레이더 처리 및 광학위성 영상 처리

6차 핵실험과 같이 변위량이 큰 경우 위성레이더 간섭기법은 비상관(decorrelation)으로 인해 긴밀도 생성이 불가능함을 확인하였다. 산사태, 빙하의 이동 관측 등 지표 변위가 큰 경우 레이더 간섭도 생성이 힘든 단점을 극복하기 위해 오프셋 트래킹 방법을 많이 사용한다(Chae et al., 2017). 오프셋 트래킹은 영상레이더의 강도 영상으로부터 교차상관기법을 사용하여 시선방향(LOS:Line Of Sight)과 비행방향(Along Track)의 2차원 변위를 관측할 수 있다. 하지만 기존의 오프셋 트래킹 방법은 변위 추정 윈도우 크기에 따라 정밀도가 떨어지는 단점이 있다(Jung et al., 2013). 이를 극복하기 위해 다중변위 커널(Multi-Kernel)을 활용한 오프셋 트래킹 방법이 제안되었다(Chae et al., 2017).

Fig. 3(a)~(d)는 시선 방향의 결과이며 Fig. 3(e)~()는비행 방향의 윈도우 크기에 따른 오프셋 트래킹 결과의 지표변위 계산값이다.

Fig. 3. Pixel Tracking Results, (a)~(d) are range deformation, (e)~(h) ae azimuth deformation. Window sizes are various: (a), (e): 96, (b), (f): 128, (c), (g):160, (d), (h): 256.

윈도우 크기가 작을 때는 세밀한 지표변위 산정이 가능하지만 상대적으로 노이즈가 심하고 두 영상간의 상관관계가 낮아져서 결과 산출이 어려운 지역이 많아졌다. 이와 반대로 윈도우 크기가 큰 경우 해상도가 낮아지기 때문에 국지적 지표변위를 표현하지 못하고 있다. 본 연구에서는 지표변위 관측 정확도를 높이기 위해 윈도우 크기를 달리하는 향상된 오프셋 트래킹 방법(Baek et al.,2018)을 이용하였다. 윈도우 크기를 32부터 224까지 16단계로 나누어 오프셋 트래킹을 적용하였으며 각 윈도우 크기에 따라 추정된 시선방향과 비행방향의 추정된 변위를 중첩하여 신뢰구간을 설정하고 이상치를 제거한 후 평균을 이용하여 각 방향에 대한 최종 지표변위를 생성하였다(Fig. 4(a), Fig. 4(b)). 2차원 지표변위 생성을 위해서는 비행방향에 대한 정확한 지표변위 산정이 필요하지만 대각선 방향으로 선(streak) 형태의 오차가 관측되었으며 이는 전리층에 의한 오차로 판단된다. 전리층에 의한 왜곡은 파장에 비례하며 ALOS-2의 경우 파장의 길이가 X-, C- 밴드에 비해 길기 때문에 전리층에 의한 오차가 더욱 크게 발생한다(Lee et al., 2015b). 따라서 비행방향에 발생하는 전리층 오차는 선(streak)을 따라 방향성을 지니며 이러한 방향성을 추출하는 방향 필터를 이용하여 제거하였으며 최종적으로 오차가 감소된 비행방향의 변위를 산출하였다(Fig. 4(b), Fig. 4(c)).

Fig. 4. Multi-kernel pixel tracking results, (a) range, (b) before ionosphere effect correction in azimuth direction, (c) after ionosphere effect correction in azimuth direction.

지표변위뿐만 아니라 육안 비교를 위해 고해상도 광학영상을 활용하였다. 광학영상의 경우 위성 기하차이로 인해 6차 핵실험 전·후 영상이 매칭되지 않는다. 정밀한 지표변화 비교를 위해 바위 또는 도로의 꺾이는 부분 등 시각적으로 6차 핵실험 전·후 영상에서 뚜렷한 특징을 지니는 지점에 대해 Matching Point 198점을 선정하였다. Delaunay Triangulation을 수행하여 두 영상을 Matching한 후 지표변화를 관측하였다.

2) 연구결과

(1) 영상레이더를 이용 지표변위 산정

Fig. 5는 위성 영상레이더를 이용하여 최종적으로 산출한 지표변위이다. Fig. 5(a)는 비행방향의 변위로 6차 핵실험을 기준으로 남-북방향의 서로 반대되는 변위가 발생하였음을 확인할 수 있다. 또한 Fig. 5(b)에 나타난 시선방향의 경우 만탑산 좌측부근에 약 3미터의 지표 변위가 관측되었다. 비행방향과 시선방향을 이용하여 변위의 크기와 방향을 2차원 화살표로 나타내면 Fig.5(b)와 같이 6차 핵실험을 기준으로 방사상으로 지표의 변위가 발생하였음을 확인할 수 있다. Fig. 5(c)는 지형에 따른 지표변위량을 확인하기 위해 수치표고모델(DEM: Digital Elevation Model)을 이용하여 3차원으로 지형과 같이 표현하였다. 변위 크기만을 확인하기 위해 시선방향과 비행방향의 각각 제곱의 합에 제곱근을 계산하여 절대값으로 표현하였다. 주로 경사면을 따라 지표변위가 발생하였음을 확인할 수 있으며 특히 경사가 심한 만탑산 서쪽 지역에 지표변위량이 큰 것을 확인할 수 있다.

Fig. 5. Multi-kernel pixel tracking deformation. (a) azimuth deformation, (b) range deformation, arrow represents 2D direction of deformation (c) 2D magnitude of deformation with digital elevation model.

6차 핵실험에 의해 발생한 인공지진 파형을 이용하여 핵실험의 수평 위치를 파악할 수 있지만 대부분의 지진관측소가 지표에 위치하기 때문에 깊이 추정의 정확도는 매우 낮다. 하지만 계산된 지표변위량으로부터 모델링을 통해 6차 핵실험의 깊이를 추정할 수 있으며 추가적인 정보로 활용할 수 있다. 본 연구에서 관측된 변위 값은 산 정상 부분을 기준으로 왼쪽은 (–)값을 지니며 오른쪽은 (+)값이 산출되었다. 만약 6차 핵실험에 의해 함몰이 발생하였다면 전체적으로 (–)값이 계산되지만 본 연구에서 관측된 결과는 산 정상을 기준으로 서로 상반되는 값을 지니고 있다. 따라서 힘의 방향이 방사형으로 퍼진다는 가정하에 Point 모델을 사용하는Mogi 모델(Mogi, 1958)을 적용하였다(Fig. 6). 시선방향의 지표변위를 Mogi 모델의 입력 값으로 이용하였으며(Fig. 6(a)) 모델 적용 결과(Fig. 6(b)) 핵실험이 수행된 깊이는 약 450 m로 추정되었다. 다만 잔차값을 확인하였을 때 RMSE (Root Mean Square Error)는 0.13으로 작지만 최대값이 1.2 m로 계산되었으며 이러한 잔차값은 6차 핵실험에 의해 방사형으로 에너지가 전달되었을 뿐만 아니라 산사태와 같은 다른 요소도 작용하여 나타난 것으로 판단된다. 하지만 본 연구에서 사용된 모델링 기법은 한 방향의 시선벡터 지표변위만을 사용하여 간단한 Mogi 모델을 사용한 결과로 이에 대한 추가 연구가 필요할 것으로 판된다.

Fig. 6. Mogi Modeling using range deformation (a) input, (b) model result, (c) residual.

(2) 고해상도 광학영상을 이용한 지표변화 측

위성 영상레이더의 경우 정밀 지표변위 관측은 가능하지만 직관적 판단이 어렵다. 이에 반해 고해상도 광학영상은 실제 비교를 통해 직관적 판단이 가능한 장점이 있다. Fig. 7은 분광영상(Fig. 7(a), 7(b), 7(c))과 흑백영상(Fig. 7(d), 7(e))의 핵실험 전·후 비교영상이다. 분광영상은 수목과 암석지역 구별이 용이하여 변화 관측에 활용하였으며 15곳 이상의 산사태가 발생하였음을 확인할 수 있었다. 지진파를 이용하여 6차 핵실험 약 8분 후 함몰지진으로 추정되는 지진이 관측되었으며 6차 핵실험 지점에서 남쪽으로 700 m 떨어진 것으로 추정된다고 보고되었다(Wang et al., 2018). 고해상도 광학영상을 이용하여 함몰지진 발생지점 주변을 관찰한 결과 Fig. 7의 A 지역에서 6차 핵실험 지점을 기준으로 남서쪽 약 850 m 지점에 산사태 흔적과 모양이 다른 형태의 원형 돌 무더기가 발견되었으며 그 반지름은 약 15 m, 그 면적은 740 m²으로 계산된다. 이 지역은 약 2.5 m 이상의 지표변위가 관측된 곳으로 핵실험으로 인해 함몰이 발생하여 지표면에 그 일부가 드러난 것으로 추정된다. 이는 6차 핵실험 이후 발생한 지표변위가 대부분 만탑산의 지형에 따른 산사태에 의한 지표변위이지만 북한 핵실험에 의해 발생한 함몰에 의한 지표변위도 포함하고 있을 것으로 분석할 수 있다.

Fig. 7. Optical image comparison before and after of 6th nuclear test. (a), (d): 2017-08-27, (b), (e): 2017-09-23. (c) zoom of white rectangle area in Fig. (b). (d) and (e) are zoom of area A in panchromatic image.

4. 결론

본 연구에서는 2017년 9월 3일에 발생한 북한 6차 핵실험에 대해 고해상도 위성자료를 이용하여 핵실험에 의한 지표변화를 추정하였다. 위성 레이더간섭기법의 긴밀도 지도를 이용하여 대규모 지표변위와 지표변화가 발생하였으며 그 지역의 크기는 약 9.2 km² 임을 확인하였다. 낮은 긴밀도로 인해 레이더 간섭기법 적용이 힘들 것으로 판단되며 영상레이더 강도의 상관관계를 이용하는 오프셋 트래킹 방법을 적용하였다. 사용된 윈도우 크기에 따라 계산된 변위가 달라지는 단점을 극복하기 위해 윈도우 크기를 순차적으로 증가시켜 그 결과의 통계적 처리를 수행하는 향상된 오프셋 트래킹 방법을 이용하였다. 또한 비행방향으로 생기는 전리층 왜곡을 감소시키기 위해 향 필터를 적용하였으며 최종적으로 2차원 지표변위 지도를 생성하였다. 또한 고해상도 광학영상을 이용하여 6차 핵실험 주변 여러 곳에서 산사태 흔적을 찾아냈다. 6차 핵실험 장소의 남서쪽 850m 지점에 산사태 흔적과 다른 형태의 반경 약 15 m 원형 돌무더기가 관측되었으며 이는 함로 인해 지표 일부분이 무너져 내린 곳으로 추정된다.

이와 같이 서로 다른 센서의 고해상도 위성자료를 이용하여 기존의 음파·지진파에 의한 분석을 보조할 수 있는 결과를 산출할 수 있었으며 이는 대규모 이벤트가 발생하였을 때 위성자료를 이용하여 다양한 부가정보생산이 가능함을 보여준다.