Teleseismic inversion of the 2004 Sumatra-Andaman earthquake rupture process using complete Green’s functions

Abstract

Spatial and temporal variations in coseismic slip distribution are often obtained by rupture process analyses using teleseismic body waves. Many analyses using teleseismic body waves were based on the ray theory because of the very efficiently computable direct P-, S-, and major reflected waves near the source. The 2004 Sumatra-Andaman earthquake was one of the largest earthquakes recorded in history, and the data that are required for the entire rupture process analysis include later phases such as PP waves and a very long period phase called a W phase. However, calculating these later phases using the conventional ray theoretical method is difficult. Here we investigate the rupture process of the 2004 Sumatra-Andaman earthquake using complete Green’s functions, including all later phases such as PP waves and W phase. We use the direct solution method, which computes complete synthetic seismograms up to 2 Hz for transversely isotropic spherically symmetric media, to calculate the Green’s functions. The obtained synthetic seismograms generated fit the observed seismograms quite well from a short period to a long period. The estimated slip distribution consists of four large slip areas: the largest slip occurred in the shallow part off the Sumatra west coast with a maximum slip of approximately 29 m, the second and third largest slips occurred in the shallow and deep parts of the Nicobar region with maximum slips of approximately 8 m and 7 m, respectively, and the fourth slip occurred in the middle Andaman region with a maximum slip of approximately 6 m. The estimated average rupture velocity is 2.8 km/s, but the rupture may have slowed between the Sumatra and the shallow Nicobar slip areas, and between the Nicobar and the middle Andaman slip areas. The delayed initiation of the shallow slips in the Nicobar region may possibly have been triggered by the deeper slip in the Nicobar region. There were no distinct depth-varying properties for the shallow and deep slips in the Nicobar region, as were reported for the 2011 Tohoku-Oki and 2010 Chile earthquake.