harmonica.ellipsoid_gravity

harmonica.ellipsoid_gravity(coordinates, ellipsoids)

Forward model gravity fields of ellipsoids.

Compute the gravity acceleration components for an ellipsoidal body at specified observation points.

Warning

The vertical direction points upwards, i.e. positive and negative values of upward represent points above and below the surface, respectively. But, g_z is the downward component of the gravitational acceleration so that positive density contrasts produce positive anomalies.

Important

The gravity acceleration components are returned in mGal (\(\text{m}/\text{s}^2\)).

Parameters:

coordinateslist of array

List of arrays containing the easting, northing and upward coordinates of the computation points defined on a Cartesian coordinate system. All coordinates should be in meters.

ellipsoidellipsoid or list of ellipsoids

Ellipsoidal body represented by an instance of harmonica.TriaxialEllipsoid, harmonica.ProlateEllipsoid, or harmonica.OblateEllipsoid, or a list of them.

Returns:

g_e, g_n, g_zarray

Easting, northing and downward component of the gravity acceleration in mGal.

References

[Clark1986] [Takahashi2018]

For derivations of the equations, and methods used in this code.