1. INTRODUCTION
The Korea Pathfinder Lunar Orbiter (KPLO) project has a robotic lunar orbiting mission to collect high-resolution images of the lunar surface, and various scientific data. This data collection requires a stable satellite communication between the lunar orbiter and ground stations, and accurate tracking information. Accurate tracking information enables precise orbit determination (OD) of a spacecraft not only in the phasing loop transfer (maneuver burn) and lunar orbit injection phase, but also in the nominal mission operation phase (Bae et al. 2017). This OD accuracy affects attitude control accuracy, and accurate attitude control allows the payloads to capture data of the desired locations (Song et al. 2014). Except for a few contingency cases of the KPLO mission, ground stations exchange signals with the spacecraft and keep tracking the location of the spacecraft via the Doppler measurement and several ranging techniques similar to the approach used in many other lunar missions such as in the Chang’e 3 mission (Wenlin et al. 2017) and the Lunar Reconnaissance Orbiter (LRO) mission (Mazarico et al. 2012). In this paper, we focused on the ranging system and found the best ranging technique that can serve the KPLO mission successfully.
A ranging system utilizes radio frequency (RF) waveform propagation characteristics to determine the location of a spacecraft. The ranging signal being received at the ground stations, particularly the phase of the ranging signal, enables the calculation of the propagation time and the corresponding distance between the ground stations and the spacecraft. The performance of the range measurement varies depending on the structure of the ranging signal and the level of thermal noise. There are three commonly used criteria for range measurement performance. First is the measurement error of the time delay, which determines how accurate the ranging measurement is. Second is the probability of the ranging signal acquisition, while the third criterion is the measurement time that is required to measure the ranging signal. A longer measurement time renders higher accuracy and acquisition probability, but it also yields less location measurements in a given time. Thus, a system needs to be designed to meet the accuracy and acquisition probability requirements while maintaining a short measurement time.
Assuming that the thermal noise level is the same, the accuracy of a ranging system differs depending on the ranging technique used in the system. Radiometric ranging techniques such as tone ranging and pseudo-noise (PN) ranging have been the most widely used methods for tracking spacecrafts. Recently, satellite laser ranging (SLR) gained attention owing to its high precision (Lim et al. 2011; Choi et al. 2014). The Science and Technology Satellite (STSAT)-2C OD was also successfully launched and enhanced orbit predictions can be generated owing to its cm-level range precision (Kim et al. 2015). Although SLR is a subject of active research, owing to the technical assistance agreement (TAA) between Korea Aerospace Research Institute (KARI) and National Aeronautics and Space Administration (NASA), radiometric ranging techniques, which are supported by the Deep Space Network (DSN) and Near Earth Network (NEN), will be used for the KPLO mission. Since the radiometric ranging measurement accuracy for the DSN and NEN can be 1–2 m (JPL 2000), in this paper, it is demonstrated that the radiometric ranging is sufficiently accurate to satisfy the KPLO accuracy requirement.
Sequential ranging and PN ranging with T2B and T4B (weighted-voting balanced Tausworthe, voting v = 4) code, which are radiometric ranging techniques, can be supported as ranging methods. Thus, it is necessary to understand each ranging technique’s characteristics and find the best ranging technique for the KPLO mission. There have been several attempts to compare the performance of sequential ranging and PN ranging systems (Berner & Bryant 2002; Berner et al. 2007). It is known that PN ranging has the advantage of modifying its integration time in real time, while sequential ranging integration time is fixed by the code definition (Berner & Bryant 2002). This gives PN ranging more freedom to adjust its measurement performance in real time. A comparison between various PN codes was carried out in the work of Boscagli et al. 2007. The authors demonstrated that T4B is a valid choice for their mission considering the jitter requirement. However, an appropriate choice of a ranging method for the KPLO mission will be different since the requirements and specifications are different.
Aside from the ranging method being used in the system, the spacecraft communication system and the influence of ground stations on the measurement accuracy is also important. This naturally raises the problem of finding the end-to-end ranging accuracy. Most studies in the literature focused on the accuracy error due to thermal noise and the ranging methodology, while other error factors were often neglected (Ren et al. 2011; Jensen et al. 2013; Sheng 2014). In the CCSDS (2014), the ranging accuracy performance was addressed in terms of the received signal-to-noise ratio (SNR) regime, and the effect of the ranging clock frequency was evaluated. However, when other error sources exist, the decrease in the accuracy is not inversely proportional to the square root of the SNR. For instance, if the dominant factor of the end-to-end ranging accuracy is on the transponder, increasing the clock frequency or SNR will not significantly increase the accuracy.
In this paper, we first compared the accuracy of sequential ranging and PN ranging due to thermal noise. The endto-end ranging accuracy includes the ground sub-system instrumentation range error and the spacecraft’s communication RF circuitry error. The end-to-end ranging system accuracy analysis can help in the design of an overall efficient system specification. We further investigated the dominant factors of the ranging accuracy error and determined appropriate ranging parameters for the KPLO mission. Finally, we determined the optimal ranging method in terms of the measurement time that meets the requirements of accuracy and acquisition probability. Finding the ranging method with the shortest measurement time allows as many range measurements as possible, and this could help in conducting a successful mission.
The rest of this paper is organized as follows. The system model is described in Section 2, where we also analyzed the accuracy of the PN and sequential ranging methods. In Section 3, the end-to-end ranging accuracy and the optimal ranging methods suitable for the KPLO mission are presented. Finally, the conclusions are given in Section 4.
2. SYSTEM AND SIGNAL MODEL
2.1 System Model
The pulse-coded modulation/phase-shift keyed/phasemodulated (PCM/PSK/PM) modulation has been widely used in space communication systems and this modulation will also be utilized for the KPLO mission. This modulation process enables simultaneous transmission of telecommand/ telemetry and ranging signal. For instance, in uplink transmission, the non-return to zero (NRZ) telecommand data is PSK modulated, which is then PM modulated along with the ranging signal. Depending on the bandwidth constraint, a squarewave or sinewave subcarrier can be employed. Details of the characteristics of the modulation schemes are well addressed in the DSN handbook (JPL 2000). A sienwave subcarrier requires spectrum bandwidth than a squarewave, which is suitable for the KPLO mission since it will be conducted in a bandwidth limited condition. Therefore, in this study, we considered a sinewave modulated signal. The uplink PCM/PSK/PM modulated signal, x(t), can be expressed as follows
\(x(t)=\sqrt{2} \sin \left(2 \pi f_{c} t-\mathrm{m}_{\mathrm{TC}} d_{\mathrm{TC}}(t) \sin 2 \pi f_{T C} t+\mathrm{m}_{\mathrm{R}} R(t) \sin 2 \pi f_{R} t\right)\) (1)
where fc is the carrier frequency. The second phase term is the telemetry signal, where mTC is the modulation index, dd( t ) is the data, and fTC is the telecommand frequency. The third term represents the ranging signal, where mR is the modulation index, and f R is the range signal frequency. The ranging signal R ( t ) is configured in the ground station, and differs depending on the ranging method adopted in the system. For instance, R ( t ) is a constant one when sequential ranging is used. When PN ranging is used, R(t) is ±1, where its period is - fR /2. We have are assumed that the ground stations that support the KPLO mission are capable of configuring sequential and PN ranging. This can possibly enhance the ranging performance by allowing flexibility in selecting a ranging method interchangeably.
In terms of the ranging method, the KPLO communication system will use transparent ranging. Transparent ranging does not require decoding and re-encoding of the ranging signal, which renders a simplified transponder design at the expense of uplink noise accumulation. In contrast, regenerative ranging eliminates the uplink noise term; thus, regenerative PN ranging is widely used in SNR limited Category B missions beyond Mars. Since the KPLO communication system will be operating in the Category A mission, transparent ranging is appropriate, and the following analysis is based on transparent ranging.
2.2 PN Ranging Signal Model
NASA’s DSN ranging system supports PN ranging in addition to the conventional sequential ranging, and it has been widely used in recent space missions. The PN code is a composite code generated from a logical operation of periodic binary component codes. Details of the code construction are well described in the CCSDS (2014). In this section, the ranging error due to thermal noise is addressed, while the composite end-to-end ranging error will be addressed in Section 3.1.
For a clock frequency f , the standard deviation of the ranging measurement error in rms meters, - σ T , is defined as follows (JPL 2000)
\(\sigma_{T}=\frac{c}{f \cdot A_{C} \cdot R_{1} \sqrt{32 \pi^{2} \cdot T_{P N} \cdot P_{R} / N_{0}}}\) (2)
where the parameters represent the following
c : velocity of light
Ac : fractional loss of correlation amplitude due to frequency mismatch
R1 : cross-correlation factor of the first component code
TPN : range measurement duration
P R / N 0 downlink SNR.
The cross-correlation factor of the first component code influences the accuracy, while the rest of the component codes ( R n , where n≠1) determine the acquisition probability, where they are defined by the code construction (Massey at el. 2007). The range measurement duration, T PN , affects both the accuracy and acquisition probability. Since the duration, T PN , can be adjusted in real-time during operation, it should be carefully analyzed to carry out appropriate adjustment during the mission. It should also be noted that PN ranging enables the measurement of all ranging components at once; thus, the integration time is the measurement time. Assuming that the measurement duration that satisfies the accuracy requirement is Tacc , and the duration that satisfies the acquisition probability requirement is Tacq , the measurement duration that satisfies both conditions should be the maximum of the two measure-ment time durations, as given below (JPL 2000)
\(T_{P N}=\max \left(T_{a c q}, T_{a c c}\right)\) (3)
Since the measurement duration takes the maximum of the required duration time, there exists a margin on one measurement criterion, and vice versa. For instance, if Tacq is greater than Tacc , there exists an accuracy margin corresponding to 10 log(T PN - T acc ) dB. Similarly, if T acc is greater than T acq , the system achieves higher acquisition probability corresponding to 10 log( T PN - T acq ) dB.
2.3 Sequential Ranging Signal Model
Sequential ranging signal structure is composed of a sequence of tone components, in which there are one major tone and n-1 number of minor tones. The major tone has the highest frequency and it determines the accuracy of ranging. The minor tones, which are often called ambiguity resolving tones, are followed by the major tone. They determine the acquisition probability and ranging ambiguity. Given a clock frequency f , the measurement error in rms meters, σ T , can be expressed as follows (JPL 2000)
\(\sigma_{T}=\frac{c}{f \cdot A_{C} \sqrt{32 \pi^{2} \cdot T_{1} \cdot P_{R} / N_{0}}}\) (4)
where T1 is the integration time of the major tone. We have assumed that all the parameters are defined the same as in PN ranging, except the integration time. The range measurement time of sequential ranging includes not only the integration time of the major tone, T1 , but also the integration time of minor tones denoted as T2 . Assuming that the total number of tones is n, the measurement time of sequential ranging Tseq is defined as follows (JPL 2000)
\(T_{s e q}=\left(2+T_{1}\right)+(n-1)\left(1+T_{2}\right)+1\) (5)
where n+1 number of one second dead times are located at the start of the tone, in between each tone, and at the end of the tone. Dead times are inserted as a margin for error, but they are not used in the signal integration. This is considered as the drawback of a sequential ranging system.
3.
PERFORMANCE ANALYSIS AND COMPARISON
3.1 End-to-End Ranging Accuracy Analysis -
The system performance is limited by systematic errors and random errors. The former could be eliminated through calibration or compensated by adding offsets in the system. The latter, on the other hand, are usually unavoidable and influence the reliability of the measurement. We have assumed that systematic errors are predictable and controlled through calibration, while random errors are uncontrollable and influence the performance of ranging accuracy. In the end-to-end ranging system, the major random error sources are the ground stations, spacecraft communication system, and thermal noise. We have assumed that the random variables that can capture each of the random delay error sources are D G , D S , and D T , and their standard deviations are σ G , σ S , and σ T , respectively.
The ground systems that support the KPLO mission will be NASA’s DSN and KDSA (Korea Deep Space Antenna). Since KDSA is still under development, we will consider the specifications of NASA’s DSN as a preliminary study. The ground system error sources include oscillator noise, group delay variation, and oscillator calibration. Details of the error sources and specifications are available in the DSN Handbook (JPL 2000).
The simplified spacecraft communication system archi-tecture is illustrated in Fig. 1 . The transponder is composed of the transmitter and the receiver, where modulation and demodulation are executed, respectively. In the radio fre-quency distribution unit (RFDU) block, there are diplexers that allow the transmitted and received signals to be conveyed on the same antenna. Furthermore, there are spurious rejection filters and switches that enable routing of the primary and secondary signal. The random delay error in the communication system is due to the signal routed from the RFDU and transponder backward and forward.
Fig. 1. Spacecraft communication system block diagram relevant for ranging performance.
We have assumed that the total delay is represented by a random variable D, which is the summation of each delay error source as follows
\(D=D_{G}+D_{S}+D_{T}\) (6)
Since the ground systems and the spacecraft communication system are designed independently, it is safe to assume that the random variables D G , D S , and D T are mutually independent. Thus, the total delay error variance is the summation of each variance, and its standard variation can be expressed as follows
\(\sigma=\sqrt{\sigma_{G}^{2}+\sigma_{S}^{2}+\sigma_{T}^{2}}\) (7)
In the remainder of this section, analytical results of the accuracy are presented in terms of the received SNR and the measurement time. In the analysis, we have assumed that the fractional loss - A c is 1, and we used the T2B code structure as PN ranging and compared its accuracy performance with sequential ranging. The outperformance of T2B over T4B in the KPLO mission is presented in Section 3.2. We used T to simplify the notation of integration time for PN ( - T PN ) and sequential ranging ( T 1 ).
Fig. 2 shows the analytical results of the ranging measure-ment error due to thermal noise ( σT ). The figure shows that the accuracy increases as the clock frequency increases. For instance, if the clock frequency is doubled, the accuracy is reduced by half. The figure also shows that the sequential ranging accuracy is slightly higher than that of PN ranging. This is an expected result since the PN ranging peak power is spread out beside the clock frequency; thus, PN ranging has lower peak power than sequential ranging. However, the difference decreases as the SNR increases, and since KPLO will mostly operate in a high SNR regime, it is safe to assume that the difference is negligible. Note that the integration time of PN ranging is the measurement time, and PN ranging requires less measurement time for each range measurement than sequential ranging. Thus, PN ranging has more margin to increase its measurement time and accuracy than sequential ranging.
Fig. 2. Measurement error profile of sequential and PN ranging due to thermal noise.
Fig. 3 shows the end-to-end range measurement error. The figure shows that in the KPLO communication system, the estimated standard deviation of the time delay is approxi-mately 64 ns. The corresponding one-way distance error variation, σ S , is 9.6 m. In terms of the ground station error, we used the specification defined in module 203 of the DSN Handbook (JPL 2000), where the distance error variation, - σ G , is 1.89 m. The figure shows the same overall tendency depicted in Fig. 2 . However, the major difference is that the accuracy does not increase as much as in Fig. 2 in the high SNR regime. This is because the dominant error source in the high SNR regime is from the spacecraft. Thus, the accuracy improvement by SNR is insignificant in the high SNR regime. Similarly, accuracy improvement by clock frequency is also negligible since the accuracy error due to thermal noise is much smaller than the error from the spacecraft. This finding can be useful in selecting the ranging clock frequency for the KPLO mission. A typical value for the ranging clock frequency is 1 MHz; however, the KPLO transponder under development only operates on 500 kHz ranging clock frequency and it requires additional efforts and resources to have the capability to handle 1 MHz ranging tone. These results show that 1 MHz clock frequency does not render accuracy improvement considering the current end-to-end system specifications.
Fig. 3. End-to-end measurement error profile of sequential and PN ranging incorporating the KPLO communication system error of 9.6 m and the ground system error of 1.89 m.
In Fig. 4 , the ranging measurement error due to thermal noise and the end-to-end error are compared for clock frequencies of 500 kHz and 1 MHz. The figure demonstrates that increasing the clock frequency from 500 kHz to 1 MHz results in up to 1 m accuracy improvement for the KPLO communication operation region, which is approximately 25 to 30 dB SNR. However, since the end-to-end error incorporates the KPLO communication system error, which is the major error source, the benefit from increasing the clock frequency is meaningless since the error due to thermal noise is much smaller than the one caused by the KPLO communication system (σ T ≪ σ S )
fig. 4. Comparison of measurement error due to thermal noise and endto- end error.
3.2 Comparisons on Ranging Accuracy
This section presents a comparison of the performance of the ranging techniques in terms of the measurement time.
The ranging techniques considered were sequential ranging, T2B code, and T4B code ranging, which are all CCSDS defined formats and are available for the KPLO mission. Table 1 presents the ranging requirements and specifications, where the values have been practically selected, but could differ from the actual values used in the KPLO mission. The accuracy requirement of 1 m (rms) is for the error due to thermal noise only, and we omitted other error sources from the specification since thermal noise is the only error source that is manageable by adjusting parameters such as the transmit power and the measurement time.
Table 1. Ranging requirements and specifications
Fig. 5 plots the duration times of each ranging measurement in terms of the received SNR. The figure shows that in the low SNR regime where the SNR is below 17.5 dB, sequential ranging has the shortest measurement time, while in the high SNR regime where the SNR is above 17.5 dB, the T2B code has the shortest measurement time. Although the T4B code out performed sequential ranging above SNR of 22 dB, its measurement time is greater than that of T2B in the SNR span considered in this analysis. The figure also demonstrates that selecting the best ranging technique can help in saving power. For instance, sequential ranging required 25 dB received SNR and 25 sec measurement time to achieve an accuracy of 1 m, but T2B requires only approximately 20 dB received SNR with the same measurement time of 25 sec, which renders approximately 5 dB gain over sequential ranging. This received SNR gain allows the possibility to increase the data rate or reduce the error rate of the telecommand, and this can be achieved by increasing the telecommand modulation index mTC in Eq. (1).
Fig. 5. Comparison of measurement time duration of each ranging under the constraint of accuracy requirement of 1 m (rms) and clock frequency of 500 kHz.
Fig. 6 shows a detailed analysis of the measurement time. In general, the T4B code provides a better accuracy than the T2B code, whereas the T2B code has better acquisition probability than the T4B code. It can be observed in Fig. 6 that Tacc of T4B is shorter than that of T2B, and Tacq of T2B is shorter than that of T4B. The figure also reveals that when the accuracy requirement is relaxed, such that Tacc of T2B is shorter than Tacq of T4B, T2B performs better than T4B.
Fig. 6. T2B/T4B PN ranging measurement duration required to satisfy the requirement of 1 m (rms) accuracy with 500 kHz clock frequency
Fig. 7 shows another comparison of results under the same constraints as in Fig. 6. The figure shows the ratio of the measurement time of sequential ranging over PN ranging. This result also demonstrates that sequential ranging had the best performance in the low SNR regime while T2B had the best performance in the high SNR regime. A notable analytical result is that the ratio of the measurement time exponentially increases as the received SNR increases.
Fig. 7. Comparison of measurement duration of sequential ranging and T2B/T4B PN ranging that satisfies the requirement of 1 m (rms) accuracy with 500 KHz clock frequency.
4. CONCLUSIONS
In this paper, the analytical results of the end-to-end ranging measurement performance for the KPLO mission were presented. Error sources such as in the ground stations and spacecraft communication system are often neglected when it comes to selecting a ranging technique and finding the corresponding ranging parameters. However, when these error sources are incorporated, results presented in this paper demonstrate that the selection of a ranging technique and its clock frequency (T2B with 500 kHz clock frequency) gives a different performance from the ones that solely considered the thermal noise error (T4B with 1 MHz clock frequency). We also showed that this is because the dominant error source is independent of the error due to thermal noise.
We further analyzed and compared ranging techniques using the criterion of measurement duration, while the accuracy requirements and specifications were met. We found that the ranging technique with the shortest measurement duration differ in the received SNR, i.e., the best ranging technique differ in the received SNR. Given the KPLO requirements, it was shown that the performance threshold between sequential and T2B PN ranging technique was a received SNR of 17.5 dB. Considering the KPLO operating received SNR region of 25–30 dB, we found that the T2B code performed better than the T4B code. We further demonstrated that by selecting the best ranging technique, there was a power gain of 5 dB for the received SNR of 25 dB. This is because the T4B code does not take advantage of the short measurement duration required to achieve accuracy in the KPLO mission. These findings reveal that PN ranging with the T2B code is the best ranging technique among the ranging techniques available for the KPLO mission.
ACKNOWLEDGMENTS
This work was supported by the Korea Aerospace Research Institute (KARI) under a contract with the Ministry of Science and ICT (MSIT) through the “Development of Pathfinder Lunar Orbiter and Key Technologies” project (No. SR18040).
References
- Bae J, Song YJ, Kim YK, Kim B, Burn delay analysis of the lunar orbit insertion for Korea Pathfinder Lunar Orbiter, J. Astron. Space Sci. 34, 281-287 (2017). https://doi.org/10.5140/JASS.2017.34.4.281
- Berner JB, Bryant SH, Operations comparison of deep space ranging types: Sequential tone vs. pseudo-noise, Proceedings of the 2002 IEEE Aerospace Conference, Big Sky, MT, 9-16 Mar 2002. https://doi.org/10.1109/AERO.2002.1035264
- Berner JB, Bryant SH, Kinman PW, Range measurement as practiced in the deep space network, Proc. IEEE 95, 2202-2214 (2007). https://doi.org/10.1109/JPROC.2007.905128
- Boscagli G, Holsters P, Vassallo E, Visintin M, PN regenerative ranging and its compatibility with telecommand and telemetry signals, Proc. IEEE 95, 2224-2234 (2007). https://doi.org/10.1109/JPROC.2007.905186
- CCSDS (The Consultative Committee for Space Data Systems), Pseudo-Noise (PN) ranging systems, CCSDS 414.1-B-2, Blue Book (2014). https://public.ccsds.org/Pubs/414x1b2.pdf
- Choi MS, Lim HC, Choi EJ, Park E, Yu SY, et al., Performance analysis of the first Korean satellite laser ranging system, J. Astron. Space Sci. 31, 225-233 (2014). https://doi.org/10.5140/JASS.2014.31.3.225
- Jensen JR, Haskins CB, DeBoy CC, Regenerative PN ranging experience with New Horizons during 2012, Proceedings of the 2013 IEEE Aerospace Conference, Big Sky, MT, 2-9 Mar 2013. https://doi.org/10.1109/AERO.2013.6496983
- JPL (Jet Propulsion Laboratory), DSN Telecommunication Link Design Handbook, 810-005, Rev. E (2000).
- Kim YR, Park E, Kucharski D, Lim HC, Orbit determination using SLR data for STSAT-2C: Short-arc analysis, J. Astron. Space Sci. 32, 189-200 (2015). https://doi.org/10.5140/JASS.2015.32.3.189
- Lim HC, Bang SC, Yu SY, Seo YK, Park ES, et al., Study on the optoelectronic design for Korean mobile satellite laser ranging system, J. Astron. Space Sci. 28, 155-162 (2011). https://doi.org/10.5140/JASS.2011.28.2.155
- Massey JL, Boscagli G, Vassallo E, Regenerative pseudo-noise (PN) ranging sequences for deep-space missions. Int. J. Satell. Commun. Netw. 25, 285-304 (2007). https://doi.org/10.1002/sat.877
- Mazarico E, Rowlands DD, Neumann GA, Smith DE, Torrence MH, et al., Orbit determination of the Lunar Reconnaissance Orbiter, J. Geod. 86, 193-207 (2012). https://doi.org/10.1007/s00190-011-0509-4
- Ren G, Sun C, Ni H, Bai Y, OFDM-based precise ranging technique in space applications, IEEE Trans. Aerosp. Electron. Syst. 47, 2217-2221 (2011). https://doi.org/10.1109/TAES.2011.5937294
- Sheng B, Enhanced OFDM-based ranging method for space applications, IEEE Trans. Aerosp. Electron. Syst. 50, 1606-1609 (2014). https://doi.org/10.1109/TAES.2014.120583
- Song YJ, Ahn SI, Sim ES, Development strategy of orbit determination system for Korea's lunar mission: Lessons from ESA, JAXA, ISRO and CNSA's experiences, J. Astron. Space Sci. 31, 247-264 (2014). https://doi.org/10.5140/JASS.2014.31.3.247
- Wenlin T, Peng X, Songjie H, Jianfeng C, Peng D, et al., Chang'e 3 lunar mission and upper limit on stochastic background of gravitational wave around the 0.01 Hz band. Adv. Space Res. 60, 1315-1326 (2017). https://doi.org/10.1016/j.asr.2017.06.008