Chasing the Magma Chamber: MT meets Geodynamics and Seismology – A numerical case study of magmatic plumbing at Oldoinyo Lengai Volcano

How well can geophysical methods image magmatic systems? Geophysical methods are commonly used to image magmatic systems; however, synthetic studies which give insights into the resolution of such methods and their interpretational scope are rare. Gravity anomalies, magnetotelluric, seismological and geodynamical modelling all have a different sensitivity to the rock parameters and are thus likely complementary methods. Our study aims to better understand their interplay by performing joint modelling of a synthetic magmatic system.  Our model setup of a magma chamber is inspired by seismological observations at the Natron plumbing system including active volcano Oldoinyo Lengai within the East African Rift system. The geodynamic modelling is guided by shear-wave velocity anomalies and it is constrained by a large Bouguer gravity anomaly which is modelled by a voxel-based gravity code. It yields the 3D distribution of several geological parameters (pressure, temperature, stress, density, rock type). The parameters are converted into a 3D resistivity distribution. By 3D forward modelling including the topography, synthetic MT transfer functions (phase tensor, induction vectors) are calculated for a rectangular grid of 441 sites covering the area. The variation of geodynamic parameters and/or petrological relations alters the related resistivity distribution and thus yields the sensitivity of MT responses to geodynamic parameters. In turn, MT observations may constrain geodynamic modelling by inverting MT transfer functions. The inversion is performed allowing for the recent seismicity distribution beneath the Natron plumbing system, assuming that active seismic areas are related to enhanced resistivity. The inversion is performed for a realistic distribution (in view of logistic accessibility) of about 40 MT sites.

By combining multiple forward models, this study yields insights into the sensitivity of different observables and thus provides a valuable base on how MT, gravity and seismological observations can help imaging a complex geological setting.