Soil moisture (SM) is a critical physical parameter in land surface processes that affects the atmospheric and hydrological cycles. Soil Moisture Active Passive (SMAP) mission produces high-quality global SM data, but its low spatial resolution limits the applications at a regional or local scale. In this study, we developed a novel method to select optimal factors from 34 candidate downscaling factors for each month using three machine learning methods (Back Propagation Neural Network (BPNN), Support Vector Machine (SVM) and Random Forest (RF)) to produce monthly time series SM products at a spatial resolution of 1 km for the entire Loess Plateau in China from 2015 to 2023. 11 in-situ SM measurements distributed in Loess Plateau and precipitation data were used to evaluate the downscaling performances of the three machine learning approaches. The results show that the spatial trends of the three downscaled SM were essentially the same as the SMAP SM, with similar spatial patterns, and all three downscaled SM can provide more spatial information and texture features than SMAP SM. Among the three types of downscaled SM, the RF downscaled SM reached the highest R (0.654) and the smallest unbiased Root Mean Square Error (ubRMSE) (0.044 m3m−3), better than the BPNN and SVM downscaled SM, and most closely matched the in-situ measurements. The dynamics of SM were successfully captured by RF downscaled SM, which exhibited strong temporal consistency with the in-situ SM and responded well to precipitation events, with significant increases in SM values in the month of high precipitation and subsequent months. In conclusion, the RF model has the best downscaling effect and it can be used to provide high spatial resolution SM data for applications at a regional or local scale across the Loess Plateau.