GPS Satellite Research Articles

Prospects of GENESIS and Galileo joint orbit and clock determination

The European Space Agency (ESA) is preparing a satellite mission called GENESIS to be launched in 2027 as part of the FutureNAV program. GENESIS co-locates, for the first time, all four space geodetic techniques on one satellite platform. The main objectives of the mission are the realization of the International Terrestrial Reference Frames and the mitigation of biases in geodetic measurements; however, GENESIS will remarkably contribute to the determination of the geodetic parameters. The precise GENESIS orbits will be determined through satellite-to-satellite tracking, employing two GNSS antennas to observe GPS and Galileo satellites in both nadir and zenith directions. In this research, we show results from simulations of GENESIS and Galileo-like constellations with joint orbit and clock determination. We assess the orbit quality of GENESIS based on nadir-only, zenith-only, and combined nadir–zenith GNSS observations. The results prove that GENESIS and Galileo joint orbit and clock determination substantially improves Galileo orbits, satellite clocks, and even ground-based clocks of GNSS receivers tracking Galileo satellites. Although zenith and nadir GNSS antennas favor different orbital planes in terms of the number of collected observations, the mean results for each Galileo orbital plane are improved to a similar extent. The 3D orbit error of Galileo is improved from 27 mm (Galileo-only), 23 mm (Galileo + zenith), 16 mm (Galileo + nadir), to 14 mm (Galileo + zenith + nadir GENESIS observations), i.e., almost by a factor of two in the joint GENESIS + Galileo orbit and clock solutions.

A simplified GNSS/LEO joint orbit determination method

Joint orbit determination (JOD) of Global Navigation Satellite Systems (GNSS) satellites and low Earth orbit (LEO) satellites has been substantiated as an efficacious approach to compensate for the ground station geometry. In the traditional JOD, the precise orbit determination (POD) of LEO satellite is mainly processed by the reduced-dynamic approach, however, this approach involves complex calculations and the credibility of the determined positions diminishes when LEO satellites undergo orbital maneuvers. Therefore, a simplified JOD method is designed that employs kinematic approach to determine the LEO satellites orbit. To verify the effectiveness of the proposed method, the orbit of GPS satellites and LEO satellites are jointly estimated utilizing the regional and global networks. 8 LEO satellites, including GRACE-C/D, SWARM-A/B/C, SENTINEL-3A/B, and JASON-3, are chosen for JOD. The comparative analysis between the proposed method and traditional method are achieved in terms of GPS orbit accuracy, LEO orbit accuracy, computation time and the JOD performance during LEO maneuvers. Under regional station scenario, the GPS orbit accuracy determined using the proposed method and the traditional method is 3.64 cm and 2.52 cm, respectively. In the case of global station scenario, the accuracies are 1.71 cm and 1.64 cm. Additionally, the traditional method yields superior enhancement and higher accuracy of the LEO orbits. However, it exhibits a noticeable increase in computation time compared to the proposed method and the performance of JOD declines significantly when LEO satellites undergo orbital maneuvers. Alternatively, although the accuracy of the LEO orbits using the proposed method is comparatively lower, it offers a substantial reduction in the overall network computation time compared to traditional method. Moreover, the proposed method based on LEO kinematic precise orbit determination (KPOD) is nearly unaffected by orbital maneuvers of LEO satellites, presenting unique advantages in practical data processing.

The effect of satellite code bias on GPS PPP time transfer

ABSTRACT Time transfer technology using PPP relies onhigh-precision clockproducts of GPS satellites. When the observation we use is different that of the clock products, a SCB will be incorporated into the calculation of clock.For analyze the effect of SCB, we first derive the principles of correcting SCB. Thenthe effect of SCB on GPS PPP time transfer under three different conditionsis analyzed. The results indicate that the SCB can significantly improve the frequency stability, precision and PPP time transfer convergence speed.The precision of time transfer is improved by approximately 7% and the 5000s-200000s frequency stabilityis improved significantly, up to 20%.

Power Spectral Characteristics of In‐Situ Irregularities and Topside GPS Signal Intensity at Low Latitudes Using High‐Sample‐Rate Swarm Echo (e‐POP) Measurements

AbstractIonospheric density structures at low latitudes range in size from thousands of kilometers down to a few meters. Radio frequency (RF) signals, such as those from global navigation satellite systems, that propagate through irregularities suffer from rapid fluctuations in phase and intensity, known as scintillations. In this study, we use the high‐sample‐rate measurements of the Swarm Echo (CASSIOPE/e‐POP) satellite's GPS Occultation (GAP‐O) receiver taken after its antenna was re‐oriented to vertical‐pointing, simultaneously with e‐POP Ion Mass Spectrometer surface current observations as a proxy for plasma density, to obtain the spectral characteristics of GPS signal intensity and in‐situ irregularities at altitudes from 350 to 1,280 km. We show that the power spectra of both measurements can generally be characterized by a power law. In the case of density irregularities, the spectral index with the highest occurrence rate is around 1.7, which is consistent with previous studies. Also, all the power spectra of GPS signal intensity in this study show a single spectral index near 2. Moreover, roll‐off frequencies estimated in this work range from 0.4 to 2.5 Hz, which is significantly higher than Fresnel frequencies calculated from ground GPS receivers at low latitudes (between 0.2 and 0.45 Hz). Part of this increase is due to the 8 km/s orbital velocity of Swarm Echo near perigee. Another key difference is that variations in the GPS signals in this study are dominated by the topside ionosphere, whereas GPS signals received from ground are affected mostly by the relatively dense F‐region plasma in the 250–350 km altitudinal range.

An Autonomous Navigation Method for Planetary Rover Based on Multi-modal Fusion and Multi-factor Graph Optimization

Planetary surface is a complex unstructured environment with strong light, strong shadows, weak textures and numerous obstacles. There is no GPS signal support and satellite networking support on the planetary surface. The limited computer and sensor capabilities of planetary rover have a significant gap from the ground autonomous driving configuration, which presents a great challenge to the autonomous navigation methods of planetary rover. Extreme lighting conditions will increase the camera mismatching rate, and poor texture will increase the lidar matching error, which will reduce the estimation accuracy. Moreover, camera and lidar produce great motion distortion in the rugged terrain and bumpy environment, which leads to a large increase in the cumulative error of odometer. Camera and lidar are not suitable as the main sensor of SLAM algorithm in unstructured environment because they are sensitive to the change of environment. Therefore, an autonomous navigation method for planetary rover based on multi-modal fusion and multi-factor graph optimization is presented in this paper. The IMU (Inertial Measurement Unit) odometer node is added to the lidar odometer and visual odometer thread, and the IMU is taken as the central node to build a multi-factor graph optimization model. A multi-factor graph optimization model and a strong adaptive constraint strategy with reasonable weights are constructed. The pose estimation of other sensors is used to constrain IMU bias, meanwhile historical rut tracking, obstacle contour matching and skyline feature matching are introduced. Finally, motion prediction is achieved in the IMU odometer node. Simulation results and field tests show that this method can effectively cope with the unstructured environment on planetary surface, which achieves adaptive and robust autonomous navigation of planetary rover.

GPS-Based Precise Orbit Determination of LEO Satellites Using Space-Based Double-Differenced Observations

One of the important GPS applications in space is precise orbit determination (POD) of Low-Earth Orbit (LEO) satellites. Thousands of LEO satellites are currently in orbit. One of the challenges is how to efficiently and precisely determine the orbits to satisfy the relative and absolute accuracy needs of missions. Currently, GPS-based POD of LEO satellites can be performed using either un-differenced (UD) or ground-based double-differenced (DD) observations. The UD POD needs both precise GPS satellite orbits and clocks; the DD POD needs not only the precise orbits, but also global ground reference receivers. Therefore, the GPS-based LEO POD is based on either global ground stations or precise GNSS clocks, which are not convenient for near real-time or real-time data processing. The GPS orbits can be precisely predicted in certain time; the clocks are not. For some formation flying satellite missions (two or more LEO satellites), the absolute orbit accuracy requirements are not as stringent as the relative requirements. The problem is how to perform LEO POD without global station data and precise GPS clocks in near real-time or real-time to achieve the mission orbit accuracy requirement. Based on this motivation, we investigated the GPS-based LEO POD using space-based DD observations without using global ground station GPS data and precise GPS clock products. For this study, we processed the real GPS observations from two LEO satellites. The absolute and relative orbit accuracy is assessed using several tests, including analysis of the orbit fits, external orbit comparisons, Satellite Laser Ranging (SLR) and K-band Ranging (KBR) residuals.

Introduction of the Runge-Kutta method in GPS orbit computation

In all Global Navigation Satellite Systems (GNSS) applications, the determination of the satellite orbits is an important task. In this study, we present the equations given in the Interface Specification Document of GPS and the Runge-Kutta method in the computation of the position P, velocity V, and acceleration A of the GPS satellites using the broadcast ephemeris. The definition of the differential equation describing the GPS satellite's motion has enabled us to introduce the Runge-Kutta method in the GPS orbit computation; this method uses the initial conditions determined in this study from the Keplerian elements provided in the broadcast ephemeris files. The Lagrange interpolation method is used for comparison of the results, where the vectors P, V, and A are estimated using the precise ephemeris. The difference not exceeding 2.4 m was obtained in the X, Y, and Z axes during seven days on the position of the GPS satellite number 9 tested in this study. In velocity and acceleration, the difference is about a few mm/s and mm/s2, respectively.

Wheat height and phenology retrieval using GPS/BDS interferometric reflectometry technology

Monitoring vegetation growth utilizing the Global Navigation Satellite System-interferometric reflectometry (GNSS-IR) technology has obtained widespread attention in recent years. GNSS-IR based approaches demonstrate the advantageous characteristics of high temporal resolution while maintaining cost-effectiveness. Most existing studies generally utilize reflected signals to retrieve the plant height and vegetation water content, while few studies focus on vegetation phenology and yield. In this study, both phenology metric dates and plant height of wheat were extracted using reflected signals data from GPS and BeiDou Navigation Satellite System (BDS) in a wheat farmland, located at northeast of Henan, China. The results showed that there was reasonable consistency between the retrieval heights and in-situ measurements, with the correlation coefficient of 0.836 for GPS and 0.868 for BDS. The damping coefficient calculated by signal-to-noise ratio (SNR) presented obvious characteristics during different wheat growth stages and the trends of the damping coefficient were basically consistent with that of the normalized difference vegetation index (NDVI). When taking a lag between the leaf area index (LAI) and the damping coefficient into consideration, the correlation was greatly increased. Moreover, the phenology metrics dates extracted by the damping coefficient were close to those estimated by NDVI with errors less than 8 days. Hence, this research preliminarily proves the feasibility of obtaining the wheat phenological period by using GNSS-IR technology.

Coseismic and Early Postseismic Deformation Mechanism Following the 2021 Mw 7.4 Maduo Earthquake: Insights from Satellite Radar Interferometry and GPS

Exploring the deformation mechanism of the 2021 Mw 7.4 Maduo Earthquake is crucial for better understanding the seismic hazard of the faults with low strain rates inside the Bayan Har block. This study leverages deformation information derived from Sentient-1 A/B images and GPS data to investigate in detail the co- and postseismic deformation mechanisms using multiple methods. The main results are as follows. First, the postseismic InSAR time series robustly identified the reactivation of the Changmahe fault, indicating the impact of the Maduo event on surrounding active faults. Second, the joint inversion of Interferometric Synthetic Aperture Radar and GPS revealed that (1) there was a complementary and partially overlapping relationship between the coseismic slip and postseismic afterslip of the main rupture; and (2) the Changmahe fault exhibited thrust compression dislocation in the early stage and experienced a sustained compressive effect from afterslip in the one year after the mainshock. Third, modeling the processes of viscoelastic relaxation and poroelastic rebound revealed that the postseismic deformation was probably caused by a combination of afterslip (near-field) and viscoelastic relaxation (near and far field). Fourth, the stress changes driven by the Maduo event revealed that the seismic gaps inside the Maqin-Maqu segment and the Kunlun Pass-Jiangcuo fault will be potential seismic risks in the future.

Precise Positioning of Primary System of Geodetic Points by GNSS Technology in Railway Operating Conditions

This article deals with the analysis of the accuracy of the geodetic real-time GNSS measurement procedure used in railway operating conditions in the Czech Republic. The purpose was to determine to what extent the operating conditions affect the accuracy of the measurement result and whether an accuracy of standard deviation σx,y = 5 mm in the horizontal plane could be achieved. The use of geodetic GNSS equipment with an IMU unit was also tested. The accuracy obtained in operational conditions is compared with the accuracy obtained on a calibration base using the same measurement procedure. The consistency between the accuracy of the primary system (satellite-based) and the secondary system (terrestrially measured by the traverse method) is also discussed. The analysis includes the issue of residual inhomogeneities of the uniform transformation key when converted to the Czech national coordinate system S-JTSK. It is shown that a homogeneous accuracy in coordinate standard deviation better than σx,y = 5 mm can be achieved. The results indicate that the accuracy under operational conditions is two–three times worse than the accuracy achieved by the same procedure under ideal conditions on a calibration base. This is due to the non-ideal observing conditions, i.e., horizon occlusion by overlays, surrounding vegetation and multipath effects. It has been shown that the effect of multipath can be reduced by repeating short observations 3–4 h apart. Older GNSS instruments using an IMU unit in combination with an electronic compass (eCompass) are at risk of a systematic bias of up to several tens of millimeters, which can be detected by rotating the antenna by 180°. The current uniform transformation key used in the Czech Republic for the conversion of GNSS coordinates into the national system has residual geometric inhomogeneities (p = 0.90 to 10 mm/km, sporadically up to 20 mm/km), which metrologically deteriorate the results of the calculation of the terrestrially measured secondary system inserted into the GNSS measured primary system. Achieving homogeneous accuracy in coordinate standard deviation in a horizontal plane better than σx,y = 5 mm has been demonstrated in non-ideal railway operating conditions with increased risk of multipath. The innovative aspect of the approach used is that it simplifies and thus increases the efficiency of the measurement with respect to the availability of GPS, GLONASS, Galileo and BeiDou satellites, as well as reducing the effect of multipath on the noise by repeating the measurement procedure.

Determination of total electron content in the ionosphere over the territory of the Republic of Belarus based on global navigation satellite systems data

We present the results of experimental studies of electron content in the ionosphere over the territory of the Republic of Belarus based on data from global navigation satellite systems. The results of measurements of the precise positioning system of the Republic of Belarus and navigation data of GPS satellites in RINEX format were used as input data. Expressions for calculation of the total electron content using the two-frequency method and a combination of measurements by phase and code delays are given. Algorithms for eliminating cycle slip and determining differential code biases are used. Examples of calculating the vertical electron content over the Republic of Belarus at different moments of time are demonstrated. The obtained results are reasonable to use in monitoring of the ionosphere in order to provide reliable operation of radio systems, detection of ionospheric anomalies of natural and artificial origin, as well as forecasting of natural phenomena on their basis.

Performance comparison of orbit and clock augmentation corrections from PPP-B2b, HAS and CLAS

Precise Point Positioning (PPP) augmentation corrections can enhance the accuracy of Global Navigation Satellite System (GNSS) positioning. Three prominent PPP augmentation services: the third generation of BeiDou Navigation Satellite System (BDS-3) PPP-B2b, Europe's Galileo High Accuracy Service (HAS) and Quasi-Zenith Satellite System (QZSS) Centimeter Level Augmentation Service (CLAS), are evaluated and compared, contributing to the optimal utilization of these services. Based on a 25-day dataset in March 2023, all the three services have made significant improvements in GPS orbit accuracy compared to the broadcast ephemeris. Galileo HAS achieved the highest three-dimension (3D) orbit accuracy for both GPS and Galileo (10.5 and 10.4 cm), followed by CLAS (48.1 and 16.3 cm) and PPP-B2b (GPS: 73.0 cm, BDS-3: 44.6 cm). This distinct performance of PPP-B2b orbits is found for the first time to be related to the large constellation rotation errors, which accounts for roughly half of the orbit errors of along and cross components. Regarding clock errors, PPP-B2b and HAS enhanced GPS clock accuracy to 0.20 ns (70 %) and 0.17 ns (74 %), respectively, over the broadcast clock (0.66 ns), while CLAS shows a much smaller improvement. Moderate correlations of 0.3–0.8 between orbit and clock errors are observed for the three augmentation corrections. For signal in space ranging error (SISRE), HAS (GPS:5.5 cm, Galileo:4.9 cm) outperforms PPP-B2b (GPS:14.3 cm, BDS-3:8.9 cm) and CLAS (GPS:23.6 cm, Galileo:13.7 cm). The 2-hourly kinematic positioning results with an accuracy of better than 20 cm (95th percentile) for the horizontal and vertical component, could be achieved for both PPP-B2b and HAS within 20 min within their service area. It is not recommended rely solely on the orbit/clock/code bias of CLAS for precise positioning, as the positioning accuracy is influenced by the SISRE. PPP-RTK of CLAS achieves horizontal 10 cm (95th percentile) and rapid convergence of mostly less than a minute by incorporating phase bias and atmospheric corrections.

A method to assess the quality of GNSS satellite phase bias products

As part of the International GNSS Service (IGS), several analysis centers provide GPS and Galileo satellite phase bias products to support precise point positioning with ambiguity resolution (PPP-AR). Due to the high correlation with satellite orbits and clock offsets, it is difficult to assess directly the precision of satellite phase bias products. Once outliers exist in satellite phase biases, PPP-AR results are no longer reliable and the combination of satellite phase bias products from IGS analysis centers also gets difficult. In this contribution, we propose a method independent of ground measurements to detect outliers in satellite phase biases by computing the total Difference of satellite Orbits, Clock offsets and narrow-lane Biases at the midnight epoch between two consecutive days. Results over 180 days show that about 0.2, 1.1, 2.0 and 0.1% of the total DOCB values for GPS satellites exceed 0.15 narrow-lane cycles for CODE final, CODE rapid, CNES/CLS final and WUHN rapid satellite products, respectively, while the same outlier-ratios for Galileo satellites are 0.1, 0.9, 0.4 and 0.1%, respectively. As an important contribution to the orbit, clock and bias combination task, we check the consistency of satellite phase bias products between two analysis centers before and after removing these detected outliers from individual analysis centers. It is convincing that the number of large differences of satellite phase biases between two analysis centers is notably reduced.

Golden Eagle Populations, Movements, and Landscape Barriers: Insights from Scotland

GPS satellite tracking allows novel investigations of how golden eagles Aquila chrysaetos use the landscape at several scales and at different life history stages, including research on geographical barriers which may prevent or limit range expansion or create population/sub-population isolation. If there are significant barriers to golden eagle movements, there could be demographic and genetic consequences. Genetic studies have led investigations on the identification of sub-species, populations, and sub-populations but should be conjoined with demographic studies and dispersal movements to understand fully such designations and their geographic delimitation. Scottish eagles are genetically differentiated from continental European birds, with thousands of years of separation creating a distinct population, though without sub-species assignation. They present unique research opportunities to examine barriers to movements illustrated by satellite tracking under Scotland’s highly variable geography. We primarily examined two features, using more than seven million dispersal records from satellite tags fitted to 152 nestlings. The first was the presence of unsuitable terrestrial habitat. We found few movements across a region of largely unsuitable lowland habitat between upland regions substantially generated by geological features over 70 km apart (Highland Boundary Fault and Southern Uplands Fault). This was expected from the Golden Eagle Topography model, and presumed isolation was the premise for an ongoing reinforcement project in the south of Scotland, translocating eagles from the north (South Scotland Golden Eagle Project: SSGEP). Second was that larger expanses of water can be a barrier. We found that, for a northwestern archipelago (Outer Hebrides), isolated by ≥24 km of sea (and with prior assignation of genetical and historical separation), there were no tagged bird movements with the Inner Hebrides and/or the Highlands mainland (the main sub-population), confirming their characterisation as a second sub-population. Results on the willingness of eagles to cross open sea or sea lochs (fjords) elsewhere in Scotland were consistent on distance. While apparently weaker than the Outer Hebrides in terms of separation, the designation of a third sub-population in the south of Scotland seems appropriate. Our results validate the SSGEP, as we also observed no movement of birds across closer sea crossings from abundant Highland sources to the Southern Uplands. Based on telemetric results, we also identified where any re-colonisation of England, due to the SSGEP, is most likely to occur. We emphasise, nevertheless, that our study’s records during dispersal will be greater than the natal dispersal distances (NDDs), when birds settle to breed after dispersal, and NDDs are the better shorter arbiter for connectivity.

Using Reinforcement Learning to Iteratively Construct Road Networks from Satellite Images and GPS Data

Constructing road networks manually is a time consuming and labor-intensive process. This paper proposes a new method to iteratively construct road networks using reinforcement learning from a combined tensor-based representation of satellite image and GPS trajectory data.

Disturbances of GLONASS and GPS signals during magnetic storm on March 23—24, 2023, according to observations on the Kola peninsula

The growth of scintillations of GLONASS and GPS satellite signals using the Septentrio GNSS receiver installed in the city of Apatity during a strong magnetic storm on March 23—24, 2023, is analyzed. According to the ionosonde data at the Lovozero station and the data of the EISCAT radar in Tromsø, it was shown that the growth of phase scintillations is caused by an increase in the plasma concentration mainly in the E-layer of the ionosphere. The growth of phase scintillations is accompanied by the appearance of discrete forms of auroras.

Simultaneous monitoring of the limited area ionosphere with the use of GPS and ionosonde

The importance of this study lies in its contribution to the field of ionospheric science, particularly in the monitoring and prediction of ionospheric parameters. This study focuses on comparing two techniques for spatio-temporal kriging of critical ionospheric parameters, namely the F2 layer critical frequency (foF2) and Total Electron Content (TEC), under both quiet and disturbed ionospheric conditions. The first technique, Ordinary Kriging (OK), relies solely on foF2 data, while the second technique, Kriging with External Drift (KED), incorporates the dynamic relationship between foF2 and TEC data to improve the interpolation process. This study demonstrates how KED, which utilizes dynamic secondary information, offers advantages over OK. By simultaneously monitoring the ionosphere from a single ionosonde station and collecting data from visible GPS satellites, KED enhances the accuracy of spatial structural analysis, enabling more precise estimations of ionospheric characteristics across interpolation grid nodes. The research highlights that the non-uniform distribution of ionosonde stations, especially in regions with limited foF2 data, can benefit from KED's incorporation of TEC data. Additionally, the study reveals the potential for improving TEC maps by assimilating foF2 data located at the periphery of the observation area. Furthermore, the paper discusses the determination of semivariogram structures for both foF2 and TEC data, emphasizing the spatial correlations between data points. These findings contribute to a better understanding of ionospheric behavior and the potential for enhanced ionospheric modeling.

Mathematical simulation of radar scattering characteristics of ultra-short-pulse signals on complex shape objects

The article discusses the method of mathematical modeling of characteristics of scattering of ultrashort-pulse signals of complex shape objects, taking into account the shape of the envelope of the ultrashort-pulse radio signal and the features of radiation of an ultrawideband antenna. To estimate the reflecting properties of objects, the determination of the ultra-short-pulse radar cross section of complex shape objects with their radar location by non-harmonic signals is given. The mathematical model of radar scattering characteristics of location objects under ultra-short-pulse signals takes into account the shape of the envelope of the signal generated by a transmitter of a single-chip microchip NVA 6100 with duration of 300 ps and obtained through direct measurements. The electric scattering field of a location object in the frequency domain is determined on the use of a polygonal model of the object and asymptotic methods of diffraction theory. The transfer function of the Vivaldi antenna to radiation obtained by numerical modeling using the method of moments was used in the electric field calculations. Verification of the developed mathematical model was carried out on the basis of determining the ultra-short pulse radar cross section of perfect electric conductive sphere. A diagram of the ultra-short-pulse radar cross section of the GPS satellite was obtained, taking into account the radiation of the ultra-short-pulse signal by the Vivaldi antenna. The diagram of ultra-short pulse radar cross section of GPS satellite does not contain interference maxima and minima, and the averaged values ultra-short radar cross section of GPS satellite from lateral angles are on average 3..8 times less than radar cross section at the center frequency of the spectrum. The mathematical model of radar scattering characteristics of location objects under ultra-short-pulse signals is intended for preliminary estimates of the ultra-short-pulse radar cross section at an early step of designing ultra-wideband radars. The obtained estimates of the ultra-short-pulse radar cross section of the GPS satellite make it possible to identify the main features in the scattering of ultra-short-pulse signals at the location of extended targets.

Thermospheric Density Estimation Method Using a First-Order Gauss–Markov Process

Low-Earth-orbit (LEO) spacecraft are significantly influenced by atmospheric drag. Accurately estimating thermospheric density is pivotal for the precise calculation of drag acceleration. However, thermospheric density along a specific orbit, computed using existing thermospheric models, has certain inaccuracies. In this work, a first-order Gauss–Markov process is used to model the deviation of atmospheric drag acceleration. With the Markov parameter of the initial state iteratively computed through sequential estimation and the smoothing method, the thermospheric density is derived from high-precision GPS measurements. In simulation scenarios, the root-mean-square error and relative error of the estimated thermospheric density reduce by about 45 and 50% relative to the prior density, respectively. Using the estimated density for orbit propagation, satellite trajectories’ one-day position and velocity error are, respectively, within 100 m and 0.1 m/s, and an average improvement in orbit precision is over 80%. The proposed method has been applied to the real Tsinghua Science Satellite (Q-SAT) GPS measurements for effectiveness verification. It shows strong adaptability under extreme space weather and during the occurrence of geomagnetic storms. Due to the estimated Markov parameter of the initial state obeying the Langevin dynamics properties, the proposed method also offers short-term thermospheric density forecasting potential.

Yere Eşlenik Uydularda Tümleşik Uzay Tabanlı Alıcılar Kullanarak Yörünge Belirleme

The use of space-born GNSS receiver onboard geosyncronous communication satellite is an effective alternative orbit determination technique compared to traditional methods, although the communication satellite operator and manufacturers need some proof. In this study, GPS-based range measurement via space-born receivers’ onboard communication satellite platforms was proposed for precise orbit determination of geostationary satellites. Thus this research investigates how proposed method can be applied to obtain reliable and sustainable results. The analysis was performed to find out the accuracy and validity of the proposed method. Moreover, the results were compared with the simulation of four satellites’ orbital parameters at different orbital slots. The sequential process (SP) Kalman filter method was chosen to process the measured data. The processing method was kept the same for all analyses to investigate the gathered measurement data effects. The results show that the average position differences between the determined and actual orbit are 17.613 m, 13.637 m and 14.998 m, for in-track, cross-track, and radial direction. The average in-track, cross-track, and radial velocity differences are 0.339 cm/s, 0.160 cm/s, and 0.066 cm/s respectively. The results show that the proposed method is approved in communication satellites orbit determination. The utilization of GPS satellites is a promising approach for geostationary satellite orbit determination. The communication satellites operators can utilize the GPS-based orbit determination with its benefits.