The goal of this study is to constrain the origins layering in seismic velocity structure within cratonic mantle lithosphere (i.e. mid-lithospheric discontinuities [MLDs]). For long-lived stations cratons worldwide, we calculated S-to-P converted phase receiver function stacks using time domain deconvolution and a k-means algorithm select robust, consistent functions. Negative MLDs appear only ...