Orbital elements are the parameters required to uniquely identify a specific orbit. In celestial mechanics these elements are generally considered in classical two-body systems, where a Kepler orbit is used (derived from Newton's laws of motion and Newton's law of universal gravitation). There are many different ways to mathematically describe the same orbit, but certain schemes each consisting of a set of six parameters are commonly used in astronomy and orbital mechanics.

A real orbit (and its elements) changes over time due to gravitational perturbations by other objects and the effects of relativity. A Keplerian orbit is merely an idealized, mathematical approximation at a particular time.

Keplerian elements[link]

In this diagram, the orbital plane (yellow) intersects a reference plane (gray). For earth-orbiting satellites, the reference plane is usually the Earth's equatorial plane, and for satellites in solar orbits it is the ecliptic plane. The intersection is called the line of nodes, as it connects the center of mass with the ascending and descending nodes. This plane, together with the Vernal Point, (♈) establishes a reference frame.

The traditional orbital elements are the six Keplerian elements, after Johannes Kepler and his laws of planetary motion.
When viewed from an inertial frame, two orbiting bodies trace out distinct trajectories. Each of these trajectories has its focus at the common center of mass. When viewed from the non-inertial frame of one body only the trajectory of the opposite body is apparent; Keplerian elements describe these non-inertial trajectories. An orbit has two sets of Keplerian elements depending on which body is used as the point of reference. The reference body is called the primary, the other body is called the secondary. The primary is not necessarily more massive than the secondary, and even when the bodies are of equal mass, the orbital elements depend on the choice of the primary.

The main two elements that define the shape and size of the ellipse:

• Eccentricity (Failed to parse (Missing texvc executable; please see math/README to configure.): e\,\!

) - shape of the ellipse, describing how flattened it is compared with a circle. (not marked in diagram)

• Semimajor axis (Failed to parse (Missing texvc executable; please see math/README to configure.): a\,\!

) - the sum of the periapsis and apoapsis distances divided by two. For circular orbits the semimajor axis is the distance between the bodies, not the distance of the bodies to the center of mass.

Two elements define the orientation of the orbital plane in which the ellipse is embedded:

• Inclination - vertical tilt of the ellipse with respect to the reference plane, measured at the ascending node (where the orbit passes upward through the reference plane) (green angle i in diagram).
• Longitude of the ascending node - horizontally orients the ascending node of the ellipse (where the orbit passes upward through the reference plane) with respect to the reference frame's vernal point (green angle Ω in diagram).

And finally:

• Argument of periapsis defines the orientation of the ellipse (in which direction it is flattened compared to a circle) in the orbital plane, as an angle measured from the ascending node to the semimajor axis. (violet angle Failed to parse (Missing texvc executable; please see math/README to configure.): \omega\,\!

in diagram)

• Mean anomaly at epoch (Failed to parse (Missing texvc executable; please see math/README to configure.): M_o\,\!

) defines the position of the orbiting body along the ellipse at a specific time (the "epoch").

The mean anomaly is a mathematically convenient "angle" which varies linearly with time, but which does not correspond to a real geometric angle. It can be converted into the true anomaly Failed to parse (Missing texvc executable; please see math/README to configure.): \nu\,\! , which does represent the real geometric angle in the plane of the ellipse, between periapsis (closest approach to the central body) and the position of the orbiting object at any given time. Thus, the true anomaly is shown as the red angle Failed to parse (Missing texvc executable; please see math/README to configure.): \nu\,\!

in the diagram, and the mean anomaly is not shown.

The angles of inclination, longitude of the ascending node, and argument of periapsis can also be described as the Euler angles defining the orientation of the orbit relative to the reference coordinate system.

Note that non-elliptic orbits also exist; If the eccentricity is greater than one, the orbit is a hyperbola. If the eccentricity is equal to one and the angular momentum is zero, the orbit is radial. If the eccentricity is one and there is angular momentum, the orbit is a parabola.

Required parameters[link]

Given an inertial frame of reference and an arbitrary epoch (a specified point in time), exactly six parameters are necessary to unambiguously define an arbitrary and unperturbed orbit.

This is because the problem contains six degrees of freedom. These correspond to the three spatial dimensions which define position (the x, y, z in a Cartesian coordinate system), plus the velocity in each of these dimensions. These can be described as orbital state vectors, but this is often an inconvenient way to represent an orbit, which is why Keplerian elements are commonly used instead.

Sometimes the epoch is considered a "seventh" orbital parameter, rather than part of the reference frame.

If the epoch is defined to be at the moment when one of the elements is zero, the number of unspecified elements is reduced to five. (The sixth parameter is still necessary to define the orbit; it is merely numerically set to zero by convention or "moved" into the definition of the epoch with respect to real-world clock time.)

Alternative parametrizations[link]

Keplerian elements can be obtained from orbital state vectors (x-y-z coordinates for position and velocity) by manual transformations or with computer software.^[1]

Other orbital parameters can be computed from the Keplerian elements such as the period, apoapsis, and periapsis. (When orbiting the earth, the last two terms are known as the apogee and perigee.) It is common to specify the period instead of the semi-major axis in Keplerian element sets, as each can be computed from the other provided the standard gravitational parameter, GM, is given for the central body.

Instead of the mean anomaly at epoch, the mean anomaly Failed to parse (Missing texvc executable; please see math/README to configure.): M\,\! , mean longitude, true anomaly Failed to parse (Missing texvc executable; please see math/README to configure.): \nu_o\,\! , or (rarely) the eccentric anomaly might be used.

Using, for example, the "mean anomaly" instead of "mean anomaly at epoch" means that time Failed to parse (Missing texvc executable; please see math/README to configure.): t

must be specified as a "seventh" orbital element. Sometimes it is assumed that mean anomaly is zero at the epoch (by choosing the appropriate definition of the epoch), leaving only the five other orbital elements to be specified.

Different sets of elements are used for various astronomical bodies. The eccentricity, e, and either the semi-major axis, a, or the distance of periapsis, q, are used to specify the shape and size of an orbit. The angle of the ascending node, Ω, the inclination, i, and the argument of periapsis, ω, or the longitude of periapsis, ϖ, specify the orientation of the orbit in its plane. Either the longitude at epoch, L₀, the mean anomaly at epoch, M₀, or the time of perihelion passage, T₀, are used to specify a known point in the orbit. The choices made depend whether the vernal equinox or the node are used as the primary reference. The semi-major axis is known if the mean motion and the gravitational mass are known.^[2][3]

Sets of Elements for Astronomical Objects

Object	Elements used
Major Planet	e,a, i, Ω, ϖ, L0
Comet	e, q,i,Ω, ω,T0
Asteroid	e,a,i,Ω, ω, M0
TLE	e,i,Ω,ω, n,M0

Euler angle transformations[link]

The angles Failed to parse (Missing texvc executable; please see math/README to configure.): \Omega , i, \omega

are the Euler angles (Failed to parse (Missing texvc executable; please see math/README to configure.): \alpha , \beta, \gamma
with the notations of that article) characterizing the orientation of the coordinate system

Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{x},\hat{y},\hat{z}

from the inertial coordinate frame Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{I},\hat{J},\hat{K}

where:

Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{I},\hat{J}

is in the equatorial plane of the central body. Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{I}
is in the direction of the vernal equinox. Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{J}
is perpendicular to Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{I}
and with Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{I}
defines the reference plane. Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{K}
is perpendicular to the reference plane.

Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{x},\hat{y}

are in the orbital plane and with Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{x}
in the direction to the pericenter (periapsis). Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{z}
is perpendicular to the plane of the orbit. Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{y}
is mutually perpendicular to Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{x}
and Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{z}

.

Then, the transformation from the Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{I},\hat{J},\hat{K}

coordinate frame to the Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{x},\hat{y},\hat{z}
frame with the Euler angles Failed to parse (Missing texvc executable; please see math/README to configure.): \Omega , i, \omega
is:

Failed to parse (Missing texvc executable; please see math/README to configure.): x_1= \cos \Omega \cdot \cos \omega - \sin \Omega \cdot \cos i \cdot \sin \omega

Failed to parse (Missing texvc executable; please see math/README to configure.): x_2= \sin \Omega \cdot \cos \omega + \cos \Omega \cdot \cos i \cdot \sin \omega

Failed to parse (Missing texvc executable; please see math/README to configure.): x_3= \sin i \cdot \sin \omega

Failed to parse (Missing texvc executable; please see math/README to configure.): y_1=-\cos \Omega \cdot \sin \omega - \sin \Omega \cdot \cos i \cdot \cos \omega

Failed to parse (Missing texvc executable; please see math/README to configure.): y_2=-\sin \Omega \cdot \sin \omega + \cos \Omega \cdot \cos i \cdot \cos \omega

Failed to parse (Missing texvc executable; please see math/README to configure.): y_3= \sin i \cdot \cos \omega

Failed to parse (Missing texvc executable; please see math/README to configure.): z_1= \sin i \cdot \sin \Omega

Failed to parse (Missing texvc executable; please see math/README to configure.): z_2=-\sin i \cdot \cos \Omega

Failed to parse (Missing texvc executable; please see math/README to configure.): z_3= \cos i\,

where

Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{x}= x_1\hat{I} + x_2\hat{J} + x_3\hat{K}

Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{y}= y_1\hat{I} + y_2\hat{J} + y_3\hat{K}

Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{z}= z_1\hat{I} + z_2\hat{J} + z_3\hat{K}

The transformation from Failed to parse (Missing texvc executable; please see math/README to configure.): \hat{x},\hat{y},\hat{z}

to Euler angles Failed to parse (Missing texvc executable; please see math/README to configure.): \Omega , i, \omega
is:

Failed to parse (Missing texvc executable; please see math/README to configure.): \Omega= \operatorname{arg}(\ -z_2\ ,\ z_1\ )

Failed to parse (Missing texvc executable; please see math/README to configure.): i = \operatorname{arg}(\ z_3\ ,\ \sqrt{{z_1}^2 + {z_2}^2}\ )

Failed to parse (Missing texvc executable; please see math/README to configure.): \omega= \operatorname{arg}(\ y_3\ ,\ x_3\ )

where Failed to parse (Missing texvc executable; please see math/README to configure.): \operatorname{arg}(x\ ,\ y)

signifies the polar argument that can be computed with

the standard function ATAN2(y,x) (or in double precision DATAN2(y,x)) available in for example the programming language FORTRAN.

Orbit prediction[link]

Under ideal conditions of a perfectly spherical central body, and zero perturbations, all orbital elements, with the exception of the Mean anomaly are constants, and Mean anomaly changes linearly with time, scaled by the Mean motion, Failed to parse (Missing texvc executable; please see math/README to configure.): n=\sqrt{\frac{\mu } {a^3}} .^[2] Hence if at any instant Failed to parse (Missing texvc executable; please see math/README to configure.): t_0

the orbital parameters are Failed to parse (Missing texvc executable; please see math/README to configure.): [e_0,a_0,i_0,\Omega_0,\omega_0,M_0]

, then the elements at time Failed to parse (Missing texvc executable; please see math/README to configure.): t_0+\delta t

is given by Failed to parse (Missing texvc executable; please see math/README to configure.): [e_0,a_0,i_0,\Omega_0,\omega_0,M_0+n\delta t]

Perturbations and elemental variance[link]

Unperturbed, two-body, Newtonian orbits are always conic sections, so the Keplerian elements define an ellipse, parabola, or hyperbola. Real orbits have perturbations, so a given set of Keplerian elements accurately describes an orbit only at the epoch. Evolution of the orbital elements takes place due to the gravitational pull of bodies other than the primary, the nonsphericity of the primary, atmospheric drag, relativistic effects, radiation pressure, electromagnetic forces, and so on.

Keplerian elements can often be used to produce useful predictions at times near the epoch. Alternatively, real trajectories can be modeled as a sequence of Keplerian orbits that osculate ("kiss" or touch) the real trajectory. They can also be described by the so-called planetary equations, differential equations which come in different forms developed by Lagrange, Gauss, Delaunay, Poincaré, or Hill.

Two-line elements[link]

Keplerian elements parameters can be encoded as text in a number of formats. The most common of them is the NASA/NORAD "two-line elements"(TLE) format [1], originally designed for use with 80-column punched cards, but still in use because it is the most common format, and can be handled easily by all modern data storages as well.

Depending on the application and object orbit, the data derived from TLEs older than 30 days can become unreliable. Orbital positions can be calculated from TLEs through the SGP/SGP4/SDP4/SGP8/SDP8 algorithms.^[4]

Example of a two line element:^[5]

 1 27651U 03004A   07083.49636287  .00000119  00000-0  30706-4 0  2692
 2 27651 039.9951 132.2059 0025931 073.4582 286.9047 14.81909376225249

See also[link]

• Beta Angle
• Ephemeris
• Orbital state vectors
• Proper orbital elements
• Osculating orbit

References[link]

External links[link]

• Keplerian Elements tutorial
• another tutorial
• Spacetrack Report No. 3, a really serious treatment of orbital elements from NORAD (in pdf format)
• Celestrak Two-Line Elements FAQ
• The JPL HORIZONS online ephemeris. Also furnishes orbital elements for a large number of solar system objects.
• NASA Planetary Satellite Mean Orbital Parameters.
• Introduction to exporting JPL planetary and lunar ephemerides
• State vectors: VEC2TLE Access to VEC2TLE software

Wikibooks has a book on the topic of Astrodynamics/Classical Orbit Elements