Estimation of postseismic deformation parameters from continuous GPS data in northern Sumatra after the 2004 Sumatra-Andaman earthquake

Continuous global positioning system (GPS) in northern Sumatra detected signal of the ongoing physical process of postseismic deformation after the M 9.2 2004 Sumatra-Andaman earthquake. We analyze the characteristics of postseismic deformation of the 2004 earthquake based on GPS networks operated by BIG, and the others named AGNeSS, and SuGAr networks located in northern Sumatra. We use a simple analytical logarithmic and exponential function to evaluate the postseismic deformation parameters of the 2004 earthquake. We ﬁnd that GPS data in northern Sumatra during time periods of 2005–2012 are ﬁt better using the logarithmic function with s log of 104.2 ± 0.1 than using the exponential function. Our result clearly indicates that other physical mechanisms of postseismic deformation should be taken into account rather than a single physical mechanism of afterslip only.


Introduction
Indonesia is surrounded by four major plates; they are Eurasian plate, Indo-Australian plate, Pacific plate, and Philippine plate. This complex tectonic condition may cause earthquakes. These earthquakes possibly occur in the area of subduction zones, such as the 2005 M8.7 Nias earthquake (e.g., Kreemer et al. 2006), the 2010 M7.8 Mentawai earthquake (e.g., Hill et al. 2012), the 1994 M7.8 and 2006 M7.8 Java tsunami earthquake (e.g., Bilek and Engdahl 2007), and along the inland faults, such as the 1987 M6.6 Toru earthquake (Sieh and Natawidjaja 2000) and the 2013 M6.1 Aceh earthquake (Ito et al., manuscript in preparation, 2015). One of the active inland fault zones is the Great Sumatran Fault (GSF). GSF, one of the world's greatest fault zones located inland of Sumatra, consists of 20 segments (Sieh and Natawidjaja 2000) starting from the southern Sumatra with a plate movement of 12 mm/a (Meilano, Susilo and Gunawan, manuscript in preparation, 2015) to the northern Sumatra with the increasingly rapid movements of 20 mm/a (Ito et al. 2012). Historically, magnitude 7 class earthquakes have been occured along GSF (Sieh and Natawidjaja 2000).
One of the large megathrust earthquakes occured due to subduction process between Indo-Australian plate and Eurasian plate was the 26 December 2004 M9.2 Sumatra-Andaman earthquake ( Fig. 1) (e.g., Subarya et al. 2006). According to the geodetic and geological estimation, the continental shelf of Southeast Asia moves independently than Eurasian plate defined as Sundaland block (e.g., Bock et al. 2003). The 2004 Sumatra-Andaman earthquake resulted from devastating tsunami with the lost of more than 280,000 inhabitants (e.g., Sieh 2006). This earthquake is considered as a megathrust earthquake with a longest rupture compared with other earthquakes (Shearer and Bürgmann 2010). Rupture along the *1600 km starts from the epicenter near the island of Simeulue (3.3160°N, 95.8540°E) and continuously up to the Andaman Islands with a speed of 1 km/s (Fujii and Satake 2007). That rupture completed within about 10 min and be a source of nuisance of sea water volume in which became a devastating tsunami (Chlieh et al. 2007).
After the 2004 Sumatra-Andaman earthquake, governments and institutions of many countries have collaborated together to understand deeply this earthquake. Many researchers have been trying to understand the characteristics of the 2004 Sumatra-Andaman earthquake so that they might determine the further earthquakes occurrence potential as well as analyze and evaluate the disaster mitigation in northern Sumatra. For instance, the cooperation between the universities from Indonesia and Japan resulted in the project of global positioning system (GPS) observation spread across several regions in northern Sumatra. This GPS network named AGNeSS (Aceh GPS Network for the Sumatran Fault System) (Ito et al. 2012;Gunawan et al. 2014). Previous studies suggested that there are ongoing postseismic activities due to the 2004 Sumatra-Andaman earthquake with the movement of 5-15 cm per year observed in the Andaman Islands (e.g., Som et al. 2011Som et al. , 2014, and also in Thailand (e.g., Panet et al. 2010). Furthermore, Gunawan et al. (2014) suggested that postseismic deformation in the long-term period after earthquake occurrences corresponds to viscous flow in the asthenosphere. They showed that contribution of afterslip during one year periods from end of 2005 to end of 2006 corresponds to *60% of GPS data displacements.
This study analyzes the continuous GPS data available in northern Sumatra which clearly detect the significant signal of postseismic deformation after the 2004 Sumatra-Andaman earthquake. We use more new continuous GPS data observed by the network that is installed and maintained by Geospatial Information Agency of Indonesia (BIG) in northern Sumatra. This GPS data is verẙ Inset shows larger regional setting important because it has never been published by any previous researches before. Additionally, this continuous GPS data may capture the ongoing process of postseismic deformation in northern Sumatra after the 2004 Sumatra-Andaman earthquake. Here, we focus on determining the postseismic deformation parameters of the GPS data time series using simple analytical logarithmic and exponential functions.

Data and method
Data used in this study are continuous GPS data observed in northern Sumatra. One GPS network is operated by BIG, and the other GPS networks are named AGNeSS, and SuGAr (Sumatran GPS Array). Figure 1 shows the distribution of GPS sites used in this study. GPS sites of CSAB, LHMI, MEUL, and SAMP were installed and maintained by BIG. Among these sites, SAMP is the oldest one which has been established since 1997. Then, MEUL was constructed in 2008, followed by LHMI in 2010 and CSAB in 2011. ACEH site is part of AGNeSS, has been installed and maintained since 2005 by the collaboration of Nagoya University, Kochi University, Tohoku University, Bandung Institute of Technology, and Syiah Kuala University (Ito et al. 2012). BSIM, LEWK, and UMLH sites are part of SuGAr installed and maintained since 2005 by Earth Observatory of Singapore (EOS) and Indonesian Institute of Sciences (LIPI) .
Time series of ACEH was published by previous studies (Ito et al. 2012;Gunawan et al. 2014), so we used their result for further analysis. For other sites, we performed GPS data processing using GAMIT 10.5 to obtain daily solutions (Herring et al. 2010). We used three steps to process the data with GAMIT 10.5. First step is data preparation. It means to prepare raw phase and pseudorange data in the form of ASCII X-files, station coordinates in the form of L-file, receiver and antenna information for each site, satellite list and scenario, initial conditions for the satellites orbits in a G-file, satellite and station clock values, control files for the analysis, standard tables to provide lunar/solar ephemerides, the Earth's rotation, geodetic datums, and instrumentation information. Second, creating model of each parameter to solve the phase ambiguity to obtain a preliminary coordinate of each site observation. Third, we fix the preliminary coordinate result to IGS sites such as KIT3, HYDE, COCO, DGAR, IISC, KUNM, CUSV, NTUS, BAKO, XMIS, TNML, TWTF, PIMO, DARW, KARR, YAR2, CNMR, GUAM, TOW2, KOUC, and obtain daily solutions in International Reference Frame (ITRF) 2008 (Altamimi et al. 2011).
In order to make the GPS data much more appropriate for tectonics discussion, we need to use the time series data in the Sundaland block reference frame. So after obtaining daily solutions in ITRF2008, we transformed these time series data to obtain daily solutions in the Sundaland block. We use rotation pole proposed by Simons et al. (2007) with parameters of 49.0°N, -94.2°E, 0.336°/Ma, which is calculated using 28 GPS sites located in South East Asia region. In order to use these parameters correctly, first, we transformed the daily solutions from ITRF2008 into ITRF2000 (Altamimi et al. 2011). Second, we transformed daily solutions from ITRF2000 reference frame into Sundaland block reference frame.
We matched these data based on characteristics of afterslip using logarithmic function (Marone et al. 1991), defined as u(t) = c ? a ln (1 ? t/slog). We compare the results with other definition of data fit using exponential function (Savage and Prescott 1978), defined as u(t) = c ? a (1-e -t/sexp ). In these functions, u(t) is horizontal component position, c is data offset and a is amplitudes, while slog and sexp are decay time of afterslip for logarithmic function and exponential function. We fit the time series horizontal data to these two functions, and search for the unknown parameters of c, a, slog, and sexp, respectively.

Time series analysis
The daily solutions of GPS data processing were used to analyze and determine the characteristics of each GPS observation site and the quality of observational data distribution. The displacement characteristics of each observation station were then determined by analyzing the time series of the horizontal components in each GPS site that refer to the Sundaland block reference frame. Figure 2 shows the daily solutions of each GPS site in northern Sumatra.
Previous studies suggested that northern Sumatra experience significant coseismic and postseismic deformation of the 12 April 2012 M8.6 Indian Ocean earthquake (e.g., Yadav et al. 2012;Wei et al. 2013;Hill et al. 2015;Han et al. 2015). In addition, GPS data in Java also experienced crustal deformation due to this earthquake (Gunawan et al., manuscript in preparation, 2015). Since our analysis focused on the 2004 Sumatra-Andaman earthquake, we did not use the GPS data after the 2012 earthquake for further analysis as it was clearly influenced by the 2012 earthquake.
The significant signal of postseismic deformation of the 2005 Nias earthquake was detected by BSIM and LEWK. Previous studies suggested that significant postseismic deformation of afterslip after the 2005 Nias earthquake occurred at shallow portion of the main rupture at 97°E, 1°N (Hsu et al. 2006;Kreemer et al. 2006). They suggested that afterslip influenced deformation observed by GPS surrounding Nias Island. On the other hand, due to location of GPS sites in northern Sumatra much closer to the rupture of the 2004 Sumatra-Andaman earthquake and much larger magnitude size of this earthquake than the 2005 Nias earthquake, previous studies suggested that GPS data in northern Sumatra are most significantly influenced by the 2004 earthquake (Subarya et al. 2006;Gunawan et al. 2014). That was why BSIM and LEWK would not be used for further analysis.
Similarly to the cases of the 2008 Simeulue earthquake and the 2010 Northern Sumatra earthquake, where location of GPS sites in northern Sumatra was much closer to the rupture of the 2004 Sumatra-Andaman earthquake and much larger magnitude size of the 2004 earthquake than those of the 2008 and 2010 earthquakes. So for MEUL site, we removed only the coseismic displacements of the 2010 earthquake with -0.0184 m easting and 0.0385 m northing, as postseismic deformation, was insignificant compared to the 2004 earthquake. We concluded that in this study, we use only GPS data from ACEH, UMLH, CSAB, SAMP, MEUL, and LHMI.

Postseismic deformation parameters
Using the corrected daily solutions, we matched the GPS data to logarithmic and exponential functions. The postseismic parameters of logarithmic and exponential functions for each GPS site were shown in Tables 1 and 2. In our calculation, we found that slog is 104.2 ± 0.1 days while sexp is 704.5 ± 0.1 days, respectively. Based on the best matching component of logarithmic and exponential functions above, our results indicated that postseismic parameters of logarithmic and exponential functions for each GPS site have different results. We found that misfit was obtained for GPS data using the models of logarithmic functions, that was 5.9, smaller than the exponential ones, that was 16.9 (Tables 1, 2). This result indicates that misfit of fitting time series data to logarithmic function obtained *2.8 times smaller than that to the exponential function.
This showed that for 104.2 days after the earthquake, deformation mechanism was dominated by afterslip deformation due to Sumatra-Andaman earthquake in 2004. Since the time series, GPS data captured postseismic deformation for more than 100 days, so the other postseismic mechanisms should be taken into account. This could indicate that the postseismic deformation of the Sumatra-Andaman earthquake was not only caused by a single mechanism of afterslip, but also influenced by the other mechanisms of postseismic deformation, such as viscoelastic relaxation, for GPS data of 100 days forward.
Our finding supported previous research, which suggested that afterslip corresponds to early deformation after earthquake occurrences (Gunawan et al. 2014). In this study, we showed a detail result of postseismic parameters of the 2004 Sumatra-Andaman earthquake. Using a simple analysis model of logarithmic function, we could obtain the ongoing physical mechanisms process after earthquake occurrences exhaustively.
The difference between data fitting to logarithmic and exponential functions was significantly observed and visible during early periods after the 2004 earthquake, especially at ACEH and UMLH sites where early time series data are available (Fig. 2). This indicates that ACEH and UMLH play major role in our analysis. In the longer time series data, discrepancy between logarithmic and exponential functions to fit the data is invisible. Our analysis suggests the importance of maintaining continuous GPS data in prior, during, and after earthquake occurrences for further hazard analysis.

Conclusions
We show that GPS sites in northern Sumatra clearly detect significant postseismic deformation after the 2004 Sumatra-Andaman earthquake. We conclude that GPS data in northern Sumatra during the period of 2005-2012 are very well matched using logarithmic function with smaller average misfit of *2.8 times compared to matching using exponential function. The decay time of postseismic Earthq Sci (2015) 28(5-6):347-352 351 deformation of 104.2 days indicates that other mechanisms of postseismic deformation should be taken into account rather than a single physical mechanism of afterslip only.