The linear sampling method seeks to localize the unknown source of an observed, time-dependent field. The unknown source could be, for example, a scatterer embedded within a medium, or an impulsive excitation such as an earthquake or explosion. The source of the observed field is localized by means of solving the so-called near-field equation and mapping the obtained solutions through an indicator functional over a test region assumed to contain the source. In its current formulation, however, the linear sampling method suffers from an ambiguous time parameter that strongly influences its ability to localize the unknown source. Our paper consists of two fundamental results central to the theoretical understanding of the linear sampling method and its numerical implementation. First, we prove the so-called blowup behavior of solutions to the near-field equation for a general source function that is separable in space and time. Second, we show that the linear sampling method can be formulated such that the ambiguous time parameter is irrelevant. We demonstrate that a dependence of the linear sampling method on the time parameter arises from an incorrect implementation of a convolution-type operator found in the near-field equation. When the operator is implemented correctly, the dependence on the time parameter vanishes. We provide detailed algorithms for efficient and proper implementations of the convolutional operator in both the time and frequency domains. The crucial result of the improved implementations is that they allow the linear sampling method to be completely automated, as one does not need to know the space-time dependence of the unknown source. We demonstrate the effectiveness of the improved time-and frequency-domain implementations using several numerical examples applied to imaging scatterers.