Quadrature errors in the partical derivatives required for the direct recovery of gravity anomalies from satellite observations
The equations of motion of a geodetic satellite in the earth's gravitational field expressed by gravity anomalies require the evaluation, amongst others, of the partial derivatives of the disturbing force with respect to individual gravity anomalies. Data are discussed on how anomaly blocks should be subdivided so that the partial derivatives may be numerically evaluated for each subdivision, and then finally meaned to give the value representative of the whole blocks, with accuracies better than 2 to 3 percent for all blocks. The number of subdivisions is large for the blocks nearest to the satellite subpoint and decreases away from it. The actual values of this spherical distance and the actual subdivision of the mean gravity anomaly blocks was determined numerically for 184 15 deg x 15 deg equal area blocks. Satellite heights above the earth of 400 km, 800 km and 1600 km were considered. The computer times for the suggested scheme were compared with alternative solutions.