Three-dimensional gravity inversion for Moho depth at rifted continental margins incorporating a lithosphere thermal gravity anomaly correction

This paper describes a method for determining Moho depth, lithosphere thinning factor (gamma = 1 - 1/beta) and the location of the ocean-continent transition at rifted continental margins using 3-D gravity inversion which includes a correction for the large negative lithosphere thermal gravity anomaly within continental margin lithosphere. The lateral density changes caused by the elevated geotherm in thinned continental margin and adjacent ocean basin lithosphere produce a significant lithosphere thermal gravity anomaly which may be in excess of -100 mGal, and for which a correction must be made in order to determine Moho depth accurately from gravity inversion. We describe a method of iteratively calculating the lithosphere thermal gravity anomaly using a lithosphere thermal model to give the present-day temperature field from which we calculate the lithosphere thermal density and gravity anomalies. For continental margin lithosphere, the lithosphere thermal perturbation is calculated from the lithosphere thinning factor (gamma = 1 - 1/beta) obtained from crustal thinning determined by gravity inversion and breakup age for thermal re-equilibration time. For oceanic lithosphere, the lithosphere thermal model used to predict the lithosphere thermal gravity anomaly may be conditioned using ocean isochrons from plate reconstruction models to provide the age and location of oceanic lithosphere. A correction is made for crustal melt addition due to decompression melting during continental breakup and seafloor spreading. We investigate the sensitivity of the lithosphere thermal gravity anomaly and the predicted Moho depth from gravity inversion at continental rifted margins to the methods used to calculate and condition the lithosphere thermal model using both synthetic models and examples from the North Atlantic.