🔭 astrophysics

Perturbation-Aware Admissible Regions: First-Order J2 Corrections to Too-Short-Arc Boundary Construction

This paper derives, validates, and scopes the domain of first-order $J_2$ corrections to the energy and perigee boundaries of the admissible region for too-short-arc initial orbit determination, demonstrating that these perturbations significantly improve boundary accuracy in LEO, MEO, and GEO regimes while reclassifying a small but non-negligible fraction of the solution space.

Original authors: Praise Nesvinga

Published 2026-08-21

📖 6 min read🧠 Deep dive

CC BY 4.0

Original authors: Praise Nesvinga

Original paper licensed under CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Every time a telescope spots a piece of space debris or a satellite, it captures a fleeting glimpse: a tiny dot moving across the sky for just a few seconds. From this brief snapshot, astronomers know exactly where the object is pointing and how fast it is moving across the sky, but they do not know how far away it is or how fast it is moving toward or away from them. To predict where that object will go next, they must fill in these missing pieces. For decades, the standard method for doing this has relied on a simplified model of gravity, treating the Earth as a perfect sphere and the object as moving in a smooth, unchanging path. This approach works well enough to generate a list of possible orbits, but it contains a subtle flaw: the initial list is built on an idealized world, while the tools used to track the object later account for the messy reality of Earth's actual shape.

This inconsistency creates a gap between the starting point and the destination. The Earth is not a perfect sphere; it bulges slightly at the equator, a feature known as oblateness. This bulge exerts a gentle, uneven pull on orbiting objects, tugging them in ways that a perfect sphere would not. While astronomers have long known how to calculate these tugs for objects already in orbit, they have not previously adjusted the very first step of the process—the creation of the initial list of possibilities—to account for them. The result is that the starting list of potential orbits is generated using a map that ignores a major geographical feature of the terrain.

A new study by Praise Nesvinga addresses this disconnect by rewriting the rules for how these initial lists are constructed. The researcher developed a way to build the list of possible orbits directly from a model that includes the Earth's bulge, rather than adding the bulge's effects later. The study focuses on two specific boundaries that define the list: the energy limit, which ensures the object stays in orbit rather than flying off into space, and the perigee limit, which ensures the object does not crash into the Earth or its atmosphere. By deriving new mathematical corrections for these boundaries that account for the Earth's oblateness, the study provides a more consistent and physically accurate starting point for tracking space objects.

The findings reveal that the correction for the energy limit is straightforward. Because the Earth's bulge depends only on where an object is located and not on how fast it is moving, the adjustment simply shifts the energy boundary slightly. This shift is small but consistent, and it remains valid across almost the entire range of possible orbits, except right at the very edges where the calculations become unstable. The correction for the perigee limit is more complex. It requires understanding how the Earth's bulge causes the shape of the orbit to wiggle slightly over time. The study shows that these wiggles depend on the object's speed and the tilt of its orbit. For objects with a moderate amount of elliptical shape in their path, the correction is reliable and can be calculated with high precision.

However, the study also identifies a clear limit to where this new method works. The mathematical approach relies on the object's orbit being sufficiently stretched out. When an object's orbit is nearly circular, the wiggles caused by the Earth's bulge become comparable in size to the orbit itself, breaking the assumptions used in the calculation. In these cases, the new correction becomes inaccurate. The researchers determined that the method is trustworthy for orbits where the elliptical shape is larger than a specific threshold, which corresponds to a value of about 0.01. For orbits smaller than this, a different mathematical approach would be needed, a task the author notes is left for future work.

Within the range where the method is valid, the results offer a practical safety margin. For objects in low Earth orbit, the correction suggests that the safe distance from the Earth's atmosphere should be increased by roughly 10 to 12 kilometers to account for the Earth's bulge. This margin is not a guess; it is a calculated buffer that ensures the object will not dip too low during its orbit. As the orbit gets higher, this safety margin shrinks, dropping to just a few kilometers for objects in geostationary orbit, because the Earth's bulge has less influence the farther away one goes. The study also examined how the orientation of the orbit affects this margin, finding that the worst-case scenario varies by a factor of about 1.7 depending on the angle, but the standard calculation remains a safe, conservative estimate.

To test these ideas, the researcher ran detailed computer simulations that tracked objects under the influence of the Earth's bulge and compared the results against the new formulas. The simulations confirmed that the new corrections match the reality of the motion very closely, with errors of only a few percent for the valid range of orbits. The study also checked whether ignoring other, smaller gravitational effects would change the outcome, finding that the Earth's bulge is the dominant factor and that higher-order gravitational details do not significantly alter the results. When applied to a synthetic example of a real-world observation, the new method shifted the boundary of the possible orbits by a tiny amount, reclassifying less than one hundredth of one percent of the total area. While this seems small, it represents a fundamental improvement in consistency: the list of possibilities is now generated using the same physical rules that govern the object's actual motion.

This work does not claim to solve every problem in tracking space debris, nor does it suggest that the old methods were useless. Instead, it removes a long-standing inconsistency in the process. By embedding the reality of the Earth's shape directly into the initial construction of the orbit list, the study ensures that the starting point and the subsequent tracking are aligned. The result is a more robust foundation for predicting where space objects will go, offering a clearer, more reliable picture of the traffic in our skies without requiring massive increases in computing power. The study provides a blueprint for how to handle other forces, such as air resistance or the gravity of the moon, in the same way, pointing the way toward a future where the initial guess is as physically grounded as the final prediction.

Technical Summary: Perturbation-Aware Admissible Regions

Problem Statement
The admissible region (AR) method is a standard approach for initial orbit determination of "too-short-arc" (TSA) tracklets, where angular data (right ascension, declination, and rates) are known, but range ( $\rho$ ) and range-rate ( $\dot{\rho}$ ) are unconstrained. Current AR formulations define the set of dynamically consistent $(\rho, \dot{\rho})$ pairs using two boundary curves: an energy constraint (ensuring a bound orbit) and a perigee constraint (ensuring the object does not intersect the Earth). However, these boundaries are constructed exclusively using unperturbed two-body (Keplerian) dynamics. This creates a foundational inconsistency: candidate orbits are generated from an idealized model, while subsequent propagation and refinement routinely employ perturbed dynamics (e.g., Earth's oblateness, drag). No existing formulation embeds perturbations directly into the boundary definitions themselves; perturbations are typically applied only downstream during orbit refinement.

Methodology
This paper derives, verifies, and scopes the domain of validity for first-order corrections to both AR boundaries arising from the Earth's $J_2$ oblateness perturbation.

Energy-Constraint Correction:
The author derives that the $J_2$ potential adds a purely position-dependent term to the specific orbital energy. Since this term depends only on $\rho$ and not $\dot{\rho}$ , it modifies the $C(\rho)$ component of the energy equation without altering the $\dot{\rho}^2$ or $\dot{\rho}$ terms. The correction is expressed as a shift to the boundary curve, derived via a first-order Taylor expansion of the square-root discriminant. The author verifies this symbolically using SymPy and numerically by comparing the expanded solution against exact $J_2$ -perturbed integration. They identify that the correction diverges near the AR edge where the unperturbed discriminant approaches zero, scaling as the inverse 3/2 power of the discriminant.
Perigee-Constraint Correction:
Correcting the perigee constraint is structurally more complex because the osculating perigee radius ( $r_p = a(1-e)$ ) depends on short-period oscillations of the semi-major axis ( $a$ ) and eccentricity ( $e$ ) driven by $J_2$ .
- The author derives the short-period response of $a$ and $e$ to $J_2$ using Lagrange planetary equations, expanding to first order in eccentricity ( $e$ ).
- They obtain closed-form expressions for the oscillations $\delta a(u)$ and $\delta e(u)$ as functions of the argument of latitude $u$ .
- These are combined to form a closed-form expression for the perigee radius perturbation $\delta r_p(u)$ .
- The derivation is extended to a general argument of perigee ( $\omega$ ), resulting in a five-term trigonometric polynomial. The author demonstrates that for general $\omega$ , finding the worst-case extremum requires solving a degree-six polynomial (via Weierstrass substitution), whereas the $\omega=0$ case reduces to a cubic.
Validation:
The analytical results are validated against direct numerical integration of the $J_2$ -perturbed equations of motion (RK45) across Low Earth Orbit (LEO), Medium Earth Orbit (MEO), and Geostationary Orbit (GEO) regimes. The study also validates the results against a zonal gravity field extending through $J_6$ to ensure higher-order harmonics do not significantly alter the residual error.

Key Results

Energy Boundary: The correction is a position-dependent shift valid throughout the AR except near the boundary edge. Numerical tests confirm second-order convergence of the residual as the perturbation strength is scaled.
Perigee Boundary Validity: The first-order correction is accurate to within a few percent for eccentricities satisfying $e \gtrsim 10 J_2 (R_\oplus/a)^2$ . Below this threshold (approximately $e \gtrsim 9 \times 10^{-3}$ for LEO), the correction breaks down because the $J_2$ -forced oscillation in eccentricity becomes comparable to the nominal eccentricity, violating the perturbative assumption.
Safety Margins: Within the valid domain, the correction yields a near-constant safety margin for the perigee constraint. At LEO altitudes, this margin is approximately 10–12 km. The margin scales as $J_2 R_\oplus^2 / a$ at higher altitudes (e.g., ~2–3 km at GEO).
Dependence on $\omega$ : The worst-case safety margin varies by a factor of $\sim 1.7$ depending on the argument of perigee $\omega$ . The $\omega=0$ case provides the maximum margin, making it a conservative bound for general application.
Impact on Admissible Region: In a synthetic example, applying the correction shifted the perigee boundary by $0.0007–0.0047$ km/s in $\dot{\rho}$ . This resulted in a reclassification of only 0.0029% of the admissible region, a fraction too small to be resolved by typical coarse AR sampling grids (which often have spacings of 0.01–0.1 km/s).
Higher-Order Harmonics: Validation against $J_2$ through $J_6$ showed that the residual discrepancy between the analytical correction and numerical integration is dominated by the $O(e^2)$ truncation error and the small-eccentricity floor mechanism, not by neglected higher-order zonal harmonics.

Significance and Claims
The paper claims to remove a long-standing modeling inconsistency in TSA orbit determination by embedding perturbation theory directly into the admissible-region boundary construction. The primary contribution is not necessarily the reclassification of existing coarse-grid samples (which is minimal for the tested geometry), but the establishment of a methodology where candidate orbits are generated from a dynamical model consistent with the one used for their subsequent propagation.

The author positions this work as a foundational extension of the AR methodology, demonstrating that boundaries can be constructed from perturbed dynamics rather than exclusively Keplerian dynamics. They provide a closed-form, computationally negligible safety margin for the perigee constraint that improves physical consistency. The paper explicitly notes that for eccentricities below the derived threshold, a different formalism (likely using quasi-nonsingular elements) is required, which is left for future work. The framework is presented as generalizable to other perturbations, including atmospheric drag, third-body effects, and solar radiation pressure.

More like this