Abstract-During the past decade, compressed sensing (CS) has proven to be extremely useful for sparse signal reconstruction. Here the method of CS is applied to fast mass distribution determination based on cold atom gravity gradiometer measurements. Specifically we consider an array of M gradiometers placed around a 2D target area in order to determine the interior mass distribution in cases where the set of (M) sensor measurements under samples the distribution. This was done by assuming that the system's sparsity comes from the mass distributions only having K non-zero masses, which led to a non-orthogonal basis dictionary. This lack of orthogonality caused interesting behaviors, including weakened noise performance, and breaking of the typical CS motivated logarithmic scaling of required M values for larger K values. However, for low K values M scaled as expected. Modifications to the gravity sensor model to promote orthogonality and test its impact on signal recovery postponed the onset of anomalous scaling; suggesting that lack of orthogonality is the primary cause. While CS works for this sparse, but intrinsically ill-posed problem, this sensor system displayed increased noise sensitivity and a smaller upper bound on the size of recoverable K sets. However, while these limitations decrease CS performance, significant improvements over traditional sensing approaches are still possible.