The accuracy assessment of lithospheric density models

Abstract

The Earth’s synthetic gravitational and density models can be used to validate numerical procedures applied for global (or large-scale regional) gravimetric forward and inverse modeling. Since the Earth’s lithospheric structure is better constrained by tomographic surveys than a deep mantle, most existing 3D density models describe only a lithospheric density structure, while 1D density models are typically used to describe a deep mantle density structure below the lithosphere-asthenosphere boundary. The accuracy of currently available lithospheric density models is examined in this study. The error analysis is established to assess the accuracy of modeling the sub-lithospheric mantle geoid while focusing on the largest errors (according to our estimates) that are attributed to lithospheric thickness and lithospheric mantle density uncertainties. Since a forward modeling of the sub-lithospheric mantle geoid also comprises numerical procedures of adding and subtracting gravitational contributions of similar density structures, the error propagation is derived for actual rather than random errors (that are described by the Gauss’ error propagation law). Possible systematic errors then either lessen or sum up after applying particular corrections to a geoidal geometry that are attributed to individual lithospheric density structures (such as sediments) or density interfaces (such as a Moho density contrast). The analysis indicates that errors in modeling of the sub-lithospheric mantle geoid attributed to lithospheric thickness and lithospheric mantle density uncertainties could reach several hundreds of meters, particularly at locations with the largest lithospheric thickness under cratonic formations. This numerical finding is important for the calibration and further development of synthetic density models of which mass equals the Earth’s total mass (excluding the atmosphere). Consequently, the (long-to-medium wavelength) gravitational field generated by a synthetic density model should closely agree with the Earth’s gravitational field.