A B S T R A C TIn seismic interpretation and seismic data analysis, it is of critical importance to effectively identify certain geologic formations from very large seismic data sets. In particular, the problem of salt characterization from seismic data can lead to important savings in time during the interpretation process if solved efficiently and in an automatic manner. In this work, we present a novel numerical approach that is able to automatically segmenting or identifying salt structures from a post-stack seismic data set with a minimum intervention from the interpreter. The proposed methodology is based on the recent theory of sparse representation and consists in three major steps: first, a supervised learning assisted by the user which is performed only once, second a segmentation process via unconstrained 1 optimization, and finally a post-processing step based on signal separation. Furthermore, since the second step only depends upon local information at each time, the whole process greatly benefits from parallel computing platforms. We conduct numerical experiments in a synthetic 3D seismic data set demonstrating the viability of our method. More specifically, we found that the proposed approach matches up to 98.53% with respect to the corresponding 3D velocity model available in advance. Finally, in appendixes A and B, we present a convergence analysis providing theoretical guarantees for the proposed method.