Retrieving earthquake signature in grace gravity solutions

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.     

Improved ocean-geoid resolution from retracked ERS-1 satellite altimeter waveforms

The ocean geoid can be inferred from the topography of the mean sea surface. Satellite altimeters transmit radar pulses and determine the return traveltime to measure sea-surface height. The ERS-1 altimeter stacks 51 consecutive radar reflections on board the satellite to a single waveform. Tracking the time shift of the waveform gives an estimate of the distance to the sea surface. We retrack the ERS-1 radar traveltimes using a model that is focused on the leading edge of the waveforms. While earlier methods regarded adjacent waveforms as independent statistical events, we invert a whole sequence of waveforms simultaneously for a spline geoid solution. Smoothness is controlled by spectral constraints on the spline coefficients. Our geoid solutions have an average spectral density equal to the expected power spectrum of the true geoid. The coherence of repeat track solutions indicates a spatial resolution of 31  km, as compared to 41  km resolution for the ERS-1 Ocean Product. While the resolution of the latter deteriorates to 47  km for wave heights above 2  m, our geoid solution maintains its resolution of 31  km for rough sea. Retracking altimeter waveform data and constraining the solution by a spectral model leads to a realistic geoid solution with significantly improved along-track resolution.     

Coseismic and post-seismic signatures of the Sumatra 2004 December and 2005 March earthquakes in GRACE satellite gravity

The GRACE satellite mission has been measuring the Earth's gravity field and its temporal variations since 2002 April. Although these variations are mainly due to mass transfer within the geofluid envelops, they also result from mass displacements associated with phenomena including glacial isostatic adjustment and earthquakes. However, these last contributions are difficult to isolate because of the presence of noise and of geofluid signals, and because of GRACE's coarse spatial resolution (>400 km half-wavelength). In this paper, we show that a wavelet analysis on the sphere helps to retrieve earthquake signatures from GRACE geoid products. Using a wavelet analysis of GRACE geoids products, we show that the geoid variations caused by the 2004 December ( M _w= 9.2) and 2005 March ( M _w= 8.7) Sumatra earthquakes can be detected. At GRACE resolution, the 2004 December earthquake produced a strong coseismic decrease of the gravity field in the Andaman Sea, followed by relaxation in the area affected by both the Andaman 2004 and the Nias 2005 earthquakes. We find two characteristic timescales for the relaxation, with a fast variation occurring in the vicinity of the Central Andaman ridge. We discuss our coseismic observations in terms of density changes of crustal and upper-mantle rocks, and of the vertical displacements in the Andaman Sea. We interpret the post-seismic signal in terms of the viscoelastic response of the Earth's mantle. The transient component of the relaxation may indicate the presence of hot, viscous material beneath the active Central Andaman Basin.     

Numerical aspects of a spline-based multiresolution recovery of the harmonic mass density out of gravity functionals

We show the numerical applicability of a multiresolution method based on harmonic splines on the 3-D ball which allows the regularized recovery of the harmonic part of the Earth's mass density distribution out of different types of gravity data, for example, different radial derivatives of the potential, at various positions which need not be located on a common sphere. This approximated harmonic density can be combined with its orthogonal anharmonic complement, for example, determined out of the splitting function of free oscillations, to an approximation of the whole mass density function. The applicability of the presented tool is demonstrated by several test calculations based on simulated gravity values derived from EGM96. The method yields a multiresolution in the sense that the localization of the constructed spline basis functions can be increased which yields in combination with more data a higher resolution of the resulting spline. Moreover, we show that a locally improved data situation allows a highly resolved recovery in this particular area in combination with a coarse approximation elsewhere which is an essential advantage of this method, for example, compared to polynomial approximation.