In this paper, we consider how near-term solar sails can be used in formation above the ecliptic plane to provide platforms for accurate and continuous remote sensing of the polar regions of the Earth. The dynamics of the solar sail elliptical restricted three-body problem are exploited for formation flying by identifying a family of periodic orbits above the ecliptic plane. Moreover, we find a family of 1-year periodic orbits in which each orbit corresponds to a unique solar sail orientation using a numerical continuation method. It is found through a number of example numerical simulations that this family of orbits can be used for solar sail formation flying. Furthermore, it is illustrated numerically that solar sails can provide stable formation-keeping platforms that are robust to injection errors. In addition, practical trajectories that pass close to the Earth and wind onto these periodic orbits above the ecliptic are identified. Nomenclature a = semimajor axis a s = solar sail acceleration, m=s 2 a x , a y , a z = components of solar sail acceleration such that a s a x ; a y ; a z T , m=s 2 e = eccentricity of the Earth's orbit about the sun e 0:0167 f = true anomaly of the Earth about the sun, rad G = universal gravitational constant, 6:673 10 11 m 3 kg 1 s 2 m 1 = mass of the sun, 1:98892 10 30 kg m 2 = mass of the Earth, 5:9742 10 24 kĝ n = unit normal to the sail p = semilatus rectum t = dimensionless time r = position vector of the solar sail in the rotating frame, astronomical units (AU) r 1 = position of the solar sail with respect to the sun, AU r 2 = position of the solar sail with respect to the Earth, AU V = velocity vector of the solar sail X, Y, Z = rotating coordinate frame, AU x, y, z = rotating-pulsating coordinate frame, AU = solar sail lightness-number ratio of solar sail radiation pressure acceleration to solar gravitational acceleration 0:05 = the angle that the sail normal makes with the y axis, rad = the angle that the sail normal makes with the x axis, rad = distance between the sun and the Earth, AU = dimensionless mass of the Earth, 3 10 6 ! = angular velocity vector of the rotating frame