Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. Join them; it only takes a minute:

Sign up
Here's how it works:
  1. Anybody can ask a question
  2. Anybody can answer
  3. The best answers are voted up and rise to the top

I am designing an orbit around Mercury. I know the values I want for the semi-major axis, eccentricity, inclination, and RAAN. I want the altitude of closest approach (periapse) to occur at $60.0^{\circ}$ N latitude above Mercury's surface.

My Question: How do I convert this constraint for the latitude of closest approach into an argument of periapsis, $\omega$?

Specifically, this is modeled after the orbit below:

enter image description here

share|cite|improve this question
up vote 1 down vote accepted

The answer is:
$ \sin(\omega) = \frac{\sin(lat)}{\sin(i)} $

this expression has has two solutions:
$ \omega = \arcsin \left( \frac{\sin(lat)}{\sin(i)} \right), \ \ \text{and} \ \ \pi - \arcsin \left( \frac{\sin(lat)}{\sin(i)} \right) $

The second solution matches the diagram.

So if the latitude is 60° and the inclination is 82.5°, then the argument of periapsis (ω) is 119.13225456935°.

(The term ''perigee'' applies only to orbits around the Earth; the correct term for Mercury is ''perihermion''. The general term is ''periapsis''.)

The perihelion direction vector is: $ \vec{q} = [1 \ 0 \ 0] \cdot \operatorname{rot}_z(\omega) \cdot \operatorname{rot}_x(i) $.
where rotz is the rotation matrix about the z axis, and rotx is the rotation matrix about the x axis.

Convert the perihelion direction vector to spherical coordinates, set the latitude equal to 60°, then solve for ω.

share|cite|improve this answer
Great answer. Thank you! – dinkelk Dec 18 '12 at 22:57

Your Answer


By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.