I would like to know the position and orientation of the Earth, relative to the
ID: 1374330 • Letter: I
Question
I would like to know the position and orientation of the Earth, relative to the sun, as a function of time.
I.e. select any coordinate system with the Sun at the origin so that the earth lies in the x-y coordinate plane. Provide a function which takes as input the time and produces
The x-y position of the center of the Earth
The unit vector from the center of the Earth to the North Pole.
The orientation of the Earth around this vector
This function is chaotic in general. I would be happy with a decent approximation.
This problem is solved as part of the "what is the position of the sun in the sky" problem. Solutions to this problem tend not to be well separable (at least I have trouble separating my question out of the complete position-of-sun-in-sky solution.)
Explanation / Answer
For a precise measurement of the movement of the Earth around the Sun as a function of time you will need to correct for the changes in the orbit of the Earth caused by the gravitational pull of the different objects on the usual two-body problem.
These perturbations are made to the usual Kepler equation (of course, if these perturbations are long range, i.e., months or years). If all of this is unfamiliar to you, please read Chapter 2 of Murray and Dermott's, 1999 book.
On the following days (if I can fix my old hard drive that seemed to die), I can post a numerical solution to Kepler's equation in Python that I did on an undergraduate course on Extrasolar Planets. However, I must tell that it's not that hard and you can understand easily how to do it if you read Chapter 2 of the book I quoted. If you get stuck, please come back to the post and we'll try to help :-).