The seislet transform uses a prediction operator that is connected to the local slope or frequency of seismic events. We have combined the 1D nonstationary seislet transform with empirical-mode decomposition (EMD) in the - domain. We used EMD to decompose data into smoothly variable frequency components for the following 1D seislet transform. The resultant representation showed remarkable sparsity. We developed a detailed algorithm and used a field example to demonstrate the application of the new seislet transform for sparsity-promoting seismic data processing.