strike slip tectonics in the sw barents sea during north atlantic rifting swaen graben northern norway this study uses high quality three dimensional 3 d seismic data to investigate the occurrence of strike slip faults in the swaen graben sw barents sea the swaen graben is divided into two principal sub basins ssb1 and ssb2 the along strike and along dip displacement variations and scale relationships are analyzed for forty two 42 faults the displacement profiles for these faults are complex in the swaen graben showing clear evidence for polycyclic fault growth and marked syn sedimentary activity the observed variations in the displacement profiles indicate complex along strike segmentation linkage and mechanical interactions at distinct structural levels along dip displacement minima indicate fault reactivation by dip linkage importantly geometric evidence for strike slip faulting in the swaen graben includes the presence of extensional horsetail splay faults positive flower structures and minor transfer faults this study shows that the faults in the swaen graben developed under extensional regimes during late jurassic to early cretaceous rifting and were reactivated by regional stresses during the late cretaceous the two principal strike slip faults in the swaen graben reveal sinistral movement and are linked at a shallow depth by minor transfer faults at a relay zone our work further demonstrates the occurrence of late mesozoic strike slip movements in the sw barents sea which were induced by regional tectonics halokinesis and fault block rotation importantly strike slip faulting in the region extends perhaps into the cenozoic interacting with extension during the north atlantic rifting

grace detects coseismic and postseismic deformation from the sumatra andaman earthquake 1 we show that spherical harmonic sh solutions of the gravity recovery and climate experiment grace are now of sufficient quality to observe effects of co seismic and post seismic deformation due to the rupture from the mw 9 3 sumatra andaman earthquake on december 26 2004 and its companion nias earthquake mw 8 7 on march 28 2005 the improved ggm 03 sh level 2 solutions and improved filtering methods provide estimates with spatial resolution comparable to earlier estimates from range rate level 1 grace data the gravityfield disturbance extends over 1800 km along andaman and sunda subduction zones and changes with time following events gravity changes may be due to afterslip viscoelastic relaxation or other processes associated with dilatation satellite gravity measurements from grace provide a unique new measure of deformation and post seismic processes associated with major earthquakes especially in areas which are primarily oceanic citation chen j l c r wilson b d tapley and s grand 2007 grace detects coseismic and postseismic deformation from the sumatraandaman earthquake geophys res lett 34 l13302

postseismic rebound in fault step overs caused by pore fluid flow near field strain induced by large crustal earthquakes results in changes in pore fluid pressure that dissipate with time and produce surface deformation synthetic aperture radar sar interferometry revealed several centimeters of postseismic uplift in pull apart structures and subsidence in a compressive jog along the landers california 1992 earthquake surface rupture with a relaxation time of 270 45 days such a postseismic rebound may be explained by the transition of the poisson s ratio of the deformed volumes of rock from undrained to drained conditions as pore fluid flow allows pore pressure to return to hydrostatic equilibrium

far reaching transient motions after mojave earthquakes require broad mantle flow beneath a strong crust 1 geodeticallyobservedpostseismicsurfacedisplacements in the 7 years following the 1999 hector mine earthquake demonstrate a previously unrecognized broad pattern of transient deformation throughout southern california and into nevada more than 200 km from the epicenter unlike previous postseismic observations in which trade offs between postseismic mechanisms and the depth of flow lead to non unique solutions this deformation pattern can only be explained by viscoelastic flow in a region of the mantle 100s of km wide and below a depth of 40 km this result enables two robust conclusions regarding the nature of lithospheric strength in this region the mantle is weaker than the lower crust and flow occurs over a wide region of mantle as opposed to within a narrow shear zone beneath

coseismic slip and afterslip of the great mw 9 15 sumatra andaman earthquake of 2004 we determine coseismic and the first month postseismic deformation associated with the sumatra andaman earthquake of 26 december 2004 from near field global positioning system gps surveys in northwestern sumatra and along the nicobar andaman islands continuous and campaign gps measurements from thailand and malaysia and in situ and remotely sensed observations of the vertical motion of coral reefs the coseismic model shows that the sunda subduction mega thrust ruptured over a distance of about 1500 km and a width of less than 150 km releasing a total moment of 6 7 7 0 10 22 n m equivalent to a magnitude m w 9 15 the latitudinal distribution of released moment in our model has three distinct peaks at about 4 n 7 n and 9 n which compares well to the latitudinal variations seen in the seismic inversion and of the analysis of radiated t waves our coseismic model is also consistent with interpretation of normal modes and with the amplitude of very long period surface waves the tsunami predicted from this model fits rela tively well the altimetric measurements made by the jason and topex satellites neither slow nor delayed slip is needed to explain the normal modes and the tsunami wave the near field geodetic data that encompass both coseismic deformation and up to 40 days of postseismic deformation require that slip must have continued on the plate interface after the 500 sec long seismic rupture the postseismic geodetic moment of about 2 4 10 22 nm mw 8 8 is equal to about 30 5 of the coseismic moment release evolution of postseismic deformation is consistent with rate strengthening frictional afterslip online material summary of geodetic data used in this study

rupture process of the 2004 sumatra andaman earthquake the 26 december 2004 sumatra andaman earthquake initiated slowly with small slip and a slow rupture speed for the first 40 to 60 seconds then the rupture expanded at a speed of about 2 5 kilometers per second toward the north northwest extending 1200 to 1300 kilometers along the andaman trough peak displacements reached 15 meters along a 600 kilometer segment of the plate boundary offshore of northwestern sumatra and the southern nicobar islands slip was less in the northern 400 to 500 kilometers of the aftershock zone and at least some slip in that region may have occurred on a time scale beyond the seismic band

viscoelastic flow in the lower crust after the 1992 landers california earthquake space geodesy showed that broad scale postseismic deformation occurred after the 1992 landers earthquake three dimensional modeling shows that afterslip can only explain one horizontal component of the postseismic deformation whereas viscoelastic flow can explain the horizontal and near vertical displacements the viscosity of a weak about 10 km thick layer in the lower crust beneath the rupture zone that controls the rebound is about 10 18 pascal seconds the viscoelastic behavior of the lower crust may help to explain the extensional structures observed in the basin and range province and it may be used for the analysis of earthquake hazard

implications of deformation following the 2002 denali alaska earthquake for postseismic relaxation processes and lithospheric rheology 1 during the first 2 years following the 2002 mw 7 9 denali alaska strike slip earthquake a large array of global positioning system gps receivers recorded rapid postseismic surface motions extending at least 300 km from the rupture and at rates of more than 100 mm yr in the near field here we use three dimensional 3 d viscoelastic finite element models to infer the mechanisms responsible for these postseismic observations we consider afterslip both from an inversion of gps displacements and from stress driven forward models poroelastic rebound and viscoelastic flow in the lower crust and upper mantle several conclusions can be drawn 1 no single mechanism can explain the postseismic observations 2 significant postseismic flow below a depth of 60 km is required to explain observed far field motions best explained by a weak upper mantle with a depth dependent effective viscosity that ranges from 1019 pa s at the moho 50 km depth to 3 4 1018 pa s at 100 km depth 3 shallow afterslip within the upper crust occurs adjacent to and beneath the regions of largest coseismic slip 4 there is a contribution from deformation in the middle and lower crust from either lower crustal flow or stress driven slip afterslip is preferred over broad viscoelastic flow owing to the existence of seismic velocity discontinuities across the fault at depth though our modeling does not favor either mechanism if the process is viscoelastic relaxation the viscosity is a factor of 3 greater than the inferred mantle viscosity 5 poroelastic rebound probably contributed to the observed postseismic deformation in the immediate vicinity of the denali totschunda junction these conclusions lead us to infer an alaskan mechanical lithosphere that is about 60 km thick overlying a weak asthenosphere and a denali fault that cuts through the entire lithosphere with shear accommodated by faulting in the top 20 km and time dependent aseismic shear below

post seismic relaxation following the great 2004 sumatra andaman earthquake on a compressible self gravitating earth summary the mw 9 0 2004 december 26 sumatra andaman and mw 8 7 2005 march 28 nias earthquakes which collectively ruptured approximately 1800 km of the andaman and sunda subduction zones are expected to be followed by vigorous viscoelastic relaxation involving both the upper and lower mantle because of these large spatial dimensions it is desirable to fully account for gravitational coupling effects in the relaxation process we present a stable method of computing relaxation of a spherically stratified compressible and self gravitating viscoelastic earth following an impulsive moment release event the solution is cast in terms of a spherical harmonic expansion of viscoelastic normal modes for simple layered viscoelastic models which include a low viscosity oceanic asthenosphere we predict substantial post seismic effects over a region several 100s of km wide surrounding the eastern indian ocean we compare observed gps time series from ten regional sites mostly in thailand and indonesia beginning in 2004 december with synthetic time series that include the coseismic and post seismic effects of the 2004 december 26 and 2005 march 28 earthquakes a viscosity structure involving a biviscous burgers body rheology in the asthenosphere explains the pattern and amplitude of post seismic offsets remarkably well

co and postseismic characteristics of indian sub continent in response to the 2004 sumatra earthquake abstract the m w 9 3 sumatra earthquake of december 26 2004 caused extensive coseismic displacements globally measurements of which were made essentially using modern geodetic techniques this earthquake induced considerable perturbation in stress distribution as far as 8000 km away from the epicenteral region which is tending to relax to its normal rates as seen from postseismic transient deformation the monitoring of crustal displacements from strategically located sites using gps provides coseismic as well as postseismic deformation that facilitates the understanding of the fault geometry elastic thickness postseismic relaxation mechanisms rheology and earthquake recurrence time interval we investigated coseismic and postseismic gps derived displacements in indian region together with the gps data collected from andaman and sumatra region it is found that while ew displacements are significantly large in peninsular india those in the region to the north of central india tectonic zone citz are relatively small we could delineate the postseismic transients from position time series and interpreted them in terms of viscoelastic relaxation it is inferred that the postseismic deformation is characterized by a power law viscoelastic flow in the mantle in indian peninsula region the timescale parameter of the exponential decay 250 days would require an extremely low viscosity for the upper mantle relying on the prevailing coseismic and postseismic displacement fields the present study also reflects upon the contemporary litho tectonics of the indian sub continent

subcrustal low velocity layers in central india and their implications a 2 d subcrustal velocity model for the central indian continental lithosphere has been derived by travel time and relative amplitude modeling of a digitally normalized analog seismic record section of the hirapur mandla dss profile using a ray tracing technique some prominent wave groups with apparent velocities slightly higher than the moho reflection phase pmp are identified on the normalized record sections assembled with a reduction velocity of 6 km s 1 we interpret these phases as the wide angle reflections from subcrustal lithospheric boundaries comparison of synthetic seismograms with the observed record section shows that the observed phases cannot be explained either by multiples or by the p to s converted phase pms from the moho subcrustal velocity models either with a velocity increase or with a single low velocity layer lvl also do not provide a satisfactory fit we infer that a subcrustal velocity model with two alternate lvls velocity 7 2 km s 1 separated by a 6 km thick high velocity layer velocity 8 1 km s 1 can satisfy both the observed travel times and amplitudes the prominent reflection phases are modeled at depths of 49 51 57 and 60 km it is inferred that the subcrustal lithosphere in the central indian region has a lamellar structure with varying structural and mechanical properties the alternating lvls occurring at relatively shallow depths below moho may be associated with the zones of weakness and lower viscosity suggesting continued mobility with a possible thermal source in the upper mantle this explains the source of observed high heat flow values in the central indian region

interpretation of shallow electromagnetic instruments resistivity and magnetic susceptibility measurements using rapid 1d 3d inversion we have developed an inversion process of electromagnetic induction emi data based on a two step approach with 1d inversion of the entire studied surface and a fast 3d inversion applied over limited areas this process is similar to that formerly used in resistivity prospection for the study of soil environmental engineering or archaeological explorations low frequency electromagnetic instruments referred to as slingram emi have highly useful specificities they are light are easy to move in the field and can simultaneously measure the ground s electric conductivity and magnetic susceptibility they have thus been used to map these properties over large surface areas within relatively short periods of time and at reasonable expense the possibility of combining several coil geometries has opened up the potential for multidepth techniques and systematic 1d inversion which are found to be sufficiently revealing to allow larger portions of surveyed areas to be analyzed in the targeted areas selected for 3d inversion the geometries of the 3d features and the resistivity and or susceptibility contrasts are determined this step is based on the method of moments where only 3d heterogeneities are meshed and only a small number of major characteristics such as contrast thickness width etc are sought we first applied this process to synthetic data then to data acquired at an experimental test site and finally to field cases the rapid 3d inversion complements the 1d inversion by solving a series of issues correction for the apparent anisotropy generated by the instrument configuration multiarched anomalies precise location of lateral changes and determination of the properties contrasts our inversion results highlighted the importance of the instrument geometry we also have determined that apparent magnetic susceptibility data can be more appropriate for the determination of the volume of man made features and can be highly complementary to conductivity data

recognition of pre grenvillian and grenvillian tectonothermal events in the central indian tectonic zones implications on rodinian crustal assembly

plate boundary deformation associated with the great sumatra andaman earthquake the sumatra andaman earthquake of 26 december 2004 is the first giant earthquake moment magnitude m w 9 0 to have occurred since the advent of modern space based geodesy and broadband seismology it therefore provides an unprecedented opportunity to investigate the characteristics of one of these enormous and rare events here we report estimates of the ground displacement associated with this event using near field global positioning system gps surveys in northwestern sumatra combined with in situ and remote observations of the vertical motion of coral reefs these data show that the earthquake was generated by rupture of the sunda subduction megathrust over a distance of 1 500 kilometres and a width of 150 kilometres megathrust slip exceeded 20 metres offshore northern sumatra mostly at depths shallower than 30 kilometres comparison of the geodetically and seismically inferred slip distribution indicates that 30 per cent additional fault slip accrued in the 1 5 months following the 500 second long seismic rupture both seismic and aseismic slip before our re occupation of gps sites occurred on the shallow portion of the megathrust where the large aceh tsunami originated slip tapers off abruptly along strike beneath simeulue island at the southeastern edge of the rupture where the earthquake nucleated and where an m w 7 2 earthquake occurred in late 2002 this edge also abuts the northern limit of slip in the 28 march 2005 m w 8 7 nias simeulue earthquake

retrieving earthquake signature in grace gravity solutions summary the grace satellites have been orbiting the earth since 2002 monitoring the time variable gravity field some of the observed fluctuations are due to geodynamic causes but they are often hidden in the complex signal composed of hydrology ocean atmosphere and geodynamics the signal of geodynamic origin being usually the smallest in addition dealiasing residuals and noise make the separation of the signal from the different causes more difficult we proposed a method based on the empirical orthogonal function decomposition to extract the signal of physical origin under the hypothesis that the physical signal is spatially more consistent than the noise and aliasing incomplete correction we used synthetic geoid variations associated with earthquakes located at nearly 2000 positions at the earth surface based on several examples of large actual subduction events we show that with the present day accuracy we can retrieve the geoid variations associated with more than 98 per cent of the earthquakes of magnitude 9 or above around 60 per cent for magnitude 8 8 40 per cent for magnitude 8 6 and 33 per cent for magnitude 8 3 some events with the right properties and location can be detected with magnitude as low as 8 we then applied the method to the grace solutions and retrieved the hokkaido event 2003 and the sumatra event 2004 which is in agreement with the retrieval rates mentioned here above

mobility of continental mantle evidence from postseismic geodetic observations following the 1992 landers earthquake the crust around the rupture zone of the 1992 landers earthquake has continued to deform in the years following the earthquake at rates 3 times greater than pre earthquake rates we use a combination of global positioning system gps and synthetic aperture radar insar data collected during a 3 year epoch following the earthquake in order to investigate postseismic mechanisms responsible for the high transient velocities in order to maximize the potential signal from viscoelastic relaxation we evaluate and model postseismic relaxation following the first few months of documented accelerated deformation the combination of gps and insar data allows us to establish viscoelastic relaxation of the lower crust and upper mantle as the dominant postseismic process and to discriminate among possible viscoelastic models the data particularly require the presence of a highly ductile uppermost mantle beneath the central mojave domain with temperature between the wet and dry basalt solidus this is consistent with independent seismic and geochemical inferences of a regionally warm uppermost mantle further consideration of seismic velocity variations in conjunction with faulting patterns within the mojave desert suggests that the primary faulting characteristics of the mojave desert namely the pervasive late cenozoic deformation within the eastern california shear zone versus the near absence of faults in the western mojave domain are controlled by the rheology of the uppermost mantle

postseismic relaxation theory on the spherical earth abstract postseismic displacements and strains observed at the earth9s surface can be explained through the relaxation of a viscoelastic asthenosphere underlying a purely elastic crust a theory is presented for modeling the effects of postseismic relaxation on a spherically symmetric earth for an earthquake point source located in the upper elastic layer the displacement field is decomposed into its toroidal and spheroidal components a linear maxwell rheology in the viscoelastic layers is assumed enabling the use of the correspondence principle for the solutions of the equations of static equilibrium displacement and strain fields are then calculated using normal mode summation computational results are presented for two simple earthquake sources a strike slip fault and a uniaxial thrust fault the patterns of postseismic displacements and strains are found to depend strongly on both the earth model and earthquake source geometry in particular the elastic plate thickness asthenosphere thickness fault type and fault length each play a major role in determining the spatial pattern of postseismic relaxation effects asthenospheric viscosity controls the temporal pattern of relaxation

constraints on the viscosity of the continental crust and mantle from gps measurements and postseismic deformation models in western mongolia 1 we use gps measurements and models of postseismic deformation caused by seven m6 8 to 8 4 earthquakes that occurred in the past 100 years in mongolia to assess the viscosity of the lower crust and upper mantle we find an upper mantle viscosity between 1 1018 and 4 1018 pa s the presence of such a weak mantle is consistent with results from independent seismological and petrological studies that show an abnormally hot upper mantle beneath mongolia the viscosity of the lower crust is less well constrained but a weak lower crust 3 1016 to 2 1017 pa s is preferred by the data using our best fit upper mantle and lower crust viscosities we find that the postseismic effects of viscoelastic relaxation on present day horizontal gps velocities are small 2 mm yr 1 but still persist 100 years after the 1905 m8 4 bolnay earthquake this study shows that the gps velocity field in the baikal mongolia area can be modeled as the sum of 1 a rigid translation and rotation of the whole network 2 a 3 5 mm yr 1 simple shear velocity gradient between the siberian platform to the north and northern china to the south and 3 the contribution of postseismic deformation mostly caused by the 1905 bolnay tsetserleg sequence and by the smaller but more recent 1957 bogd earthquake

insight into the 2004 sumatra andaman earthquake from gps measurements in southeast asia data collected at 60 global positioning system gps sites in southeast asia show the crustal deformation caused by the 26 december 2004 sumatra andaman earthquake at an unprecedented large scale small but significant co seismic jumps are clearly detected more than 3 000 km from the earthquake epicentre the nearest sites still more than 400 km away show displacements of 10 cm or more here we show that the rupture plane for this earthquake must have been at least 1 000 km long and that non homogeneous slip is required to fit the large displacement gradients revealed by the gps measurements our kinematic analysis of the gps recordings indicates that the centroid of released deformation is located at least 200 km north of the seismological epicentre it also provides evidence that the rupture propagated northward sufficiently fast for stations in northern thailand to have reached their final positions less than 10 min after the earthquake hence ruling out the hypothesis of a silent slow aseismic rupture data culled from gps stations in indonesia malaysia singapore thailand india and china have answered more questions relating to the 26 december sumatra andaman megathrust earthquake permanent displacements were seen 3 000 km away final displacements occurred within 10 minutes of the earthquake ruling out slow aseismic slip on the fault

uplift and subsidence associated with the great aceh andaman earthquake of 2004 rupture of the sunda megathrust on 26 december 2004 produced broad regions of uplift and subsidence we define the pivot line separating these regions as a first step in defining the lateral extent and the downdip limit of rupture during that great m w 9 2 earthquake in the region of the andaman and nicobar islands we rely exclusively on the interpretation of satellite imagery and a tidal model at the southern limit of the great rupture we rely principally on field measurements of emerged coral microatolls uplift extends from the middle of simeulue island sumatra at 2 5 n to preparis island myanmar burma at 14 9 n thus the rupture is 1600 km long the distance from the pivot line to the trench varies appreciably the northern and western andaman islands rose whereas the southern and eastern portion of the islands subsided the nicobar islands and the west coast of aceh province sumatra subsided tilt at the southern end of the rupture is steep the distance from 1 5 m of uplift to the pivot line is just 60 km our method of using satellite imagery to recognize changes in elevation relative to sea surface height and of using a tidal model to place quantitative bounds on coseismic uplift or subsidence is a novel approach that can be adapted to other forms of remote sensing and can be applied to other subduction zones in tropical regions

global positioning system constraints on plate kinematics and dynamics in the eastern mediterranean and caucasus we present and interpret global positioning system gps measurements of crustal motions for the period 1988 1997 at 189 sites extending east west from the caucasus mountains to the adriatic sea and north south from the southern edge of the eurasian plate to the northern edge of the african plate sites on the northern arabian platform move 18 2 mm yr at n25 5 w relative to eurasia less than the nuvel 1a circuit closure rate 25 1 mm yr at n21 7 w preliminary motion estimates 1994 1997 for stations located in egypt on the northeastern part of africa show northward motion at 5 6 2 mm yr also slower than nuvel ia estimates 10 1 mm yr at n2 4 e eastern turkey is characterized by distributed deformation while central turkey is characterized by coherent plate motion internal deformation of 2 mm yr involving westward displacement and counterclockwise rotation of the anatolian plate the anatolian plate is de coupled from eurasia along the right lateral strike slip north anatolian fault naf we derive a best fitting euler vector for anatolia eurasia motion of 30 7 0 8 n 32 6 0 4 e 1 2 0 1 myr the euler vector gives an upper bound for naf slip rate of 24 1 mm yr we determine a preliminary gps arabia anatolia euler vector of 32 9 1 2 n 40 3 1 1 e 0 8 0 2 myr and an upper bound on left lateral slip on the east anatolian fault eaf of 9 1 mm yr the central and southern aegean is characterized by coherent motion internal deformation of 2 mm yr toward the sw at 30 1 mm yr relative to eurasia stations in the se aegean deviate significantly from the overall motion of the southern aegean showing increasing velocities toward the trench and reaching 10 1 mm yr relative to the southern aegean as a whole

regional conductance map of andaman and nicobar region abstract to image the electrical conductivity distribution fluxgate magnetometers are operated at five sites in andaman and nicobar region transfer functions are estimated for the period range 8 128 min from nighttime transient geomagnetic variations using robust regression analysis the observed induction arrows in andaman islands are found to point towards east despite deep sea located towards its west this indicates that fore arc basin andaman nicobar deep is more conducting than the region of outer non volcanic island arc thin sheet model requires the conductance of 10 000 35 000 s with increase conductivity towards the south for explaining the observed induction pattern the observed induction pattern at andaman nicobar stations can be explained in terms of high conducting cretaceous tertiary sediments filling the andaman nicobar deep high conductivity over invisible bank has been attributed to the partial melts volatile fluids derived from the subducting indian plate that are intruding into the eastern margin of fore arc basin through the west andaman fault waf the induction pattern at great nicobar station campbell bay may be related to the highly conducting sediments filling the mergui basin along with mafic intrusions also crustal transition occurs below the mergui terrace at the malayan coast contributing to the enhanced conductivity anomaly

the size and duration of the sumatra andaman earthquake from far field static offsets the 26 december 2004 sumatra earthquake produced static offsets at continuously operating gps stations at distances of up to 4500 kilometers from the epicenter we used these displacements to model the earthquake and include consideration of the earth s shape and depth varying rigidity the results imply that the average slip was 5 meters along the full length of the rupture including the approximately 650 kilometer long andaman segment comparison of the source derived from the far field static offsets with seismically derived estimates suggests that 25 to 35 of the total moment release occurred at periods greater than 1 hour taking into consideration the strong dip dependence of moment estimates the magnitude of the earthquake did not exceed mw 9 2

implications of postseismic gravity change following the great 2004 sumatra andaman earthquake from the regional harmonic analysis of grace intersatellite tracking data 1 we report gravity recovery and climate experiment grace satellite observations of coseismic displacements and postseismic transients from the great sumatra andaman islands thrust event mw 9 2 earthquake in december 2004 instead of using global spherical harmonic solutions of monthly gravity fields we estimated the gravity changes directly using intersatellite range rate data with regionally concentrated spherical slepian basis functions every 15 day interval we found significant step like coseismic and exponential like postseismic behavior in the time series of estimated coefficients from may 2003 to april 2007 for the spherical slepian functions after deriving coseismic slip estimates from seismic and geodetic data that spanned different time intervals we estimated and evaluated postseismic relaxation mechanisms with alternate asthenosphere viscosity models the large spatial coverage and uniform accuracy of our grace solution enabled us to clearly delineate a postseismic transient signal in the first 2 years of postearthquake grace data our preferred interpretation of the long wavelength components of the postseismic gravity change is biviscous viscoelastic flow we estimated a transient viscosity of 5 1017 pa s and a steady state viscosity of 5 1018 1019 pa s additional years of the grace observations should provide improved steady state viscosity estimates in contrast to our interpretation of coseismic gravity change the prominent postearthquake positive gravity change around the nicobar islands is accounted for by seafloor uplift with less postseismic perturbation in intrinsic density in the region surrounding the earthquake