Correspondence should be addressed to Rodolpho Vilhena de Moraes,

under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Low-altitude, near-polar orbits are very desirable as science orbits for missions to planetary satellites, such as the Earth’s Moon. In this paper, we present an analytical theory with numerical simulations to study the orbital motion of lunar low-altitude artificial satellite. We consider the problem of an artificial satellite perturbed by the nonuniform distribution of the mass of the Moon J2–J5, J7, and C22. The conditions to get frozen orbits are presented. Using an approach that considers the single-averaged problem, we found families of periodic orbits for the problem of an orbiter travelling around the Moon, where frozen orbits valid for long periods of time are found. A comparison between the models for the zonal and tesseral harmonics coefficients is presented. 1