An algorithm to calculate the gravity anomaly of sedimentary basins with exponential density-depth relationships

We derive wavenumber domain expressions to calculate the gravity anomaly of a body with irregular bounding surfaces and an exponential density-depth relationship. We apply the method to sedimentary basins, which commonly have this type of geometry and density distribution. The mathematical formulation also allows the exponential density-depth relationship to be measured from an arbitrary irregular surface rather than the top surface. Using this arrangement, the gravity anomaly of exhumed sedimentary basins can be predicted if the amount of eroded section can be estimated. The corresponding inverse algorithms are also derived. Examples of the use of the forward algorithms, from the Galicia Interior Basin and the Central Irish Sea Basin, are used to illustrate these methods.