1. Introduction
The Global Geodetic Observing System (GGOS) has been established by the International Association of Geodesy (IAG) in order to satisfy the expected future requirements of science and society, which are facing increasing challenges on a changing planet. GGOS aims to establish a terrestrial reference frame with an accuracy on the level of 1 mm or better on a global scale [
1]. This can only be achieved by a rigorous combination and integration of the different ground and space geodetic techniques. One of the most accurate realizations of a reference system aiming at this goal is the current International Terrestrial Reference Frame (ITRF2014) [
2]. As a combination product, the ITRF makes use of the full history of observation time series of the four space geodetic techniques Very Long Baseline Interferometry (VLBI), Global Navigation Satellite Systems (GNSS), Satellite Laser Ranging (SLR), and Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS). In addition, the ITRF determination fundamentally relies on differential coordinates—called local ties—connecting the reference points of the geodetic instruments that are typically located within distances of a few hundred meters at co-location sites (for a definition of local ties, we refer to [
3]). These local tie vectors are usually determined by classical terrestrial measurements (angles plus distances and leveling) or by GPS measurements. Although local ties play a critical role in the ITRF combination by linking the contributing techniques (along with Earth rotation parameters), comparisons of the local ties and differential coordinates determined by various space geodetic techniques reveal significant discrepancies with residuals on the centimeter level in many cases. This issue was discussed in detail among others by Altamimi et al. [
4], Altamimi et al. [
2], Seitz et al. [
5], and Thaller et al. [
6].
As proposed, e.g., by Rothacher et al. [
7], an additional way to combine space geodetic techniques—complementary to local ties at co-located sites on the Earth’s surface—is the co-location of geodetic instrumentation on board Earth orbiting satellites. Such so-called “space ties” (in accordance with the term local tie) would provide unprecedented possibilities to connect the coordinate frames of different space geodetic techniques, and they would help to reveal unknown technique-specific biases. Alternative approaches to establish frame ties are the co-location of VLBI with other techniques on the moon [
8], and the application of so-called global ties [
9], i.e., of commonly estimated parameters in different technique solutions of space geodesy.
So far, there have been several successful attempts to combine the three satellite techniques GNSS, SLR, and DORIS via co-location on Earth-orbiting satellite platforms. Thaller et al. [
6] demonstrated that it is feasible to combine GNSS and SLR by utilizing GPS and GLONASS satellites equipped with SLR retro-reflectors as co-location platforms in space. In a more recent study by Sosnica et al. [
10], satellites of other GNSS constellations (BeiDou and Galileo) were used to establish the SLR-to-GNSS frame connection. Zoulida et al. [
11] studied the combined analysis of GNSS, SLR, and DORIS on board the Jason-2 satellite with the goal to set a basis for the integration of a low Earth orbit (LEO) multi-technique satellite such as Jason-2 as space tie in the TRF determination.
Contrary to that, space ties with VLBI have not been established so far and have solely been able to be studied in simulations: Plank et al. [
12] investigated general aspects such as suitable satellite orbits and antenna networks, finding that station position estimates at the level of a few millimeters in weekly solutions are achievable based on VLBI satellite observations only. Initial scheduling strategies for VLBI observations of the GNSS constellations were simulated and investigated by Plank et al. [
13]. Anderson et al. [
14] simulated VLBI observation, time-of-flight measurements using a time-encoded signal in the satellite transmission, and differential VLBI observations using angularly nearby quasars as calibration sources, and investigated their suitability for the estimation of frame tie parameters.
However, the actual link between the satellite techniques GNSS, SLR, and DORIS and VLBI has not yet been established. The essential problem in this respect is the absence of real satellite observations with VLBI radio telescopes. While GNSS, SLR, and DORIS are operationally applied for observations of Earth satellites, geodetic VLBI is intrinsically different, with its roots being settled in astronomy. Present applications of VLBI satellite tracking can mainly be found in the field of planetary sciences, with VLBI being used for the navigation and positioning of spacecrafts (for example: [
15,
16,
17]). The requirements on observation schemes and data processing essentially differ between VLBI observations being carried out for the purpose of interplanetary navigation and for establishing frame ties via Earth satellites. Hence, observation and processing routines established for planetary science cannot simply be adopted for use in geodetic VLBI.
In the past few years, several groups within the geodetic community reported an effort on observing L-band GNSS satellite signals. Initial experiments were carried out to investigate basic signal chain characteristics and tracking capabilities mainly involving the European observatories in Wettzell (Germany), Onsala (Sweden), and Medicina (Italy) (described in, e.g., [
18,
19,
20,
21]). In an extensive series of experiments with observations of GNSS satellites (L1 and L2 signals) on the Australian baseline Hobart–Ceduna the complete process chain—from scheduling, actual observations, correlation and fringe fitting, to the final analysis of the derived group delays—could be closed for the first time [
22].
One main challenge in these experiments was the standard equipment of legacy geodetic VLBI antennas, being restricted to observations in the S- and X-bands and not covering the L-band frequency range. However, this may change in the near future with the advent of broadband antenna feeds and matching receivers following the VGOS design concepts [
23]. VGOS-compatible antenna equipment is expected to cover the frequency domain from 2 to 14 GHz. This will enable more flexibility in the selection of observation targets, though the frequency range again does not cover the L-Band of GNSS satellites. Hence, scientists had to fall back on antenna equipment that is not operationally used for geodetic applications involving the risk of introducing unexpected biases and unknown systematic effects caused by the non-standard signal chain.
Two dedicated co-location satellite missions were proposed by the geodetic community in recent years that would employ all four techniques used in the TRF combination: the Geodetic Reference Antenna in Space (GRASP [
24]) and the European Geodetic Reference Antenna in Space (E-GRASP/Eratosthenes [
25].) The GRASP satellite was planned to be equipped with a VLBI beacon emitting signals in the S- and X-bands that are compatible with frequency ranges of legacy and future VLBI antennas. The signals should have been encoded using Binary Offset Carrier (BOC) modulation to enable range measurements. The E-GRASP mission was proposed to carry a broadband noise beacon emitting signals in the S-, X- and Ku-bands.
Since a GRASP-like satellite has not yet been realized, alternative observation targets must be employed for hands-on observation experiments. Actual VLBI satellite tracking experiments are considered valuable for future co-location satellite missions as they provide plenty of hands-on expertise, complementing the information derived from theoretical simulation studies. Only actual tests are capable of revealing certain deficiencies on the data acquisition level, such as tracking issues, and allow for the experience required to adopt existing VLBI infrastructure components—hardware as well as software—according to the needs of satellite observations. Hence, today’s experiments—although they may not deliver highly accurate results—provide a great test bed for future missions with potentially a high impact on geodesy.
In this paper, we discuss VLBI observations of the APOD-A nano satellite (hereinafter abbreviated as APOD). APOD can be considered as the first realization of a co-location LEO satellite combining the techniques SLR, GNSS, and VLBI on one platform (details are provided in
Section 2.2).
We describe a set of intensive tracking sessions to APOD using the AuScope VLBI array in Australia [
26]. We focus on the data acquisition scheme for observations of a fast moving LEO satellite on regional baselines and discuss the resulting challenges for the tracking procedures and the data processing.
The established process chain for our experiments is outlined in
Figure 1. All tasks and processing steps are described in detail in the following sections.
Section 2 provides an overview of all experiments discussed in this paper, and introduces the station network and the observation target. The observations, including all steps from scheduling with the Vienna VLBI and Satellite Software (VieVS [
27]) to the actual data acquisition, are discussed in
Section 3. In
Section 4, the correlation and post-processing steps, based on the software correlator DiFX [
28] and the Haystack Observatory Post-processing System (HOPS,
http://www.haystack.mit.edu/tech/vlbi/hops.html), are described. The analysis of the derived delay observables using VieVS is presented in
Section 5. The paper is concluded with a discussion and future prospects in
Section 6.
3. Scheduling and Observations
The starting point of each VLBI experiment is the observation planning, referred to as scheduling. In general, VLBI schedules define which antennas observe which source at what time. The scheduling task itself is complicated due to the large number of different observation criteria that have to be considered, such as common visibility of the target from remote sites, slew times between consecutive scans, and the determination of required on-source time to reach the target SNR. For planning the discussed experiments, the VieVS satellite scheduling module [
32] was used. It allows for the scheduling of VLBI observations of near-field targets along with observations of natural extragalactic radio sources, which are routinely observed in geodetic VLBI campaigns. The natural sources are selected automatically by optimizing the sky coverage at all stations as common for geodetic VLBI sessions [
33].
The observation geometry that was determined by continental-wide baselines and the very low orbit of APOD was a major limiting factor for scheduling. On average, the APOD satellite was visible by individual AuScope telescopes four times a day for a couple of minutes only. The projected fields of view are displayed in terms of shaded red circles in
Figure 2. Only in the intersecting areas was APOD simultaneously visible by two stations. As indicated in
Figure 2, common visibility from all three AuScope telescopes was restricted to very low elevation angles and scan durations as short as a couple of minutes at most. In general, APOD was simultaneously visible by only two of the three antennas, limiting the experiment design to single-baseline scans. Due to these restrictions, it was not possible to observe more than two single-baseline scans shortly after another during a flyover, as exemplarily indicated in
Figure 2 for two consecutive scans (168 and 169) in Experiment a332. In this example, APOD was observed first on the baseline Hb–Yg for 86 s followed by a second scan on Ke–Yg lasting for 256 s. The next occasion for an APOD scan only existed about 10.5 h later that day. Despite these circumstances, we took every chance from November 11 to 14, 2016, and observed APOD whenever common visibility from two AuScope antennas was given. All observed experiments are listed in
Table 1, which also provides details on the observation durations. During observations, APOD was tracked continuously by the antennas as long as common visibility was given and the signal was recorded continuously. (Each of these recorded APOD tracks is referred to as
scan in this paper, which must not be mixed up with the term
observation. Due to the high SNR, delays could be derived in a 1 s interval, i.e., a single scan is chopped up into multiple observations.)
Due to the observation geometry APOD scans were restricted to low elevation angles in general with most experiments being observed at elevations between 5 and 20
, at some rare occasions at higher elevations up to 37
at most. In general, the data analysis becomes more challenging due to the consistently low elevations, as the effects of the neutral atmosphere and the ionosphere on the observed signals increase with decreasing elevation. Furthermore, estimates of station heights, zenith wet delays, and station clocks do not decorrelate properly without observations at varying elevations (e.g., [
34]). On the other side, the change rates of the topocentric antenna pointing directions (azimuth and elevation rates) decrease as the distance between antenna and satellite increases with lower elevation angles. Hence, the demands on the tracking and on the pointing accuracy of the antenna decrease at lower elevation angles. Therefore, the APOD satellite was easier to keep within the field of view via the observing antennas at low elevations than at high elevations.
We were able to assess the received signal live during APOD tracks by a spectrum analyzer connected to the intermediate frequency (IF) channels at the recorder racks. Our experience was that the signal amplitude became increasingly unstable at higher elevation angles, especially in the X-band due to the narrower field of view. While the S-band signals were rather stable throughout all experiments, we even experienced signal loss in the X-band occasionally at elevations above ∼30
. We think that these tracking issues were mainly caused by the low accuracy of the predicted APOD orbits used for tracking, which show offsets to the final orbits of up to 1000 m (see
Section 2.2). At an elevation of 30
the max. pointing error caused by an offset of 1000 m in the APOD orbit was about 4.2
. With a beam-width of about 10.2
in the X-band, pointing errors of a few arc-minutes can already be critical. Furthermore, the internal interpolation of the tracking data in the ACU might not be accurate enough to precisely follow such a fast satellite. If there is already a pointing offset caused by the tracking data, additional inaccuracies due to a non-optimal interpolation of the data in the ACU may be enough to cause severe pointing problems in the X-band—even at moderate elevation angles.
On average five quasars were observed before and after the APOD track(s) in Experiments 316a to 319a. (For the sake of brevity, we only mention quasars when we talk about natural radio sources, as the most common type of active galactic nucleus (AGN) observed in the geodetic VLBI. To be precise, we do not exclude the possibility that other types of AGN were observed as well in the described experiments.) The main reason for the inclusion of quasars in the observation plan was to use them as calibrator sources in order to establish an initial clock model in the correlation process as outlined in
Section 4.2).
Experiment a332 was designed differently: Basically, it is a geodetic 24 h VLBI session consisting of 761 three-station scans to strong quasars with a minimum flux density of 0.65 Jy. For the scheduling we optimized the sky coverage at each site in order to decorrelate estimates of station clocks, station heights, and troposphere delays. Whenever visible, APOD was observed in between the quasars, resulting in four single-baseline scans to APOD in 24 h. The first two scans (168 and 169), which are illustrated in
Figure 2, are used as a generic example in this paper. This particular session design with satellite scans being embedded in quasar scans allows for additional analysis options, such as using Zenith Wet Delays (ZWDs) estimated in a quasar-only solution, to correct the APOD observations a priori (for more details, we refer to
Section 5).
The observing mode used for all described experiments was designed pursuing two purposes: (1) to record the full APOD S- and X-band signals and (2) to derive reasonable geodetic results from standard observations of quasars at the same time. Changing the observation mode when switching between quasars and APOD was not an option as it might have introduced unknown systematic biases. The observation mode has been slightly modified from the mode used for the AUSTRAL sessions described in [
35]. Data was recorded in 16 channels with 16 MHz bandwidth each, 10 in the X-band and 6 in the S-band, applying 2-bit sampling, which results in a recording rate of 64 Mbps per channel. The frequency allocation is illustrated in
Figure 3. Due to strong Radio Frequency Interference (RFI) in the lower S-band at Hobart, the S-band channels were allocated contiguously, yielding continuous frequency coverage. Only 4 of 16 channels were dedicated to record the APOD tones. The purpose of the others was to increase the frequency coverage for quasar observations to enable the calculation of reasonable multi-band delays.
Publicly available two-line element (TLE) data sets (e.g., at
https://celestrak.com/) were used for the initial orbit determination by VieVS and enable one to plan the experiments up to several weeks in advance. However, a final schedule iteration is recommended shortly before a session, as the observation times may slightly change with updated orbit elements. VieVS was updated for the APOD observations to generate the station-specific AZEL tracking files enabling the AuScope antennas to continuously track satellites (see
Section 2.1). The data points for the tracking files were determined in a 1 s interval with the calculations being based on the latest APOD orbit predictions provided by BACC about 12 h before an experiment. The observation plans were written to VEX formatted schedule files (see
http://www.vlbi.org/vex/), which are standard in VLBI. They contain a full description of the experiments, including all receiver and recorder settings for the whole antenna network, and the observation time schedule. The participating stations extract all relevant information from the VEX file using the Field System program drudg. Drudg creates station-specific control files which contain all required commands to run the experiment fully automated by the Field System, including antenna motion control, calibration routines, and the setup of the signal chain.
Only one major modification from this standard procedure was required to enable satellite tracking: Prior to satellite scans, the tracking mode of the ACUs had to be manually switched over from the star tracking mode, which is used per default for astronomical sources, to the AZEL tracking mode described in
Section 2.1, and the prepared AZEL tracking files had to be loaded. These setup changes had to be done individually for all observing antennas. Manually changing the tracking mode only affected the antenna motion control, while the signal chain was still controlled by the Field System based on the setup parameters defined in the VEX file. To incorporate these additional steps, about 5 min of idling time was defined in the observation schedules prior and after APOD tracking.
The raw data were recorded with Mark5B+ recorders and were shipped to the computation cluster in Hobart for correlation and post-processing.
5. Data Analysis
The multi-band delay observables derived by the processing scheme outlined in
Section 4 were analyzed with VieVS 3.0., which is able to handle quasar observations, as well as observations of satellites. All target parameters can be estimated in a least-squares adjustment (Gauss–Markov model, e.g., [
39]) based on reduced observations, which are calculated by subtracting the computed (theoretical) delays from the observed delays (observed minus computed, o-c).Theoretical delays for APOD observations were calculated with the near-field delay model outlined in
Section 4.1 using the final orbit solution by BACC (see
Section 2.2).
The analysis is based on the so-called total delays, which were computed as the sum of multi-band delay residuals calculated in fourfit and the a priori delays of the correlator input model. The observation reference epochs are the times of the signal reception at Station 1 (of the baseline), calculated at integer seconds, as common in geodetic processing. Basically, we have three sets of delay results for each experiment: (1) S-band delays, (2) X-band delays, and (3) the ionosphere free linear combination. To minimize the effect of the ionosphere, we decided to use Option 3, although the quality of the S-band observations—in terms of amplitude stability and SNR (see
Section 4)—seems to be higher.
Multi-band delays were also obtained from the quasars observed in Experiment a332 by standard VLBI data processing with DiFX and fourfit. These quasar delays were written to an Mk3 database, processed in
Solve [
40] and exported as a Level 4 NGS file. They were then analyzed in VieVS by estimating a set of standard parameters: station coordinates, source coordinates, EOP, ZWD, and station clock parameters. The WRMS of the post-fit residuals was 38 ps. The obtained ZWD estimates shown in
Figure 11 were then applied a priori on the APOD observations to correct for effects caused by the wet fraction of the atmosphere. The Vienna Mapping Function (VMF1 [
41]) was used to determine the slant delay corrections.
Investigating the o-c delay residuals provide a first impression of how well the observed delays match the theoretical delays modeled in VieVS. The top panel in
Figure 12 depicts o-c values for two consecutive APOD scans, 168 and 169, in Experiment a332. The observed delays were corrected for the station clock offsets and rates applied in the correlation step to keep the absolute o-c values small.
The bottom panel in
Figure 12 depicts the major delay corrections applied on the purely geometric near-field delays (see
Section 4.1), which target the effects of the hydrostatic (trop_h) and the wet (trop_w) fraction of the troposphere. The hydrostatic part was modeled according to [
42] using in situ pressure measurements and the VMF1
mapping function. The wet fraction was modeled by mapping the ZWD estimates derived by the quasar observations in Experiment a332 down to slant directions using the VMF1
mapping function. The total troposphere corrections are in the range between ∼−61 ns and ∼+51 ns. Other corrections not being illustrated here due to their comparably small magnitudes are thermal antenna deformation (max. 1.8 ps [
43]), gravitational effects on the near-field delays according to [
37] (max. ±24 ps) and the effect of antenna axes offsets (max. 44 ps).
In Scan 168, the residual delays show little variation of ∼1.5 ns over the scan duration. In contrast, the delay residuals in Scan 169 are noisier and they show a systematically curved variation with a magnitude of ∼12 ns over the scan. Looking at the elevation angles at which APOD was observed during these scans (
Figure 12, third panel), we suppose that the more pronounced scatter in the second scan can be explained by tracking issues that are more likely to occur at higher elevations. At first glance, the distinct bent shape in Scan 169 can lead to the assumption that this is caused by elevation-dependent effects, such as troposphere delays. However, investigations showed that this signature can be explained by a constant offset in the satellite’s along-track position present in the a priori orbit data that were used to determine the theoretical delays. The second panel in
Figure 12 depicts the o-c time series derived by the same observation data, but with different constant along-track offsets (from +10 to −10 m) applied when modeling the computed delays. Typically, the largest uncertainty in orbit determination is present in the along-track direction, especially in the case of LEO satellites due to the predominant effect of the atmospheric drag, which is constantly decelerating the satellite (e.g., [
44,
45]). Considering the low quality of the used orbit data (see
Section 2.2), deviations of 10 m (or even more) along-track can be expected. When applying a constant along-track offset of about −8 m on the a priori orbit, the bent shape disappears. The remaining constant offsets between the two scans can be explained by uncorrected clocks.
Such o-c values could in principle be used to estimate all parameters that are usually determined by the geodetic VLBI technique, plus satellite orbit parameters. Due to the severe inaccuracies in the a priori orbit data, it is mandatory to estimate orbit parameters. Otherwise, unmodeled orbit errors would propagate into other estimates. Basically, VieVS provides features to estimate satellite position offsets in terms of piece-wise linear functions, which refers to a kinematic orbit modeling approach. The presented APOD observations were used to validate and test these features. However, for kinematic orbit modeling, the available VLBI observations are insufficient. Firstly, the total number of available tracks is by far too low. Secondly, a global tracking network would be required, rather than a regional one, to acquire observations distributed over whole orbit arcs. A global network would also yield more observations in total.
6. Summary and Outlook
The Chinese APOD-A nano satellite can be considered as the first prototype of a co-location mission enabling VLBI with SLR and GNSS on an LEO satellite.Unfortunately, it has hardly been observed by VLBI, so major studies towards the actual frame ties could not be performed. The main reason for the lack of observations was that VLBI observations of satellites are non-standard, and suitable observing strategies were not in place for this mission. This work now presents the first serious attempt to observe the satellite with a VLBI network over multiple passes. We discuss a series of APOD observations with the AuScope geodetic VLBI array carried out in November 2016. We describe all steps—from the initial experiment design, to the observation and data acquisition procedures, to the correlation and post-processing scheme, and to the analysis of the derived delay observables in VieVS. By the example of Experiment a332, we discussed all applied processing steps, yielding o-c residuals on the level of a few nanoseconds.
One major limitation was the bad quality of the APOD orbit solutions due to the absence of GNSS raw data caused by malfunction of the on-board receiver. The predicted orbit data used for tracking show offsets of up to 1000 m compared to the final orbit solutions, which were used to calculate theoretical near-field delays for correlation and analysis. We assume that the low quality of the orbit data, in combination with potential deficiencies of the AZEL tracking mode provided by the antenna controllers, lead to pointing inaccuracies on the level of a few arc-minutes. The mis-pointing caused wide amplitude variations in the received X-band signals, resulting in a lower and more variable SNR of the X-band data compared to the S-band (where the antenna field of view is about four times wider). On the data processing level, the low quality of the a priori orbit limited the accuracy of the computed theoretical delays. Hence, the correlator input model was not accurate enough to fully stop phase wrapping. However, the delay model was sufficient to decrease the phase rate to such an extent that a practicable integration time could be used without introducing decorrelation.
Facing the problems caused by the inaccurate orbit solutions, it would be beneficial to use SLR observations (in addition to the available GNSS measurements) to improve the orbit determination. SLR measurements could not only be used to calibrate orbit errors in radial direction but also along-track errors due to the low elevation angles at which most of the LEO tracking data are collected. According to Arnold et al. [
46], SLR can contribute along-track corrections on the millimeter level and thus substantially improve the overall orbit accuracy. Unfortunately, SLR measurements were sparsely available at the times VLBI observations were performed, precluding this option. It is highly advisable that, for future experiments, VLBI observations are coordinated with SLR tracking campaigns.
For the sake of brevity, only one experiment is discussed in detail. Experiment a332 was chosen as a representative example. However, delay observables could be obtained from all experiments using the presented process chain. (Results in terms of total multi-band delays for all sessions listed in
Table 1 are publicly available in the IVS working group 7 Wiki at
http://auscope.phys.utas.edu.au/opswiki/doku.php?id=wg7:apod:apod_obs_auscope. In addition to the results, all session control files, e.g., VEX schedules, are provided along with operator notes.) All residual multi-band delays show similar characteristics in terms of magnitude, scatter (S-band is always smoother than X-band), and SNR (higher and more stable in the S-band). Presumably, systematic signatures in the delay residuals were mainly caused by errors in the a priori orbits used to model the computed delays. Additionally, the described tracking issues could have caused small systematic effects, as the satellite beacon was observed with slight off-axis angles. In general, the delay observables become noisier with increasing elevation, especially in the X-band, indicating that improved tracking features utilizing more accurate orbit data would be beneficial.
Considering the experience gained through the presented experiments, we provide some suggestions for future experiments. In general, the common visibility of APOD was limited to elevations lower than 40 due to the observation geometry, which was determined by continental-wide baselines (between 2360 and 3432 km long) and the low satellite orbit (about 450 km). Shorter baselines would enable one to observe longer tracks with a larger variation in the elevation angle. Nevertheless, a global distribution of observation sites is a prerequisite for the estimation of satellite orbits. Further (simulation) studies are recommended to investigate suitable observation networks for LEO satellites.
The VLBI signal emitted by APOD, consisting of narrow band carriers plus symmetric DOR tones at redundant spacing, is not ideal for obtaining delay observables as common in the geodetic VLBI, as it yields a very narrow ambiguity spacing. For classical VLBI processing, wide band noise would be preferable, although the present signal is suitable for other processing schemes, such as Doppler ranging.
Furthermore, our experiments clearly showed that adequate satellite tracking features, which enable the antennas to track continuously and accurately are a basic requirement for observations of LEO satellites. It is important to consider this when designing the future global VLBI network in the view of co-location satellite missions.
The next step towards fully automated satellite observations is the integration of the AZEL tracking mode provided by the ACU in the Field System (station-specific code). In combination with the full support of the new VEX2.0 format (
https://safe.nrao.edu/wiki/bin/view/VLBA/Vex2), which will enable one to integrate satellites in the observation schedules and control files, satellites could then be observed like quasars in geodetic sessions.
Our work showed that a sophisticated dynamic orbit model is highly beneficial for the analysis of future experiments. Considering a potentially low number of observed tracks—caused by the restricting observation geometry for LEO satellites—and a non-ideal distribution of stations, a kinematic orbit model is not suitable. To be prepared, we plan to upgrade VieVS with a suitable orbit model.
Although the results of this study are not suitable yet to study actual frame ties, this work is valuable due to the gained experience in terms of observation, data processing, and analysis strategies. The developed observation and data processing procedures can now serve as guidelines for other observatories, hopefully better preparing the global VLBI network for the next co-location satellite mission.