Abstract. Profile soil moisture (SM) in mountainous areas are significant in water resources management and ecohydrological studies of downstream arid watersheds. Satellite products are useful in providing spatially distributed SM information, but only have limited penetration depth (e.g. top 5 cm). In contrast, in situ observations can provide multi-depth measurements, but only with limited spatial coverage. Spatially continuous estimates of subsurface SM can be obtained from surface observations using statistical methods, but this requires sufficient coupling strength among surface and subsurface SM. This study evaluates methods to calculate subsurface SM from surface SM and an application to the satellite SM product based on a SM observation network in the Qilian Mountains (China) established since 2013. First, we used cross-correlation to analyze the coupling strength among surface (0–10 cm) and subsurface (10–20, 20–30, 30–50, 50–70 cm, and profile of 0–70 cm) SM. Our results indicated an overall strong coupling among surface and subsurface SM in this study area. Afterwards, three different methods were tested to estimate subsurface SM from in-situ surface SM: the exponential filter (ExpF), artificial neural networks (ANN) and cumulative distribution function matching (CDF) methods. The results showed that both ANN and ExpF methods were able to provide accurate estimates of subsurface soil moisture at 10–20 cm, 20–30 cm, and for the profile of 0–70 cm using surface (0–10 cm) soil moisture only. Specifically, the ANN method had the lowest estimation error (RSR) of 0.42, 0.62 and 0.49 for depths of 15 and 25 cm and profile SM, respectively, while the ExpF method best captured the temporal variation of subsurface soil moisture. Furthermore, it could be shown that the performance of the profile SM estimation was not significantly lower with using an area-generalized Topt (optimum T) compared to the station-specific Topt. In a final step, the ExpF method was applied to the satellite SM product (Soil Moisture Active Passive Level 3: SMAP_L3) to estimate profile SM, and the resulting profile SM was compared to in situ observations. The results showed that the ExpF method was able to estimate profile SM from SMAP_L3 surface products with reasonable accuracy (median R of 0.718). It was also found that the combination of ExpF method and SMAP_L3 surface product can significantly improve the estimation of profile SM in the mountainous areas in comparison to the SMAP_L4 root zone product. Overall, it was concluded that the ExpF method is able to estimate profile SM using SMAP surface products in the Qilian Mountains.