For decades, how to improve the resolution limit of seismic data has been a concern for seismologists. For this reason, several geoscientists have proposed various methods of improving the bandwidth of the data. Discovering an easy and scientifically reliable means of improving seismic data resolution would undoubtedly help geophysicists interpret more complex details of the subsurface geology. In this study, we transformed the bandlimited time-domain seismic data to the frequency domain using the Fourier analysis method, and a basis pursuit atomic algorithm was applied to decompose the real and imaginary parts of the spectrum into summations of cosines and sines. The resultant reflectivity spectrum (in the frequency domain) was deconvolved by a pre-estimated wavelet spectrum to obtain the true earth�s reflectivity data spectrum and was subsequently extrapolated to beyond the original band limit. The result shows an extended bandwidth from 68 Hz to 161 Hz and 80 Hz to 170 Hz for both synthetic trace model and the main seismic data, respectively. Consequently, this improved vertical resolution of sub-seismic geologic features, such as crevasse splay, levee, barrier bar complex, lagoon inlet channels, alluvial fans, and fluvial channels, and shows subtle facies variations in Inas field. © 2022 by the authors. Licensee MDPI, Basel, Switzerland.