Estimating surface optical properties and thermal thrust for Galileo satellite body and solar panels

Precise orbit determination of GNSS (Global Navigation Satellite System) satellites requires accurate models of perturbing forces acting on the spacecraft, i.e., solar radiation pressure (SRP) and thermal radiation forces. With the officially published satellite metadata, the analytical box-wing model is usually used to describe most of the SRP accelerations and the rest is assumed to be compensated by estimating ECOM/ECOM2 (Empirical CODE Orbit Model) parameters. However, we find that the precision of Galileo satellite orbits shows notable degradation inside eclipse seasons for 3-day-arc solutions and 24-h predictions. For instance, the RMS (root-mean-square) of orbit misclosures increases by about a factor of two in the eclipse season when using the box-wing model as the a priori and the 5-parameter ECOM model on top. The reason is proven to be mostly due to ignoring imbalanced thermal radiation forces (i.e., radiator emission and thermal radiation of solar panels) as satellite thermal properties are unknown. These imbalanced thermal effects cannot be fully absorbed by the ECOM/ECOM2 parameters inside eclipse seasons because the earth’s shadowing of a satellite in orbit causes periodic changes of the thermal environment. To cope with this problem, we first estimate satellite optical and thermal parameters as part of orbit determination based on Galileo tracking data covering 1 year. Then, we add physical thermal radiation models for radiators and solar panels as part of the a priori model and evaluate the performance of different ECOM models in Galileo satellite orbit determination. As shown by orbit misclosures, 24-h orbit predictions and SLR (Satellite Laser Ranging) residuals, the 7-parameter ECOM2 model performs better than the 5-parameter ECOM and the 9-parameter ECOM2 model for Galileo satellites. When using the 7-parameter ECOM2 model on top, the impact of the radiator emission and the thermal radiation of solar panels on Galileo satellite orbits is about 1 and 2 cm, respectively, inside eclipse seasons for 3-day-arc solutions.


Introduction
Galileo is the European Global Satellite Navigation System (GNSS), providing precise positioning, navigation and timing (PNT) services. In 2005 and 2006, the first two Galileo In-Orbit Validation Element (GIOVE) satellites were launched to secure the Galileo frequencies and to flight-test some onboard satellite components (Montenbruck et al. 2006;Steigenberger et al. 2011and 2012 In-Orbit Validation (IOV) satellites were launched for the assessment of the navigation payload and the infrastructure on ground. The first pair of Full Operational Capability (FOC) satellites were launched in August 2014 and the initial service capability of Galileo satellites was announced in December 2016. Currently, the operational Galileo constellation consists of 4 IOV and 22 FOC satellites. On 5 December 2021, another two Galileo satellites of the new generation were launched to reinforce Galileo satellite PNT services further.
Compared to other GNSS satellites (e.g., GPS), Galileo satellites have relatively high area-to-mass ratio, making the spacecraft more sensitive to non-gravitational forces. Also, satellite body of Galileo satellites is elongated in shape, SRP accelerations thus show a periodic pattern due to the rotation of the satellite body. Therefore, the radiation forces of Galileo satellites have to be more carefully considered in satellite orbit determination. For satellites with very simple spacecraft structures, Milani et al. (1987) formulate the analytical SRP model for each satellite surface based on satellite attitude law, dimensions, total mass and optical properties. For satellites with a high level of complexity in the spacecraft structure, the precise ray-tracing model is preferred if the detailed CAD model of the satellite is known (Li et al. 2018;Ziebart 2001;Ziebart and Dare 2001). In the absence of a perfect analytical model, the ECOM or hybrid ECOM model is widely used (Arnold et al. 2015;Beutler et al. 1994;Tang et al. 2021;Tseng 2021) in GNSS satellite orbit determination. Prange et al. (2017) confirm that the ECOM2 model performs much better than the 5-parameter ECOM model (Springer et al. 1999) for Galileo satellites. The reason is that higher order Fourier coefficients in the ECOM2 D direction (pointing toward the Sun) can partially absorb periodic SRP accelerations caused by the satellite body attitude motion. When combining the ECOM model with an a priori SRP model, (Montenbruck et al. 2015b;Steigenberger and Montenbruck 2017;Steigenberger et al. 2015) show that Galileo satellite orbits are clearly improved compared to the ECOM-only results.
In 2017, Galileo satellite metadata including satellite attitude law, dimensions, total mass and surface optical properties was published by the European Global Navigation Satellite Systems Agency (GSA, https:// www. gsc-europa. eu/). With this metadata information, an analytical boxwing model can be formed to compensate for most of the SRP effect (Li et al. 2019). Bury et al. (2020) show that the overall accuracy of the Galileo-FOC orbits at the level of 23 mm can be obtained when considering the box-wing model as the a priori model. However, Galileo satellite orbits still show notable degradation inside eclipse seasons. This is actually the same as we have seen for other constellations, for instance for GLONASS and BeiDou satellites (Duan et al. 2020(Duan et al. , 2021a. The reason is proven to be mostly due to ignoring the radiator emission and the thermal radiation of the solar panels as these perturbations also contribute inside the shadow, while in the shadow all the ECOM parameters are deactivated (Duan et al. 2021a). Also, thermal accelerations caused by radiators and solar panels cannot be fully absorbed by ECOM parameters.
The impact of radiator emission for Galileo satellite orbits is assessed by Sidorov et al. (2020). The emission power of radiators in the − X surface is given by accumulating the power consumption of all the clocks onboard a Galileo satellite. They show that the orbit accuracy in the radial component during eclipse seasons is improved by 14% by modeling additionally the radiator thrust. Apart from the radiator thrust, Bhattarai et al. (2022) use precise ray-tracing SRP models for Galileo satellites and set up specific models to account for the SRP force caused by LRA (Laser Retroreflector Array) surfaces and thermal force for navigation antennas. The accuracy of Galileo satellite orbits is further improved to better than 2 cm in the radial component. Thermal radiation of solar panels is assessed by Vigue et al. (1994) for GPS Block II satellites based on the thermal properties of solar panels, i.e., thickness, density, heat capacity, conductivity and emissivity. The force caused by imbalanced thermal radiation of the solar panels (front and back sides) for GPS Block II satellites produces accelerations on the order of 1.0 nm/s 2 and may cause orbit errors of more than 10 m after 7 days.
In absence of officially published thermal properties of satellite surfaces, we present physical thermal radiation force models for Galileo satellites based on the estimated thermal parameters, such as radiator emission power and scaling factor of solar panel thermal thrust acceleration. We first determine satellite optical and thermal parameters using Galileo tracking data covering 1 year. Satellite attitude errors, such as solar sensor bias (causing yaw bias) and solar panel rotation lag are considered as additional parameters in the adjustment. Then, we show the impact of radiator emission and thermal radiation of solar panels on Galileo satellite orbits. Finally, we assess the performance of different ECOM models in Galileo satellite orbit determination when combining with the analytical a priori models.

Analytical radiation and thermal models
The a priori analytical radiation and thermal forces considered in this contribution consist of SRP, earth radiation (ERP), radiator emission and imbalanced thermal force of the solar panels. Milani et al. (1987) formulate the physical SRP acceleration for solar panels and satellite body surfaces (covered by Multilayer Insulation, MLI) separately. By introducing a thermal reradiation factor , we combine these two equations for a single surface into one equation  where SRP denotes the SRP acceleration, A the surface area, M the total mass of the satellite, S 0 the solar flux constant, c the speed of light in vacuum, the thermal reradiation factor (0 for solar panels and 1 for satellite body surfaces covered with MLI). , , represent the fractions of absorbed, diffusely scattered and specularly reflected photons. D and N denote the incident unit vector and the surface normal vector, the angle between both vectors. Earth radiation, a combination of albedo and infrared radiation is modeled using the same equation. The amount of earth-reflected visible solar radiation and thermal radiation is calculated from the CERES (Clouds and Earth's Radiant Energy System) (Priestley et al. 2011). A default set of (1) satellite infrared properties is used as given by Rodriguez-Solano et al. (2012b). In order to dump excess heat generated by electronic devices, e.g., clocks inside the satellite body, radiators are used onboard Galileo satellites. As published by GSA, optical surface radiators are located on − X and ± Y surfaces for IOV satellites and on − X , ± Y and − Z surfaces for FOC satellites (https:// www. gsc-europa. eu/). The Galileo satellite body-fixed reference frame described in this study follows the IGS (International GNSS Service) definition (Montenbruck et al. 2015a). As a Galileo satellite keeps transmitting signals continuously, we assume that heat generated by the equipment is constant over time. The radiator emission thrust can be modeled as where R denotes the radiator acceleration vector and R the corresponding radiator acceleration, which is unknown for Galileo satellites. We will estimate these radiator accelerations for X and Y surfaces of Galileo satellites based on tracking data covering 1 year, as described in the next section.
Radiative heat transfer between a satellite and its environment may result in imbalanced thermal forces. As shown by Cook (1989), the dominant source for imbalanced thermal force on a GPS-like satellite is the solar panels, because of their large exposed area and low heat capacity. The density of the earth's atmosphere in space at GNSS's altitude is extremely low, the only transfer mechanism between solar panels and the environment is thus thermal radiation (Walter 2018). According to Stefan Boltzmann's law the satellite solar panels emit energy at a rate that is proportional to the fourth power of the temperature and to the emissivity on both sides of the solar panels. The temperature variation of solar panels may be formulated by balancing the input energy flux from the solar radiation and the thermal reradiation on both sides of the array where C A is the heat capacity per unit area depending on the thermal properties of solar array material, i.e., thickness, density, and specific heat (Vigue et al. 1994). represents the Sun illumination factor, the electric efficiency of the solar panel, the fractions of absorbed photons, f and b the emissivity on the front and back sides of the solar panels, the Stefan-Boltzmann constant. models the temperature difference of the front and back sides of the solar panel from thermal conductivity through the panel.
The solution to this differential equation requires the specification of boundary conditions. In our computation, the boundary conditions are defined such that the temperature of the solar panels is the same after one satellite revolution. With the given temperature the imbalanced thermal acceleration can then be computed from the difference of the thermal thrust on both sides of the panel.
For Galileo satellites, we assume a total thickness of the solar cells of 2 mm, an electric efficiency of 12%, a temperature difference of 0.91 (corresponding to a temperature difference of 0.02 T), and an emissivity of 0.75 on the front side and 0.90 on the back side of the solar panels. Solar cell layers and thermal properties are assumed to be similar as for GPS Block IIR satellites (Fliegel and Gallini 1996). Figure 1 shows the thermal thrust acceleration caused by thermal radiation of the solar panels. The acceleration is constant over time if the satellite is continuously illuminated, whereas it shows periodic variations when the satellite crosses the earth's shadow. As the assumptions of the thermal properties for Galileo solar panels are not accurate enough, a scaling factor S will be estimated to compensate for modeling errors, as described in the next section.

Parameter set up
To set up precise physical radiation and thermal models, we developed a method to determine satellite optical and thermal parameters as part of orbit determination (Duan  (1)) are estimated. The respective partial derivatives can be written as In the adjustment, surface normal vector ( ), dimension ( A ) and total mass ( M ) of the satellite are fixed to the published values. As described by Rodriguez-Solano et al. (2012a), solar panels may follow the Sun with a small lag ( lag ). We estimate this lag parameter whose partial derivative can be written as where ̇ denotes the rate of the elongation angle, the angle between the normal vector of the solar panels and the Z-axis, GNSS satellite attitude is partly controlled by the data from solar sensors. A sensor bias may cause a yaw bias which is a function of the elongation angle where sbias denotes the sensor bias, the elongation angle, ybias the yaw bias. By rotating the satellite body by the ybias angle around the Z-axis, X and Y change accordingly. From the SRP differences caused by a ybias using (1), the numerical partial derivative can be computed to estimate the sensor bias parameter.
As radiator information in (2) for Galileo satellites is not publicly known, we estimate constant radiator accelerations in − X and + Y surfaces. A radiator acceleration in the − Z-axis is not estimated because it is fully correlated with the semi-major axis parameter in satellite orbit determination. The reason is the insensitivity of pseudorange observations on the orbit height. The constant thermal radiation force of solar panels outside eclipses is compensated by estimating the optical parameter of the solar panels. The periodic thermal accelerations (cooling down and heating up parts) are corrected by estimating a scaling factor S . For the sake of simple usage, we also provide the function coefficients to direct compute the cooling down (analytical solution) and heating up (approximate solution) accelerations without the need for a deep understanding of the physical background. The analytical solution of the acceleration during cooling down is where acc shd denotes the thermal acceleration inside the shadow, p and c are function coefficients, t the current time epoch and t 0 the shadow entering epoch. S represents the scaling factor that compensates for model errors caused by assumptions of solar panel thermal properties. The empirical exponential function for the acceleration during heating up is where acc sun denotes the thermal acceleration outside the shadow, t h and h are function coefficients. t sft represents a time shift assuring optimal fit to the solution of the differential equation after shadow exit where t shd denotes the shadow duration in seconds. All the function coefficients are given in Table 1. The scaling factor S is given in Tables 3 and 4 together with other adjusted parameters. In summary, there are six optical parameters for satellite body surfaces, one optical parameters for solar panels, one solar panel rotation lag, one sensor bias, two radiator accelerations, and one scaling factor of the solar panel imbalanced thermal acceleration, i.e., in total 12 parameters for each satellite.

Parameter estimation
We estimate all the parameters using 90 globally distributed IGS MGEX ) tracking stations. The time span starts from day of year (doy) 275, 2020 to doy 274, 2021. In the first step, we compute Galileo satellite orbits using the a priori box-wing model (based on the published metadata) and the 5-parameter ECOM model on top. The geodetic datum is defined by a minimum constraint (no net rotation condition) of all the stations to the IGS14 reference frame. Earth rotation parameters are fixed to the values of the IERS Bulletin A. Earth radiation and satellite antenna thrust (Steigenberger et al. 2018) effects are modeled. Satellite antenna phase center offsets and variations (PCOs, PCVs) are corrected using the igs14.atx file. Station PCO and PCV corrections for Galileo E1 and E5a frequencies are assumed to be the same as those for GPS L1 and L2 frequencies. Detailed settings are given in Table 2.
In the second step, instead of estimating ECOM parameters, we estimate 12 radiation force parameters together with the 6 keplerian elements for each satellite in the orbit determination procedure. Integer ambiguities and stationrelated parameters that were determined in the first step are fixed as known. In a daily arc solution, the 6 orbital parameters and 12 radiation force parameters of one satellite are highly correlated, it is impossible to obtain reasonable daily solutions. To reduce high correlations between parameters, we pre-eliminate orbital parameters from daily normal equations leaving only radiation force parameters, and we stack daily normal equations over 1 year (corresponding to two times the full beta angle range) to compute reasonable radiation force parameters. Because, on the one hand, radiation force parameters do not change within 1 year (ignoring aging effects of materials), on the other hand, we have proven that correlations between parameters are greatly reduced when including full beta angle range measurements in the adjustment . Tables 3 and 4 show the estimated radiation force parameters and the published optical properties of Galileo IOV and FOC satellites. Formal errors of all the estimated optical properties are smaller than 0.02. Note that the published optical properties of different materials in the same satellite body surface are combined into one set of values according to the proportion of the area of each element. The estimated optical properties show some differences of about 5-20% compared to the published values, i.e., in the − Z surface of IOV satellites and in the + X surface of FOC satellites. The reason could be partially due to the correlation between parameters, in particular between optical parameters and radiator parameters. However, it is also possible that optical properties change with the long-term exposure in space. Galileo IOV satellites show a radiator acceleration of − 1.1 nm/ s 2 in the − X direction, corresponding to an emission power of about 360 W. The radiator effect in the Y surface for IOV satellites either does not exist or is balanced by both sides. The radiator acceleration for FOC satellites is − 0.9 nm/s 2 and 0.68 nm/s 2 in − X and + Y direction, respectively, corresponding to a total emission power of about 500 W. The solar panel ρ estimates represent the constant thermal radiation of solar panels outside eclipse seasons. The scaling factor (S) estimates shows that the computed thermal acceleration in Fig. 1 is not accurate and has to be compensated by  the scaling factor. Galileo satellites seem to have no solar sensor bias according to our estimation but the solar panel rotation lag exists with a value smaller than 1 deg.

Galileo satellite orbit determination using different radiation force models
We use all the estimated radiation force parameters in a physical a priori model where apr denotes the total a priori acceleration, SRP the solar radiation pressure acceleration, sbias the ybias effect, lag the solar panel rotation lag effect, ERP the earth radiation acceleration, R the radiator acceleration, sun and shd are thermal thrust accelerations of the solar panels. We use ECOM and ECOM2 models on top to assess the performance of different SRP models in Galileo satellite orbit determination. All the modeling options are given in Table 5. We use the same tracking stations and time period (1 year) as in the last section. Orbit arc length is 3 days and we extract the middle day solution as the final daily solution. Pseudo-stochastic pulses are not considered. Orbit misclosures between two consecutive arcs, 24-h orbit predictions and SLR residuals are analyzed to assess the quality of the determined Galileo satellite orbits.
We compute orbit differences at the day boundaries between consecutive arcs and define the orbit misclosure as where i and i + 1 refer to days, s,i+1 denotes the orbit position vector of satellite s of day i + 1 (Lutz et al. 2016). Figure 2 illustrates the 3D (position distance) mean RMS of orbit misclosures for each modeling option. The use of the gsa model based on the published metadata shows reasonable improvement for all the ECOM models compared to those without any a priori model. For instance, the improvement of gsa + 7ECOM2 compared to 7ECOM2-only is about 40% for FOC satellites outside the eclipse season. However, the ECOM-only solution inside the eclipse season is about two times worse than the ECOM2-only solution, and the use of the gsa model does not help much. When having considered thermal force models in est, the ECOM model solution inside eclipse seasons is improved by a factor of two. Thus, the imbalanced thermal forces are not negligible in satellite orbit determination, especially during eclipse seasons. The mean RMS values of all the radiation pressure models are shown in Table 6. The est + 7ECOM2 model shows the best RMS value of about 2.0 and 3.5 cm for Galileo satellites outside and inside eclipse seasons.
Satellite orbit predictions are usually used in real-time applications, the accuracy of the predicted orbits depends tightly on the radiation force models . As a further step to assess the performance of different radiation force models, we predict Galileo satellite orbits over 24 h based on the 3-day-arc solution. The time periods cover the same 1 year. All the predicted Galileo satellite orbits are compared to the est + 7ECOM2 daily solutions of the same epoch. Figures 3 and 4  and 10% compared to the use of the gsa model for ECOM, 7ECOM2 and ECOM2 model. The est + 7ECOM2 model provides the best orbit predictions. RMS values of IOV satellites are 2.0, 7.8 and 2.5 cm in radial, along-track and cross-track directions outside the eclipse season, and are 3.4, 12.2 and 4.9 cm inside the eclipse season. RMS values of FOC satellites are slightly better, with 1.1, 5.0 and 2.1 cm in radial, along-track and cross-track directions outside the eclipse season, and 3.9, 14.5 and 2.9 cm inside the eclipse season.

show mean RMS values
All the Galileo satellites are equipped with laser retroreflector arrays (LRAs). Satellite laser ranging (SLR) measurements collected by the International Laser Ranging Service (ILRS) can be used as external reference to evaluate satellite orbit results (Pearlman and Degnan 2002;Pearlman et al. 2019). Table 7 shows the mean offset and STD (standard deviation) of SLR residuals for different radiation pressure models. SLR stations that show systematic effects are excluded (Bury et al. 2021). SLR residuals exceeding ± 25 cm are not used for the statistic. SLR residuals from ECOM2 model are about two times better than those from ECOM and 7ECOM2 models without the use of any a priori model. The use of the a priori gsa or est model results in a clear improvement by combining with any ECOM model on top. For instance, the STD value is reduced by about 60% for the ECOM and 7ECOM2 model, and by about 25% for the ECOM2 model. However, the use of the gsa model increases the mean residual value from less than 1 cm to about 3 cm for FOC satellites, whereas, it decreases to almost zero when using the est model. In general, est + 7ECOM2 model shows the best SLR residuals.

Impact of individual radiation forces on satellite orbits
In order to show the impact of each term in (12) on Galileo satellite orbits, we take the est + 7ECOM2 solution as the reference and evaluate satellite orbits that are determined by ignoring one of the radiation forces in (12). All the settings are the same as in the previous sections. The ERP effect is confirmed by Rodriguez-Solano et al. (2012b) to have an impact of about 1-2 cm in the radial direction, and therefore, we do not repeat the analysis. The yaw bias angle for Galileo satellites is very close to zero and does not need to be considered. The importance of the a priori SRP model for Galileo satellites is clear, we can already draw conclusions from our aforementioned results. The interesting point is whether our estimated Galileo satellite optical properties perform better than the published values in Galileo satellite orbit determination. So, in this section we first assess both sets of optical parameters in satellite orbit determination, considering no other physical effects but only the a priori SRP box-wing model. Table 8 shows RMS values of orbit misclosures and statistics of SLR residuals using the published (GSA) and the estimated surface optical parameters. The a priori box-wing model and the 7ECOM2 model are used, thermal radiation forces and other physical effects are not considered. As seen from orbit misclosures. a minor improvement (1 mm) is observed when using the estimated optical values. The STD values of SLR residuals are almost the same when using both parameters, but an improvement of the mean offset for FOC satellites is clearly visible when using our estimated optical parameters. Therefore, the estimated optical parameters of Galileo satellites are preferred.
Then, we show the impact of solar panel rotation lag, radiator emission and thermal radiation of the solar panels on Galileo satellite orbits. Figure 5 shows the daily RMS of orbit differences by ignoring solar panel rotation lag, radiator emission and thermal radiation of solar panels separately (Galileo E03 satellite is taken as an example). All the solutions employ the 7ECOM2 model on top and are compared  to the reference solution (est + 7ECOM2). Solar panel rotation lag and thermal radiation of solar panels are correlated and have a notable impact of 1-3 cm in the eclipse season. In the non-eclipse season, the impact of these two perturbations is close to zero because the lag and thermal effects of solar panels in this case can be well absorbed by the ECOM BS (sine term in the B direction) and D0 (constant term in the D direction) parameter, respectively. We can also find that when the absolute value of the β angle (Sun elevation above the orbital plane) is larger than 60 deg the mean impact of solar panel rotation lag and thermal radiation of solar panels is about 0.5 cm. This is because one of the other two orbital planes during this time period is in the eclipse season, satellite orbit solutions of the E03 orbital plane are affected through other common parameters, i.e., receiver clock offsets and station coordinates. The radiator thrust shows an overall impact of about 0.5 and 1.0 cm on satellite orbits outside and inside the eclipse season. The mean RMS of all satellites for all these three perturbations are given in Table 9.

Summary and outlook
This contribution aims to set up precise thermal radiation force models for Galileo satellites. In the absence of published thermal properties, we estimate thermal modeling parameters together with optical parameters and attitude bias parameters as part of orbit determination. The current result shows that the estimated surface optical properties of Galileo satellites have differences of about 5-20% in different satellite surfaces (+ X and ± Z) compared to the published values.
These differences cause 1 mm differences in the RMS of orbit misclosures and in the STD of SLR residuals, but cause about 3 cm differences in the mean offset of SLR residuals for the FOC satellites. All the Galileo satellites have no yaw bias according to our estimation, but have a rotation lag of about 0.7 and 0.5 deg for IOV and FOC satellite, respectively. Without any yaw bias, the Y-bias estimates of the FOC satellites are dominated by the radiator thrust in the Y surface. The emission power of radiators on the − X surface for IOV satellites is about 360 W, while the total emission power of radiators on the − X and Y surfaces for FOC satellites is about 500 W. Thermal radiation of solar panels is first calibrated by using a physical thermal radiation model based on some assumed material compositions and thermal properties. Then, thermal radiation modeling errors are compensated by estimating a scaling factor using tracking data covering 1 year. With the box-wing SRP model and all the thermal radiation models, it is evident that Galileo satellite orbits are further improved compared to solutions using only the boxwing model and ECOM models (ECOM, 7ECOM2 and ECOM2). In particular, the improvement for the 5-parameter ECOM model is more than a factor of two in the eclipse season. The est + 7ECOM2 model shows the best orbit misclosure of about 2.0 and 3.5 cm outside and inside eclipse seasons. When predicting orbit solutions over 24 h the use of the a priori est model shows an improvement of about 30%, 15% and 10% over the gsa model for ECOM, 7ECOM2 and ECOM2 model, respectively. However, we observe minor improvement in SLR residuals (STD value) by considering thermal radiation effects since imbalanced thermal perturbations affect satellite orbits mainly in the along-and crosstrack directions.
To assess the impact of solar panel rotation lag, radiator effects and thermal radiation of solar panels, we use the 7ECOM2 model on top and ignore each of these perturbations one by one. Radiator effects of Galileo satellites have an impact of about 0.5 cm and 1.0 cm on satellite orbits outside and inside eclipse seasons. Solar panel rotation lag and thermal radiation of solar panels have minor impact on satellite orbit outside eclipses, because the lag and thermal effects of solar panels can be very well absorbed by the ECOM BS and D0 parameter, respectively. Although we  compute the impact of lag and thermal effects of solar panels inside the eclipse season individually, correlations between these two effects cannot be avoided. It makes more sense to summarize the total impact of both perturbations, which is 3.2 cm for both types of Galileo satellites inside the eclipse season if the 7ECOM2 is used on top.
Acknowledgments This research is based on the analysis of Galileo measurements provided by the IGS MGEX (publicly available). SLR measurements are provided by ILRS. The effort of all the agencies and organizations as well as all the data centers and analysis centers is acknowledged. Calculations are based on the Bernese GNSS Software (license available, University of Bern) and its specific modifications by the authors (not public).
Funding Open Access funding enabled and organized by Projekt DEAL.
Data availability Galileo satellite RINEX and SLR observations are downloaded from (ftp:// gdc. cddis. eosdis. nasa. gov). The Galileo satellite metadata is made available by the European Global Navigation Satellite Systems Agency (GSA, https:// www. gsc-europa. eu/).
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.