Process diagnostics of snowmelt runoff in global hydrological and land surface models – Part 1: A systematic evaluation across basins of increasing complexity
Abstract. Accurate simulation of snowmelt runoff (SMR) is critical for water resource management. However, despite the abundance of global hydrological models, little is known about their SMR performance. This study presents a comprehensive evaluation of SMR across 15 state-of-the-art large-scale models and runoff products by focusing on their biases in first-order indices, i.e., the total volume (Qsum), peak flow (Qmax), and centroid timing (CTQ) of runoff in the snowmelt period. Then by introducing 1455 snow-dominated basins with diverse topography and vegetation complexities, we further proposed a novel model robustness metric to test how different models perform under increasing basin complexity, thereby allowing for a quantification on how they adapt to complex environmental conditions. Our results reveal that (1) most models exhibit underestimated Qsum and Qmax and predict CTQ too early. These biases are particularly pronounced in regions such as the western United States, northern Europe, and northeastern China. (2) Model biases systematically increase with basin complexity, with CTQ exhibiting strong sensitivity to mean elevation and topographic variability, while Qsum and Qmax being shaped more by mean elevation and the diversity of vegetation types in the basin. (3) The robustness assessment further shows that observation-constrained runoff products exhibit the most outstanding performance (i.e., low biases and strong adaptability to stern conditions), followed by the hydrological and land surface models. Notably, while global hydrological models generally exhibit stronger robustness in simulating Qsum and Qmax, land surface models show a clear advantage in simulating CTQ, highlighting their structural strength in capturing melt timing rather than runoff magnitude. This study provides a large-sample benchmark for SMR evaluation and complements existing model assessment approaches by examining model performance across basin complexity gradients, offering useful insights for future model development and uncertainty reduction.