A simple method to evaluate the uncertainty of magnetotelluric forward modeling for practical three-dimensional conductivity structure models

Earth, Planets and Space

Technical Report

Kiyoshi Baba¹ 

[1] Earthquake Research Institute, The University of Tokyo, 1-1-1, Yayoi, Bunkyo-Ku, 113-0032, Tokyo, Japan;

关键词: Geomagnetic induction;    Magnetotellurics;    Three-dimensional structure;    Forward modeling;    Uncertainty;   

DOI  :  10.1186/s40623-023-01832-5

 received in 2022-11-15, accepted in 2023-04-22,  发布年份 2023

来源: Springer

【 摘 要 】

The forward calculation of magnetotelluric (MT) responses is generally assumed to be sufficiently accurate compared with typical observational errors in practical modeling and inversion studies. Although the uncertainty of the forward calculation may be examined by comparison with analytical or other numerical solutions for some simple models, such an examination does not guarantee that the uncertainty is similar for more realistic complex structures. In this study, I propose a simple method to evaluate the uncertainty of MT forward modeling for practical three-dimensional (3D) conductivity structure models in a Cartesian coordinate system. The method is based on the idea that the horizontal coordinate system can be selected arbitrarily for a general 3D structure. The synthesized MT responses are ideally identical irrespective of the selection but are different because of the difference in discretization angles, boundary values, and numerical errors. By synthesizing MT responses to the model in several different coordinate systems, the mean, standard deviation, and coefficient of variation can be calculated. These statistics provide quantitative information on how stably the forward calculations synthesize MT responses under the given conditions of the structure model, observation array, periods, numerical algorithm for the forward modeling, and mesh design. The proposed method was applied to two practical situations of seafloor MT arrays in the northwestern Pacific and southern Atlantic and a land MT array in Hokkaido, Japan. The results show that the uncertainty is comparable to real observation errors and is significantly dependent on the MT impedance element, period, site, structure model, and horizontal coordinate system. The uncertainty of the forward calculation should be considered for each element, period, and site to quantitatively evaluate how well a given model explains the data. I propose a new root-mean-square in which the residuals are normalized by both the standard errors of the MT responses observed and synthesized. This would help avoid overfitting data in the inversion analysis by ignoring the uncertainty of the forward calculation. This method is also useful for testing the appropriate selection of the coordinate system and mesh design.Graphical Abstract

【 授权许可】

CC BY   
© The Author(s) 2023

【 图 表 】

41116_2023_36_Article_IEq473.gif

Fig. 2

41116_2023_36_Article_IEq471.gif

41116_2023_36_Article_IEq470.gif

Fig. 1

41116_2023_36_Article_IEq469.gif

41116_2023_36_Article_IEq468.gif

41116_2023_36_Article_IEq467.gif

41116_2023_36_Article_IEq466.gif

Fig. 4

41116_2023_36_Article_IEq465.gif

41116_2023_36_Article_IEq464.gif

Fig. 2

Fig. 1

41116_2023_36_Article_IEq461.gif

41116_2023_36_Article_IEq459.gif

41116_2023_36_Article_IEq458.gif

41116_2023_36_Article_IEq457.gif

41116_2023_36_Article_IEq456.gif

41116_2023_36_Article_IEq455.gif

41116_2023_36_Article_IEq454.gif

41116_2023_36_Article_IEq453.gif

41116_2023_36_Article_IEq452.gif

41116_2023_36_Article_IEq451.gif

41116_2023_36_Article_IEq450.gif

Fig. 2

41116_2023_36_Article_IEq449.gif

41116_2023_36_Article_IEq448.gif

41116_2023_36_Article_IEq447.gif

41116_2023_36_Article_IEq446.gif

40517_2023_255_Article_IEq1.gif

41116_2023_36_Article_IEq445.gif

41116_2023_36_Article_IEq444.gif

Fig. 3

41116_2023_36_Article_IEq443.gif

41116_2023_36_Article_IEq441.gif

Fig. 3

41116_2023_36_Article_IEq439.gif

41116_2023_36_Article_IEq438.gif

41116_2023_36_Article_IEq437.gif

41116_2023_36_Article_IEq436.gif

41116_2023_36_Article_IEq435.gif

41116_2023_36_Article_IEq434.gif

Fig. 3

Fig. 2

41116_2023_36_Article_IEq431.gif

41116_2023_36_Article_IEq430.gif

41116_2023_36_Article_IEq429.gif

Fig. 2

41116_2023_36_Article_IEq427.gif

41116_2023_36_Article_IEq426.gif

41116_2023_36_Article_IEq425.gif

Fig. 1

41116_2023_36_Article_IEq424.gif

41116_2023_36_Article_IEq423.gif

41116_2023_36_Article_IEq422.gif

41116_2023_36_Article_IEq421.gif

41116_2023_36_Article_IEq420.gif

41116_2023_36_Article_IEq419.gif

41116_2023_36_Article_IEq418.gif

41116_2023_36_Article_IEq417.gif

41116_2023_36_Article_IEq416.gif

Fig. 1

41116_2023_36_Article_IEq414.gif

41116_2023_36_Article_IEq413.gif

41116_2023_36_Article_IEq412.gif

41116_2023_36_Article_IEq410.gif

41116_2023_36_Article_IEq408.gif

41116_2023_36_Article_IEq407.gif

41116_2023_36_Article_IEq404.gif

Fig. 3

41116_2023_36_Article_IEq401.gif

Fig. 2

41116_2023_36_Article_IEq398.gif

41116_2023_36_Article_IEq397.gif

41116_2023_36_Article_IEq396.gif

Fig. 1

41116_2023_36_Article_IEq394.gif

41116_2023_36_Article_IEq393.gif

41116_2023_36_Article_IEq392.gif

41116_2023_36_Article_IEq391.gif

41116_2023_36_Article_IEq390.gif

41116_2023_36_Article_IEq389.gif

41116_2023_36_Article_IEq388.gif

41116_2023_36_Article_IEq387.gif

41116_2023_36_Article_IEq386.gif

41116_2023_36_Article_IEq385.gif

41116_2023_36_Article_IEq384.gif

41116_2023_36_Article_IEq383.gif

41116_2023_36_Article_IEq382.gif

Fig. 1

41116_2023_36_Article_IEq380.gif

41116_2023_36_Article_IEq379.gif

41116_2023_36_Article_IEq378.gif

12938_2023_1111_Article_IEq24.gif

41116_2023_36_Article_IEq377.gif

41116_2023_36_Article_IEq376.gif

41116_2023_36_Article_IEq375.gif

41116_2023_36_Article_IEq374.gif

Fig. 5

41116_2023_36_Article_IEq372.gif

41116_2023_36_Article_IEq370.gif

41116_2023_36_Article_IEq369.gif

41116_2023_36_Article_IEq367.gif

41116_2023_36_Article_IEq366.gif

41116_2023_36_Article_IEq365.gif

41116_2023_36_Article_IEq364.gif

41116_2023_36_Article_IEq362.gif

Fig. 4

41116_2023_36_Article_IEq361.gif

41116_2023_36_Article_IEq360.gif

41116_2023_36_Article_IEq359.gif

41116_2023_36_Article_IEq358.gif

Fig. 3

Fig. 2

Fig. 1

41116_2023_36_Article_IEq354.gif

41116_2023_36_Article_IEq353.gif

41116_2023_36_Article_IEq352.gif

41116_2023_36_Article_IEq351.gif

41116_2023_36_Article_IEq350.gif

Fig. 3

41116_2023_36_Article_IEq347.gif

41116_2023_36_Article_IEq346.gif

41116_2023_36_Article_IEq345.gif

41116_2023_36_Article_IEq344.gif

Fig. 5

41116_2023_36_Article_IEq342.gif

41116_2023_36_Article_IEq340.gif

Fig. 1

41116_2023_36_Article_IEq338.gif

40517_2023_256_Article_IEq256.gif

41116_2023_36_Article_IEq337.gif

41116_2023_36_Article_IEq336.gif

41116_2023_36_Article_IEq335.gif

41116_2023_36_Article_IEq334.gif

41116_2023_36_Article_IEq333.gif

41116_2023_36_Article_IEq332.gif

Fig. 4

41116_2023_36_Article_IEq330.gif

Fig. 2

41116_2023_36_Article_IEq328.gif

41116_2023_36_Article_IEq327.gif

41116_2023_36_Article_IEq326.gif

40517_2023_256_Article_IEq244.gif

41116_2023_36_Article_IEq325.gif

41116_2023_36_Article_IEq324.gif

41116_2023_36_Article_IEq323.gif

41116_2023_36_Article_IEq322.gif

41116_2023_36_Article_IEq321.gif

Fig. 3

Fig. 2

41116_2023_36_Article_IEq318.gif

41116_2023_36_Article_IEq317.gif

41116_2023_36_Article_IEq316.gif

41116_2023_36_Article_IEq315.gif

41116_2023_36_Article_IEq314.gif

41116_2023_36_Article_IEq313.gif

41116_2023_36_Article_IEq312.gif

41116_2023_36_Article_IEq311.gif

41116_2023_36_Article_IEq310.gif

41116_2023_36_Article_IEq309.gif

41116_2023_36_Article_IEq308.gif

41116_2023_36_Article_IEq307.gif

Fig. 2

40517_2023_256_Article_IEq225.gif

Fig. 1

41116_2023_36_Article_IEq304.gif

41116_2023_36_Article_IEq302.gif

41116_2023_36_Article_IEq301.gif

Fig. 1

41116_2023_36_Article_IEq297.gif

41116_2023_36_Article_IEq296.gif

41116_2023_36_Article_IEq295.gif

41116_2023_36_Article_IEq294.gif

40517_2023_256_Article_IEq213.gif

41116_2023_36_Article_IEq293.gif

41116_2023_36_Article_IEq292.gif

Fig. 1

41116_2023_36_Article_IEq289.gif

Fig. 1

41116_2023_36_Article_IEq287.gif

41116_2023_36_Article_IEq286.gif

41116_2023_36_Article_IEq285.gif

41116_2023_36_Article_IEq284.gif

41116_2023_36_Article_IEq282.gif

41116_2023_36_Article_IEq281.gif

40517_2023_256_Article_IEq199.gif

41116_2023_36_Article_IEq280.gif

Fig. 2

41116_2023_36_Article_IEq278.gif

41116_2023_36_Article_IEq277.gif

41116_2023_36_Article_IEq276.gif

Fig. 1

41116_2023_36_Article_IEq274.gif

41116_2023_36_Article_IEq273.gif

41116_2023_36_Article_IEq272.gif

41116_2023_36_Article_IEq270.gif

41116_2023_36_Article_IEq268.gif

41116_2023_36_Article_IEq266.gif

41116_2023_36_Article_IEq258.gif

Fig. 1

41116_2023_36_Article_IEq246.gif

41116_2023_36_Article_IEq242.gif

41116_2023_36_Article_IEq237.gif

41116_2023_36_Article_IEq232.gif

41116_2023_36_Article_IEq214.gif

41116_2023_36_Article_IEq205.gif