The topological derivative method is used to solve a pollution sources reconstruction problem governed by a steady-state convection-diffusion equation. To be more precise, we are dealing with a shape optimization problem which consists of reconstruction of a set of pollution sources in a fluid medium by measuring the concentration of the pollutants within some subregion of the reference domain. The shape functional measuring the misfit between the known data and solution of the state equation is minimized with respect to a set of ball-shaped geometrical subdomains representing the pollution sources. The necessary conditions for optimality are derived with the help of the topological derivative method which consists in expanding the shape functional asymptotically and then truncate it up to the second order term. The resulting expression is trivially minimized with respect to the parameters under consideration which leads to a noniterative second-order reconstruction algorithm. Two different cases are considered. Firstly, when the velocity of the leakages is given and we reconstruct the support of the unknown sources, including their locations and sizes. In the second case, we consider the size of the pollution sources to be known and find out the mean velocity of the leakages and their locations. Numerical examples are presented showing the capability of the proposed algorithm in reconstructing multiple pollution sources in both cases. [ABSTRACT FROM AUTHOR]