Soil moisture (SM) is a critical component of the land surface hydrological cycle, significantly impacting various sectors such as hydrology, meteorology, and agriculture. Accurate, high-resolution SM data are essential for effective flood forecasting, water resource management, and understanding the soil freeze-thaw processes in cold regions. This study aims to generate 1 km resolution liquid surface SM (SSM) data with a twice-daily update frequency by downscaling SMAP Level-4 SSM data using random forest (RF) and multiple linear regression (MLR) in the source region of the Yellow River (SRYR), by considering the differences in SM changes between freezing and thawing periods. To obtain the SSM data, 16 downscaling schemes of both RF and MLR were designed for each of the three scenarios. In each downscaling process, both land surface temperature (LST) and normalized difference vegetation index (NDVI) were utilized in MLR and RF models, alongside various combinations of additional variables such as albedo, elevation, leaf area index (LAI), soil texture. Results showed that during the freezing period, RF produced superior SSM estimates when supplemented with NDVI, LST, albedo, elevation, LAI, and soil texture. MLR was more effective during the thawing period when paired with NDVI, LST, elevation, LAI, and soil texture. During the freezing period, the downscaled SMAP SSM exhibited average R, RMSE, ubRMSE of 0.76, 0.029 m3·m-3, and 0.023 m3·m-3, respectively, when compared with in-situ measurements. During the thawing period, the average R, MAE, RMSE, and ubRMSE between the downscaled SMAP SSM and in-situ measurements were 0.52, 0.057 m3·m-3, 0.067 m3·m-3, and 0.054 m3·m-3, respectively, compared to 0.45, 0.070 m3·m-3, 0.083 m3·m-3, and 0.060 m3·m-3 for the original SMAP SSM. Thus, the research significantly enhances both the accuracy and spatial resolution of SMAP SSM estimations, underscoring its vital role in advancing hydrological studies within the SRYR.