이 연구에서는 기존의 탄성파 모델링 알고리즘에 지진 송신원을 적용하고, 음향-탄성파 결합 매질을 구현하여 남극대륙 주변과 같은 극지해역에서 발생할 수 있는 지진파의 거동을 모사한다. 기존의 변위근사 유한차분법과 달리 속도-응력 식으로 구성되는 엇격자 유한차분법의 경우 다양한 송신원을 구현하는데 적합하므로 변위-속도-응력 식에 기초하여 개발된 3차원 엇격자 유한차분법 알고리즘과 이중 우력(Double Couple Forces)을 이용하여 구현한 지진 송신원을 접목시켜 지진파의 거동을 모사한다. 좌수향 주향이동단층, 정단층, 역단층 형태의 지진 송신원에 대해서 개발된 알고리즘을 검증한 결과 이론적으로 예측되는 P파의 초동을 정확히 모사할 수 있었고, 섭입대 모델에 대한 수치모형실험 결과 섭입대에서 역단층에 의해 발생된 후 대륙지각, 해양지각 및 해양에서 전파되는 지진파의 거동양상이 정확하게 모사되는 것을 확인할 수 있었다.
We simulate the propagation of earthquake waves in the continental margin of Antarctica using the elastic wave modeling algorithm, which is modified to be suitable for acoustic-elastic coupled media and earthquake source. To simulate the various types of earthquake source, the staggered-grid finite-difference method, which is composed of velocity-stress formulae, can be more appropriate to use than the conventional, displacement-based, finite-difference method. We simulate the elastic wave propagation generated by earthquakes combining 3D staggered-grid finite-difference algorithm composed of displacement-velocity-stress formulae with double couple mechanisms for earthquake source. Through numerical tests for left-lateral strike-slip fault, normal fault and reverse fault, we could confirm that the first arrival of P waves at the surface is in a good agreement with the theoretically-predicted results based on the focal mechanism of an earthquake. Numerical results for a model made after the subduction zone in the continental margin of Antarctica showed that earthquake waves, generated by the reverse fault and propagating through the continental crust, the oceanic crust and the ocean, are accurately described.