  This work focuses on providing closed form analytical expressions to define
frozen orbits under the effects of the zonal harmonics of an Earth-like planet.
Particularly, the perturbation effects from the terms J2, J3, J4, J5, J6, and
J7 are considered in this work. This is done using a power series expansion in
the small parameter that allow not only to provide an approximate solution to
the system, but also to determine the analytical expressions that define the
initial osculating conditions that generate frozen orbits. As a result of that,
the proposed methodology allows to study the bifurcation of frozen orbits close
to the critical inclination by purely analytical methods. Additionally, the
derivation of the secular variation of the orbital elements as well as the
transformation from osculating to mean elements is provided based on the second
order analytical solution proposed in this work. Examples of application are
also provided to show the error performance of the results included in this
document.
