Periodic orbits around areostationary points in the Martian gravity field

This study investigates the problem of areostationary orbits around Mars in the three-dimensional space. Areostationary orbits are expected to be used to establish a future telecommunication network for the exploration of Mars. However, no artificial satellites have been placed in these orbits thus far. In this paper, the characteristics of the Martian gravity field are presented, and areostationary points and their linear stability are calculated. By taking linearized solutions in the planar case as the initial guesses and utilizing the Levenberg-Marquardt method, families of periodic orbits around areostationary points are shown to exist. Short-period orbits and long-period orbits are found around linearly stable areostationary points, and only short-period orbits are found around unstable areostationary points. Vertical periodic orbits around both linearly stable and unstable areostationary points are also examined. Satellites in these periodic orbits could depart from areostationary points by a few degrees in longitude, which would facilitate observation of the Martian topography. Based on the eigenvalues of the monodromy matrix, the evolution of the stability index of periodic orbits is determined. Finally, heteroclinic orbits connecting the two unstable areostationary points are found, providing the possibility for orbital transfer with minimal energy consumption.

💡 Research Summary

This paper presents a comprehensive dynamical study of areostationary orbits—circular, equatorial orbits that remain fixed over a point on the rotating surface of Mars. Using the latest Mars gravity field model (e.g., MGS‑GRS‑1200), the authors expand the planetary potential in spherical harmonics up to the required degree and order, and formulate the effective potential in a rotating reference frame. By solving the equilibrium condition ∇U = 0, four areostationary points (E₁–E₄) are located, two on the longitudinal axis (E₁, E₃) and two on the transverse axis (E₂, E₄).

Linearization of the equations of motion about each equilibrium yields a state‑transition matrix whose eigenvalues determine local stability. E₁ and E₃ possess only purely imaginary eigenvalue pairs, indicating linear stability, whereas E₂ and E₄ exhibit a real eigenvalue pair, making them linearly unstable. This classification reflects the asymmetry of Mars’s gravity field and the influence of higher‑order zonal and tesseral terms.

To explore the nonlinear regime, the authors first obtain planar periodic solutions as initial guesses and then apply the Levenberg‑Marquardt (LM) algorithm to enforce the full three‑dimensional periodicity condition (state vector repeats after a period T). The LM method simultaneously minimizes the deviation from periodicity and satisfies any imposed geometric constraints (e.g., limiting latitude excursions).

The numerical continuation reveals three families of periodic orbits around each equilibrium:

1. Short‑period orbits (periods of a few tens of minutes to a couple of hours). These exist around both stable and unstable points and keep the satellite within a few degrees of longitude from the nominal areostationary longitude, making them attractive for continuous communication coverage with modest station‑keeping.

2. Long‑period orbits (periods of several days). They are found only around the linearly stable points (E₁, E₃) and involve larger longitudinal excursions (tens of degrees). Their existence is linked to a nonlinear resonance between the natural frequencies of the linearized system.

3. Vertical periodic orbits that oscillate primarily in the out‑of‑plane direction while maintaining a near‑constant longitude. Such orbits are possible around all four points and could be used to vary the sub‑satellite latitude without large fuel penalties.

For each periodic solution the monodromy matrix (state transition over one period) is computed, and its eigenvalues are used to define a stability index (the product of the absolute values of the eigenvalues). Index values close to unity indicate marginal stability; deviations signal exponential divergence. The authors show that short‑period orbits around stable points retain an index ≈ 1, whereas long‑period orbits display index growth in certain parameter ranges, implying limited long‑term stability.

A particularly noteworthy result is the identification of heteroclinic connections between the two unstable equilibria (E₂ ↔ E₄). By integrating trajectories along the unstable manifolds of one point and the stable manifolds of the other, the authors construct low‑energy transfer paths that require minimal thrust. Such heteroclinic orbits could serve as “highways” for moving a spacecraft between different areostationary longitudes with negligible propellant consumption.

The paper concludes that a rich set of periodic dynamics exists around Mars’s areostationary points, and that these dynamics can be harnessed for practical mission design: (i) short‑period orbits for quasi‑stationary communication platforms, (ii) long‑period orbits for periodic coverage of a broader surface area, (iii) vertical orbits for latitude‑adjustable observations, and (iv) heteroclinic transfers for fuel‑efficient repositioning. The authors’ methodology—combining high‑fidelity gravity modeling, linear stability analysis, and robust nonlinear continuation—provides a solid foundation for future Mars telecommunication constellations and scientific observation campaigns.