1. Introduction
Missions to Lagrange points have been effective because they provide opportunities for flexible, low-energy trajectories. Lagrange points are equilibrium points in the space between two celestial bodies where the net gravitational force on an object at that point is zero. The five stable Lagrange points are L1, L2, L3, L4, and L5. In particular, L1 and L2 points are of interest because they offer a unique location in space where a spacecraft can maintain a stable position relative to the two celestial bodies. These points can be used for a variety of scientific and practical applications, such as satellite constellations, thermal observation, and studying the environment of celestial bodies.
Halo orbits and Near Rectilinear Halo Orbits (NRHO) are types of orbits that are related to Lagrange points. Halo orbits are defined by their proximity to a Lagrange point and their ability to remain in the vicinity of that point for an extended period of time, allowing for a broad range of scientific observations and technological applications. NRHO orbits are a specific type of Halo orbits that do not coincide with resonance with the motions of celestial bodies. These orbits are called "near-rectilinear" because they are close to being a straight line. NRHO orbits have the advantage of being less affected by the perturbations caused by other celestial bodies, making them useful for long-term missions.
One example of the application of this type of orbit is the James Webb Space Telescope (JWST or Webb), a joint NASA/ESA/CSA mission that launched an infrared-optimized 6.5 m space telescope in December 2021. The JWST will allow for a range of astronomical and cosmological studies, such as the observation of the first stars, the creation of the first galaxies, and a detailed analysis of potentially habitable exoplanets [
1]. Another example is the Near-Earth Object Surveillance Mission (NEOSM), which is scheduled to be launched by NASA in 2025. The NEOSM will use an infrared telescope to survey the solar system for potentially hazardous asteroids and will be conducted by the NEO Surveyor spacecraft, which will be stationed at the Sun–Earth L1 Lagrange point. This location allows it to look close to the Sun and detect objects inside Earth’s orbit that would otherwise be difficult to spot from Earth or low Earth orbit [
2].
Moreover, in recent years, several proposals have been made for using libration point orbits in the Earth–Moon system for communication and navigation services, particularly around the collinear L1 and L2 points. For example, Grebow et al. [
3] have investigated the use of libration points for lunar South Pole coverage. Hamera et al. [
4] analyzed the feasibility of using small satellites on Halo orbits for a lunar constellation. A lunar positioning system has also been designed using Lissajous [
5] or Halo orbits [
6] around both L1 and L2 points. Last, but not least, a future mission will deploy in a Moon Halo orbit, the Lunar Gateway: a space station that is a key part of NASA’s Artemis program. The Lunar Gateway will serve as a science laboratory, a testing ground for new technologies, and a stepping stone for deep space exploration, including toward Mars.
For all these missions, and for future ones, Orbit Determination (OD) is a crucial element in order to meet the requirements and the goals of the missions. Conventionally, Earth-based tracking stations, such as the Deep Space Network (DSN), are used for tracking spacecraft. However, this practice has several limitations, especially if the number of simultaneous missions to be tracked increases. In fact, it has a capacity limit in terms of spacecraft that can be served simultaneously; it is costly and cannot support direct full coverage of a celestial object (for example, the case of spacecraft on the dark side of the Moon).
In recent years, a study introduced the Linked, Autonomous, Interplanetary Satellite Orbit Navigation (LiAISON) architecture as an alternative way for tracking spacecraft deployed in orbits around libration points, such as the Halo orbits. LiAISON is a concept of OD that uses relative inter-satellite measurements to estimate the absolute states of two or more spacecraft simultaneously [
7].
An example of the use of this technique is the Cislunar Autonomous Positioning System Technology Operations and Navigation Experiment (CAPSTONE) mission [
8]. It has started its operational activities, and it will serve as a pathfinder for navigation and operations in the same NRHO that will be utilized by NASA’s Lunar Gateway. In addition to this objective, CAPSTONE contains a dedicated flight board for the Cislunar Autonomous Positioning System (CAPS), a framework for autonomous or near-autonomous navigation in cislunar space. The navigation architecture is based on inter-satellite radiometric data transferred between the CAPSTONE spacecraft, which is deployed in the NRHO orbit, and the Lunar Reconnaissance Orbiter (LRO), which is deployed in a low lunar orbit [
9]. The combination of an Earth–Moon L2 Halo orbit and a low-lunar orbit has been shown to yield enough geometric and dynamic asymmetry such that the states of both spacecraft are observable.
The basic working principle proposed by the LiAISON study is based on the assumption that the Halo orbit has a unique size and shape because of the strong asymmetry of the three-body force field. This means that a spacecraft on these types of orbits can track a second spacecraft using inter-satellite measurements and determine their absolute positions and velocities without any Earth-based tracking or any other mathematical constraints.
In this work, we propose and evaluate, for a given case study, a method for Orbit Determination and Time Synchronization (OD&T) of spacecraft deployed in the same Halo orbit exploiting an inter-satellite link (ISL) for dual one-way measurements of the range and range rate between spacecraft. The proposed approach not only gives the possibility to have a fully autonomous orbit determination but also a total independent time synchronization of the spacecraft, avoiding the use of the DSN, which can be used only for periodic checks and calibrations.
For example, in the case of a single spacecraft mission, such as the JWST or any other scientific observation mission, exploiting only one spacecraft, additional small and low-cost orbiters (called “floating anchors” hereafter) can be deployed from the same vector in the same orbit to perform the inter-satellite measurements to be used for autonomous OD. The mission’s main spacecraft, exploiting the measurements (to and from the floating anchors), will be able to solve the OD problem autonomously and without any operation from the ground.
Moreover, for other types of missions, such as for navigation systems like the ones proposed in [
10,
11] or telecommunication systems in which, typically, more spacecraft are needed to cover a given region, the additional spacecraft can also be used to provide these services, being telecommunication relays to/from the Earth or other celestial bodies, or providing ranging signals for navigation purposes to other spacecraft or rovers. Finally, they can also be equipped, if needed, with instruments for remote sensing such as SAR, laser altimeters, or others, exploiting the precise orbit of the spacecraft.
After the description of the proposed OD and Timing procedure, the possible performances are evaluated for application in a NAV/COM constellation. In particular, in the case study, the possibility of adding two floating anchors to the Lunar Gateway, and to exploit this system of three spacecraft to deliver navigation and telecommunication signals to the Moon’s South Pole, is investigated.
The paper is arranged as follows: in
Section 2, the system concept and the inter-satellite measurements approach are introduced. In
Section 3 and
Section 4, the Orbit Determination and Timing algorithms are explained in detail.
Section 5 introduces and evaluates the case study, and discussions and conclusions follow in
Section 6.
2. System Concept and Dual One-Way Measurements
As stated in the Introduction, there have been many missions that exploited Lagrange orbits because they offer a stable, convenient, and cost-effective location for missions that require a long-term presence in space or a fixed view of a particular region of the celestial sphere. In these particular orbits, it is possible to exploit relative tracking measurements to estimate relative and absolute orbit states for all participating spacecraft simultaneously if at least two spacecraft are deployed [
7]. Usually, these measurements are conducted exploiting different communication bands such as X-band or K-band, or exploiting inter-satellite optical links.
In general, if the additional spacecraft are not used for other services, it can be assumed that only the main spacecraft have an onboard system able to process the measurements, allowing small, low complexity and cheap floating anchors, used for inter-satellite measurements.
Figure 1 illustrates an example in which two spacecraft are deployed in a Lagrange orbit.
Different inter-satellite measurements can be performed by exploiting an ISL. The most used technique is the conventional two-way ranging that is also used by ground stations: it consists of sending an up-link signal coming from the TX-satellite (which could be a sequence of discrete tones or a pseudo-noise sequence); this signal is then received and re-transmitted back by the RX-satellite and received by the TX satellite, and this exchange allows the computation of the two-way propagation time [
12]. Another technique is the telemetry-based ranging that provides a round-trip time derived from the telemetry stream [
13,
14]. This method does not require a return ranging signal, but it provides, as part of the telemetry stream, the delay between the acquired ranging signal and the start of the next telemetry frame. Both techniques provide a round-trip delay time solution at the signal source. Finally, an alternative technique is based on time transfer, as in the CCSDS Proximity-1 Space Link Protocol [
15]: users can exchange epochs between satellites and derive the round-trip time using dedicated algorithms [
16]. All of these techniques can be implemented, in principle, in any communication band, including the optical one. Finally, another effective way for two-way ranging can be the use of a satellite laser ranging technique exploiting retro-reflectors on the spacecraft [
17].
In this work, the use of a different technique is proposed: the so-called dual one-way. It has the advantage of producing not only an estimation of the satellite range but also an estimation of the satellite clock biases.
Referring to
Figure 2, each spacecraft is equipped with a transmitter and a receiver to exchange ranging signals and with a reference clock to measure the Time of Transmission and the Time of Arrival of these signals. The working principle is mutated from the one used in GNSS: pseudo-random codes are transmitted from one spacecraft and received in the other exploiting a Delay Lock Loop (DLL) to measure the signal Time of Arrival (TOA) and a Phase Lock Loop (PLL) to measure the signal Doppler shift (Frequency Of Arrival-FOA). These measurements are used to estimate the range, the range rate, and the clock bias at regular time intervals (
); the measurements are then accumulated for a given time span (
T) and used to solve the OD&T problem in the main spacecraft.
The orbiting spacecraft exchange data and measurements exploiting the ISL and, if they have a precise OD&T, can be used to transmit to the celestial body Navigation and Telecommunication signals.
In more detail, the range between two satellites is defined as the linear distance between the spacecraft and can be measured by exploiting the Time Of Arrival (TOA), if the TOA measurements taken on the anchors are transmitted back to the main spacecraft. In detail, let us have two spacecraft
i and
j in position
and
and let us assume a common reference time
t. The spacecraft can have different clock biases with respect to the reference time;
,
, and the time of flight of signal can be measured by the use of the TOAs taken by both of the satellites, in the following way:
where
is the Time of Transmission on satellite
j (encoded in the transmitted signal), and
is called “pseudorange” from satellite
i to satellite
j.
Finally,
is the measurement error. The pseudo-ranges,
and
do not coincide with the real range, but, if both of the spacecraft produce a measurement, it is possible to compute the range making the average of the two pseudo-ranges because, by definition,
:
Moreover, if the difference between the pseudo ranges also computed, an estimation of the clock bias between the two spacecraft can be obtained:
As mentioned before, the time of arrival is measured in the receiver exploiting a Delay Locked Loop as in the case of GNSS [
18].
Exploiting the same signals, it is also possible to measure the range rate between two spacecraft, with
i and
j exploring the Doppler effect. Let us recall that the Doppler shift for a radio frequency signal is given by:
where
is the radial velocity of the spacecraft
j with respect to the spacecraft
i,
c is the speed of light, and
is the carrier frequency of the data link. It follows that
can be seen as the derivative of the geometrical distance of the satellites and that the Doppler shift can be expressed as follows:
where
represents the measurement error in measuring the Doppler shift of satellite
i from satellite
j. The Doppler shift can be measured by the Phase Lock Loop implemented in the receivers [
18].
Concerning the measurement errors, measuring the TOA of a given signal is a well-known problem, and various methods have been proposed in the literature to find the Cramer–Rao limit, and the most used expression is given by [
19]:
where
E is the signal energy,
is the noise power spectral density, and
is the Gabor bandwidth, which is defined as:
where
is the power spectrum of the received signal.
In the same manner, a similar expression can be formulated for the measurement of the Frequency of Arrival (FOA) of the signal [
19]:
where
is the effective duration, which is defined as [
19]:
where
is the received signal.
In the case of pseudo-random code and Delay Look Loop (DLL), the TOA error model can also be refined by taking into account many of the receiver and signal parameters, by the following expression [
18]:
where
is the code loop noise bandwidth,
is the double-sided front-end bandwidth,
is the predetection integration time,
is the chip period,
is the spreading code rate,
is the carrier to noise power ratio, and
D is the early-to-late correlator spacing.
Moreover, Ref. [
20] demonstrates that using the PLL filter to estimate the Doppler frequency and the range rate gives the following error due to thermal noise (for the range rate):
where
is the wavelength of the signal,
is the coherent integration time,
is the carrier to noise ratio, and
is the Doppler bandwidth defined as:
where
is the frequency noise impulse responses transformed into the Z-domain. The Doppler bandwidth
summarizes in a single parameter the ability of the tracking loop to produce smooth frequency estimates.
Having these formulations, it is possible to find the range and range-rate measurement errors knowing the receiver and signal parameters. Referring to the literature, in [
21], the range measurement error between GPS satellites is assumed to be 1 m by using ranging measurementswith a maximum ranging distance of 50,000 km; Ref. [
22] proposes a method for relative navigation for microsatellites using inter-satellite measurements, and the range error is assumed to be at a maximum of 10 cm for a maximum ranging distance of 10 km. Another example can be found in [
23], where a precise ranging and time synchronization system was implemented based on dual one-way ranging with a precision of the ranging measures of 1.038 m. Ref. [
24] proposes a solution for relative OD using inter-satellite range measurements, and the ranging measurement error is assumed to be Gaussian with a standard deviation of 1 m.
Refs. [
18,
25,
26] explain in more detail the range measurement assumptions and the relative models, and it is shown that the measurement error always depends on the carrier-to-noise ratio,
, remaining below 1 m for most cases.
Concerning the range rate, in [
27], the measurement accuracy is assumed to be 0.01 mm/s, and in [
28], an autonomous navigation system for satellites deployed in Martian orbits was proposed, and the worst accuracy in range rate was assumed to be below 4
m/s. The inter-satellite link, in this case, was in the X-band, and the Carrier-to-Noise ratio was assumed to be equal to 30 db-Hz.
In this work, we assume the use of the Ka-band.
(
GHz) for the inter-satellite link. It is compliant with the standards defined in [
29], allowing a range error of 1 m, obtained, for example, with a typical DLL code noise bandwidth
of 0.2 Hz, an integration time
T of 1 s and
dB-Hz, with the following other typical parameters:
MHz,
Mcps, and
.
The range-rate error is assumed to be below 0.06 mm/s, and it can be obtained with the PLL with a Doppler bandwidth of 1 Hz, an integration time of less than five seconds, and a carrier-to-noise ratio of 35 dB-Hz.
In the following sections, the proposed method to use these measurements for a precise OD&T is reported.
3. Orbit Determination with ISL Measurements
Let us assume to have a system of orbiters in a Halo orbit. This system can be easily described by exploiting the Circular Restricted Three-Body Problem (CRTBP) where two massive bodies are considered orbiting around their mutual center of mass, and the spacecraft with infinitesimal masses experience forces due to the gravitational influence of both bodies simultaneously. All the bodies are approximated as point masses. Moreover, to simplify the notation, let us assume to use a reference frame centered in the two celestial body center of gravity and keep the x-axis pointed towards the second celestial body, obtaining a frame rotating with constant angular velocity, equal to the mean motion of the system.
Finally, let us scale the system in such a way that the distance between the two main bodies, as well as the angular velocity, are set to 1 (see
Figure 3).
In this coordinate frame, it is possible to define the state vector for any spacecraft (containing its position and velocity) as:
or, in the general case of
N spacecraft as:
Moreover, in this reference frame for the generic spacecraft
i, the equations take the following form: [
30]:
where
is the potential function:
and
and
and
are the celestial body masses.
In this particular frame, a batch-weighted least square orbit determination algorithm [
30] can be used to estimate the initial state of all the spacecraft.
The differential equation of the satellite motion, , can be derived from the motion law and can be used to propagate the initial state to any possible time t by integration. It follows that, knowing the initial state of the system, it is possible to derive the state, , for any possible time t.
In the same way, having a reference solution (that is a guess with an error of the initial state), it is possible to also compute a reference state for the time t, .
Finally, for any time
t, it is also possible to define the expected measurement vector
for the reference solution and the real measurement vector for
. The vector
contains the set of the range and range-rate measurements taken at time
t as defined in Equation (
3) and in Equation (
6). Making the hypothesis conduct the measurement process with an update time of
for a Time span of
T, it is possible to compute the state deviation,
, and the observation residual,
, for each measurement instant
:
The state deviations,
, can be mapped back to the initial epoch,
, with the state transition matrix,
, as follows:
The state transition matrix is obtained by numerical integration, and its differential function is:
where:
where []* means computed on the reference trajectory.
Considering the case of
N spacecraft,
has the following formulation
, where:
The same approach can be used to map back all the measurements to the initial time
, exploiting the observation residuals, defining the
matrix as:
where
is the function that relates the observations to the state vector in the following way:
as given in Equations (
3) and (
6).
Therefore, it is possible to move the set of observations at time
back to the initial epoch
defining the matrix
:
When the total number of observations is larger than the number of unknowns, a solution for
can be found by exploiting the normal equation:
where
is the information matrix, and
is the reduced residual vector, where
is the block matrix containing the observation–state relationship matrices
for all the epochs, and is given by:
and
is the weighting matrix:
assuming that the observation errors can be modeled as white Gaussian noise, and the standard deviation in the observation noise can be denoted as
,
.
4. Satellite Synchronization
As introduced before, by exploiting the dual one-way technique, it is possible to compute the clock bias between the two spacecraft at the time of the measurement.
Unfortunately, computing only the time bias for a given specific time is not sufficient to model the deviation between two clocks. A more complex model that considers deterministic and random components of the clock deviation, such as bias, drift, and drift rate, is usually needed. A common solution is the use of stochastic differential equations.
The clock time and frequency errors for a time epoch
k can be derived from the previous one by exploiting, for example, the model proposed in [
31] where two main clock noises are considered: White Frequency Modulation (WFM) and Random Walk Frequency Modulation (RWFM). These two noises generate on the clock’s phase a Wiener noise (integral of the WFM) and an integrated Wiener noise (integral of the RWFM):
where
represents the clock phase deviation (time Bias),
represents the clock frequency deviation, and
and
are the clock deterministic drift terms. The vector
is a bivariate random variable with zero mean and variance that depends on two parameters:
, which is the diffusion coefficient of the phase noise, and
, which is the diffusion coefficient of the frequency noise; it represents the covariance matrix of the error related to the specific clock. By setting the
,
,
, and
, it is possible to set the noise model representing the desired type of clock to be represented for a long time.
For short-time modeling, this model can be replaced by a simpler one that takes into account only the deterministic components of the clock deviation, which are the time bias, the drift, and the drift rate using the following polynomial representation for the bias of each satellite clock at given epoch (
) [
32,
33]:
where
is the clock offset,
is the clock drift, and
is the clock drift rate. This simpler model is commonly used to represent the clock error in the GNSS system and can be accurate enough to represent the clock error if it is frequently updated. For example, in the GPS system, the parameters of the clock are updated every two hours [
34].
In this work, this second approach is proposed and, considering the main spacecraft clock as the reference one, the clock parameters to compute the bias with respect to the main spacecraft for the anchor satellite i can be collected in the vector .
Recalling that, for each time epoch
, the TOA measurements are able to produce an estimation of the time bias
between the spacecraft clocks as derived in Equation (
4), for each spacecraft, it is possible to estimate
by linear regression or a Least Square solution, inverting the system of equations:
where
K is the number of epochs used for the estimation.
It must be noted that, as mentioned before, the proposed model works for a short time, and it is not possible to use the estimated parameters for long periods: the parameter estimation must be repeated, after a given time, for example, every 2 or 3 h.
It follows that the Time Synchronization and Orbit Determination processes, also if based on the same set of measurements, may have the same renewal time () but different time spans.
In
Figure 4, a schematic representation of the whole algorithm is reported. Every
, a set of TOA and FOA measurements, is taken and collected from the main spacecraft. The measurements are used to derive the range, range rate, and clock biases between spacecraft. At the time
, the range and range rate accumulated from the last epochs in
are used to derive the state deviation and the measurement residuals. Having all the epochs of the selected time span, the normal equation can be solved and the state of the spacecraft at time
can be estimated. In the same way, the clock parameters for the epoch
are estimated with the last
epochs. After that, the orbits and clocks parameters can be used to propagate the satellite’s position and time to the future epochs, and the sliding windows can be moved to the epoch
.
5. Evaluation for Lunar NHRO Orbiters for NAV/COM Applications
One of the next missions that will be deployed in a Moon Halo orbit is the Lunar Gateway. The Lunar Gateway will be a science laboratory, a test bed for innovative technologies, and represents a stepping stone to deep space exploration, especially toward Mars [
9]. The current baseline orbit for the Gateway is an NRHO with a period of approximately 6.56 days. It exhibits a 9:2 resonance with the lunar synodic period, completing nine revolutions around the Moon every two synodic months. The mean radius at perilune is approximately 3400 km, while the orbit extends to a radius of approximately 71,000 km at apolune [
35].
The case study proposed and evaluated here is a system of three spacecraft: the main one, which can be the Gateway, and two additional anchor orbiters for a totally autonomous and precise OD&T and for providing Navigation and Telecommunication services on the Moon’s surface, and in particular on the Moon South Pole, where the future Artemis missions are planned [
36].
The idea is to test the possibility to add a reduced number of small and inexpensive spacecraft in the same orbit as the Lunar Gateway to be used as relays for communication with the Moon’s Surface and for the transmission of ranging signals (as in the case of terrestrial GNSS system) to provide a navigation service. The additional orbiters can be equipped only with a telecommunication payload for communication with the Moon’s surface (and for transmitting ranging signals) and with ISLs able to perform the TOA and FOA measurements with the Gateway and, in principle, can be deployed from the same vector used for the Gateway components or from the Gateway itself. All of the OD&T operations are performed on the Gateway.
The constellation for the test is composed of the Gateway and two additional spacecraft and is represented in
Figure 5; all the details about the constellation are reported in
Table 1.
As can be seen from
Figure 5, the spacecraft are not equally distributed in the orbit; this configuration was chosen because it has good performance in terms of Navigation and Telecommunication service availability on the South Pole of the Moon.
Figure 6 reports the elevation angle for each spacecraft for all the orbiting period, for different epochs, as seen from the South Pole. It is possible to note that there is always at least one satellite in the view from the South Pole that can be used for telecommunication purposes with an elevation angle higher than 40 degrees.
Moreover, three satellites are in view for more than 95% of the time with an elevation angle higher than 10 degrees and can be used as ranging sources for navigation purposes. Having three synchronized satellites in view, it is possible to have an estimation of the user position if it is imposed to be on the Moon surface exploiting a Digital Elevation Model (DEM) or, for flying users, exploiting the DEM together with an altimeter as shown, for example, in [
37].
Simulations of the OD&T process are performed using the previously introduced orbit determination algorithm and considering the errors reported in
Table 2, exploding different time spans (between one and two orbital periods) and measuring intervals (between 15 and 60 min) for two system conditions: with (RR) and without (NORR) range-rate measurements. These two different configurations are compared to understand if the simpler configuration exploiting anchors without the ability to perform PLL-based range-rate measurements is also able to allow OD&T performances compatible with the stringent requirements to provide a navigation service.
Regarding the measurement errors reported in
Table 2, they are the typical measurement performances that can be obtained using a receiver implementing DLL and PLL and referring to Equations (
11) and (
12) with the previously mentioned parameters.
Monte Carlo runs were carried out with an initial error on the reference orbit generated for each run with a Gaussian distribution, zero mean, and a standard deviation of 100 m for the position error and 0.1 mm/s for the velocity error. The error on the inter-satellite range measurements was considered Gaussian distributed with 1 m of standard deviation for each measurement session, and the error on the range rate is considered with a standard deviation of 0.06 mm/s.
To evaluate the performance of the Orbit Determination process, the Distance Root Mean Square (DRMS) of the estimated trajectories with respect to the nominal one was computed: the estimated,
, and the nominal,
, trajectories were propagated for three hours, and the norm of the displacement vector (distance) for each epoch is used to compute the DRMS for the given spacecraft:
where
m indicates the
m-th satellite, and
l indicates the epoch of the propagation. Finally, the average of the DRMS error for all the spacecraft (Gateway and anchors) was considered:
For comparison, a classic orbit determination process exploiting an Earth-based measuring station performing the same measurements (range and range rate) was also simulated for the same time spans and measurement rate (measurements are taken only from the ground, without the use of any ISL).
Examples of results for a time span of 9.89 days (1.5 T) and measurement time interval of 30 min are reported in
Figure 7 and
Figure 8, where the histograms for the DRMS orbit error (in meters) for 500 Monte Carlo runs are reported for the two cases: with and without the range-rate measurements, respectively, and compared with the ground-based approach. It is possible to note that, with this configuration, the obtained performance is similar to the one obtained with the classical Earth-based OD method. Moreover, the results show that sub-meter positional DRMS error can be obtained by the use of range and range-rate measurements with a totally autonomous approach.
In
Table 3, the 90th percentile for the DRMS errors computed for different simulations, varying the time span and measurement intervals, are reported. Apart from the first value of the time span, the proposed Autonomous OD gives similar or better results than the OD from the Earth and with a sub-meter DRMS error for a large part of the tested cases.
Concerning the synchronization, the performances are calculated with the same approach: the clock was propagated with the nominal and estimated parameters, and the RMS clock bias was computed for different epochs in a time span of three hours.
The obtained results are compared with the a posteriori synchronization method based on GPS satellites proposed in [
38] in which, after an independent orbit determination, the GPS satellites are used for Time Transfer to the Moon orbiting spacecraft.
To mitigate the simplified clock model errors, the time span for estimating the clock parameters was chosen to be equal to 3 h.
In addition, in this case, the results are encouraging and always better than the ones obtained in [
38].
Figure 9 and
Figure 10 report the histograms of the RMS error for 500 Monte Carlo runs with and without RR measurements, compared with the method proposed in [
38].
In
Table 4, more results for the different measurement rates are reported. As before, the 90th percentile of the RMS error obtained for each run on 500 runs is reported for both cases and using the GNSS satellite or the, here proposed, autonomous solution.
These results show that, without the need for an additional payload for deep-space reception of the GPS signal, it is possible to obtain high synchronization performance for the constellation of NHRO satellites with an RMS error of the order of one nanosecond on a time span of three hours.
6. Conclusions
The results reported in the previous section open the way for various considerations.
First of all, it was demonstrated that the proposed method could be successfully implemented and used in NHRO orbits, leading to precise orbit determination and timing with performances comparable to or better than those obtained by exploiting Earth-based measurement stations. In particular, concerning the case without the use of the range-rate measurements, the results for the autonomous orbit determination are similar to the Earth-based ones for a measurement time span of 9.89 days (1.5 T) and obtain better than the Earth-based if the time span is increased, bounding the 90th percentile of the DRMS errors between 0.32 m and 0.62 m.
On the other hand, in the case of the range and range rate measurements, the results for the autonomous orbit determination are similar to the Earth-based ones starting from a measurement time span of 9.89 days (1.5 T). They remain comparable also for longer time spans, bounding the 90th percentile of the DRMS errors between 0.21 m and 0.43 m.
Last but not least, in many cases, the 90th percentile of the DRMS error of the proposed method is lower than one meter, obtaining ultra-precise orbit parameters for the spacecraft. This consideration is also true in the case of the very simple configuration that does not use the range-rate measurements, allowing for a simplified HW in the anchors. In addition, in this last case, for Time Spans longer than 1.5 T, the performances reach the meter or sub-meter levels.
Concerning the synchronization performance, the proposed method always gives better results than the one proposed in [
38] by the same authors. One of the reasons is that, in the proposed method, the timing accuracy does not depend on the spacecraft position accuracy and on the GNSS satellite’s position accuracy, moreover, the propagation errors are lower (the range measurement from the GNSS satellite used for time synchronization in [
18] depends, for example, also on the ionosphere delay).
Coming to the application, the Lunar orbiting constellation in NHRO, it is possible to assume that the constellation can be used for navigation purposes due to its high OD&T performance also with no or minimal contact with the Earth, assuring an autonomous navigation service to other spacecraft and rovers. This is an important result because, as mentioned before, in the future, the number of missions to be controlled from the Earth will grow, leading to a possible capacity limit of the terrestrial Deep Space Network. Having the possibility to provide an autonomous and independent navigation service on the Moon by exploiting a small constellation of spacecraft will overcome this limit.
Moreover, these autonomous properties can be much more useful for future missions to other planets, such as those for Mars exploration and colonization, just because, as the spacecraft gets further away, it becomes more difficult to perform the OD&T from Earth.
Last but not least is some consideration about the limit of the proposed approach. It must be recalled that at least two important hypotheses carried out in this work may be considered as potential drawbacks.
First, the main spacecraft must have the ability to solve the orbit determination problem, accumulating measurements for a long period (comparable with the orbiting period on the Halo orbit) and with a high rate of measurements (on the order of hours or fraction of hour depending on the selected configuration). In the author’s opinion, for the proposed test case and many other future missions, this is not a real limit just because it is related to the HW and SW capability (i.e., storage memory and processing power on board that are always improving with time following, for example, the Moore law [
39]). Moreover, in the specific case of the here proposed application, it is based on the Gateway, which will be able to carry all the processing and storing facilities needed because it will be a space-station designed for humans and space laboratories.
The second possible limit is due to the gravitational model used. The Circular Restricted Three-Body used here does not take into account perturbations such as, for example, the perturbation of gravitational fields for the two celestial bodies, or the effects of other deep space celestial bodies, or solar radiation pressure. These perturbations were not added to keep the method formulation as general as possible and applicable to any possible system of three bodies, postponing the introduction of the perturbations for more precise evaluation for specific cases and for the different possible solutions or applications. However, the used model is the same for the two scenarios compared in the paper, leading to the same limitation in both the methods and, in the author’s opinion, it is enough for a first comparison of the performances. Moreover, if the specific applications or spacecraft limitations will not allow the introduction of a more complex model that also include the most important perturbations, another approach can also be attempted: a time-by-time Earth-based calibration procedure, by exploiting the Earth-based Deep Space Network that, in any case, must be used for periodical Command and Control purposes.