Sumatra-Andaman Megathrust Earthquake Slip: Insights From Mechanical Modeling of ICESat Surface Deformation Measurements

Modeling the kinematics of the 2004 Great Sumatra-Andaman earthquake is limited in the northern two-thirds of the rupture zone by a scarcity of near-rupture geodetic deformation measurements. Precisely repeated Ice, Cloud, and Land Elevation Satellite (ICESat) profiles across the Andaman and Nicobar Islands provide a means to more fully document the spatial pattern of surface vertical displacements and thus better constrain geomechanical modeling of the slip distribution. ICESat profiles that total ~45 km in length cross Car Nicobar, Kamorta, and Katchall in the Nicobar chain. Within the Andamans, the coverage includes ~350 km on North, Central, and South Andaman Islands along two NNE and NNW-trending profiles that provide elevations on both the east and west coasts of the island chain. Two profiles totaling ~80 km in length cross South Sentinel Island, and one profile ~10 km long crosses North Sentinel Island. With an average laser footprint spacing of 175 m, the total coverage provides over 2700 georeferenced surface elevations measurements for each operations period. Laser backscatter waveforms recorded for each footprint enable detection of forest canopy top and underlying ground elevations with decimeter vertical precision. Surface elevation change is determined from elevation profiles, acquired before and after the earthquake, that are repeated with a cross-track separation of less than 100 m by precision pointing of the ICESat spacecraft. Apparent elevation changes associated with cross-track offsets are corrected according to local slopes calculated from multiple post-earthquake repeated profiles. The surface deformation measurements recorded by ICESat are generally consistent with the spatial distribution of uplift predicted by a preliminary slip distribution model. To predict co-seismic surface deformation, we apply a slip distribution, derived from the released energy distribution computed by Ishii et al. (2005), as the displacement discontinuity boundary condition on the Sumatra-Andaman subduction interface fault. The direction of slip on the fault surface is derived from the slip directions computed by Tsai et al. (in review) for centroid moment tensor focal mechanisms spatially distributed along the rupture. The slip model will be refined to better correspond to the observed surface deformation as additional results from the ICESat profiles become available.