Satellite orbits around a central body with arbitrary zonal harmonics are considered in a relativistic framework. Our starting point is the relativistic Celestial Mechanics based upon the first post-Newtonian approximation to Einstein’s theory of gravity as it has been formulated by Damour et al. (Phys Rev D 43:3273–3307, 1991; 45:1017–1044, 1992; 47:3124–3135, 1993; 49:618–635, 1994). Since effects of order \((\mathrm{GM}/c^2R) \times J_k\) with \(k \ge 2\) for the Earth are very small (of order \( 7 \times 10^{-10}\,\times \,J_k\)) we consider an axially symmetric body with arbitrary zonal harmonics and a static external gravitational field. In such a field the explicit \(J_k/c^2\)-terms (direct terms) in the equations of motion for the coordinate acceleration of a satellite are treated first with first-order perturbation theory. The derived perturbation theoretical results of first order have been checked by purely numerical integrations of the equations of motion. Additional terms of the same order result from the interaction of the Newtonian \(J_k\)-terms with the post-Newtonian Schwarzschild terms (relativistic terms related to the mass of the central body). These ‘mixed terms’ are treated by means of second-order perturbation theory based on the Lie-series method (Hori–Deprit method). Here we concentrate on the secular drifts of the ascending node \( \) and argument of the pericenter \( \). Finally orders of magnitude are given and discussed.