Structure constrained joint inversion of magnetotelluric data and receiver function

Using the quasi Newton method, we jointly invert the magnetotelluric and receiver function data to obtain resistivity model and velocity model with similarity constraint, and a new joint constraint operator is defined to control spatial gradients. Our method accommodates the two models with different properties to evolve towards coupled and similar structure, but still provides enough freedom to model all essential features. The evolution of joint operator and spatial gradient distribution indicates the coupling degree of two different models. The synthetic data and field measurement results show that the joint inversion can improve the coupling degree between resistivity and velocity models compared with the single inversion, and maintain a good misfit at the same time. We also find that the joint inversion significantly improves the resolution of magnetotelluric and reduces the non-uniqueness of magnetotelluric.