Computing osculating orbital elements

In summary, Kepler calculates Earth's orbital elements in heliocentric coordinates, while Horizons calculates them in barycentric coordinates. This causes a difference in the results.
  • #1
cptolemy
48
1
Good afternoon,

I wonder if someone can help me in a small doubt.

I'm trying to calculate the orbital parameters of some solar system bodies. It's quite easy knowing their positions and velocities.

But my question is this: these two components, position and velocity, must they be heliocentric or barycentric?

And as far the equinox and ecliptic of reference, must it be the J2000 or the mean of date?

I'm asking this, because I get different ascending nodes, for instance for Earth, compared to the NASA release values. Actually, their value is for J2000 about -11 degrees, which is strange; it should be a positive number (adding 360).

Does anyone has expertise in this field?

Cheers,

Kepler
 
Astronomy news on Phys.org
  • #2
Position and velocity can be either heliocentric or barycentric. If you use heliocentric position and velocity, you'll get heliocentric orbital elements.

For computing Earth's orbital elements, it doesn't make much sense to use barycentric since Earth orbits inside the orbit of Jupiter, the Sun's biggest perturber.

-11 and 349 are just 2 ways of saying the same thing. Where is the NASA page that states -11 degrees? Using Horizons, I get 160 degrees.

Regardless, it's not surprising to get vastly different numbers when computing the ascending node of an orbit with no inclination. Since Earth's orbit defines the ecliptic plane, Earth has no inclination. Only the z-component tugs from the Moon and other planets barely give Earth an instantaneous non-zero inclination, but it changes all the time, causing Earth's calculated Longitude of Ascending Node to dance all over the place.

You'll get similar poorly-defined results if you try to compute Earth's argument of perihelion, since it's measured from the poorly-defined LAN.

To a lesser degree, Venus' argument of Perihelion dances all over the place too, since its orbit is very circular. A perfectly-circular orbit has no perihelion.

How do your computed values compare to NASA's values for other orbital elements?
 
  • #3
cptolemy said:
Good afternoon,
But my question is this: these two components, position and velocity, must they be heliocentric or barycentric?
And as far the equinox and ecliptic of reference, must it be the J2000 or the mean of date?
Depends on what you want out of it. If you start with heliocentric, J2000 pos/vel, you will get heliocentric, J2000 orbital elements. Just be careful in doing your comparisons with someone else's results (even NASA's) - if they started with, for instance, barycentric, J2000 pos/vel, their results will be different. Or, if they started with some other data set - pos/vel as measured by a spacecraft rather than pos/vel computed and published in an almanac - again, different results if you happen to be comparing yours to the almanac elements.

The best thing to do, if you're verifying a set of equations or something, is to find some output data with known input data and put that through your program.
 
  • Like
Likes cptolemy
  • #5
Hi everyone,

Thanks for trying to help. Oh, and of course, thank you very much for mention the "great Jean Meeus" - one of my idols regarding time fit astronomical alghorythms. I have - I think - all his books. Willman Bell is a great provider of specific astronomy technical books that we do not find elsewhere...

I did some digging - and reversed the problem. Calculated the state vectors of Horizons from their orbital elements. Conclusion: they were diferent from mine... thus the difference.

Horizons calculates the simple vectors choice in barycentric coordinates. As for the orbital elements, it's calculation is made in heliocentric ones, referred to the equinox and equator plane of J2000, TDB scale (or similar - they use Coordinate Time), with true coordinates - no light effect, aberration or whatever. Not astrometric - just true.

I'm getting now very similar results since my state vectors differ only in some decimals (10^-9 aprox.). But the results are quite accurate, I must say :)

So, up until next time.

Cheers,

CPtolemy
 

1. What are osculating orbital elements?

Osculating orbital elements are a set of six parameters used to describe the orbit of a celestial object, such as a planet or satellite, around another larger object, such as a star or planet. These elements include the semi-major axis, eccentricity, inclination, longitude of the ascending node, argument of periapsis, and true anomaly.

2. How are osculating orbital elements calculated?

Osculating orbital elements are calculated using mathematical models and calculations based on observations of the celestial object's position and velocity. These elements can be updated and refined over time as more precise observations are made.

3. Why are osculating orbital elements important?

Osculating orbital elements are important because they provide a standardized way to describe the orbit of a celestial object. They allow for predictions of future positions and movements of the object, as well as comparisons and analyses of different orbits.

4. What is the difference between osculating orbital elements and mean orbital elements?

The main difference between osculating and mean orbital elements is that osculating elements are based on the current position and velocity of the celestial object, while mean elements are averaged over a longer period of time. Osculating elements are more accurate for predicting short-term orbital behavior, while mean elements are better for long-term predictions.

5. Are osculating orbital elements constant?

No, osculating orbital elements are not constant. They can change over time due to various factors, such as gravitational influences from other objects, atmospheric drag, and the effects of relativity. However, they are usually relatively stable and can be updated as needed based on new observations.

Similar threads

  • Astronomy and Astrophysics
Replies
3
Views
3K
  • Astronomy and Astrophysics
Replies
4
Views
4K
  • Astronomy and Astrophysics
Replies
2
Views
3K
  • Astronomy and Astrophysics
Replies
7
Views
4K
  • Astronomy and Astrophysics
Replies
4
Views
2K
Replies
1
Views
2K
Replies
13
Views
2K
  • Astronomy and Astrophysics
Replies
2
Views
3K
  • Astronomy and Astrophysics
Replies
1
Views
5K
  • Astronomy and Astrophysics
2
Replies
45
Views
81K
Back
Top