No abstract

We have developed the open-source software Tesseroids, a set of command-line programs to perform forward modeling of gravitational fields in spherical coordinates. The software is implemented in the C programming language and uses tesseroids (spherical prisms) for the discretization of the subsurface mass distribution. The gravitational fields of tesseroids are calculated numerically using the Gauss-Legendre quadrature (GLQ). We have improved upon an adaptive discretization algorithm to guarantee the accuracy of the GLQ integration. Our implementation of adaptive discretization uses a “stack-based” algorithm instead of recursion to achieve more control over execution errors and corner cases. The algorithm is controlled by a scalar value called the distance-size ratio ([Formula: see text]) that determines the accuracy of the integration as well as the computation time. We have determined optimal values of [Formula: see text] for the gravitational potential, gravitational acceleration, and gravity gradient tensor by comparing the computed tesseroids effects with those of a homogenous spherical shell. The values required for a maximum relative error of 0.1% of the shell effects are [Formula: see text] for the gravitational potential, [Formula: see text] for the gravitational acceleration, and [Formula: see text] for the gravity gradients. Contrary to previous assumptions, our results show that the potential and its first and second derivatives require different values of [Formula: see text] to achieve the same accuracy. These values were incorporated as defaults in the software.

At periods greater than 1000 seconds, Earth's seismic free oscillations have anomalously large amplitude when referenced to the Harvard Centroid Moment Tensor fault mechanism, which is estimated from 300- to 500-second surface waves. By using more realistic rupture models on a steeper fault derived from seismic body and surface waves, we approximated free oscillation amplitudes with a seismic moment (6.5 x 10(22) Newton.meters) that corresponds to a moment magnitude of 9.15. With a rupture duration of 600 seconds, the fault-rupture models represent seismic observations adequately but underpredict geodetic displacements that argue for slow fault motion beneath the Nicobar and Andaman islands.