Vallado (2010): ICATT Astrodynamics Tutorial

Source

Vallado, David A. “Fundamentals of Astrodynamics and Applications.” Tutorial Lectures at the 4th International Conference on Astrodynamics Tools and Techniques (ICATT), Madrid, Spain, May 2010. 120 slides. From ESA trajectory server.

Download PDF

Teaching Presentation, Not Research Paper

This is a 120-slide tutorial condensing Vallado’s textbook into a single session. It provides an accessible overview of the full SGP4 pipeline — from time systems through observation geometry — without the mathematical depth of the individual foundation papers. Use this as a map of the territory before diving into the specifics.

Why This Matters for SGP4

The tutorial is structured around a single end-to-end problem: “Determine when you can see a satellite from a ground site.” This mirrors Astrolock’s exact pipeline:

1. Time systems — converting between UTC, UT1, TAI, TT, and GPS time
2. Coordinate frames — TEME, J2000, ITRF, geodetic, topocentric
3. Propagation — Kepler’s equation, perturbation theory, SGP4
4. Orbit determination — observations to orbital elements
5. Observation geometry — altitude, azimuth, range, visibility constraints

Where the individual papers in this archive each address one piece of that pipeline, the tutorial shows how the pieces connect. It is a pedagogical bridge between the foundation papers and the textbook.

Tutorial Structure

The 120 slides map directly to chapters in the textbook:

Slides	Topic	Textbook Chapter
1—7	Introduction	—
8—20	Time and coordinate systems	Ch. 3
21—49	Kepler’s equation and two-body	Ch. 1—2
50—80	Perturbations and propagation	Ch. 8—9
81—108	Orbit determination and estimation	Ch. 10
109—120	Applications	Ch. 11

The weighting is telling: perturbations and propagation (31 slides) and orbit determination (28 slides) together consume half the tutorial. These are the two domains where SGP4 operates.

Perturbation Coverage

Slides 50—80 cover the same forces that SGP4 models:

• Gravity harmonics ($J_2$ through $J_5$): zonal, sectoral, and tesseral harmonics with clear visualizations (slides 53—56)
• Nodal regression and apsidal rotation: visual demonstrations of how oblateness rotates the orbital plane and the line of apsides (slides 57—59)
• Atmospheric drag: force model and altitude-dependent effects
• Third-body perturbations: Sun, Moon, and planetary gravitational influences
• Solar radiation pressure: area-to-mass ratio effects on orbit evolution

The Perturbation Decomposition

Slide 52 presents the decomposition of perturbation effects into three categories:

$\delta\text{(element)} = \delta_\text{secular} + \delta_\text{long-periodic} + \delta_\text{short-periodic}$

This is the same three-stage structure that Brouwer (1959) derived analytically and that the SGP4 FORTRAN implements as three sequential code blocks:

Perturbation Type	Period	SGP4 FORTRAN Lines
Secular	Cumulative	29—52
Long-periodic	Months	53—72
Short-periodic	One orbit	120—175

The tutorial’s visual treatment of this decomposition is the clearest available diagram of how these terms combine — more accessible than the formal derivation in the original papers.

The Applicability Chart

Slide 51 shows which perturbation forces matter at which orbital altitudes. This is the physical justification for two key design decisions in the Spacetrack propagators:

• SGP4’s drag model dominates at LEO (<2000 km) where atmospheric density is significant, but is negligible for higher orbits
• SDP4’s third-body perturbations (Sun and Moon) become the dominant non-$J_2$ force for deep-space orbits (period $\geq$ 225 minutes)

The 225-minute threshold that the STR#3 DRIVER uses to select between SGP4 and SDP4 sits at the crossover point on this chart.

Orbit Determination Overview

Slides 81—108 cover the statistical orbit determination pipeline, progressing from raw observations to refined orbital elements:

Slides	Topic	Key Concept
81—99	Observation types and errors	Bias, noise, drift in radar and optical data
100	Observation modeling	Transforming measurements to state-space quantities
101	Statistical foundations	Multidimensional probability distributions
102	Covariance matrices	$\mathbf{P} = (\mathbf{A}^T \mathbf{W} \mathbf{A})^{-1}$
104	Least squares	Step-by-step matrix inverse approach
105	Sequential batch least squares	Bayesian update for incorporating new observations
106—107	Extended Kalman filter	Full prediction-update cycle
108	Fit spans and averaging	Connection to mean elements (what TLEs contain)

This OD material complements the AIAA 2008-6770 paper by providing the theoretical foundation in an accessible format. Where the 2008 paper gives the full mathematical treatment with operational results, these slides provide the conceptual framework needed to follow that paper’s derivations.

Slide 108’s treatment of fit spans and mean element averaging is particularly relevant for understanding why TLE accuracy degrades with age — the mean elements are a best-fit average over the observation span, not an instantaneous state.

Edition Mismatch

The slides reference the 3rd edition textbook (2007), but the tutorial was delivered in 2010. Our archive has the 4th edition (2013). The SGP4-relevant content — perturbation theory, propagation models, and orbit determination — is substantially the same between editions. The 4th edition adds updated numerical examples and expanded treatment of newer propagation methods, but the foundational material that maps to SGP4 is unchanged.

---

Vallado (2013): Astrodynamics Textbook The comprehensive textbook this tutorial condenses

Vallado (2008): SGP4 Orbit Determination Full treatment of the OD pipeline covered in slides 100-108

Rev-1 Overview The SGP4 modernization referenced throughout the tutorial

Brouwer (1959): Mean Element Theory The perturbation decomposition shown in slide 52