Mahalanobis Distance Method for Enhancing Magnetotelluric Data Processing: A Case Study from the Rochechouart Impact Structure, France
نویسندگان
1 Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran.
2 Department of Earth Sciences, Geosciences Paris-Saclay Laboratory (GEOPS), Université Paris-Saclay, Orsay, France.
3 CIRIR – Center for International Research and Restitution on Impacts and on Rochechouart, Rochechouart, France.
4 Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran.
5 Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran.
6 Department of Earth Physics, Institute of Geophysics, University of Tehran, Tehran, Iran.
7 CEREGE, Aix-Marseille Université, CNRS, IRD, INRAE, Aix-en-Provence, France.
doi
10.22059/jesphys.2026.409669.1007750چکیده
One of the primary challenges in magnetotelluric data processing is the presence of noise and outliers (anomalous values). These disturbances often come from human-made sources such as power lines, electronic devices, and nearby infrastructure. They can significantly affect the results related to apparent resistivity and phase, leading to unreliable models of subsurface electrical resistivity.To identify and remove these outlier and noisy components, the Mahalanobis distance method is proposed as an effective solution. This approach — applied here in a four-dimensional feature space comprising the real and imaginary parts of two components of the impedance transfer function — calculates the distance of each data point from the dataset, meanwhile accounting for variances, covariances, and correlations between variables, thereby enabling the detection of anomalous points more accurately than simpler (2D) approaches.In this study, to identify outliers, we applied the Mahalanobis distance method to real data from a single MT station located at the Chassenon Forage site within the Rochechouart impact structure, France. The results demonstrate that this approach not only improves the accuracy of subsequent analyses and enables the extraction of more precise information from subsurface structures, but also reduces processing time by efficiently eliminating contaminated windows before final impedance estimation.
