Inferences of Mantle Viscosity Based on Ice Age Data Sets: The Bias in Radial Viscosity Profiles Due to the Neglect of Laterally Heterogeneous Viscosity Structure. Issue 9 (15th September 2018)

• Record Type: Journal Article
• Authors: Lau, H. C. P.; Austermann, J.; Mitrovica, J. X.; Crawford, O.; Al‐Attar, D.; Latychev, K.
• Is Part Of: Journal of geophysical research. Volume 123:Issue 9(2018)
• Page Start: 7237
• Page End: 7252
• Publication Date: 2018-09-15
• Subjects: Geomagnetism -- Periodicals; Geochemistry -- Periodicals; Geophysics -- Periodicals; Earth sciences -- Periodicals; 551.1
• DOI: 10.1029/2018JB015740
• Languages: English
• ISSNs: 2169-9313

Abstract

Abstract: Inferences of mantle viscosity using glacial isostatic adjustment (GIA) data are hampered by data sensitivity to the space‐time geometry of ice cover. A subset of GIA data is relatively insensitive to this ice history: the Fennoscandian relaxation spectrum (FRS), postglacial decay times in Canada and Scandinavia, and the rate of change of the degree‐2 zonal harmonic of the geopotential ( J ̇ 2 ). These geographically limited data have been inverted to constrain the radial (one‐dimensional [1D]) mantle viscosity profile. We explore potential biases in these 1D inversions introduced by neglecting a three‐dimensional (3D) viscosity structure. We perform 1D Bayesian inversions of synthetic GIA data generated from Earth models with realistic 3D variations in mantle viscosity and lithospheric thickness and compare results to the 1D viscosity profile associated with the 3D model used to generate the synthetics. Differences between these two 1D profiles reflect GIA data resolution and biasing introduced by neglecting, in the inversions, a 3D viscosity structure. We focus on the second issue, demonstrating that the largest bias occurs within the upper mantle (in particular, the transition zone). This remains consistent when varying inversion parameters (e.g., prior/starting models) and the 1D/3D viscosity fields adopted in generating the synthetics. Inversions of individual data sets show 3D biasing increases for data exhibiting shallower (thus more localized) sensitivity to viscosity. Of the data considered herein, inversions of the FRS are subject to the largest bias followed by decay time data. The bias is minimal for J ̇ 2, as its deeper sensitivity is accompanied by broader averaging of structure in radial and lateral directions. Key Points: Synthetic glacial isostatic adjustment data produced using realistic models of 3D mantle viscosity are inverted for a 1D viscosity structure The bias introduced in 1D inversions due to the neglect of 3D viscosity structure is isolated from that due to limited radial resolution The bias peaks within the lower half of the upper mantle and originates from data with shallow sensitivity to mantle viscosity