Variable source depth beneath the Indian Ocean geoid low area: Insights from L1 and L2 norm-based scaling power spectrum inversion

Indian Ocean geoid low, situated at the south of the Indian peninsula, is an intriguing feature. Hitherto, several interpretations have been put forth to explain the genesis and causative depths of the anomaly. In this work, we investigate the scaling spectral inversion approach utilizing gravity data to estimate the causative depth and nature of the source from this region. In practice, crustal or sub-surface heterogeneities follow scaling or fractal characteristics, and comparative studies across various regions suggest that the scaling power spectrum method provides better insights into the sub-surface source distribution. An optimum fitting and comparison between the modeled and observed power spectrum have been attained by utilizing two norm penalties e.g. L1 and L2 norms for various combinations of depth and scaling exponents. Equally important, the L1 norm misfit function yields better results when compared to the L2 norm function. We propose that estimated source depths (̴80 km and 120 km, respectively) are perfect examples of a half-space model with a source following scaling or fractal characteristics (scaling exponent value of ̴1.78) along the upper mantle boundary. However, the scaling spectrum windowing analysis of the studied region reveals a variable depth to long wavelength sources as 1014 km, 431 km, and 94 km, which conform well with the recently published depth range. Further, the weight of low and high wavenumbers indicates that the causative source depths might be due to the above-stated depth range rather than a single model fit. The study has demonstrated a good case for wider use of L1 norm-based power spectrum inversion for source depth estimation from potential field data. We further suggest that the causative source for the Indian Ocean geoid low may extend deeper than upper mantle discontinuities.