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 R2 of 0.98 and 0.87, and average RMSE of 0.61℃ and 0.03 m3·m−3, respectively. At the basin scale, DHTC outperformed (Daily: R2 = 0.66, RMSE = 0.75 mm; Monthly: R2 = 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.
Read full abstract