Reverse simulation was carried out to specify the point source of pollutants in closed waters. If pollutants flow out into waters, their source must be specified as soon as possible to take a quick measure against the pollution problem. Actually, a wide variety of pollutants have been released into the aquatic environment, some with catastrophic consequences for aquatic life or man. Numerical simulation is one of the powerful tools to specify the point source of pollutants. An inverse trajectory analysis has been used in the case when an advective process of pollutants is dominant. However, the inverse trajectory analysis does not take the diffusion of pollutants into account. Several analytical techniques are applicable only to simple current fields. The present study proposes reverse simulation to specify the point source of pollutants. The basic equation of the reverse simulation is given by changing the positive time derivative term of the advection-diffusion equation of pollutants into a negative term. First, numerical simulation was executed in a forward direction. The results of water current velocities were preserved for the following reverse simulation. Assuming that the pollutants are subject to the surrounding water current and turbulence like a tracer, the advection and diffusion processes of the tracer could be solved for obtaining the initial condition of reverse simulation, and for the following comparison with the results of reverse simulation. The predicted result of water current velocities in the forward simulation was given to calculate the advection-diffusion equation of the tracer. One of the major problems of reverse simulation is the instability of numerical simulation. A Gaussian filter was used for the spatial distribution of the tracer or of the flux of the tracer to eliminate the numerical instability, and the optimum filter width was discussed. As a result, the instability of reverse simulation was suppressed by applying the Gaussian filter for the spatial distribution of the tracer. However, the concentration of the tracer was not condensed in comparison with the result of the tracer in the forward simulation. When the Gaussian filter was used for the special distribution of the flux of the tracer, the accuracy of the prediction of the point source was improved. This is because the high frequency variation was suppressed, keeping the low frequency variation. The concentration of the tracer was more condensed with smaller filter width in both cases. The future studies are to determine the filter width adequately and to modify the difference equation of the Gaussian filter.