Relativistic positioning: four-dimensional numerical approach in Minkowski space-time

We simulate the satellite constellations of two Global Navigation Satellite Systems: Galileo (EU) and GPS (USA). Satellite motions are described in the Schwarzschild space-time produced by an idealized spherically symmetric non rotating Earth. The trajectories are then circumferences centered at the same point as Earth. Photon motions are described in Minkowski space-time, where there is a well known relation, Coll, Ferrando & Morales-Lladosa (2010), between the emission and inertial coordinates of any event. Here, this relation is implemented in a numerical code, which is tested and applied. The first application is a detailed numerical four-dimensional analysis of the so-called emission coordinate region and co-region. In a second application, a GPS (Galileo) satellite is considered as the receiver and its emission coordinates are given by four Galileo (GPS) satellites. The bifurcation problem (double localization) in the positioning of the receiver satellite is then pointed out and discussed in detail.Comment: 16 pages, 9 figures, published (online) in Astrophys. Space Sc

Closed-Form transformation between geodetic and ellipsoidal coordinates

We present formulas for direct closed-form transformation between geodetic coordinates(Φ, λ, h) and ellipsoidal coordinates (β, λ, u) for any oblate ellipsoid of revolution.These will be useful for those dealing with ellipsoidal representations of the Earth's gravityfield or other oblate ellipsoidal figures. The numerical stability of the transformations for nearpolarand near-equatorial regions is also considered

Ellipsoidal area mean gravity anomalies - precise computation of gravity anomaly reference fields for remove-compute-restore geoid determination

Gravity anomaly reference fields, required e.g. in remove-compute-restore (RCR) geoid computation, are obtained from global geopotential models (GGM) through harmonic synthesis. Usually, the gravity anomalies are computed as point values or area mean values in spherical approximation, or point values in ellipsoidal approximation. The present study proposes a method for computation of area mean gravity anomalies in ellipsoidal approximation ('ellipsoidal area means') by applying a simple ellipsoidal correction to area means in spherical approximation. Ellipsoidal area means offer better consistency with GGM quasi/geoid heights. The method is numerically validated with ellipsoidal area mean gravity derived from very fine grids of gravity point values in ellipsoidal approximation. Signal strengths of (i) the ellipsoidal effect (i.e., difference ellipsoidal vs. spherical approximation), (ii) the area mean effect (i.e., difference area mean vs. point gravity) and (iii) the ellipsoidal area mean effect (i.e., differences between ellipsoidal area means and point gravity in spherical approximation) are investigated in test areas in New Zealand and the Himalaya mountains. The impact of both the area mean and the ellipsoidal effect on quasigeoid heights is in the order of several centimetres. The proposed new gravity data type not only allows more accurate RCR-based geoid computation, but may also be of some value for the GGM validation using terrestrial gravity anomalies that are available as area mean values

The zero gravity curve and surface and radii for geostationary and geosynchronous satellite orbits

A geosynchronous satellite orbits the Earth along a constant longitude. A special case is the geostationary satellite that is located at a constant position above the equator. The ideal position of a geostationary satellite is at the level of zero gravity, i.e. at the geocentric radius where the gravitational force of the Earth equals the centrifugal force. These forces must be compensated for several perturbing forces, in particular for the lunisolar tides. Considering that the gravity field of the Earth varies not only radially but also laterally, this study focuses on the variations of zero gravity not only on the equator (for geostationary satellites) but also for various latitudes. It is found that the radius of a geostationary satellite deviates from its mean value of 42164.2 km only within ±2 m, mainly due to the spherical harmonic coefficient J22, which is related with the equatorial flattening of the Earth. Away from the equator the zero gravity surface deviates from the ideal radius of a geosynchronous satellite, and more so for higher latitudes. While the radius of the former surface increases towards infinity towards the poles, the latter decreases about 520 m from the equator to the pole. Tidal effects vary these radii within ±2.3 km

Groebner Basis in Geodesy and Geoinformatics

In geodesy and geoinformatics, most problems are nonlinear in nature and often require the solution of systems of polynomial equations. Before 2002, solutions of such systems of polynomial equations, especially of higher degree remained a bottleneck, with iterative solutions being the preferred approach. With the entry of Groebner basis as algebraic solution to nonlinear systems of equations in geodesy and geoinformatics in the pioneering work “Gröbner bases, multipolynomial resultants and the Gauss Jacobi combinatorial algorithms : adjustment of nonlinear GPS/LPS observations", the playing field changed. Most of the hitherto unsolved nonlinear problems, e.g., coordinate transformation problems, global navigation satellite systems (GNSS)'s pseudoranges, resection-intersection problems in photogrammetry, and most recently, plane fitting in point clouds in laser scanning have been solved. A comprehensive overview of such applications are captured in the first and second editions of our book Algebraic Geodesy and Geoinformatics published by Springer. In the coming third edition, an updated summary of the newest techniques and methods of combination of Groebner basis with symbolic as well as numeric methods will be treated. To quench the appetite of the reader, this presentation considers an illustrative example of a two-dimension coordinate transformation problem solved through the combination of symbolic regression and Groebner basis