The paper presents a numerical method to deconvolve complex body waves into a multiple shock sequence, assuming that all constituent events have identical fault geometry and depth. The far-field source time function is derived as a superposition of ramp functions, with the height and onset time determined by matching synthetic waveforms with observed ones. The method is applied to single and multi-station data, using the 1976 Guatemala earthquake as a test case. The results show that the method can systematically analyze complex event sequences, providing insights into the source mechanisms and heterogeneity of the fault zone. The validity of the model is assessed through comparisons with local strong-motion data, surface break distributions, and macroseismic data. The study highlights the importance of asperities in various seismological problems and offers a useful tool for interpreting complex observed records.The paper presents a numerical method to deconvolve complex body waves into a multiple shock sequence, assuming that all constituent events have identical fault geometry and depth. The far-field source time function is derived as a superposition of ramp functions, with the height and onset time determined by matching synthetic waveforms with observed ones. The method is applied to single and multi-station data, using the 1976 Guatemala earthquake as a test case. The results show that the method can systematically analyze complex event sequences, providing insights into the source mechanisms and heterogeneity of the fault zone. The validity of the model is assessed through comparisons with local strong-motion data, surface break distributions, and macroseismic data. The study highlights the importance of asperities in various seismological problems and offers a useful tool for interpreting complex observed records.