Introduction of TROPS ionospheric TEC products for FORMOSAT-7/COSMIC-2 mission

Six low Earth orbit (LEO) satellites were launched on June 25th, 2019 for a radio occultation (RO) mission for the FORMOSAT-7/COSMIC-2 (F7/C2) program. The GPS and GLONASS RO signals received by these F7/C2 satellites can be used to retrieve atmospheric and ionospheric parameter profiles for atmospheric and ionospheric research. In order to process the received RO signal, the processing system named Taiwan Radio Occultation Processing System (TROPS) is built. TROPS is developed by National Space Organization, Taiwan Analysis Center for COSMIC, and GPS Science and Application Research Center in Taiwan. The ionospheric products of TROPS are electron density profile, ionospheric scintillation index (S4 index), and absolute total electron content (TEC). S4 index has been calculated on board the satellites and other two products are retrieved by TROPS after the observation data downlink to ground. TEC is the linear integration of electron density along the signal propagation path. The electron density profile is retrieved from the relative TEC when the elevation angle of GNSS satellite is negative from F7/C2 satellite. The absolute TEC is the TEC from GNSS satellite to F7/C2 satellite. The difference between absolute and relative TEC is the TEC with/without differential code bias (DCB) correction. Currently, the data for the electron density profile and absolute TEC are provided by TROPS. Users can obtain the products freely from the internet. In this study, the retrieval method and the preliminary F7/C2 ionospheric TEC products retrieved by TROPS are presented in detail. Electron density profiles data service of FORMOSAT-7/COSMIC-2. Absolute TEC data service of FORMOSAT-7/COSMIC-2. Validation results confirms the data is ready to use. Electron density profiles data service of FORMOSAT-7/COSMIC-2. Absolute TEC data service of FORMOSAT-7/COSMIC-2. Validation results confirms the data is ready to use.


Introduction
The Taiwan Radio Occultation Processing System (TROPS) was developed for data processing and retrieval for the FORMOsa SATellite Mission-7/Constellation Observing System for Meteorology, Ionosphere, and Climate-2 (FORMOSAT-7/COSMIC-2, F7/C2) program, which is an international collaboration between the National Space Organization (NSPO) of Taiwan and the National Oceanic and Atmospheric Administration

Open Access
Terrestrial, Atmospheric and Oceanic Sciences  (Fong et al. 2008), which is also an international collaborative effort between NSPO and the University Corporation for Atmospheric Research (UCAR) of the United States. F3/C was the first weather satellite mission in Taiwan with six low Earth orbit (LEO) satellites and was launched in 2006. The F3/C products allowed the generation of atmospheric and ionospheric parameter profiles by using radio occultation (RO) technique (Schreiner et al. 1999;Wickert et al. 2001).
The mission payload of F3/C, the GPS Occultation Experiment (GOX), was used for RO observation and to receive Global Positioning System (GPS) signals for retrieval of the atmospheric and ionospheric parameter profiles and can provide around 2000 pressure, temperature, and electron density profiles per day. Due to the success of F3/C, the F7/C2 program was later developed. The mission payload of F7/C2, Tri-GNSS Radio Occultation System (TGRS) can not only receive GPS signals but also GLObal NAvigation Satellite System (GLONASS) for RO observation. The F7/C2 data products are similar to F3/C's but the generated data volume is about twice that of F3/C. The ionospheric products of the F3/C GOX payload have been used for the study of ionospheric activity Liu et al. 2010;Chang et al. 2014), space weather nowcasting (Lin et al. 2015(Lin et al. , 2017 and forecasting (Lee et al. 2012;Hsu et al. 2014), as well as detection of ionospheric irregularities (Yeh et al. 2012(Yeh et al. , 2014Liu et al. 2016;Chen et al. 2017), and plasmaspheric activity (Chen 2008;Cherniak et al. 2012). Due to the difference in the orbital inclination angle, from 72° for the F3/C satellites to 24° for the F7/C2 satellite, the data distribution of the F7/C2 TGRS is roughly between 45° N to 45° S. Although the data distribution of F7/C2 only includes the mid-and lower-latitudes, the density is around 4 times that of F3/C, because of having double the data volume and half the distribution area. The ionospheric products of F7/C2 allow more detailed ionospheric research studies.
In order to thoroughly understand the process of the data retrieval and achieve the data latency required by the F7/C2 mission, it is necessary to obtain as much information as possible for data applications. TROPS was developed as a joint effort between NSPO, Taiwan Analysis Center for COSMIC (TACC) of Central Weather Bureau (CWB), and GPS Science and Application Research Center (GPSARC) at National Central University (NCU) in Taiwan. The architecture of TROPS is shown in Fig. 1. The architecture is divided into three main segments for data collection, data processing, and data service, shown on the left, top right, and bottom right sides of Fig. 1, respectively. In data processing portion, the blue ovals indicate the modules for precision orbit determination (POD) and green ones are for atmospheric and ionospheric data retrieving. The POD modules are not only for Global Navigation Satellite System (GNSS) and LEO satellites' POD but also for the calculation of intermediate data, such as the satellite clock offset (Tseng et al. 2018). The main products of the atmospheric data retrieving module are the pressure, temperature, and water vapor profiles. The main products of the ionospheric data retrieving module are the electron density profile, ionospheric scintillation S4 index (S4 index), and absolute total electron content (TEC). Since the S4 index is calculated on board the satellite and does not need to be retrieved, the S4 index of TROPS has been compared with the results of COSMIC Data Analysis and Archive Center (CDAAC) retrieval system operated in TACC to confirm the offset and scaling factors for file transferring by using ~ O(10) observations. Only the details of electron density profile, which is retrieved from relative TEC, and absolute TEC retrieval process module are presented in this study. The difference between relative and absolute TEC is the TEC without/with difference code bias (DCB) correction.
The retrieval method used in the module is described in Sect. 2. The retrieval results are displayed, followed by the conclusions.

Highlights of FORMOSAT-3/COSMIC retrieval methods
The ionospheric RO observations are illustrated in Fig. 2. The signal is transmitted from the left side GNSS satellites (G 1 -G 4 ) and propagates through Earth's ionosphere, then it is received by the TGRS of F7/C2 satellite (point F 1 -F 4 ). The influence of the received GNSS signal by the ionosphere is presented in the phase and amplitude of the received signal. The dual band of the GNSS signal and the different frequency causes a difference in the refractive index when the signal propagates through a medium with the same electron density. The total electron content (TEC) can be calculated from observed pseudorange (TEC p ) and carrier phase (TEC c ) by where λ 1 , λ 2 , N 1 , N 2 , k t , k r , f 1 , f 2 , P 1 , P 2 , S 1 , and S 2 , are wavelength of L1 band, wavelength of L2 band, (1a) Page 3 of 9 Yeh et al. Terrestrial, Atmospheric and Oceanic Sciences (2022) 33:23 ambiguity of L1 band, ambiguity of L2 band, differential code bias (DCB) of transmitter, DCB of receiver, the frequency of the L1 band, the frequency of the L2 band, the pseudorange of L1 band, the pseudorange of L2 band, the carrier L1 phase, and the carrier L2 phase, respectively (Liu et al. 1996). The TEC from Fig. 1 Architecture of TROPS. The left side, upper right side, and lower left side show the sections for data collection, data processing, and data service, respectively. In the data processing portion, the blue oval indicates the POD module and the green ovals are the modules for atmospheric and ionospheric data retrieval

Fig. 2 Illustration of RO ionospheric observations
Page 4 of 9 Yeh et al. Terrestrial, Atmospheric and Oceanic Sciences (2022) 33:23 point G 1 to point F 1 (TEC G1F1 ) and point G 1 to point F 1 ʹ (TEC G1F1ʹ ) can be calculated by Eq. (1). The TEC from point F 1 ʹ to point F 1 (TEC F1ʹF1 ) can be calculated by (TEC F1ʹF1 ) = (TEC G1F1 ) − (TEC G1F1ʹ ). Note that TEC G1F1 and TEC G1F1ʹ are relative to TEC, which does not remove the influence of the signal ambiguity and DCB. Since the ambiguity and DCB in one arc observation (from point F 1 ʹ to point F 1 ) is the same, (TEC G1F1 ) − (TEC G1F1ʹ ) can remove the influence of the signal ambiguity and DCB. In Fig. 2, points F 1 ʹ to F 4 ʹ are the reference points for F 1 to F 4 , respectively. Furthermore, F 1 to F 4 are in the same LEO satellite orbit and TEC F1ʹF1 , TEC F2ʹF2 , TEC F3ʹF3 , and TEC F4ʹF4 can be derived from one arc observation. Although only 4 TECs between the position of the LEO satellite and its reference point are derived with the signal propagating through different altitudes in Fig. 2, several hundred TECs can be derived from real RO observations. Under the assumption that the ionosphere is spherically symmetrical, the Abel transform (Phinney and Anderson 1968), which is often used to analyze spherically symmetric functions, is applied to transfer the TEC profile to the electron density profile. The Abel transform equation for electron density derived from the TEC is (Schreiner et al. 1999) where N, T, p, r, and r LEO are the electron density, TEC (the TEC between the LEO position and its reference point), altitude, the altitude of the retrieval electron density, and the altitude of the LEO satellite, respectively. The GNSS satellites (G1 to G4) on the left side in Fig. 2 are on the RO side. The elevation angle of the GNSS satellites from F7/C2 is negative, and their signals can be used to retrieve the electron density profiles. The F7/C2 TGRS not only receives the GNSS signals from the RO side but also from POD side, where the elevation angle of the GNSS satellites from F7/C2 is positive; Satellites G5 to G7 are shown on the right hand side in Fig. 2. The signal received on the POD side of the GNSS satellites is used for POD and also can be used for the calculation of the TEC along the signal propagation path. The TEC can be calculated by using Eq. (1). However, it is necessary to remove the influence of ambiguity, which is the integer refers to the first epoch of observations and (2) N(r)=-1 π r LEO r dT /dp p 2 − r 2 dp remains constant during the period of observation and indicated by N 1 and N 2 in Eq. (1b), and the DCB, which is the systematic bias between two GNSS code observations and indicated by k t and k r in Eq. (1a), before calculating the TEC. The DCB is calculated from the POD side of the GNSS signal from the last 3 days. Due to the signal of GPS being based on code division multiple access (CDMA) and GLONASS being based on frequency division multiple access ( where θ and mapf are the elevation angle of the GPS satellite from F7/C2 and the mapping function, respectively. The mapping function is from Foelsche and Kirchengast (2002) where r ion and r orb are the ionospheric and satellite altitude from the Earth's center. r ion can be selected to be several hundreds or thousands of kilometres above r orb (Yue et al. 2011) and r orb + 200 km is used in TROPS. For GLONASS signal, Eqs. (3) and (4) are also used for calculation of DCB, but TEC 1 and TEC 2 are calculated from the same GLONASS satellite at different times with mapf(θ 1 ) and mapf(θ 2 ) as the corresponding mapping function. The calculation results of GLONASS DCB is the summation of k t and k r . On the other hand, the ambiguity of L1 and L2 band in Eq. (1b) is the least square results of N i = (P i − λ i S i − 2d ion )/λ i , i = 1, 2 implies L1 and L2 band, in one arc observation and d ion = 40.3TEC c /f i 2 is ionospheric delay. After calculating ambiguity and DCB, the absolute TEC can be derived by Eq. (1). ( Page 5 of 9 Yeh et al. Terrestrial, Atmospheric and Oceanic Sciences (2022)

FORMOSAT-7/COSMIC-2 mission status retrieval results and discussion
Six examples of electron density profiles in the day of year 2020.039 are shown in Fig. 3. The altitude range of the retrieval electron density profile is from several tens of kilometers to the satellite altitude. The profiles in panel (a) to (f ) are from satellites FS701 to FS706. At day of year 2020.039, FS701 and FS704 are at the mission altitude and the others are at the original altitudes, which are 550 km and 720 km, respectively. So, the highest altitude of the profiles in panels (a) and (d) is around 550 km and the others are around 720 km. The profiles in Fig. 3 show the obvious F region which F2-layer peak density (NmF2) and height (HmF2) are 10 5 to 7 × 10 5 /cm 3 in density and around 200 to 400 km in altitude. In panels (b), (c), (d), and (e) of Fig. 3, the sporadic E (Es) layer occur in the profiles around 100 km in altitude. The electron density of Es layer in panels (c) and (d) is up to around half and full NmF2, respectively. However, due to the assumption of Abel transform, the spherically symmetric electron  Page 6 of 9 Yeh et al. Terrestrial, Atmospheric and Oceanic Sciences (2022) 33:23 density in ionosphere, the negative value occurs in profiles in low altitude, like the profile in panel (a). The variation in the daily amount of electron density profile in 2019 and 2020 is shown in Fig. 4. As can be seen in Fig. 4, F7/C2 can provide around four thousand and sometimes up to five thousand electron density profiles per day if all six satellites are working. A comparison between the F7/C2 electron density profiles and ionosonde observations from day of year 001 to 100 in 2020 is shown in Fig. 5. The NmF2 and HmF2 are used for the comparisons. The ionosonde observations are from the Digital Ionogram Data Base (DIDBase). The time and spatial difference between the F7/C2 and ionosonde observations are 5 min and 3°, respectively, and a total of around 3250 pairs of F7/C2 and ionosonde observations are compared. The blue and green lines in Fig. 5 indicate the lines where the slope is equal to 1 and the linear regression line, respectively. The green line in Fig. 5a is very close to the blue line and the NmF2 of F7/ C2 and ionosonde is in good agreement. The green line in Fig. 5b is not as close as to the blue line in Fig. 5a and HmF2 of F7/C2 and ionosonde is not in good agreement as NmF2. The disagreement of HmF2 between F7/C2 and ionosonde is due to ionosonde cannot directly scale the true HmF2 from ionograms and only provide the measurement of the virtual height. The HmF2 provided by ionosonde is derived by the variation of electron density with height is parabolic with plasma frequency under the assumption of a simple ionospheric layer, or by the conversion of a plasma frequency versus virtual height curve to an electron density true height profile (Tsai et al. 2009). Although the agreement of HmF2 is not as good as NmF2 between F7/C2 and ionosonde, the positive slope of regression line and the linear distribution of the red dots can validate the TROPS retrieval electron density profile.
Two examples of the GNSS signal DCB are shown in Figs. 6 and 7. The green, blue, purple, and red curves in Fig. 6 show the variations of DCB for the L2C signal of antenna 1, L2P signal of antenna 1, L2C signal of antenna 2, and L2P signal of antenna 2, respectively, of the F7/C2 FS701 satellite, from day of year 001 to 365 in 2020. As can be seen in Fig. 6, the variation of the GPS L2C and L2P signal DCB of F7/C2 FS701 antennas 1 and 2 is stable, with a maximum of around 2 NS perturbations and the trends are around − 1 NS in 1 year. The colored curves in Fig. 7a and b indicate the DCB variations of GLONASS R01 to R26 satellite signals of antenna 1 and 2, respectively, for the F7/C2 FS706 satellite, from day of year 001 to 080 in 2020. As can be seen in Fig. 7, the variation of the GLONASS signal DCB, like of the GPS, is stable with perturbations. The colored curves in Fig. 7b  show the maximum variation of the DCB of the GLO-NASS satellite signal for antenna 2 is around 2 NS. However, as shown in Fig. 7a, the perturbation of the variation in the DCB of the GLONASS satellite signal is larger for antenna 1 than antenna 2, up to a maximum of around 4 NS. The higher perturbation shown by the colored curves in Fig. 7a than in Fig. 7b is due to fewer observations for analysis for F7/C2 FS706 antenna 1 than antenna 2. In early 2020, before the TGRS software was updated, the ability of TGRS to receive GLONASS signals was not as good as it is today, so most of the observations for antenna 1 were for rising RO, which is more difficult to receive GLONASS signal than setting RO for antenna 2.
The gaps in Figs. 6 and 7 are due to no observations for analysis when the TGRS power was off.
Six examples of retrieved absolute TEC profiles are shown in Fig. 8. The black and red curves are derived by using Eqs. (1a) and (1b), respectively. Due to the low accuracy of pseudorange, the black curves have larger perturbations than red curves. The profiles in panels  and GLONASS signals at the same time and the angle between the GPS to F7/C2 and GLONASS to F7/C2 is smaller than 3° it is retrieved as the absolute TEC product, and the GPS and GLONASS absolute TEC pair is used for comparison. From day of year 001 to day of year 100 in 2020, around 36 such pairs are found. The numerical distribution of the difference between GPS and GLONASS absolute TEC is shown in Fig. 9. The mean difference and the standard deviation of the TEC difference between the GPS and GLONASS signals are around 0.7 TEC units (TECu) and 2.7 TECu. The small TEC difference between GPS and GLONASS can validate the absolute TEC of the TROPS retrieval.

Conclusions
TROPS is a Taiwan built data processing system for RO observation. After the launch of the six LEO satellites on June 25th, 2019, named F7/C2, TROPS began to provide F7/C2 ionospheric retrieval products in March, 2020. The F7/C2 mission provided three products for the ionosphere. Except for the S4 index, which was calculated on board the satellite, the details of the retrieval process for the electron density profiles and absolute TEC are described in this study. The retrieval method used in TROPS is almost similar to the method of CDAAC in TACC. However, the number of electron density profile of TROPS is usually more than CDAAC in TACC for around several hundred after 2020. The reason perhaps due to the electron density profile need at least around continuous 10 min observations for retrieving and sometimes one profile is separated into two continuous dumps and TROPS connects the two dump separated profile well. Users can freely obtain the F7/C2 ionospheric products retrieved by TROPS from the TACC (https:// tacc. cwb. gov. tw/ v2/ en/ trops_ downl oad. html) and Space Weather Operational Office (https:// swoo. cwb. gov. tw/) website. Fig. 9 Number distribution of TEC difference between GPS and GLONASS satellites received by TGRS. The angle between GPS to F7/ C2 and GLONASS to F7/C2 is less than 3° at the same time from day of year 001 to 100 in 2020. A total of 36 pairs of absolute TECs are found and the mean and standard deviations of the difference are around 0.7 TECu and 2.7 TECu, respectively