Geodynamic tomography: constraining upper-mantle deformation patterns from Bayesian inversion of surface waves J. K. Magali, T Bodin, N Hedjazian, H Samuel, S Atkins
To cite this version:
J. K. Magali, T Bodin, N Hedjazian, H Samuel, S Atkins. Geodynamic tomography: constraining upper-mantle deformation patterns from Bayesian inversion of surface waves. Geophysical Journal International, Oxford University Press (OUP), 2021, 224 (3), pp.2077 - 2099. 10.1093/gji/ggaa577. hal-03189035
HAL Id: hal-03189035 https://hal.archives-ouvertes.fr/hal-03189035 Submitted on 2 Apr 2021
HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Geophys. J. Int. (2021) 224, 2077–2099 doi: 10.1093/gji/ggaa577 Advance Access publication 2020 December 03 GJI Seismology
Geodynamic tomography: constraining upper-mantle deformation patterns from Bayesian inversion of surface waves
J. K. Magali,1 T. Bodin ,1 N. Hedjazian ,1 H. Samuel 2 and S. Atkins3 1UCBL, CNRS, LGL-TPE, Universite´ de Lyon, 69622 Villeurbanne, France. E-mail: [email protected] 2Institut de Physique du Globe de Paris, Universite´ de Paris, CNRS, F-75005 Paris, France 3Laboratoire de Geologie,´ Ecole Normale Superieure,´ PSL Res. Univ, 75005 Paris, France Downloaded from https://academic.oup.com/gji/article/224/3/2077/6019874 by INFU BIBLIO PLANETS user on 02 April 2021
Accepted 2020 December 1. Received 2020 November 21; in original form 2020 September 23
SUMMARY In the Earth’s upper mantle, seismic anisotropy mainly originates from the crystallographic preferred orientation (CPO) of olivine due to mantle deformation. Large-scale observation of anisotropy in surface wave tomography models provides unique constraints on present- day mantle flow. However, surface waves are not sensitive to the 21 coefficients of the elastic tensor, and therefore the complete anisotropic tensor cannot be resolved independently at every location. This large number of parameters may be reduced by imposing spatial smoothness and symmetry constraints to the elastic tensor. In this work, we propose to regularize the tomographic problem by using constraints from geodynamic modelling to reduce the number of model parameters. Instead of inverting for seismic velocities, we parametrize our inverse problem directly in terms of physical quantities governing mantle flow: a temperature field, and a temperature-dependent viscosity. The forward problem consists of three steps: (1) calculation of mantle flow induced by thermal anomalies, (2) calculation of the induced CPO and elastic properties using a micromechanical model, and (3) computation of azimuthally varying surface wave dispersion curves. We demonstrate how a fully nonlinear Bayesian inversion of surface wave dispersion curves can retrieve the temperature and viscosity fields, without having to explicitly parametrize the elastic tensor. Here, we consider simple flow models generated by spherical temperature anomalies. The results show that incorporating geodynamic constraints in surface wave inversion help to retrieve patterns of mantle deformation. The solution to our inversion problem is an ensemble of models (i.e. thermal structures) representing a posterior probability, therefore providing uncertainties for each model parameter. Key words: Inverse theory; Probability distributions; Seismic anisotropy; Seismic tomog- raphy; Surface wave and free oscillations.
seismic data, tomographers have produced detailed models of az- 1 INTRODUCTION imuthal anisotropy (e.g. Debayle et al. 2005; Deschamps et al. Seismic anisotropy reveals key insights into the Earth’s interior 2008; Adam & Lebedev 2012; Yuan & Beghein 2013, 2014), and structure and dynamics. In the upper mantle, the propagation of radial anisotropy (e.g., Plomerova´ et al. 2002;Lebedevet al. 2006; seismic waves appears to be anisotropic, which has generally been Nettles & Dziewonski´ 2008; Chang et al. 2014, 2015). Numerous associated with the preferred alignment of mantle minerals (Nico- studies have inverted dispersion curves by minimizing the differ- las & Christensen 1987; Montagner 1994). This so-called intrinsic ence between observed and synthetic phase and/or group velocities, anisotropy relates to the strain history induced by regional-scale proving that they can effectively constrain the depth dependence convection and is observable with various seismological tools, in- of anisotropy (e.g., Montagner & Tanimoto 1990; Ritzwoller et al. cluding surface waves. 2002). Seismic anisotropy can be described with 21 independent com- ponents of the elastic tensor. In practice however, the full tensor 1.1 Surface wave tomography studies cannot be resolved by the seismic data independently at every lo- Surface wave tomography offers a powerful technique to constrain cation, and generally only a restricted number of parameters are seismic anisotropy and to image the structure of the upper man- inverted for. This is done by assuming specific symmetry classes, tle at both regional and global scales. With growing amounts of or by using petrological constraints to impose relations between