The ground tracking support is a critical factor for the navigation performance of spacecraft orbiting around the Moon. Because of the tracking limit of antennas, only a small number of facilities can support lunar missions. Therefore, case studies for various ground tracking support conditions are needed for lunar missions on the stage of preliminary mission analysis. This study analyzes the ground supporting condition effect on orbit determination (OD) of Korea Pathfinder Lunar Orbiter (KPLO) in the lunar orbit. For the assumption of ground support conditions, daily tracking frequency, cut-off angle for low elevation, tracking measurement accuracy, and tracking failure situations were considered. Two antennas of deep space network (DSN) and Korea Deep Space Antenna (KDSA) are utilized for various tracking conditions configuration. For the investigation of the daily tracking frequency effect, three cases (full support, DSN 4 pass/day and KDSA 4 pass/day, and DSN 2 pass/day and KDSA 2 pass/day) are prepared. For the elevation cut-off angle effect, two situations, which are 5 deg and 10 deg, are assumed. Three cases (0%, 30%, and 50% of degradation) were considered for the tracking measurement accuracy effect. Three cases such as no missing, 1-day KDSA missing, and 2-day KDSA missing are assumed for tracking failure effect. For OD, a sequential estimation algorithm was used, and for the OD performance evaluation, position uncertainty, position differences between true and estimated orbits, and orbit overlap precision according to various ground supporting conditions were investigated. Orbit prediction accuracy variations due to ground tracking conditions were also demonstrated. This study provides a guideline for selecting ground tracking support levels and preparing a backup plan for the KPLO lunar mission phase.

The Korean Institute of Technology Satellite (KITSAT-1) is the first satellite developed by the Satellite Technology Research Center and the University of Surrey. KITSAT-1 is orbiting the Earth’s orbit as space debris with a 1,320 km altitude after the planned mission. Due to its relatively small size and altitude, tracking the KITSAT-1 was a difficult task. In this research, we analyzed the tracking results of KITSAT-1 for one year using the Midland Space Radar (MSR) in Texas and the Poker Flat Incoherent Scatter Radar (PFISR) in Alaska operated by LeoLabs, Inc. The tracking results were analyzed on a weekly basis for MSR and PFISR. The observation was conducted by using both stations at an average frequency of 10 times per week. The overall corrected range measurements for MSR and PFISR by LeoLabs were under 50 m and 25 m, respectively. The ionospheric delay, the dominant error source, was confirmed with the International Reference of Ionosphere-16 model and Global Navigation Satellite System data. The weekly basis orbit determination results were compared with two-line element data. The comparison results were used to confirm the orbital consistency of the estimated orbits.

In this study, the observational arc-length effect on orbit determination (OD) for the Korea Pathfinder Lunar Orbiter (KPLO) in the Earth-Moon Transfer phase was investigated. For the OD, we employed a sequential estimation using the extended Kalman filter and a fixed-point smoother. The mission periods, comprised between the perigee maneuvers (PM) and the lunar orbit insertion (LOI) maneuver in a 3.5 phasing loop of the KPLO, was the primary target. The total period was divided into three phases: launch–PM1, PM1–PM3, and PM3–LOI. The Doppler and range data obtained from three tracking stations [included in the deep space network (DSN) and Korea Deep Space Antenna (KDSA)] were utilized for the OD. Six arc-length cases (24 hrs, 48 hrs, 60 hrs, 3 days, 4 days, and 5 days) were considered for the arc-length effect investigation. In order to evaluate the OD accuracy, we analyzed the position uncertainties, the precision of orbit overlaps, and the position differences between true and estimated trajectories. The maximum performance of 3-day OD approach was observed in the case of stable flight dynamics operations and robust navigation capability. This study provides a guideline for the flight dynamics operations of the KPLO in the trans-lunar phase.

In this study, the precise orbit determination (POD) software is developed for optical observation. To improve the performance of the estimation algorithm, a nonlinear batch filter, based on the unscented transform (UT) that overcomes the disadvantages of the least-squares (LS) batch filter, is utilized. The LS and UT batch filter algorithms are verified through numerical simulation analysis using artificial optical measurements. We use the real optical observation data of a low Earth orbit (LEO) satellite, Cryosat-2, observed from optical wide-field patrol network (OWL-Net), to verify the performance of the POD software developed. The effects of light travel time, annual aberration, and diurnal aberration are considered as error models to correct OWL-Net data. As a result of POD, measurement residual and estimated state vector of the LS batch filter converge to the local minimum when the initial orbit error is large or the initial covariance matrix is smaller than the initial error level. However, UT batch filter converges to the global minimum, irrespective of the initial orbit error and the initial covariance matrix.

The Optical Wide-field patroL-Network (OWL-Net) is a Korean optical space surveillance system used to track and monitor objects in space. In this study, the characteristics of four Initial Orbit Determination (IOD) methods were analyzed using artificial observational data from Low Earth Orbit satellites, and an appropriate IOD method was selected for use as the initial value of Precise Orbit Determination using OWL-Net data. Various simulations were performed according to the properties of observational data, such as noise level and observational time interval, to confirm the characteristics of the IOD methods. The IOD results produced via the OWL-Net observational data were then compared with Two Line Elements data to verify the accuracy of each IOD method. This paper, thus, suggests the best method for IOD, according to the properties of angles-only data, for use even when the ephemeris of a satellite is unknown.

In this study, orbit determination (OD) simulation for the Korea Pathfinder Lunar Orbiter (KPLO) was accomplished for investigation of the observational arc-length effect using a sequential estimation algorithm. A lunar polar orbit located at 100 km altitude and 90° inclination was mainly considered for the KPLO mission operation phase. For measurement simulation and OD for KPLO, the Analytical Graphics Inc. Systems Tool Kit 11 and Orbit Determination Tool Kit 6 software were utilized. Three deep-space ground stations, including two deep space network (DSN) antennas and the Korea Deep Space Antenna, were configured for the OD simulation. To investigate the arc-length effect on OD, 60-hr, 48-hr, 24-hr, and 12-hr tracking data were prepared. Position uncertainty by error covariance and orbit overlap precision were used for OD performance evaluation. Additionally, orbit prediction (OP) accuracy was also assessed by the position difference between the estimated and true orbits. Finally, we concluded that the 48-hr-based OD strategy is suitable for effective flight dynamics operation of KPLO. This work suggests a useful guideline for the OD strategy of KPLO mission planning and operation during the nominal lunar orbits phase.

This study presents the application of satellite laser ranging (SLR) to orbit determination (OD) of high-Earth-orbit (HEO) satellites. Two HEO satellites are considered: the Quasi-Zenith Satellite-1 (QZS-1), a Japanese elliptical-inclinedgeosynchronous- orbit (EIGSO) satellite, and the Compass-G1, a Chinese geostationary-orbit (GEO) satellite. One week of normal point (NP) data were collected for each satellite to perform the OD based on the batch least-square process. Five SLR tracking stations successfully obtained 374 NPs for QZS-1 in eight days, whereas only two ground tracking stations could track Compass-G1, yielding 68 NPs in ten days. Two types of station bias estimation and a station data weighting strategy were utilized for the OD of QZS-1. The post-fit root-mean-square (RMS) residuals of the two week-long arcs were 11.98 cm and 10.77 cm when estimating the biases once in an arc (MBIAS). These residuals were decreased significantly to 2.40 cm and 3.60 cm by estimating the biases every pass (PBIAS). Then, the resultant OD precision was evaluated by the orbit overlap method, yielding three-dimensional errors of 55.013 m with MBIAS and 1.962 m with PBIAS for the overlap period of six days. For the OD of Compass-G1, no station weighting strategy was applied, and only MBIAS was utilized due to the lack of NPs. The post-fit RMS residuals of OD were 8.81 cm and 12.00 cm with 49 NPs and 47 NPs, respectively, and the corresponding threedimensional orbit overlap error for four days was 160.564 m. These results indicate that the amount of SLR tracking data is critical for obtaining precise OD of HEO satellites using SLR because additional parameters, such as station bias, are available for estimation with sufficient tracking data. Furthermore, the stand-alone SLR-based orbit solution is consistently attainable for HEO satellites if a target satellite is continuously trackable for a specific period.

The deep space orbit determination software (DSODS) is a part of a flight dynamic subsystem (FDS) for the Korean Pathfinder Lunar Orbiter (KPLO), a lunar exploration mission expected to launch after 2018. The DSODS consists of several sub modules, of which the orbit determination (OD) module employs a weighted least squares algorithm for estimating the parameters related to the motion and the tracking system of the spacecraft, and subroutines for performance improvement and detailed analysis of the orbit solution. In this research, DSODS is demonstrated and validated at lunar orbit at an altitude of 100 km using actual Lunar Prospector tracking data. A set of a priori states are generated, and the robustness of DSODS to the a priori error is confirmed by the NASA planetary data system (PDS) orbit solutions. Furthermore, the accuracy of the orbit solutions is determined by solution comparison and overlap analysis as about tens of meters. Through these analyses, the ability of the DSODS to provide proper orbit solutions for the KPLO are proved.

This paper presents an overview of deep space orbit determination software (DSODS), as well as validation and verification results on its event prediction capabilities. DSODS was developed in the MATLAB object-oriented programming environment to support the Korea Pathfinder Lunar Orbiter (KPLO) mission. DSODS has three major capabilities: celestial event prediction for spacecraft, orbit determination with deep space network (DSN) tracking data, and DSN tracking data simulation. To achieve its functionality requirements, DSODS consists of four modules: orbit propagation (OP), event prediction (EP), data simulation (DS), and orbit determination (OD) modules. This paper explains the highest-level data flows between modules in event prediction, orbit determination, and tracking data simulation processes. Furthermore, to address the event prediction capability of DSODS, this paper introduces OP and EP modules. The role of the OP module is to handle time and coordinate system conversions, to propagate spacecraft trajectories, and to handle the ephemerides of spacecraft and celestial bodies. Currently, the OP module utilizes the General Mission Analysis Tool (GMAT) as a third-party software component for highfidelity deep space propagation, as well as time and coordinate system conversions. The role of the EP module is to predict celestial events, including eclipses, and ground station visibilities, and this paper presents the functionality requirements of the EP module. The validation and verification results show that, for most cases, event prediction errors were less than 10 millisec when compared with flight proven mission analysis tools such as GMAT and Systems Tool Kit (STK). Thus, we conclude that DSODS is capable of predicting events for the KPLO in real mission applications.

The optical wide-field patrol network (OWL-Net) is a Korean optical surveillance system that tracks and monitors domestic satellites. In this study, a batch least squares algorithm was developed for optical measurements and verified by Monte Carlo simulation and covariance analysis. Potential error sources of OWL-Net, such as noise, bias, and clock errors, were analyzed. There is a linear relation between the estimation accuracy and the noise level, and the accuracy significantly depends on the declination bias. In addition, the time-tagging error significantly degrades the observation accuracy, while the time-synchronization offset corresponds to the orbital motion. The Cartesian state vector and measurement bias were determined using the OWL-Net tracking data of the KOMPSAT-1 and Cryosat-2 satellites. The comparison with known orbital information based on two-line elements (TLE) and the consolidated prediction format (CPF) shows that the orbit determination accuracy is similar to that of TLE. Furthermore, the precision and accuracy of OWL-Net observation data were determined to be tens of arcsec and sub-degree level, respectively.

A strategy is needed for a regional survey of geosynchronous orbits (GSOs) to monitor known space objects and detect uncataloged space objects. On the basis of the Inter-Agency Debris Committee’s recommendation regarding the protected region of geosynchronous Earth orbit (GEO), target satellites with perigee and apogee of GEO ± 200 km and various inclinations are selected for analysis. The status of the GSO region was analyzed using the satellite distribution based on the orbital characteristics in publicly available two-line element data. Natural perturbation effects cause inactive satellites to drift to two stable longitudinal points. Active satellites usually maintain the designed positions as a result of regular or irregular maneuver operations against their natural drift. To analyze the detection rate of a single optical telescope, 152 out of 412 active satellites and 135 out of 288 inactive satellites in the GSO region were selected on the basis of their visibility at the observation site in Daejeon, Korea. By using various vertical view ranges and various numbers of observations of the GSO region, the detection efficiencies were analyzed for a single night, and the numbers of follow-up observations were determined. The orbital estimation accuracies were also checked using the arc length and number of observed data points to maintain the GSO satellite catalog.

We estimated the orbit of the Communication, Ocean and Meteorological Satellite (COMS), a Geostationary Earth Orbit (GEO) satellite, through data from actual optical observations using telescopes at the Sobaeksan Optical Astronomy Observatory (SOAO) of the Korea Astronomy and Space Science Institute (KASI), Optical Wide field Patrol (OWL) at KASI, and the Chungbuk National University Observatory (CNUO) from August 1, 2014, to January 13, 2015. The astrometric data of the satellite were extracted from the World Coordinate System (WCS) in the obtained images, and geometrically distorted errors were corrected. To handle the optically observed data, corrections were made for the observation time, light-travel time delay, shutter speed delay, and aberration. For final product, the sequential filter within the Orbit Determination Tool Kit (ODTK) was used for orbit estimation based on the results of optical observation. In addition, a comparative analysis was conducted between the precise orbit from the ephemeris of the COMS maintained by the satellite operator and the results of orbit estimation using optical observation. The orbits estimated in simulation agree with those estimated with actual optical observation data. The error in the results using optical observation data decreased with increasing number of observatories. Our results are useful for optimizing observation data for orbit estimation.

To protect and manage the Korean space assets including satellites, it is important to have precise positions and orbit information of each space objects. While Korea currently lacks optical observatories dedicated to satellite tracking, the Korea Astronomy and Space Science Institute (KASI) is planning to establish an optical observatory for the active generation of space information. However, due to geopolitical reasons, it is difficult to acquire an adequately sufficient number of optical satellite observatories in Korea. Against this backdrop, this study examined the possible locations for such observatories, and performed simulations to determine the differences in precision of optical orbit estimation results in relation to the relative baseline distance between observatories. To simulate more realistic conditions of optical observation, white noise was introduced to generate observation data, which was then used to investigate the effects of baseline distance between optical observatories and the simulated white noise. We generated the optical observations with white noise to simulate the actual observation, estimated the orbits with several combinations of observation data from the observatories of various baseline differences, and compared the estimated orbits to check the improvement of precision. As a result, the effect of the baseline distance in combined optical GEO satellite observation is obvious but small compared to the observation resolution limit of optical GEO observation.

In this study, we present the results of orbit determination (OD) using satellite laser ranging (SLR) data for the Science and Technology Satellite (STSAT)-2C by a short-arc analysis. For SLR data processing, the NASA/GSFC GEODYN II software with one year (2013/04 – 2014/04) of normal point observations is used. As there is only an extremely small quantity of SLR observations of STSAT-2C and they are sparsely distribution, the selection of the arc length and the estimation intervals for the atmospheric drag coefficients and the empirical acceleration parameters was made on an arc-to-arc basis. For orbit quality assessment, the post-fit residuals of each short-arc and orbit overlaps of arcs are investigated. The OD results show that the weighted root mean square post-fit residuals of short-arcs are less than 1 cm, and the average 1-day orbit overlaps are superior to 50/600/900 m for the radial/cross-track/along-track components. These results demonstrate that OD for STSAT-2C was successfully achieved with cm-level range precision. However its orbit quality did not reach the same level due to the availability of few and sparse measurement conditions. From a mission analysis viewpoint, obtaining the results of OD for STSAT-2C is significant for generating enhanced orbit predictions for more frequent tracking.

In this paper, a brief but essential development strategy for the lunar orbit determination system is discussed to prepare for the future Korea’s lunar missions. Prior to the discussion of this preliminary development strategy, technical models of foreign agencies for the lunar orbit determination system, tracking networks to measure the orbit, and collaborative efforts to verify system performance are reviewed in detail with a short summary of their lunar mission history. Covered foreign agencies are European Space Agency, Japan Aerospace Exploration Agency, Indian Space Research Organization and China National Space Administration. Based on the lessons from their experiences, the preliminary development strategy for Korea’s future lunar orbit determination system is discussed with regard to the core technical issues of dynamic modeling, numerical integration, measurement modeling, estimation method, measurement system as well as appropriate data formatting for the interoperability among foreign agencies. Although only the preliminary development strategy has been discussed through this work, the proposed strategy will aid the Korean astronautical society while on the development phase of the future Korea’s own lunar orbit determination system. Also, it is expected that further detailed system requirements or technical development strategies could be designed or established based on the current discussions.

In this study, we present preliminary results of precise orbit determination (POD) using satellite laser ranging (SLR) observations for International Laser Ranging Service (ILRS) Associate Analysis Center (AAC). Using SLR normal point observations of LAGEOS-1, LAGEOS-2, ETALON-1, and ETALON-2, the NASA/GSFC GEODYN II software are utilized for POD. Weekly-based orbit determination strategy is applied to process SLR observations and the post-fit residuals check, and external orbit comparison are performed for orbit accuracy assessment. The root mean square (RMS) value of differences between observations and computations after final iteration of estimation process is used for post-fit residuals check. The result of ILRS consolidated prediction format (CPF) is used for external orbit comparison. Additionally, we performed the precision analysis of each ILRS station by post-fit residuals. The post-fit residuals results show that the precisions of the orbits of LAGEOS-1 and LAGEOS-2 are 0.9 and 1.3 cm, and those of ETALON-1 and ETALON-2 are 2.5 and 1.9 cm, respectively. The orbit assessment results by ILRS CPF show that the radial accuracies of LAGEOS-1 and LAGEOS-2 are 4.0 cm and 5.3 cm, and the radial accuracies of ETALON-1 and ETALON-2 are 30.7 cm and 7.2 cm. These results of station precision analysis confirm that the result of this study is reasonable to have implications as preliminary results for administrating ILRS AAC.