Simulation of aftershock using a wavelet-based description of the time frequency characteristics of a class of earthquake accelerograms and proposes three-parameter models for the modulation of wavelet coefficient envelopes of different frequency bands.