Accurately quantifying the impact of permafrost degradation and soil freeze-thaw cycles on hydrological processes while minimizing the reliance on observational data are challenging issues in hydrological modeling in cold regions. In this study, we developed a modular distributed hydro-thermal coupled hydrological model for cold regions (DHTC) that features a flexible structure. The DHTC model couples heat-water transport processes by employing the conduction-advection heat transport equation and Richard equation considering ice-water phase change. Additionally, the DHTC model integrates the influence of organic matter into the hydrothermal parameterization scheme and includes a subpermafrost module based on the flow duration curve analysis to estimate cold-season streamflow sustained by subpermafrost groundwater. Moreover, we incorporated energy consumption due to ice phase changes to the available energy, enhancing the accuracy of evaporation estimation in cold regions. A comprehensive evaluation of the DHTC model was conducted. At the point scale, the DHTC model accurately replicates daily soil temperature and moisture dynamics at various depths, achieving average R-2 of 0.98 and 0.87, and average RMSE of 0.61degree celsius and 0.03 m(3)m(-3), respectively. At the basin scale, DHTC outperformed (Daily: R-2 = 0.66, RMSE = 0.75 mm; Monthly: R-2 = 0.90, RMSE = 15.7 mm) the GLDAS/FLDAS Noah, GLDAS/VIC, and PML-V2 models in evapotranspiration simulation. The DHTC model also demonstrated reasonable performance in simulating daily (NSE = 0.70, KGE = 0.84), monthly (NSE = 0.86, KGE = 0.90), and multi-year monthly (NSE = 0.97, KGE = 0.93) streamflow in the Source Regions of Yangtze River. DHTC also successfully reproduced the snow depth in basin-averaged time series and spatial distributions (RMSE = 0.86 cm). The DHTC model provides a robust tool for exploring the interactions between permafrost and hydrological processes, and their responses to climate change.

来源平台：JOURNAL OF HYDROLOGY