Conventional Back-Projection Imaging (BPI) methods determine the origin time of an event from the peak of the Maximum Brightness Function curve, however, this assumption might not be completely valid for long duration seismic events composed of continuous bursts of energy, which may produce non-pronounced peaks on the curve due to varying levels of noise and waveform complexity.For such adverse conditions, an adaptive centroid of the energy could be considered a better alternative to describe the origin (x0, y0, z0, t0) of these events.Here we propose the Brightness PDF-based Amplitude Stacking (PbAS) method, a 4D-Centroid location method derived from the Back-Projection Imaging family, able to compute simultaneously the spatial and temporal location of a wide variety of seismic sources.We generate 100 model realizations by perturbing an initial layer-cake velocity model.For each model we perform the event location using a conventional BPI method and PbAS method.The PbAS method produced lower dispersion results and exhibited a better capability to handle higher levels of velocity uncertainty compared to the conventional BPI method.Furthermore, this centroid based method is not as dependent on discretization as most conventional BPI methods, which enable the PbAS method to achieve enough resolution without excessive mesh refinement.