This paper considers an inverse problem for a partial differential equation to identify a pollution point source in a watershed. The mathematical model of the problem is a weakly coupled system of two linear parabolic equations for the concentrations u(x, t) and v(x, t) with an unknown point source F(x, t) = A( t)δ(x- s) related to the concentration u(x, t), where s is the point source location and A(t) is the amplitude of the pollution point source. Assuming that source F becomes inactive after time T*, it is proved that it can be uniquely determined by the indirect measurements { v(0, t), v( a, t), v( b, t), v( l, t), 0 〈 t ≤ T, T* 〈 T}, and, thus, the local Lipschitz stability for this inverse source problem is obtained. Based on the proof of its uniqueness, an inversion scheme is presented to determine the point source. Finally, two numerical examples are given to show the feasibility of the inversion scheme.