The second-order analytical state propagation under an arbitrary continuous thrust is proposed for relative motion in <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><tex-math notation="LaTeX">$J_{2}$</tex-math></inline-formula> -perturbed elliptic chief orbits. The mean relative orbital elements are used to describe the relative motion. The second-order relative motion equations are obtained by performing the Taylor expansion of the nonlinear relative motion equations about the chief orbit. The original thrust acceleration is approximated by a truncated Fourier series in terms of eccentric anomaly of the chief orbit. Based on the first-order state propagation and the thrust acceleration Fourier-series expansion, the second-order propagation of the mean relative orbital elements is obtained. Numerical simulations show that the state propagation errors of the second-order solution are reduced by 1 <inline-formula xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink"><tex-math notation="LaTeX">$\sim$</tex-math></inline-formula> 2 orders of magnitude when compared with the first-order solution.