Assessment of Swarm Kinematic Orbit Determination Using Two Different Double-Difference Methods

: The Swarm mission aims to study the principle and change regularities of the Earth’s magnetic ﬁeld. Precise orbit determination is essential to the successful implementation of the mission and relevant scientiﬁc research. This article focuses on using two different double-difference methods to improve the accuracy of Swarm kinematic orbit determination. The accuracy of the kinematic orbit determination relies entirely on the space-borne observation data, independent of any dynamic parameters. The article analyzes the data quality of the Swarm space-borne global positioning system (GPS) receiver and presents a detailed introduction to the data pre-processing algorithms. The double-difference observation gathering and the applied orbit determination strategy using two different double-difference methods are discussed. The results of the kinematic orbits under different solar cycle conditions are presented, along with an evaluation based on analysis of GPS carrier phase residuals, subtracting from the post-processed orbits, and assessment with satellite laser ranging (SLR) measurements. The results show that the accuracy of the kinematic orbit determination is at the centimeter level for the three Swarm satellites’ orbit solutions. The daily root mean square (RMS) values of the three satellites’ phase residuals remain at around the 6 mm level. The RMS values of the position residuals between the kinematic orbits and the reduced dynamic orbits released by the European Space Agency (ESA) are at about the 2–3 cm level. The external evaluation with SLR measurements shows a good agreement with the ESA level, with the RMS values of the SLR residuals for kinematic orbits around 2 cm.


Introduction
The computation of precise orbits for satellites in the low earth orbit (LEO) is an essential precondition to complete important scientific tasks such as ocean monitoring, earth gravity field model recovery, electromagnetic field exploration, topographic mapping, and atmospheric and space environment detection [1].LEO satellite precise orbit determination (POD) can be divided into three approaches: kinematic orbit determination, dynamic orbit determination, and reduced dynamic orbit determination.During the past few years, since the continuing evolution of space-borne global positioning system (GPS) technology, the accuracy of kinematic precision orbit determination can reach the centimeter level [2][3][4][5][6][7][8][9].It has been demonstrated that using kinematic positions is better suited for independent gravity field recovery.This is because kinematic positions are obtained from GPS measurements taken at different times, which allows for a larger separation between the positions and reduces the correlation between the observations.This makes the data more appropriate for independent analysis.In addition, kinematic orbit determination is a purely geometrical algorithm to study the motion of satellites.This algorithm does not use any information from the gravity field or LEO dynamic parameters.By contrast, the anomalous dynamic parameters will assimilate some gravitational orbital disturbances and impact the kinematic positions [10,11].
The European Space Agency (ESA) 's Earth Explorer Swarm mission consists of three satellites with different orbital altitudes [12].The constellation was successfully launched on 22 November 2013 and provided Earth observations.The mission was granted an extension until the end of 2021.Swarm is dedicated to producing a comprehensive and precise examination of the magnetic field, including its variation over time.This is achieved through a constellation of satellites equipped with advanced magnetometers and other cutting-edge instruments [13].Swarm A and Swarm C were initially positioned to fly in close proximity to each other at the equator, separated by a distance of 150 km.Both satellites were placed in orbit at an altitude of approximately 460 km.Meanwhile, Swarm B was placed in a higher orbit, with an initial altitude of about 510 km.In fact, as the extension of time, all three satellites' orbits will undergo a natural decay after experiencing drag, the decay will be much more pronounced for the lower pair satellites [14].Each satellite has a set of seven core instruments installed on it.Theseinclude the vector and scalar magnetometers for measuring the Earth's magnetic field's direction and strength, the Langmuir probe and thermal ion imager to measure the Earth's electric field, the accelerometer for measuring the non-gravitational accelerations affecting the Swarm satellites, and the GPS receiver and laser retro-reflectors for orbit determination and validation [15].
Each of the Swarm satellites is loaded with two dual-frequency eight-channel GPS receivers manufactured by RUAG Space, a main and a redundant one [14].The ESA provides Level 1b and Level 2 data products related to orbits with high accuracy.The Level 1b products contain space-borne GPS observation data and satellite attitude data.The Level 2 products contain corrected and calibrated accelerometer observations and time series of precise science orbits (PSOs) of each satellite.The PSOs consist of kinematic orbit and reduced dynamic orbit solution, and both orbits are calculated by the GPS high-precision orbit determination software tools (GHOST) [16] combined with observation files and external data from IGS, IERS, and CODE.The reduced-dynamic orbit products offer both position and velocity solutions with a smooth 10-s integer sampling, whereas the kinematic orbit products solely offer position information with non-integer time stamps [14].The kinematic orbits have a sampling rate of 1 Hz and depend on the frequency of the GPS receiver.Additionally, a brief description detailing the accuracy of the two orbit types is contained in the Level 2 product.The accuracy of the final reduced dynamic orbit is less than 3 cm but is released at least 13 days later (generally 21 days) due to the use of precise ephemeris (ultrafast-3 h delay, fast-17 h delay, and final-13 days delay).
After the launch of Swarm satellites, Swarm POD has gained a great deal of attention.Zangerl et al. carried out a performance evaluation of the Swarm GPS receiver, presented the in-orbit measurement data, and compared it with ground-test results [17].Several modifications have been made to the onboard receiver settings to further boost the tracking performance, thereby possibly increasing the quality of the results for both kinematic orbits and reduced-dynamic orbits [18].Van den IJssel et al. explained the adopted POD strategies for the reduced dynamic approach and kinematic approach and presented the first precise science orbits based on over a year of GPS observations [14].Zehentner developed an approach for kinematic POD with raw GPS measurements and validated that the proposed method can produce orbit positions as accurate or better than existing approaches [10,11].Allende-Alba et al. provided a brief overview of the half-cycle ambiguity resolution process utilized in generating carrier phase observations and implemented the strategies for integer ambiguities fixing and the space baseline determination between Swarm A and Swarm C [19,20].Zhang proposed kinematic and reduced dynamic POD methods for the Swarm satellites based on space-borne GPS precise point positioning measurements, and then the acceleration approach was adopted to recover the gravity field model [21,22].Montenbruck et al. showed the PSO results from four years of GPS observation data by using a refined model and reprocessing technique.Following an assessment from the analysis of SLR observation, and overall, a 50% performance improvement can be inferred for the kinematic orbits and a 30% performance improvement can be inferred for the reduced-dynamic orbits [23].Kermarrec et al. proposed a new way using the Matérn covariance family to reduce the influence of ionospheric scintillations [24] on Swarm orbit determination [25].Ren and Schön developed software using precise point positioning batch least-squares estimation for kinematic POD and computed the 1-year orbit results of Swarm satellites with the software [26].Li et al. processed the space-borne GPS data from Swarm and Sentinel-3 satellites via the undifferenced ambiguities resolution method for real-time LEO kinematic orbit determination [27].Bezděk et al. showed that GPS data quality and calculation results from the Swarm mission have similar accuracy as the GRACE mission, thus the Swarm mission could fill the blank in gravity field observation before the launch of GRACE Follow-On in August 2017 [28][29][30][31][32]. Encarnação et al. demonstrated the best combination method for generating an optimal model based on the products provided by three institutes (ASU, Bezděk et al., 2016 [28]; AIUB, Jäggi et al., 2016 [30]; IfG, Zehentner et al., 2015 [10]).The quality of the Swarm-derived gravity field model is proven to be good [33].Lück et al. utilized Swarm data from 37 months of kinematic orbit to derive and investigate time-variable gravity field and ocean mass variations [34].
In GPS data processing, the double-difference model is one of the most commonly used basic methods [35].However, with the increase in multi-frequency, multi-mode, and station network scale, the calculation time of the double-difference model increases exponentially.
To break the computational bottleneck of the double-difference model, some scholars have proposed precise point positioning (PPP) based on the zero-difference model [36].PPP can obtain the absolute coordinates of the target by single-station processing, which significantly reduces the calculation time of data processing.For most of the research mentioned above, the POD strategies are based on zero-difference approaches.With the rapid improvement of computer performance in recent years, the double-difference model is preferred for solving in the absence of extra-long baselines.Hwang et al. used double-difference methods for processing KOMPSAT-2 and SAC-C satellite GPS observations, which satisfied the requirement of 1-m accuracy [37].Švehla et al. used the double-difference method for the CHAMP satellite [3].Jäggi et al. discussed the impact of fixed ambiguities on the double-difference baselines and obtained good orbital results for GRACE through different double-difference baseline solutions [4].Kang et al. used zero-difference and doubledifference observation values for reduced dynamic orbit determination for GRACE and GRACE-FO satellites and analyzed the results, the three-dimensional accuracy of the two orbits was better than 2 cm [38,39].
In this article, we focus on increasing the accuracy of Swarm kinematic precise orbit determination using two different double-difference methods.The orbit determination process with this method is a bit complicated and depends on the quality of the observation data.The satellite-borne GPS data in 2021 were processed to validate the proposed methods under different solar cycle conditions.The year 2021 is currently in the upward period of the solar cycle 25 [40], and the level of solar activity is gradually increasing, with a consequent increase in ionospheric activity.The deterioration of signal reception caused by the ionosphere can significantly reduce the signal-to-noise ratio, making it difficult to capture and track signals, thereby reducing positioning accuracy [25,41].In severe cases, strong ionospheric activity can cause the navigation receiver to lose lock and signal interruption.The first ground enhancement of the 25th solar cycle occurred on 28 October 2021 [42].Therefore, a special analysis was conducted on the impact of ionospheric activity.This article is divided into the following sections to explain.First, Section 2 introduces the double-difference observation gathering and satellite-borne GPS data pre-processing.Then, the applied POD strategy using the double-difference method is discussed.Afterward, the results of the orbit calculations are presented in Section 3, and the accuracy and reliability of the obtained results are evaluated through three methods: analysis of space-borne GPS phase observation residuals, subtracting from the ESA's PSOs, and assessment with satellite laser ranging measurements.Section 4 is about an analysis and discussion of the results.Finally, a summary is given in Section 5.

Observation Analysis
This article analyses the quality of GPS observations offered by space-borne GPS receivers which provides a basis for data pre-processing.At the beginning of the mission, the Swarm constellation provided 10-s GPS observation, but no simultaneous GPS measurement epoch, which complicates the baseline determination between the Swarm constellation.Therefore, on 15 July 2014, the instrument setting was modified, and the three Swarm satellites delivered 1-s GPS observation.Each Swarm satellite has eight channels which means that the Swarm satellite supports tracking up to 8 GPS satellites simultaneously.For analyzing daily observation data, we calculated the number of tracked satellites and plotted them as a histogram.The total number of daily epochs approached to or was equal to 86,400.Figure 1 shows the GPS tracking performance for all three Swarm satellites from 1 January to 31 December 2021.On 1 November 2021, which is an exception, a malfunction occurred onboard of Swarm C. The quantity of available epochs was only 59,237 due to the anomaly.As indicated from the figure, the GPS receiver tracked eight satellites more than 80% of the time on most days of 2021.The average percent of tracking eight GPS satellites for the three Swarm satellites are 83.52%,82.38%, and 83.67%, respectively.The reduction in the number of visible satellites was accompanied by a decline in the quality of GPS observation data.We observed a decline in data quality during one day in late November and the second half of December.This could account for the medium-high solar activity.A specific analysis will be conducted in the following section.
borne GPS phase observation residuals, subtracting from the ESA's PSOs, and assessment with satellite laser ranging measurements.Section 4 is about an analysis and discussion of the results.Finally, a summary is given in Section 5.

Observation Analysis
This article analyses the quality of GPS observations offered by space-borne GPS receivers which provides a basis for data pre-processing.At the beginning of the mission, the Swarm constellation provided 10-s GPS observation, but no simultaneous GPS measurement epoch, which complicates the baseline determination between the Swarm constellation.Therefore, on 15 July 2014, the instrument setting was modified, and the three Swarm satellites delivered 1-s GPS observation.Each Swarm satellite has eight channels which means that the Swarm satellite supports tracking up to 8 GPS satellites simultaneously.For analyzing daily observation data, we calculated the number of tracked satellites and plotted them as a histogram.The total number of daily epochs approached to or was equal to 86,400.Figure 1 shows the GPS tracking performance for all three Swarm satellites from 1 January to 31 December 2021.On 1 November 2021, which is an exception, a malfunction occurred onboard of Swarm C. The quantity of available epochs was only 59,237 due to the anomaly.As indicated from the figure, the GPS receiver tracked eight satellites more than 80% of the time on most days of 2021.The average percent of tracking eight GPS satellites for the three Swarm satellites are 83.52%,82.38%, and 83.67%, respectively.The reduction in the number of visible satellites was accompanied by a decline in the quality of GPS observation data.We observed a decline in data quality during one day in late November and the second half of December.This could account for the medium-high solar activity.A specific analysis will be conducted in the following section.The true observation epoch of the Swarm constellation is expressed in GPS time and pretty close to an integer second.Therefore, the timestamp of the PSOs published by ESA is also pretty close to an integer second.In order to process single-and double-difference GPS data, an appropriate interpolation of IGS precise ephemeris and clocks is required.The high-order polynomial interpolation suffers from a potentially random noise process that is not appropriate for the clock parameters [43].We choose to implement a simple linear interpolation method.

Data Pre-Processing
Suppose that there are Swarm satellite receiver r and GPS satellite p, then, neglecting part of the errors, the simplified pseudo and phase basic observation equations are as follows where ρ p r stands for the distance between the GPS satellite p and the Swarm receiver r, measured by pseudo-range observations; t r is the receiver clock at the time of signal reception; t p is the satellite clock at the time of signal transmission; b r and b p refer to the code hardware delay of the receiver and satellite, whereas B r and B p are the receiverdependent and satellite-dependent phase biases, respectively; λ is the wavelength; N p r is the phase ambiguity, and I p r is the ionospheric delay; ε P and ε L stand for the sum of random noise.Since the altitudes of the three Swarm satellites are around 500 km, the tropospheric delay is not considered in the observation equation.
The ionosphere-free (IF) combination of pseudo-range and phase observations at a specific epoch can be expressed as  3) is assimilated into the clock offsets of the satellite and receiver, respectively, then eliminated.In addition, the phase ambiguity in Equation (4) will absorb both phase bias and the code hardware delay.
The Melbourne-Wübbena (MW) linear combination can eradicate the impact of geometric distance, ionospheric effects, GNSS satellite, and the receiver clock, including the wide-lane part and the narrow-lane part, and finally represented as where L p r,MW represents the MW linear combination; λ W represents the wavelength; N p r,W represents the ambiguity; ε MW is the sum of MW measurement noise and multipath errors.At two continuous observation epochs t i and t i+1 , the epoch-side phase single-difference observable is defined by where ∆ is the operator for epoch difference in corresponding parameters.The singledifference between epochs eliminates the ambiguities but ignores the correlation between the phase observations.If the number of epochs of the observed arc segment is less; then, the coordinate correction ∆X ∆Y ∆Z cannot be calculated.In addition, Swarm orbit is required to have a higher precision priori orbit.
With the extension of the observation period, the possibility of cycle slip formation is increased.First of all, the methods of cycle slip detection are briefly introduced below.According to the results of kinematic point positioning using the ionosphere-free (IF) combination pseudo range observation, we can calculate the position and receiver clock bias of the Swarm satellite X = X Y Z ct r .Then, the residual information is used to judge whether the pseudo-range observations are grossly inaccurate.Similarly, the various parameters of the position and receiver clock offset X = ∆X ∆Y ∆Z ∆ct r are estimated using a least squares algorithm, the cycle slips can be precisely detected in the current epoch through the IF combination epoch carrier phase single-difference observations.The observation equations of pseudo-range gross error and phase cycle slip can be modeled as follows where P is the IF combination pseudo-range, ∆L is the epoch carrier phase single-difference observations.X Y Z ct r and ∆X ∆Y ∆Z ∆ct r are the parameters to be estimated.In addition to identifying gross errors and detecting cycle slips in the data pre-processing, the difference observation data must be optimized and corrected.
Possible configurations for the double-difference approaches are illustrated in Figure 2. The double-difference observation data involves two GPS satellites, p and q, and two receiver stations.First, receiver-side single-difference between receiver stations is employed to eliminate the satellite clock error and satellite-dependent phase biases.Then, satelliteside double-difference is used to eradicate the receiver clock bias and receiver-dependent phase bias as shown in Figure 2

( ) ( )
, , , The moving baseline mode is used in Equation ( 9) because the positions of one Swarm satellite are not fixed.When both Swarm satellites are moving, one of the satellites is selected as the base station, and its position is determined through the single-point positioning method.The zero-difference IF combination phase observations can be expressed as Equation ( 4).The estimated parameters are positions and clock offsets.In order The moving baseline mode is used in Equation ( 9) because the positions of one Swarm satellite are not fixed.When both Swarm satellites are moving, one of the satellites is selected as the base station, and its position is determined through the single-point positioning method.The zero-difference IF combination phase observations can be expressed as Equation ( 4).The estimated parameters are positions and clock offsets.In order to ensure the accuracy of the reference orbit, numerical integration and iteration are required [21].
Compare reference positions with the final ESA orbits and eliminate the points with a difference greater than 10 cm.Then, the position of another Swarm satellite compared to the base station can be computed on an epoch-by-epoch basis.
When Equation ( 8) is used, double-difference baselines need to be formed between the Swarm satellite and all IGS ground stations.The positioning accuracy depends not only on the GPS satellite orbit but also on the coordinate accuracy of the ground station.Tropospheric delays must also be considered once ground stations are involved.Precise ephemeris for GPS satellite coordinates should be employed to reduce the broadcast ephemeris errors for the long-distance estimation between Swarm satellites and IGS ground stations.For the corrections not incorporated into the equation, the multipath effect is absorbed into the ambiguities and position coordinates.The IF linear combination is formed to remove the ionosphere terms.The Earth tidal effect and relativistic effect also need to be considered for the long-distance estimation between ground stations and satellites.
Fixing the phase ambiguities to integer values will increase the accuracy of the satellite orbit determination solution [44].The long baseline processing represents a challenging situation for integer ambiguity resolution.We first used the traditional Kalman filtering algorithm to compute the ambiguity parameters, and the results showed that the fixing rate is very low.It may be caused by the modeling error and divergence of values in the filtering process.Therefore, the ambiguities must be determined epoch-by-epoch for kinematic positioning.The modified least-squares ambiguity decorrelation adjustment (LAMBDA) method has been used for single epoch integer ambiguities least-squares estimation [45].Suppose the float ambiguity solutions vector N in Equations ( 8) and ( 9) have been solved, Q N is its variance-covariance matrix.The direct purpose is to obtain the integer parameter vector N, make the least squares residual minimize, The LAMBDA algorithm changes the ambiguity parameter and its variance-covariance matrix from the original space to a new space by Gauss transformation (so-called ztransformation) and aims to realize the reduction process of ambiguity, Q N and Q N are symmetric positive definite.For the convenience of calculation, Cholesky factorization technique is usually used with them, where L and L are lower triangular matrices, and D and D are diagonal matrices.
There are two main steps in reduction process: first, to make the absolute values of the off-diagonal elements in the L matrix less than or equal to 0.5, and second, to sort the diagonal elements of the D matrix in descending order.The detailed methods are presented in [45].After reducing the correlation between ambiguity parameters, the efficiency of the ambiguity search can be greatly accelerated.By resolving the float IF combination ambiguities, the positioning solution can subsequently be realized as long as the number of synchronous observation satellites is more than or equal to 4.

Kinematic Orbit Determination Strategy
For Swarm B, 50 IGS ground stations were selected based on station performance and geographic distribution to generate the double-difference model.The coordinates of ground stations calculated by CODE and corresponding tropospheric delays are introduced.Nevertheless, for the low-pair Swarm satellites, we take a different approach by forming a baseline of variable length to compute their kinematic position coordinates with Equation ( 9).These baselines were constrained to a maximum cross-track spacing of around 160 km at the equator.The distances from Swarm B to the ground stations are more than 500 km.These two types of double-difference methods represent a complex situation for integer ambiguity fixing and orbit determination due to the potential ionospheric and tropospheric delays in GPS observations.The Melbourne-Wübbena (MW) linear combination algorithm is used to determine carrier phase integer ambiguities.The algorithm works by first solving the wide-lane ambiguities of the MW combination.The wide-lane (WL) ambiguities are the difference in the number of whole cycles between the transmitted and received signals.Once the WL ambiguities are resolved, the narrow-lane (NL) ambiguities are solved by substituting the WL ambiguities into the equation of IF combination phase observations.In the process of ambiguity determination, the ambiguity should be fixed as an integer as possible.This means that the ambiguity should be resolved in such a way that it is as close to a whole number as possible.The resolved integer ambiguities are then employed as known parameters in kinematic orbit determination formulas using IF observation models.This allows for the computation of kinematic position coordinates using only carrier phase observations with fixed integer ambiguities.This significantly improves the accuracy of the positioning information obtained from GPS technology.
By utilizing the data processing algorithms mentioned in Section 2.3, the first-order ionospheric effect which accounts for over 99% of the whole ionospheric influence [46] is eliminated through the IF linear combination of GPS code and carrier observations.Receiver coordinates and phase ambiguities are then estimated to determine the doubledifference parameters.The flow chart of orbit determination is shown in Figure 3.In this article, we utilize the final GPS orbits and high-rate 5-s GPS satellite clock corrections provided by the IGS analysis center.Additionally, we account for the absolute GPS transmit antenna phase center bias and variation in accordance with the igs14.atxfile.The attitude data of the Swarm constellation need to be transformed from the satellite frame to the Earth center Earth fixed coordinates system.Throughout the mission lifetime, the satellite center of mass will experience slight movement with fuel consumption.The actual values of the center of mass position are also contained in the ESA PSO products.Finally, the elevation cut-off angle is set to zero for making full use of all observations.Table 1 showed the POD strategy, which details the measurement model and reference frame used for the Swarm kinematic orbit determination.

Results
Based on the processing standards described above, the time series of three kinds of kinematic orbits for the three Swarm satellites have been generated from 1 January to 31 December 2021: 1.
double-difference kinematic POD with ambiguity resolution.
The zero-difference orbit that meets the accuracy requirements is mainly used as an a priori orbit for a double-difference baseline solution.In theory, the accuracy of fixed ambiguity solutions is higher than that of float ambiguity solutions.The double-difference orbits with ambiguity resolution are the final results.As part of this section, the accuracy of the calculations is evaluated and contrasted with precise science orbits from ESA and satellite laser ranging (SLR) measurements.According to ESA documents, the kinematic and reduced dynamics orbit determination are carried out using the Bernese software.The Bernese software performs orbit determination operations based on zero-difference ionosphere-free observations.Further evidence of the quality of kinematic orbits was obtained by analyzing phase residuals as well.The contrast with reduced-dynamic orbits and independent SLR measurements is an external assessment of the orbits.The residuals of GPS observations provide another internal validation of the calculations.

Ambiguity Resolutions
In order to complete the task of improving the success rate of fixing ambiguities, we adopted a method of selecting an optimal subset of ambiguities and attempted to fix ambiguities multiple times within an epoch.Specifically, when the radio test verification has failed, we removed the satellite with the smallest weight (weighted according to elevation angle) and attempted to fix the ambiguity again.If the ambiguity was still not resolved, we continued removing the satellite with the lowest weight and repeated the process until there were no more satellites available or the ambiguity was successfully resolved.This section analyzes the ambiguity-fixing effect of Swarm satellites from the fixing rate and the changes in phase residuals before and after the fixed ambiguity.Figure 4 shows the results of the fixing ambiguity rate in 2021, which is the proportion of successfully fixed double difference ambiguities to all double difference ambiguities.Through applying this method, the WL ambiguity fixing rates reached 99% for the majority of days in 2021.The average fixing rates of NL ambiguity are 94-95%.Figure 5 demonstrates the mean value of the daily phase residuals before and after the fixed ambiguity.After fixing the ambiguity, the phase residuals of all three satellites increased significantly.The float ambiguities can absorb some other unmodeled errors, whereas the fixed ambiguities make the ambiguity parameters constrained and unable to absorb other errors.Therefore, the phase residuals of the float solutions are slightly larger than that of the fixed solutions.can absorb some other unmodeled errors, whereas the fixed ambiguities make the ambiguity parameters constrained and unable to absorb other errors.Therefore, the phase residuals of the float solutions are slightly larger than that of the fixed solutions.can absorb some other unmodeled errors, whereas the fixed ambiguities make the ambiguity parameters constrained and unable to absorb other errors.Therefore, the phase residuals of the float solutions are slightly larger than that of the fixed solutions.

Carrier Phase Residuals
The carrier phase residual of GPS observations is the difference between the predicted value and the observed value.It can be used to measure the accuracy of the observations and evaluate the internal consistency of the kinematic orbit solutions.The phase residuals were similar to the noise level of the observation when the observations are of satisfactory quality and pre-processed properly [22].The mean values of daily RMS error of the phase residuals for three kinds of kinematic orbits are shown in Table 2.There are significant differences between the three kinds of orbits.The RMS values of zero-difference phase residuals remain at around 4 mm, which is smaller than that of the double-difference phase residuals.This indicates that the accuracy of the double-difference orbits is higher than that of the zero-difference orbits.The detailed daily RMS values of the IF carrier phase observation residuals for final orbits in 2021 are shown in Figure 6.On 1 November at 10:27 UTC, an anomaly has been registered on board Swarm C. As a result, the carrier phase residuals have been affected (the daily RMS is 20 mm).For the rest, the average of the daily RMS for Swarm A, B, and C are 6.2 mm, 6.0 mm, and 6.4 mm, respectively.As expected, the RMS of three satellite GPS phase residuals are kept at around the 6 mm level and indicate that the observations are accurate.
vations and evaluate the internal consistency of the kinematic orbit solutions.Th residuals were similar to the noise level of the observation when the observation satisfactory quality and pre-processed properly [22].The mean values of daily RM of the phase residuals for three kinds of kinematic orbits are shown in Table 2. T significant differences between the three kinds of orbits.The RMS values of zer ence phase residuals remain at around 4 mm, which is smaller than that of the difference phase residuals.This indicates that the accuracy of the double-differen is higher than that of the zero-difference orbits.The detailed daily RMS values of the IF carrier phase observation residuals orbits in 2021 are shown in Figure 6.On 1 November at 10:27 UTC, an anomaly h registered on board Swarm C. As a result, the carrier phase residuals have been (the daily RMS is 20 mm).For the rest, the average of the daily RMS for Swarm A C are 6.2 mm, 6.0 mm, and 6.4 mm, respectively.As expected, the RMS of three GPS phase residuals are kept at around the 6 mm level and indicate that the obse are accurate.

Differences with Reduced Dynamic Orbit Products
The Faculty of Aerospace Engineering at Delft University of Technology is r ble for Swarm satellite precise science orbit determination.The reduced-dynami calculated by GHOST software, and the orbit is checked by the observation da SLR station on the ground, with a precision of less than 2 cm.The products of determination have been published on the ESA website: http://swarm-diss.eo.esa cessed on 23 March 2023).According to the orbit determination strategy given ab results of three kinds of kinematic orbits are compared with those of ESA-redu namic orbits.The statistical results for daily RMS values are presented in Tabl

Differences with Reduced Dynamic Orbit Products
The Faculty of Aerospace Engineering at Delft University of Technology is responsible for Swarm satellite precise science orbit determination.The reduced-dynamic orbit is calculated by GHOST software, and the orbit is checked by the observation data of the SLR station on the ground, with a precision of less than 2 cm.The products of the orbit determination have been published on the ESA website: http://swarm-diss.eo.esa.int(accessed on 23 March 2023).According to the orbit determination strategy given above, the results of three kinds of kinematic orbits are compared with those of ESA-reduced dynamic orbits.The statistical results for daily RMS values are presented in Table 3.The differences of the final double-difference fixed resolutions in all three directions are significantly reduced.The accuracies of the final orbits have been improved by 0.8-1.1 cm and 0.4-0.6 cm in three directions compared to the ZD and DD-float results.In the 3D direction, the accuracy of Swarm A has been improved by 35.2% and 17.4%, respectively, Swarm B has been improved by 37.5% and 20.7%, respectively, and Swarm C has been improved by 36.3% and 16.8%, respectively.Figures 7-9 show the detailed daily RMS of final position residuals for Swarm satellites with respect to ESA-reduced dynamic orbits in three directions.During 2021, the daily RMS values for Swarm A and Swarm C are prone to larger fluctuations, especially in late December, whereas Swarm B is relatively smoother.For the low pair Swarm A and Swarm C, the mean daily RMS of the difference between our kinematic orbit and reduced-dynamic orbit released by ESA in radial, along-track, and cross-track directions are around 2.9 cm, 2.5 cm, and 2.6 cm, respectively.The average of the daily RMS values for Swarm B is about 2.60 cm, 2.13 cm, and 2.45 cm in three directions.The computed Swarm B kinematic orbits are better than Swarm A and C due to their higher quality of observations.Furthermore, compared with the along-track and cross-track directions, the mean RMS values for the radial direction have lower accuracy.This may be caused by the weak observation geometry.The radial direction is also usually not well determined in kinematic orbit determination because the receiver clock corrections must be estimated at the same time as the kinematic positions [14].On the whole, the mean 3D RMS values of the position residuals between kinematic calculations and reduced-dynamic orbits are below 4.7 cm for Swarm A and C and below 4.2 cm for Swarm B, which can satisfy the requirement for precision (cm level) and have application for gravity field recovery.
direction, the accuracy of Swarm A has been improved by 35.2% and 17.4%, respect Swarm B has been improved by 37.5% and 20.7%, respectively, and Swarm C has improved by 36.3% and 16.8%, respectively.9 show the detailed daily RMS of final position residuals for Swarm lites with respect to ESA-reduced dynamic orbits in three directions.During 202 daily RMS values for Swarm A and Swarm C are prone to larger fluctuations, espe in late December, whereas Swarm B is relatively smoother.For the low pair Swarm A Swarm C, the mean daily RMS of the difference between our kinematic orbit and red dynamic orbit released by ESA in radial, along-track, and cross-track direction around 2.9 cm, 2.5 cm, and 2.6 cm, respectively.The average of the daily RMS valu Swarm B is about 2.60 cm, 2.13 cm, and 2.45 cm in three directions.The computed S B kinematic orbits are better than Swarm A and C due to their higher quality of obs tions.Furthermore, compared with the along-track and cross-track directions, the RMS values for the radial direction have lower accuracy.This may be caused by the observation geometry.The radial direction is also usually not well determined in matic orbit determination because the receiver clock corrections must be estimated same time as the kinematic positions [14].On the whole, the mean 3D RMS values position residuals between kinematic calculations and reduced-dynamic orbits are b 4.7 cm for Swarm A and C and below 4.2 cm for Swarm B, which can satisfy the req ment for precision (cm level) and have application for gravity field recovery.

SLR Residuals
The calculated kinematic orbits are also compared against SLR measurements precision of millimeters to centimeters can be achieved with fully independent op distance observations between the Swarm satellite and tracking station by using the system [47,48].In this way, the differences between measured and modeled values, other words, the SLR residuals, can be used to validate not only the precision but als accuracy of satellite orbits.SLR observations are typically received from the Internat Laser Ranging Service (ILRS) tracking stations.
During the study period, there are 28 ground tracking stations that provided av ble SLR data.The global distribution of these 28 stations is given in Figure 10.The from Wuhan station are all unavailable.The total number of normal points for the Swarm satellites are 28,516, 82,126, and 26,952, respectively.More normal points are Swarm B since it flies in a higher orbit.A moderate separation between Swarm A Swarm C results in common station visibility.The computed kinematic solutions are polynomially interpolated to the observed time of the SLR observations.Next, with sideration of the tropospheric delay, relativistic effects, solid Earth tides, and oc

SLR Residuals
The calculated kinematic orbits are also compared against SLR measurements precision of millimeters to centimeters can be achieved with fully independent o distance observations between the Swarm satellite and tracking station by using the system [47,48].In this way, the differences between measured and modeled values, other words, the SLR residuals, can be used to validate not only the precision but als accuracy of satellite orbits.SLR observations are typically received from the Interna Laser Ranging Service (ILRS) tracking stations.
During the study period, there are 28 ground tracking stations that provided a ble SLR data.The global distribution of these 28 stations is given in Figure 10.The from Wuhan station are all unavailable.The total number of normal points for the Swarm satellites are 28,516, 82,126, and 26,952, respectively.More normal points are Swarm B since it flies in a higher orbit.A moderate separation between Swarm A Swarm C results in common station visibility.The computed kinematic solutions ar polynomially interpolated to the observed time of the SLR observations.Next, with sideration of the tropospheric delay, relativistic effects, solid Earth tides, and oc loads, the distances between the kinematic orbits and the ILRS stations are calculate

SLR Residuals
The calculated kinematic orbits are also compared against SLR measurements.The precision of millimeters to centimeters can be achieved with fully independent optical distance observations between the Swarm satellite and tracking station by using the SLR system [47,48].In this way, the differences between measured and modeled values, or in other words, the SLR residuals, can be used to validate not only the precision but also the accuracy of satellite orbits.SLR observations are typically received from the International Laser Ranging Service (ILRS) tracking stations.
During the study period, there are 28 ground tracking stations that provided available SLR data.The global distribution of these 28 stations is given in Figure 10.The data from Wuhan station are all unavailable.The total number of normal points for the three Swarm satellites are 28,516, 82,126, and 26,952, respectively.More normal points are from Swarm B since it flies in a higher orbit.A moderate separation between Swarm A and Swarm C results in common station visibility.The computed kinematic solutions are first polynomially interpolated to the observed time of the SLR observations.Next, with consideration of the tropospheric delay, relativistic effects, solid Earth tides, and oceanic loads, the distances between the kinematic orbits and the ILRS stations are calculated and compared with the laser ranges.In order to detect outliers that may still exist in the SLR measurements, a screening threshold of 20 cm was used for the residuals.The SLR residuals for the three Swarm satellites during 2021 are shown in Figures 11-13.The mean RMS values of SLR residuals are summarized in Table 4.It is worth mentioning that the residual values solved based on the SLR stations located in China are relatively large.For most SLR observations, the residuals between the kinematic orbits and SLR measurements are under the 5 cm level.The RMS value for Swarm A's kinematic orbits is 1.8 cm, Swarm B is 2.1 cm, and Swarm C is 1.7 cm, respectively.4. It is worth mentioning that the residual values solved based on the SLR stations located in China are relatively large.For most SLR observations, the residuals between the kinematic orbits and SLR measurements are under the 5 cm level.The RMS value for Swarm A's kinematic orbits is 1.8 cm, Swarm B is 2.1 cm, and Swarm C is 1.7 cm, respectively.4. It is worth mentioning that the residual values solved based on the SLR stations located in China are relatively large.For most SLR observations, the residuals between the kinematic orbits and SLR measurements are under the 5 cm level.The RMS value for Swarm A's kinematic orbits is 1.8 cm, Swarm B is 2.1 cm, and Swarm C is 1.7 cm, respectively.4. It is worth mentioning that the residual values solved based on the SLR stations located in China are relatively large.For most SLR observations, the residuals between the kinematic orbits and SLR measurements are under the 5 cm level.The RMS value for Swarm A's kinematic orbits is 1.8 cm, Swarm B is 2.1 cm, and Swarm C is 1.7 cm, respectively.

Discussion
For the assessment of the proposed double-difference methods, the three kinds of Swarm kinematic orbits were listed from 1 January to 31 December 2021.The zero-difference model is efficient in solving the problem but requires an accurate error correction model.On the contrary, the double-difference model can effectively weaken or eliminate the same kind of error sources, such as the GPS satellite clock offset and receiver clock offset.expected, the accuracy of the final double-difference model with ambiguity resolution is higher than the zero-difference model and float resolution.However, it is inevitable that solving double-difference observations requires a significant amount of computational time and storage space.Analysis of the final ionosphere-free carrier phase residuals suggests that the mean daily RMS values for all three Swarm satellites are at around the 6 mm level.The carrier phase residuals under different solar cycle conditions remained below 9 mm and stayed at a satisfactory level, showing strong internal consistency.A comparison with the reduced-dynamic orbits released by ESA shows the mean daily RMS of the orbit differences in radial, along-track, and cross-track directions are around 2.9 cm, 2.5 cm, and 2.6 cm, respectively, for Swarm A and Swarm C. The average of the daily RMS values are about 2.60 cm, 2.13 cm, and 2.45 cm in three directions for Swarm B. For most SLR observations, the residuals between the kinematic orbits and SLR measurements are under the 5 cm level.The RMS values of SLR residuals for the three Swarm satellites' kinematic orbits are 1.8 cm, 2.1 cm, and 1.7 cm, respectively, showing a very good agreement with the ESA level.For Swarm B, an improved orbit condition is apparent.It is also worth mentioning that the daily RMS of the position residuals for the low pair satellites has comparatively great fluctuation than Swarm B.
Solar activity has a significant impact on signal communication and satellite navigation.Strong solar activity increases the concentration of electrons and ions in the ionosphere, which increases thermospheric density and atmospheric drag, thus affecting satellite operations.The solar radiation index F10.7 is often used to indicate the intensity of solar activity.F10.7 is the radio radiation flow of the sun at 10.7 cm wavelength (2800 MHz), and its unit is sfu (1 sfu = 1 × 10 −22 W•m −2 •Hz −1 ).The value of the solar radiation index F10.7

Discussion
For the assessment of the proposed double-difference methods, the three kinds of Swarm kinematic orbits were listed from 1 January to 31 December 2021.The zerodifference model is efficient in solving the problem but requires an accurate error correction model.On the contrary, the double-difference model can effectively weaken or eliminate the same kind of error sources, such as the GPS satellite clock offset and receiver clock offset.As expected, the accuracy of the final double-difference model with ambiguity resolution is higher than the zero-difference model and float resolution.However, it is inevitable that solving double-difference observations requires a significant amount of computational time and storage space.Analysis of the final ionosphere-free carrier phase residuals suggests that the mean daily RMS values for all three Swarm satellites are at around the 6 mm level.The carrier phase residuals under different solar cycle conditions remained below 9 mm and stayed at a satisfactory level, showing strong internal consistency.A comparison with the reduced-dynamic orbits released by ESA shows the mean daily RMS of the orbit differences in radial, along-track, and cross-track directions are around 2.9 cm, 2.5 cm, and 2.6 cm, respectively, for Swarm A and Swarm C. The average of the daily RMS values are about 2.60 cm, 2.13 cm, and 2.45 cm in three directions for Swarm B. For most SLR observations, the residuals between the kinematic orbits and SLR measurements are under the 5 cm level.The RMS values of SLR residuals for the three Swarm satellites' kinematic orbits are 1.8 cm, 2.1 cm, and 1.7 cm, respectively, showing a very good agreement with the ESA level.For Swarm B, an improved orbit condition is apparent.It is also worth mentioning that the daily RMS of the position residuals for the low pair satellites has comparatively great fluctuation than Swarm B.
Solar activity has a significant impact on signal communication and satellite navigation.Strong solar activity increases the concentration of electrons and ions in the ionosphere, which increases thermospheric density and atmospheric drag, thus affecting satellite operations.The solar radiation index F 10.7 is often used to indicate the intensity of solar activity.F 10.7 is the radio radiation flow of the sun at 10.7 cm wavelength (2800 MHz), and its unit is sfu (1 sfu = 1 × 10 −22 W•m −2 •Hz −1 ).The value of the solar radiation index F 10.7 can properly reflect the level of ionospheric conditions.As the value of F 10.7 increases, the ionospheric scintillations which are disturbed by solar activity will be more intense [21].In order to apprehend whether the kinematic orbits are impacted by ionospheric activities and further verify the effectiveness of the kinematic orbit determination strategy proposed above, we did an analysis of the F 10.7 index.The data of F 10.7 used in this article come from the National Centers for Environmental Information (NCEI) of the National Oceanic and Atmospheric Administration (NOAA): https://www.ngdc.noaa.gov(accessed on 23 March 2023).The average daily F 10.7 index in 2019 is 69.7 sfu, in 2020, it was 73.8 sfu, and in 2021 it was 81.8 sfu. Figure 14 represents the sequence diagram of the solar radiation index F 10.7 in 2021, and the data show that the overall level of solar activity is gradually rising.We can find that index F 10.7 fluctuated between 110 sfu and 150 sfu from 16 December to 29 December (DoY 350-363), showing that the ionospheric activity during this time is relatively intense.This may explain why the daily RMS values for the low pair satellites are highly volatile in late December.can properly reflect the level of ionospheric conditions.As the value of F10.7 increases, the ionospheric scintillations which are disturbed by solar activity will be more intense [21].
In order to apprehend whether the kinematic orbits are impacted by ionospheric activities and further verify the effectiveness of the kinematic orbit determination strategy proposed above, we did an analysis of the F10.  Figure 14 also shows the difference between the final double-difference orbit and the ESA-reduced dynamic orbit in 3D.We discussed the statistical method of the position difference to demonstrate the impact of the F10.7 index on the daily RMS variations.One of the statistical indicators of the strength of a monotonic relationship between paired data is Spearman's correlation coefficient rs [49].Spearman correlation analysis describes the correlation of two variables by a linear equation and does not require a priori information of the original data.It is a non-parametric statistical program that has a wide range of applications in mathematical statistics.The correlation coefficient rs and the p-value are computed and entered in Table 5.The significance test results are referred to as the pvalue, if the p-value is small enough, we can be more certain that the observed correlation is not due to chance or random variation.A smaller p-value suggests a higher level of confidence that there is a correlation between kinematic orbits and ionospheric scintillations.The bigger absolute values of rs reflects the more relevant the correlation.As shown in Table 5, the rs of three directions range from 0.2 to 0.3 for the low pair satellites, which indicates the F10.7 index is in weak relativity to orbit residuals.However, the rs values for Swarm B are approximately equal to zero and show there is almost no correlation between the radio radiation flow of the sun and Swarm B kinematic orbit residuals.Figure 14 also shows the difference between the final double-difference orbit and the ESA-reduced dynamic orbit in 3D.We discussed the statistical method of the position difference to demonstrate the impact of the F 10.7 index on the daily RMS variations.One of the statistical indicators of the strength of a monotonic relationship between paired data is Spearman's correlation coefficient r s [49].Spearman correlation analysis describes the correlation of two variables by a linear equation and does not require a priori information of the original data.It is a non-parametric statistical program that has a wide range of applications in mathematical statistics.The correlation coefficient r s and the p-value are computed and entered in Table 5.The significance test results are referred to as the p-value, if the p-value is small enough, we can be more certain that the observed correlation is not due to chance or random variation.A smaller p-value suggests a higher level of confidence that there is a correlation between kinematic orbits and ionospheric scintillations.The bigger absolute values of r s reflects the more relevant the correlation.As shown in Table 5, the r s of three directions range from 0.2 to 0.3 for the low pair satellites, which indicates the F 10.7 index is in weak relativity to orbit residuals.However, the r s values for Swarm B are approximately equal to zero and show there is almost no correlation between the radio radiation flow of the sun and Swarm B kinematic orbit residuals.Furthermore, we calculated the position residuals between the obtained kinematic orbits and the PSO products provided by ESA in intense ionospheric activity (sfu ≥ 110) and in relatively peaceful periods (sfu < 110).The results are shown in Table 6.In the high ionosphere active period, the mean daily RMS values of the difference between two orbit types in the radial direction are 3.38 cm, 2.61 cm, and 3.39 cm for three Swarm satellites, and 2.94 cm, 2.17 cm, and 2.90 cm, respectively, in the along-track direction.In the cross-track direction, the mean daily RMS are 3.09 cm, 2.47 cm, and 2.98 cm, respectively.The mean values of daily RMS in three directions are all larger than in a peaceful period, whereas the changes for Swarm B were not significant.From the above studies, it is worth mentioning that the daily RMS values of the difference with the PSO products change little for Swarm B, with comparatively great fluctuation for Swarm A and Swarm C, especially in late December.This can be interpreted that Swarm B has a higher flight orbit, and the effect of the ionosphere is weaker than the low pair satellites.

Conclusions
In this study, we discussed approach for determining the kinematic precise orbits with two different double-difference methods, one applied to the low pair satellites and the other applied to Swarm B. Swarm kinematic orbits are generated with moving base and ground stations, respectively.The accuracy and reliability of kinematic orbit rely on the quality of GPS data; therefore, an extensive analysis of GPS data was first performed.Although only eight channels are available for the Swarm constellation, the average percent of tracking eight GPS satellites for three Swarm satellites are 83.52%,82.38%, and 83.67%, respectively.After processing all the GPS observations, we found that the ambiguities could not be fixed effectively by the Kalman filtering approach.The final WL and NL ambiguity fixing rates reached over 99% and 94% for the Swarm constellation by adopting a program of selecting an optimal subset of ambiguities and attempting to fix ambiguities multiple times within an epoch.
After evaluating the orbit quality, the overall results of the Swarm satellites' kinematic orbits determination accord with the requirement precision (cm level).Compared to the zero-difference orbits, the accuracy of Swarm satellites was improved by 35.2%, 37.5%, and 36.3%,For the double-difference orbits without ambiguity resolutions, the accuracy was improved by 17.4%, 20.7%, and 16.8%, respectively.The daily RMS values of the IF phase residuals for Swarm satellites were maintained at about the 6 mm level.In contrast, the RMS values of the position residuals with the reduced-dynamic orbits offered by ESA are at about the 2-3 cm level in three directions, respectively.Compared to the SLR measurements, the RMS values of the SLR residuals for kinematic orbits are about 2 cm.
Ionospheric scintillations occur when ionospheric activity becomes intense.GPS receiver operations can be interrupted or degraded, contributing to the position difference between kinematic orbits and reduced-dynamic orbits showing a larger fluctuation.The changes in ionospheric activity are rapid in late December, and for Swarm A and Swarm C, the orbit residuals are strongly perturbed.Swarm B has a higher flight orbit than the low-pair satellites, so it is less affected by the influence of the ionosphere.In addition, ithere were weak monotonic relationships between the F 10.7 index and the low pair satellites' kinematic orbit residuals that were identified through statistical analysis.In a peaceful period, the mean daily RMS of positions residuals between the obtained kinematic orbits and ESA reduced-dynamic orbits in the along-track direction is 2.86 cm and 2.91 cm for Swarm A and Swarm C. Thecross-track directions are 2.45 cm and 2.41 cm, and the radial directions are 2.62 cm and 2.61 cm, respectively.However, in periods with strong ionosphere solicitations, the mean daily RMS are 3.38 cm, 2.94 cm, and 3.09 cm for Swarm and 3.39 cm, 2.90 cm, and 2.98 cm for Swarm C in three directions.Notwithstanding, ionospheric scintillations do not disturb the carrier phase residuals much.
As expected, using both double-difference methods under different solar cycle conditions satisfies the requirement of centimeter accuracy in the three Swarm satellite orbit solutions and has application value for gravity field determination.Although this article has made a systematic study of Swarm satellite observation data quality and double-difference kinematic orbit determination, some useful conclusions have been obtained.However, there are still some shortcomings that need further research and improvement.For example, in phase residual analysis, it is also necessary to study the inherent effects of phase residual change, multipath error, and observation noise.For zero-difference kinematics orbit determination, the performance of the error processing model needs to be improved.Consider adding dynamic models to achieve reduced dynamic orbit determination.

Figure 1 .
Figure 1.Swarm GPS data tracking observation of the year 2021.Figure 1. Swarm GPS data tracking observation of the year 2021.

Figure 1 .
Figure 1.Swarm GPS data tracking observation of the year 2021.Figure 1. Swarm GPS data tracking observation of the year 2021.

Figure 2 .
Figure 2. (a) Configurations of double-difference methods for the ground station; (b) Configurations of double-difference methods for satellite-to-satellite.

Figure 2 .
Figure 2. (a) Configurations of double-difference methods for the ground station; (b) Configurations of double-difference methods for satellite-to-satellite.

Figure 3 .
Figure 3.The flow chart of orbit determination.

Figure 3 .
Figure 3.The flow chart of orbit determination.

Figure 4 .
Figure 4.The time series of WL and NL ambiguities fixing rates of the Swarm satellites.

Figure 5 .
Figure 5.The mean daily RMS of the phase residuals of float and fixed ambiguity results.

Figure 4 .
Figure 4.The time series of WL and NL ambiguities fixing rates of the Swarm satellites.

Figure 4 .
Figure 4.The time series of WL and NL ambiguities fixing rates of the Swarm satellites.

Figure 5 .
Figure 5.The mean daily RMS of the phase residuals of float and fixed ambiguity results.

Figure 5 .
Figure 5.The mean daily RMS of the phase residuals of float and fixed ambiguity results.

Figure 6 .
Figure 6.Daily RMS error of the carrier phase residuals of kinematic orbit determination fo Swarm satellites during 2021.

Figure 6 .
Figure 6.Daily RMS error of the carrier phase residuals of kinematic orbit determination for all three Swarm satellites during 2021.

Figure 7 .
Figure 7.The daily RMS of positions residuals for Swarm A with respect to ESA reduced-dynamic orbits in the along-track, cross-track, and radial component for 2021.

Figure 7 .
Figure 7.The daily RMS of positions residuals for Swarm A with respect to ESA reduced-dyn orbits in the along-track, cross-track, and radial component for 2021.

Figure 8 .
Figure 8.The daily RMS of positions residuals for Swarm B with respect to ESA reduced-dyn orbits in the along-track, cross-track, and radial component for 2021.

Figure 9 .
Figure 9.The daily RMS of positions residuals for Swarm C with respect to ESA reduced-dyn orbits in the along-track, cross-track, and radial component for 2021.

Figure 8 .
Figure 8.The daily RMS of positions residuals for Swarm B with respect to ESA reduced-dynamic orbits in the along-track, cross-track, and radial component for 2021.

Figure 7 .
Figure 7.The daily RMS of positions residuals for Swarm A with respect to ESA reduced-dy orbits in the along-track, cross-track, and radial component for 2021.

Figure 8 .
Figure 8.The daily RMS of positions residuals for Swarm B with respect to ESA reduced-dy orbits in the along-track, cross-track, and radial component for 2021.

Figure 9 .
Figure 9.The daily RMS of positions residuals for Swarm C with respect to ESA reduced-dy orbits in the along-track, cross-track, and radial component for 2021.

Figure 9 .
Figure 9.The daily RMS of positions residuals for Swarm C with respect to ESA reduced-dynamic orbits in the along-track, cross-track, and radial component for 2021.
Remote Sens. 2023, 15, x FOR PEER REVIEW 16 of 22 measurements, a screening threshold of 20 cm was used for the residuals.The SLR residuals for the three Swarm satellites during 2021 are shown in Figures 11-13.The mean RMS values of SLR residuals are summarized in Table

Figure 10 .
Figure 10.The distribution of the SLR tracking stations on the ground.

Figure 11 .
Figure 11.The SLR residuals for the Swarm A kinematic orbits.

Figure 12 .
Figure 12.The SLR residuals for the Swarm B kinematic orbits.

Figure 10 .
Figure 10.The distribution of the SLR tracking stations on the ground.

Figure 10 .
Figure 10.The distribution of the SLR tracking stations on the ground.

Figure 11 .
Figure 11.The SLR residuals for the Swarm A kinematic orbits.

Figure 12 .
Figure 12.The SLR residuals for the Swarm B kinematic orbits.

Figure 11 .
Figure 11.The SLR residuals for the Swarm A kinematic orbits.

Figure 10 .
Figure 10.The distribution of the SLR tracking stations on the ground.

Figure 11 .
Figure 11.The SLR residuals for the Swarm A kinematic orbits.

Figure 12 .
Figure 12.The SLR residuals for the Swarm B kinematic orbits.

Figure 12 .
Figure 12.The SLR residuals for the Swarm B kinematic orbits.

Figure 13 .
Figure 13.The SLR residuals for the Swarm C kinematic orbits.

Figure 13 .
Figure 13.The SLR residuals for the Swarm C kinematic orbits.
7 index.The data of F10.7 used in this article come from the National Centers for Environmental Information (NCEI) of the National Oceanic and Atmospheric Administration (NOAA): https://www.ngdc.noaa.gov(accessed on 23 March 2023).The average daily F10.7 index in 2019 is 69.7 sfu, in 2020, it was 73.8 sfu, and in 2021 it was 81.8 sfu.Figure14represents the sequence diagram of the solar radiation index F10.7 in 2021, and the data show that the overall level of solar activity is gradually rising.We can find that index F10.7 fluctuated between 110 sfu and 150 sfu from 16 December to 29 December (DoY 350-363), showing that the ionospheric activity during this time is relatively intense.This may explain why the daily RMS values for the low pair satellites are highly volatile in late December.

Figure 14 .
Figure 14.Solar radiation index F10.7 (black) and 3D orbit difference RMS (blue for Swarm A, red for Swarm B, and green for Swarm C) of year 2021.

Figure 14 .
Figure 14.Solar radiation index F 10.7 (black) and 3D orbit difference RMS (blue for Swarm A, red for Swarm B, and green for Swarm C) of year 2021.

2 are the pseudo range and carrier phase observation, respectively. Since the t r and t p are estimated by IF combination observations, the code hardware delay of the GPS satellite and Swarm receiver in Equation (
P,IF and ε L,IF are the sum of IF random noise and multipath effects; f 1 and f 2 are the carrier frequencies; P 1 /P 2 and L 1 /L

Table 1 .
Summary of the measurement model and reference frame used for the Swarm kinema orbit determination.

Table 1 .
Summary of the measurement model and reference frame used for the Swarm kinematic orbit determination.

Table 2 .
Mean of daily RMS error of the carrier phase residuals for all three Swarm satellites during 2021.

Table 2 .
Mean of daily RMS error of the carrier phase residuals for all three Swarm satellit 2021.

Table 3 .
The mean daily RMS of the difference between the kinematic orbit and ESA reduced dynamic orbit during 2021.

Table 3 .
The mean daily RMS of the difference between the kinematic orbit and ESA reduce namic orbit during 2021.

Table 4 .
The mean and RMS values of SLR residuals for the Swarm satellites using observations from SLR stations.

Table 4 .
The mean and RMS values of SLR residuals for the Swarm satellites using observations from SLR stations.

Table 5 .
Spearman's correlation between the F10.7 index and orbit difference for the Swarm satellites.

Table 5 .
Spearman's correlation between the F 10.7 index and orbit difference for the Swarm satellites.

Table 6 .
The mean daily RMS of the difference between the kinematic orbit and reduced-dynamic orbit under different solar cycle conditions.