1. Introduction
Satellites can fly in proximity and in formation as a team in order to outperform traditional single satellite missions. Such a strategy has motivations from the animal kingdom in which, for example, bird species form V-shapes to improve the aerodynamic efficiency of flight, communicate, and coordinate better within the flock, as well as more easily sense their prey [
1]. Several benefits can be expected from satellite formation flying (SFF), including scalability in terms of the fleet configuration, multidimensional flexibility for scientific observations, lower cost budget by launching a group of small satellites, and sustainability of mission operation. The mission architecture of SFF has been applied to many space science missions, such as the Magnetospheric Multiscale (MMS) mission [
2], the Gravity Recovery and Climate Experiment (GRACE) mission [
3], and the European Space Agency (ESA) Swarm mission [
4].
The performance of guidance, navigation, and control (GNC) systems is very critical in SFF missions. As the most ubiquitous modern positioning, navigation, and timing (PNT) technology, the global navigation satellite system (GNSS) is commonly used in absolute navigation and relative navigation for SFF. The navigation accuracy and reliability continuously increase rapidly using multiple global GNSS systems including the American Global Positioning System (GPS), the Russian GLONASS constellation, the European Union’s Galileo constellation and the Chinese BeiDou (or COMPASS) constellation. The modernization of these GNSS constellations plus the development of newer generation GNSS receivers with more advanced software algorithms have brought GNSS users into a golden era of multi-constellation GNSS with unprecedented quality of PNT services. Precise relative navigation can be accomplished by using differential GNSS, where centimeter or sub-centimeter level accuracy of short baseline (∼1 km) relative position determination in low Earth orbit (LEO) can be achieved by utilizing the single or double differential carrier phase technique [
5,
6].
Beside navigation, GNSS is also widely used for ionospheric remote sensing as the Earth’s ionosphere dynamically impacts GNSS signal propagation. Multi-frequency GNSS receivers are widely used for ground-based ionospheric remote sensing, such as total electron content (TEC) measurements and GNSS scintillation observations. MIT Madrigal database gathers thousands of GPS stations to create a global TEC map database [
7], which greatly benefits the community for ionospheric space weather monitoring and studies. Space-compliant GNSS receivers are also utilized to conduct space-based GNSS sounding projects, such as the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) mission [
8]. The key technique implemented by the COSMIC mission is radio occultation, where the bending effect of the GNSS signals propagating through the ionosphere is utilized to retrieve the electron density (
Ne) of the ionosphere, measure GNSS scintillation, and produce other atmospheric sounding data. Some recent CubeSat missions also apply this radio occultation technique to measure ionospheric scintillation and sense ionospheric irregularities as well, such as the Compact Total Electron Content Sensor (CTECS) [
9] and the Scintillation Prediction Observations Research Task (SPORT) [
10]. The Coherent Electromagnetic Radio Tomography (CERTO) constellation of radio beacons were developed to fly on LEO satellites to measure TEC, scintillations, and plasma irregularities below the satellite orbit together with ground-based beacon receivers [
11]. A pioneering project, called the Ionospheric Observation Nanosatellite Formation (ION-F), proposed and discussed the potential application of spacecraft formation flying to ionospheric measurement [
12]. GNSS-based LEO SFF opens new doors for ionospheric remote sensing.
For satellite missions, a simulated platform to test the functionality of all the key hardware and software systems is required to validate the mission feasibility before launch. Hardware-in-the-loop (HIL) simulation testbeds are an ideal platform during research and development phase to prototype GNSS algorithms and assess GNSS receiver(s) performance in various emulated scenarios of GNSS-based SFF. The first of this kind of simulation testbed was established by United States’ National Aeronautics and Space Administration (NASA) at 2001 [
13], and later on being used to support several formation flying missions (e.g., MMS) development. The German Aerospace Center (DRL) [
14] and a few other universities (e.g., Stanford University [
15], Massachusetts Institute of Technology [
5], University of Toronto [
16], and Yonsei University [
17]) have also been involved in this type of GNSS-based SFF simulation testbed development. However, most of these simulation testbeds are primarily applied to GNSS algorithm development rather than ionospheric remote sensing missions.
A GPS-based simulation testbed for two-satellite formation flying, the Virginia Tech Formation Flying Testbed (VTFFTB), was recently developed with applications to ionospheric remote sensing [
18]. An ionospheric plasma bubble observation scenario was incubated and verified on the VTFFTB by running HIL simulations. Vertical
Ne profiles can be retrieved with differential GPS vertical TEC, and simulated ionospheric irregularities can be investigated by analyzing
Ne, TEC, and amplitude scintillation measurements. The purpose of this study is to describe the extension of the VTFFTB into a multi-constellation simulation testbed by adding Galileo. Various benefits come with the addition of Galileo, which includes, but is not limited to, more signal and data types, better spatial and temporal coverage, higher navigation, and remote sensing accuracy. This is advantageous to multi-scale observation and multi-constellation GNSS comparisons. Compared to a GPS-only scenario, applying multi-constellation GNSS will enhance the capability of both ionospheric remote sensing and navigation with the additional Galileo satellites and signals. Twenty-two (22) Galileo satellites were operational as of March 2019 with more to be launched. The historical improvement on ionospheric remote sensing due to the growing Galileo constellation can be quantitatively evaluated by simulating the Galileo almanac at different times.
The hardware and software infrastructure of the new VTFFTB are overviewed in
Section 2. The ionospheric measurement and navigation results from HIL simulations are presented, analyzed, and discussed in
Section 3. Finally, the conclusions, current issues and future work are summarized in
Section 4.
2. Methods and Infrastructure
An overview of the experimental infrastructure and methods is presented here. The hardware configuration (as shown in
Figure 1) and software algorithms are developed based on the previously established VTFFTB [
18].
The previous version of VTFFTB includes a GPS radio frequency (RF) signal simulator (Spirent GSS-8000) system, a multi-constellation multi-frequency GNSS receiver (NovAtel OEM628), a GPS-based navigation and control system, a GPS-based ionospheric remote sensing system, and an STK visualization system. A GPS-based HIL simulation testbed for two-spacecraft formation flight was considered. The GSS-8000 GPS simulator is capable of simulating GPS (L1, L2, and L5) RF signals and emulating ionospheric impacts (e.g., TEC and amplitude scintillation) on GPS signals. In a satellite scenario, the GPS signals received by the satellite antenna are output from the RF simulator and then fed to the OEM628 GNSS receiver through a coaxial cable.
The GNSS measurement module in the navigation and control system extracts real time data from the GNSS receiver via an Universal Serial Bus (USB) communication interface. A single-differential carrier phase measurement model [
5] is implemented with an extended Kalman filter (EKF) in order to correct the relative states predicted by a dynamic propagator when the measurement and estimation are available. The relative position and velocity between the chief and deputy satellites serve as the input of the controller to compute the required thrust to maneuver the deputy spacecraft given the desired relative orbit. The controller implements the state dependent Riccati equation technique based on the Hill–Clohessy–Wiltshire (HCW) relative motion model [
19]. Finally a remote control module generates and transfers the motion command to the GPS simulator by TCP/IP to propagate the satellite orbit in the simulation. These iterated tasks form a closed-loop real-time feedback system with a default looping rate of 1 Hz. Note that, only one GNSS receiver (i.e., NovAtel OEM628) was available to track the LEO satellite, which is similar to the circumstance in [
13]. Therefore, the chief satellite scenario was simulated first without any active control, and then the deputy satellite was simulated after loading the previously recorded chief data.
The ionospheric remote sensing system processes the GNSS data from the simulation to generate TEC, effective amplitude scintillation index (S4), and vertical
Ne. The TEC processing algorithm is summarized as a flowchart in
Figure 2. Pseudorange, carrier-phase, and ephemeris are extracted from the observation and navigation data collected from the OEM628 receiver. After loading the LEO receiver trajectory, pseudorange, and ephemeris data, the raw pseudoranage relative slant TECs are computed given the constellation and frequency combination selections. The raw carrier-phase relative slant TECs are then computed using the LEO receiver trajectory, carrier-phase, and ephemeris. Concurrently, the elevation angles at each time step are calculated as well. Next, the raw pseudoranage relative slant TECs are utilized to level or fit the raw carrier-phase relative slant TECs and produce the fitted relative carrier-phase TEC. A differential code bias (DCB) exists in the fitted relative slant TEC for each PRN. Therefore, a differential linear least-squares method [
20] is implemented to estimate the DCB and correct the slant TEC bias for each PRN. In this way, the bias-free slant TECs are obtained. Finally, the vertical TEC (VTEC) are generated using the elevation data, and used to retrieve vertical
Ne.
The effective S4 is calculated by taking the ratio of the standard deviation of
to the mean of
over a 1-minute period as:
This algorithm is different from the standard routines to compute
, where signal intensity
I is used in Equation (
1) instead of
[
21]. However, a benchmarking indicated the effective
values computed by Equation (
1) are highly correlated with the
generated from the NovAtel GPStation-6 receiver, a commercial off-the-shell (COTS) ionospheric TEC and scintillation monitor. The STK visualization system can visualize the satellite trajectory in either real-time or replay mode. More details of these systems of VTFFTB have been introduced in [
18].
The new version of the VTFFTB (or VTFFTB 2.0) is an extension of the VTFFTB by adding the Galileo relevant components, including a Spirent GSS-7800 Galileo signal generator, an upgraded version of navigation and control system, and a multi-constellation version of ionospheric remote sensing system. The GSS-7800 signal generator is capable of simulating Galileo (E1, E5a, E5b, E5 AltBOC) RF signals. The GSS-8000 GPS signal generator and GSS-7800 Galileo signal generator have been synchronized together to simulate multi-constellation (GPS + Galileo) scenarios. The new navigation and control system is designed to handle both GPS and Galileo data from the OEM628 receiver. The previous ionospheric remote sensing is expanded to process and generate Galileo TEC, Ne, and effective S4.
3. Results and Discussion
3.1. Scenario Overview
The Equatorial Spread F (ESF) scenario simulated in [
18] is chosen as the baseline simulation for multi-constellation GNSS comparisons. ESF is a type of ionospheric irregularity that often occurs in the post-sunset equatorial region of the ionosphere and it is typically seen as plasma bubbles or plumes within unstable plasma depletions. Studies show that ESF is generated by plasma instabilities; primarily the (generalized) Rayleigh–Taylor instability along with other secondary instabilities that impact the development of the plasma bubble/plume [
22]. ESF can negatively affect communication signals and disturb satellite operation. In terms of the impacts on GNSS, it causes signal scintillation (i.e., rapid temporal fluctuations in amplitude and phase), signal delay due to TEC gradient, and eventually degradation in the navigation reliability and accuracy. Therefore, it is significant to study ESF to mitigate its potential impacts on satellite communication and operation systems.
According to the phase screen scintillation theory discussed in [
23], amplitude scintillation is mainly caused by plasma density irregularities with the scale size corresponding to the dimension of the first Fresnel zone of the RF signals. The Fresnel length is expressed as
, where
is the signal wavelength and
h is the altitude of the ionospheric irregularities. The ionospheric irregularities (assuming
h = 350 km) mostly contribute to the amplitude scintillation of GPS L1 band (19 cm) is approximately 400 m [
24]. Phase screen theory indicates amplitude scintillation is positively correlated with the plasma density fluctuations, i.e., higher electron density deviation will induce larger S4. Past ESF observations results by [
25] are consistent with these scintillation characteristics produced by equatorial plasma bubbles. Plasma bubbles can occur for hours after sunset. Smaller-scale ionospheric irregularities measured by radar appear to decay faster.
The ionospheric impact on GNSS signals can be customized in the Spirent GNSS simulator by modeling TEC and amplitude scintillation. By default, TEC variation can be modeled vertically only, and the S4 can only by modeled in a region lower than 350 km defined by a horizontal grid with the minimum resolution of 10 by 15 degrees. A cuboid region with a few plasma bubbles was simulated in the VTFFTB by setting the one-dimensional vertical TEC profile and two-dimensional horizontal S4 grid. Using a
Ne measurement from the PLUMEX I sounding rocket (SR) campaign [
26], a vertical
Ne profile including three obvious plasma bubbles is used to derive the vertical TEC profile. Horizontally, the S4 was set as 0.4 in the cuboid region ranges from 290 to 350 km vertically, 20:00 to 21:00 local time (LT) longitudinally, and 0° to 10°S latitudinally, which is above the Jicamarca Radio Observatory (JRO) in Peru. Modeling the S4 value as 0.4 is consistent with the S4 value simulated by [
27] in a similar simplified plasma bubble scenario.
When a formation of two small satellites with on-board GNSS receivers fly through this ESF region, TEC and S4 between the GNSS constellation(s) and the two spacecraft can be measured. This scenario is illustrated in
Figure 3. If the dual spacecraft are flying in proximity at different altitudes, the vertical
Ne in between can be estimated by dividing the vertical TEC difference from the vertical separation. Compared to some single satellite missions (e.g., LEO GNSS tomography, radio occultation) or ground based remote sensing techniques, this LEO formation flight concept is more advantageous to observing the global and microscale morphology (e.g., precise structure and location) of irregularities that potentially induce GNSS scintillations. Also, ionospheric disturbances associated with Tsunamis or earthquakes can potentially be detected using such a multi-point ionospheric measurement technique, which is applicable to natural hazard monitoring or prediction. The vertical boundaries of irregularities can be estimated by examining the TEC,
Ne, and S4, and the horizontal boundaries can be estimated by analyzing the TEC gradients and S4. To cater the observation needs for this scenario, a low inclination (∼10°) orbit was selected to continuously monitor the equatorial region. A small-eccentricity (0.044) orbit is considered to have both spacecraft altitudes vary with time. The initial orbit states of the chief and deputy spacecraft are listed in
Table 1. A constant in-track offset of 15 km and radial offset of 1 km was predetermined as the default desired relative orbit configuration during formation keeping, however, elliptic natural orbits should be implemented in the future for a better fuel budget.
An HIL simulation of this ESF scenario was run using the VTFFTB 2.0. The history of relative orbit (radial and in-track) and thrust (radial, in-track, and cross-track) in the Hill’s frame are presented in
Figure 4. As shown in
Figure 4a, it took 125 s to reach formation keeping in the radial direction (offset within ±2 m); and it took 60 s for the in-track relative distance to reach less than ±2 m offset in
Figure 4b. As shown in
Figure 4c, thrust in radial and in-track directions were “on duty” to perform formation acquisition (correct the relative orbit from the initial state into the desired state). The first available EKF-based estimation was available at 276 s, where the EKF initialization led to a transient orbit deviation and redundant thrust as well. The filter converged quickly after that, therefore, inconsequential thrust transients are observed for the rest of the simulation. When the satellite fleet entered the ionospheric irregularity (ESF) region around 2200 s, scintillation impacted the navigation performance and disturbed the relative orbit. Throughout this one hour HIL simulation with a maximum thrust limit (
) of 0.5 m/s
, the total thrust used in radial, in-track, and cross-track directions are 22.4665 m/s, 7.6677 m/s, and 1.1758 m/s, respectively.
3.2. Gaileo TEC and Ne Measurements
When the chief satellite soars higher than the deputy satellite, the vertical
Ne between the two satellites retrieved from a specific PRN can be computed by the following formula:
where VTEC
is the vertical TEC measured by the chief satellite GNSS receiver from a PRN, VTEC
is the vertical TEC measured by the deputy satellite GNSS receiver from the same PRN, and ∆
h is the height difference between the chief and deputy. Using this method, space-born GNSS receivers can be utilized to retrieve localized
Ne given a formation flying orbit with non-zero radial offset. If the measurement accuracy is comparable to or outperforms other measurement techniques (e.g., in-situ Langmuir probe or radio-occultation), the design of a satellite mission (e.g., cost, size, and power budget) can be significantly reduced while the same
Ne product can still be generated.
The vertical
Ne retrieved from Galileo PRN 30 is plotted against the central height between the two LEO satellites as shown in
Figure 5. The TEC result using the E1 and E5a frequency combination is plotted on the left (a), while the result from the E1 and E5b frequency combination is plotted on the right (b). The SR
Ne profile is plotted in red, and the measurements using Equation (
2) are plotted in blue. Note that, these are not localized (same latitude and longitude) vertical
Ne profiles but a vertical projection of the three-dimensional trace of the sounding rocket vertical
Ne. Due to a high frequency scattering effect from the raw signals, a low pass filter (critical frequency = 0.01 Hz) is applied to process the raw measurements and the filtered results are plotted in black. After filtering, the
Ne measurement results from the HIL simulation become much more consistent with the SR model. In
Figure 5b, the measurement results above 500 km are relatively deviated from the model. This is due to the inaccurate DCB estimation with respect to PRN 30 of the E1 and E5b frequency combination. Further analysis indicates the scattering features of the
Ne raw measurement is a heritage from TEC, which is associated with the particular receiver model and frequency combination of a specific GNSS constellation. For the NovAtel OEM628 receiver, GPS L2P and L2C tracking is aided by GPS L1. While Galileo E5a, E5b, and E5 AltBOC tracking is unaided, as is GPS L5. The unaided carrier-phase TEC combinations (e.g., E1 and E5a and E1 and E5b) are noisier than the aided ones (e.g., L1 and L2). Therefore, the vertical
Ne retrieved using GPS L1 and L2 TEC by the OEM628 receiver is less fluctuated than the vertical
Ne retrieved by Galileo (E1 and E5a or E1 and E5b) TEC in the baseline ESF scenario.
3.3. Multi-Constellation Data Fusion
An accuracy improvement on
Ne retrieval is found by combining both GPS and Galileo measurements. As shown in
Figure 6a, the
Ne retrieved from selected GPS (L1 and L2) is combined with the Galileo (E1 and E5b)
Ne from a simulation using the 13 July 2018 almanac, when there were 13 operational Galileo satellites. After taking an average of those selected GPS and Galileo PRNs (with outlier-free
Ne measurements), a 0.01 Hz low pass filter was applied to plot the final values in black. The SR
Ne profile is plotted in red to represent the “true” values. The measurement errors (discrepancy between filtered
Ne and “true” value) are plotted in green. It is found that the mean and standard deviation of measurement errors were decreased by 32.83% and 46.12%, respectively, compared to just using the GPS constellation in this simulated scenario.
As expected, an accuracy improvement can also be seen by simulating the GPS and Galileo almanac in a more recent time. As shown in
Figure 6b, the
Ne retrieved from selected GPS (L1 and L2) is combined with the Galileo (E1 and E5b)
Ne from a simulation using a more recent 8 March 2019 almanac, where there were 22 operational Galileo satellites. It is found that the mean and standard deviation of measurement errors were decreased further by 43.34% and 49.92%, respectively, compared to using only the GPS data in the simulated scenario. This accuracy improvement is consistent with the growing number of Galileo satellites and indicates a benefit of applying multi-constellation GNSS SFF for ionospheric remote sensing. Adding a new constellation can improve the geometry and increase the number of line-of-sight measurements, which is beneficial to smooth out the bias in each measurement.
Wavenumber spectra of the
Ne measurements were generated to analyze the retrieval resolution on multiple spatial scales of ionospheric plasma structure, based on the 8 March 2019 Almanac and the SR
Ne profile. As shown in
Figure 7, the wavenumber spectrum of the SR
Ne model is plotted in red, the wavenumber spectrum of
Ne measurements using L1 and L2 TEC from all visible GPS satellites is plotted in blue, and the wavenumber spectrum of
Ne measurements using TEC from all visible GPS satellites plus eight Galileo satellites is plotted in green. The wavenumber
k is a function of spatial scale
as
k = 2
/
. Ne(
k) is the spatial Fourier transform of the
Ne profile. The three plasma bubbles can be clearly seen in the spectra at 0.1 < k < 0.2, which corresponds to 30 km <
< 60 km. Compared to the blue spectrum, the green spectrum is closer to the red spectrum across the whole spatial range except a noisy spike near
k = 6. This, once again, demonstrates an improvement of ionospheric electron density retrieval resolution using more satellites from multi-constellation GNSS.
3.4. Formation Configuration Sensitivity Study
A sensitivity study was undertaken to investigate the impact of the variation of the LEO satellite radial offset on the vertical Ne measurements. Based on the 8 March 2019 almanac, three formation flying configurations were simulated respectively with different radial separations during formation keeping: (i) 100 m, (ii) 1000 m, and (iii) 3000 m. The Ne retrieval using Galileo PRN 1 and the E1 and E5b TEC are compared between the three configurations.
By simulating configuration (i) with a radial offset of 100 m in the ESF scenario, very noisy
Ne measurements are obtained as shown in
Figure 8i. The measurement errors (in green) were computed by differencing the SR modeled
Ne (in red) and raw measured
Ne (in blue). The average absolute error is
cm
−3, and standard deviation of errors is
cm
−3. The raw
Ne measurements here are too noisy to distinguish the three plasma bubbles (located approximately at 340–380 km, 395–435 km, and 480–520 km, respectively). The main reason causing such scattering features is the relatively short altitude separation (∼100 m) between the LEO satellites. The VTEC measured at each altitude contains a certain amount of noise and such TEC noise levels are comparable to the real TEC difference between 100 m. Therefore, 100-meter radial separation is not large enough to resolve the E1 and E5b TEC noise from the OEM628 receiver.
By simulating configuration (ii) with a radial offset of 1000 m in the ESF scenario,
Ne measurements with a limited amount of noise are retrieved as shown in
Figure 8ii. This is the same formation configuration implemented in
Section 3.2 and
Section 3.3. The average absolute error is
cm
−3, and standard deviation of error is
cm
−3. The raw
Ne measurements look much better than the result of configuration (a), however, applying a low pass filter can further smooth the results and help visualize the plasma bubbles as demonstrated in the earlier two subsections. By simulating configuration (iii) with a radial offset of 3000 m in the ESF scenario,
Ne measurements are retrieved without the data scattering as shown in
Figure 8iii. The average absolute error is
cm
−3, and standard deviation of error is
cm
−3. The raw
Ne measurements are “clean” enough to distinguish the three plasma bubbles without further applying filtering.
A wavenumber spectrum analysis was performed to compare the resolution results between the three different configurations. As shown in
Figure 9, the wavenumber spectra of
Ne retrieval and SR model are plotted together to determine the multi-scale resolution using each formation configuration in this scenario. Clearly, the 100-m case exhibit the worst performance across the whole spatial range. The 1-km and 3-km cases show comparable agreements with the spectrum of SR
Ne model. The agreement for the 1-km case outperforms the 3-km case for smaller scales (
k > 1), vice versa.
The two GNSS receivers vertically separated by 3000 m in LEO allows a near optimal Ne retrieval for the larger scale plasma bubbles using Galileo TEC, however, this will decrease the resolution of GPS Ne retrieval because the GPS’s TEC noise level is lower than Galileo’s. Therefore, the 1000 m radial separation was chosen due to a (measurement noise) trade-off between the GPS and Galileo given the OEM628 receiver. The results indicate there is an optimal range of vertical separation (between the satellites in formation) to resolve the spatial scale of a specific ionospheric structure. In future mission designs, the VTFFTB can offer a HIL simulation testbed to determine an optimal geometry of formation flight for ionospheric remote sensing with respect to specific GNSS receivers.
3.5. C/N Level Sensitivity Study
Another sensitivity study was conducted to investigate the impact of the GNSS signal C/N level on the vertical Ne measurements. The GNSS antenna gain for the incoming signals to the receiver on a small satellite can be unstable or below the anticipated level in some scenarios, such as, satellite attitude control error or failure, spacecraft (power system) in a low power mode, and degraded satellite antenna. Any anisotropy in the GNSS antenna gain will imply the tracking PRN C/N will be different because the GNSS satellite signals come from different directions with different relative distances. Also, the power levels of GNSS signals transmitted by different GNSS constellations or different frequency bands from the same constellation are different. Moreover, ionospheric scintillations can impact the C/N level as well. Therefore, it is important to characterize the effect of C/N level on vertical Ne retrieval.
By adjusting the global power level offset in the GNSS simulator, different C/N
level tracked by the receiver can be effectively simulated with the isolation of other controlled variables. The result in
Figure 8ii is chosen as the reference (power level). Two scenarios of global power offset (applied to all PRNs) on the deputy satellite’s antenna were simulated: (a) −3 dB and (b) −8 dB. The HIL simulation comparison used the default formation configuration based on the 8 March 2019 almanac.
The simulation results show that the lower the C/N
is, the more the TEC measurement fluctuates. As a result, more measurement noise is introduced into the vertical
Ne retrieval.
Figure 10a shows the vertical
Ne retrieval for PRN 1 using Galileo E1andE5b TEC with a power offset of −3 dB, where the average absolute error is
cm
−3 and the standard deviation of error is
cm
−3.
Figure 10b shows the same vertical
Ne retrieval except with a power offset of −8 dB, where the average absolute error is
cm
−3 and the standard deviation of errors is
cm
−3. By comparing to the reference result in
Figure 8ii, the absolute mean and standard deviation of
Ne measurement errors becomes larger when the signal power (C/N
) decreases.
A comparison using the GPS constellation indicates a similar trend of decreasing accuracy that is in line with the Galileo comparison described above. Choosing the vertical Ne retrieval results using GPS L1 and L2 TEC from PRN 7, the absolute mean and standard deviation of raw measurement errors for the reference power level scenario are cm−3 and cm−3, respectively. For the −3 dB power offset scenario, the absolute mean and standard deviation of Ne measurement errors increase to cm−3 and cm−3, respectively. For the −8 dB power offset case, the absolute mean and standard deviation of measurement errors further increase to cm−3 and cm−3, respectively.
This sensitivity study demonstrates a correlation between the C/N level and the vertical Ne retrieval characteristics: When the C/N is reduced, the Ne retrieval accuracy is reduced and associated with more scattering in the Ne measurements. During the antenna selection process for a LEO satellite with a GNSS receiver, a link budget calculation including the C/N is important in characterizing the ionospheric remote sensing capability.
3.6. Galileo Scintillation Measurements
Other than TEC and electron density, adding the Galileo system will produce more types of GNSS scintillation observations (e.g., different line-of-sights, frequency bands, signal power levels, and modulation schemes), which has advantages for multi-scale ionospheric remote sensing investigations. Several Galileo scintillation observations from a multi-constellation HIL simulation of the baseline ESF scenario based on the 13 July 2018 almanac are demonstrated here.
Figure 11a–c present the vertical distribution of Galileo E1, E5a, and E5b effective 1-Hz S4 observed by the chief satellite receiver from all eight visible PRNs.
Figure 12a–c show the horizontal view of Galileo E1, E5a, and E5b effective S4 level observed by the chief receiver following the satellite ground track, as the “snapshots” of 1-Hz S4 taken every minute from eight different visible PRNs. Similar observations from the deputy receiver of the fleet (not shown) can be obtained with similar macroscale behavior as the chief. In total, 10 Galileo PRNs were tracked by the deputy LEO receiver during the 1 h ESF scenario. Since the amplitude scintillation modeled in the GNSS simulator is confined in a cuboid region, the high S4 observations from Galileo mostly occur in the corresponding region (290–350 km, 20-21 LT, 0–10°S) near the JRO as anticipated. Therefore, the observation of S4 by SFF can potentially be designed as a tracer or new indicator of ionospheric irregularities.
Because the GNSS simulator currently is only capable of modeling a simplified ionosphere, the amplitude scintillation affects all frequency bands of Galileo and GPS similarly. A multi-constellation, multi-frequency software GNSS receiver was developed and applied to ground-based ionospheric scintillation studies by [
28]. Reference [
29] utilized a ground-based, multi-constellation, multi-band GNSS data collection system to observe equatorial scintillation. Reference [
30] analyzed the multi-frequency responses of GNSS receivers to ionospheric scintillations and characterized the difference in receiver behavior under scintillations across multiple frequency bands. This is currently an important area of investigation and, therefore, space-based, multi-constellation, multi-frequency-band GNSS scintillation measurements will facilitate studying multi-scale ionospheric irregularity impacts on GNSS signals from the perspective of different constellations, frequency bands, modulation schemes, chipping rates, powers, etc. [
31].
3.7. Relative Navigation Improvement
The relative navigation robustness is very important when an SFF mission rigorously demands the separation accuracy among their distributed space systems. For instance, the precision of formation keep and relative state estimation for a 1000-m baseline SFF mission is sensitive to the spatial resolution of small-scale (e.g., ≤100 m) ionospheric irregularity or large-scale ionospheric morphology (e.g., size or boundary of plasma structure) observations.
Besides ionospheric remote sensing, the new version of VTFFTB after incorporating Galileo also shows an improvement on relative state estimation. The baseline ESF scenario described in
Section 3.1 is used for HIL simulations, in order to compare the EKF-based relative state estimations between a GPS scenario and a multi-constellation scenario. A number of qualitative comparisons were performed, and one comparison based on the 8 March 2019 almanac is shown here as an example. The relative state estimation errors presented in
Table 2 are based on a simulation using the GPS (L1) constellation only for relative navigation. The relative state estimation errors in
Table 3 is based on a simulation using both the GPS (L1) and Galileo (E1) constellations. The error is defined as the difference between the EKF estimated value and the simulator recorded “truth” after each simulation, and it is calculated including the initialization period before the EKF converges. All the error values in
Table 2 and
Table 3 are rounded to three decimal places. Note that, the scintillation effect is not fully simulated and the induced cycle slips were not fully addressed in the HIL simulation for simplicity. This is likely the main factor degrading the overall relative estimation performance. However, several non-control simulations on the VTFFTB with more robust algorithms (to handle cycle slips) were used for formation flying relative navigation and their relative position errors were less than 10-cm after the EKF converges.
Table 4 gives the computed error decrease percentages of
Table 3 relative to
Table 2, which quantifies how much relative state estimate errors were decreased by using multi-constellation GNSS (GPS L1 and Galileo E1) against using GPS only. The positive values indicate error decrease, while a negative value indicates an error increase. After a large ensemble of simulations, it is found that the average absolute errors in all ECEF (Earth-centered, Earth-fixed) components of relative position and relative velocity typically decreased. The standard deviation of errors decreased in most cases as well. The simulations based on a newer almanac usually showed better improvement as well due to a larger number of Galileo satellites. Also, the EKF-based relative estimation converge time is typically faster using multi-constellation, compared to using GPS-only. Therefore, relative navigation with both GPS and Galileo constellations can offer better estimation performance, as the dimension of differential GNSS measurement matrix becomes larger with more PRNs.
4. Conclusions
The GPS-based VTFFTB has been successfully upgraded into a multi-GNSS (GPS + Galileo) version. By running HIL simulations of a baseline ESF scenario, the ionospheric remote sensing and EKF estimation capabilities were found to be improved with the addition of Galileo. By comparing to just using the GPS constellation for ionospheric sounding, using GPS (L1, L2) and Galileo (E1, E5b) together decreased the mean and standard deviation of vertical Ne retrieval errors by 32.83% and 46.12% when simulating the 13 July 2018 almanac. A more recent simulation of the 8 March 2019 almanac shows the mean and standard deviation of vertical Ne retrieval errors were decreased further by 43.34% and 49.92% when multi-constellation (GPS and Galileo) data were utilized instead of implementing GPS-only. Using the NovAtel OEM628 receiver, vertical Ne retrieved by the Galileo constellation using E1 and E5a or E1 and E5b TEC are nosier than the Ne retrieved by the GPS constellation using L1 and L2 TEC. The characteristics of retrieved vertical Ne are sensitive to the altitude separation between the two satellites in formation, and level. Sufficiently small altitude separation increases measurement noise, while sufficiently large altitude separation reduces spatial resolution. Lower level decreases the vertical Ne retrieval accuracy and increases the measurement variance. Ionospheric scintillations from more frequency bands, different power levels and modulation schemes can be observed with the new Galileo system. This offers more opportunities for future GNSS-based SFF missions to detect multi-scale ionospheric irregularities. Also, the relative state estimation performance is increased by using both GPS L1 and Galileo E1 data for differential carrier-phase measurement compared to using GPS L1 only.
In order to place onboard GNSS receivers at different altitudes, the formation configurations simulated in the ESF scenario feature a constant radial offset during formation keeping, which is not fuel-efficient especially for small satellite formation flight. Natural relative orbits with bounded motion will be implemented to optimize the fuel budget. To realize an altitude difference between satellites, elliptic or circular relative orbits should be considered and accessed to characterize the Ne retrieval capability while varying the altitude difference between two LEO satellites. Also, the simulation fidelity of the ionosphere requires some improvement to emulate the ionospheric impacts on GNSS signals more realistically. This will facilitate the design and development of new applications to a wider variety of ionospheric phenomena, such as sub-auroral polarization streams (SAPS), polar tongue of ionization, polar cap patches, ULF waves, etc. Last but not least, more GNSS receivers and simulators are being incorporated into the VTFFTB in order to simulate and study multiple-satellite (≥3) real-time formation flight and its applications to multi-scale space weather observation. This will be reported on in the near future.