Soil MoistureRetrieval from TM-1 GNSS-R Reflections with Auxiliary Geophysical Variables: A Multi-Cluster and Seasonal Evaluation

Current passive microwave satellites like SMAP still face limitations in observational frequency and responsiveness in regions with frequent cloud cover, dense vegetation, or complex terrain, making it difficult to achieve continuous global monitoring with high spatio-temporal resolution. To enhance global high-frequency monitoring capabilities, this study utilizes global reflectivity data provided by the Tianmu-1 (TM-1) constellation since 2023, combined with multiple auxiliary variables, including NDVI, VWC, precipitation, and elevation, to develop a 9 km resolution soil moisture retrieval model. Several spatial clustering and temporal partitioning strategies are incorporated for systematic evaluation. Additionally, since the publicly available TM-1 L1 reflectivity data does not provide separable polarization channels, this study uses DDM/specular point reflectivity as the primary observable quantity for modeling and mitigates non-soil factor interference by introducing multi-source priors such as NDVI, VWC, precipitation, terrain, and roughness. Unlike SMAP’s “single orbit daily fixed local time” observation mode, TM-1, leveraging multi-constellation and multi-orbit reflection geometry, offers more balanced temporal sampling and availability in cloudy, rainy, and mid-to-high latitude regions. This enables temporal gap filling and rapid event response (such as moisture transitions within hours after precipitation events) during periods of SMAP’s quality masking or intermittent data loss. Results indicate that the model combining LC-cluster with seasonal partitioning delivers the best performance at the cluster level, achieving a correlation coefficient (R) of 0.8155 and an unbiased RMSE (ubRMSE) of 0.0689 cm3/cm3, with a particularly strong performance in barren and shrub ecosystems. Comparisons with SMAP and ISMN datasets show that TM-1 is consistent with mainstream products in trend tracking and systematic error control, providing valuable support for global and high-latitude studies of dynamic hydrothermal processes due to its more balanced mid- and high-latitude orbital coverage.