Unknown future precipitation is the dominant source of uncertainty for many streamflow forecasts. Numerical weather prediction (NWP) models can be used to generate quantitative precipitation forecasts (QPF) to reduce this uncertainty. The usability and usefulness of NWP model outputs depend on the application time and space scales as well as forecast lead time. For streamflow nowcasting (very short lead times; e.g., 12 h), many applications are based on measured in situ or radar-based real-time precipitation and/or the extrapolation of recent precipitation patterns. QPF based on NWP model output may be more useful in extending forecast lead time, particularly in the range of a few days to a week, although low NWP model skill remains a major obstacle. Ensemble outputs from NWP models are used to articulate QPF uncertainty, improve forecast skill, and extend forecast lead times. Hydrologic prediction driven by these ensembles has been an active research field, although operational adoption has lagged behind. Conversely, relatively little study has been done on the hydrologic component (i.e., model, parameter, and initial condition) of uncertainty in the streamflow prediction system. Four domains of research are identified: selection and evaluation of NWP model–based QPF products, improved QPF products, appropriate hydrologic modeling, and integrated applications.

Skillful and reliable forecasts of seasonal streamflows are highly valuable to water management. In a previous study, we developed a Bayesian joint probability (BJP) modeling approach for seasonal forecasting of streamflows at multiple sites. The approach has been adopted by the Australian Bureau of Meteorology for seasonal streamflow forecasting in Australia. This study extends the applicability of the BJP modeling approach to streams with zero flow occurrences. The aim is to produce forecasts for these streams in the form of probabilities of zero seasonal flows and probability distributions of above zero seasonal flows. We turn a difficult mathematical problem of mixed discrete‐continuous multivariate probability distribution modeling into one of continuous multivariate probability distribution modeling by treating zero flow occurrences as censored data. This paper presents the mathematical formulation and implementation of the modeling approach, methods for forecast verification, and results of a test application to the Burdekin river catchment in northern Queensland, Australia.