Chinese BeiDou Navigation Satellite System (BDS) as a mixed constellation, consists of not only Medium Earth Orbit (MEO) satellites but also Inclined Geosynchronous Orbit (IGSO) and Geostationary Earth Orbit (GEO) satellites, of which the GEO satellites remain visible and stationary in the Asia- Pacific region almost all the time so as to better ensure the reflected signals' stability and continuity, compared to using only MEO satellites in previous GNSS-R remote sensing research. This paper demonstrates ocean altimetry capability using reflected signals from BDS GEO satellite with coastal experiments conducted in Weihai, Shandong, China on July 6 and July 9, 2019. By using only one GEO satellite, the precision of DDM (Delay Doppler Map) and phase altimetry is up to decimeter (close to millimeter) level. And the results show that BDS GEO- R technique can be used to monitor small changes in ocean altimetry with high temporal resolution.

The Yunyao Aerospace Constellation Program is the core project being developed by Yunyao Aerospace Technology Co., Ltd., Tianjin, China. It aims to provide scientific data for weather forecasting, as well as research on the ionosphere and neutral atmosphere. It is expected to launch 90 high time resolution weather satellites. Currently, the Yunyao space constellation provides nearly 16,000 BDS, GPS, GLONASS, and Galileo multi-system occultation profile products on a daily basis. This study initially calculates the precise orbits of Yunyao LEO satellites independently using each GNSS constellation, allowing the derivation of the neutral atmospheric refractive index profile. The precision of the orbit product was evaluated by comparing carrier-phase residuals (ranging from 1.48 cm to 1.68 cm) and overlapping orbits. Specifically, for GPS-based POD, the average 3D overlap accuracy was 4.93 cm, while for BDS-based POD, the average 3D overlap accuracy was 5.18 cm. Simultaneously, the global distribution, the local time distribution, and penetration depth of the constellation were statistically analyzed. BDS demonstrates superior performance with 21,093 daily occultation profiles, significantly exceeding GPS and GLONASS by 15.9% and 121%, respectively. Its detection capability is evidenced by 79.75% of profiles penetrating below a 2 km altitude, outperforming both GPS (78.79%) and GLONASS (71.75%) during the 7-day analysis period (DOY 169–175, 2023). The refractive index profile product was also compared with the ECWMF ERA5 product. At 35 km, the standard deviation of atmospheric refractivity for BDS remains below 1%, while for GPS and GLONASS it is found at around 1.5%. BDS also outperforms GPS and GLONASS in terms of the standard deviation in the atmospheric refractive index. These results indicate that Yunyao satellites can provide high-quality occultation product services, like for weather forecasting. With the successful establishment of the global BDS-3 network, the space signal accuracy has been significantly enhanced, with BDS-3 achieving a Signal-in-Space Ranging Error (SISRE) of 0.4 m, outperforming GPS (0.6 m) and GLONASS (1.7 m). This enables superior full-link occultation products for BDS.

Global Navigation Satellite System (GNSS) Reflectometry (GNSS-R) is the use and study of GNSS signals’ distortion and delay reflected or scattered by the earth’s surface to tell us the sea level, wind speed and direction, humidity, ice density, and even moving targets. The space-borne GNSS-Reflectometer from Tianjin Yunyao Aerospace Technology Co., Ltd. is China’s first commercial GNSS-R instrument. On July 3, 2021, this payload was equipped on the Jilin-1 wide-range 01B satellite (J1-01B) developed by Changguang Satellite Technology Co., Ltd. and entered the polar orbit with the height about 500km from the Taiyuan Satellite Launch Center.The attitude of J1-01B is stable to the ground within 0.1 degrees, which can guarantee the high accuracy of GNSS-R antenna gain model on the detection of sea surface wind speed and direction, soil moisture and etc. Meanwhile, its orbital characteristics ensure that its trajectory can cover the north and south poles of the earth, providing opportunities for researches on sea-ice detection and sea-ice thickness retrieval. The payload supports the mutual switching between dual system (GPS L1 or BDS B1) signal mode and dual working mode (generate the doppler delay maps (DDM) in manual-adjustable resolution on the satellite, or collect intermediate frequency (IF) data with 58Mb/s rate), which guarantees the verification and continuous development of multi-system compatibility, data quality analysis and application scenario expansion.After two weeks of satellite’s launch, the satellite has entered normal working mode. In two periods from 6:0 to 6:15 and from 7:32 to 7:42 on July 16, 2021, GNSS-R data were collected in DDMs as well as IFs formats, when J1-01B was over the Pacific Ocean. And L1a DDMs with the horseshoe shape were generated with correction from raw DDMs, as well as from IFs.The GNSS-R instrument on the Jilin-1 wide 01B satellite is an important component of the Yunyao constellation plan, which would provide rich and stable data, and be an important role in global and regional exploration in aerospace, meteorology, ocean, earthquake and other fields.

Through Global Navigation Satellite System (GNSS) occultation measurement, the global ionosphere and atmosphere can be observed. When the navigation satellites’ signal passes through the lower atmosphere, the rapid change of the atmospheric refractive index gradient will cause serious multipath phenomena in radio wave propagation. Atmospheric doppler frequency shift and amplitude signal fluctuations increase drastically. Due to the attenuation of signal amplitude and the rapid change of the Doppler frequency, the general phase locked loop (PLL) cannot work properly. Hence, a more stable tracking technology is needed to track the occultation signal passing through the lower atmosphere. In this paper, a mountain-top based radio occultation experiment is performed, where we employ an open-loop receiver and remove the navigation bits by the internal demodulation. In the process of the experiment, we adopt the open-loop tracking technique and there is no feedback between the observed signal and the control model. Specifically, taking the pseudo-range and doppler information from models as input, three key parameters, i.e., accurate code phase, carrier doppler and code doppler, can be obtained, and furthermore, the accurate accumulation is determined by them. For the full open-loop occultation data, a closed-loop observation assisted strategy is presented to compare the tracking results between open-loop and closed-loop occultation data. Through the compared results, we can determine whether the initial phase has been reversed or not, and obtain the high consistency corrected open-loop data that can be directly used for subsequent atmospheric parameters inversion. To verify the effect of open-loop tracking and open-loop inversion, we used the company’s self-developed occult receiver system for verification. The company’s self-developed occult receiver system supports Global Position System (GPS)/Beidou satellites constellation (BD, the 2nd and 3rd generations) dual systems. We have verified GPS and BD open-loop tracking and inversion, carried out in a three-week mountain-based experiment. We used closed-loop and open-loop strategies to track and capture the same navigation star to detect its acquisition effect. Finally, we counted the results of a week (we only listed the GPS data; BD’s effect is similar). The experimental results show that the open-loop has expanded the signal-cut-off angle by nearly 20% under the condition of counting all angles, while the open-loop has increased the signal-cut-off angle value by nearly 89% when only calculating the negative angle. Finally, the atmosphere profiles retrieved from observations in open-loop tracking mode are evaluated with the local observations of temperature, humidity and pressure provided by the Beijing Meteorological Bureau, and it is concluded that the error of open-loop tracking method is within ~4% in MSER (mean square error of relative error), which meets the accuracy of its applications (<5%, in MSER).

This article proposes a new sea surface wind speed (SSWS) retrieval modeling algorithm based on the empirical orthogonal function (EOF) analysis for observations acquired by the global navigation satellite system reflectometry (GNSS-R). As a nonparametric modeling algorithm, it is simpler compared with the nonlinear methods. The influence of wind speed and incident angle on the modeling error is analyzed for the first time using a spectrum analysis. Three types of data from 80% CYGNSS 2019–2020 observations [delay Doppler map average (DDMA) and leading edge slope (LES)], signal incident angle, and the European Centre for Medium-Range Weather Forecasts Reanalysis V5 (ERA5) reference wind speed are used in the EOF analysis to establish two retrieval models. The remaining 20% of the data are used for accuracy evaluation after getting the final wind speed by the minimum variance (MV) estimator. As a result, when using three 0–20-m/s wind speeds of ERA5, Advanced Scatterometer (ASCAT), and the Modern-Era Retrospective Analysis for Research and Applications V2 (MERRA2) as contrasts, the root mean squared errors (RMSEs) are 1.51, 1.45, and 1.43 m/s, respectively. Compared with CYGNSS wind product, the performance of this algorithm is closer to the L2 Climate Data Record (CDR) V1.1 product than V1.0. The results demonstrate that the EOF algorithm has a good performance in retrieving SSWS and can better retain the influence of the incident angle on the observations.

This study introduced a Kalman filtering assimilation model that considers the DCB errors of GPS/LEO satellites and GNSS stations. The assimilation results and reliability were verified by various types of data, such as ionMap, ionosonde, ISR, and the EDP of ionPrf from COSMIC. The following analyses were carried out. Assimilating the measured ground-based/spaceborne ionospheric observation data from DOY 010, 2008 and DOY 089, 2012 revealed that the introduction of GPS/LEO satellite and GPS station DCB errors can effectively suppress the STEC observation errors caused by the single-layer hypothesis. Furthermore, the top of the ionosphere contributes 2.8 TECU (approximately 10–20% of the STEC) of electrons during the ionospheric quiet period, greatly influencing the ionospheric assimilation at altitudes of 100–800 km. The assimilation results also show that, after subtracting the influence of the top of the ionosphere, the ionospheric deviation during the quiet period improved from 1.645 TECU to 1.464 TECU; when the ionosphere was active, the standard deviation was improved from 4.408 TECU to 3.536 TECU. The IRI-Imp model introduced by Wu et al. and the IRI (2007) model were used as background fields to compare the effects of COSMIC occultation observation data on the ionospheric assimilation process. Upon comparison, the occultation data introduced by the improved model showed the greatest improvement in the vertical structure of the ionosphere; additionally, the assimilation process reused the horizontal structure information of the occultation data, and the assimilation result (IRI-Imp-Assi) was the most ideal. Due to the lack of an occultation data correction, the IRI2007 model was relatively more prone to errors. With the strategy of the IRI-Imp-Assi model, the introduction of occultation data caused a more significant reduction in the error between the assimilation model with the IRI model as the background field and the ionMap.