Abstract.Poorly monitored river flows in many regions of the world have been hindering our ability to accurately estimate global water usage as well as the budgets and variability of the global water cycle. In-situ gauging sites, as well as a number of satellite-based systems, make observations of river discharge throughout the globe; however, these observations are often sparse due to, e.g., the sampling frequencies of sensors or a lack of reporting. Recently, efforts have been made to develop methods to integrate these discrete observations to gain a better understanding of the underlying processes. This paper 10 presents an application of a fixed interval Kalman smoother based model, called Inverse Streamflow Routing (ISR), to generate spatially and temporally continuous river discharge fields from discrete observations. The method propagates the observed information across all reachable parts of the river network (up/downstream from gauging point) and all reachable times (before/after observation time) using a two-sweep procedure that first propagates information backward in time to the furthest upstream locations (inverse routing) and then propagates it forward in time to the furthest downstream locations (forward 15 routing). The key advantages of this approach are that it (1) maintains all the physical consistencies embodied by a diffusive wave routing model (flow confluence relationships on the river network and the resulting mass balance, wave velocity and diffusivity), (2) updates the lateral influx (runoff) at pixel level (furthest upstream) to guarantee exhaustive propagation of observed information, (3) works both with a first guess of initial river discharge conditions from a routing model (assimilation) and without a first guess (pure interpolation of observations). Two sets of experiments are carried out under idealized 20 conditions as well as real-world conditions provided by U.S. Geological Survey observations. Results show that the method can effectively reproduce the spatial and temporal dynamics of river discharge in each of the experiments presented. The performance is driven by the density of the gauge network as well as the quality of the data being assimilated. We find that when assimilating the actual USGS observations, the performance decreases relative to our idealized scenario; however, we are still able to produce an improved discharge product at each validation site. With further testing as well as global application, 25 ISR may prove to be a useful method for extending our current network of global river discharge observations.