China’s BeiDou satellite navigation system (BDS) has completed its full constellation in orbit since June 2020. Services have been evolved from regional (BDS-2) to global (BDS-3). This contribution evaluates the impact of solar radiation pressure (SRP) modeling on satellite orbits and geodetic parameters. To that end, we process 2 years of BDS observations (2019-2021), collected by a network of 100 ground stations. A physical a priori box-wing (bw) model based on the estimated optical properties is introduced. Various physical effects, such as yaw bias, self-shadowing, radiator emission and thermal radiation of solar panels are considered. The ECOM (Empirical CODE orbit Model, 5 parameters), ECOM+along-track and ECOM2 (both 7 and 9 parameters) models are employed on top of the a priori box-wing model in the experiment. We show that without the use of the a priori box-wing model, the ECOM+along-track model shows clear better orbit solutions during eclipse seasons for BDS-3 satellites. This is proven to be mainly due to the thermal radiation of solar panels. However, the along-track acceleration is highly correlated with LOD (length of day) and ECOM parameters. LOD estimates in this case are contaminated. The STD (standard deviation) of daily LOD estimates with respect to IERS-C04-14 series increases from 40 us (ECOM) to 85 us (ECOM+along-track). After the consideration of the a priori boxwing model, satellite orbital errors are greatly reduced for all the ECOM models. For instance, orbit misclosures of BDS-3 CAST (China Academy of Space Technology) satellites improve by a factor of two for the ECOM model during eclipse seasons; dependencies of SLR (satellite laser ranging) residuals on the sun elongation angle almost vanish for BDS-3 satellites. Furthermore, the use of the a priori box-wing model mitigates a great majority of the spurious signals in the geodetic parameters. In particular, the total amplitude of the 1, 3, 5, 7 cpy signals for the geocenter Z component has been reduced by a factor of 4.5 for the ECOM model. In general, the combination of the introduced physical a priori box-wing model and the ECOM model is preferred for BDS satellites.

Different Earth orientation parameter (EOP) time series are publicly available that typically arise from the combination of individual space geodetic technique solutions. The applied processing strategies and choices lead to systematically differing signal and noise characteristics particularly at the shortest periods between 2 and 8 days. We investigate the consequences of typical choices by introducing new experimental EOP solutions obtained from combinations at either normal equation level processed by DGFI-TUM and BKG, or observation level processed by ESA. All those experiments contribute to an effort initiated by ESA to develop an independent capacity for routine EOP processing and prediction in Europe. Results are benchmarked against geophysical model-based effective angular momentum functions processed by ESMGFZ. We find, that a multi-technique combination at normal equation level that explicitly aligns a priori station coordinates to the ITRF2014 frequently outperforms the current IERS standard solution 14C04. A multi-GNSS-only solution already provides very competitive accuracies for the equatorial components. Quite similar results are also obtained from a short combination at observation level experiment using multi-GNSS solutions and SLR from Sentinel-3A and -3B to realize space links. For ΔUT1, however, VLBI information is known to be critically important so that experiments combining only GNSS and possibly SLR at observation level perform worse than combinations of all techniques at normal equation level. The low noise floor and smooth spectra obtained from the multi-GNSS solution nevertheless illustrates the potential of this most rigorous combination approach so that further efforts to include in particular VLBI are strongly recommended.