Bayesian Inference of Mantle Viscosity From Whole‐Mantle Density Models

Abstract The long‐wavelength geoid is sensitive to Earth's mantle density structure as well as radial variations in mantle viscosity. We present a suite of inversions for the radial viscosity profile using whole‐mantle models that jointly constrain the variations in density, shear‐ and compress...

Full description

Bibliographic Details
Main Authors: M. L. Rudolph, P. Moulik, V. Lekić
Format: Article
Language:English
Published: Wiley 2020-11-01
Series:Geochemistry, Geophysics, Geosystems
Online Access:https://doi.org/10.1029/2020GC009335
_version_ 1797637919785091072
author M. L. Rudolph
P. Moulik
V. Lekić
author_facet M. L. Rudolph
P. Moulik
V. Lekić
author_sort M. L. Rudolph
collection DOAJ
description Abstract The long‐wavelength geoid is sensitive to Earth's mantle density structure as well as radial variations in mantle viscosity. We present a suite of inversions for the radial viscosity profile using whole‐mantle models that jointly constrain the variations in density, shear‐ and compressional‐wavespeeds using full‐spectrum tomography. We use a Bayesian approach to identify a collection of viscosity profiles compatible with the geoid, while enabling uncertainties to be quantified. Depending on tomographic model parameterization and data weighting, it is possible to obtain models with either positive‐ or negative‐buoyancy in the large low shear velocity provinces. We demonstrate that whole‐mantle density models in which density and VS variations are correlated imply an increase in viscosity below the transition zone, often near 1,000 km. Many solutions also contain a low‐viscosity channel below 650 km. Alternatively, models in which density is less‐correlated with VS—which better fit normal mode data—require a reduced viscosity region in the lower mantle. This feature appears in solutions because it reduces the sensitivity of the geoid to buoyancy variations in the lowermost mantle. The variability among the viscosity profiles obtained using different density models is indicative of the strong nonlinearities in modeling the geoid and the limited resolving power of the geoid kernels. We demonstrate that linearized analyses of model resolution do not adequately capture the posterior uncertainty on viscosity. Joint and iterative inversions of viscosity, wavespeeds, and density using seismic and geodynamic observations are required to reduce bias from prior assumptions on viscosity variation and scalings between material properties.
first_indexed 2024-03-11T12:55:10Z
format Article
id doaj.art-078d13e4d53e4444adfd0c52429726ea
institution Directory Open Access Journal
issn 1525-2027
language English
last_indexed 2024-03-11T12:55:10Z
publishDate 2020-11-01
publisher Wiley
record_format Article
series Geochemistry, Geophysics, Geosystems
spelling doaj.art-078d13e4d53e4444adfd0c52429726ea2023-11-03T17:01:17ZengWileyGeochemistry, Geophysics, Geosystems1525-20272020-11-012111n/an/a10.1029/2020GC009335Bayesian Inference of Mantle Viscosity From Whole‐Mantle Density ModelsM. L. Rudolph0P. Moulik1V. Lekić2Department of Earth and Planetary Sciences University of California Davis CA USADepartment of Geology University of Maryland College Park MD USADepartment of Geology University of Maryland College Park MD USAAbstract The long‐wavelength geoid is sensitive to Earth's mantle density structure as well as radial variations in mantle viscosity. We present a suite of inversions for the radial viscosity profile using whole‐mantle models that jointly constrain the variations in density, shear‐ and compressional‐wavespeeds using full‐spectrum tomography. We use a Bayesian approach to identify a collection of viscosity profiles compatible with the geoid, while enabling uncertainties to be quantified. Depending on tomographic model parameterization and data weighting, it is possible to obtain models with either positive‐ or negative‐buoyancy in the large low shear velocity provinces. We demonstrate that whole‐mantle density models in which density and VS variations are correlated imply an increase in viscosity below the transition zone, often near 1,000 km. Many solutions also contain a low‐viscosity channel below 650 km. Alternatively, models in which density is less‐correlated with VS—which better fit normal mode data—require a reduced viscosity region in the lower mantle. This feature appears in solutions because it reduces the sensitivity of the geoid to buoyancy variations in the lowermost mantle. The variability among the viscosity profiles obtained using different density models is indicative of the strong nonlinearities in modeling the geoid and the limited resolving power of the geoid kernels. We demonstrate that linearized analyses of model resolution do not adequately capture the posterior uncertainty on viscosity. Joint and iterative inversions of viscosity, wavespeeds, and density using seismic and geodynamic observations are required to reduce bias from prior assumptions on viscosity variation and scalings between material properties.https://doi.org/10.1029/2020GC009335
spellingShingle M. L. Rudolph
P. Moulik
V. Lekić
Bayesian Inference of Mantle Viscosity From Whole‐Mantle Density Models
Geochemistry, Geophysics, Geosystems
title Bayesian Inference of Mantle Viscosity From Whole‐Mantle Density Models
title_full Bayesian Inference of Mantle Viscosity From Whole‐Mantle Density Models
title_fullStr Bayesian Inference of Mantle Viscosity From Whole‐Mantle Density Models
title_full_unstemmed Bayesian Inference of Mantle Viscosity From Whole‐Mantle Density Models
title_short Bayesian Inference of Mantle Viscosity From Whole‐Mantle Density Models
title_sort bayesian inference of mantle viscosity from whole mantle density models
url https://doi.org/10.1029/2020GC009335
work_keys_str_mv AT mlrudolph bayesianinferenceofmantleviscosityfromwholemantledensitymodels
AT pmoulik bayesianinferenceofmantleviscosityfromwholemantledensitymodels
AT vlekic bayesianinferenceofmantleviscosityfromwholemantledensitymodels