Unraveling the upper mantle heterogeneity from integrated multi-observable inversions

The lateral and vertical thermochemical heterogeneity in the mantle is a long standing question in geodynamics. The forces that control mantle flow and therefore Plate Tectonics arise from the density and viscosity lateral and vertical variations. A common approach to estimate the density field for geodynamical purposes is to simply convert seismic tomography anomalies sometimes assuming constraints from mineral physics. Such converted density field does not match in general with the observed gravity field, typically predicting anomalies the amplitudes of which are too large. Knowledge on the lateral variations in lithospheric density is essential to understand the dynamic/residual isostatic components of the Earth’s topography linking deep and surface processes. The cooling of oceanic lithosphere, the bathymetry of mid oceanic ridges, the buoyancy and stability of continental cratons or the thermochemical structure of mantle plumes are all features central to Plate Tectonics that are dramatically related to mantle temperature and composition.

Conventional methods of seismic tomography, topography and gravity data analysis constrain distributions of seismic velocity and density at depth, all depending on temperature and composition of the rocks within the Earth. However, modelling and interpretation of multiple data sets provide a multifaceted image of the true thermochemical structure of the Earth that needs to be appropriately and consistently integrated. A simple combination of gravity, petrological and seismic models alone is insufficient due to the non-uniqueness and different sensitivities of these models, and the internal consistency relationships that must connect all the intermediate parameters describing the Earth involved. In fact, global Earth models based on different observables often lead to rather different, even contradictory images of the Earth.

 Here we present a new global thermochemical model of the lithosphere-upper mantle (WINTERC-grav) constrained by state-of-the-art global waveform tomography, satellite gravity (geoid and gravity anomalies and gradiometric measurements from ESA's GOCE mission), surface elevation and heat flow data. WINTERC-grav is based upon an integrated geophysical-petrological approach where all relevant rock physical properties modelled (seismic velocities and density) are computed within a thermodynamically self-consistent framework allowing for a direct parameterization of the temperature and composition variables.