Identifying and Interpreting Hydrological Model Structural Nonstationarity Using the Bayesian Model Averaging Method

Understanding hydrological nonstationarity under climate change is important for runoff prediction and it enables more robust decisions. Regarding the multiple structural hypotheses, this study aims to identify and interpret hydrological structural nonstationarity using the Bayesian Model Averaging (BMA) method by (i) constructing a nonstationary model through the Bayesian weighted averaging of two lumped conceptual rainfall-runoff (RR) models (the Xinanjiang and GR4J model) with time-varying weights; and (ii) detecting the temporal variation in the optimized Bayesian weights under climate change conditions. By combining the BMA method with period partition and time sliding windows, the efficacy of adopting time-varying model structures is investigated over three basins located in the U.S. and Australia. The results show that (i) the nonstationary ensemble-averaged model with time-varying weights surpasses both individual models and the ensemble-averaged model with time-invariant weights, improving NSE[Q] from 0.04 to 0.15; (ii) the optimized weights of Xinanjiang model increase and that of GR4J declines with larger precipitation, and vice versa; (iii) the change in the optimized weights is proportional to that of precipitation under monotonic climate change, as otherwise the mechanism changes significantly. Overall, it is recommended to adopt nonstationary structures in hydrological modeling.