1 Introduction

Sentinel-6 is a multi-national satellite mission aiming to provide continuity of service for sea surface altimetry in the third decade of the current century (Scharroo et al. 2016; Donlon et al. 2021). It builds on the legacy of the TOPEX/Poseidon and Jason-1/-2/-3 missions, and is therefore also known as Jason Continuity of Service (CS) mission. The prime instrument is a Poseidon-4 Ku/C-band altimeter (Cullen and Francis 2014; Donlon et al. 2021), which offers both low- and high-resolution modes and represents a major evolution over the Poseidon-3 instruments of Jason-2 and -3. Using two frequencies, ionospheric path delays in the altimeter measurements can be compensated and the total electron content below the spacecraft can be measured. The payload is complemented by the AMR-C (advanced microwave radiometer—climate quality; Maiwald et al. 2020) instrument and a GNSS radio occultation (RO) receiver (Young 2017). For precise orbit determination (POD), Sentinel-6 carries a Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS; Auriol and Tourain 2010) receiver, a redundant pair of PODRIX GPS/Galileo receivers, and a laser retro-reflector array (LRA).

Sentinel-6A “Michael Freilich,” the first satellite of the mission, was launched on November 21, 2020 into a low Earth orbit (LEO) of 66\(^\circ \) inclination. At a nominal altitude of 1336 km, it completes one revolution in 112.5 min and its ground track repeats after 10 days (127 orbits). The spacecraft has a house-shaped body with approximate dimensions of 3.5 m (length), 2.4 m (width), and 1.8 m (height). Two solar panels with an individual size of about \({3.6}\,{\hbox {m}} \times {2.4}\,{\hbox {m}}\) are mounted in a roof-like configuration on the portside and starboard side of the spacecraft and extend well beyond the main body (Fig. 1).

Fig. 1
figure 1

Artist’s drawing of the Sentinel-6A spacecraft with the dish antenna of the advanced microwave radar and three star camera baffles on the front panel. The top panel carries three GNSS choke ring antennas for the POD and RO receivers. The three axes indicate the orientation of the spacecraft body system, with \(+x\) pointing in flight direction and \(+z\) pointing in nadir direction during nominal operations. Image credit: ESA/ATG Medialab, CC BY-SA IGO 3.0

This study aims to characterize the POD accuracy achievable for Sentinel-6A with dual-constellation GNSS tracking from the PODRIX receiver. Following an overview of the receiver, calibrations of the raw observations are discussed that are required for seamless single-receiver ambiguity resolution in the POD process. Subsequently, the need for inflight calibration of the GNSS antenna is demonstrated and the resulting phase center estimates and phase patterns for GPS and Galileo frequency bands are presented. Thereafter, an overview of measurements models and dynamical models for the POD process is given with focus on the aspects of bias handling for single-receiver ambiguity fixing as well as the modeling of non-gravitational surface forces. POD results obtained for the first 6 months of the Sentinel-6A mission are presented in the following section and discussed along with an assessment of the achieved performance using various forms of self-consistency tests and external validation with satellite laser ranging (SLR) measurements.

2 GNSS instrument and calibration

The Sentinel-6A spacecraft is equipped with a redundant pair of triple-frequency, dual-constellation Precise Orbit Determination Receivers (PODRIX, Fig. 2) developed by RUAG, Austria. They are based on a combination of second- and fourth-generation versions of the Advanced GPS/GLONASS ASIC (AGGA), which were developed under the lead of the European Space Agency (ESA) for use in spaceborne GNSS receivers (Roselló et al. 2012).

Fig. 2
figure 2

PODRIX precise orbit determination receiver (right) and patch excited cup antenna with choke ring (left). Image credit: RUAG

The PODRIX combines a single AGGA-4 chip for tracking of open GPS and Galileo signals with two legacy AGGA-2 chips for GPS C/A- and semi-codeless P(Y)-code tracking. Supported frequencies cover the L1/E1 band, the L2 band and the L5/E5a band. All signals are jointly processed on a LEON-FT micro-processor embedded in the AGGA-4 chip. In total, dual-frequency tracking of up to 18 GPS or Galileo satellites is supported, including P(Y) tracking for up to 8 GPS satellites.

Table 1 PODRIX observation types and corresponding RINEX (Romero 2020) identifiers

An overview of supported signals and tracking modes is provided in Table 1 along with the corresponding designations of the receiver independent exchange format (RINEX; Romero 2020) for GNSS measurements. In case of the modernized signals, PODRIX observations are exclusively based on tracking of the pilot component, while the navigation data bits are obtained with a slaved correlator for the data component. For the Galileo E1 open service signal, only the BOC(1,1) binary offset carrier is considered in the code replica generation, while the BOC(6,1) contribution is ignored (Roselló et al. 2012). This results in a small loss of carrier-to-noise density ratio, \(C/N_0\), compared to a tracking of the actual composite BOC (CBOC) modulation, which amounts to roughly 1 dB and can well be tolerated in the present application.

In view of various resource limitations, the set of jointly tracked GPS signals in the PODRIX receiver is kept to the necessary minimum for dual-frequency navigation, using a configurable choice of signals for the individual blocks of GPS satellites. In the default scheme adopted for the nominal receiver during the commissioning phase, C/A and L1/L2 P(Y)-code tracking is used for Block IIR satellites, while the C/A and L2C signals are tracked for IIR-M, IIF, and GPS III satellites. No use is made of L5 observations despite their superior characteristics, since these are only available for half of the current GPS constellation and would have resulted in different frequency pairs (L1/L5 vs L1/L2) across the set of tracked satellites. For Galileo, the pilot signals on E1 and E5a are tracked for all healthy satellites in nominal orbit, but tracking of space vehicle numbers SVN E201 and E202 in eccentric orbits is currently not supported due the lack of Galileo almanac data for visibility prediction. On average, GPS P(Y), GPS C/A+L2C, and Galileo E1+E5a observations were available for approximately 1.6, 6.2, and 6.3 satellites in early 2021 at a 10\(^\circ \) elevation mask.

The use of different observation types within the GPS constellation poses specific challenges for a proper POD processing and requires careful consideration of inter-signal biases. While biases related to the GPS satellites themselves can be considered based on corresponding products made available through the International GNSS Service (IGS; Johnston et al. 2017), user-related equipment biases need to be corrected based on dedicated calibrations or treated as unknown solve-for parameters. Code biases in the PODRIX receivers and antenna were found to cause a bias of roughly 4 m in the ionosphere-free (IF) combination of C/A+L2C observations as compared to L1/L2 P(Y) on Sentinel-6A, which reflects itself in a corresponding inconsistency of the estimated receiver clock offsets or pseudorange residuals. Likewise, these biases affect the Hatch–Melbourne–Wuebbena combination (HMW; Hatch 1983; Melbourne 1985; Wübbena 1985) and therefore impact the fixing of integer ambiguities in the processing of GPS carrier phase measurements.

Table 2 Signal-specific group delay corrections of raw pseudoranges (in [m]) for the nominal PODRIX receiver of Sentinel-6A as a function of the frontend temperature

Calibrations of temperature-dependent code biases for all signal types have been performed by the manufacturer for the individual flight units of the PODRIX receivers and are made available as part of the instrument calibration database (ICDB). Along with measurements of the front-end temperature provided in the receiver telemetry, these data were used by the authors to correct the code measurements prior to the generation of RINEX observation files for the POD processing from the raw GNSS receiver data. Frontend temperatures in the nominal receiver were found to vary between daily means of about 15–20\(^{\circ }\) that partly relate to the Sun’s elevation angle above the orbital plane (\(\beta \)-angle) but also show step-wise variations attributed to operational changes of onboard equipment. In addition to the long-term changes, superimposed orbital variations of about ±1\(^{\circ }\) may be observed in response to the varying Sun illumination. An overview of ICDB corrections for the temperature range of interest is provided in Table 2. Application of the ICDB corrections greatly reduces the inconsistencies in the IF and HMW combinations for the two signal groups, but does not provide a full correction. Based on an analysis of flight data from the first month, additional empirical corrections of 0.22 m and −0.12 m are therefore applied to the L1 and L2 P(Y) code observations, respectively. It is presently unclear whether these additional biases relate to imperfections of the pre-flight receiver calibration, incompatible temperature data, or biases outside the actual receiver, i.e., cables, low noise amplifier, and antenna. Further monitoring and adjustment of the empirical biases throughout the mission is therefore suggested to facilitate Sentinel-6A POD based on PODRIX GNSS measurements.

Concerning carrier phase observations, L1 phase measurements in the PODRIX receiver are consistently obtained from C/A-code tracking in the AGGA-4 correlator. For L2, in contrast, phase measurements are based on either L2C tracking in the AGGA-4 or semi-codeless P(Y) tracking in the AGGA-2 chip. Both signals are transmitted in phase-quadrature by the modernized GPS satellites and a −1/4 cy adjustment of the L2C phase is required for alignment to the P(Y) phase reference in RINEX observations files (Romero 2020). However, additional fractional phase biases may be observed in the PODRIX L2 phase measurements which can most likely be attributed to different delays in the internal clock signals of the two chipsets. While such biases do not affect a POD processing based on float-ambiguity estimation, they inhibit a seamless use of mixed signal types in carrier phase ambiguity resolution and need to be corrected as part of the measurement generation process.

Since L2 P(Y) and L2C carrier phase observation are never obtained concurrently for a given GPS satellite, we adjusted the bias of the L2C phase measurements relative to those of P(Y) tracking in a short baseline test with the redundant receiver of Sentinel-6A while that receiver was temporarily operated in a P(Y)-only configuration on November 28, 2020. After accounting for the known antenna baseline and line-of-sight directions, double-difference observations for pairs of satellites with mixed carrier phase types can directly be used to estimate the L2C−P(Y) phase bias of the nominal receiver. While potentially affected by antenna phase pattern differences and local multipath, the resulting estimate of 0.075 cy was found to agree well with that of a laboratory test (0.073 cy) of a PODRIX engineering model and a geodetic reference receiver. Compared to an expected quarter cycle offset, the actual L2C phase bias amounts to roughly 1/3 cy and clearly needs to be taken into account for all forms of carrier phase-based positioning with ambiguity resolution.

3 Antenna calibration

The PODRIX receiver is operated with a triple-band (L1/E1, L2, L5/E5a) patch excited cup antenna (PEC; Öhgren et al. 2011), with choke ring for improved multipath suppression (Fig. 2).

Table 3 Phase center offsets from antenna reference point of Sentinel-6A antenna including manufacturer calibration and adjusted reference values for precise orbit determination
Fig. 3
figure 3

Zero-mean antenna patterns of the nominal Sentinal-6A GNSS antenna for ionosphere-free L1/L2 (top) and E1/E5a (bottom) carrier phase observations as obtained from the pre-flight chamber calibration (left) and from carrier phase residuals in the POD process (right)

Nominal values of the phase center relative to the mechanical antenna reference point (ARP) can be obtained from pre-flight anechoic chamber calibrations of the phase pattern for individual frequencies relative to a reference phase center offset (PCO) and are provided by the manufacturer as part of the ICDB. Based on these data, mean PCOs have been adjusted that minimize the deviations from a spherical wavefront over the \(+z\) hemisphere. The resulting PCOs exhibit horizontal offsets of less than 1 mm from the geometric center, but notable differences may be observed in boresight direction for the individual frequencies (Table 3). For the ionosphere-free combination of L1/L2 and E1/E5a carrier phases, the effective PCOs differ by 16 mm in z-direction. Phase variations (PVs) with respect to the mean PCOs as measured in the pre-flight calibration exhibit a high level of rotational symmetry around the boresight axis and are confined to less than ±2 mm for elevations above \(10^\circ \) (Fig. 3).

Similar to many other LEO missions with dual-frequency GNSS receivers, pronounced phase pattern distortions may, however, be observed for the GNSS antennas after integration into the spacecraft (Montenbruck et al. 2008, 2009) as a result of local multipath and other near-field effects. In case of Sentinel-6A, these may exhibit up to 10-times higher amplitudes than phase variations of the standalone antenna and necessitate a complementary in-flight calibration for precise modeling of the carrier phase observation in the POD process. Following Jäggi et al. (2009), we made use of a residual stacking approach to obtain distinct L1/L2 and E1/E5a antenna phase patterns from GPS and Galileo observations collected over a 35 days period from day of year (DOY) 2020/354 to 2021/022. The resulting patterns are shown in Fig. 3 and reveal a patchy structure with RMS amplitudes of about 4 mm and peaks up to 5 times this value over the \(80^\circ \) boresight angle cone considered for tracking by the receiver.

The patterns refer to conventional values of the PCOs listed in Table 3, which were initially adjusted to remove the mean vertical offset between kinematic and reduced-dynamic orbit determination solutions within a 1 mm tolerance. We note that the resulting calibration is implicitly tied to the employed force models and slightly different results may, for example, be obtained when using different macro-models for solar and Earth radiation pressure. While the 18 mm difference of the adjusted L1/L2 and E1/E5a z-PCOs is close to that of the pre-flight chamber calibration results, the adjusted z-PCOs themselves are roughly 25 mm smaller than suggested by the factory calibration.

Considering possible sources of this discrepancy, it may be noted that the radial component of the position vector in a dynamic or reduced-dynamic orbit determination is essentially determined by the employed force models. For near-circular LEO satellite orbits, a change of \(\Delta a_\text {R}\) in the modeled radial acceleration within the POD process causes a difference of

$$\begin{aligned} \Delta r_\text {R} = - \frac{1}{3}\left( \frac{T}{2\pi }\right) ^2 \Delta a_\text {R} \end{aligned}$$
(1)

in the estimated orbital radius, where T denotes the period of period of revolution (Montenbruck et al. 2018). For the specific case of Sentinel-6A, a 10 \(\hbox {nm}/\hbox {s}^2\) increase in the radial acceleration lowers the resulting orbit by about 3.8 mm. The kinematic position, in contrast, is determined by the orbital radius of the GNSS satellites, the assumed transmit antenna phase center offset and pattern, and, finally, the assumed offset of the receiver ARP from the LEO center of mass (CoM).

Models and standards used in the present analysis are summarized in Table 4. Even though the modeling of non-gravitational forces introduces some uncertainties in the dynamical leveling of the CoM trajectory, experience from past LEO missions and cross comparison with other POD centers shows that associated uncertainties in mean height are typically less than 5 mm (Kobel et al. 2019; Fernández et al. 2019). As such, force model deficiencies are unlikely to explain the observed \(\approx 25\) mm PCO difference between manufacturer and in-flight calibration. GPS satellite orbits and transmit antenna data for the period of interest are aligned to the IGS14b reference frame and provide consistent kinematic positions of the receiver antenna phase center in this frame. In case of Galileo, possible frame height inconsistencies of 6–7 mm may arise from the use of manufacturer-calibrated transmit antenna patterns in the generation of Galileo orbit and clock products (Villiger et al. 2020; Huang et al. 2021) but would only explain potential discrepancies in the relative location of the GPS L1/L2 and Galileo E1/E5a PCOs. However, these amount to 16 mm and 18 mm for chamber and inflight calibrations (Table 3), respectively, and are thus consistent within 2 mm.

Table 4 Models and data sets for Sentinel-6A precise orbit determination and SLR validation

4 POD concepts and models

For Sentinel-6A orbit determination, we make use of a reduced-dynamic approach (Wu et al. 1991), which combines a state-of-the-art a priori model of the spacecraft dynamics with the estimation of empirical accelerations to overcome prevailing uncertainties in the modeling of non-gravitational forces. The processing with DLR’s GNSS High precision Orbit determination Software Tools (GHOST) builds on basic concepts and algorithms for reduced-dynamic orbit determination and single-receiver ambiguity fixing as described in Montenbruck et al. (2005) and Montenbruck et al. (2018). In preparation for Sentinel-6 and other upcoming space missions with advanced GNSS receivers, the software has been extended to support multiple constellations as well as mixed signal types within an individual constellation. Among others, this includes support for signal-specific code biases as well as the capability to adjust epoch-wise inter-system biases (ISBs).

An overview of the employed models and data is given in Table 4. Satellite-specific data including GNSS antenna and LRA positions and orientations, body dimensions and surface properties, time-dependent mass and CoM locations, as well as the spacecraft attitude were provided by the manufacturer and the Sentinel-6A operations team. Auxiliary GNSS data including precise orbits and clocks, as well as code and phase biases were obtained from the IGS or its individual analysis centers. For comparison of multi-constellation processing and ambiguity resolution performance, independent POD solutions were generated using products from the Centre National d’Études Spatiales/Collecte Localisation Satellites (CNES/CLS) and the Center for Orbit Determination in Europe (CODE).

Table 5 Sentinel-6A macro-model. Adapted from Cullen (2021)

The observation model takes into account conventional GNSS transmit antenna offsets and patterns consistent with the employed orbit and clock products. In accord with established IGS standards, both CNES/CLS and CODE clock products refer to a ionosphere-free combination of L1/L2 P(Y)-code measurements for GPS and E1/E5a measurements for Galileo. While this enables direct processing of the respective PODRIX measurements, satellite-specific code biases need to be applied in the modeling of GPS L1 C/A and L2C observations. When working with CODE orbit and clock products, we made use of observable-specific biases from the associated CODE bias product, while differential code biases (DCBs) provided by the Chinese Academy of Science were used with the CNES/CLS clock products.

Single-receiver carrier phase ambiguity resolution is performed using a wide-lane/narrow-lane fixing strategy as described in Laurichesse et al. (2009) and Montenbruck et al. (2018). The processing is based on individual passes, i.e., arcs of continuous and cycle-slip-free carrier phase tracking. In a first step, an average value of the HMW combination is determined for each pass from the dual-frequency pseudorange and phase observations. Aside from measurement noise, it represents the sum of the pass-specific wide-lane ambiguity as well as the fractional receiver and satellite wide-lane biases. Taking into account the known satellite biases, integer rounding is used to separate the fractional bias, which is common to all passes, and to determine the integer-valued wide-lane ambiguities for all passes.

As part of the orbit determination a float-valued estimate \(b_i\) of the ambiguity of the ionosphere-free, dual-frequency combination is obtained for each pass i. Making use of the fixed wide-lane ambiguities, \(b_i\) can be expressed as the sum of a receiver-specific fractional phase bias, which is common to all passes, a satellite-specific fractional phase bias, and a pass-specific, integer-valued \(N_{1,i}\) ambiguity on the L1 (or E1) frequency. The integer ambiguities of all passes are then obtained through rounding of between-pass differences of the ambiguities \(b_i\) after accounting for the known satellite biases. Subsequently, the fixed wide- and narrow-lane ambiguities are introduced as soft constraints for the estimation of \(b_i\) in subsequent iterations of the orbit determination process. This is accomplished by adding the between-passes differences of the ionosphere-free carrier phase ambiguities \(b_i-b_j\) computed from the integer-valued wide- and narrow-lane ambiguities as pseudo-observations with small (0.1 mm) a priori standard deviation to the overall normal equations.

Depending on the orbit/clock product supplier, wide-lane biases for the GNSS satellites are either provided directly as an independent product (CNES/CLS) or can be computed from the observable specific code and phase biases in the generic bias product (CODE). Fractional phase biases of the GNSS satellites may be lumped into the clock estimates (CNES/CLS) or provided as part of the observable specific bias product. As discussed in Banville et al. (2020), both representations are fully equivalent and can be used inter-changeably provided that a consistent set of clock and bias products is used. In case of GPS, receiver-specific biases between L1 C/A and L1 P(Y) observations as well as L2C and L2 P(Y) observations have been removed as part of the PODRIX measurement generation to the extent possible. As such, ambiguity resolution can be performed across all satellites of the GPS constellation irrespective of the particular signals tracked.

The trajectory model for Sentinel-6A orbit determination makes use of established models for Earth gravity including tidal contributions, luni-solar gravitation and relativistic corrections in a post-Newtonian approximation. Surface forces due to solar and Earth radiation pressure (SRP, ERP) as well as drag are described through a macro-model taking into account the principal body shape and dimensions (Table 5). Optical surface properties comprise the fractions of absorbed (\(\alpha \)) as well as diffusely (\(\delta \)) and specularly (\(\rho \)) reflected photons in the visible and infrared region, which are used to compute the net radiation pressure force for each surface element as a function of the incident radiation (Milani et al. 1987). The respective values are based on manufacturer data for mylar insulation, optical surface reflectors and solar cells and take into account the share of individual materials for each of the panel (Cullen 2021). In our macro-model, diffuse re-emission of absorbed photons in both the visible and infrared region is assumed for all surfaces including the solar panels. This is different from common practice (Cerri et al. 2010), where diffuse re-emission is only considered for surfaces covered by thermal insulation layers but not assumed for solar panels. However, the SRP contribution of the large solar panels in the a priori macro-model without considering such re-emission is found to be roughly one-third too small compared to the observed accelerations. This empirical evidence suggests that the standard model is mainly applicable for solar panel “wings” where thermal re-emission from the front and backside is largely balanced. Sentinel-6, in contrast, makes use of surface mounted panels where additional thermal insulation is likely to create an asymmetric re-radiation toward the sun-lit hemisphere.

The employed macro-model represents an obvious simplification of the actual spacecraft geometry, but, nevertheless, captures the dominating variation of radiation pressure forces under the varying Sun-incidence angles. Among others, contributions of the AMR-C dish antenna with a cross section of about 1 \(\hbox {m}^2\) as well as the concave regions between the satellite body and the back-side of the outer solar panels are not considered. A coarse approximation of radiation pressure due to structures on the nadir-looking side (including the conic cover of the altimeter antenna, the DORIS antenna, and other elements) is obtained by considering two “panels” with \(\pm y\) surface normals and a 1 \(\hbox {m}^2\) cross section in the macro-model. The model simplifications are largely motivated by the computational complexity of a full-featured ray-tracing model that would be needed to account for self-shadowing or multiple reflections, but is also justified by the use of approximate optical coefficients with a presumed uncertainty of 0.1. Within the POD, these modeling deficiencies are compensated by the adjustment of piece-wise constant empirical accelerations in radial (R), along-track (T) and cross-track (N) direction over intervals of ten minutes. For Sentinel-6A, a priori values of zero and standard deviations of 4 \(\hbox {nm}/\hbox {s}^2\), 8 \(\hbox {nm}/\hbox {s}^2\), and 15 \(\hbox {nm}/\hbox {s}^2\) are adopted for the three axes. These constraints are chosen in accord with the observed amplitude of estimated accelerations, except for the radial component, which shows a notable coupling with the along-track accelerations. At the given magnitude of the a priori constraints, the estimated empirical accelerations may absorb systematic errors in the assumed antenna offsets from the CoM at a level of about 1.5 mm in radial direction (Eqn. (1)) and 15 mm in cross-track direction (Sect. 5.4).

A batch least-squares estimator is used for the Sentinel-6A orbit determination. Next to the receiver clock offsets and Galileo-GPS inter-system biases, the adjusted parameters comprise the initial state vector, force model scale factors, and empirical accelerations, as well as carrier phase ambiguities. Similar to the receiver clock offsets, ISBs are also estimated as epoch-wise parameters without a priori constraints to account for temperature related, orbital variations as well as irregular phase slips in the L1/E1, L2, and L5/E5a frontends of the PODRIX receiver. At a 30 s sampling roughly 5800 clock and ISB parameters need to be adjusted for a 24 h data arc, which drives the overall dimension of the least squares system. Using pre-elimination, the resulting normal equations can, however, be reduced to the dimension of orbit and ambiguity parameters. Considering 10 min intervals for the empirical accelerations and about 550 passes for roughly 50 GPS and Galileo satellites observed in each revolution, a total of 1000 parameters are determined in the reduced normal equations.

5 Results and validation

Making use of the methodology described above, precise orbits for the Sentinel-6A satellite were computed over a 6-month time span starting with the first transmission of science telemetry on November 27, 2020. Throughout this period, GNSS data were collected with the “nominal” (primary) PODRIX receiver, which was connected to the DORIS ultrastable oscillator (USO) and configured for “mixed,” block-specific GPS signal types.

Observation data in RINEX format were generated from the raw receiver data using code and phase bias corrections as described in Sect. 2. The receiver time scale used in the generation of the pseudorange and phase measurements was related to the instrument measurement time (IMT), which measures time based on the nominal oscillator frequency, by a third-order polynomial in IMT that minimizes the difference of receiver time and GPS time over the daily data arc. In this way, clock offsets are minimized while retaining a traceable link to the DORIS USO frequency.

As part of the measurement generation, corrections for occasional half-cycle phase slips occurring concurrently for all observations in a specific frequency band were corrected manually after initial inspection and analysis. In total, 30 such events were encountered in the 6-month data arc with an increased number of events found between end of December 2020 and mid-January 2021 as well as the second half of March 2021. All such phase slips occurred in the vicinity of the South Atlantic Anomaly (SAA), where the inner Van Allen radiation belt comes close to the surface of the Earth (Hess 1968). They may thus be interpreted as radiation-induced single-event effects (SEEs) in the receiver frontend, even though the detailed mechanism of their occurrence still needs further investigation. While uncritical for float-ambiguity POD solutions, the phase slips impact the wide- and narrow-lane ambiguity and need to be repaired for proper ambiguity resolution on the affected days.

5.1 Measurement-related quality indicators

Independent POD solutions were computed using combined GPS/Galileo orbit and clock products from the CODE and CNES/CLS analysis centers of the IGS. Internal quality checks of these solutions are provided by code and carrier phase residuals, as well as the success rate of the ambiguity solution, which provide a measure of the data quality and the overall goodness of fit. Results for the combined GPS/Galileo solutions as obtained with the two sets of GNSS products are summarized in Tables 6 and 7. The Galileo measurements clearly outperform those of GPS with RMS residuals that are roughly 30–50% smaller. Observations based on the civil C/A and L2C signals of GPS show slightly higher residuals than those of semi-codeless P(Y)-code tracking for the tracking loop bandwidth adopted in the PODRIX receiver, but can be expected to offer increased robustness and reduced sensitivity to ionospheric scintillation (van den IJssel et al. 2016).

Table 6 Residuals statistics

In comparing PODs based on CODE products with those based on CNES/CLS products, an obvious increase of the C1C/C2L code residuals level may be observed. This degradation is not related to the measurements themselves, but can be attributed to an inferior quality of CODE pseudorange biases for the L2C signal. Other than the CAS DCBs, which originate from a large L2C-capable receiver network, CODE’s observable specific biases (Villiger et al. 2019) for the L2C signals are presently derived from only a small subset of IGS stations still delivering their observations in the old RINEX v2 format. However, the reduced quality of the CODE OSBs has only a negligible effect on the wide-lane ambiguity determination from the HMW combination and does not impact the overall success rate of wide-lane and \(N_1\)-ambiguity resolution. Overall, a better than 92% fixing rate is obtained in 95% of all days over the data analysis period (Table 7).

Table 7 Ambiguity resolution statistics for wide-lane (\(N_\text {WL}\)) and \(N_1\) ambiguities
Fig. 4
figure 4

Estimated Sentinel-6A clock offsets (top) and Galileo-GPS inter-system biases (bottom) for January 4 2021 based on CNES/CLS clock products

5.2 Clock performance

The DORIS USO is a highly accurate and stable oven-controlled quartz oscillator (Debaisieux et al. 1985; Auriol and Tourain 2010; Jayles et al. 2016), which serves as the common frequency reference for the DORIS instrument as well as the main PODRIX receiver. Starting from an initial value of about \(-3.2\) ns/s, the USO drift relative to GPS time decreased to roughly \(-0.1\) ns over the first 6 months of operation.

Based on a “one-way carrier phase” (OWCP) analysis (Gonzalez and Waller 2007) of detrended 1 Hz PODRIX carrier phase observations an upper limit of the USO’s Allan deviation ranging from \(5.3 \times 10^{-12}\)s/s at time intervals of 1 s down to \(1.2 \times 10^{-13}\)s/s at 100 s can be derived. These values are governed by thermal noise of the GNSS measurements and possibly short-term clock stability of the tracked GNSS satellites. As such, they can only provide an upper limit on the actual USO performance, which is considered better by up to one order of magnitude. Nevertheless, the OWCP analysis shows that no cross talk between the internal PODRIX oscillator and the external DORIS USO can be recognized, which was previously found to cause a few-cm oscillation of the clock signal with a 1.5–3 Hz beat frequency in the GPS receiver of the Sentinel-3A mission (Jalabert and Mercier 2018).

Clock variations over time scales of up to one day are illustrated Fig. 4, which shows the offset of the detrended USO time with respect to GPS time as obtained in the POD for a selected day in January 2021. Here, periodic relativistic variations (Kouba 2019) related to the eccentricity and Earth oblateness have been removed from the apparent clock. These vary with one and two cycles per revolution, respectively, and exhibit amplitudes at the level of one decimeter for the Sentinel-6A orbit. Local peaks in the clock offset relate to SAA passages where the oscillator experiences small frequency jumps in response to ionizing radiation. In between those peaks, a largely parabolic variation of the clock offset may be noticed. As discussed in Štěpánek et al. (2020), GNSS-based observations of the radiation-induced USO frequency variation will allow to correct observations of the DORIS instrument on Sentinel-6A and obtain improved position estimates of DORIS beacons.

Concerning absolute onboard time synchronization, it may be noted that the estimated clock offsets depend on the GPS time realization in the respective GNSS clock product. Here, differences at the 0.3 m (1 ns) level may typically be observed between solutions obtained with either CNES/CLS or CODE products.

Figure 4 also illustrates the temporal variation of the Galileo-GPS inter-system biases as estimated in the POD process. While the absolute value and the overall drift of the estimated ISB depends again on the employed GNSS clock product, short-term variations are consistently reflected in all POD solutions. Aside from orbit-periodic variations attributed to temperature-related bias changes, irregular ISB jumps at the level of about 3 cm may be noticed. These have previously been identified in PODRIX ground tests and are related to phase jumps in the local oscillators for signal down-conversion in the three frontends of the receiver. Depending on the specific frequency band (L1/E1, L2, L5/E5A) affected by a jump, discontinuities will arise in the estimated receiver clock offset and/or the ISB. While these glitches put a natural limit to the observability of the DORIS frequency stability, they do not affect the DORIS operation itself in the coupled GNSS/DORIS configuration.

5.3 Orbit consistency

For an assessment of the precision, i.e., repeatability, of the orbit determination solutions we assess the consistency of different Sentinel POD solutions. This includes the comparison of results based on different GNSS orbit and clock products, the comparison of single-constellation solutions with combined GPS/Galileo solutions, and, finally, the comparison of kinematic and reduced-dynamic orbits.

Fig. 5
figure 5

Comparison of reduced-dynamic orbits based on CNES/CLS products with those based on CODE products for a sample 3-day arc

Table 8 Performance of reduced-dynamic and kinematic Sentinel-6A orbit determination solutions relative to a reduced-dynamic POD using GPS and Galileo observations and CODE GNSS products. For each axis, the mean ± standard deviation over a 6-month period are provided. All values in [mm]

For illustration, Fig. 5 shows the difference of reduced-dynamic Sentinel-6A orbits based on GPS/Galileo observations and orbit and clock products from CNES/CLS with those based on CODE products for a representative three-day interval. 3D RMS position difference over this interval amount to 1 cm. The errors are dominated by repeated peaks with typical duration of about one hour and amplitudes at the 5 cm level, which can be attributed to inconsistent ambiguity fixing for individual passes in the two solutions. Long-term statistics over the a 6-month test period show a matching 3D RMS error and confirm the overall error distribution (Table 8). Overall, 3D RMS differences between the two solutions exceed a 1 cm threshold in 4.3% and a 2 cm threshold in 1.2% of all epochs. These values appear generally compatible with the ambiguity fixing rates reported in Table 7, but further investigations may be required to obtain an optimum trade-off between fixing rates and the occurrence of false ambiguity fixes.

The summary statistics Table 8 also show a good consistency of reduced-dynamic orbit solutions for the individual constellations with the dual-constellation processing. GPS-only and Galileo-only results are essentially unbiased and only show a slightly larger scatter, which may be related to the lower number of processed satellites and observations.

Kinematic solutions obtained with previously fixed carrier phase ambiguities from the reduced-dynamic processing exhibit a 1.3 cm (3D RMS) consistency with the CODE-based reference solution. Errors are most pronounced in the radial direction which simply reflects the higher dilution of precision (DOP) in vertical direction as compared to the horizontal DOP (Langley et al. 2017). Closer inspection of the temporal evolution shows that the radial difference of kinematic and reduced-dynamic positions reveals an almost exponential decrease by roughly 10 mm over the first month of the Sentinel-6A mission before converging to a stable value near the start of 2021 (Fig. 6).

Fig. 6
figure 6

Variation of radial offset between kinematic and reduced-dynamic Sentinel-6A orbits over the first 6 months of the Sentinel-6A mission

In the reduced-dynamic POD, the estimated CoM trajectory is constrained in radial direction by the dynamical orbit model. The kinematic POD, in contrast yields an estimate of the antenna position from which the CoM location is obtained based on the assumed CoM offset of the antenna in the spacecraft frame. Varying differences between kinematic and reduced-dynamic orbits might be understood in terms of changes in the CoM location within the spacecraft that are not reflected by the manufacturer’s model for CoM shifts in response to fuel consumption and varying tank loads. Alternatively, a mismodeling of the radial acceleration may be considered, in which case the observed differences could be explained (see Eqn. (1)) by an exponentially decreasing radial acceleration with an initial value of about \(-25\) \(\hbox {nm}/\hbox {s}^2\) that is not considered in the POD process. Given the inbound orientation of the empirical acceleration and its variation in time, outgassing of the solar cells appears as a plausible cause of the observed phenomenon.

5.4 Empirical accelerations

As discussed in Sect. 4, the reduced-dynamic approach for Sentinel-6A POD makes use of empirical accelerations in addition to an a priori force model. The magnitude of these accelerations provides an indicator for the quality of the macro-model for the non-gravitational forces and can also hint at possible errors in the adopted antenna offsets from the center of mass.

Fig. 7
figure 7

Daily mean values of empirical acceleration in radial (R), along-track (T), and cross-track (N) direction

Daily mean values of the estimated empirical accelerations in radial, along-track, and cross-track direction are shown in Fig. 7. Comparison with the Sun’s elevation \(\beta \) above the orbital plane shows essentially no seasonal variation and thus confirms that the macro-model of Table 5 is able to capture the overall variation of solar radiation pressure under a wide range of Sun-orbit geometries and incidence angles.

During the first weeks after launch, gradual variations of the mean accelerations may be observed in all three axes before settling at near-constant values. As discussed above, these variations may be understood by the assumption of outgassing of satellite surfaces during the early orbit phase. Following the initial convergence, near-zero mean values are obtained for accelerations in the radial and along-track direction, whereas a systematic acceleration bias of about 8 \(\hbox {nm}/\hbox {s}^2\) may be recognized in cross-track direction. Following Montenbruck et al. (2018), a constant cross-track acceleration \(\Delta a_\text {N}\) causes a lateral shift

$$\begin{aligned} \Delta r_\text {N} = \left( \frac{T}{2\pi }\right) ^2 \Delta a_\text {N} \end{aligned}$$
(2)

of the orbital plane and may be required to align the modeled center-of-mass trajectory with the observed GNSS antenna, if the assumed GNSS antenna offset from the center-of-mass does not match the actual value. In case of Sentinel-6A, the estimated cross-track acceleration suggests that the assumed antenna reference point offset \(y_\text {ARP}-y_\text {CoM}= 8\,\hbox {mm}\) (Table 4) is too large by roughly 8 mm.

5.5 Satellite laser ranging

Satellite laser ranging constitutes a key technique for independent validation of the GNSS-based POD results of LEO satellites (Arnold et al. 2019). Since the early mission phase, Sentinel-6A has been regularly tracked by stations of the International Laser Ranging Service (ILRS; Pearlman et al. 2019). For the present analysis, a set of 13 high-performance ILRS stations (Arequipa, Graz, Greenbelt, Haleakala, Hartebeesthoek, Herstmonceaux, Matera, Monument Peak, Mount Stromlo, Papete, Potsdam, Yarragadee, and Zimmerwald) has been selected, which offer normal points with precision at the 5–10 mm level. The modeled station positions are based on the SLRF2014 reference frame and take into account time-dependent variations as summarized in Table 4.

The Sentinel-6A LRA is made up of a central prism surrounded by a ring of eight prisms with a 50\(^\circ \) tilt of their viewing directions. The centers of the prisms’ front surfaces are located on a common sphere of 82.55 mm around the optical reference point. For a group refractive index of 1.4853 at 532 nm and a prism vertex height of 22.5 mm, a range correction of about 49 mm is obtained for rays perpendicular to any of the prism surfaces. Range corrections for arbitrary incidence angles have been computed for this study based on the LRA design data using a nearest-prism approximation (Montenbruck and Neubert 2011), and are illustrated in Fig. 8.

Fig. 8
figure 8

Range correction of Sentinel-6 LRA relative to the optical reference point. The bold line at 55\(^\circ \) off-boresight angle marks the LRA field of view for stations tracking the Sentinel-A spacecraft above an elevation of 10\(^\circ \)

Fig. 9
figure 9

Sentinal-6A satellite laser ranging residuals before (top) and after (bottom) adjustment of orbit and station offsets. The solid gray lines show the variation of the Sun’s angle above the orbit plane (\(\beta \)-angle)

Residuals of roughly 37000 normal points collected up to end of May 2021 by the 13 ILRS stations are shown in Fig. 9. Compared to the ambiguity-fixed reduced-dynamic orbits obtained with CODE GNSS products, the SLR observations exhibit differences with a mean value of 1.0 mm and a standard deviation of 9.7 mm. Virtually identical results are obtained with the CNES/CLS-based POD solution.

While these statistics already demonstrate a high overall performance of the Sentinel-6 orbit determination and compatibility with the specification of a better than 1.5 cm radial error, they are still affected by systematic, cm-level errors in the SLR site coordinates and range biases that mask the actual orbit quality. For calibrating these biases, we made complementary use of ambiguity-fixed orbits of five additional LEO satellites (Swarm-A/B/C, Sentinel-3A/B) based on the same processing standards and CODE GNSS products as the Sentinel-6A POD solutions. Given the large number of observations with good hemispherical coverage of the individual stations, the offsets can be determined in a least-squares adjustment (Arnold et al. 2019) with good confidence. In addition, systematic offsets of the tracked LEO satellite in radial, along-track, and cross-track direction may be adjusted that can be caused by errors in LRA positions, GNSS antenna coordinates, or CoM location. For jointly estimating both station-related corrections and orbit errors, we applied a zero-mean constraint for both the range biases and the up-component of site positions over all contributing stations. Overall, the adjustment almost halves the residuals (Fig. 9) and yields a standard deviation of 5.7 mm which characterizes the actual precision of the SLR normal points and the GNSS-based POD solutions. Aside from the availability of normal points of individual stations, no obvious seasonal variations can be recognized in the S6A SLR residuals. In particular, no \(\beta \)-angle dependence is discernible, which might hint at systematic deficiencies in the solar radiation pressure modeling.

From the SLR residuals analysis, mean orbit errors at the few-mm level in radial, along-track, and cross-track direction, or, equivalently, the \(-z\), \(+x\), and \(-y\) body axis can be inferred for the Sentinel-6A POD solutions (Table 9). Better than 1 mm consistency of SLR observations and GNSS-based orbit determination is obtained for the radial orbit component relative to the height system defined by the 13 ILRS stations. Aside from the specific choice of stations, this favorable result depends also on the validity of the employed non-gravitational force models, as well as the proper knowledge of CoM location, LRA mount point, and LRA range corrections. Uncertainties in these factors are likely to limit the actual orbit height accuracy to a level of about 5 mm, but further monitoring and comparison of independent POD solutions will be required to verify and consolidate this error budget.

Table 9 Systematic errors of GNSS-based Sentinel-6A orbits as inferred from SLR observations over a 3-month period

Concerning the along-track component, multiple causes of the observed 2.9 mm difference between the GNSS-based orbit determination and the SLR observations can be imagined. These include GNSS antenna phase center offsets or CoM model errors in the longitudinal spacecraft axis but also biases in the timing of GNSS carrier phase observations (Montenbruck et al. 2006), which are essentially unobservable in the POD for the specific S6A attitude. SLR timing biases can likewise contribute to the apparent SLR-GNSS difference. Following Arnold et al. (2019), systematic timing errors at the sub-microsecond level are present in various ILRS stations, which might then map into an apparent along-track error at the few-mm level. However, no unambiguous explanation can be given at this stage without independent calibration of the various error sources.

The most pronounced mismatch arises in the cross-track direction (or \(-y\) body axis), where the SLR analysis suggests a systematic \(-6\) mm error in the LRA trajectory as described by the GNSS-based POD solution. This offset is likewise seen in the ambiguity-fixed kinematic and reduced-dynamic solutions, which eliminates the modeled CoM location within the spacecraft body as a potential error source. Instead, the observed offset indicates that the true separation of the LRA and the GNSS antenna in y-direction is about 6 mm smaller than modeled. While no evidence for such a modeling error is available from the manufacturer documentation, it is instructive to note that a similar (8 mm; see Sect. 5.4) discrepancy can be inferred for the CoM offset of the GNSS antenna offset from the analysis of empirical accelerations estimated in the POD process. Overall, both SLR residuals and empirical accelerations can thus be explained by assuming a modified GNSS antenna reference point location \(y_\text {ARP}\) of about \(-6\) to \(-8\) mm instead of the nominal position \(y_\text {ARP}=0\,\hbox {mm}\).

6 Summary and conclusions

Sentinel-6A is the first science mission in low Earth orbit using a dual-constellation, GPS-Galileo receiver for precise orbit determination. Based on data from the first 6 months of the Sentinel-6A mission, the achievable orbit precision and accuracy were assessed. Compared to GPS tracking, Galileo offers reduced measurement errors and yields an improved success rate in single-receiver ambiguity fixing. On the other hand, the Galileo constellation still offers a notably lower number of operational satellites placed in only three rather than six orbital planes. Despite these differences, reduced-dynamic POD solutions using only GPS or Galileo observations exhibit very good (better than \({1}\,\mathrm{cm}\) 3D RMS) consistency with the solutions based on a combined set of GPS and Galileo observations.

Processing of undifferenced GPS and Galileo observations with ambiguity fixing is presently supported by precise GNSS orbit clock, and bias products of two IGS analysis centers (CNES/CLS and CODE), which are made available on a weekly basis with 4–7 days latency. POD results obtained with the two types of GNSS products show 3D RMS differences of well below 1 cm over the entire data period. However, peak errors at the 5 cm level can be noted over representative intervals of 0.5–1 h that are attributed to erroneously fixed ambiguities for individual passes in one of the two solutions. Further adaptation of acceptance criteria for the wide-/narrow-lane ambiguity fixing may be required to reduce the occurrence of such false fixes at the possible expense of a reduced overall fixing rate.

For a consistent processing of all measurements, various calibrations of the Sentinel-6A GNSS receiver and antenna are required. These include a pre-mission or inflight calibration of hardware related inter-signal code and phase biases between semi-codeless tracking of the P(Y) signals and tracking of the civil L2C signal of GPS as a prerequisite for ambiguity fixing with mixed signal sets. Secondly, proper antenna phase center calibration in the frequency bands used for GPS (L1/L2) and Galileo (E1/E5a) tracking is required. In view of obvious (\(\Delta z \approx 25\,\hbox {mm}\)) discrepancies between manufacturer-specified and observed phase center locations relative to the spacecraft center-of-mass, an inflight calibration is indispensable for proper modeling of carrier phase observations in the POD process. Unfortunately, GNSS can no longer provide independent height information after adjustment of the antenna coordinates, and the orbital radius determined in the POD process is ultimately tied to the dynamic orbit models. Errors in these models are mainly related to the non-gravitational forces with presumed uncertainties at the 10–20% level. In radial direction, accelerations due to solar and Earth radiation pressure as predicted by the current macro-model exhibit representative values of 50 \(\hbox {nm}/\hbox {s}^2\) and 35 \(\hbox {nm}/\hbox {s}^2\), respectively. This translates into a 9–17 \(\hbox {nm}/\hbox {s}^2\) uncertainty of the radial acceleration and finally a 3–7 mm uncertainty in the height of the resulting center-of-mass trajectory.

Independent verification of the resulting POD accuracy is achieved through comparison with satellite laser ranging where 1-cm consistency of measured normal points with the GNSS-based orbits is achieved and a mean height error of less than 2 mm is identified. Following adjustment of station-specific position corrections and range biases as well as average orbit/sensor offsets, a 6 mm standard deviation of the SLR normal point residuals is achieved, which further evidences an excellent precision of the POD results.