Atmospheric moisture holding capacity increases with temperature by about 7% per °C according to the Clausius‐Clapeyron relationship. Thermodynamically then, precipitation intensity should exponentially intensify and thus worsen flood conditions as the climate warms. However, regional and global analyses often report a nonmonotonic (hook) scaling of precipitation and runoff, in which extremes strengthen with rising temperature up to a maximum or peak point (T_pp) and decline thereafter. The underlying cause of this hook structure is not yet well‐understood, and whether it may shift and/or regulate storm runoff extremes under anthropogenic climate warming remains unknown. Here, we examine temperature scaling of precipitation and storm runoff extremes under different climate conditions using observations and large ensemble hydroclimatic simulations over mainland China. In situ observations suggest a spatially homogeneous, negative response of relative humidity to warming climates over 34.6% of the land area, and the remaining hook‐dominated regions usually show a colder T_pp than that of precipitation or storm runoff extremes. The precipitation and streamflow series over mainland China's catchments throughout the 21st century are projected by a model cascade chain under a high‐end emission scenario (RCP 8.5), which involves 31 CMIP5 climate models, 11 CMIP6 climate members, a daily bias correction method, and four lumped conceptual hydrological models. The CMIP5 ensemble projects that the hook structures shift toward warmer temperature bins, resulting in 10%–30% increases in storm runoff extremes over mainland China, while the CMIP6 ensemble projects more severe flood conditions in future warming climates.

Quantification of the inherent uncertainty in hydrologic forecasting is essential for flood control and water resources management. The existing approaches, such as Bayesian model averaging (BMA), hydrologic uncertainty processor (HUP), copula-BMA (CBMA), aim at developing reliable probabilistic forecasts to characterize the uncertainty induced by model structures. In the probability forecast framework, these approaches either assume the probability density function (PDF) to follow a certain distribution, or are unable to reduce bias effectively for complex hydrological forecasts. To overcome these limitations, a copula Bayesian processor associated with BMA (CBP-BMA) method is proposed with ensemble lumped hydrological models. Comparing with the BMA and CBMA methods, the CBP-BMA method relaxes any assumption on the distribution of conditional PDFs. Several evaluation criteria, such as containing ratio, average bandwidth and average deviation amplitude of probabilistic application, are utilized to evaluate the model performance. The case study results demonstrate that the CBP-BMA method can improve hydrological forecasting precision with higher cover ratios more than 90%, which are increased by 4.4% and 3.2%, 2.2% and 1.7% over those of BMA and CBMA during the calibration and validation periods, respectively. The proposed CBP-BMA method provides an alternative approach for uncertainty estimation of hydrological multi-model forecasts.