Pith. sign in

REVIEW 3 major objections 4 minor 71 references

A dynamical coherency gate for state recovery: Statistical requiem for the long arc of cislunar orbital mis-prediction

T0 review · 3 major / 4 minor · reviewed 2026-08-06 · deepseek-v4-flash

Pith's one-line read A statistical coherency gate on historical TLE data recovers OGO-1's orbit well enough to predict its 2020 reentry within five hours.

desk verdict A plausible and well-written TLE recovery method whose headline 18-year reentry match is in-sample — the 2002 state is built from TLEs spanning through 2020 — but the UKF windows offer partial out-of-sample support and the paper deserves a serious referee. read the letter →

arxiv 2506.22748 v1 pith:3T2QS7PO submitted 2025-06-28 astro-ph.EP

classification astro-ph.EP
keywords CislunarspaceDynamicalevolutionandstabilityGaussianmixturemodeldebrissituationalawarenessUnscentedKalmanfilterTwo-lineelementsetsOrbitdetermination
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper claims that a statistically filtered batch of historical two-line element sets (TLEs) can stand in for precise orbit determination for a satellite in the dynamically sensitive cislunar region. Applying the method to OGO-1, a 1960s science satellite that was tracked intermittently and untracked for roughly thirty years, the authors recover an initial state at a 2002 epoch and show that propagating it reproduces the satellite's full multi-decade trajectory, including its August 2020 atmospheric reentry within about five hours. The practical significance is that accurate long-arc reconstruction and reentry forecasting for high-altitude debris and legacy missions could be done from public TLE data alone, without the original tracking and covariance records.

What carries the argument

The load-bearing mechanism is the coherency gate: filtering in mean-element space rather than osculating-element space. Mean elements are obtained by FFT-based numerical averaging of short propagated arcs, which removes short-period variations so that inter-TLE scatter reflects genuine dynamical consistency rather than epoch artifacts. A Gaussian mixture model, with median-absolute-deviation or quantile fallback, selects the dominant statistical core of mean elements; the consensus osculating state is then rebuilt from the inlier osculating elements. This gated state seeds an unscented Kalman filter for intermediate-epoch estimation, and all long-arc propagations are done with a high-fidelity Cowell-type integrator that handles the strong lunar and solar perturbations of the cislunar regime.

What would settle it

Run the coherency-gated recovery on the OGO-1 TLEs with all TLEs after 2007 removed, or with the gate epoch moved to a different date, and check whether the predicted reentry still lands within five hours of MJD 59090.86389. Because the paper identifies two distinct reentry-prediction clusters in the post-2007 data that may reflect systematic generation changes, a reentry match that depends on including one of those clusters would show the gate is selecting a biased mode rather than a dynamically coherent core.

Watch

Extended reading notes

Core claim

The central discovery is that the coherent statistical core of a TLE ensemble, once mapped into mean-element space, carries enough dynamical information to define a high-fidelity initial condition. For each TLE, the paper evaluates SGP4 at its native epoch, propagates the resulting Cartesian state to a common reference epoch with an ephemeris-quality integrator, converts to osculating elements, and applies FFT-based numerical averaging to obtain mean elements. Outlier detection with a Gaussian mixture model and median-absolute-deviation fallback isolates the dynamically coherent inlier subset; the element-wise mean or median of the inlier osculating elements then forms the recovered state. At the gate epoch MJD 52461.31528 the two independent filter runs produce states that, propagated forward, yield reentry at MJD 59090.75562 and 59090.67808, within roughly five hours of the true decay at MJD 59090.86389. The same recovered state also captures the long-period von Zeipel-Lidov-Kozai oscillations in eccentricity and inclination over the full arc, and the method is extended to a windowed unscented Kalman filter setting that produces consistent state estimates at intermediate epochs from sparse, irregular TLEs.

Load-bearing premise

After SGP4 evaluation and numerical averaging, the dominant cluster of mean elements in the TLE batch is an unbiased estimate of the true orbital state at the gate epoch; if systematic TLE-generation biases shift that cluster, the recovered state and the close reentry match would be an artifact of the data reduction rather than an independent physical prediction.

Editorial extensions

If this is right

  • Long-arc trajectory reconstruction for objects without precise ephemerides becomes possible from public TLE archives, provided enough TLEs span the object's dynamical evolution.
  • Reentry forecasts for high-eccentricity cislunar objects can improve from spreads of nearly a year, as seen with raw TLE propagation of OGO-1, to agreement within hours of the observed decay.
  • The recovered initial condition reproduces not just the decay date but the full secular evolution, including von Zeipel-Lidov-Kozai-driven oscillations in eccentricity and inclination, so the method can serve as a dynamical-consistency check on historical catalogs.
  • The windowed unscented Kalman filter variant yields consistent state estimates and covariances at epochs not directly represented in the TLE record, bridging sparse observation gaps.
  • Operationally, the same gating logic could flag and down-weight inconsistent TLEs in catalog and conjunction-analysis pipelines without requiring new tracking data.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • A natural test, not pursued in the paper, is to shift the gate epoch or truncate the TLE batch before and after the two post-2007 clusters the paper identifies as systematically biased; if the five-hour reentry agreement persists only when the biased cluster is included, the match is a property of the data reduction rather than of the recovered physics.
  • The method implies that the long-period dynamical content of even a sparse, decades-long TLE record is dominated by a single coherent secular mode; if true, similar gates should recover accurate decay epochs for other long-lived high-eccentricity objects with known reentry dates, turning historical decay records into validation benchmarks.
  • Because the gate operates on mean elements, it could be combined with resonance analysis to decide whether a TLE batch spans a single secular regime or multiple regimes separated by a close lunar encounter or resonance crossing; the paper's OGO-1 arc appears to stay in one such regime, and objects that transition regimes would challenge the element-wise median-of-inliers reconstruction.
  • The approach treats TLE noise as a statistical ensemble property rather than as per-object measurement error; a practical consequence left implicit by the authors is that archives of old TLEs could be reprocessed in bulk to produce candidate initial conditions for conjunction screening of the entire cislunar population.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper proposes a 'coherency-gated filtering' method to recover a statistically meaningful initial condition from historical TLEs for highly eccentric, resonance-dominated cislunar objects. The method evaluates each TLE with SGP4, propagates all states to a common epoch with a high-fidelity integrator (ASSIST/IAS15), converts to osculating elements, forms mean elements by FFT-based numerical averaging, and then applies MAD or GMM-based outlier rejection with fallbacks to select a dynamically coherent subset. The consensus state is then used as an initial condition for forward propagation or as the seed for a UKF. The central demonstration is OGO-1: a state recovered at MJD 52461.31528 (6 July 2002) is propagated with ASSIST and reported to reproduce the satellite's August 2020 reentry within about five hours, and windowed UKF experiments around 2010-2016 give reentry predictions within about 1.5-2 days. The paper argues that filtered TLE statistics can serve as a proxy for precise orbit determination in regimes where resonances and long-period perturbations dominate.

Significance. If the central claim were validated out-of-sample, the method would be a useful contribution to space situational awareness for xGEO objects, where public TLE data are often the only historical record. The paper has real strengths: it uses a transparent, reproducible algorithmic pipeline; it propagates with a modern, validated high-fidelity integrator; it filters in mean-element space rather than on raw osculating residuals; and it explicitly targets a genuinely hard dynamical regime (vZLK-driven high-eccentricity motion). The windowed UKF experiments are a partial out-of-sample check because each selected window ends years before the 2020 reentry. However, the headline 18-year reentry 'prediction' is not out-of-sample: the recovered 2002 state is built from the same 2000-2020 TLE set that includes the final reentry arc. The absence of an independent precise ephemeris comparison, the lack of uncertainty quantification, and the presence of several hand-tuned thresholds mean that the current evidence does not establish that the method recovers an unbiased physical state rather than a statistically filtered average that is consistent with the known outcome.

major comments (3)
  1. [§4.1, Algorithm 1, Figure 5 caption] The headline 18-year reentry 'prediction' is not out-of-sample. Algorithm 1 Step 2 propagates every TLE in the input set -- which, according to the Figure 5 caption, spans 23 June 2000 to 20 August 2020 -- backward to the gate epoch MJD 52461.31528. The recovered state therefore encodes information from post-2002 TLEs, including the final reentry arc that is later used as the prediction target. The 5-hour agreement in Table 1 is consequently a consistency check on the ensemble-averaging filter, not an independent prediction. The paper should re-run the gate-epoch recovery using only TLEs with epochs before the gate epoch (or before a training cutoff) and report the resulting recovered state and reentry prediction; alternatively, the current experiment should be explicitly relabeled as a fit/reconstruction rather than a prediction.
  2. [Table 1 and Figure 4] Reentry epoch is too weak a diagnostic to validate the recovered six-dimensional state. The two Table 1 runs differ by about 4.15 degrees in mean anomaly (193.22310 vs 197.37701) yet agree in reentry epoch to about two hours, so the recovered state is highly non-unique with respect to the stated success metric. Moreover, Figure 4 shows two systematic reentry-prediction clusters in the post-2007 data, which the paper attributes to TLE generation or sensor-coverage changes; if the GMM selects one of these biased modes at the gate epoch, the recovered state inherits that bias. The paper should compare the recovered state to an independent precise ephemeris, or to pseudo-observations from a strictly held-out TLE subset, and should report per-element uncertainties or a covariance rather than a single scalar reentry time.
  3. [Algorithm 1 and Section 4 footnote 1] No sensitivity analysis is provided for the free parameters of the method (MAD scale factor, GMM thresholds, percentile fallback, arc mode), and the 'reentry' prediction is based on the criterion that osculating perigee falls below 50 km altitude, not on a full drag-perturbed decay propagation. The 5-hour agreement could in principle be a consequence of tuning the thresholds or of the specific 50 km cutoff rather than of the recovered state being physically correct. The authors should report how the recovered state and predicted reentry epoch vary under reasonable perturbations of the threshold parameters, and should validate the 50 km crossing-time approximation against a drag-included decay simulation for at least one trajectory.
minor comments (4)
  1. [Figure 3] The left-hand schematic box says 'Element-wise medium or mean'; 'medium' should be 'median'.
  2. [Table 2] Table 2 is misaligned: each row appears to contain eight numeric entries while the header has seven columns, and the second MJD-like value (e.g., 55246.20907) is not identified. Please reformat the table and label all columns.
  3. [Abstract and Section 4.1] The abstract states that the recovered state reproduces reentry 'to within one day,' while Table 1 and Section 4.1 claim agreement within five hours; the more precise statement should be used consistently throughout.
  4. [Figure 6] Figure 6 compares the propagation to the 'full TLE time history' including the 1964-1971 era, but the recovered state is from 2002 and the pre-2000 TLEs are not used in the recovery. The caption should clarify that the early TLEs are shown for context only.

Circularity Check

1 steps flagged · score 6.0 of 10

Headline 18-year reentry 'prediction' is in-sample: the 2002 state is reconstructed from the same 2000–2020 TLE batch that includes the final reentry arc, so the match is a smoothing consistency check, not an independent forecast.

  1. fitted input called prediction [Section 4.1 / Figure 5 caption; Algorithm 1 Steps 1–6; Table 1]
    "Application of Algorithm 1 to the approximately 3,200 TLEs of OGO-1 spanning from 23 June 2000 to 20 August 2020, with the coherency gate epoch set to 6 July 2002. (Figure 5 caption); Algorithm 1 Step 2: Propagate to t0: x_i(t0)=P(t_i→t0; x_i)."

    The recovered state at MJD 52461.31528 is constructed by propagating every TLE in the batch — including TLEs from 2002–2020, i.e., after the gate epoch and through the final pre-reentry arc — backward to the gate epoch, filtering the resulting mean elements, and averaging the inlier osculating elements. The subsequent ASSIST forward run and the reported reentry epoch (Table 1) therefore do not test the state against unseen future data; they replay information already contained in the batch used to build the state. The paper itself shows in Figure 4 that raw post-2007 TLE-based predictions already cluster near the true reentry, so selecting a coherent core from the full 2000–2020 dataset and recovering that same reentry is an in-sample consistency check, not an out-of-sample prediction.

full rationale

The central circularity is in-sample validation rather than self-citation. Algorithm 1 takes the full OGO-1 TLE batch spanning 23 June 2000 to 20 August 2020, propagates every TLE — including post-2002 and post-2007 epochs — backward to the 6 July 2002 gate epoch, filters the resulting mean elements, and averages the inliers (Figure 5 caption; Algorithm 1 Steps 1–6). The resulting state is then propagated forward and the 2020 reentry epoch is presented as a prediction (Table 1, Section 4.1). Because the batch already contains the final pre-reentry arc and the paper's own Figure 4 shows raw TLE-based reentry predictions clustering near the true decay, the reported 5-hour agreement is a consistency check, not an independent forecast. The UKF window runs (Section 4.2, Table 2) are less in-sample because their selected windows end before 2020, but the headline claim rests on the 2002 coherency-gate state. No load-bearing self-citation chain was found: the averaging and propagation tools (Schubart, Ely, ASSIST/IAS15, GMAT) are external, and the GMM/MAD machinery is described in the paper itself. The result could be made non-circular by training only on TLEs before the gate epoch or otherwise withholding all data after 2002; as written, the central 'reentry prediction' reduces to an in-sample reconstruction.

Assumptions & free parameters 5 free parameters · 5 assumptions · 0 invented entities

The central result depends on several hand-set parameters (gate epoch, thresholds, UKF noise) and on standard but untested-in-this-paper assumptions about TLE/SGP4 fidelity, FFT averaging, and the ASSIST integrator. No new physical entities are introduced.

free parameters (5)
  • Coherency gate epoch t0 = MJD 52461.31528 (6 July 2002)
    Chosen for the demonstration; the method does not specify how to select it, and the recovered state depends on it.
  • MAD scale factors = e.g., 1.5-sigma, 3-sigma, 1.2-sigma (element-dependent, Figure 5)
    Hand-tuned for each orbital element; not derived from TLE error statistics.
  • GMM configuration and thresholds = modes 'adaptive', 'knee', 'percentile'; mean anomaly de-weighted at 15th percentile
    User-selected fallback logic; no objective criterion given.
  • UKF noise parameters = P0, Q, R user-defined
    Not specified in the paper; needed to reproduce the windowed UKF results.
  • Arc mode = 'long' for main run, 'short' for UKF windows
    Choice affects whether GMM or MAD is used; no rule given.
assumptions (5)
  • domain assumption TLEs, when evaluated with SGP4 at their native epochs, provide osculating states accurate enough for mean-element reconstruction.
    Used in Algorithm 1 Step 1; TLE accuracy in HEO is tens of kilometers, but the method does not incorporate this uncertainty.
  • domain assumption FFT-based numerical averaging over a short arc removes short-period variations and yields mean elements that are comparable across TLEs.
    Step 4; relies on Schubart (1964), Uphoff (1973), Ely (2015).
  • domain assumption The GMM/MAD outlier rejection identifies dynamically incoherent TLEs, and the remaining inlier set is statistically representative of the true state.
    Step 5; assumes the coherent core is the truth and that systematic biases are absent or removable.
  • domain assumption ASSIST/IAS15 provides an ephemeris-quality propagator that captures lunisolar perturbations and drag for OGO-1.
    Used in Steps 2 and 7; validated against GMAT per the text, but no numerical validation shown.
  • domain assumption Reentry occurs when perigee falls below 50 km altitude.
    Footnote 1, based on Colombo et al. (2015); affects all reentry predictions.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A dynamical coherency gate for state recovery: Statistical requiem for the long arc of cislunar orbital mis-prediction." pith.science (2026). https://pith.science/paper/3T2QS7PO

@misc{pith2026250622748,
  author       = {Pith},
  title        = {Pith review of: A dynamical coherency gate for state recovery: Statistical requiem for the long arc of cislunar orbital mis-prediction},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/3T2QS7PO}},
  note         = {Machine review of arXiv:2506.22748}
}
read the original abstract

We present a statistically grounded methodology for recovering physically consistent initial conditions from historical two-line element sets (TLEs), enabling accurate long-arc trajectory reconstruction for distant and highly eccentric Earth satellites. The approach combines numerical averaging of osculating orbital elements with Gaussian-mixture-model (GMM) filtering with robust fallback strategies to isolate a dynamically coherent subset of mean elements at a common reference epoch. From this filtered ensemble, a representative osculating element is reconstructed, yielding a recovered Cartesian state vector with predictive capability far exceeding that of raw TLE-based propagations and existing approaches. We apply this method to the case of OGO-1 (1964-054A), a spacecraft launched into a cislunar orbit and tracked intermittently over five decades. Despite large observational gaps and significant secular evolution, our recovered state, used within an unscented Kalman filter (UKF), accurately reproduces the full trajectory, including its atmospheric reentry in August 2020, to within one day of the true decay time. This result demonstrates the viability of filtered TLE statistics as a proxy for precise orbit determination, particularly in dynamical regimes where resonances and long-period perturbations dominate. The techniques presented here provide a framework for trajectory reconstruction and prediction using only publicly available data and are broadly applicable to the study of high-altitude debris objects and legacy space missions whose original tracking and covariance data are unavailable.

Figures

Figures reproduced from arXiv: 2506.22748 by the authors.

Figure 1
Figure 1. A snapshot of the historic and current cataloged xGEO space objects in the planes of semimajor axis-ecliptic inclination (a, i) (top panel) and semimajor axis-eccentricity (a, e) (bottom panel), where the colored circles correspond to the osculating elements obtained from the latest TLE of each object. The geocentric semi-major axis of GEO is indicated by the vertical line at a„6.6RC. Objects that reach the Earth-gr… view at source ↗
Figure 2
Figure 2. Time history of the osculating Keplerian elements of NASA’s Orbiting Geophysical Obser￾vatory 1 (OGO 1, 1964-054A), derived from TLEs using the SGP4 algorithm at each native epoch. The orbital elements are measured with respect to the ecliptic frame. This xGEO satellite was launched on 5 September 1964 01:23 UTC (MJD 38643.05764) into a distant, highly eccentric orbit (a „ 12.7 R‘, e „ 0.918), exhibiting pronounced … view at source ↗
Figure 3
Figure 3. Schematic of the coherency-gated initialization process ( [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: Reentry epoch predictions for OGO-1, derived from forward propagation of individual TLEs under high-fidelity dynamics. The distribution reveals substantial variance in predicted decay dates, a manifestation of both observational sparsity and structural inconsistencies …
Figure 5
Figure 5. Figure 5: Application of Algorithm 1 to the approximately 3,200 TLEs of OGO-1 spanning from 23 June 2000 to 20 August 2020, with the coherency gate epoch set to 6 July 2002. A Gaussian mixture model (GMM) with fallback logic was first applied (top panel), followed by a median ab…
Figure 6
Figure 6. Figure 6: The multi-decadal orbital evolution of NASA’s Orbiting Geophysical Observatory 1 (OGO 1, 1964-054A), obtained using an ASSIST prediction from the recovered state at the coherency gate epoch of 6 July 2002 7:34 UTC (MJD 52461.31528), compared to the historical TLE time …
Figure 7
Figure 7. Figure 7: Reentry epoch predictions for OGO-1, showing the ensemble of raw TLE-based predictions (background, as in [PITH_FULL_IMAGE:figures/full_fig_p016_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

71 extracted references · 68 canonical work pages

  1. [1]

    J., 2020

    Amato, D., Malhotra, R., Sidorenko, V., Rosengren, A. J., 2020. Lunar close encounters compete with the circumterrestrial Lidov--Kozai effect: The dynamical demise of Luna 3 . Celestial Mechanics and Dynamical Astronomy 132, 35 (18 pp.)

  2. [2]

    F., 1963

    Arenstorf, R. F., 1963. Existence of periodic solutions passing near both masses of the restricted three-body problem. AIAA Journal 1, 238--240

  3. [3]

    E., 1976

    Ash, M. E., 1976. Doubly averaged effect of the Moon and Sun on a high altitude Earth satellite orbit. Celestial Mechanics 14, 209--238

  4. [4]

    A comprehensive review on cislunar expansion and space domain awareness

    Baker-McEvilly , B., Bhadauria, S., Canales, D., Frueh, C., 2024. A comprehensive review on cislunar expansion and space domain awareness. Progress in Aerospace Sciences 147, 101019 (16 pp.)

  5. [5]

    Reliable and repeatable transit through cislunar space using 2:1 resonant spatial orbits

    Binder, D., Arnas, D., 2024. Reliable and repeatable transit through cislunar space using 2:1 resonant spatial orbits. Journal of Guidance, Control, and Dynamics 47, 1973--1979

  6. [6]

    Fragmentation characterization in the circular restricted three body problem for cislunar space domain awareness

    Black, A., Frueh, C., 2025. Fragmentation characterization in the circular restricted three body problem for cislunar space domain awareness. Advances in Space Research 75, 1177--1204

  7. [7]

    Solar-wind observations with satellite ESRO HEOS-1 in December 1968

    Bonetti, A., Moreno, G., Cantarano, S., et al., 1969. Solar-wind observations with satellite ESRO HEOS-1 in December 1968. Il Nuovo Cimento 64, 307--323

  8. [8]

    R., Bettinger, R

    Boone, N. R., Bettinger, R. A., 2021. Debris collision risk analysis following simulated cislunar spacecraft explosions. Journal of Spacecraft and Rockets 60, 668--684

Show all 71 references
  1. [9]

    A., Gernet, E

    Boyarchuk, A. A., Gernet, E. D., Gershberg, R. E., et al., 1988. Orbital astrophysical observatory `` Astron '': Results from 5 years of operation. Kosmicheskie Issledovaniya 26, 917--933

  2. [10]

    Interplanetary monitoring platform

    Butler, P., 1980. Interplanetary monitoring platform. Tech. Rep. TM-80758, NASA

  3. [11]

    C., Hoots, F., 2018

    Chao, C. C., Hoots, F., 2018. Applied Orbit Perturbation and Maintenances, 2nd Edition. Aerospace Press, El Segundo, CA

  4. [12]

    M., van der Weg , W., et al., 2015

    Colombo, C., Alessi, E. M., van der Weg , W., et al., 2015. End-of-life disposal concepts for libration point orbit and highly elliptical orbit missions. Acta Astronautica 110, 298--312

  5. [13]

    E., Scott, D

    Cook, G. E., Scott, D. W., 1967. Lifetimes of satellites in large-eccentricity orbits. Planetary and Space Science 15, 1549--1556

  6. [14]

    A., 1959

    Ehricke, K. A., 1959. Cislunar orbits. In: Birkhoff, G., Langer, R. E. (Eds.), Orbit Theory. American Mathematical Society, Providence, RI, pp. 48--74

  7. [15]

    A., Ditrikh, A

    Eismont, N. A., Ditrikh, A. V., Janin, G., et al., 2003. Orbit design for launching INTEGRAL on the Proton/Block-DM launcher. Astronomy and Astrophysics 411, L37--L41

  8. [16]

    A., 2015

    Ely, T. A., 2015. Transforming mean and osculating elements using numerical methods. The Journal of the Astronautical Sciences 62, 21--43

  9. [17]

    Dynamics of satellites with multi-day periods

    \'Erdi, B., 1999. Dynamics of satellites with multi-day periods. In: Steves, B., Roy, A. (Eds.), The Dynamics of Small Bodies in the Solar System. Kluwer Academic Publishers, Dordrecht, pp. 303--307

  10. [18]

    Assessment and categorization of TLE orbit errors for the US SSN catalogue

    Flohrer, T., Krag, H., Klinkrad, H., 2008. Assessment and categorization of TLE orbit errors for the US SSN catalogue. In: Advanced Maui Optical and Space Surveillance Technologies Conference (AMOS). Maui, Hawaii

  11. [19]

    Accuracy of two-line-element data for geostationary and high-eccentricity orbits

    Fr \"u h, C., Schildknecht, T., 2012. Accuracy of two-line-element data for geostationary and high-eccentricity orbits. Journal of Guidance Control and Dynamics 35, 1483--1491

  12. [20]

    A., Gal'Perin, Y

    Galeev, A. A., Gal'Perin, Y. I., Zelenyi, L. M., 1996. The INTERBALL project to study solar-terrestrial physics. Kosmicheskie Issledovaniya 34, 339--362

  13. [21]

    K., Bucci, L., et al., 2023

    Guardabasso, P., Skoulidou, D. K., Bucci, L., et al., 2023. Analysis of accidental spacecraft break-up events in cislunar space. Advances in Space Research 72, 1550--1569

  14. [22]

    I., Amato, D., 2025

    Hallgarten La Casta , M. I., Amato, D., 2025. Debiasing of two-line element sets for batch least squares pseudo-orbit determination in MEO and GEO . Advances in Space Research 75, 7259--7289

  15. [23]

    Predicting cislunar orbit lifetimes from initial orbital elements, arXiv:2503.03892

    Higgins, D., Yeager, T., McGill, P., et al., 2025. Predicting cislunar orbit lifetimes from initial orbital elements, arXiv:2503.03892

  16. [24]

    J., Akmal, A., Farnocchia, D., Rein, H., Payne, M

    Holman, M. J., Akmal, A., Farnocchia, D., Rein, H., Payne, M. J., Weryk, R., Tamayo, D., Hernandez, D. M., 2023. ASSIST : An ephemeris-quality test-particle integrator. The Planetary Science Journal 4, 69 (9 pp.)

  17. [25]

    J., Chow, C

    Holzinger, M. J., Chow, C. C., Garretson, P., 2021. A primer on cislunar space. Tech. Rep. 2021-1271, AFRL

  18. [26]

    R., Glover, R

    Hoots, F. R., Glover, R. A., Schumacher Jr , P. W., 2004. History of analytical orbit modeling in the U. S. Space Surveillance System . Journal of Guidance Control and Dynamics 27, 174--185

  19. [27]

    S., 1962

    Huang, S. S., 1962. Preliminary study of orbits of interest for Moon probes. The Astronomical Journal 67, 304--310

  20. [28]

    P., Qureshi, R

    Hughes, S. P., Qureshi, R. H., Cooley, D. S., Parker, J. J. K., Grubb, T. G., 2014. Verification and validation of the General Mission Analysis Tool (GMAT) . In: AIAA/AAS Astrodynamics Specialist Conference. San Diego, California, Paper AIAA 2014-4151

  21. [29]

    A., 1976

    Janin, G., Roth, E. A., 1976. Decay of a highly eccentric satellite. Celestial Mechanics 14, 141--149

  22. [30]

    XMM-Newton observatory: I

    Jansen, F., Lumb, D., Altieri, B., et al., 2001. XMM-Newton observatory: I . The spacecraft and operations. Astronomy and Astrophysics 365, L1--L6

  23. [31]

    S., Kreisman, B

    Kardashev, N. S., Kreisman, B. B., Pogodin, A. V., et al., 2014. Orbit design for the Spektr-R spacecraft of the ground-space interferometer. Kosmicheskie Issledovaniya 52, 366--375

  24. [32]

    W., Strong, I

    Klebesadel, R. W., Strong, I. B., Olson, R. A., 1973. Observations of gamma-ray bursts of cosmic origin. The Astrophysical Journal 182, L85--L88

  25. [33]

    Perturbations of the orbits of artificial satellites by an attraction of external bodies

    Kopal, Z., 1967. Perturbations of the orbits of artificial satellites by an attraction of external bodies. Icarus 6, 298--314

  26. [34]

    A., 2025

    Koplow, D. A., 2025. Pave outer space and put up a parking lot: Lagrange points should be the common heritage of mankind. Michigan Journal of International Law 46, 403--461

  27. [35]

    A., 1959

    Kozai, Y., Whitney, C. A., 1959. Anticipated orbital perturbations of satellite 1959 Delta Two . Smithsonian Astrophysical Observatory Special Report 30, 15 (pp.)

  28. [36]

    Improved orbit prediction using two-line elements

    Levit, C., Marshall, W., 2011. Improved orbit prediction using two-line elements. Advances in Space Research 47, 1107--1115

  29. [37]

    L., 1963

    Lidov, M. L., 1963. On the approximated analysis of the orbit evolution of artificial satellites. In: Roy, M. (Ed.), Dynamics of Satellites. Springer--Verlag, Berlin, pp. 168--179

  30. [38]

    A., Gondelach, D

    Lidtke, A. A., Gondelach, D. J., Armellin, R., 2019. Optimising filtering of two-line element sets to increase re-entry prediction accuracy for GTO objects. Advances in Space Research 63, 1289--1317

  31. [39]

    TLE outlier detection based on expectation maximization algorithm

    Liu, J., Liu, L., Du, J., Sang, J., 2021. TLE outlier detection based on expectation maximization algorithm. Advances in Space Research 68, 2695--2712

  32. [40]

    Liu, J. J. F., Segrest, J., Szebehely, V., 1986. Orbit mechanics of deep space probes. The Journal of the Astronautical Sciences 34, 171--187

  33. [41]

    E., 1972

    Lowrey, B. E., 1972. Ephemeris of a highly eccentric orbit: Explorer 28. Celestial Mechanics 5, 107--125

  34. [42]

    H., 1963

    Ludwig, G. H., 1963. The orbiting geophysical observatories. Space Science Reviews 2, 175--218

  35. [43]

    Correcting TLEs at epoch: Application to the GPS constellation

    Ly, D., Lucken, R., Giolito, D., 2020. Correcting TLEs at epoch: Application to the GPS constellation. Journal of Space Safety Engineering 61, 302--306

  36. [44]

    H., 1964

    Milstead, A. H., 1964. Launch windows for orbital missions. Tech. Rep. SSD-TDR-64-31, USSF

  37. [45]

    E., 1963

    Montgomery, H. E., 1963. EGO launch window study. Tech. Rep. TM X-55413, NASA

  38. [46]

    K., 1968

    Musen, P., Squires, R. K., 1968. Orbital mechanics. In: Hess, W., Mead, G. (Eds.), Introduction to Space Science, 2nd Edition. Gordon and Breach, New York, pp. 529--554

  39. [47]

    Computational exploration of the cislunar region and implications for debris mitigation

    Namazyfard, H., 2019. Computational exploration of the cislunar region and implications for debris mitigation. Master's thesis, University of Arizona

  40. [48]

    R., Prokhorenko, V

    Nazirov, R. R., Prokhorenko, V. I., Sheikhet, A. I., 2002. A retrospective geometric analysis of the long-term evolution of orbits and the ballistic lifetimes for satellites of the Prognoz series. Kosmicheskie Issledovaniya 40, 538--554

  41. [49]

    R., 1959

    Newton, R. R., 1959. Periodic orbits of a planetoid passing close to two gravitating masses. Smithsonian Contribution to Astrophysics 3, 69--78

  42. [50]

    J., Shute, B

    Paddack, S. J., Shute, B. E., 1965. IMP-C orbit and launch time analysis. Tech. Rep. TM X-55128, NASA

  43. [51]

    3:1/3:2 resonant orbits touring L_3 -- L_5 in cislunar space

    Peng, C., Zhang, Y., He, S., 2024. 3:1/3:2 resonant orbits touring L_3 -- L_5 in cislunar space. Advances in Space Research 73, 2499--2514

  44. [52]

    J., Ross, S

    Rawat, A., Kumar, B., Rosengren, A. J., Ross, S. D., 2025. Cislunar mean-motion resonances: Definitions , widths, and comparisons with resonant satellites, arXiv:2505.10138

  45. [53]

    M., Tamayo, D., 2019

    Rein, H., Hernandez, D. M., Tamayo, D., 2019. Hybrid symplectic integrators for planetary dynamics. Monthly Notices of the Royal Astronomical Society 485, 5490--5497

  46. [54]

    L., 1970

    Renard, M. L., 1970. Practical stability of high-eccentricity orbits quasi-normal to the ecliptic. Journal of Spacecraft and Rockets 7, 1208--1214

  47. [55]

    R., Winn, J

    Ricker, G. R., Winn, J. N., Vanderspek, R., et al., 2015. Transiting exoplanet survey satellite. Journal of Astronomical Telescopes, Instruments, and Systems 1, 014003 (10 pp.)

  48. [56]

    J., Scheeres, D

    Rosengren, A. J., Scheeres, D. J., 2021. Cis-lunar trajectories. In: Cudnik, B. (Ed.), Encyclopedia of Lunar Science. Springer, Cham, p. (9 pp.)

  49. [57]

    A., 1970

    Roth, E. A., 1970. Launch window study for the highly eccentric orbit satellite HEOS-1 . Celestial Mechanics 2, 369--381

  50. [58]

    J., Shute, B

    Sandifer, R. J., Shute, B. E., 1962. Effect of solar-lunar perturbations on the lifetime of Explorer XII (abstract). The Astronomical Journal 67, 282

  51. [59]

    C., Smith, C

    Sang, J., Bennett, J. C., Smith, C. H., 2013. Estimation of ballistic coefficients of low altitude debris objects from historical two line elements. Advances in Space Research 52, 117--124

  52. [60]

    Long-period effects in nearly commensurable cases of the restricted three-body problem

    Schubart, J., 1964. Long-period effects in nearly commensurable cases of the restricted three-body problem. Smithsonian Astrophysical Observatory Special Report 149, (36 pp.)

  53. [61]

    I., 2017

    Shevchenko, I. I., 2017. The Lidov-Kozai Effect -- Applications in Exoplanet Research and Dynamical Astronomy. Springer, Berlin

  54. [62]

    P., Humphries, M., Maclean, J., et al., 2024

    Shorten, D. P., Humphries, M., Maclean, J., et al., 2024. Optimal proposal particle filters for detecting anomalies and manoeuvres from two line element data. Acta Astronautica 228, 709--723

  55. [63]

    E., Chiville, J., 1966

    Shute, B. E., Chiville, J., 1966. The lunar-solar effect on the orbital lifetimes of artificial satellites with highly eccentric orbits. Planetary and Space Science 14, 361--369

  56. [64]

    The Lidov-Kozai oscillation and Hugo von Zeipel

    Takashi, I., Ohtsuka, K., 2019. The Lidov-Kozai oscillation and Hugo von Zeipel . Monographs on Environment, Earth and Planets 7, 1--113

  57. [65]

    Numerical averaging in orbit prediction

    Uphoff, C., 1973. Numerical averaging in orbit prediction. AIAA Journal 11, 1512--1516

  58. [66]

    Lunar and solar perturbations on satellite orbits

    Upton, E., Bailie, A., Musen, P., 1959. Lunar and solar perturbations on satellite orbits. Science 130, 1710--1711

  59. [67]

    C., 2014

    Vaquero, M., Howell, K. C., 2014. Design of transfer trajectories between resonant orbits in the Earth--Moon restricted problem. Acta Astronautica 94, 302--317

  60. [68]

    A., Shockley, L

    Wilmer, A., Bettinger, R. A., Shockley, L. M., Holzinger, M. J., 2025. Preliminary investigation and proposal of periodic orbits and their utilization for logistics in the cislunar regime. Space Policy, 101635 (17 pp.)

  61. [69]

    The quest to conquer Earth's space junk problem

    Witze, A., 2018. The quest to conquer Earth's space junk problem. Nature 24, 25--26

  62. [70]

    , " * write output.state after.block = add.period write newline

    ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type url volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.senten...

  63. [71]

    write newline

    " write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 gl...

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.