We propose a methodology for local gravity field modelling from gravity data using spherical radial basis functions. The methodology comprises two steps: in step 1, gravity data (gravity anomalies and/or gravity disturbances) are used to estimate the disturbing potential using least-squares techniques. The latter is represented as a linear combination of spherical radial basis functions (SRBFs). A data-adaptive strategy is used to select the optimal number, location, and depths of the SRBFs using generalized cross validation. Variance component estimation is used to determine the optimal regularization parameter and to properly weight the different data sets. In the second step, the gravimetric height anomalies are combined with observed differences between global positioning system (GPS) ellipsoidal heights and normal heights. The data combination is written as the solution of a Cauchy boundary-value problem for the Laplace equation. This allows removal of the non-uniqueness of the problem of local gravity field modelling from terrestrial gravity data. At the same time, existing systematic distortions in the gravimetric and geometric height anomalies are also absorbed into the combination. The approach is used to compute a height reference surface for the Netherlands. The solution is compared with NLGEO2004, the official Dutch height reference surface, which has been computed using the same data but a Stokes-based approach with kernel modification and a geometric six-parameter "corrector surface" to fit the gravimetric solution to the GPS-levelling points. A direct comparison of both height reference surfaces shows an RMS difference of 0.6 cm; the maximum difference is 2.1 cm. A test at R. Klees (B) · R. Tenzer · I. Prutkin · T. Wittwer

We consider the problem of local (quasi-)geoid modelling from terrestrial gravity anomalies. Whereas this problem is uniquely solvable (up to spherical harmonic degree one) if gravity anomalies are globally available, the problem is non-unique if gravity anomalies are only available within a local area, which is the typical situation in local/regional gravity field modelling. We derive a mathematical description of the kernel of the gravity anomaly operator. The non-uniqueness can be removed using external height anomaly information, e.g., provided by GPS-levelling. The corresponding problem is formulated as a Cauchy problem for the Laplace equation. The existence and uniqueness of the solution of the Cauchy problem is guaranteed by the Cauchy-Kowalevskaya theorem. We propose several numerical procedures to compute the solution of the Cauchy problem from given differences between gravimetric and geometric height anomalies. We apply the numerical techniques to real data over the Netherlands and Germany. We show that we can compute a unique quasi-geoid from observed gravimetric and geometric height anomalies, which agree with the data within the expected noise level. We conclude that observed differences between gravimetric height anomalies and geometric height anomalies derived from GPS and levelling cannot only be attributed to systematic errors in the data sets, but are also caused by the intrinsic non-uniqueness of the problem of local quasi-geoid modelling from gravity anomalies. Hence, GPS-levelling data are necessary to get a unique solution, which also implies that they should not be I. Prutkin · R. Klees (B) Delft Institute of Earth Observation and Space Systems (DEOS), Delft University of Technology, Kluyverweg 1, 2629 HS Delft, The Netherlands e-mail: r.klees@tudelft.nl used to validate local quasi-geoid solutions computed on the basis of gravity anomalies.

The second international comparison of absolute gravimeters was held in Walferdange, Grand Duchy of Luxembourg, in November 2007, in which twenty absolute gravimeters took part. A short description of the data processing and adjustments will be presented here and will be followed by the presentation of the results. Two different methods were applied to estimate the relative offsets between the gravimeters. We show that the results are equivalent as the uncertainties of both adjustments overlap. The absolute gravity meters agree with one another with a standard deviation of 2 μgal (1 gal = 1 cm/s 2 ). In 1999, a laboratory ( Fig. 5.1) dedicated to the comparison of absolute gravimeters was built within the WULG. The laboratory lies 100 m below the surface at a distance of 300 m from the entrance of the mine. The WULG is environmentally stable (i.e., constant temperature and humidity within the lab), and is extremely well isolated from anthropogenic noise. It has the power and space requirements to be able to accommodate up 16 instruments operating simultaneously. IntroductionMultiple absolute gravimeter comparisons are regularly carried out. Being absolute instruments, these gravimeters cannot really be calibrated. Only some of their components (such as the atomic clock and the laser) can be calibrated by comparison with known standards. The only way one currently has to verify their good working order is via a simultaneous comparison with other absolute gravimeters of the same and/or if possible even of a different model, to detect possible systematic errors.During a comparison, we cannot estimate how accurate the meters are: in fact, as we have no way to know the true value of g, we can only investigate the relative offsets between instruments. This means that all instruments can suffer from the same unknown and undetectable systematic error. However, differences larger than the uncertainty of the measurements, is usually indicative of a possible systematic error.For the second comparison in Walferdange, a few new procedures have been introduced. First, some of the participants accepted to take part in a

No abstract

Please cite this article as: Prutkin, I., Saleh, A., Gravity and magnetic data inversion for 3D topography of the Moho discontinuity in the northern Red Sea area, Egypt, Journal of Geodynamics (2008Geodynamics ( ), doi:10.1016Geodynamics ( /j.jog.2008 This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain. Page 1 of 19A c c e p t e d M a n u s c r i p t AbstractThe main goal of our study is to investigate 3D topography of the Moho boundary for the area of the northern Red Sea including Gulf of Suez and Gulf of Aqaba. For potential field data inversion we apply a new method of local corrections. The method is efficient and does not require trial-and-error forward modeling. To separate sources of gravity and magnetic field in depth, a method is suggested, based on upward and downward continuation. Both new methods are applied to isolate the contribution of the Moho interface to the total field and to find its 3D topography. At the first stage, we separate near-surface and deeper sources. According to the obtained field of shallow sources a model of the horizontal layer above the depth of 7 km is suggested, which includes a density interface between light sediments and crystalline basement. Its depressions and uplifts correspond to known geological structures. At the next stage, we isolate the effect of very deep sources (below 100 km) and sources outside the area of investigation. After subtracting this field from the total effect of deeper sources, we obtain the contribution of the Moho interface. We make inversion separately for the area of rifts (Red Sea, Gulf of Suez and Gulf of Aqaba) and for the rest of the area. In the rift area we look for the upper boundary of low-density, heated anomalous upper mantle. In the rest of the area the field is satisfied by means of topography for the interface between lower crust and normal upper mantle. Both algorithms are applied also to the magnetic field. The magnetic model of the Moho boundary is in agreement with the gravitational one.