Solving the bound dependence and staircasing problems in 3D potential-field sparse inversions using a mixed Lp-norm regularization

Sparse inversions have proven to be useful for interpreting potential-field data because the recovered models are characterized by sharp boundaries, compact features and elevated values, compared with conventional smoothness-based inversion results. However, several open problems remain to be addressed, including the boundary dependence and staircasing problems. The former results in recovered anomalous values being close to the upper bound, while the latter leads to recovered anomalous bodies with either horizontal or vertical boundaries. We have developed a mixed Lp norm regularization strategy to address these two problems. Inversion results based on two synthetic examples confirm the validity of our approach. The data and codes in this upload allow users to reproduce our results in the manuscript with the same title submitted to Geophysical Prospecting.