The simultaneous identification of location and source release history in complex river networks is a very complicated ill-posed problem, particularly in a case of multiple unknown pollutant sources with time-varying release pattern. This study presents an innovative method for simultaneous identification of the number, locations and release histories of multiple pollutant point sources in a river network using minimum observation data. Considering two different type of monitoring stations with an adaptive arrangement as well as real-time data collection at those stations and using a reliable numerical flow and transport model, at first the number and suspected reach of presence of pollutant sources are determined. Then the source location and its intensity function is calculated by solving inverse source problem using a geostatistical approach. A case study with three different scenarios in terms of the number, release time and location of pollutant sources are discussed, concerning a river network with unsteady and non-uniform flow. Results showed the capability of the proposed method in identifying of sought source characteristics even in complicated cases with simultaneous activity of multiple pollutant sources.