On Sun-Synchronous Orbits and Associated Constellations Daniele Mortari, Matthew P. Wilkins, and Christian Bruccoleri Department of Aerospace Engineering, Texas A&M University, College Station, TX 77843, USA Abstract This paper contains the preliminary ideas on how to design orbits and constellations associated with extended Sun-Synchronicity concepts. The Multi Sun-Synchronous Orbits are introduced first. These orbits could be adopted by those Earth science missions whose aim is to study the local hour effects on the observed object. In particular, these orbits allow us to build the Multi Sun-Synchronous Constellations, whose mechanism is briefly introduced here. By applying the Sun-Synchronicity constraint to the Flower Constellations theory, the Sun-Synchronous Flower Constellations have been derived. For these constellations, however, the synchronization is referring to the relative trajectory in a reference frame that rotates with Earth s mean motion rather than to the Earth-Centered Inertial frame. Finally, by merging the concepts of Multi Sun-Synchronicity and the Flower Constellations, the Multi Sun-Synchronous Flower Constellations are obtained. Introduction This paper defines a set of satellite orbits associated with classic and novel Sun-Synchronicity concepts. The intent of this paper is to summarize and present the basic ideas. The complete analysis on the proposed orbits and constellations will belong to a full paper (Ref. [?]) that will contain the associated exhaustive study. In particular, the Multi Sun-Synchronous Orbits (MSSO), which generalize the Sun-Synchronous Orbits, are introduced. These orbits are built such that the difference between the orbit nodal precession rate and the Earth s mean motion is proportional to the orbit nodal precession rate or to the Earth s mean motion through a rational number. Based on these orbits, the idea of Multi Sun-Synchronous Constellations is presented. In this constellation, after every orbit nodal period, the constellation is shifted by a constant angle, such that the constellation appears identical with respect to the Sun direction. Then, within the recently proposed Flower Constellations theory (see Ref. [1]), it is also shown how to build the Sun-Synchronous Flower Constellations, with an adaptation of the definition for the Sun-Synchronous concept. The orbits of a Flower Constellations are compatible orbits, that is, all the spacecraft follow the same closed loop repeating trajectory in an Earth-Centered Earth-Fixed (ECEF) reference frame. In the Sun-Synchronous Flower Constellations, is the spacecraft relative trajectory in ECEF that track the Sun, rather than the inertial spacecraft orbit in Earth-Centered Inertial (ECI). Finally, by merging
the ideas of the Multi Sun-Synchronous Orbits and the Flower Constellations, the Multi Sun-Synchronous Flower Constellations have been devised. In order to proceed with the analysis some reference frames (see Fig. 1) have to be introduced. Figure 1: Reference Frames The first reference frame is the ECI reference frame [γ, τ, N], where γ is the vernal direction and N is the Earth s spin axis. The second frame is the Earth s orbit reference frame [x E, y E, z E ], where the x E -axis is aligned with the Sun-Earth direction and the z E -axis is the Earth s orbit normal (celestial pole). Finally, the third frame is the spacecraft orbit reference frame [x, y, z], that has the x-axis aligned with Earth-Spacecraft direction and the z E axis is the S/C orbit normal. The orientation of the [x E, y E, z E ] reference frame with respect to the ECI frame is identified by the orthogonal matrix A E = cos ϕ E sin ϕ E cos i E sin ϕ E sin i E sin ϕ E cos ϕ E cos i E cos ϕ E sin i E (1) 0 sin i E cos i E where C cos and S sin, and the orientation of the [x, y, z] reference frame in ECI has the form C ϕ C Ω S ϕ S Ω C i C ϕ S Ω + S ϕ C Ω C i S ϕ S i A S/C = S ϕ C Ω C ϕ Ω C i S ϕ S Ω + C ϕ C Ω C i C ϕ S i (2) S Ω S i C Ω S i C i Sun-Synchronous Orbits A Sun-Synchronous Orbit (SSO) is defined as the orbit whose normal (z axis) makes a constant angle with the Sun-Earth direction (x E axis), hence, the key property (z T x E = constant) of being synchronized with the Sun, as the name
implies, is derived. For circular orbits, this property implies that the local hour behavior is periodic in one orbit nodal period. The Sun-Earth direction rotates about z E with an angular velocity equal to the Earth s orbit mean motion n E = 2π T E 1.991 10 7 rad/sec (3) where the Earth orbit is T E 365.25 days. Therefore, the necessary and sufficient condition for a SSO is that it must rotate (precede) about z E with n E angular velocity. This cannot be achieved by natural motion. However, it is possible to approximate it by taking advantage of the J 2 effect that results in a precession of the inertial orbit about N. The nodal rate due to the Earth s gravitational coefficient, J 2, is Ω = 3 2 R 2 E a 2 (1 e 2 ) 2 n J 2 cos i (4) where R E = 6, 378.1363 Km is the Earth radius, a is the orbit semi-major axis, e is the orbit eccentricity, i is the orbit inclination, J 2 = 1.0826269 10 3 is a coefficient of an expansion representing the Earth s oblateness effect, while the expression of the orbit mean motion is n = µ a 3 where µ = 398, 600.4415 Km 3 /s 2 is the Earth gravitational parameter. Therefore, the Sun-synchronicity approximated condition is achieved by setting Ω = n E (5) Substituting Eqs. (3-4) into Eq. (5), we obtain the relationship cos i = 4.774 10 15 a 7/2 (1 e 2 ) 2 (6) Equation (5) does not fully comply with the Sun-Synchronicity definition (z T x E = constant) and the effect of this approximation are described in the following. While z rotates about N, the x E direction rotates about z E, and these two rotation axes are displaced by i E 23.5 (ecliptic obliquity). This implies, for instance, that the angle between z and x E (which should be constant), has a behavior that repeats in one year during which this angle experiences variations depending on the initial orbital condition in term of right ascension and inclination. Hence, the angle between z and x E would experience a minimum variation that should be eliminated by orbit maintenance control. However, this would imply continuous plane change maneuvers that are, in turn, unacceptable because they prove to be fuel cost prohibitive. For this reason, Eq. (5) is accepted as the (approximated) definition of the SSO, and the described unavoidable variations must, therefore, be accepted and tolerated. These variations affect the sub-satellite local hour. Figure 2 shows a typical local hour variation experienced in one year for a generic orbit (Ω 0 = 83.2 and i = 98 ).
Another important parameter for SSO is the angle λ between two successive ascending nodes. This angle, after each orbit nodal period, is represented by λ = ( Ω α E ) T Ω, where α E = 7.2921 10 5 rad/s is the Earth s spin rate and T Ω is the nodal period, that is defined[2] as 2π T Ω = n + Ṁ0 + ω where the rates of change for the initial mean anomaly and the perigee argument due to J 2 perturbation, have the expressions and (7) Ṁ 0 = 3R2 E J 2 4p 2 n 1 e 2 (3 cos 2 i 1) (8) ω = 3R2 E J 2 4p 2 n (5 cos 2 i 1) (9) respectively. Equations (4), (8), and (9), come from geopotential perturbation theory[2] by considering second order zonal effects, only. Therefore, the distance along the equator between successive node crossing points is ξ = R E λ. 4.5 12 4 Ω = 83.2099 11 min(max(lh)) 3.5 i = 98 10 Local Hour Variation 3 2.5 2 Local Hour Variations (in one year) 9 8 7 6 1.5 1 5 4 max(min(lh)) 0.5 0 50 100 150 200 250 300 350 400 Earth Longitude (deg) 3 0 20 40 60 80 100 120 140 160 180 Orbit Inclination (deg) Figure 2: Local Hour Example Figure 3: Local Hour Variations Finally, in accordance with the notation adopted in Fig. 1, the range of the sub-satellite point local hour h is [ 6 3 2 ] [ π cos 1 ( z T x E ) h 6 3 + 2 ] π cos 1 ( z T x E ) (10) where x E = cos ϕ E sin ϕ E cos i E sin ϕ E sin i E and z = sin Ω sin i cos Ω sin i cos i (11)
In particular, the maximum variation of h is h max = 24 π cos 1 ( z T x E ) (12) Figure 2 shows a typical local hour variation during a year for a generic orbit (Ω 0 = 83.2 and i = 98 ), while Fig. 3 summarizes the results obtained by using Eq. (12), in term of minimum of maximum and maximum of minimum values of h max experienced in the entire range of orbit inclinations. The SSO are typical of Earth observation missions that desire to minimize the effects of different solar ray inclinations (local hour) on the observed object (clouds, ice, specific atmosphere component or layer, etc). A complementary Earth Science can be accomplished when the science target is precisely the measurement of the effect(s) of the local hour on the specific observed object. The Multi Sun-Synchronous Orbits, introduced in the next section, have been developed specifically to address this complementary Earth science mission. Multi Sun-Synchronous Orbits Equation (5) characterizes the SSO: the Sun-Earth direction x E and the orbit normal direction z rotate about z E and N, respectively, by the same angular speed ( Ω = n E ). This is a continuous property. By extension, the Multi Sun-Synchronous Orbits (MSSO), can be defined based on a discrete time characterization. This paper introduces two distinct definition of the MSSO. First, consider multiplying Eq. (5) through by the orbit nodal period T Ω. We obtain T Ω Ω = TΩ n E, which allows one to view the Sun-Synchronicity concept from a different perspective but in an equivalent way. This implies that, in the time T Ω, the directions x E and z rotate by the same angle. Based upon this view, the MSSO are then defined as satisfying the equation where ξ = I n I d T Ω ( Ω n E ) = 2π ξ (13) is the MSSO rational parameter, and I n and I d are two independent integer parameters (positive or negative). In fact, this equation states that, in the orbit nodal period T Ω, the relative rotation angle between x E and z is changed by a rational amount of the 2π angle. Equation (13) generalizes Eq. (5) which, in turn, can be obtained by setting I n = 0. For any rational value of ξ, Eq. (13) represents the following relationship among a, e, and i. A cos 2 i + B cos i + C = 0 (14) where, setting χ 3R2 E J 2, the expression for A, B, and C, become 4 A = (5 + 3 1 e 2 ) ξ B = 2 [ a 2 (1 e 2 ) 2 C = ] 1 e χ 2 1 ξ + n E a 7/2 (1 e 2 ) 2 χ µ (15)
Equation (14) is a quadratic function in cos i that can be numerically solved. For the case of circular orbits, the expressions for A, B, and C reduce to A = 8 ξ B = 2 ( ) a 2 C = χ 2 ξ + n E a 7/2 χ µ (16) 180 ξ=1/10000 160 140 ξ=1/1000 ξ=1/2000 ξ= 1/10000 ξ= 1/3000 ξ= 1/2000 120 Inclination (deg) 100 80 60 ξ= 1/1000 40 ξ= 1/1000 20 0 6 8 10 12 14 16 18 20 22 24 Semi major axis ( 10,000 Km) Figure 4: Multi Sun-Synchronous Circular Orbits: Type I Figure 4 shows the relationship between a and i for various values of the integer parameters I d and I n, for the case of multi Sun-Synchronous circular orbits of this kind. The correct solution is cos i = ( B B 2 4AC)/(2A). Alternatively, a second way to introduce the MSSO is by the relationship K d Ω = Kn n E (17) where K n is a positive integer parameter and K d is an independent (positive or negative) integer parameter. This equation states that, in K d Earth orbit periods, the orbit nodal line precesses (clockwise or counter-clockwise) K n
times. Equation (17) also generalizes Eq. (5) that, in turn, can be obtained by setting K n = K d = 1. Equation (6) can be rewritten for the MSSO obtaining cos i = 4.774 10 15 K n K d a 7/2 (1 e 2 ) 2 (18) which enlarges the possibilities given by Eq. (6). Notice that the parameter η = K n orbit. K d can be negative, which avoids the SSO constraint of a retrograde 180 η=2/1 η=1/1 η=3/1 η=3/2 η=2/3 η=1/2 η=1/3 160 140 120 Inclination (deg) 100 80 60 40 20 η= 3/1 η= 2/1 η= 3/2 η= 1/1 η= 2/3 η= 1/2 η= 1/3 0 6 8 10 12 14 16 18 Semi major axis ( 10,000 Km) Figure 5: Multi Sun-Synchronous Circular Orbits: Type II Figure 5 shows the relationship between a and i for various values of the integer parameters K d and K n, for the case of multi Sun-Synchronous circular orbits of the second kind. Multi Sun-Synchronous Constellations Let us consider a satellite constellation consisting of N s satellites whose orbits have identical inclination, argument of perigee, eccentricity, and semi-major axis. As for the right ascensions, these orbits are evenly and equally separated
by the angle Ω = 2 π (19) N s In consideration of the time-varying parameters, say the mean anomalies for each spacecraft, there are no particular constraints. However, particular attention is given to those values complying with the theory of Flower Constellations[1]. Based upon these assumptions, we can introduce the Multi Sun-Synchronous Constellations (MSSC), as follows. In the time (J d T Ω ), where T Ω is the orbit nodal period and J d is a given integer number, the relative angle between node and Sun direction increases/decreases by a given angle. The MSSC are defined such that this angle is exactly m times the angle Ω given in Eq. (19), that is, is (J n Ω), where J n is another integer parameter. Now, since the relative velocity between the node and Sun directions is (n E Ω), then we can write the relationship J d T Ω (n E Ω) = J n Ω = J n 2 π N s (20) that defines the Multi Sun-Synchronous Constellations. By introducing the MSSC rational parameter η = J d J n (21) By comparing Eq. (13) with Eq. (20) we can derive the relationship between the integer parameters I d = N s J d and I n = J n (22) Equation (20) constitutes a relationship between inclination, semi-major axis, eccentricity, and the MSSC rational parameter η. Sun-Synchronous Flower Constellations The Flower Constellations (FCs), introduced in Ref. [1], are built using compatible orbits. A compatible orbit (often called repeated ground track orbit) has the orbital period such that the relative trajectory in an Earth-Centered Earth-Fixed (ECEF) system of coordinates constitutes a repeated closed-loop trajectory 1. The relative trajectory can be designed such that it repeats after the satellite completes N p revolutions (orbit nodal periods) over N d sidereal days. If T r is the period of repetition, then T r = N p T Ω = N d T ΩG (23) where T ΩG is the nodal period of Greenwich, which is defined by Carter in Ref. [?] as T ΩG = 2π α E Ω (24) 1 The name Flower Constellation is due to the petal-like shape of the closed-loop relative trajectories as viewed from an ECEF rotating frame.
In the theory of the FCs, which includes the J 2 gravitational perturbation, it is shown that there are a set of admissible orbital positions (admissible mean anomalies) that allow each spacecraft to follow the same identical closed-loop relative trajectory in ECEF. The FCs are built on this property and on a particular procedure (phasing rule) to distribute the satellites on the admissible positions. In particular, in an axial-symmetric FCs, the node lines of each satellite are equally displaced along the equatorial plane. The Sun-Synchronous Flower Constellations (SSFCs) that are introduced here are FCs whose orbits are compatible with respect to a fictitious Earth that spin with angular velocity equal to the Earth s mean motion n E. This is the description of an Earth that is gravity gradient stabilized with respect to the Sun. In other words, this fictitious body is fixed on the [x E, y E, z E ] reference frame. This implies a new way to see the Sun-synchronicity concept. The inertial spacecraft orbit is not synchronized with the rotation of the solar ray direction (Sun-Synchronous classical definition), but rather it is the spacecraft relative trajectory (with respect to a ECEF rotating with angular velocity n E ) that is synchronized with the Sun. In order to build such Sun-Synchronous relative trajectories, the Earth spin rate α E must be substituted by the Earth mean motion n E in Eq. (24). Apogee t = 8 months t = 4 months N x E xe xe n E xe t = 0 months t = 12 months Figure 6: SSFC Figure 7: MSSFC Figure 6 shows an isometric view of a 3-1 Flower Constellation where one apogee has been aligned with the Earth-Sun axis (x E ). The entire relative path is set to rotate about the Earth s polar axis (N) such that that particular apogee stays aligned with x E. Hence, this is a Sun-Synchronous Flower Constellation. Multi Sun-Synchronous Flower Constellations The concept of Multi Sun-Synchronicity can be applied to the Flower Constellations theory, thus obtaining the Multi Sun-Synchronous Flower Constellations. The synchronization is that the FC relative trajectory rotates in the
inertial space with an angular velocity such that, every time a petal (relative apogee) becomes aligned with the Sun rays direction then a spacecraft of the constellation passes to the aligned relative apogee. Figure 7 shows a polar view of a 3-1 Flower Constellation where one apogee has been initially aligned with the Earth-Sun axis (x E ). The entire relative path rotates about the Earth s polar axis (N) such that the petals (relative apogees) will sequentially align with x E at a prescribed times (0, 4, and 8 months). Hence, this is a Multi Sun-Synchronous Flower Constellation. Conclusion In this paper, we have presented the summary of the basic ideas on how to design orbits and constellations associated with extended Sun-Synchronicity concepts. In particular, the Multi Sun-Synchronous Orbits have been introduced. These orbits have been devised for Earth science missions that want to study the local hour effects on the observations. Built on these orbits, the Multi Sun- Synchronous Constellations have been proposed. Then it has been shown how to apply the Sun-Synchronicity constraint to the Flower Constellations theory, and the Sun-Synchronous Flower Constellations have been derived. The main characteristics of these constellations is that the Sun-synchronization is associated with the motion (in ECI) of the spacecraft relative trajectory that appear fixed into a reference frame that rotates with Earth s mean motion, while the classic Sun-Synchronicity concept is a particular precession of the inertial orbit. The paper then presents the Multi Sun-Synchronous Flower Constellations that are obtained by exporting the Multi Sun-Synchronicity concept into the Flower Constellations theory. All the ideas presented in this paper are preliminary. A detailed and complete analysis of these orbits and constellations is presently under preparation and it will appear in a future paper. References [1] Mortari, D., Wilkins, M. P. and Bruccoleri, C., The Flower Constellations, in Advances in the Astronautical Sciences, Texas A&M University, College Station, TX: Paper AAS 03 274 of the John L. Junkins Astrodynamics Symposium, 2003 [2] Vallado, D. A., Fundamentals of Astrodynamics and Applications, vol. 2, McGraw Hill, New York, 2001