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.     

Validation and fusion of different databases in preparation of high-resolution geoid determination

An up to date determination of a high-resolution geoid requires the use of best available databases concerning digital terrain model (DTM), bathymetry, global geopotential model and gravity field. The occasion to revisit methods to validate and merge different data sets has been created by a new project for the determination of a new European Geoid.
Since the computation of the latest European geoid and quasi-geoid model (EGG97), significant new or improved data sets have become available, such as new global geopotential models from CHAMP and GRACE missions, new national and global DTMs and new or upgraded gravity data sets.
In the context of the new European Gravity and Geoid Project (EGGP), within the IAG Commission 2, some data validation tests have been performed in the Italian zone.
In the area 19°× 17° wide, covering Italy, three kinds of tests have been performed: comparison among different DTMs in order to choose the best one to be used; comparisons in terms of geoid computation in some coastal areas, to evaluate bathymetry effects, and the validation of the EIGEN-CG01C and EIGEN-CG03C new global models up to degree and order 360.
These preliminary tests lead to the choice of SRTM DTM (integrated in no-data holes), with an added bathymetry derived by the Italian 1:25 000 official cartography near the coasts and the NOAA bathymetry in high seas. The validation of the new global models and the comparison with EGM96 model show that, in terms of geoid computation, the EGM96 yields better results. Moreover, the validation of new available land gravity data and the cross-validation of two sets of gravity data on sea have been completed.     

The reduction of aliasing in gravity anomalies and geoid heights using digital terrain data

Observations of gravity can be aliased by virtue of the logistics involved in collecting these data in the field. For instance, gravity measurements are often made in more accessible lowland areas where there are roads and tracks, thus omitting areas of higher relief in between. The gravimetric determination of the geoid requires mean terrain-corrected free-air anomalies; however, anomalies based only on the observations in lowland regions are not necessarily representative of the true mean value over the topography. A five-stage approach is taken that uses a digital elevation model, which provides a more accurate representation of the topography than the gravity observation elevations, to reduce the unrepresentative sampling in the gravity observations. When using this approach with the Australian digital elevation model, the terrain-corrected free-air anomalies generated from the Australian gravity data base change by between 77.075 and −84.335 mgal (−0.193 mgal mean and 2.687 mgal standard deviation). Subsequent gravimetric geoid computations are used to illustrate the effect of aliasing in the Australian gravity data upon the geoid. The difference between 'aliased' and 'non-aliased' gravimetric geoid solutions varies by between 0.732 and −1.816 m (−0.058 m mean and 0.122 m standard deviation). Based on these conceptual arguments and numerical results, it is recommended that supplementary digital elevation information be included during the estimation of mean gravity anomalies prior to the computation of a gravimetric geoid model.     

Instantaneous deformation from continuous GPS: contributions from quasi-periodic loads

Continuous GPS (CGPS) coordinate time-series are known to experience repeating deformation signals with seasonal and other periods. It is unlikely that these signals represent perfect sinusoids with temporally constant amplitude. We develop an analysis method that accommodates temporal variations in the amplitudes of sinusoidal signals. We apply the method to simulated coordinate time-series to numerically explore the potential consequences of neglecting decadal variation in amplitude of annual motions on the residual-error spectra of CGPS measurements, as well as potential bias in estimates for secular site velocity. We find that secular velocity bias can be appreciable for shorter time-series, and that residual-error time-series of longer duration may contain significant power in a broad band centred on semi-annual frequency if temporal variation in the amplitude of annual motions is not accounted for in the model used to reduce the observations to residuals. It may be difficult to differentiate the bandpass filtered signature of mismodelled loading signals from power-law noise, using residual-error spectra for shorter time-series. We provide an example application to a ∼9-yr coordinate time-series for a CGPS station located in southern California at Carbon Creek Control Structure (CCCS), which is known to experience large amplitude seasonal motions associated with the Santa Ana aquifer system.