It is important to identify source information after a river chemical spill incident occurs. Among various source inversion approaches, a Bayesian-based framework is able to directly characterize inverse uncertainty using a probability distribution and has recently become of interest. However, the literature has not reported its application to actual spill incidents, and many aspects in practical use have not yet been clearly illustrated, e.g., feasibility for large scale pollution incidents, algorithm parameters, and likelihood functions. This work deduced a complete modular-Bayesian approach for river chemical spills, which combined variance assumptions on a pollutant concentration time series with Adaptive-Metropolis sampling. A retrospective case study was conducted based on the ‘landmark’ spill incident in China, the Songhua River nitrobenzene spill of 2005. The results show that release mass, place, and moment were identified with biases of −26.9%, −7.9%, and 16.9%, respectively. Inverse uncertainty statistics were also quantified for each source parameter. Performance, uncertainty sources, and future work are discussed. This study provides an important real-life case to demonstrate the usefulness of the modular-Bayesian approach in practice and provides valuable references for the setting of parameters for the sampling algorithm and variance assumptions.