[Paper Review] Relativistic Positioning Systems: Numerical Simulations

N. Puchades, D. Sáez|arXiv (Cornell University)|Apr 3, 2014

Geophysics and Gravity Measurements8 references3 citations

TL;DR

This paper numerically simulates relativistic positioning systems using GPS and Galileo satellite constellations in Schwarzschild spacetime, modeling satellite orbits and emission coordinates to identify regions of bifurcation (double positioning) and low Jacobian magnitude, which correlate with large positioning errors. The key contribution is the identification of orbital altitudes above ~10⁴ km where bifurcation and high errors occur, with recommendations to select 4-satellite configurations maximizing |J| to minimize errors for deep-space users.

ABSTRACT

The motion of satellite constellations similar to GPS and Galileo is numerically simulated and, then, the region where bifurcation (double positioning) occurs is appropriately represented. In the cases of double positioning, the true location may be found using additional information (angles or times). The zone where the Jacobian, J, of the transformation from inertial to emission coordinates vanishes is also represented and interpreted. It is shown that the uncertainties in the satellite world lines produce positioning errors, which depend on the value of |J|. The smaller this quantity the greater the expected positioning errors. Among all the available 4-tuples of satellites, the most appropriate one -for a given location- should minimize positioning errors (large enough |J| values) avoiding bifurcation. Our study is particularly important to locate objects which are far away from Earth, e.g., satellites.

Motivation & Objective

• To investigate the occurrence of bifurcation (two possible user positions) in relativistic positioning systems for users far from Earth.
• To analyze positioning errors arising from uncertainties in satellite world lines, particularly near regions where the Jacobian |J| of the inertial-to-emission coordinate transformation vanishes.
• To identify and represent regions in spacetime where positioning becomes ambiguous or inaccurate due to low |J| or bifurcation.
• To guide the selection of optimal 4-satellite configurations that minimize positioning errors by maximizing |J| and avoiding bifurcation.

Proposed method

• Numerical simulation of satellite constellations (GPS and Galileo) in Schwarzschild spacetime with circular orbits at altitudes ~20,000–23,000 km.
• Use of the Newton-Raphson method to solve the system of four equations for emission coordinates (τ¹, τ², τ³, τ⁴) from known inertial coordinates.
• Employment of an analytical transformation from emission to inertial coordinates (from Ref. [4]) to avoid numerical inversion and enable systematic study of bifurcation and Jacobian behavior.
• Computation of the Jacobian determinant J of the transformation from inertial to emission coordinates to identify regions of vanishing J.
• Representation of bifurcation zones and low-|J| regions via cross-sections of the 4D emission region using color-coded time-ordered points.
• Use of multiple 1800-point segments along user world lines to track transitions between single and bifurcated positioning, with color sequences indicating evolving solutions.

Experimental results

Research questions

• RQ1 At what orbital altitudes do bifurcation regions emerge in relativistic positioning systems for distant users?
• RQ2 How do uncertainties in satellite world lines affect positioning accuracy, and what role does the Jacobian |J| play in this error amplification?
• RQ3 Which 4-satellite configurations minimize positioning errors by maximizing |J| and avoiding bifurcation for a given user location?
• RQ4 How do the dynamics of satellite motion influence the transition between single-positioning and bifurcated positioning states?
• RQ5 What is the spatial extent and structure of regions where the Jacobian J vanishes, leading to singularities in the positioning transformation?

Key findings

• Bifurcation occurs in regions where the Jacobian determinant J of the inertial-to-emission coordinate transformation vanishes, leading to two possible user positions.
• Users at altitudes above approximately 10⁴ km may cross zones with bifurcation, experiencing transitions between single and double positioning.
• Positioning errors increase significantly in regions where |J| is small, due to amplification of satellite world line uncertainties.
• The most accurate positioning is achieved by selecting a 4-satellite configuration that maximizes |J| and avoids bifurcation, especially for deep-space users.
• For users on Earth’s surface, positioning remains single-valued (χ² ≤ 0) and |J| does not vanish, confirming the robustness of current GNSS for terrestrial use.
• The simulation results show that satellite motion causes continuous transitions between bifurcation states, with one solution tending to infinity at transition boundaries.