First Galileo Single-Frequency Occultation Process and Precision Analysis of FengYun3E

This article proposes a single-frequency occultation method whose core is the reconstruction of a second frequency measurement. We process the actual received Galileo E1 single-frequency occultation observation data of FengYun3E to meet the urgent need for single-frequency Galileo occultation inversion of FengYun3 E/F/G/H satellites. Galileo single-frequency occultation events are globally distributed evenly and have stable quantities. The refractive index products and dry temperature products inverted from the single-frequency occultation data are reliable at altitudes of 5–30 km. The Galileo E1 single-frequency occultation process can effectively improve the output of FengYun3E/GNOS occultation products. These results validate the feasibility and correctness of using FengYun3/GNOS for the actual Galileo single-frequency occultation process. The atmospheric occultation products of Galileo single-frequency occultation increase the quantity of global occultation products and serve as a beneficial supplement to global numerical weather prediction data sources.


Introduction
Radio occultation (RO) technology was initially developed to explore the vertical profile parameters of the ionospheres and atmospheres of the planets in the solar system [1,2]. Subsequently, Fishbach proposed using a global navigation satellite system (GNSS) RO to probe the atmospheric parameters of the Earth [3]. The United States' Global Positioning System (GPS) became more mature in the late 1980s, and GNSS RO began to develop gradually and flourish. Currently, mature GNSS RO utilizes established GNSS systems such as BDS, GPS, Galileo, and GLONASS as signal transmitters, with GNSS receivers mounted on low-Earth-orbit satellites (LEOs) as observation platforms. By receiving and analyzing signal delay caused by the bending of the signal path, atmospheric parameters such as refractive index, pressure, temperature, humidity, and other related products can be retrieved through inversion, which are of great significance in improving numerical weather models, monitoring climate change, and studying atmospheric dynamic processes [4]. GNSS RO, with the capable of self-calibration and lack of bias, provides an effective complement to other data types, such as radiosondes and weather balloons [5]. Therefore, GNSS RO products from LEOs are one of the most significant factors in reducing errors in numerical weather prediction models [6]. If GNSS RO events are uniformly distributed globally, the greater the number of GNSS RO events is, the larger the volume of raw observational data available for inclusion in numerical weather prediction, leading to of the GNOS clock bias is high, it can improve the inversion precision of the zero-difference method and achieve similar inversion precision as the single differencing method [14].
During the algorithm verification stage, we have already processed GPS singlefrequency occultation data and BDS single-frequency occultation data from FY3E/GNOS, utilizing the single-frequency occultation process of reconstructing the second frequency with the zero-difference method. We have compared the processed products with the products derived from the dual-frequency occultation process and ERA5 reanalysis to assess their precision [15,16]. Compared with the dual-frequency occultation inversion, the correlation coefficient of the relative total electron content (relTEC) for single-frequency GPS and single-frequency BDS occultation inversion is greater than 0.95. The average deviation of the excess phase Doppler for single-frequency inversion is mostly less than 0.2%. The refractive index product obtained from single-frequency inversion has slightly lower precision compared to the dual-frequency inversion product, but they show good consistency in terms of deviation and root mean square (RMS) indicators compared to ERA5 reanalysis. This verifies the feasibility of the proposed method. In this study, actual FY3E Galileo single-frequency occultation observations were used for product inversion and precision analysis, with the aim of providing more occultation data for numerical weather prediction. Section 2 of this paper describes the FY3E/GNOS Galileo E1 single-frequency occultation process method. In Section 3, the global distribution of occultation events during the operational period of FY3E/GNOS is analyzed. Subsequent Sections 3.2 and 3.3 evaluate the precision of the refractive index product and the dry temperature product inverted from the Galileo single-frequency occultation process. Section 4 summarizes the research findings, and Section 5 discusses the potential applications of the single-frequency occultation process. Figure 1 shows the Galileo E1 single-frequency occultation process flow, which includes the following steps:
FY3E/GNOS Level 0 data are segmented, converted into an appropriate format, and preprocessed to positioning observation data and occultation observation data.

2.
FY3E precise orbit determination (POD) using FY3E positioning observation data and accurate GNSS products from IGS, resulting in precise position, velocity, and receiver clock bias. 3.
Reconstruction of the carrier phase occultation observation data at the second frequency using pseudorange and carrier phase occultation observation data at Galileo E1 single-frequency.

4.
Computation of the neutral atmosphere excess phase using the results from the previous three steps. 5.
Inversion of the neutral atmosphere bending angle and refractive index using Abel integral transformation. 6.
Retrieval of atmospheric temperature based on refractive index.
FY3E/GNOS includes two sets of positioning antennas, each with a corresponding RF unit, as well as two sets of atmospheric occultation antennas, each also equipped with an RF unit, in addition to a data processing unit. The antennas utilize low noise RF front end technology to ensure effective capture and processing of GNSS signals. The central frequencies of the receiver channels are shown in Table 1. FY3E/GNOS simultaneously receives positioning signals and occultation signals: GPS L1&L2, BDS B1&B3, and Galileo E1 [17]. The GNSS RO data record rule in GNOS is that when a rise or set GNSS RO event occurs, the carrier phase and amplitude of the radio signals are recorded at a high sampling rate of 50 Hz. The sampling rate in the pseudorange is recorded at a range of 1-50 Hz [18]. FY3E/GNOS includes two sets of positioning antennas, each with a correspond RF unit, as well as two sets of atmospheric occultation antennas, each also equipped w an RF unit, in addition to a data processing unit. The antennas utilize low noise RF fr end technology to ensure effective capture and processing of GNSS signals. The cen frequencies of the receiver channels are shown in Table 1. FY3E/GNOS simultaneou receives positioning signals and occultation signals: GPS L1&L2, BDS B1&B3, and Gali E1 [17]. The GNSS RO data record rule in GNOS is that when a rise or set GNSS RO ev occurs, the carrier phase and amplitude of the radio signals are recorded at a high sa pling rate of 50 Hz. The sampling rate in the pseudorange is recorded at a range of 1 Hz [18]. Because Galileo positioning data have a single-frequency and Galileo position channels are significantly fewer than those of GPS, we utilize GPS L1&L2 positioning servation data for ground POD. The satellite position and velocity coordinates are within the framework of the Earth-centered inertial coordinate system in the FY3E PO FY3E is equipped with two positioning antennas, and the simplified dynamics metho employed by merging the observation data from the half field of view of the two anten [19]. The POD is computed at intervals of 30 s, and all GPS positioning observation d are processed in batches with a duration of 30 h. The fundamental observations used ionosphere-free pseudorange and phase combinations. State vectors, including positi velocity, and clock bias, are estimated, along with empirical segment constants for ac eration (with 90 min intervals) and atmospheric drag (with 6-h intervals). The estima results from each iteration are used as the initial orbit for the next iteration until the e  Because Galileo positioning data have a single-frequency and Galileo positioning channels are significantly fewer than those of GPS, we utilize GPS L1&L2 positioning observation data for ground POD. The satellite position and velocity coordinates are all within the framework of the Earth-centered inertial coordinate system in the FY3E POD. FY3E is equipped with two positioning antennas, and the simplified dynamics method is employed by merging the observation data from the half field of view of the two antennas [19]. The POD is computed at intervals of 30 s, and all GPS positioning observation data are processed in batches with a duration of 30 h. The fundamental observations used are ionosphere-free pseudorange and phase combinations. State vectors, including position, velocity, and clock bias, are estimated, along with empirical segment constants for acceleration (with 90 min intervals) and atmospheric drag (with 6-h intervals). The estimated results from each iteration are used as the initial orbit for the next iteration until the estimation residuals are sufficiently small.
In the radial-tangential-normal (RTN) reference coordinate system, POD achieves precision in position and velocity, using the RMS errors from overlapping arc segments as indicators. Figure 2 illustrates the RMS errors for a 14-day FY3E POD analysis. Subfigure (A) represents the positional precision, which is lower than 2.5 cm along the T direction, and the average 3D RMS position precision is 1.91 cm. Subfigure (B) represents the velocity precision, which is lower than 25 µm/s along the R direction and, the average 3D RMS velocity precision is 15.96 µm/s. The computed position and velocity precision are sufficient for the inversion of atmospheric occultation products.
Remote Sens. 2023, 15, x FOR PEER REVIEW 5 o velocity precision, which is lower than 25 µm/s along the R direction and, the average RMS velocity precision is 15.96 µm/s. The computed position and velocity precision sufficient for the inversion of atmospheric occultation products.

Reconstruct a Second Frequency Measurement
The previous section provided an overview of the characteristics of GNOS Galileo occultation observation data and the FY3E POD. This section will discuss the reconstr tion of the second frequency for the Galileo E1 single-frequency occultation process. S cifically, the reconstruction of the carrier phase measurement of the second frequenc based on the Galileo E1 single-frequency pseudorange measurement E and car phase measurement C . In theory, the choice of the second frequency is arbitrary, but the versatility of the occultation Level 2 product process, the second frequency in this per is based on the E5a in Galileo's dual-frequency measurement, with a frequency 1176.45 MHz.
The reconstructed occultation carrier phase measurement E * is greatly affected occultation pseudorange measurement noise. Therefore, a regularization method t minimizes the second derivative is commonly used to filter the noise in the occultat pseudorange measurements to reduce its impact on E * [20][21][22]. The first step is to c struct a filtering matrix F. This is done by interpolating and patching the data gaps in (sampled at a rate of 50 Hz) and C (sampled at a rate of 1-50 Hz) using a regularizat parameter γ, thereby reducing the noise influence from the observations. A value o equal to 10 is chosen for a smoothing strength that corresponds to a low-pass filter w a cutoff frequency of approximately 0.05 Hz. To unify the sampling rates of E and C 50 Hz, the size of the filtering matrix should be 500 × 500. Equation (1) represents calculation formula for the filtering matrix F in the regularization method that minimi the second derivative [13]: where I is a diagonal matrix, in which the corresponding diagonal elements are set t at the positions where data exist in both E and C , and the remaining elements are

Reconstruct a Second Frequency Measurement
The previous section provided an overview of the characteristics of GNOS Galileo E1 occultation observation data and the FY3E POD. This section will discuss the reconstruction of the second frequency for the Galileo E1 single-frequency occultation process. Specifically, the reconstruction of the carrier phase measurement of the second frequency is based on the Galileo E1 single-frequency pseudorange measurement E 1 and carrier phase measurement C 1 . In theory, the choice of the second frequency is arbitrary, but for the versatility of the occultation Level 2 product process, the second frequency in this paper is based on the E5a in Galileo's dual-frequency measurement, with a frequency of 1176.45 MHz.
The reconstructed occultation carrier phase measurement E * 5a is greatly affected by occultation pseudorange measurement noise. Therefore, a regularization method that minimizes the second derivative is commonly used to filter the noise in the occultation pseudorange measurements to reduce its impact on E * 5a [20][21][22]. The first step is to construct a filtering matrix F. This is done by interpolating and patching the data gaps in E 1 (sampled at a rate of 50 Hz) and C 1 (sampled at a rate of 1-50 Hz) using a regularization parameter γ, thereby reducing the noise influence from the observations. A value of γ equal to 10 6 is chosen for a smoothing strength that corresponds to a low-pass filter with a cutoff frequency of approximately 0.05 Hz. To unify the sampling rates of E 1 and C 1 to 50 Hz, the size of the filtering matrix should be 500 × 500. Equation (1) represents the calculation formula for the filtering matrix F in the regularization method that minimizes the second derivative [13]: where I is a diagonal matrix, in which the corresponding diagonal elements are set to 1 at the positions where data exist in both E 1 and C 1 , and the remaining elements are set to 0. Its size is 500 × 500. The matrix S is a finite difference operator for the second derivative, which is represented by S T as its transpose.
The second step is to calculate the reconstructed carrier phase measurements, denoted as E * 5a [13].
As shown in Equation (2), the E5a frequency carrier phase measurement E * 5a is calculated by subtracting the carrier phase measurement E 1 from the filtered E 1 − C 1 . Here, f 1 represents the frequency of the Galileo E1 signal, and f 2 represents the frequency of the Galileo E5a signal, both measured in Hz. Then, the term E * 5a − C 1 can be used to invert the reconstructed relTEC of the second frequency, measured in total electron content units (TECUs). It has been statistically observed that the range of relTEC falls within 20 TECUs. The relTEC variation rate is involved in the excess phase Doppler calculation [23]. Therefore, Figure 3 illustrates the reconstructed relTEC variation rate for two representative occultation events, as shown in subplots (A) and (B). In particular, the calculation of the rate has nothing to do with the type of occultation rising or setting, and is based on the tangent point height from low to high. As can be seen from the figure, the relTEC variation rate is always positive, which indicates that relTEC is always increasing within 120 km.
As shown in Equation (2), the E5a frequency carrier phase m culated by subtracting the carrier phase measurement E from the f represents the frequency of the Galileo E1 signal, and f repr the Galileo E5a signal, both measured in Hz. Then, the term E * the reconstructed relTEC of the second frequency, measured in tota (TECUs). It has been statistically observed that the range of relTEC The relTEC variation rate is involved in the excess phase Doppler fore, Figure 3 illustrates the reconstructed relTEC variation rate occultation events, as shown in subplots (A) and (B). In particula rate has nothing to do with the type of occultation rising or settin tangent point height from low to high. As can be seen from the figu rate is always positive, which indicates that relTEC is always incre

Excess Phase Calculation
The initial purpose of reconstructing the second frequency is single-frequency occultation inversion algorithm to general GNSS

Excess Phase Calculation
The initial purpose of reconstructing the second frequency is to adapt the Galileo E1 single-frequency occultation inversion algorithm to general GNSS occultation Level 2 inversion software (ROPP 6.0 is used in this paper). Figure 4 shows an illustration of a GNSS RO. FY3E (denoted by A) simultaneously receives signals from the GNSS satellite (denoted by B) undergoing an occultation event and signals from a high elevation angle GNSS satellite (denoted by C). The red path represents the bending of the signal in the atmosphere. In this case, AB represents the occultation-LEO link, and AC represents the reference-LEO link. In the single difference method, the AB observation equation is subtracted from the AC observation equation to eliminate the receiver clock bias. However, since the reference satellite C measurement is also a single frequency, the solution is affected by the noise in the reference satellite pseudorange measurement. Fortunately, the FY3E/GNOS clock has stability reaching the order of 10 −12 in one second, and the clock drift precision is sufficient for calculating the excess phase with relatively high precision using the zero-difference method.
FOR PEER REVIEW 7 of 17 since the reference satellite C measurement is also a single frequency, the solution is affected by the noise in the reference satellite pseudorange measurement. Fortunately, the FY3E/GNOS clock has stability reaching the order of 10 −12 in one second, and the clock drift precision is sufficient for calculating the excess phase with relatively high precision using the zero-difference method. The calculation of the excess phase using the zero-difference method only involves link AB and does not include link AC. This is expressed in Equation (3): In the equation, t represents the time of GNOS signal reception. i represents the frequency of either E or E * . δ , , r , E , and ∆φ t represent the excess phase, geometric distance AB, carrier phase observation of E1 or reconstructed carrier phase observation of E5a, and corresponding phase ambiguity correction, respectively. Their units are all meters. c is the speed of light, taken as a constant value of 2.99792458 × 10 8 m/s. For the time-related parameters, δ and δ represent the clock biases of satellites A and B, respectively; δ , and δ represent the relativistic effects on the clock biases of satellites A and B, respectively; and τ represents the radio wave propagation time. δρ , represents the relativistic effects, excluding clock biases. All time-related parameters are expressed in units of seconds. Information regarding the positions, velocities, and clock biases is obtained through POD. Additionally, ε represents measurement noise.
The excess phase Doppler is the input for refractive index product inversion. Therefore, it is necessary to convert the excess phase into excess phase Doppler [24]. The equation for calculating the excess phase Doppler corresponding to the E1 and E5a frequencies is as follows: The calculation of the excess phase using the zero-difference method only involves link AB and does not include link AC. This is expressed in Equation (3): In the equation, t r represents the time of GNOS signal reception. i represents the frequency of either E 1 or E * 5a . δ ex,i , r B A , E i , and ∆ϕ B A (t r ) represent the excess phase, geometric distance AB, carrier phase observation of E1 or reconstructed carrier phase observation of E5a, and corresponding phase ambiguity correction, respectively. Their units are all meters. c is the speed of light, taken as a constant value of 2.99792458 × 10 8 m/s. For the time-related parameters, δ A and δ B represent the clock biases of satellites A and B, respectively; δ A,rel and δ B rel represent the relativistic effects on the clock biases of satellites A and B, respectively; and τ B A represents the radio wave propagation time. δρ B A,rel represents the relativistic effects, excluding clock biases. All time-related parameters are expressed in units of seconds. Information regarding the positions, velocities, and clock biases is obtained through POD. Additionally, ε represents measurement noise.
The excess phase Doppler is the input for refractive index product inversion. Therefore, it is necessary to convert the excess phase into excess phase Doppler [24]. The equation for calculating the excess phase Doppler corresponding to the E1 and E5a frequencies is as follows: The equation above, where λ i represents the wavelength of frequency i, provides the excess phase Doppler calculated from the E1 frequency occultation data and the reconstructed excess phase Doppler calculated from the E5a frequency observation data. Figure 5A,B illustrate two representative occultation events. The x-axis represents the excess phase Doppler in Hz, while the y-axis represents the height in km in the form of a logarithmic axis.
ote Sens. 2023, 15, x FOR PEER REVIEW descends through the atmosphere, the refractive index becomes stronger to bend more, thus leading to a noticeable change in the excess Doppler d Figure 5. Excess phase Doppler corresponding to heights ranging from 0 to 140 km represent two representative occultation events. The excess phase Doppler calcu frequency occultation data is represented by the blue curve, while the excess pha lated from the reconstructed E5a observation data is represented by the orange cu  From the figures, during the initial stages of occultation, the excess phase Doppler is close to zero. As the height decreases, the excess Doppler delay undergoes significant changes. This is because when the ray's tangent point crosses the ionosphere and upper troposphere, where the medium is relatively thin, the signal amplitude remains relatively constant, resulting in the excess phase Doppler being close to zero. However, as the ray descends through the atmosphere, the refractive index becomes stronger, causing the ray to bend more, thus leading to a noticeable change in the excess Doppler delay [25].    With the aid of the quantities illustrated in Figure 6, the excess phase Doppler can be expressed as:

Inversion of Refractive Index Product and Temperature Product
Thus, the above equation has only two unknowns, δ A and δ B . Using Snell's law in a spherically symmetric medium, we have the following equation: n(r A ) and n(r B ) represent the refractive indices at the positions of the FY3E and GNSS satellites, respectively, and can be approximated as 1. Utilizing Equations (5) and (6), the bending angle α i can be calculated as follows: Eliminating the ionospheric influence is crucial in inverting atmospheric products. Due to the dispersive effects of the ionosphere, the bending angle profiles calculated using Galileo E1 and reconstructed E5a exhibit differences. It can be assumed that the first-order term of the ionospheric bending angle can be eliminated through a linear combination using Equation (8) and a represents the impact parameter: In the upper atmosphere above 40 km, the bending angle is relatively small, and the received signal is lower than the receiver's thermal noise level. Therefore, it is necessary to employ a statistical optimization method to further eliminate ionospheric errors [26]. Figure 7 shows the bending angle profiles originally and optimized by statistical optimization method, which can effectively correct the scene of sharp changes in the bending angle of the original data. This is because that the statistical optimization method combines measurement and bending angle of models. The coefficients are statistically optimized based on the noise characteristics of the bending angles, extrapolating from reliable lowaltitude data to unreliable high-altitude data. Additionally, it incorporates bending angle intervention obtained from the a priori climatological model in mass spectrometer and incoherent scatter (MSIS) data [27,28]. low-altitude data to unreliable high-altitude data. Additionally angle intervention obtained from the a priori climatological mo and incoherent scatter (MSIS) data [27,28]. Under atmospheric symmetry, the atmospheric refractive index can be obtained through the Abel integral transform using Equation (9): In the equation, the additional parameters represent the following: N represents the atmospheric refractive index, a 0 corresponds to the impact parameter specific to the current occultation. Both the bending angle and refractive index can be directly assimilated into numerical weather prediction models. With the assistance of numerical weather prediction models as prior information on the atmospheric state, the refractive index can be further inverted into meteorological parameters such as temperature using Equation (10): N = 77.6 P T + 3.73 · 10 5 P W T 2 (10) In the equation, P represents atmospheric pressure in units of hPa, T represents the absolute temperature of the atmosphere in Kelvin (K), and P W represents water vapor pressure. The first term on the right side represents the contribution of dry air to the refractive index, while the second term represents the contribution of water vapor. By neglecting the influence of moisture in the atmosphere, we can obtain the dry temperature: When inverting for single-frequency occultation products, the largest error comes from the delay caused by the radio signal passing through the ionosphere. In this paper, a reconstruction of the second frequency is performed to calculate the excess phase, which is then used to validate the data through excess phase Doppler. In the bending angle calculation, the error associated with the first-order term of the ionosphere is eliminated, and the upper atmospheric ionospheric error is mitigated with the assistance of climatological models. This approach enables the inversion of occultation products with higher precision.

Distribution and Number of Galileo E1 Single-Frequency Occultation Events
As mentioned above, if the occultation events are more evenly distributed globally, the consistency of the precision of global numerical weather prediction will be improved in terms of geographical distribution. Figure 8A-C illustrate the global distribution of atmospheric occultation events received by FY3E/GNOS for Galileo, GPS, and BDS on 1 January 2022. The red upright triangles represent rising occultation events, while the blue inverted triangles represent setting occultation events. From the figures, the number of rising and setting occultation events is comparable for a day. Similar to GPS and BDS, Galileo occultation events are uniformly distributed worldwide.
When inverting for single-frequency occultation products, the largest error comes from the delay caused by the radio signal passing through the ionosphere. In this paper, a reconstruction of the second frequency is performed to calculate the excess phase, which is then used to validate the data through excess phase Doppler. In the bending angle calculation, the error associated with the first-order term of the ionosphere is eliminated, and the upper atmospheric ionospheric error is mitigated with the assistance of climatological models. This approach enables the inversion of occultation products with higher precision.

Distribution and Number of Galileo E1 Single-Frequency Occultation Events
As mentioned above, if the occultation events are more evenly distributed globally, the consistency of the precision of global numerical weather prediction will be improved in terms of geographical distribution. Figure 8A-C illustrate the global distribution of atmospheric occultation events received by FY3E/GNOS for Galileo, GPS, and BDS on 1 January 2022. The red upright triangles represent rising occultation events, while the blue inverted triangles represent setting occultation events. From the figures, the number of rising and setting occultation events is comparable for a day. Similar to GPS and BDS, Galileo occultation events are uniformly distributed worldwide. If the Galileo single-frequency occultation process is successful, it will significantly increase the number of FY3E occultation inversion products. Figure 9 shows statistics of FY3E/GNOS occultation events received from 1 January to 14 January 2022. On average, FY3E/GNOS received 480 Galileo occultation events per day, 790 BDS occultation events, and 630 GPS occultation events. Conducting occultation events with all three systems increased the number by 33.8% compared to observations with BDS and GPS. As shown in the figure, the continuous 14-day occultation reception quantity remains relatively stable, indicating that the instrument could consistently receive occultation events from various navigation satellite systems. If the Galileo single-frequency occultation process is successful, it will significantly increase the number of FY3E occultation inversion products. Figure 9 shows statistics of FY3E/GNOS occultation events received from 1 January to 14 January 2022. On average, FY3E/GNOS received 480 Galileo occultation events per day, 790 BDS occultation events, and 630 GPS occultation events. Conducting occultation events with all three systems increased the number by 33.8% compared to observations with BDS and GPS. As shown in the figure, the continuous 14-day occultation reception quantity remains relatively stable, indicating that the instrument could consistently receive occultation events from various navigation satellite systems.

Precision Analysis of the Refractive Index Product
In the preliminary validation process, we inverted the refractive index product based on GPS L1 and BDS B1 single-frequency RO data, and the precision compared to the ERA5

Precision Analysis of the Refractive Index Product
In the preliminary validation process, we inverted the refractive index product based on GPS L1 and BDS B1 single-frequency RO data, and the precision compared to the ERA5 reanalysis is satisfactory. Furthermore, the precision of the single/dual-frequency occultation process for refractive index products is within the same order of magnitude, indicating the correctness of the single-frequency occultation process in Section 2. In this section, we process the actual FY3E/GNOS Galileo E1 single-frequency occultation observation data.
The reference standard for the refractive index product of a single-frequency occultation process is ERA5 reanalysis. ERA5 reanalysis is a reliable compilation of extensive observed data from satellites, weather stations, buoys, and ships, as well as predictive data from climate models. It provides comprehensive meteorological data covering global and all weather conditions, including temperature, humidity, precipitation, wind speed, air pressure, cloud cover, and surface radiation. These data can be downloaded for free from the website [29]. ERA5 reanalysis does not directly provide the refractive index. Instead, Equation (12) is used to convert specific humidity and air pressure from the analysis data into water vapor pressure. Then, employing the aforementioned Equation (10), the water vapor pressure and atmospheric temperature are converted into the refractive index.
where q represents the specific wetness, g/kg; ε = R d R v , R d and R v represent the dry air and wet air gas constants, respectively, ε ≈ 0.622. The remaining parameters have the same meaning as the parameters in Section 2.
The refractive index profiles obtained from Galileo single-frequency occultation data are compared with those calculated from ERA5 reanalysis. The average deviation index and RMS error index are calculated to analyze the level of deviation between the two. The experiment statistically analyzed the FY3E actual measured Galileo E1 single-frequency RO data on 3 January and 10 January 2022. The precision of the inverted refractive index is shown in Figure 10. The solid line reflects the average deviation percentage between processed by Galileo single-frequency occultation and calculated using ERA5 reanalysis. The dashed line represents the distribution of the RMS deviation between the two, with an auxiliary line indicating the 2% level. The two subfigures show similar trends in their variations.
From the solid line in the figure, the average deviation of the refractive index is less than 1% and is distributed near zero between heights of 5 km and 40 km. The consistency between the Galileo E1 single-frequency RO and the ERA5 reanalysis is very high, with the RMS error below 2%.
At lower altitudes near the Earth's surface, the average deviation index of the refractive index is relatively larger. This deviation may be due to imperfect attenuation of noise during the measurement caused by reconstructing the second frequency, even with the regularization method of minimizing the second-order derivative [30]. Furthermore, inverting the refractive index assumes atmospheric symmetry and utilizes the Abel integral transform of the bending angle. The Abel transform suppresses high-frequency noise to some extent. However, at lower altitudes of the atmosphere, there can be interference from ground reflections and the influence of turbulence, laminar flow, and inhomogeneities in the atmospheric layer. This leads to significant errors in the refractive index calculated using the Abel transform. When calculating the bending angle, errors are integrated and compounded, resulting in the final deviation.
processed by Galileo single-frequency occultation and calculated usi The dashed line represents the distribution of the RMS deviation be an auxiliary line indicating the 2% level. The two subfigures show si variations. From the solid line in the figure, the average deviation of the re than 1% and is distributed near zero between heights of 5 km and 40 between the Galileo E1 single-frequency RO and the ERA5 reanalys the RMS error below 2%.
At lower altitudes near the Earth's surface, the average deviation tive index is relatively larger. This deviation may be due to imperfect during the measurement caused by reconstructing the second frequ regularization method of minimizing the second-order derivative [3 verting the refractive index assumes atmospheric symmetry and utili transform of the bending angle. The Abel transform suppresses high some extent. However, at lower altitudes of the atmosphere, there can ground reflections and the influence of turbulence, laminar flow, and the atmospheric layer. This leads to significant errors in the refract At altitudes above 30 km, the overall RMS is less than 4%, but the deviation of the refractive index increases. The sources of error could be attributed to the following three reasons:

1.
When reconstructing the second frequency carrier phase, pseudorange measurement noise is introduced. Additionally, even after smoothing the relTEC, there can still be residual noise from the ionosphere.

2.
Errors in the prior climatological model near the transition layer have a certain impact on the precision of refractive index inversion in the upper atmosphere [31].

3.
Because the actual background refractive index is smaller above 30 km, this leads to larger relative errors in the calculations.
We also statistics the lowest penetration height rate at which the occultation profile penetrated vertically through the atmosphere from 10 to 12 January 2022 in the Figure 11. Among the 984 Galileo RO events that can be inverted by refractive index products, nearly 49% of the lowest penetration heights are between 5 and 8 km deep, and the frequency of RO events below 5 km and above 8 km is 29% and 22%, respectively. Due to resource limitations, Galileo RO does not employ an open loop tracking mode in the lower atmospheric layers. Instead, only occultation events with closed loop tracking in the lower atmosphere are used. This leads to a number reduction in the effective occultation observations below 5 km.
limitations, Galileo RO does not employ an open loop tracking mode in the pheric layers. Instead, only occultation events with closed loop tracking in mosphere are used. This leads to a number reduction in the effective occulta tions below 5 km.

Precision Analysis of Dry Temperature Product
To further understand the impact of different geographical locations o single-frequency occultation temperature product, Figure 12 presents the dry precision profiles for various latitudes from 1 to 14 January 2022. The solid lin the average deviation between the temperature inverted from Galileo sing RO and ERA5 reanalysis. The dashed line represents the RMS error of the dev plots (A-C) correspond to statistical comparisons of low latitude (up to 30°N itude (between 30° and 60°), and high latitude (above 60°) regions, respectiv Regardless of the low-, mid-, or high-latitude Galileo single-frequenc temperature product, the precision characteristics are quite similar. In term deviation, the temperature products show almost no bias compared to the E sis. Analyzing the RMS, the average temperature differences within the 5-35 Figure 11. The lowest penetration height rate on 10-12 January 2022.

Precision Analysis of Dry Temperature Product
To further understand the impact of different geographical locations on the Galileo single-frequency occultation temperature product, Figure 12 presents the dry temperature precision profiles for various latitudes from 1 to 14 January 2022. The solid line represents the average deviation between the temperature inverted from Galileo single-frequency RO and ERA5 reanalysis. The dashed line represents the RMS error of the deviations. Subplots range have an RMS below 3 K. Above 30 km altitude, the RMS rapidly increases and remains at approximately 3 K beyond 35 km. The temperature precision profiles exhibit good consistency and reliability within the 5-30 km range. From Figure 12A, at altitudes of 5 km and below, there are fewer effective data available. This is because, compared to high latitudes, lower latitudes have higher atmospheric humidity, and the higher humidity environment poses significant challenges for signal closure tracking [32,33]. It is worth noting due to resource constraints, GNOS receives the single-frequency Galileo RO signals in closed-loop mode, the tracking height is higher than that of the middle and high latitudes.

Summary
The Galileo occultation signal reception of the GNOS payload designed for the Feng-yun3 series E/F/G/H satellites is based on the E1 single frequency. In this study, the actual Regardless of the low-, mid-, or high-latitude Galileo single-frequency occultation temperature product, the precision characteristics are quite similar. In terms of average deviation, the temperature products show almost no bias compared to the ERA5 reanalysis. Analyzing the RMS, the average temperature differences within the 5-35 km altitude range have an RMS below 3 K. Above 30 km altitude, the RMS rapidly increases and remains at approximately 3 K beyond 35 km. The temperature precision profiles exhibit good consistency and reliability within the 5-30 km range.
From Figure 12A, at altitudes of 5 km and below, there are fewer effective data available. This is because, compared to high latitudes, lower latitudes have higher atmospheric humidity, and the higher humidity environment poses significant challenges for signal closure tracking [32,33]. It is worth noting due to resource constraints, GNOS receives the single-frequency Galileo RO signals in closed-loop mode, the tracking height is higher than that of the middle and high latitudes.

Summary
The Galileo occultation signal reception of the GNOS payload designed for the Fengyun3 series E/F/G/H satellites is based on the E1 single frequency. In this study, the actual Galileo single-frequency RO data obtained from the FY3E/GNOS receiver are used for the single-frequency occultation process. In terms of the number of occultation observations, the Galileo atmospheric occultation events received by FY3E/GNOS are distributed globally, similar to BDS and GPS occultation observations. The distribution is uniform, and the quantity is stable, which effectively improves the output of FY3E/GNOS occultation products and provides more data sources for global numerical weather prediction.
This study evaluates the precision of the products inverted from Galileo singlefrequency RO data. The refractive index product and dry temperature product show good precision in the altitude range of 5-30 km. However, at lower altitudes, there are influences from ground reflections and the presence of turbulence, laminar flow, and inhomogeneities, which result in slightly larger biases in the refractive index and dry temperature average deviation. In the upper atmosphere, the inversion products are affected by high-frequency ionospheric and residual pseudorange noise, leading to increased errors.
In summary, the actual Galileo E1 single-frequency occultation data received by GNOS are accurate. The inversion products of single-frequency Galileo occultation serve as a valuable supplement to the FY3 occultation products, thereby enhancing the quantity of global occultation products.

Discussion
The GNSS RO remote sensing payload GNOS carried on FY3E for measuring atmospheric occultation data has the unique feature of receiving occultation signals from the GPS, BDS, and Galileo constellations simultaneously. However, the Galileo occultation signals received by FY3E/GNOS are single frequency. This study demonstrates the precise refractive index product and dry temperature product, indicating the correctness of Galileo single-frequency GNSS RO data. Therefore, it can contribute to increasing the number of global occultation events for numerical weather prediction.
The precision of single-frequency occultation process products above 30 km is significantly affected by factors such as high-frequency ionospheric and pseudorange noise residuals. However, the single-frequency occultation process products in the range of 5-30 km possess a considerable level of reliability. In terms of numerical assimilation and related applications, it is recommended to decrease the weighting of single-frequency process products below 5 km and above 30 km. The single-frequency occultation process employed in this study is derived from the Ørsted satellite and has been successfully validated in the FY3E GPS and BDS single-frequency occultation processes. It has also been applied to Galileo's single-frequency occultation process. It is believed that this method has the potential to be applied to single-frequency occultation processes in other satellite occultation missions to increase the number of available occultation events.
Not only for the single-frequency occultation process but also for the L2P band of the GPS dual-frequency occultation process, which exhibits poor tracking near the bottom of the atmosphere, the method of second frequency reconstruction introduced in this paper can be considered to eliminate ionospheric errors. This means that a single-frequency occultation process can to some extent enhance the number of occultation data sources for global numerical weather prediction, thereby improving its precision. If there is more support from GNSS RO inversion data for research on numerical weather prediction and global climate change, it will greatly enhance the application capabilities of meteorological services, regardless of whether these GNSS RO inversion products come from dual-frequency or single-frequency occultation data.