Land Surface Reflection Differences Observed by Spaceborne Multi-Satellite GNSS-R Systems

With the accelerated launch of spaceborne Global Navigation Satellite System Reflectometry (GNSS-R) satellites, GNSS-R has gradually emerged as an important technique for remote sensing. However, due to its pseudo-random observation mode, the use of a single system makes it difficult to provide continuous spatiotemporal coverage over a specific area within the short term. Although interpolation methods can partially alleviate the coverage gaps, their application is limited by accuracy and reliability constraints, which still restrict the practical use of GNSS-R in terrestrial surface monitoring. To address this issue, conducting joint analyses and data fusion of multi-satellite GNSS-R observations has become an important approach to improving the continuity and accuracy of surface monitoring. However, systematic studies on the integration of multi-satellite GNSS-R data remain relatively limited. Moreover, differences in orbital inclination, antenna design, and signal bandwidth among various spaceborne GNSS-R systems lead to discrepancies in their land observations. Therefore, this study systematically analyzes the reflectivity differences among multiple GNSS-R satellites (e.g., the Cyclone Global Navigation Satellite System (CYGNSS), Fengyun-3 (FY-3), and Tianmu-1 (TM-1)) under consistent surface roughness and land cover conditions, with the aim of providing a theoretical and methodological foundation for the fusion and integrated application of multi-satellite GNSS-R data. The results show that, except for desert regions, the spatial distribution of the correlation coefficients from the least squares fitting of reflectivity between different spaceborne GNSS-R satellites exhibits a pattern similar to that of an established variable, i.e., the vegetation–roughness composite variable (VR), with higher inter-system correlations occurring in areas characterized by lower VR values. Significant reflectivity deviations were observed near water bodies and river networks, such as the Amazon, Paraná, Congo, Niger, Nile, Ganges, Mekong, and Yangtze, where both the fitting intercepts and biases are relatively large. In addition, the reflectivity correlations between CYGNSS–TM-1 and CYGNSS–FY-3 are both strongly influenced by surface vegetation cover type. As the correlation increases, the proportion of non-vegetated and forested areas decreases, while that of grasslands, shrublands, and cropland/vegetation mosaics increases. Analysis of inter-system reflectivity correlations across different land cover types indicates that forested areas exhibit low-to-moderate correlations but maintain stable structural characteristics, whereas wooded areas show moderate correlations slightly lower than those of forests. Grasslands, shrublands, and croplands are mainly distributed within regions of moderate surface roughness and correlation, among which croplands have the highest proportion of highly correlated grids, demonstrating the greatest potential for multi-source data fusion. Wetlands display high roughness and low correlation, largely influenced by dynamic water variations, while bare soils show low roughness (0.2–0.4) but still weak correlations.