Abstract. In many simulations of historical daily streamflow distributional bias arising from the distributional properties of residuals, however small, has been noted. This bias often presents itself as an underestimation of high streamflow and an overestimation of low streamflow. Here, 1168 streamgages across the conterminous United States having at least 14 complete water years of daily data between October 01, 1980, and September 30, 2013, are used to explore a method for rescaling simulated streamflow to correct the distributional bias. Based on an existing approach that separates the simulated streamflow into com-5 ponents of timing and magnitude, the timing component is converted into simulated nonexceedance probabilities and rescaled to new volumes using an independently estimated flow-duration curve (FDC). In this study, this method is applied to a pooled ordinary kriging simulation of daily streamflow coupled with FDCs estimated by regional regression on basin characteristics.The improvement in the representation of high and low streamflows is correlated with the accuracy and unbiasedness of the estimated FDC. The method is verified by using an idealized case, though, with the introduction of regionally regressed FDCs 10 developed for this study, the method is only useful overall for the upper tails, which are more accurately and unbiasedly estimated than the lower tails. It remains for future work to determine how accurate the estimated FDCs need to be to be useful for bias correction without unduly reducing accuracy. In addition to its potential efficacy for distributional bias correction, this methodology also represents a generalization of nonlinear spatial interpolation of daily streamflow using FDCs. Rather than relying on single index stations as is commonly done to reflect streamflow timing, this approach leverages geostatistical tools 15 to allow a region of neighbors to reflect streamflow timing.